Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
SYSTEM AND METHOD FOR QUANTUM MICROGRID STATE ESTIMATION
Document Type and Number:
WIPO Patent Application WO/2023/133189
Kind Code:
A1
Abstract:
A method for performing state estimation for a microgrid system includes: obtaining a general quantum state estimation model for swing- bus-contained microgrid systems; establishing a state estimation formulation through quantum language; assigning multiple quantum registers used during a quantum state estimation iteration; generating a quantum-circuit-based state estimation solver for performing quantum phase estimation, controlled rotation, and inverse quantum phase estimation; generating a preconditioned quantum linear sol for an illconditioned state estimation iteration; iteratively optimizing voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, the optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; obtaining a quantum state estimation model for hierarchical-based microgrid systems including droop and/or secondary control modes; establishing an enhanced quantum state estimation framework and generating an enhanced state estimation output for controlling one or more parameters of the microgrid system.

Inventors:
ZHANG PENG (US)
FENG FEI (US)
ZHOU YIFAN (US)
TANG ZEFAN (US)
Application Number:
PCT/US2023/010183
Publication Date:
July 13, 2023
Filing Date:
January 05, 2023
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
UNIV NEW YORK STATE RES FOUND (US)
International Classes:
G06N10/00; G05B15/02; G06N10/20; G06N10/60; G06N10/40
Foreign References:
US20210184495A12021-06-17
Other References:
RAGHAV L. PHANI; KUMAR R. SESHU; RAJU D. KOTESWARA; SINGH ARVIND R.: "Optimal Energy Management of Microgrids Using Quantum Teaching Learning Based Algorithm", IEEE TRANSACTIONS ON SMART GRID, IEEE, USA, vol. 12, no. 6, 24 June 2021 (2021-06-24), USA, pages 4834 - 4842, XP011884040, ISSN: 1949-3053, DOI: 10.1109/TSG.2021.3092283
AKSHAY AJAGEKAR; FENGQI YOU: "Quantum computing for energy systems optimization: Challenges and opportunities", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 29 February 2020 (2020-02-29), 201 Olin Library Cornell University Ithaca, NY 14853 , XP081610982, DOI: 10.1016/j.energy.2019.04.186
MARKO ANGJELICHINOSKI; ANNA SCAGLIONE; PETAR POPOVSKI; CEDOMIR STEFANOVIC: "Decentralized DC MicroGrid Monitoring and Optimization via Primary Control Perturbations", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 26 March 2018 (2018-03-26), 201 Olin Library Cornell University Ithaca, NY 14853 , XP081308678, DOI: 10.1109/TSP.2018.2827331
VAZQUEZ ALMUDENA CARRERA, HIPTMAIR RALF, WOERNER STEFAN: "Enhancing the Quantum Linear Systems Algorithm Using Richardson Extrapolation", ACM TRANSACTIONS ON QUANTUM COMPUTING (PRINT), CORNELL UNIVERSITY LIBRARY, ARXIV.ORG, ITHACA, 5 July 2021 (2021-07-05), Ithaca, pages 1 - 37, XP093078900, Retrieved from the Internet [retrieved on 20230905], DOI: 10.1145/3490631
Attorney, Agent or Firm:
ELLENBOGEN, Wayne et al. (US)
Download PDF:
Claims:
CLAIMS

What is claimed is:

1. A quantum computing-based method for performing enhanced state estimation for a microgrid system, the method comprising: obtaining a general quantum state estimation model for swing-bus-contained microgrid systems; establishing a state estimation formulation through quantum language; assigning a plurality of quantum registers for at least one of quantum encoding, binary representation of eigenvalues, and binary representation of electrical values during a quantum state estimation iteration; generating a quantum-circuit-based state estimation solver configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation; generating a preconditioned quantum linear solver for an ill-conditioned state estimation iteration; iteratively optimizing voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, said optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; obtaining a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; establishing an enhanced quantum state estimation framework including the quantum- circuit-based state estimation solver and the preconditioned quantum linear solver for well- conditioned and ill-conditioned hierarchical-based microgrids, respectively; and utilizing the quantum-circuit-based state estimation solver to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

2. The method according to claim 1, wherein generating a preconditioned quantum linear solver is performed using incomplete Cholesky factorization.

3. The method according to claim 1, wherein assigning a plurality of quantum registers comprises initializing three quantum registers, the three quantum registers being configured for quantum encoding, binary representation of eigenvalues, and binary representation of electrical values, respectively.

4. The method according to claim 1, wherein generating the preconditioned quantum linear solver comprises: initializing ΔV0 , r0, ω0 , and p0 , where ΔV0 represents a difference voltage vector at iteration zero, r0 represents the residual vector at iteration zero, ω0 represents preconditioned error at iteration zero, and p0 refers to search direction at iteration zero; for each iteration, updating ΔVξ+1 and rξ+1 using the expression f

1 representing a preconditioned quantum linear solver iteration number, is a coefficient, is a gain matrix; and when rξ+1 reaches a prescribed tolerance of e, outputting ΔVξ+1 , otherwise updating

\ respectively.

5. The method according to claim 1, further comprising incorporating at least one of droop and secondary controls into a vector, h(V), containing microgrid system states to be estimated, using enhanced quantum state estimation (EQSE).

6. The method according to claim 5, wherein when only droop controls are applied in the microgrid system, EQSE is configured to devise a vector of power injections at unknown- voltage buses, P(V), and corresponding elements in a Jacobian matrix, as follows: where G(diag) is a diagonal part of a nodal conductance matrix, G, of the microgrid system, KG is a vector containing a reciprocal of each power and voltage droop coefficient, and Vref is a reference- voltage error.

7. The method according to claim 5, wherein when both droop and secondary controls are applied in the microgrid system, the method comprises assisting voltage recovery by adding a dummy bus, Vd, to determine a vector of power injections at unknown-voltage buses, P(V), as where V is an unknown-voltage vector, V* is a rated-voltage vector, and denotes the dummy bus voltage vector at a previous iteration.

8. The method according to claim 7, wherein at each iteration, Vd is updated until a difference between V* and V reaches convergence.

9. The method according to claim 1, wherein the quantum-circuit-based state estimation solver comprises a Harrow-Hassidim-Lloyd (HHL)-based state estimation solver.

10. The method according to claim 1, wherein iteratively optimizing the voltage and residual vectors is performed using a gradient descent rule.

11. A quantum computing-based apparatus for performing enhanced state estimation for a microgrid system, the apparatus comprising: a control module including at least one classical processor configured to update a quantum computing model to obtain eigenvectors and corresponding coefficients of the eigenvectors of a gain matrix and cost function of the quantum computing model; and a quantum-circuit-based state estimation solver operatively coupled to the control module, the quantum-circuit-based state estimation solver comprising at least one quantum processing unit configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation, the quantum-circuit-based state estimation solver further comprising a plurality of quantum registers for at least one of quantum encoding, binary representation of eigenvalues, and binary representation of electrical values during a quantum state estimation iteration; wherein the at least one quantum processing unit is configured: to generate a preconditioned quantum linear solver for an ill-conditioned state estimation iteration; to iteratively optimize voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, said optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; to obtain a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; to establish an enhanced quantum state estimation framework including the quantum-circuit-based state estimation solver and the preconditioned quantum linear solver for well-conditioned and ill-conditioned hierarchical-based microgrids, respectively; and to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

12. The apparatus according to claim 11, wherein the preconditioned quantum linear solver is configured to performed incomplete Cholesky factorization.

13. The apparatus according to claim 11, wherein the preconditioned quantum linear solver is configured: to initialize ΔV0, r0, ω0 , and p0, where ΔV0 represents a difference voltage vector at iteration zero, r0 represents the residual vector at iteration zero, ω0 represents preconditioned error at iteration zero, and p0 refers to search direction at iteration zero; for each iteration, to update ΔVξ+1 and rξ+1 using the expression . ’ where ξ is an integer representing a preconditioned quantum linear solver iteration number, is a coefficient, is a gain matrix; and when rξ+1 reaches a prescribed tolerance of e, to output ΔVξ+1 , otherwise to update

\ respectively.

14. The apparatus according to claim 11, wherein the quantum-circuit-based state estimation solver is configured to incorporate at least one of droop and secondary controls into a vector, h(V), containing microgrid system states to be estimated, using enhanced quantum state estimation (EQSE).

15. The apparatus according to claim 14, wherein when only droop controls are applied in the microgrid system, EQSE comprises devising a vector of power injections at unknown-voltage buses, P(V), and corresponding elements in a Jacobian matrix, as follows: f where G(diag) is a diagonal part of a nodal conductance matrix, G, of the microgrid system, KG is a vector containing a reciprocal of each power and voltage droop coefficient, and Vref is a reference- voltage vector.

16. The apparatus according to claim 14, wherein when both droop and secondary controls are applied in the microgrid system, the quantum-circuit-based state estimation solver is configured to assist voltage recovery by adding a dummy bus, Vd, to determine a vector of power injections at unknown-voltage buses, P(V), as where V is an unknown-voltage vector, V* is a rated-voltage vector, and Vpd denotes the dummy bus voltage vector at a previous iteration.

17. The apparatus according to claim 16, wherein at each iteration, Vd is updated until a difference between V* and V reaches convergence.

18. The apparatus according to claim 11, wherein the quantum-circuit-based state estimation solver comprises a Harrow-Hassidim-Lloyd (HHL)-based state estimation solver.

19. The apparatus according to claim 11, wherein the at least one quantum processing unit is configured to iteratively optimize the voltage and residual vectors using a gradient descent rule.

20. An apparatus for performing quantum state estimation for a microgrid system, the apparatus comprising: at least one classical processing unit configured to obtain a general quantum state estimation model for swing-bus-contained microgrid systems and to update the quantum state estimation model to obtain eigenvectors and corresponding coefficients of the eigenvectors of a gain matrix and cost function of the quantum state estimation model; and at least one quantum processing unit operatively coupled to the at least one classical processing unit, the quantum processing unit being configured: to establish a state estimation formulation through quantum language; to assign a plurality of quantum registers for at least one of quantum encoding, binary representation of eigenvalues, and binary representation of electrical values during a quantum state estimation iteration; to generate a quantum-circuit-based state estimation solver configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation; to generate a preconditioned quantum linear solver for an ill-conditioned state estimation iteration; to iteratively optimize voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, said optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; to obtain a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; to establish an enhanced quantum state estimation framework including the quantum-circuit-based state estimation solver and the preconditioned quantum linear solver for well-conditioned and ill-conditioned hierarchical-based microgrids, respectively; and to utilize the quantum-circuit-based state estimation solver to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

Description:
SYSTEM AND METHOD FOR QUANTUM MICROGRID STATE ESTIMATION CROSS-REFERENCE TO RELATED APPLICATIONS [0001] This application claims the benefit of U.S. Provisional Patent Application No. 63/296,604, filed on January 5, 2022, entitled “System and Method for Quantum Microgrid State Estimation,” the disclosure of which is incorporated by reference herein in its entirety for all purposes. STATEMENT OF GOVERNMENT RIGHTS [0002] This invention was made with government support under contract number OIA-2040599 awarded by the National Science Foundation. The government has certain rights in the invention. BACKGROUND [0003] The present invention relates generally to the electrical, electronic and computer arts, and, more particularly, to a system and method for microgrid state estimation. [0004] Microgrid is a proven paradigm that can flexibly manage distributed energy resources (DERs) and ensure the electricity resiliency against outages. Among many microgrid functions, state estimation is of fundamental importance, as it provides valuable resources for online monitoring and control of microgrids. Basic requirements for microgrid state estimation primarily include accuracy, efficiency and resiliency against noises (see, e.g., A. Primadianto et al., “A Review on Distribution System State Estimation,” IEEE Trans. on Power Systems, vol.32, no.5, pp.3875-3883, 2017, the disclosure of which is incorporated by reference herein in its entirety). [0005] State estimation is an important element of energy management systems, at least in part because it allows a system operator to control and monitor prescribed parameters of an electric power system. The state of the power system can be estimated based on supervisory control and data acquisition (SCADA) measurements, phasor measurement unit (PMU) measurements, and other measurements known in the art. Some of the most important functions of EMS, such as, but not limited to, contingency analysis, optimal power flow, and security assessment, are based in large part on state estimation. [0006] For modern microgrids, an increasingly urgent and important demand for high-frequency state estimation is caused, at least in part, by community expansion, a high penetration of uncertain renewables, and volatile operational conditions, among other factors. However, the complexities of conventional state estimation methods scale polynomially with the problem size, which makes those methods no longer suitable for a future grid with formidable real-time operation requirements.

[0007] Quantum computing provides a promising solution for overcoming the complexity issue. Unlike classical computing, quantum computing requires fewer bits (i.e., qubits) to handle a complicated problem. For microgrid state estimation, one main bottleneck is to establish an efficient solver for a sparse linear system of equations. Currently, there are two main types of quantum linear system algorithms: hybrid quantum/classical algorithms; and quantum-circuit- based algorithms. Hybrid algorithms are developed for the noisy intermediate-scale quantum (NISQ) era. Examples include the Variational Quantum Linear Solver and quantum random walk algorithms. However, the inevitable computation and optimization processes on classical computers have made most hybrid algorithms heuristic. As a consequence, the convergence and efficiency of hybrid algorithms cannot be strictly guaranteed.

[0008] In addition to the complexity issue, high- volume state measurement and information exchange also make microgrids vulnerable to exogenous disturbances (see, e.g., S. Bi, et al., “Graphical Methods for Defense Against False-Data Injection Attacks on Power System State Estimation,” IEEE Trans, on Smart Grid, vol. 5, no. 3, pp. 1216-1227, 2014, the disclosure of which is incorporated by reference herein in its entirety). However, conventional microgrid state estimation methods typically follow classical power flow formulations, which have ignored the features of droop-based microgrids where multiple DERs, instead of a swing bus, support the entire system. For hybrid AC/DC microgrids, an attempted solution to the state estimation problem has conventionally involved decomposing the system into subsystems. In those existing methods, microgrids are modeled in the same way with the traditional distribution feeders, where a main grid or an infinite source is used to support the downstream system. Failure to represent the droop/secondary regulation may provide erroneous results when a microgrid is subjected to unenforceable disturbances. SUMMARY

[0009] Principles of the present invention, as manifested in one or more embodiments thereof, are directed to a system and method for performing enhanced quantum microgrid state estimation that addresses at least the above-noted problems with conventional (i.e., classical) state estimation approaches.

[0010] In accordance with an embodiment of the invention, a quantum computing-based method for performing enhanced state estimation for a microgrid system includes: obtaining a general quantum state estimation model for swing-bus-contained microgrid systems; establishing a state estimation formulation through quantum language; assigning a plurality of quantum registers for at least one of quantum encoding, binary representation of eigenvalues, and binary representation of electrical values during a quantum state estimation iteration; generating a quantum-circuit-based state estimation solver configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation; generating a preconditioned quantum linear solver for an ill- conditioned state estimation iteration; optimizing voltage and residual vectors iteratively, using a gradient descent rule, while iteratively updating preconditioned errors by the quantum-circuit- based state estimation solver, such optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; obtaining a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; establishing an enhanced quantum state estimation framework including the quantum-circuit-based state estimation solver and the preconditioned quantum linear solver for well-conditioned and ill-conditioned hierarchical-based microgrids, respectively; and utilizing the quantum-circuit-based state estimation solver to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

[0011] In accordance with another embodiment of the invention, an apparatus for performing quantum microgrid state estimation in conjunction with a hierarchical control-based microgrid comprises a control module including at least one classical processor, and a quantum-circuit-based state estimation solver operatively coupled to the control module and including at least one quantum processing unit. The classical processor is configured to update a quantum computing model to obtain eigenvectors and corresponding coefficients of the eigenvectors of a gain matrix and cost function of the quantum computing model. The quantum processing unit is configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation. The quantum-circuit-based state estimation solver further includes a plurality of quantum registers for quantum encoding, binary representation of eigenvalues, and/or binary representation of electrical values during a quantum state estimation iteration.

[0012] In one or more embodiments, the quantum processing unit is configured: to generate a preconditioned quantum linear solver for an ill-conditioned state estimation iteration; to iteratively optimize voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, said optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; to obtain a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; to establish an enhanced quantum state estimation framework including the quantum-circuit-based state estimation solver and the preconditioned quantum linear solver for well-conditioned and ill-conditioned hierarchical-based microgrids, respectively; and to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

[0013] In accordance with yet another embodiment of the invention, an apparatus for performing quantum state estimation for a microgrid system includes at least one classical processing unit configured to obtain a general quantum state estimation model for swing -bus-contained microgrid systems and to update the quantum state estimation model to obtain eigenvectors and corresponding coefficients of the eigenvectors of a gain matrix and cost function of the quantum state estimation model, and at least one quantum processing unit operatively coupled to the classical processing unit. The quantum processing unit is configured: to establish a state estimation formulation through quantum language; to assign a plurality of quantum registers for at least one of quantum encoding, binary representation of eigenvalues, and binary representation of electrical values during a quantum state estimation iteration; to generate a quantum-circuit-based state estimation solver configured to perform quantum phase estimation, controlled rotation, and inverse quantum phase estimation; to generate a preconditioned quantum linear solver for an ill- conditioned state estimation iteration; to iteratively optimize voltage and residual vectors while iteratively updating preconditioned errors by the quantum-circuit-based state estimation solver, said optimizing and updating iterations continuing until the preconditioned quantum linear solver reaches convergence; to obtain a quantum state estimation model for hierarchical-based microgrid systems including at least one of droop and secondary control modes; to establish an enhanced quantum state estimation framework including the quantum-circuit-based state estimation solver and the preconditioned quantum linear solver for well-conditioned and ill-conditioned hierarchical-based microgrids, respectively; and to utilize the quantum-circuit-based state estimation solver to generate an enhanced state estimation output for controlling one or more parameters of the microgrid system.

[0014] As the term may be used herein, “facilitating” an action contemplates performing the action, making the action easier, helping to carry out the action, or causing the action to be performed. Thus, by way of example only and without limitation, in the context of a processor- implemented method, instructions executing on one processor might facilitate an action carried out by instructions executing on a remote processor, by sending appropriate data or commands to cause or aid the action to be performed. For the avoidance of doubt, where an actor facilitates an action by other than directly performing the action itself, it is assumed that the action is nevertheless performed by some entity or combination of entities.

[0015] One or more embodiments of the invention or elements thereof may be implemented in the form of a computer program product including a computer readable storage medium with computer usable program code for performing the method steps indicated. Furthermore, one or more embodiments of the invention or elements thereof can be implemented in the form of a system (or apparatus) including a memory, and at least one processor that is coupled to the memory and operative to perform exemplary method steps.

[0016] Yet further, in another aspect, one or more embodiments of the invention or elements thereof can be implemented in the form of means for carrying out one or more of the method steps described herein; the means can include (i) hardware module(s), (ii) software module(s) stored in a computer readable storage medium (or multiple such media) and implemented on a hardware processor, or (iii) a combination of (i) and (ii); any of (i)-(iii) implement the specific techniques set forth herein.

[0017] Techniques of the present invention can provide substantial beneficial technical effects. By way of example only and without limitation, a quantum microgrid state estimation system and method according to one or more embodiments of the invention may provide one or more of the following features, among other advantages: utilizes quantum-circuit based algorithms that achieve exponential speedup in microgrid state estimation over traditional methodologies; includes a preconditioned quantum linear solver (PQLS) for handling ill-conditioned general quantum state estimation (GQSE) with limited quantum resources; and establishes an enhanced quantum state estimation (EQSE) algorithm for hierarchical control-based microgrids with exogenous disturbances.

[0018] These and other features and advantages of the present invention will become apparent from the following detailed description of illustrative embodiments thereof, which is to be read in connection with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Non-limiting and non-exhau stive embodiments of the present invention will be described with reference to the following drawings which are presented by way of example only, wherein like reference numerals (when used) indicate corresponding elements throughout the several views unless otherwise specified, and wherein:

[0020] FIG. 1 is a block diagram depicting at least a portion of an exemplary architecture of a quantum-circuit-based general quantum state estimation (GQSE) algorithm, according to one or more embodiments of the present invention;

[0021] FIG. 2 depicts an exemplary enhanced quantum state estimation (EQSE) methodology suitable for performing DC microgrid state estimation, according to one or more embodiments of the present invention;

[0022] FIGS. 3 A and 3B are block diagrams depicting at least portions of exemplary test system architectures for a two-bus DC microgrid system and a four-bus DC microgrid system, respectively, in which embodiments of the present invention can be utilized;

[0023] FIGS. 4A and 4B are graphs depicting exemplary voltage distributions of droop-based enhanced quantum state estimation (EQSE) against disturbances o. N (0.01, 0.01 2 ) and N(0.01 0.05 2 ), respectively, according to illustrative embodiments of the present invention; [0024] FIGS. 5 A and 5B are graphs depicting exemplary voltage distributions of secondarybased EQSE against disturbances of MO.01. 0.01 2 ) and MO.01 , 0.05 2 ), respectively, according to illustrative embodiments of the present invention;

[0025] FIG. 6 is a bar graph depicting exemplary convergence performance results of EQSE according to embodiments of the present invention compared with classical state estimation (CSE); and

[0026] FIG. 7 is a block diagram depicting at least a portion of an exemplary computer system configurable for implementing at least a portion of one or more aspects and/or elements of the present invention.

[0027] It is to be appreciated that elements in the figures are illustrated for simplicity and clarity. Common but well-understood elements that may be useful or necessary in a commercially feasible embodiment may not be shown in order to facilitate a less hindered view of the illustrated embodiments.

DETAILED DESCRIPTION

[0028] Principles of the present invention, as manifested in one or more embodiments thereof, will be described herein in the context of an illustrative system, apparatus and/or method for direct current (DC) microgrid state estimation. It is to be appreciated, however, that the invention is not limited to the specific system, apparatus and/or methods illustratively shown and described herein. Rather, it will become apparent to those skilled in the art given the teachings herein that numerous modifications can be made to the embodiments shown that are within the scope of the claimed invention. For example, although the formulations described herein are presented in the context of a DC microgrid, it is to be appreciated that these formulations can be readily extended to alternating current (AC) microgrids and other hierarchical control-based microgrids as well. That is, no limitations with respect to the embodiments shown and described herein are intended or should be inferred.

[0029] A microgrid is a grouping of distributed energy resources (DERs) with control capabilities and can work in conjunction with the main/primary grid or operate autonomously as dictated by demand and/or other considerations (e.g., cost, etc.). It is essentially a smaller version of the grid with several distinct advantages. These advantages include the capability to operate when the main grid is down, a provision of stability in strengthening the grid by reducing grid disruptions, and an increase in efficiency due to the use of local energy sources which mitigate energy losses in the processes of transmission and distribution. Microgrids often use localized renewable energy sources for power generation (e.g., solar panel arrays, wind turbines, etc.), thus making them more environmentally friendly. A microgrid is mostly automated with smart technology and intelligent controls that can anticipate problems and failures and reconfigure itself to account for them.

[0030] As previously asserted, microgrid state estimation is an important aspect of energy management systems (EMS), primarily since it allows a system operator to control and monitor one or more parameters of an electric power system. The state of the power system can be estimated based on supervisory control and data acquisition (SCADA) measurements, phasor measurement unit (PMU) measurements, and/or other measurements known in the art. Some of the most important functions of EMS, such as, but not limited to, contingency analysis, optimal power flow, and security assessment, are based in large part on state estimation.

[0031] For modern microgrids, the demand for high-frequency state estimation is attributable to community expansion, a high penetration of uncertain renewables, and volatile operational conditions, among other factors. However, the complexities of conventional state estimation methods scale polynomially with the problem size, which makes those methods unsuitable for a future grid with formidable real-time operation requirements.

[0032] In order to address one or more problems associated with conventional microgrid state estimation, aspects of the present invention as described herein are directed to quantum-circuit- based algorithms for microgrid state estimation. In one or more embodiments, a general quantum state estimation (GQSE) formulation for swing-bus-constrained microgrids is devised through a quantized Gaussian-Newton iteration. To resolve the ill-conditioned matrix issue in GQSE, a preconditioned quantum linear solver (PQLS) is developed with limited quantum resources. Further, one or more embodiments of the invention establish an enhanced quantum state estimation (EQSE) algorithm for hierarchical-control-based microgrids with exogenous disturbances.

[0033] State estimation is an indispensable functionality in modern EMS for determining power system states with raw measurement data containing heterogeneous noises. A weighted least square can be used for state estimation. It minimizes the sum of weighted squared errors between measurements and estimated system states. The following equation gives a cost function of state estimation for a DC microgrid:

7(v) = [z - h<r[z - ft(v)] , (1) where , and the initialized V is a vector where each element is one). h(V) is a vector containing the system states to be estimated. In this illustrative scenario, the estimated system states include three types of states, namely, power injections at unknown-voltage buses, branch powers at unknown-voltage buses, and unknown voltages. Mathematically, it can be expressed as follows:

GV ° V h(V) = V ( ) ° ° (/( ) - ^(Z)) (2) where g, denotes the branchpower vector, V g and V/ are from-bus-voltage and to-bus voltage vectors, respectively, G and g (g,/) are the nodal conductance matrix and branch conductance vector of the DC microgrid, respectively, and ° represents a Hadamard product.

[0034] In equation (1) above, z is a vector containing the measured results of variables in h(V). R is the weight matrix (each element in R represents the coefficient for the corresponding element in [z - h(V)]).

[0035] The Gauss-Newton algorithm can be applied to solving the nonlinear least square problem in equation (1). It performs the following iterations to obtain a minimal cost function:

H T R~\Z - , (3) where - , k or at the k th iteration containing the difference between V k and V k-1 , and H is a Jacobian matrix of A( Vzc-i). Then, the bus-voltage vector Vk can be updated as V k = i + ΔV k .

[0036] A major computation burde Newton algorithm lies in equation (3) above, which involves the calculation of the gain matrix inverse; that is, (H T R -1 H) -1 . Meanwhile, the gain matrix inherits the sparsity feature from the Jacobian matrix H and is strictly

Hermitian ws the classical state estimation to be converted to a quantum formulation. [0037] In the GQSE formulation, the Harrow-Hassidim-Lloyd (HHL) algorithm is utilized to handle equation (3) above to reduce computational complexity of the entire state estimation. For conciseness, the following derivations omit the subscript k. The classical formulation shown in equation (3) can be converted to a quantum computing model as follows: where ΔV represents a difference of bus voltage vector between current iteration and previous one iteration, and |ΔV) is a normalized quantum state of ΔV. More particularly, the i th component of p amplitude of the i th basis state of the quantum state | ΔV). Similarly, is a normalized quantum state of

[0038] A primary objective of GQSE is to prepare a quantum superposition of ΔV, i.e., |ΔV), on the quantum circuit, which satisfies equation (4). To achieve this, since is Hermitian, it has a spectral decomposition as follows: where Uj) is the j th eigenvector respective eigenvalue λj. \H T R-1( z -h(V)) can also be decomposed using the eigenvectors of as follows: where bj is the coefficient of eigenvector | u j ).

[0039] Consequently, the quantum superposition of ΔV, i.e., |ΔV), can be established from equations (4) - (6) above as follows:

|

[0040] Equation (7) above allows the expression of |ΔV) using |u j ). A primary task is to establish a quantum circuit where A is combined with bj |u j). In one or more embodiments, the quantum-circuit-based GQSE algorithm is developed to achieve this task, as described in further detail herein below. [0041] FIG. 1 is a block diagram depicting at least a portion of an exemplary architecture 100 for implementing a quantum-circuit-based GQSE algorithm, according to one or more embodiments of the invention. With reference to FIG. 1, the architecture 100 includes a control module or center 102, including at least one central processing unit (CPU), operatively coupled to an HHL quantum circuit 104, which together are configured to monitor disturbances 106 on an electrical network or microgrid. An overall idea of the GQSE algorithm is that raw data z, which is the vector containing measured results of variables in h(V), the vector containing the system states to be estimated, are measured and sent to the control center 102. The CPU in the control center 102 is configured to update using equations (5) and (6) above, to obtain | and bj, respectively. This process is accomplished in at least one classical processor or computer (e.g., CPU), in one or more embodiments.

[0042] The HHL quantum circuit 104, which preferably includes at least one quantum processing unit (QPU), is then configured to obtain | ΔV). In one or more embodiments, the HHL quantum circuit 104 includes three components: namely, a quantum phase estimation (QPE) module 108, a controlled rotation module 110, and an inverse QPE (QPE ) module 112. The input to the HHL quantum circuit 104 comprises a qubits for ancilla quantum encoding (AQE), denoted as anc, n qubits for the binary representation of λj, denoted as q, and m qubits for the binary representation of hy|uy), denoted as io, as shown in FIG. 1, where a, n and m are integers. The initialized states of the qubits in anc and q are all zero. For io, the initialized states are bj |u j ) (j = 1, 2, . . . , m). In other words, the input of the HHL quantum circuit 104 can be represented as |φ 0 ) = \ , where ® refers to a tensor product.

[0043] The initialized quantum state |φ 0 ) is first handled by QPE. In one or more embodiments, the QPE module 108 comprises n Hadamard gates configured for generating quantum superpositions, n quantum unitary gates configured for obtaining a multi-qubit quantum state containing λj , and an inverse quantum Fourier transformation (QFT) module, denoted QFT^, configured for obtaining a reciprocal of each λj , i.e., λj -1 . λj -1 and bj |u j ) are thus combined through the functions of the QPE module 108 and the controlled rotation module 110. The inverse QPE module 112 is further utilized to reset the state of each qubit in q and io, i.e., making the state of each qubit in q to be zero and that in io as bj |u j )• [0044] With continued reference to FIG. 1, the QPE module 108 aims to obtain λj . Initially, n Hadamard gates are applied on n qubits in q, respectively, to create quantum superpositions. The quantum state (with all the qubits used) after all the Hadamard gates are applied (denoted as |φ1) in FIG. 1), can be expressed as follows: i where H and I are the Hadamard and unit gates, respectively.

[0045] n quantum unitary gates, represented as (where k = 0, 1, . . . n-1), are then applied on the qubits in io in series for obtaining a multi-qubit quantum state containing λj . Specifically, Afterunitary gates arc applied, the quantum state (|φ 2 )) can be expressed as follows:

[0046] Inverse QFT is further applied to obtain each λj . The quantum states after applying inverse QFT ( | φ 3 )) can be represented as

|

[0047] With continued reference to FIG. 1, the controlled rotation module 110 is configured to apply AQE to thereby obtain the reciprocal of each λj, i.e., λ The quantum state | φ 4 )) than contains , i.e., |ΔV) as in equation (7), as follows:

I where fl is a user-defined constant value. [0048] Referring again to FIG. 1, the inverse QPE module 112 is utilized to reset the state of each qubit in q and i o . After inverse quantum phase estimation is applied, the quantum state | φ 5 )) can be expressed as follows:

[0049] As shown in equation (12) above, once the measurement on anc is | l) a , the quantum state can be obtained as: where ( see equation

(7)). The voltage vector V can thus be updated for the next iteration. The GQSE iterations continue until ΔV achieves convergence.

[0050] Although theoretically HHL offers an exponential speedup in terms of the system dimension, its efficiency can be affected by several factors. A computational complexity of HHL can be expressed as where N, η and K denote the dimension, sparsity and condition number, respectively, of the gain matrix While q does not vary much due to the sparse nature of power grids, K can largely deviate in different systems. If an ill-conditioned gain matrix (where K is extremely large) appears in GQSE, the efficiency of GQSE will be largely impacted. This is because a larger K leads to a higher difference between the largest and the smallest eigenvalues, which makes the QPE in GQSE require more qubits to maintain the estimation accuracy. This inevitably leads to an increased quantum circuit depth and an overcompensation of qubit resources, which is undesirable.

[0051] To address this challenge, the present disclosure develops and provides a preconditioned quantum linear solver (PQLS). A primary objective of PQLS, according to one or more embodiments of the invention, is to use a preconditioned iterative optimization to obtain ΔV instead of directly calculating ΔV using equation (4) above. In other words, at each iteration in the Gauss-Newton algorithm (see equations (3) and (4)), ΔV is calculated by PQLS which utilizes an optimization process. Iteration procedures according to one or more embodiments are described below.

[0052] Multiple iterations are involved in PQLS. The symbol ξ is used to distinguish each iteration in PQLS with the iteration in the Gauss-Newton algorithm (denoted as K in equation (3)), where f starts at 0, and ends when a convergence is reached. At the iteration of PQLS, ΔV (+1 and the residual vector obtained from equation (3)) can be updated by using a gradient descent rule as follows: where the initialized (H T R -1 H)ΔV 0 , respectively, p ξ refers to the search direction and can be updated as follows: represents preconditioned errors. It can be calculated through quantum computing as follows:

\ where M denotes a positive-definite, fixed preconditioner. An example of a commonly used preconditioner is the incomplete Cholesky factorization (see, e.g., P. Concus, et al., “Block Preconditioning for the Conjugate Gradient Method,” SIAM Journal on Scientific and Statistical Computing, vol. 6, no. 1, pp. 220-252, 1985, the disclosure of which is incorporated by reference herein in its entirety), is initialized as follows:

\

[0053] p ξ is initialized as In equation (13), is a coefficient, and can be obtained as ( ^

[0054] The entire PQLS algorithm can be summarized as follows:

• Step 1: Initialize ΔV 0 , r 0 , ω 0 , and p 0 .

• Step 2: For each iteration, update ΔV ξ+1 and r ξ+1 using equation (13).

• Step 3: If r ξ+1 reaches a tolerance of e, output ΔV^ +1 . Otherwise, update and

Pt+i using equations (15) and (14), respectively. [0055] In micro grids, hierarchical controls are applied on DERs to support power consumption and regulate bus voltages. A hierarchical control commonly consists of a droop control and a secondary control. For DC microgrids, an active power (P) and voltage (V), referred to as P/V, droop control can be applied to balance the power consumption. A main function of the secondary control, in one or more aspects of the invention, is to retrieve deviations of local voltages to their nominal values.

[0056] Conventional microgrid state estimation methods typically follow classical power flow formulations in equation (2), which however have ignored the features of hierarchical-control- based microgrids where multiple DERs, instead of a swing bus, support the entire system. Failure to represent the droop/secondary regulation may provide erroneous results when a microgrid is subjected to disturbances. Thus, embodiments of the invention, which utilize a novel EQSE algorithm, provide a beneficial improvement in the technical art of microgrid state estimation over conventional methodologies.

[0057] FIG. 2 depicts an exemplary EQSE methodology 200 suitable for use in a DC microgrid, according to one or more embodiments of the invention. The EQSE methodology 200 is a symbolic representation of the EQSE algorithm previously described herein, showing a correlation with equations (2, 4, and 13-18), as shown in FIG. 2. The EQSE methodology 200 incorporates either droop or (droop + secondary) control into h(V) to mitigate errors of state estimation results when the measurement z contains disturbances.

[0058] For droop-based EQSE, when only droop controls are applied in a microgrid, P(V) and the corresponding elements in the Jacobian matrix can be devised as

(

] (17)

( where G(diag) is the diagonal part of G, KG is a vector containing the reciprocal of each P/V droop coefficient, and V ref is the reference- voltage error.

[0059] For secondary-based EQSE, when both droop and secondary controls are applied in a microgrid, a dummy bus V d can be added to calculate P(V) to assist the voltage recovery. Specifically, Vd can be designed as V where V is the unknown- voltage vector, V* is the rated-voltage vector, and V denotes the dummy-bus -voltage vector at the previous iteration. To be specific, it the element within (V* — V) is negative (or positive), the corresponding element in V will be added by a positive (or negative) deviation. With the dummy bus V d , P(V) and d can be expressed as where G d is a vector containing the coefficients. At each iteration, V d will be updated until the difference between V* and V reaches convergence.

[0060] It is to be understood that the aforementioned exemplary droop/secondary-based state estimation schemes according to embodiments of the invention are described in the context of DC microgrids by way of illustration only and without limitation or loss of generality. Nonetheless, the principles of EQSE according to embodiments of the invention are not limited to DC microgrids but rather are generic and can be readily applicable for AC microgrids as well, where hierarchical P//and Q/V control schemes are adopted, as will become apparent to those skilled in the art given the teachings herein.

Numerical Tests - Case Studies

[0061] Several case studies were conducted for validating the correctness of GQSE, PQLS and EQSE in two typical illustrative microgrids; namely, a two-bus DC microgrid and a four-bus DC microgrid. Specifically, FIGS. 3 A and 3B are block diagrams depicting at least portions of exemplary test system architectures for a two-bus DC microgrid system 300 and a four-bus DC microgrid system 350, respectively, in which embodiments of the invention can be utilized. With reference to FIG. 3 A, the two-bus DC microgrid system 300 includes a first bus (Bus 1), which is a swing bus, and a second bus (Bus 2), which is a load bus, are operatively coupled together through a resistive element (R) 302. With reference to FIG. 3B, the four-bus DC microgrid system 350 includes a first bus (Bus 1), to which a first DER (DERI) is connected, a second bus (Bus 2), which is a first load bus (Loadl), a third bus (Bus 3), which is a second load bus (Load2), and a fourth bus (Bus 4), to which a second DER (DER2) is connected. The Bus 1 and Bus 2 are operatively coupled together via a first impedance element 352, Bus 2 and Bus 3 are operatively coupled together via a second impedance element 354, Bus 3 and Bus 4 are operatively coupled together via a third impedance element 356, and Bus 1 and Bus 4 are operatively coupled together via a fourth impedance element 358.

[0062] The robustness and convergence performance of EQSE are verified as well. By way of example only and without limitation, all quantum algorithms are implemented in IBM’s Qisket (version: 0.23.4), Terra (version: 0.16.3), and IBMQ provider (version: 0.11.1). The classical state estimation (CSE), i.e., the Gauss-Newton algorithm, is implemented in MATLAB® (a registered trademark of The MathWorks, Inc.) running on a 2.5 GHz computer for comparison purposes.

[0063] The base voltage, i.e., the nominal rated voltage of the system, in this example is set at 400 V, and the base power is 1 KVA. In the two-bus microgrid system 300 shown in FIG. 3A, assume Bus 1 is a swing bus whose per-unit bus voltage is fixed at 1 per-unit (p.u.), and supports the power consumption of the entire system. In the four-bus microgrid system 350 shown in FIG. 3B, DERs 1 and 2 are connected with Buses 1 and 4, respectively. Assume that the P/V droop coefficients for DERs 1 and 2 are set at 1 x 10’ 5 and 2 x 10’ 4 , respectively; each element in Ga (see equation (18)) for the secondary control is set at 20.

[0064] The correctness of GQSE in the 2-bus DC microgrid system and that of PQLS and EQSE in the 4-bus DC microgrid system will now be shown, using CSE as a comparison. Tables I and III below provide state estimation results of GQSE and PQSL (together with CSE comparisons) in the normal situation (i.e., without an ill-conditioned matrix), respectively. Table II below provides the state estimation results of GQSE with different numbers of qubits in q (see FIG. 1). Table IV presents the PQLS and CSE estimation results with an ill-conditioned matrix (i.e., the maximum eigenvalue: 9.4675 x 10 10 , and the minimum eigenvalue: 277.7778).

Table I: GQSE and CSE results in normal situation (p.u.)

Table II: GQSE results in 1 st iteration with different numbers of qubits in q (p.u.)

Table III: PQLS and CSE results in the normal situation (p.u.)

Table IV: GQSE, PQLS and CSE results with an ill-conditioned matrix (p.u.)

Table V: EQSE and ESE results (p.u.)

[0065] For EQSE, an enhanced state estimation (ESE) method is utilized for comparison. Specifically, ESE has the same principle as EQSE except that the quantum procedures in EQSE are replaced with classical operations. Table V shows the comparison results of EQSE and ESE under droop control and droop plus secondary control. From Tables I-V above, the following insights can be obtained:

(i) The estimation results from GQSE (i.e., after the third iteration) and CSE are exactly the same (see Table I), which validates the correctness of GQSE.

(ii) For GQSE, having a sufficient number of qubits in q is critical to achieve an accurate estimation result. As shown in Table II, increasing the number of qubits in q can greatly improve accuracy of the GQSE result, i.e., the result obtained from GQSE is closer to the CSE result. This is primarily because more precise eigenvalues can be obtained with more qubits in q.

(iii) Tables III and IV validate the correctness of PQLS in the normal situation and with an ill- conditioned matrix, respectively. Results show that although is ill-conditioned, the

PQLS-based QSE still possesses excellent convergence by adjusting the preconditioning errors to eventually achieve accurate results. (iv) Table V validates the correctness of EQSE, i.e., the results obtained from EQSE are exactly the same with those from ESE in both droop and secondary controls.

[0066] The robustness of EQSE against disturbances will now be discussed. Disturbances are assumed to follow a Gaussian distribution with N(μ , σ 2 ). More particularly, by way of illustration only and without limitation, 100 samples of disturbances are generated based on the distribution of M0.01 , 0.01 2 ). These disturbances are then added to the measurement z, respectively. For each disturbance, EQSE and CSE are used to conduct the state estimation using the updated z (with the disturbance), respectively. Similarly, another 100 samples with the distribution of N(0.01, 0.05 2 ) are applied. The voltage distributions in droop-based and secondary-based EQSE and CSE under the two distributions are shown in FIGS. 4A-B and 5A-B, respectively, where the outlier, maximum, upper quartile, median, lower quartile, and minimum voltages are presented. Specifically, FIGS. 4A and 4B are graphs 400 and 450 depicting exemplary voltage distributions of droop-based EQSE against disturbances of M0.01 , 0.01 2 ) and M0.01 , 0.05 2 ), respectively, according to illustrative embodiments of the invention. Likewise, FIGS. 5 A and 5B are graphs 500 and 550 depicting exemplary voltage distributions of secondary-based EQSE against disturbances of M0.01 , 0.01 2 ) and M0.01 , 0.05 2 ), respectively, according to embodiments of the invention. Table VI below presents the average voltages and variances (from all the 100 samples) of droop-based and secondary-based EQSE and CSE.

Table VI: Average EQSE and CSE results with disturbances (p.u.) [0067] The results in Table VI demonstrate the superior performance achieved using EQSE according to embodiments of the invention. Specifically, from Table VI the following observations can be perceived:

(i) The droop-based EQSE outperforms CSE in terms of the robustness against disturbances. For instance, in FIG. 4A, the voltage distributions of EQSE under the disturbance of M0.01 , 0.01 2 ) are much narrower than those of CSE. When the disturbance has a distribution of 2V(0.01, 0.05 2 ) (see FIG. 4B), the voltage variations of CSE become larger, while those of EQSE remain small. This is primarily because CSE relies on a swing bus to balance the disturbances, while EQSE according to embodiments of the invention redistributes the disturbances into multiple DERs. FIGS. 5A-B show that the secondary-based EQSE is also robust against disturbances.

(ii) Table VI further validates the robustness of both droop-based and secondary-based EQSE under disturbances. For instance, when the disturbance has a distribution of M0.01 . 0.01 2 ), for CSE under secondary control, Vi becomes 0.9986 p.u., which largely deviates from the value with no disturbance, i.e., 1.0000 p.u. under N(0, 0 2 ). By comparison, for EQSE, Vi remains 0.9998 p.u., i.e., a value much closer to 1.0000 p.u. This becomes more apparent when the disturbance has a distribution of M0.01 , 0.05 2 ).

[0068] In terms of convergence, EQSE according to embodiments of the invention also provides superior performance results. By way of example only and without limitation, for both EQSE and CSE, under each disturbance (from 100 samples with M0.01 , 0.01 2 )), the required number of iterations for producing the final result is recorded, in droop-based and secondary-based systems, respectively. FIG. 6 is a bar graph 600 visually depicting exemplary convergence performance results of EQSE according to embodiments of the invention compared with CSE. With reference to FIG. 6, the following insights are obtained:

(i) For either droop-based or secondary-based EQSE, the required number of iterations has a similar level with that of CSE. For instance, all the recorded numbers of iterations for the droop-based EQSE are 4 (see EQSE 1 in FIG. 6), while those for CSE are either 3 or 4 (see CSE 1 in FIG. 6).

(ii) The convergence performance of EQSE is less likely to be affected by disturbances than that of CSE. For instance, in FIG. 6, all the recorded numbers of iterations for the droop-based EQSE are 4 (see EQSE 1 is FIG. 6) and those for the secondary-based EQSE are 5 (see EQSE 2 in FIG. 6). However, those for CSE are distributed (see CSE 1 and CSE 2 in FIG. 6).

[0069] FIG. 7 is a block diagram depicting at least a portion of an exemplary computing system 700 configurable for implementing at least a portion of one or more methods of the invention described herein. The computing system 700 is only one example of a suitable computer system and is not intended to limit the scope of use or functionality of embodiments of the invention described herein. In the computing system 700, there are components, which are operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well-known computing systems, environments, and/or configurations that may be suitable for use with the computing system 700 include, but are not limited to, quantum computing systems, personal computer systems, server computer systems, thin clients, thick clients, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, distributed computing systems, programmable consumer electronics, network PCs, minicomputer systems, mainframe computer systems, and distributed cloud computing environments that include any of the above systems or devices, and the like.

[0070] The computing system 700 may be described in the general context of computer systemexecutable instructions, such as program modules, being executed by one or more processors. Generally, program modules may include routines, programs, objects, components, logic, data structures, and so on that perform particular tasks or implement particular abstract data types. The computing system 700 may be practiced in a distributed cloud computing environment where tasks are performed by remote processing devices that are linked through a communications network. In a distributed cloud computing environment, program modules may be located in both local and remote computer system storage media, including memory storage devices.

[0071] As shown in FIG. 7, the computing system 700 is shown in the form of a general-purpose computing device. Components of the computing system 700 may include, but are not limited to, one or more processors or processing units 702, a system memory 704, and a bus 706 that operatively couples various system components including system memory 704 to the processor 702. Bus 706 represents one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. By way of example only and without limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnects (PCI) bus. The computing system 700 typically includes a variety of computer system readable media. Such media may be any available media that is accessible by the computing system 700, and it includes both, volatile and non-volatile media, removable and non-removable media.

[0072] The system memory 704 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 708 and/or cache memory 710. The computing system 700 may further include other removable/non-removable, volatile/non-volatile computer system storage media. By way of example only, a storage system 712 may be provided for reading from and writing to a non-removable, non-volatile magnetic media (not shown and typically called a “hard drive”). Although not explicitly shown, a magnetic disk drive for reading from and writing to a removable, non-volatile magnetic disk (e.g., a “floppy disk”), an optical disk drive for reading from or writing to a removable, non-volatile optical disk such as a CD-ROM, DVD-ROM or other optical media, and/or a solid-state drive (SSD) may be provided. In such instances, each can be connected to the bus 706 by one or more data media interfaces. As will be further depicted and described below, memory 704 may include at least one program product having a set (e.g., at least one) of program modules 714 that are configured to carry out one or more functions according to embodiments of the invention.

[0073] The program/utility, having a set (at least one) of program modules 714, may be stored in memory 704 by way of example, and not limiting, as well as an operating system, one or more application programs, other program modules, and program data. Each of the operating systems, one or more application programs, other program modules, and program data or some combination thereof, may include an implementation of a networking environment. Program modules 714 generally carry out the functions and/or methodologies of embodiments of the invention, as described herein.

[0074] The computing system 700 may also communicate with one or more external devices 716 such as a keyboard, a pointing device, a display 718, etc.; one or more devices that enable a user to interact with the computing system 700; and/or any devices (e.g., network card, modem, etc.) that enable the computing system 700 to communicate with one or more other computing devices. Such communications can occur via input/output (I/O) interfaces 720. Still yet, the computing system 700 may communicate with one or more networks, such as a local area network (LAN), a general wide area network (WAN), and/or a public network (e.g., the Internet), via a network adapter 722. As depicted, the network adapter 722 may communicate with the other components of the computing system 700 via bus 706. It should be understood that although not explicitly shown, other hardware and/or software components could be used in conjunction with the computing system 700. Examples, include, but are not limited to, microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data archival storage systems, etc.

[0075] The processing unit 702 in the computing system 700 may communicate with one or more other processors, such as a general processing unit (GPU) 724, via the bus 706, such as in a distributed computing arrangement. The GPU 724 may include its own internal memory for at least temporarily storing intermediate processing results, program instructions, etc. The GPU 724 may also be configured to communicate with the one or more external devices 716 via the I/O interface 720 or with the one or more networks via the network adapter 722.

[0076] The descriptions of the various embodiments of the present invention have been presented for purposes of illustration, but are not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skills in the art given the teachings herein without departing from the scope and spirit of the described embodiments. The terminology used herein was chosen to best explain the principles of the embodiments, the practical application or technical improvement over conventional technologies, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.

[0077] The present invention may be embodied as a system, a method, and/or a computer program product. The computer program product may include a non-transitory computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.

[0078] The non-transitory storage medium may be an electronic, magnetic, optical, electromagnetic, infrared or a semi-conductor system for a propagation medium. Examples of a computer-readable medium may include a semiconductor or solid-state memory, magnetic tape, a removable computer diskette, a random access memory (RAM), a read-only memory (ROM), a rigid magnetic disk, a flash drive, and an optical disk. Current examples of optical disks include compact disk-read only memory (CD-ROM), compact disk-read/write (CD-R/W), digital versatile disk (DVD), Blu-Ray Disk, etc.

[0079] The computer readable storage medium can be any tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.

[0080] Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.

[0081] Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object-oriented programming language such as Smalltalk, C++ or the like, and conventional procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user’s computer, partly on the user's computer as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user’s computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.

[0082] Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood by those skilled in the art that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.

[0083] These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.

[0084] The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or another devices to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or another device implement the functions/acts specified in the flowchart and/or block diagram block or blocks. [0085] The flowcharts and/or block diagrams in the Figures illustrate the architecture, functionality, and/or operation of possible implementations of systems, methods, and computer program products according to one or more embodiments of the present invention. In this regard, each block in the flowchart or block diagrams, or each program code line in an illustrated algorithm, may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified function(s). In some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be appreciated that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or act or carry out combinations of special purpose hardware and computer instructions.

[0086] At least a portion of the techniques of the present invention may be implemented in an integrated circuit. In forming integrated circuits, identical die are typically fabricated in a repeated pattern on a surface of a semiconductor wafer. Each die includes a device described herein, and may include other elements, structures and/or circuits. The individual die are cut or diced from the wafer, then packaged as an integrated circuit. One skilled in the art would know how to dice wafers and package die to produce integrated circuits. Any of the exemplary systems and circuits illustrated in the accompanying figures, or portions thereof, may be part of an integrated circuit. Integrated circuits so manufactured are considered part of this invention.

[0087] Those skilled in the art will appreciate that the exemplary systems and circuits discussed above, or portions thereof, can be distributed in raw form (i.e., a single wafer having multiple unpackaged chips), as bare dies, in packaged form, or incorporated as parts of intermediate products or end products that benefit from having a microgrid controller formed therein in accordance with one or more embodiments of the invention, such as, for example, an SDC microgrid control system.

[0088] An integrated circuit in accordance with aspects of the present disclosure can be employed in essentially any microgrid control application and/or electronic system. Systems incorporating such integrated circuits are considered part of this invention. Given the teachings of the present disclosure provided herein, one of ordinary skill in the art will be able to contemplate other implementations and applications of embodiments of the invention.

[0089] The illustrations of embodiments of the invention described herein are intended to provide a general understanding of the various embodiments and are not intended to serve as a complete description of all the elements and features of apparatus and systems that might make use of the module, circuits and/or methods described herein. Many other embodiments will become apparent to those skilled in the art given the teachings herein; other embodiments are utilized and derived therefrom, such that structural and logical substitutions and changes can be made without departing from the scope of this disclosure. The drawings are also merely representational and are not drawn to scale. Accordingly, the specification and drawings are to be regarded in an illustrative rather than a restrictive sense.

[0090] Embodiments of the invention are referred to herein, individually and/or collectively, by the term “embodiment” merely for convenience and without intending to limit the scope of this application to any single embodiment or inventive concept if more than one is, in fact, shown. Thus, although specific embodiments have been illustrated and described herein, it should be understood that an arrangement achieving the same purpose can be substituted for the specific embodiment(s) shown; that is, this disclosure is intended to cover any and all adaptations or variations of various embodiments. Combinations of the above embodiments, and other embodiments not specifically described herein, will become apparent to those of skill in the art given the teachings herein.

[0091] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms “a,” “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms “comprises” and/or “comprising,” when used in this specification, specify the presence of stated features, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, elements, components, and/or groups thereof. Locational terms such as “above,” “below,” “in,” “out,” “internal,” and “external,” as may be used herein, are intended to describe the relative position of elements or structures to each other as opposed to an absolute position of the elements or structures. [0092] The corresponding structures, materials, acts, and equivalents of all means or step-plus- function elements in the claims below are intended to include any structure, material, or act for performing the function in combination with other claimed elements as specifically claimed. The description of the various embodiments has been presented for purposes of illustration and description, but is not intended to be exhaustive or limited to the forms disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the invention. The embodiments were chosen and described in order to best explain the principles of the invention and the practical application, and to enable others of ordinary skill in the art to understand the various embodiments with various modifications as are suited to the particular use contemplated.

[0093] The abstract is provided to comply with 37 C.F.R. § 1.72(b), which requires an abstract that will allow the reader to quickly ascertain the nature of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. In addition, in the foregoing Detailed Description, it can be seen that various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the appended claims reflect, inventive subject matter lies in less than all features of a single embodiment. Thus, the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as separately claimed subject matter.

[0094] Given the teachings of embodiments of the invention provided herein, one of ordinary skill in the art will be able to contemplate other implementations and applications of the techniques of embodiments of the invention. Although illustrative embodiments of the invention have been described herein with reference to the accompanying drawings, it is to be understood that embodiments of the invention are not limited to those precise embodiments, and that various other changes and modifications are made therein by one skilled in the art without departing from the scope of the appended claims.