CINXE.COM

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head><script type="text/javascript" src="/_static/js/bundle-playback.js?v=HxkREWBo" charset="utf-8"></script> <script type="text/javascript" src="/_static/js/wombat.js?v=txqj7nKC" charset="utf-8"></script> <script>window.RufflePlayer=window.RufflePlayer||{};window.RufflePlayer.config={"autoplay":"on","unmuteOverlay":"hidden"};</script> <script type="text/javascript" src="/_static/js/ruffle/ruffle.js"></script> <script type="text/javascript"> __wm.init("https://web.archive.org/web"); __wm.wombat("http://www.patentgenius.com/patent/8626353.html","20190415173352","https://web.archive.org/","web","/_static/", "1555349632"); </script> <link rel="stylesheet" type="text/css" href="/_static/css/banner-styles.css?v=S1zqJCYt" /> <link rel="stylesheet" type="text/css" href="/_static/css/iconochive.css?v=3PDvdIFv" />  <title>Integration of demand response and renewable resources for power generation management - Patent # 8626353 - PatentGenius</title> <meta name="description" content="System and method of solving, in a single-period, an optimal dispatching problem for a network of energy generators connected via multiple transmission lines, where it is sought to find the lowest operational cost of dispatching of various energy sources to satisfy demand. The model includes traditional thermal resources and renewable energy resources available generation capabilities within the grid. The method considers demand reduction as a virtual generation source that can be dispatched quickly to hedge against the risk of unforeseen shortfall in supply. Demand reduction is dispatched in response to incentive signals sent to consumers. The control options of the optimization model consist of the dispatching order and dispatching amount energy units at generators together with the rebate signals sent to end-users at each node of the network under a demand response policy. Numerical experiments based on an analysis of representative data illustrate the effectiveness of demand response as a hedging option."/><meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <link href="/web/20190415173352cs_/http://www.patentgenius.com/patentgenius.css" rel="stylesheet" type="text/css"> </head> <body leftmargin="0" topmargin="0" marginwidth="0" marginheight="0"> <table width="100%" height="100%" border="0" cellpadding="0" cellspacing="0"> <tr> <td align="left" valign="top"> <table width="100%" border="0" cellspacing="0" cellpadding="0"> <tr align="left" valign="top"> <td width="23%" background="/web/20190415173352im_/http://www.patentgenius.com/images/bg_topraw.jpg"><a href="/web/20190415173352/http://www.patentgenius.com/"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/logo.jpg" border="0" width="418" height="156"></a></td> <td width="77%" height="156" align="right" background="/web/20190415173352im_/http://www.patentgenius.com/images/bg_topraw.jpg"> <table width="385" border="0" cellspacing="0" cellpadding="0"> <tr> <td height="37" align="right" valign="middle"> </td> </tr> <tr> <td height="41" align="right" valign="top"><form action="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/searchresults.php" id="cse-search-box"> <input type="hidden" name="cx" value="partner-pub-1302519442017013:b6x8u1-id06"/> <input type="hidden" name="cof" value="FORID:10"/> <input type="hidden" name="ie" value="ISO-8859-1"/> <table width="100%" border="0" cellspacing="0" cellpadding="0"> <tr align="left" valign="middle"> <td width="22%" background="https://web.archive.org/web/20190415173352im_/http://www.patentgenius.com/images/bg_searchbox.jpg"><img src="https://web.archive.org/web/20190415173352im_/http://www.patentgenius.com/images/search.jpg" width="83" height="41"></td> <td width="64%" background="https://web.archive.org/web/20190415173352im_/http://www.patentgenius.com/images/bg_searchbox.jpg"><input name="q" type="text" size="40"></td> <td width="14%" background="https://web.archive.org/web/20190415173352im_/http://www.patentgenius.com/images/bg_searchbox.jpg"><input type="submit" name="sa" value="Search" src="https://web.archive.org/web/20190415173352im_/http://www.patentgenius.com/images/butt_go.jpg"></td> <script type="text/javascript" src="https://web.archive.org/web/20190415173352js_/http://www.google.com/cse/brand?form=cse-search-box&lang=en"></script> </tr> </table> </form></td> </tr> <tr> <td height="28" align="right" valign="top"> </td> </tr> <tr> <td height="31" align="center" valign="top" class="link1"><table width="80%" border="0" cellspacing="0" cellpadding="0"> <tr align="center" valign="middle" class="link1"> <td width="32%"><a href="/web/20190415173352/http://www.patentgenius.com/resources.html" class="link1">Resources</a></td> <td width="4%"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/linksep.gif" width="1" height="37"></td> <td width="30%"><a href="/web/20190415173352/http://www.patentgenius.com/contact.html" class="link1">Contact Us</a></td> <td width="2%"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/linksep.gif" width="1" height="37"></td> <td width="32%"><a href="/web/20190415173352/http://www.patentgenius.com/" class="link1">Home</a></td> </tr> </table></td> </tr> </table> </td> </tr> </table> <table width="100%" border="0" cellspacing="0" cellpadding="0"> <tr> <td align="left" valign="top" background="/web/20190415173352im_/http://www.patentgenius.com/images/bg_bottomlinks.gif"> <table width="600" border="0" cellspacing="0" cellpadding="0"> <tr class="link1"> <td width="96" align="right" valign="middle"><font color="#FFFFFF">Browse by:</font></td> <td width="129" align="center" valign="middle"><a href="/web/20190415173352/http://www.patentgenius.com/inventor.html" class="link1">INVENTOR</a></td> <td width="12"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/linksep2.gif" width="1" height="22"></td> <td width="113" align="center" valign="middle"><a href="/web/20190415173352/http://www.patentgenius.com/assignee.html" class="link1">PATENT HOLDER</a></td> <td width="10" align="center" valign="middle"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/linksep2.gif" width="1" height="22"></td> <td width="143" align="center" valign="middle"><a href="/web/20190415173352/http://www.patentgenius.com/number.html" class="link1">PATENT NUMBER</a></td> <td width="10" align="center" valign="middle"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/linksep2.gif" width="1" height="22"></td> <td width="77" height="35" align="center" valign="middle"><a href="/web/20190415173352/http://www.patentgenius.com/date.html" class="link1">DATE</a></td> </tr> </table></td> </tr> </table> <table width="100%" border="0" cellspacing="1" cellpadding="0"> <tr> <td width="84%" align="center" valign="top"><table width="100%" border="0" cellspacing="0" cellpadding="0"> <tr> <td width="1%" align="right" valign="bottom"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/box_leftcorner.gif" width="29" height="22"></td> <td width="98%" background="/web/20190415173352im_/http://www.patentgenius.com/images/box_topbg.gif"> </td> <td width="1%" align="left" valign="bottom"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/box_rightcorner.gif" width="29" height="22"></td> </tr> <tr> <td align="right" valign="top" background="/web/20190415173352im_/http://www.patentgenius.com/images/box_leftvrtical.gif"> </td> <td align="center" valign="top"> <table width="100%" border="0" cellspacing="0" cellpadding="0"> <tr bgcolor="#D6D6D6"> <td width="2%" align="left" valign="top" bgcolor="#D6D6D6"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/edge.gif" width="21" height="17"></td> <td width="98%" height="30" align="left" valign="middle" bgcolor="#D6D6D6" class="heading2"><strong> Integration of demand response and renewable resources for power generation management</div></strong></td> </tr> </table> <table width="100%" border="0" cellspacing="2" cellpadding="3"> <tr align="left" valign="top"> <td width="50%" align="left" class="dots"> <script type="text/javascript" src="/web/20190415173352js_/http://www.patentgenius.com/js/wrapper.js"></script> <table width="100%" border="0" cellspacing="2" cellpadding="2" style="font-size: 11px;"> <tr style="font-size: 12px;"> <td width="150"><strong>8626353</strong></td> <td><strong>Integration of demand response and renewable resources for power generation management</strong> </td> </tr> <tr> <td colspan="2"></td> </tr></td></tr><tr><td valign="top"><strong>Patent Drawings:</strong></td><td align="center" id="images"></td></tr><tr><td colspan="2"><br><center> <script type="text/javascript"> </script> <script type="text/javascript" src="https://web.archive.org/web/20190415173352js_/http://pagead2.googlesyndication.com/pagead/show_ads.js"> </script> </center> </td></tr><tr bgcolor="#F0F0F0"> <td><strong>Inventor:</strong></td> <td>Ghosh, et al.</td> </tr> <tr bgcolor="#FFFFFF"> <td><strong>Date Issued:</strong></td> <td>January 7, 2014</td> </tr> <tr bgcolor="#F0F0F0"> <td><strong>Application:</strong></td> <td></td> </tr> <tr bgcolor="#FFFFFF"> <td><strong>Filed:</strong></td> <td></td> </tr> <tr bgcolor="#F0F0F0"> <td valign="top"><strong>Inventors:</strong></td> <td></td> </tr> <tr bgcolor="#FFFFFF"> <td><strong>Assignee:</strong></td> <td></td></tr> <tr bgcolor="#F0F0F0"> <td><strong>Primary Examiner:</strong></td> <td>Padmanabhan; Kavita</td> </tr> <tr bgcolor="#FFFFFF"> <td><strong>Assistant Examiner:</strong></td> <td>Dunn; Darrin</td> </tr> <tr bgcolor="#F0F0F0"> <td><strong>Attorney Or Agent:</strong></td> <td>Scully, Scott, Murphy & Presser, P.C.Morris, Esq.; Daniel P.</td> </tr> <tr bgcolor="#FFFFFF"> <td valign="top"><strong>U.S. Class:</strong></td> <td><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/700/291.html">700/291</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/700/28.html">700/28</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/700/286.html">700/286</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/700/297.html">700/297</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/700/32.html">700/32</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/703/1.html">703/1</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/703/3.html">703/3</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/705/37.html">705/37</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/705/7.31.html">705/7.31</a></td> </tr> <tr bgcolor="#F0F0F0"> <td valign="top"><strong>Field Of Search:</strong></td> <td><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/28.html">;700/28</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/32.html">;700/32</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/33.html">;700/33</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/36.html">;700/36</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/286.html">;700/286</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/291.html">;700/291</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;700/297.html">;700/297</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;703/1.html">;703/1</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;703/3.html">;703/3</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;705/7.31.html">;705/7.31</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;705/10.html">;705/10</a>; <a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/class/;705/37.html">;705/37</a></td> </tr> <tr bgcolor="#FFFFFF"> <td valign="top"><strong>International Class:</strong></td> <td>G05D 3/12; G05D 13/02; G06F 17/50; G06Q 40/00; G06Q 10/00; G06F 7/48</td> </tr> <tr bgcolor="#F0F0F0"> <td valign="top"><strong>U.S Patent Documents:</strong></td> <td></td> </tr> <tr bgcolor="#FFFFFF"> <td valign="top"><strong>Foreign Patent Documents:</strong></td> <td></td> </tr> <tr bgcolor="#F0F0F0"> <td valign="top"><strong>Other References:</strong></td> <td>Bhattacharya et al., "Real Time Optimal Interruptible Tariff Mechanism Incorporating Utility-Customer Interactions," IEEE Transactions onPower Systems, vol. 15, No. 2, May 2000, pp. 700-706. cited by applicant.<br>Office Action dated Jun. 28, 2013 received in a related U.S. Appl. No. 12/890,119. cited by applicant.</td> </tr> <tr bgcolor="#FFFFFF"> <td colspan="2"><br><br><br><br><center> <script type="text/javascript"> </script> <script type="text/javascript" src="https://web.archive.org/web/20190415173352js_/http://pagead2.googlesyndication.com/pagead/show_ads.js"> </script> </center> <br><br><br><br><br></td> </tr> <tr bgcolor="#FFFFFF"> <td valign="top"><strong>Abstract:</strong></td> <td>System and method of solving, in a single-period, an optimal dispatching problem for a network of energy generators connected via multiple transmission lines, where it is sought to find the lowest operational cost of dispatching of various energy sources to satisfy demand. The model includes traditional thermal resources and renewable energy resources available generation capabilities within the grid. The method considers demand reduction as a virtual generation source that can be dispatched quickly to hedge against the risk of unforeseen shortfall in supply. Demand reduction is dispatched in response to incentive signals sent to consumers. The control options of the optimization model consist of the dispatching order and dispatching amount energy units at generators together with the rebate signals sent to end-users at each node of the network under a demand response policy. Numerical experiments based on an analysis of representative data illustrate the effectiveness of demand response as a hedging option.</td> </tr> <tr bgcolor="#F0F0F0"> <td valign="top"><strong>Claim:</strong></td> <td>Having thus described our invention, what we claim as new, and desire to secure by Letters Patent is:<br><br>1. A computer-implemented method for managing dispatch of energy for distribution tocustomers of an energy utility, said method comprising: receiving, at a processor device, configuration data of an energy distribution network including at least one energy generator for generating energy units, and generator busses carrying energy unitsdispatched from said generators to customers via transmission lines, receiving, at a processor device, data representing a renewable energy source injected for power generation on one of said generator busses; specifying decision variables forcontrolling a dispatching of energy units from energy sources over said network and a total power extracted at connected generator busses over said network in a defined time period; specifying a set of optimization constraints including a risk controlconstraint as a function of said decision variables, said configuration data and injected energy data; formulating, via mixed integer programming at said processor device, an optimization problem having an a cost-based objective function representing anobjective to minimize a total cost metric within an acceptable risk level; solving, by said processing device, said optimization problem subject to said risk control constraint, set of optimization constraints, the objective and the cost-based objectivefunction to jointly determine an optimal dispatching policy for generated energy units and an incentive rate for customers of the utility to extract load reduction from said customers total power at connected generator busses over said network in saiddefined time period that minimizes said total cost metric within said acceptable risk level, said dispatching policy comprising an on/off status indicator of a generator i for each of one or more generator buses i, a real power generation P.sup.g.sup.iof each generator bus i, a voltage phase angle .theta..sub.i of each bus, and an energy unit incentive rebate price r offered to customers on the one or more buses; and, outputting said dispatching policy and incentive rate for adoption by the utilityto control energy dispatch and customer demand response.<br><br>2. The computer-implemented method as claimed in claim 1, wherein said cost-based objective function represents a sum of a utility's generation costs minus a total revenue from a load position exposed to an energy spot market, said methodcomprising configuring said cost-based loss function as a nonlinear mixed integer program (MIP).<br><br>3. The computer-implemented method as claimed in claim 2, wherein said solving includes implementing a gradient descent technique to solve said cost-based objective function.<br><br>4. The computer-implemented method as claimed in claim 1, wherein said output solution indicates: a dispatching order and the dispatching amount of said energy units, and in addition, said utility offering a rebate to customers according tosaid rebate price for reduced energy load through communicated price signals sent to customer smart meters via a generator bus for said time period.<br><br>5. The computer-implemented method as claimed in claim 1, wherein said set of optimization constraints include one or more bounding constraints on the decision variables that specify minimum and maximum limits of the a real power generationP.sup.g.sup.i of each generator bus, a power (load) demand P.sup.d.sup.i at each bus; and voltage phase angles .theta..sub.i limits of all busses; and, a constraint to ensure balancing of a total generation load plus the reduced load against a totaldemand.<br><br>6. The computer-implemented method as claimed in claim 1, wherein said set of optimization constraints include one or more power flow balance constraints on the decision variables to achieve a system level power balance when aggregating: a realload power generation output at said bus P.sup.g.sup.i including injection of a renewable energy source P.sub.RN at a generator bus, a power demand response P.sub.DR as a source of virtual generation at said bus, and a real load output of P.sup.g.sup.iadjustment by an amount of load bought from (sold to) said spot market at a generator bus having access to said spot market, said cost-based objective function modeled as function of decision variables z, P.sup.g and r, and random variables P.sup.d,P.sub.DR(r) and P.sub.RN.<br><br>7. The computer-implemented method as claimed in claim 4, further comprising: determining a risk control value within a level that is traded off against offered rebates for load reduction for said time period.<br><br>8. The computer-implemented method as claimed in claim 1, wherein a risk measure imposed by said risk constraint on the cost-based function for said optimization problem is a Conditional Value-at-Risk (CVaR), said method includes receiving aspecified risk level for said risk control constraint.<br><br>9. A system for managing dispatch of energy for distribution to customers of an energy utility, said method comprising: a memory storage device; and a processor device connected to the memory storage device, said processor device configured toperform a method comprising: receiving configuration data of an energy distribution network including at least one energy generator for generating energy units, and generator busses carrying energy units dispatched from said generators to customers viatransmission lines, receiving data representing a renewable energy source injected for power generation on one of said generator busses; specifying decision variables for controlling a dispatching of energy units from energy sources over said networkand a total power extracted at connected generator busses over said network in a defined time period; specifying a set of optimization constraints including a risk control constraint as a function of said decision variables, said configuration data andinjected energy data; formulating, via mixed integer programming, an optimization problem having an a cost-based objective function representing an objective to minimize a total cost metric within an acceptable risk level; solving said optimizationproblem subject to said risk control constraint, set of optimization constraints, the objective and the cost-based objective function to jointly determine an optimal dispatching policy for generated energy units and an incentive rate for customers of theutility to extract load reduction from said customers at connected generator busses over said network in said defined time period that minimizes said total cost metric within said acceptable risk level, said dispatching policy comprising an on/off statusindicator of a generator i for each of one or more generator buses i, a real power generation P.sup.g.sup.i of each generator bus i, a voltage phase angle .theta..sub.i of each bus, and an energy unit incentive rebate price r offered to customers on theone or more buses; and, outputting said dispatching policy and incentive rate for adoption by the utility to control energy dispatch and customer demand response.<br><br>10. The system as claimed in claim 9, wherein said cost-based objective function represents a sum of a utility's generation costs minus a total revenue from a load position exposed to an energy spot market, said method comprising configuringsaid cost-based loss function as a nonlinear mixed integer program (MIP).<br><br>11. The system as claimed in claim 10, wherein said solving includes implementing a gradient descent technique to solve said cost-based objective function.<br><br>12. The system as claimed in claim 9, wherein said output solution indicates: a dispatching order and the dispatching amount of said energy units, and in addition, said utility offering a rebate to customers according to said rebate price forreduced energy load through communicated price signals sent to customer smart meters via a generator bus for said time period.<br><br>13. The system as claimed in claim 9, wherein said set of optimization constraints include one or more bounding constraints on the decision variables that specify minimum and maximum limits of the a real power generation P.sup.g.sup.i of eachgenerator bus, a power (load) demand P.sup.d.sup.i at each bus; and voltage phase angles .theta..sub.i limits of all busses; and, a constraint to ensure balancing of a total generation load plus the reduced load against a total demand.<br><br>14. The system as claimed in claim 9, wherein said set of optimization constraints include one or more power flow balance constraints on the decision variables to achieve a system level power balance when aggregating: a real load powergeneration output at said bus P.sup.g.sup.i including injection of a renewable energy source P.sub.RN at a generator bus, and a power demand response P.sub.DR as a source of virtual generation at said bus, and a real load output of P.sup.g.sup.iadjustment by an amount of load bought from (sold to) said spot market at a generator bus having access to said spot market, said cost-based objective function modeled as function of decision variables z, P.sub.g and r, and random variables P.sup.dP.sub.DR (r) and P.sub.RN.<br><br>15. The system as claimed in claim 12, further comprising: determining a risk control value within a level that is traded off against offered rebates for load reduction.<br><br>16. The system as claimed in claim 9, wherein a risk measure imposed by said risk constraint on the cost-based function for said optimization problem is a Conditional Value-at-Risk (CVaR), said method includes receiving a specified risk levelfor said risk control constraint.<br><br>17. A computer program device for managing dispatch of energy for distribution to customers of an energy utility, the computer program device comprising a storage medium, said medium not a propagating signal, said medium readable by aprocessing circuit and storing instructions run by the processing circuit for performing a method, the method comprising: receiving configuration data of an energy distribution network including at least one energy generator for generating energy units,and generator busses carrying energy units dispatched from said generators to customers via transmission lines, receiving data representing a renewable energy source injected for power generation on one of said generator busses; specifying decisionvariables for controlling a dispatching of energy units from energy sources over said network and a total power extracted at connected generator busses over said network in a defined time period; specifying a set of optimization constraints including arisk control constraint as a function of said decision variables, said configuration data and injected energy data; formulating, via mixed integer programming, an optimization problem having an a cost-based objective function representing an objectiveto minimize a total cost metric within an acceptable risk level; solving said optimization problem subject to said risk control constraint, set of optimization constraints, the objective and the cost-based objective function to jointly determine anoptimal dispatching policy for generated energy units and an incentive rate for customers of a utility to extract load reduction from said customers at connected generator busses over said network in said defined time period that minimizes said totalcost metric within said acceptable risk level, said dispatching policy comprising an on/off status indicator of a generator i for each of one or more generator buses i, a real power generation P.sup.g.sup.i of each generator bus i, a voltage phase angle.theta..sub.i of each bus, and an energy unit incentive rebate price r offered to customers on the one or more buses; and, outputting said dispatching policy and incentive rate for adoption by an utility to control energy dispatch and customer demandresponse.<br><br>18. The computer program device as claimed in claim 17, wherein said cost-based objective function represents a sum of a utility's generation costs minus a total revenue from a load position exposed to an energy spot market, said methodcomprising configuring said cost-based loss function as a nonlinear mixed integer program (MIP).<br><br>19. The computer program device as claimed in claim 18, wherein said solving includes implementing a gradient descent technique to solve said cost-based objective function.<br><br>20. The computer program device as claimed in claim 19, wherein said output solution indicates: a dispatching order and the dispatching amount of said energy units, and in addition, said utility offering a rebate to customers according to saidrebate price for reduced energy load through communicated price signals sent to customer smart meters via a generator bus for said time period.<br><br>21. The computer program device as claimed in claim 20, further comprising: determining a risk control value within a level that is traded off against offered rebates for load reduction for said time period.</td> </tr> <tr bgcolor="#FFFFFF"> <td valign="top"><strong>Description:</strong></td> <td>CROSS REFERENCE TO RELATED APPLICATIONS<br><br>This application relates to commonly-owned, co-pending U.S. patent application Ser. No. 12/890,119 entitled SYSTEM AND METHOD FOR LOWEST COST AGGREGATE ENERGY DEMAND REDUCTION, the whole contents and disclosure of which is incorporated byreference as if fully set forth herein.<br><br>The present invention relates generally to smart grid technologies as it pertains to managing the generation, distribution and usage of energy, and, more particularly, to system and method for controlling the dispatching of energy (e.g.,committed regular generation units) over an energy distribution network to customers that integrates demand response and availability of renewable resources for power generation management.<br><br>BACKGROUND<br><br>With the advent of the Smart Grid that provides an infrastructure for energy supply (e.g., electricity) generation and transmission for delivery from suppliers to consumers under control of digital communications, the ability to access real-timeinformation on supply availability and prices supported by demand offers unique opportunities to improve the overall efficiency of the energy grid in terms of both long-term supply-demand management and near-term dispatching of diverse generationfacilities to meet current demand. The responsiveness and flexibility envisioned for the Smart Grid provide additional advantages in facing the significant new challenges of integrating distributed and intermittent generation capability, such assmall-scale generators and renewable energy sources (wind, solar, etc.), at a scale that current grid technology simply cannot achieve. This is becoming more critical as renewable energy technologies are playing an increasingly important role in theportfolio mix of electricity generation.<br><br>Currently, only a few published methods exist that considers dispatching policy with distributed renewable generation (e.g., wind energy). However, the industry has yet to produce a system or method that performs integrating renewablegeneration and demand response together into power generation management and risk control.<br><br>That is, existing utility generation dispatching policies do not integrate demand response for risk management.<br><br>It would be highly desirable to provide a system and method that achieves improvements in both total financial measures and risk control over existing mechanisms that do not consider virtual generation (demand response) and distributedintermittent generation (renewable sources).<br><br>SUMMARY<br><br>System and method for determining optimal dispatching policy of aggregate energy generation capacity while integrating consideration of renewable resources and demand response program at one or more electricity network levels such as an energygeneration and distribution network.<br><br>The system and method when implemented by an electricity generation utility helps optimize their total financial performance by providing options for hedging against energy shortfalls with intelligent responsive dispatching of various energysources including renewables (e.g., wind, solar, thermal, etc.) and a "virtual" generation source of incentivized demand reduction from end-users (a "Demand response").<br><br>In one aspect, the system and method performs solving, in a time period(s), an optimal dispatching problem for a network of one or more energy utilities connected via multiple transmission lines, and in connection with a configuredmicroprocessor(s) of a computing system, generate solutions enabling choices for utilities to achieve find a lowest operational cost dispatching of various energy sources to satisfy demand. The model includes traditional thermal resources and renewableenergy resources available generation capabilities within the grid.<br><br>In a further aspect, the method considers demand reduction as a virtual generation source that can be dispatched quickly to hedge against the risk of unforeseen shortfall in supply. Demand reduction is dispatched in response to incentive rebatesignals sent to consumers.<br><br>In a further aspect, an optimization model is generated and solutions provided that includes control options including: the dispatching order and dispatching amount of the thermal generators together with the rebate signals sent to end-users(e.g., incentives to residential customers to increase the energy efficiency of their homes) at each node of the network under a simple demand response policy.<br><br>Accordingly, there is provided a system, method and computer program product for managing dispatch of energy for distribution to customers of an energy utility, the method comprising:<br><br>receiving, at a processor device, configuration data of an energy distribution network including at least one energy generator for generating energy units, and generator busses carrying energy units dispatched from the generators to customersvia transmission lines,<br><br>receiving, at a processor device, data representing a renewable energy source injected for power generation on one of the generator busses;<br><br>specifying decision variables for controlling a dispatching of energy units from energy sources over the network and a total power extracted at connected generator busses over the network in a defined time period;<br><br>specifying a set of optimization constraints including a risk control constraint as a function of the decision variables, the configuration data and injected energy data;<br><br>formulating, via mixed integer programming at the processor device, an optimization problem having an a cost-based objective function representing an objective to minimize a total cost metric within an acceptable risk level;<br><br>solving, by the processing device, the optimization problem subject to the risk control constraint and set of optimization constraints, the objective and the cost-based objective function to jointly determine an optimal dispatching policy forgenerated energy units and an incentive rate for customers of a utility to extract load reduction from the customers total power at connected generator busses over the network in the defined time period that minimizes the total cost metric within theacceptable risk level; and,<br><br>outputting the dispatching policy and incentive rate for adoption by an utility to control energy dispatch and customer demand response.<br><br>A computer program product is provided for performing operations. The computer program product includes a storage medium readable by a processing circuit and storing instructions run by the processing circuit for running a method. The methodis the same as listed above.<br><br>BRIEF DESCRIPTION OF THE DRAWINGS<br><br>The objects, features and advantages of the present invention will become apparent to one skilled in the art, in view of the following detailed description taken in combination with the attached drawings, in which:<br><br>FIG. 1 depicts a schematic of all end users on a feeder network 201 where each user has a price elasticity represented by the demand response curves;<br><br>FIG. 2 depicts an example IEEE 30-bus 6-generator system topology for which the present system and method may be implemented in accordance with one embodiment;<br><br>FIG. 3 depicts an example plot comparing Mean Loss vs. Wind with/without DR for example varying spot market price assumptions of a numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 4A shows an example plot of Cost & Rebate Load vs. Risk Level and FIG. 4B shows an example plot of load profiles vs. the tolerated risk levels at an example Market Price=3 and Wind=45 according to one example numerical experiment computedfor example system topology of FIG. 2;<br><br>FIG. 5 shows an example plot of Loss/Rebate vs. Risk Level at various states of Market Price=5 and Wind=45 according to one example numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 6 shows an example plot of Load Profiles vs. Risk Level at various states of Market Price=5 and Wind=45 according to one example numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 7 shows an example plot of Loss vs. Wind's Coefficient of Variation according to one example numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 8 shows an example plot of Rebate vs. Wind's Coefficient of Variation according to one example numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 9 illustrates the optimal cost under changing market price relations with respect to different choices for the busses that interconnect wind, DR or spot market load injections according to one example numerical experiment computed forexample system topology of FIG. 2;<br><br>FIG. 10 is a plot that compares the optimal cost versus risk relations under different wind generation injection bus locations according to one example numerical experiment computed for example system topology of FIG. 2;<br><br>FIG. 11 depicts a computer-implemented method for managing dispatch of energy over a network according to one embodiment;<br><br>FIG. 12 depicts an example solution for said optimization problem implementing a two-stage procedure in one embodiment;<br><br>FIG. 13 depicts conceptually an example high performance computing application 400 running on a multiprocessor computing system that employs the selective pairing facility in one embodiment;<br><br>DETAILED DESCRIPTION<br><br>In view of FIG. 1, employed at the Smart Grid 12 is a computer implemented system and method for determining an optimal dispatching policy with integration of renewable resources and demand response program at multiple network levels such asgeneration and distribution networks including, e.g., distribution and feeder networks 200.<br><br>In an embodiment, the computer implemented system and method employs a model including an algorithm for determining an optimal dispatching policy for electricity generation utilities. This model accounts for: 1) intermittency and volatility ofrenewable generations; 2) for heterogeneous customer flexibility in load reduction; 3) for non-discriminate incentive pricing for customers; and, 4) can account for hedging risks with demand response for generation utilities with renewable resources.<br><br>More particularly, a method includes a computer-based system applying a model that incorporates risk-based control in the short-term economic dispatching of energy generation capacity including distributed renewable energy sources and demandresponse together with physical constraints on transmission of energy through the distribution network. Examples of dispatchable energy sources include traditional thermal, distributed generation, etc.<br><br>The method jointly determines an optimal dispatching policy for generation (thermal, renewable etc.) and an incentive rate for customers of the utility to extract load reduction from them.<br><br>In additional embodiments, the method accounts for one or more of: intermittency and volatility of renewable generations; heterogeneous customer flexibility in load reduction; non-discriminate incentive pricing for customers; hedging risks withdemand response for generation utilities with renewable resources, etc.<br><br>As shown in FIG. 1, an electric utility controls the dispatching of committed (regular and peaker) generation units (thermal sources) over a network of multiple local buses interconnected via transmission lines. For example, FIG. 1 particularlydepicts a schematic of all end users on an example feeder (energy distribution) network 200 where each user has a price elasticity represented by demand curves. As represented in FIG. 2, demand users include classes of customers such as offices 202,residential 204 (e.g., homes) and commercial 206, e.g., retail stores. A network of distribution feeders 201 import the power from grid 12 in response to the aggregated reactive and potential demand, and the utility provides power to the demand users,e.g., organized according to zones 210 corresponding to a respective demand class, for example. Thus, in one aspect, the system and method employing the model, is applied at multiple levels of the energy distribution hierarchies such as shown in FIG. 1where the various levels at which the method can be applied include: across multiple classes of customers (demand) for a utility (e.g., office, residential, retail) within a grid of local generation; across multiple grids of local generation, e.g.,within a city energy demand model; across multiple local generation distribution feeders; and so on within regional grids, across regional grids, within the national grid, across national grids (e.g., US and Canada). This generation including renewableenergy sources (solar, wind, etc) may be modeled for a single period in fine or course time scale (.about.sec, Minutes, hours, days etc.) In the model described, class is an abstraction that allows the ISO/utility to differentiate among different groupsof entities.<br><br>The methodology described accounts for heterogeneous customer flexibility in load reduction; and, a demand response 225 realized via a utility's rebate signal (not shown), as shown in FIG. 2.<br><br>Accordingly, in the description of the method, and provided as inputs to a programmed computer, there is defined:<br><br>N.sub.g is number of traditional generator buses;<br><br>N.sub.d is the number of load buses;<br><br>N is the total number of buses; and,<br><br>P.sup.d.sup.i is the power demand at each (load) bus.<br><br>In one embodiment, an assumption is made that power demand P.sup.d.sup.i at each bus "i" is normally distributed, P.sup.d.sup.i:N(.mu..sub.i,.sigma..sub.i.sup.2). This assumption is typically used in stochastic Optimal Power Flow (OPF) models,and in particular fits in well with a demand response model to be used.<br><br>The N.sub.g generators are peak generators, namely active generators that are currently unused but spinning and thus can be activated within a short time span, e.g., on the order of minutes. Let P.sup.g.sup.i.sup.min and P.sup.g.sup.i.sup.maxrepresent the lower and upper bounds on the power outputs from generators i, i=1, . . . , N.sub.g. Let .theta..sub.i.sup.min and .theta..sub.i.sup.max represent the lower and upper bounds on the voltage phase angles of bus i, i=1, . . . , N (where iis representative of bus or generator).<br><br>It is further assumed that renewable generation is connected to a single bus such that the total power generated lies within P.sub.RN.sup.min and P.sub.RN.sup.max, where the forecasts are P.sub.RN: (truncated) N(.mu..sub.RN,.sigma..sub.RN.sup.2) i.e., the values of the normally distributed random variable cannot be negative. Other renewable resources share similar intermittent and volatile properties with wind and can be modeled in this manner. It is further assumed thatthe utility also integrates a demand-response policy as a source of virtual generation, where M of the buses offer rebates to customers for reduced energy load. This can be implemented by allowing the utility to interact with different householdelectric appliances (laundry, temperature control, etc.) through price signals sent to smart meters as is known in the art ([e.g., http://www.smartmeters.com/the-news/617-us-northwest-perfect-for-- smart-grid-rollout.html). The customer response tothese incentives is a reduced demand P.sub.DR.sub.j(r) with a mean and variance depending on the rebate value r, P.sub.DR.sub.j(r):N(.mu..sub.DR.sub.j(r), .sigma..sub.DR.sub.j.sup.2(r)), j=1, . . . , M with j representing a bus subject of the loadreduction. In example numerical experiments, the mean value of the demand response, .mu..sub.DR.sub.j, is linear in the rebate value: .mu..sub.DR.sub.j=b.sub.1jr, j=1, . . . , M. (1)<br><br>Finally, any excess (shortfall) of supply is sold to (bought from) the spot market at a unit price c.<br><br>In the following, the subscripts .sub.RN and .sub.DR are used to indicate connections with wind generation (renewable energy) resource and demand response, respectively. In one embodiment, it is assumed that the total demand, demand responseand wind power generation are mutually independent random processes.<br><br>In one embodiment, the model includes formulating an optimization problem with the following decision variables:<br><br>z.sub.i is a 0/1 indicator of whether traditional generator i is on or off, i=1, . . . , N.sub.g;<br><br>P.sup.g.sup.i is a total real power output extracted from generator bus i, i=1, . . . , N.sub.g;<br><br>.theta..sub.i is a voltage phase angle for each bus i, i=1, . . . , N;<br><br>r is a non-discriminate unit rebate price offered to all customers on the M participating busses.<br><br>It is noted that {z.sub.1} are discrete control variables, while all the other decision variables are continuous. As for the continuous control variable r, the utility is assumed to have the forecasting profiles of each type of load reduction.<br><br>Constraints<br><br>In one embodiment, the problem is formulated as a single-period, nonlinear Mixed Integer Program (MIP) and can be solved using Optimization Subroutine Library OSL, CPLEX.RTM. (registered trademark of CPLEX Optimization, Inc. By single period,the optimization is run once, and the period may be of any time scale, e.g., minutes. Further, the optimization is adapted/expanded to multiperiod, analogous to optimization provided in herein-incorporated Ser. No. 12/890,119. In the model employed, areal power generation P.sup.g.sup.i of each generator bus, the load demand P.sup.d.sup.i of each load bus, and the voltage phase angles .theta..sub.i of each bus lie within their minimum and maximum limits. In addition, the total generation load plusthe reduced load balance the total demand. This therefore results in the following bounding constraints: z.sub.iP.sup.g.sup.i.sup.min.ltoreq.P.sup.g.sup.i.ltoreq.z.sub.iP.sup.g.s- up.i.sup.max, i=1, . . . , N.sub.g; (2).theta..sub.i.sup.min.ltoreq..theta..sub.i.ltoreq..theta..sub.i.sup.max, i=1, . . . , N; (3) b 0.ltoreq.r<c; (4) 0.ltoreq.P.sub.DR.ltoreq.P.sub.DR.sup.d; (5)<br><br>where the last constraint restricts the total demand-response load (variable P.sub.DR) to be no more than the load of the hosting bus (variable P.sub.DR.sup.d), where P.sub.DR.sup.d refers to total demand at the bus.<br><br>Further, in one embodiment, the formulation imposes power flow balance constraints on the decision variables. For this economic dispatching problem, the DC power flow constraints commonly used in optimal power flow and economic dispatchproblems, are used as set forth in equation 6):<br><br>.function..theta..times..times..theta..theta. ##EQU00001## where x.sub.ij is the reactance of the transmission line between bus i and bus j. The generator bus injected with renewable energy source, e.g., wind power, has a real load output ofP.sup.g.sup.i+P.sub.RN. The generator bus with demand response options has a real load output of P.sup.g.sup.i+P.sub.DR, and the generator bus with access to the spot market has a real load output of P.sup.g.sup.i plus (minus) the amount of load boughtfrom (sold to) the market.<br><br>Power balance also has to be achieved at the system level:<br><br>.times..times. ##EQU00002## where S represents the load position in the spot market. This results in: S=.SIGMA..sub.i=1.sup.NP.sup.d.sup.i-.SIGMA..sub.i=1.sup.N.sup.gP.sup.g.s- up.i-P.sub.RN-P.sub.DR, which, if positive, represents a shortfallin supply and, if negative, represents an excess of supply sold to the spot market on a particular bus. Risk Control<br><br>For a loss value l(u;v), which is a cost function of decision variables u and random variables v (e.g., representing realized demand, generation from renewable energy sources), the .beta.-Value-at-Risk (the .beta.-quantile of the distributioninduced on the loss function. The risk measure .beta.-Conditional Value-at-Risk with confidence level .beta. is defined as the expected value of the cost-based loss function l(u;v) conditional on loss being beyond the .beta.-th quantile. The.beta.-VaR and .beta.-CVaR values are denoted as .alpha..sub..beta.(u) and .phi..sub..beta.(u), where .alpha..sub..beta.(u).ident.min{L.sub.0.epsilon.R:Prob{l(u;v).ltoreq.L.su- b.0}.gtoreq..beta.}, (8).phi..sub..beta.(u).ident.E[l(u;v)|l(u;v).gtoreq..alpha..sub..beta.(u)] (9)<br><br>It has been noted that CVaR is a coherent risk measure; moreover, it is a convex function of the decision variables. The model may consider CVaR as the risk measure for the power generation portfolio. For example, let L be a pre-specified risk(loss) level. Then the risk control constraint is formulated as .phi..sub..beta.(u).ltoreq.L, which can be transformed to the following equivalent formulation under the assumption, in one embodiment, that v is normally distributed: .mu.(u)+( {squareroot over (2.pi.)}(1-.beta.)e.sup.[erf.sup.-1.sup.(2.beta.-1)].sup.2).sup.-1.sigma.- (u).ltoreq.L, (10) .mu.(u)=E[l(u;v)], .sigma.(u)=Var(l(u;v)), and error function:<br><br>.function..pi..times..intg..times.e.times.d ##EQU00003##<br><br>When v has a general distribution, the CVaR constraint can be approximated by sampling a collection of K paths, with v.sub.1, . . . , v.sub.K generated from the probability distribution of v:<br><br>.function..beta..times..times..function..ltoreq. ##EQU00004##<br><br>There is further defined the cost function l(z, P.sup.g, r, .theta.; P.sup.d, P.sub.DR(r), P.sub.RN) as a function of decision variables z.epsilon.R.sup.N.sup.g, P.sup.g.epsilon.R.sup.N.sup.g and r, and random variables P.sup.d.epsilon.R,P.sub.DR(r).epsilon.R.sup.M and P.sub.RN.epsilon.R. In the model, the loss function (total cost) equals the sum of all the generation costs minus the total revenue from the load position exposed to the spot market (positive when selling to the marketand negative when purchasing from the market).<br><br>In one embodiment, it is assumed a production cost f.sub.i, i=1, . . . , N.sub.g, for those operating units, e.g., generators, to be a quadratic function of the amount dispatched, namelyf.sub.i=a.sub.2iP.sup.g.sup.i.sup.2+a.sub.1iP.sup.g.sup.i+a.sub.0i. (13)<br><br>Where total real power output extracted from generator bus i is P.sup.g.sup.i. In addition, for example, it may be supposed that the start-up cost for load reduction is zero and that this virtual generator only incurs a "generation" cost linearin the rebate r. In general, there is no operational cost for wind renewables. Thus, the following form for the loss function is had in terms of the utility's total cost:<br><br>.function..theta..function..times..function..times..times..times..times..- times..times..times..times..function..times..times..times..function..times- . ##EQU00005##<br><br>The mean and variance of the loss induced by the random variables P.sup.d, P.sub.DR(r) and P.sub.RN are functions of z, P.sup.g and r. That is:<br><br>.function..function..theta..function..mu..function..theta..times..functio- n..times..times..times..times..times..times..times..mu..function..times..t- imes..times..mu..function..mu..times..mu..function..function..theta..funct-ion..sigma..function..theta..times..times..sigma..function..times..sigma..- sigma. ##EQU00006##<br><br>The formulation of the model is now provided. An objective of the formulation is to minimize the total cost (loss function), which yields a final formulation:<br><br>.theta..times..function..function..theta..function..times..times..times..- times..times..times..times..times..times..times..times..times. ##EQU00007##<br><br>Thus, in one aspect, the model addresses the impact that the unpredictable nature of renewable generation has on the network dispatching problem; and further addresses that demand response in the form of offering energy reduction incentives(rebates) is effective to control this stochasticity.<br><br>FIG. 11 depicts a computer-implemented method 100 for managing dispatch of energy over a network according to one embodiment. As shown in FIG. 11, there is depicted the step of 103 of receiving from wired or wireless communications and/oraccessing from a memory storage device, data representing the network configuration, any renewable energy injection, and bus power demand (a random variable). Further, there is received and or accessed, one or more decision variables: z.sub.i,P.sup.g.sup.i, .theta..sub.i, and, r. Further, there is received and or accessed, data representing the power flow constraints, decision variable bounding constraints, risk control constraints, including specified risk level value for CVar riskconstraint. Then, at 106, there is performed the programming, as a single period MIP, of the optimization problem including the cost-based function E[l(z, P.sup.g, r, .theta.; P.sup.d, P.sub.DR(r), P.sub.RN)] to be minimized which is a function of thedecision and random variables. Then, at 110, there is performed solving the optimization problem to jointly determine optimal dispatching policy for generation and an incentive rate for customers of the utility to extract load reduction on one or morebusses for implementation by the utility at step 130.<br><br>Towards this end, instances of the above nonlinear mixed integer program are solved. In the solution, many combinations of the specific values for the decisions variables are tested, and each combination requires a solution e.g., a linear ornonlinear optimization problem.<br><br>In one embodiment, FIG. 12 depicts an example method for solving said optimization problem indicated at step 110, FIG. 11, that implements a two-stage procedure: In the two-stage solution implemented: a first outer stage provides a listing ofpossible joint values of the binary variables z.sub.i<br><br>.di-elect cons..times. ##EQU00008## that indicate whether generator i is dispatched, and the inner stage where each z.sub.i has been fixed and the formulation simplifies to a non-linear formulation that is solved using a gradient descentalgorithm given the expressions for the non-linear constraint of equation (10).<br><br>Thus, as shown in FIG. 12, at 113 there is depicted the step of determining sets of all possible joint values of variables z.sub.i to indicate energy generator dispatch subject to input constraints. These sets are denoted as sets Y.sub.1, . .. , Y.sub.K where K is a set limit. In the method depicted, for example, considering index q a variable set Y.sub.q represents a set of values of variables z.sub.i's for all i. For example, given a topology including 2 generators, Y.sub.1 is {z.sub.1=0,z.sub.2=0}, Y.sub.2 is {{z.sub.1=1, z.sub.2=0}, Y.sub.3 is {z.sub.1=0, z.sub.2=1} and Y.sub.4 is {z.sub.1=1, z.sub.2=1}. Then, after setting index q as q=1, the solution proceeds to step 115 where variables z.sub.i are replaced by elements of setY.sub.q. In the loop, a next step 117 represents the performing of a gradient descent algorithm wherein for the q=1 case, there is solved the<br><br>.theta..times..function..function..theta..function. ##EQU00009## subject to constraints of eq. (2)-(7), (10). At step 117, any solution result from the gradient descent technique implemented may be stored in a memory storage device. Continuing in FIG. 12, the index q is incremented (i.e., by setting q=q+1) and the process proceeds to step 120 to determine if the limit last set Y.sub.K to be processed has been reached, i.e., whether q=K+1. In the example topology 50 of FIG. 2, sixgenerators may represent K=64 (2.sup.6) possible energy dispatch indicator combinations of (decision variables z.sub.i). If a maximum number of Y.sub.K sets has not been processed, the loop processing continues by fixing a new z.sub.i (e.g., i isincremented) at 123 and returning to step 115 for performing the next gradient descent processing on the next set Y.sub.q+1. If a maximum number of Y.sub.K sets has been processed as determined at step 120, then the process proceeds to step 125 todetermine the particular z.sub.i and Y.sub.q set values providing the optimal solution result obtained over all the iterations are output/stored. That is, in one embodiment, the process includes, at each iteration, comparing each solution determinedagainst solutions of all previous iterations of the loop, to track the most optimal z.sub.i and Y.sub.q set values.<br><br>In consideration of representative numerical experiments, the distribution network of an IEEE 30-bus 6-generator system topology shown in FIG. 2 was modeled. As shown in FIG. 2, the illustrative IEEE 30-bus 6-generator system 50 includesgenerators 35, and multiple generator busses (e.g., numbered busses 1, . . . , 30) operative interconnected with generators for carrying and distributing energy generated from the generators to customers via transmission lines (e.g. feeders such asshown in FIG. 1). Example parameter values for the quadratic cost function a.sub.0, a.sub.1, a.sub.2 are provided as are the upper P.sup.g.sup.max power bound and lower (power) bounds P.sup.g.sup.min specified constraints for generator busses areprovided in Table 1. Example voltage phase bounds .theta..sup.min and .theta..sup.max for all load busses are 0.degree. and 30.degree., respectively.<br><br>TABLE-US-00001 TABLE 1 Generator Unit Parameters for IEEE 30-Bus System Bus a.sub.0 a.sub.1 a.sub.2 P.sup.g.sup.max P.sup.g.sup.min 1 0 2.00 0.00375 200 50 2 0 1.75 0.00175 80 20 5 0 1.00 0.00625 50 15 8 0 3.25 0.0083 35 10 11 0 3.00 0.025 30 1013 0 3.00 0.025 40 12<br><br>A computing apparatus is programmed to solve the cost-based function of the model. The apparatus may use the optimization methods or optimization packages for solving Mixed Integer Programming problems, e.g., CPLEX.RTM., results in solutionoutputs, i.e., decision variables.<br><br>FIG. 3 is a plot 140 depicting, from solution results for the example topology 50 of FIG. 2, how the total expected cost when allowing demand response (DR) is lower than the total expected cost without DR for varying spot market priceassumptions. That is, as shown in plot 140 of FIG. 3, the computed total cost when calculated without implementation of DR is represented as plots 141a, 142a, 143a when priced according to various spot market prices, and each is shown plotted at adecreased cost as represented by example plots 141b, 142b, 143b when calculated with implementation of DR. In addition, under relatively high spot market prices, if adopted, the utility can even realize additional profits through the use of DR, which canlead to more efficient overall management of the energy grid.<br><br>Table 2 compares the optimal dispatching policy, rebate price and spot market load position S (e.g., positive to buy and negative to sell), for different spot market price and wind generation assumptions, both with and without DR. From Table 2,it can be seen that when the market price is relatively low, wind injection into the system will not change the dispatching load nor the rebate price, and will only reduce the amount of load bought from the spot market. In one example embodiment, anoptimal rebate price is about c/2. When the spot market price is relatively high, the utility dispatches more generation to realize a profit from this excess load. A dispatching order is inherent in the solution provided in results of Table 2.<br><br>From FIG. 3 plots 140 show how more wind injection under such conditions may reduce the dispatched thermal power and rebate price, which seems to result from the power flow constraints limiting the ability of the utility to transfer the excessgenerated load to the spot-market entry node. The optimal rebate price in this case is less than c/2. Under every condition, the total dispatched thermal generation is lower when DR is allowed.<br><br>TABLE-US-00002 TABLE 2 Optimal Dispatching Without CVaR Constraint c Wind P.sup.g.sub.1 P.sup.g.sub.2 P.sup.g.sub.3 P.sup.g.sub.4 P.sup.g.sub- .5 P.sup.g.sub.6 r S Cost 3 15 133.3 35.7 16.0 0 0 0 off 83.4 700.2 3 15 133.3 35.7 16.0 0 0 0 1.568.4 677.7 3 45 133.3 35.7 16.0 0 0 0 off 53.3 610.2 3 45 133.3 35.7 16.0 0 0 0 1.5 38.4 587.7 5 15 200 58.9 23.4 35 17.8 19.4 off -86.1 593.8 5 15 196.4 53.5 22 33.4 14.3 16.2 1.81 -85.5 559.7 5 45 188.3 51.9 21.6 30.6 13.3 15.2 off -82.5 488.7 5 45178.9 50.1 21.1 27.2 12.1 14.1 1.75 -82.8 457.5 5 80 167.6 47.9 20.6 23.2 10.7 12.8 off -82.5 372.0 5 80 158.2 46.1 20.1 19.9 10 12 1.68 -79.7 343.2 10 80 126.4 53.9 25.8 35 21.4 30.0 off -88.9 -52.5 10 80 115.0 51.7 25.2 35 19.9 28.6 1.78 -89.8 -84.8<br><br>Of the numerical examples, FIG. 4A illustrates a plot 150 of the mean cost, rebate price and load profiles as the tolerated risk levels vary and the spot market price is kept relatively low (low vs. high is relative to other costs). As seenfrom the plot 155, FIG. 4B, the dispatching policy and total dispatch load remains low (low vs. high is relative to capacity) for different CVaR constraints (one example measure), due to the relatively low spot market price. From these, the utility mayprefers to increase load reduction incentive prices to purchase less from the spot market instead of dispatching more generation load due to the higher cost for this extra dispatching.<br><br>FIG. 5 is a plot 160 depicting the optimal rebate and the minimum cost as the tolerated risk levels vary, thus providing the efficient frontier for the cost-based loss function.<br><br>FIG. 6 is a plot 170 of the optimal total dispatching amount and demand response load reduction with respect to these risk tolerances. It is seen from FIG. 6 that with relatively high spot market prices, the utility generates (either directlyfrom thermal sources or virtually via load reduction) and sells extra load to the market to make more profit. With a relatively loose CVaR constraint, the utility tends to pay more load reduction incentives rather than generate more energy from thermalsources.<br><br>FIGS. 7 and 8 show in respective plots 175, 180 that as the coefficient of variation (CV) for wind generation increases, the mean loss increases while the rebate price decreases under the same risk constraint. This is because with the same riskconstraint but a higher proportion of risk being due to wind generation, the most efficient way to keep CVaR within its prescribed limit is to reduce the risk from demand response, which in turn is achieved with a decreasing rebate. Moreover, thechanges in cost and rebate are more significant with a tighter risk constraint, which follows from the risk measure in the model being quadratic in r.<br><br>FIG. 9 illustrates a plot 190 of the optimal cost under changing market price relations with respect to different choices for the busses that interconnect wind, DR or spot market load injections. The optimal costs tend to increase with the spotmarket price at first, upon reaching a peak, and then tend to decrease because the utility can realize more profits by selling load to the market when the marginal generation cost is lower than the spot market price. In addition, a lower cost isincurred when the phase angle constraints are not active for the bus accessing the spot market or demand reduction.<br><br>FIG. 10 is a plot 195 that compares the optimal cost versus risk relations under different wind generation injection bus locations. Lower costs are incurred when the power flow constraints are not active, or when the bus that has wind injectionis also attached to a lower marginal cost generator. These results serve to illustrate and quantify the influence that the bus connections for renewable generation and spot market access can have on the resulting optimal dispatch.<br><br>The solution described assumes the existence of a feasible solution (given the constraints described), as demonstrated in the example numerical experiments. In one embodiment, the approach ensures feasibility by first considering a version ofthe model within the context of a robust optimization framework. Upon solving this problem, the corresponding robust optimal solution provides a lower bound on the potential benefits achievable by the utility, which the utility can decide to directlyemploy. On the other hand, this solution can instead be exploited to properly define a feasibility region for the original optimization problem and then apply a variant of the original solution approach taken over this feasible region. By employingthis solution, the utility can further improve upon the benefits realized under the optimal solution.<br><br>The importance of capturing in the dispatching problem the effects of any renewable generation sources introduced at various nodes in the distribution network is demonstrated. Although they generate energy at essentially no operational cost,renewable sources can create issues arising from supply forecast uncertainty, where a mismatch in the planned supply and actual demand due to lower than expected renewable generation results in the utility having to purchase expensive units from the spotmarket.<br><br>The approach by implementing the model hedges this uncertainty by imposing a risk constraint on the objective function (total cost of meeting demand) of the dispatching problem, where a risk metric based on the Conditional Value-at-Risk of thecost function is considered. This risk constraint has the effect of slightly increasing the total dispatched generation capacity to hedge the renewable generation uncertainty. In addition or as an alternative strategy, adopt a policy to incentivizeend-users to reduce their load, which can be effectuated almost instantly (as compared to increasing thermal generation) but at the expense of such incentives. The formulation also considers this strategy.<br><br>Results from numerical experiments show that the introduction of renewable sources can have a significant impact on the dispatching policy, as can the location of this injection within the transmission network. Additionally, demand responseload reduction is an effective tool in controlling shortfalls created by renewable generation forecast mismatches.<br><br>In further embodiments, a variant of the formulation considers AC power flow and incorporates into the model AC power flow equations, which yields a highly non-linear stochastic optimization problem.<br><br>Thus, the system and method of the invention optimizes the utility's total financial performance by providing the ability to intelligently hedge against energy shortfalls by dispatching various energy sources, including thermal, renewables, andvirtual generation, e.g., by determining the optimal dispatching policy and non-discriminate incentive rate for the utility.<br><br>Further, the system and method achieves improvements in both total financial measures and risk control over existing mechanisms that do not consider virtual generation and renewables. Moreover, the model incorporates risk-based control ofrenewable energy sources and demand response together with physical constraints on transmission of energy through the distribution network.<br><br>In this manner, resulting benefits include: integration of demand response as a virtual generator to hedge against the volatility and intermittency of renewable generator allows utilities to decrease the "spinning reserve" level. This is theadditional reserve generation capacity from thermal sources that has to be provisioned to match any shortfall due to load or generation fluctuation (e.g. from renewables). This is a high cost saving benefit; providing a dispatching policy with riskmanagement can provide higher coordinated, responsive control on utility's generation assets as well as the load from end-users. As shown in FIG. 1, this can be used in many levels of distribution hierarchy. Further, incentive mechanism providesignificant benefits to utilities and their customers.<br><br>FIG. 13 illustrates an exemplary hardware configuration of a computing system 400 running and/or implementing the method steps described herein. The hardware configuration preferably has at least one processor or central processing unit (CPU)411. The CPUs 411 are interconnected via a system bus 412 to a random access memory (RAM) 414, read-only memory (ROM) 416, input/output (I/O) adapter 418 (for connecting peripheral devices such as disk units 421 and tape drives 440 to the bus 412), userinterface adapter 422 (for connecting a keyboard 424, mouse 426, speaker 428, microphone 432, and/or other user interface device to the bus 412), a communication adapter 434 for connecting the system 400 to a data processing network, the Internet, anIntranet, a local area network (LAN), etc., and a display adapter 436 for connecting the bus 412 to a display device 438 and/or printer 439 (e.g., a digital printer of the like).<br><br>As will be appreciated by one skilled in the art, aspects of the present invention may be embodied as a system, method or computer program product. Accordingly, aspects of the present invention may take the form of an entirely hardwareembodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining software and hardware aspects that may all generally be referred to herein as a "circuit," "module" or "system." Furthermore,aspects of the present invention may take the form of a computer program product embodied in one or more computer readable medium(s) having computer readable program code embodied thereon.<br><br>Any combination of one or more computer readable medium(s) may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example,but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storagemedium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory),an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be anytangible medium that can contain, or store a program for use by or in connection with a system, apparatus, or device running an instruction.<br><br>A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated signal may take any of a variety of forms,including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, ortransport a program for use by or in connection with a system, apparatus, or device running an instruction. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless,wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.<br><br>Computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or thelike and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may run 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).<br><br>Aspects of the present invention are described below 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 understoodthat 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 program instructions. These computer program instructions may beprovided 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 run via the processor of the computer or other programmable dataprocessing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer program instructions may also be stored in a computer readable medium that can direct a computer,other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer readable medium produce an article of manufacture including instructions which implement the function/actspecified in the flowchart and/or block diagram block or blocks.<br><br>The computer program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other devices to cause a series of operational steps to be performed on the computer, other programmable apparatus or otherdevices to produce a computer implemented process such that the instructions which run on the computer or other programmable apparatus provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.<br><br>The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more operable instructions for implementing the specified logical function(s). It should also be noted that, insome alternative implementations, the functions noted in the block may occur out of the order noted in the FIG.s. For example, two blocks shown in succession may, in fact, be run substantially concurrently, or the blocks may sometimes be run in thereverse order, depending upon the functionality involved. It will also be noted 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 byspecial purpose hardware-based systems that perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.<br><br>While there has been shown and described what is considered to be preferred embodiments of the invention, it will, of course, be understood that various modifications and changes in form or detail could readily be made without departing from thespirit of the invention. It is therefore intended that the scope of the invention not be limited to the exact forms described and illustrated, but should be construed to cover all modifications that may fall within the scope of the appended claims.<br><br><center>* * * * *</center></td> </tr> </table> </td> </tr> </table> </td> <td align="left" valign="top" background="/web/20190415173352im_/http://www.patentgenius.com/images/box_rightvertical.gif"> </td> </tr> <tr> <td align="right" valign="top"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/box_bottomleft.gif" width="29" height="22"></td> <td background="/web/20190415173352im_/http://www.patentgenius.com/images/box_bottombg.gif"> </td> <td align="left" valign="top"><img src="/web/20190415173352im_/http://www.patentgenius.com/images/box_bottomright.gif" width="29" height="22"></td> </tr> </table></td> <td width="16%" align="center" valign="top" bgcolor="#05507f"><table width="100%" border="0" cellpadding="0" cellspacing="0" background="/web/20190415173352im_/http://www.patentgenius.com/images/bg_heading.gif"> <tr> <td width="15%"> </td> <td width="85%" height="25" align="left" valign="middle" class="heading1">Recently Added Patents</td> </tr> </table> <table width="97%" border="0" cellspacing="3" cellpadding="3"> <tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8713123.html" class="textblue">Method of sending CTI messages in a communication system</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/D703638.html" class="textblue">Antenna for an automobile</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8712044.html" class="textblue">Key management system</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8709999.html" class="textblue">Substituted bicyclic HCV inhibitors</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8711030.html" class="textblue">Single-pass Barankin Estimation of scatterer height from SAR data</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/D703698.html" class="textblue">Lower set insert with a lower ball seat for a downhole plug</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8713078.html" class="textblue">Method for building taxonomy of topics and categorizing videos</a></td></tr><script src="https://web.archive.org/web/20190415173352js_/http://www.google-analytics.com/urchin.js" type="text/javascript"> </script> <script type="text/javascript"> _uacct = "UA-801222-1"; urchinTracker(); </script> </table> <table width="100%" border="0" cellpadding="0" cellspacing="0" bgcolor="000000"> <tr> <td width="100%" height="3px"></td> </tr> </table> <table width="100%" border="0" cellpadding="0" cellspacing="0" background="/web/20190415173352im_/http://www.patentgenius.com/images/bg_heading2.gif"> <tr> <td width="15%"> </td> <td width="85%" height="25" align="left" valign="middle" class="heading1">Randomly Featured Patents </td> </tr> </table> <table width="97%" border="0" cellspacing="3" cellpadding="3"> <tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/8112668.html" class="textblue">Dynamic broadcast of configuration loads supporting multiple transfer formats</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/6540532.html" class="textblue">Electrical connector assembly for connecting electrical contacts</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/4395346.html" class="textblue">Method for cleaning contact lenses</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/4173026.html" class="textblue">Self clocking speed tolerant magnetic recording method and apparatus</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/7397503.html" class="textblue">Systems and methods for reducing artifacts caused by illuminant flicker</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/7460826.html" class="textblue">Image forming apparatus</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/7010590.html" class="textblue">System and method for secure transactions over a network</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/4674623.html" class="textblue">Bucket elevator including means for modifying bucket speed in the loading zone</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/6731251.html" class="textblue">Device for supporting two interfaces including the TV image card and the LCD panel image card</a></td></tr><tr><td class="dots2"><a href="https://web.archive.org/web/20190415173352/http://www.patentgenius.com/patent/7627379.html" class="textblue">Method of speech conversion in a cochlear implant</a></td></tr> </table> </tr> <tr> <td height="30" align="center" valign="middle" bgcolor="#EFEFEF" class="copyright">© 2006-2011. PatentGenius.com All Rights Reserved.</td> <td bgcolor="#03324E"> </td> </tr> </table> </td> </tr> </table> </body> </html> <script type="text/javascript">displayImages('8626353',5,1);</script>