Systemic Simulation of Smart Grids : First Results

<div dir="ltr" style="text-align: left;" trbidi="on"><w:sdt contentlocked="t" id="89512093" sdtgroup="t"><span style="font-family: "Calibri","sans-serif"; font-size: 1.0pt; mso-ansi-language: FR; mso-ascii-theme-font: minor-latin; mso-bidi-font-family: "Times New Roman"; mso-bidi-language: AR-SA; mso-bidi-theme-font: minor-bidi; mso-fareast-font-family: "Times New Roman"; mso-fareast-language: FR; mso-fareast-theme-font: minor-fareast; mso-hansi-theme-font: minor-latin;"><w:sdtpr></w:sdtpr><w:sdt docpart="60B50BC023144821ABA7290EC34FE67B" id="89512082" storeitemid="X_66CC8B3F-D460-4797-8657-E0B3B1ECE80D" text="t" title="Titre du message" xpath="/ns0:BlogPostInfo/ns0:PostTitle"></w:sdt></span> </w:sdt><br /><div class="Publishwithline"><br /></div><div style="border-bottom: solid #4F81BD 1.0pt; border: none; mso-border-bottom-themecolor: accent1; mso-element: para-border-div; padding: 0cm 0cm 2.0pt 0cm;"> <div class="underline"><br /></div></div><div class="PadderBetweenControlandBody" style="text-align: justify;"><span style="font-family: Verdana, sans-serif;"><span style="font-size: 10pt; line-height: 12pt; text-indent: 14.2pt;">The month of August provides me with the opportunity to get back to my computer and to resume my programming projects. I have been able to complete the first step of S3G (</span><a href="http://informationsystemsbiology.blogspot.fr/2011/09/systemic-simulation-of-smart-grids.html" style="font-size: 10pt; line-height: 12pt; text-indent: 14.2pt;">Systemic Simulation of Smart Grids</a><span style="font-size: 10pt; line-height: 12pt; text-indent: 14.2pt;">) which I started last year. I had the opportunity to attend the </span><a href="http://www.naefrontiers.org/Symposia/EU-USFOE/32025/2011EU-USFOE.aspx" style="font-size: 10pt; line-height: 12pt; text-indent: 14.2pt;">EU-US Frontiers of Engineering Symposium</a><span style="font-size: 10pt; line-height: 12pt; text-indent: 14.2pt;">, where I have presented my project and received encouraging and interesting feedback.</span></span></div><div class="MHeading2" style="line-height: 12pt; margin: 0cm 0cm 0.0001pt; text-align: justify;"><span lang="EN-US" style="font-size: 10pt;"><span style="font-family: Verdana, sans-serif;">The objective of S3G is to simulate the production and consumption of electricity throughout a long period of time (15 years). The S3G model, which is illustrated by the following figure, may be summarized with five parts:<o:p></o:p></span></span></div><div class="MHeading2" style="margin: 0cm 0cm 0.0001pt;"></div><ul style="text-align: left;"><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Energy demand: for each city, energy demand is generated from an hour-by-hour and day-by-day template, adding some random variation (the extent of which is a model parameter) together with a city-specific variation. This number is then reduced by the amount of “negaWatts”, computed from the total amount invested by the city. The model uses a ratio obtained from a concave-increasing function of the investment.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Dynamic Pricing: both suppliers and operators use a simple affine pricing model, with a constant price when the demand is less than a “base power”, and a linear formula when the demand is higher.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Production: suppliers use nuclear power according to planned schedule and adjust to resulting demand with fossil plants. Operators always use their green power (store it in the “buffer” or resell it when there is too much of it). They adjust to the city demand with their own fossil plant and wholesale electricity from suppliers, at the lowest marginal cost.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;"><span style="font-size: 7pt; line-height: normal;"> </span></span><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Consumption: The actual electricity consumption for each city is the demand, minus “shaving”, which is obtained by applying an S-curve to the sale price.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Market-share: for each city, the market balance between the national supplier and the local operator is determined yearly using another S-curve.</span></li></ul><br /> <div class="MHeading2" style="line-height: 12.0pt; margin-bottom: .0001pt; margin: 0cm;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjyXOwvxedvl-ED4CIdQZNX6S3gAieOSMNzlkQhG6iRuWOCVs_cL9nDdUg1z0DzACJZizEmAPi8HgKQ_hPMLaeuqAt9ndyqo_MfABSmk2fTOHmBXwtEEkUQPaEpeit9gQDquQrEPeU3cd4/s1600/Diapositive1.JPG" imageanchor="1" style="clear: left; float: left; margin-bottom: 1em; margin-right: 1em;"><img border="0" height="480" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjyXOwvxedvl-ED4CIdQZNX6S3gAieOSMNzlkQhG6iRuWOCVs_cL9nDdUg1z0DzACJZizEmAPi8HgKQ_hPMLaeuqAt9ndyqo_MfABSmk2fTOHmBXwtEEkUQPaEpeit9gQDquQrEPeU3cd4/s640/Diapositive1.JPG" width="640" /></a><span style="font-family: Verdana, sans-serif;"></span></div><div style="text-align: justify;"><span style="font-family: Verdana, sans-serif;"><span lang="EN-US" style="font-size: 10pt; line-height: 12pt;">This S3G model is both simple and complex. It is simple because it is based on a handful of equations, resulting in a simulation code that is 500 lines long. On the other hand, it is a complex model for two reasons. On the one hand, there are multiple feedback and interaction loops that make it difficult to analyze how the system will react to perturbations. On the other hand, there are many unknown parameters in this model (such as market sensitivity, demand-response behavior, negaWatt capabilities, etc.).</span><span lang="EN-US" style="font-size: 10pt; line-height: 12pt;"> </span></span></div><br /><div class="MHeading2" style="line-height: 12.0pt; margin-bottom: .0001pt; margin: 0cm;"><span style="font-family: Verdana, sans-serif;"><span lang="EN-US" style="font-size: 10pt;"><br /></span></span></div><div class="MHeading2" style="line-height: 12.0pt; margin-bottom: .0001pt; margin: 0cm;"><span style="font-family: Verdana, sans-serif;"><span lang="EN-US" style="font-size: 10pt;"><br /></span></span></div><div class="MHeading2" style="line-height: 12.0pt; margin-bottom: .0001pt; margin: 0cm;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiBur6yqBxsJcyAB2eYRNSmsgScw8FBuTsIMEBmmAKPFEaX-Y4DuVG9UtLoTf4s3wqnA7XMuqebL4kpT8tfuHnRUlBNeHif7ano-qquIdmkWWtCOO138KDHjnMybAqj0LDTZLOF1vAQd40/s1600/S3G-FoE-Poster-Dec12.jpg" imageanchor="1" style="clear: left; float: left; margin-bottom: 1em; margin-right: 1em;"><img border="0" height="480" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiBur6yqBxsJcyAB2eYRNSmsgScw8FBuTsIMEBmmAKPFEaX-Y4DuVG9UtLoTf4s3wqnA7XMuqebL4kpT8tfuHnRUlBNeHif7ano-qquIdmkWWtCOO138KDHjnMybAqj0LDTZLOF1vAQd40/s640/S3G-FoE-Poster-Dec12.jpg" width="640" /></a><span style="font-family: Verdana, sans-serif;"><span lang="EN-US" style="font-size: 10pt;">The following figure shows, on the right part, a rough summary of the simulation loop that is run for each time period (3 hours, hence 8 times per day). On the left part, it shows the three main <a href="https://sites.google.com/site/yvesresearchagenda/home">GTES generic procedures</a>.<o:p></o:p></span></span></div><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiQoHEooBvOHZnz7dU7e-xTRPQvxwC5RQUOTf9jk3kji9Oo9pR-enFz9kSUFeYuOvkN7wRSubPBlQAo8WClqhQ45GSkFev26fHMSvXZdSCOMwQfj6FQsuhQO4VnpJdCzLF0RTCyBd-qLJ0/s1600/S3G-simulation-loop.jpg" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="480" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiQoHEooBvOHZnz7dU7e-xTRPQvxwC5RQUOTf9jk3kji9Oo9pR-enFz9kSUFeYuOvkN7wRSubPBlQAo8WClqhQ45GSkFev26fHMSvXZdSCOMwQfj6FQsuhQO4VnpJdCzLF0RTCyBd-qLJ0/s640/S3G-simulation-loop.jpg" width="640" /></a></div><div class="MText" style="line-height: 12pt; text-align: justify;"><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">I have already spent some time on computational experiments that will be presented at </span><a href="http://www.csdm2012.csdm.fr/" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">CSDM</a><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">in December. </span><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">A GTES simulation run returns the average and standard deviation of a few key business parameters, as well as some indication of the </span><a href="http://organisationarchitecture.blogspot.fr/2008/01/lapproche-gtes-et-les-quilibres-entre.html" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">Nash convergence</a><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">. Giving averages and a few deviations is a poor restitution of the rich data gathered during the computational experiment, but the goal is simply to “get a feeling for what is happening”, as opposed to producing a forecast. The current set of experiments is designed to understand the main issues that were exposed in </span><a href="http://informationsystemsbiology.blogspot.fr/2011/09/systemic-simulation-of-smart-grids.html" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">the previous post</a><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;">. I have defined a fictional country somehow similar to France decomposed into 10 regions/cities. I have run 8 experiments that may be defined as follows:</span></div><div class="MText"></div><ol style="text-align: left;"><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The “default” is a reference point, from which “what-if” sensitivity analysis is made. The economic parameters are set in such a way that alternate operators start with a 20% market-share and should be able to increase it if they demonstrate a better management of variability.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The second experience raises the variability of energy consumption (globally), while the third experience raises the local variability (each city is more different from each other).</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The fourth experiment doubles the fossil energy price (gas and coal). In the default scenario, it is randomly drawn between 20€ and 40€/MWh.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The fifth experiment imposes a 5% reduction of the nuclear assets for the supplier during the first 5 years.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The sixth experiment sets a carbon tax at 100€/t, the proceeds of which is used by the “regulator” to subsidize green energy investment.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The seventh experiment explores the impact of overall demand variation in the next 15 years. The model assumes a constant growth/decline of electricity demand which is expressed as a percentage. This experiment plays with different values to see if <i>demand impacts the profitability of smart grids</i>.</span></li><li><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">The last experiment is a small variation of the first one, where wholesale prices are more rigidly constrained. During this set of experiment, I considered that the constraints that regulate wholesale prices for the supplier are fixed. An interesting next step will be </span><i style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">to make them a strategic parameter for the regulator</i><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;"> (which is a better representation of reality and makes for a more interesting game).</span></li></ol><div style="text-indent: -23.799999237060547px;"><span style="font-family: Verdana, sans-serif; font-size: x-small;"><span style="line-height: 16px;"><br /></span></span></div><span style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt;"><div style="text-align: justify;"><span style="font-size: 10pt; line-height: 12pt;">I still need to do some additional tuning before presenting my results at CSDM. I will then make my slides available on this blog. For once, I </span><i style="font-size: 10pt; line-height: 12pt;">need to run longer simulations</i><span style="font-size: 10pt; line-height: 12pt;"> (with larger samples) to ensure that the Monte-Carlo sampling is stable. I also want </span><i style="font-size: 10pt; line-height: 12pt;">to introduce better local optimization meta-heuristics</i><span style="font-size: 10pt; line-height: 12pt;">, such as Tabu search. One of the key issues for smart grid is the pricing model, both for the operator and the supplier. The more complex the pricing model, the better the actor may play with the </span><a href="http://en.wikipedia.org/wiki/Demand_response" style="font-size: 10pt; line-height: 12pt;">demand-response</a><span style="font-size: 10pt; line-height: 12pt;">concept. However, a more complex princing model requires better optimization techniques to ensure that the automatic computation from GTES (where each actor “learns” its “best response”) is relevant.</span></div></span><br /> <div class="MText" style="line-height: 12.0pt;"><br /></div><div class="MText" style="line-height: 12.0pt; text-indent: 0cm;"><span lang="EN-US" style="font-size: 10.0pt; mso-fareast-language: EN-US;"><span style="font-family: Verdana, sans-serif;">This being said, here are a few findings that may be drawn from the hundreds of runs made with the S3G model, and that are worth sharing: <o:p></o:p></span></span></div><div class="MText" style="text-indent: 0cm;"></div><ul style="text-align: left;"><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;"><span style="font-size: 7pt; line-height: normal;"> </span></span><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -17.85pt;">There is a systemic benefit of distribution and autonomy to cope with variation. This is shown in our second and third line. It is a “subtle” variation (small effects), which means that the economy of the local operator is dominated by its capacity to operate at much lower customer management costs than the supplier.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">CO2 tax increases play a very small role, and one that is difficult to anticipate since it both favors the local operator (support green subsidies) and the supplier (raises the difference between fossil and nuclear).</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">“De-nuclearization” is a favorable scenario for smart grid operators, as are most regulations that are adverse to the supplier. The obvious limitation is the resulting price increase that reduces the total economy output (and the country’s competitiveness).</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">The “community advantage” (that is, the ability for a local operator to better manage the demand-response loop because it is “closer” to its end customer) is marginal, and it is quite unclear if the payback from demand-response management is enough to sustain the operator’s business model.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">Investing in local storage is never an interesting option (at current prices). We needed to slash the price by over an order of magnitude to see a viable payback in less than 10 years.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">There is a clear competition between local operators and suppliers. The learning component of GTES makes for “agile” players who react closely to each other signals. The pricing structure plays an important role (we have only explored a simple variable pricing scheme). A logical consequence is the importance of regulation.</span></li><li><span lang="EN-US" style="font-family: Verdana, sans-serif; font-size: 10pt; line-height: 12pt; text-indent: -18pt;">The results are sensitive to the strategies of the player. <i>A next step for S3G is to build a “strategy matrix”</i>, which is a tabular “what-if” sensitivity analysis where we see what happens if the goals of the players are changed.</span></li></ul><br /> <div class="MsoNormal"><br /></div></div>

Systemic Simulation of Smart Grids : First Results

The month of August provides me with the opportunity to get back to my computer and to resume my programming projects. I have been able t...

Read more »