arxiv: v1 [math.ds] 22 May 2012

Transcription

1 FTT:Power : A global model of the power sector wth nduced technologcal change and natural resource depleton Jean-Franços Mercure Cambrdge Centre for Clmate Change Mtgaton Research (4CMR), Department of Land Economy, Unversty of Cambrdge, 19 Slver Street, Cambrdge, CB3 1EP, Unted Kngdom arxv: v1 [math.ds] 22 May 2012 Abstract Ths work ntroduces a model of Future Technology Transformatons for the power sector (FTT:Power), a representaton of global power systems based on market competton, nduced technologcal change (ITC) and natural resource use and depleton. It s the frst component of a famly of sectoral bottom-up models of technology, desgned for ntegraton nto the global macroeconometrc model E3MG. ITC occurs as a result of technologcal learnng produced by cumulatve nvestment and leads to hghly nonlnear, rreversble and path dependent technologcal transtons. The model uses a dynamc coupled set of logstc dfferental equatons. As opposed to tradtonal bottom-up energy models based on systems optmsaton, such dfferental equatons offer an approprate treatment of the tmes and structure of change nvolved n sectoral technology transformatons, as well as a much reduced computatonal load. Resource use and depleton are represented by local cost-supply curves, whch gve rse to dfferent regonal energy landscapes. The model s explored for a sngle global regon usng two smple scenaros, a baselne and a mtgaton case where the prce of carbon s gradually ncreased. Whle a constant prce of carbon leads to a stagnant system, mtgaton produces successve technology transtons leadng towards the gradual decarbonsaton of the global power sector. Keywords: Energy technology model, Clmate change mtgaton, Induced technologcal change 1. Introducton The future level of anthropogenc greenhouse gas (GHG) emssons s one of the prmary unknowns n estmatng the rate of clmate change n the medum and long term. Dfferent assumptons about emssons pathways lead to wdely dfferng warmng temperatures, rangng between 2 C, for hgh mtgaton polces to 6-8 C for the most pessmstc non mtgaton scenaros (IPCC, 2007). Anthropogenc GHG emssons depend on human actvtes, and thus on the structure of the future economc system. They stem prmarly from energy use (57%) and land use change (17%) (IPCC, 2007, fgures for 2004). Whle the second s a complex subject to model (see for nstance Bouwman et al. (2006)), the frst s better understood and dfferent approaches have been used to model energy demand and supply (IEA, 2010b, Bouwman et al., 2006, Messner and Strubegger, 1995, Seebregts et al., 2001). The problem of GHG emssons mtgaton requres changes to be made to the structure of global energy use. Snce the latter les at the very core of the world s economy, these changes have deep mplcatons and effects felt n all aspects of socety. Therefore, a smulaton of GHG emssons cannot easly be separated from smulatons of the global economc system. One approach to explore Energy-Economy- Envronment (E3) nteractons was made wth an nput-output Correspondng author: Jean-Franços Mercure Emal address: jm801@cam.ac.uk (Jean-Franços Mercure ) structure based macroeconometrc model of the global economy E3MG (Barker et al., 2006, Dagoumas and Barker, 2010, Barker and Screcu, 2010). E3MG s a dsaggregated model that features 20 world regons, 12 energy carrers, 19 energy users, 28 energy technologes, 14 atmospherc emssons and 42 ndustral sectors. The decarbonsaton of the global power system depends frst and foremost on the rate at whch hghly emttng technologes based on fossl fuels can be substtuted for cleaner ones. Whle fossl fueled electrcty generaton technologes are mature and well establshed, cleaner systems such as renewable energes are more expensve and therefore not yet compettve. Polcy favourng green technologes are expected to play an mportant role n the transton towards low GHG emssons n the energy sector. The costs of technologes tend to decrease wth cumulatve nvestment, an effect termed technologcal learnng, and moreover, newer technologes see ther costs decrease more rapdly than those of mature technologes (IEA, 2000, Pan and Koehler, 2007, McDonald and Schrattenholzer, 2001). R & D nvestments and government subsdes to green technologes are expected to be necessary only for an ntal push, untl ther reducng costs reach compettve levels, gvng rse to technologcal change (IEA, 2000). The composton of future energy landscapes are hghly dependent on the development n tme of polcy portfolos, where some strateges can generate technologcal avalanches whch would not have occurred otherwse. Energy landscapes are, however, very local n nature and technologcal substtuton occurs accordng to what nteractons Preprnt submtted to Energy Polcy January 24, 2014

2 are possble n a partcular world regon. The energy portfolo of each regon s n large parts determned by the avalablty of natural resources. Therefore, any modellng attempt of technologcal substtuton should reflect the local nature of natural resources. Ths paper ntroduces the model FTT:Power (Future Technology Transformatons n the Power sector), a dynamc model of the global power sector based on market prce competton, technologcal substtuton and resource use and depleton. It s the frst member of a famly of bottom-up dynamc technology substtuton models, bult as components for E3MG. FTT:Power makes use of a dynamc set of coupled logstc dfferental equatons, smlar to those used to represent populaton growth (Verhulst, 1838) and competton n bologcal systems (Lotka, 1925, Volterra, 1939). Evolutonary dynamcs have been used to descrbe technologcal change theoretcally (see for nstance Dercole et al., 2008). Technology transtons have been shown numerous tmes to follow emprcally logstc functons (for a revew see Grubler et al., 1999), n partcular n energy systems (Marchett and Nakcenovc, 1978). There exsts an extensve lterature where logstc systems of equatons have been used n the analyss of growth and competton n markets (for nstance Bass, 1969, Sharf and Kabr, 1976, C., 1989, Morrs and Pratt, 2003). Whle these systems follow closely logstc functons, ther underlyng changeover tmescales have never been explaned satsfactorly by any theory. The subject has been explored more recently by Anderson et al. (2007), where they ntroduced a famly of coupled dfferental logstc equatons to represent energy technology transtons. Ther system s not dynamc, however, and leads to statc results when technology costs do not change, and furthermore lack approprate symmetry propertes. Meanwhle, current energy technology substtuton models generally fal to reproduce the sgmod (S -shaped) character of technologcal transtons that stem from ther dynamc nature. Most models use nstead procedures of energy systems optmsaton based on cost mnmsaton, whch produces optmal equlbrum solutons. They lack the ntrnsc tme constants (of the order of 50 to 100 years) assocated wth technologcal transtons, and smoothng of changes n tme are generally mposed exogenously. The FTT:Power model s the frst power sector model based on a dynamc set of coupled logstc dfferental equatons, of the Lotka-Volterra famly, whch represents gradual technologcal substtuton processes. It s parametrsed by natural tme constants whch determne the rates at whch technologcal transtons can occur. It uses a probablstc treatment of nvestor decson-makng usng local dstrbutons of cost values and resultng lkelhood of technologcal swtchng. It allows for non-ratonal nvestor behavour, assumng that the trend rather than ndvdual actons s orented towards choosng lowest cost technologes. Due to the long tme constants nvolved, whch are assocated wth lfetmes and constructon tmes, however, the model does not reach an equlbrum state as t would n an optmsaton of the system for cost mnmsaton, unless all parameters were to reman stable for perods several tmes longer than these tme constants. Due to the smplcty of evaluatng n tme steps a set of dfferental equatons, 2 FTT:Power takes a small fracton of the computatonal tme requred by optmsaton models; t s thus lkely to be one of the most compact and fastest models of technology around. Addtonally, whle technology learnng curves generate multple optmal ponts to optmsaton models and assocated computatonal dffcultes, they do not requre addtonal computatonal tme n ths model. The model s desgned to calculate n parallel the dynamc evoluton of local energy landscapes n 20 world regons based on local cost values. These are constructed usng a standard framework whch ncludes components representng polcy decsons such as taxes and carbon markets, experence curves and the avalablty of natural resources and ther cost. The treatment of natural resource use and depleton s nspred by that featured n the TIMER model, whch defnes global sets of local cost-supply curves (Bouwman et al., 2006, Hoogwjk, 2004, Hoogwjk et al., 2004, 2005, 2009, Rogner, 1997), but expanded to nclude a representaton of uncertanty and evoluton n resource assessments. As an ntegral part of E3MG, FTT:Power s desgned to reproduce feedback nteractons between the global economy and the energy sector, whch occurs through the energy demand, tself hghly nfluenced by energy prces. Electrcty prces are derved through FTT:Power wth the margnal costs of each technology, whch depend on those of natural resources and technology. Ths paper s dvded nto three sectons. The logstc model of technology substtuton n dfferental equaton form and the assocated probablstc approach to local cost dstrbutons are frst ntroduced. We dscuss the effects of technologcal learnng onto cost values, and thus to decson makng, gven local polcy portfolos. We then descrbe how ths equaton s constraned usng a smple representaton of technologcal lmts mposed by local grd networks. The second secton ntroduces our framework of natural resource use and depleton featurng an explct representaton of uncertanty n global resource dstrbutons and avalablty. Ths s fnally followed by the presentaton of one set of global results n order to explore and understand the propertes of FTT:Power. For ths purpose, the model was run for ths paper by tself wthout feedback wth E3MG, n order to avod the sgnfcant ncrease n complexty of the model when feedback effects are present. The results were obtaned usng two sets of smplfed global polcy portfolos, a baselne and a mtgaton scenaro. We dscuss the appearance of a phenomenon we name the energy technology ladder, whch emerges from the complexty of the set of equatons at the root of FTT:Power. Ths effect s a clear representaton of rreversble and path dependent nduced technologcal change and the way n whch the energy sector mght evolve n hgh mtgaton scenaros n order to acheve hgh levels of decarbonsaton effcently at low costs. Ths result stems from our dynamc approach rather than beng mposed through a cost mnmsaton procedure and assocated assumpton of market equlbrum, and s consstent wth the non-equlbrum Post- Keynesan approach of E3MG. Addtonal nformaton regardng FTT:Power can be found n Mercure (2011).

3 2. Dynamcs as a set of dfferental equatons 2.1. The shares equaton Fgure 1: Top Probablty dstrbutons f (C, C ) and f j (C, C j ) of cost values for two technologes and j based on real recently experenced costs. The number of unts of technology j that come out at a cost cheaper than the medan value of the dstrbuton of technology corresponds to red shaded area, a value much smaller than that of the reverse, whch s the blue area. Bot. Cumulatve probablty dstrbuton functons for and j. We ntroduce the central assumpton of our model by presentng a probablstc framework desgned to represent the results of local decson-makng by nvestors n the electrcty sector. The electrcty sector sees a local varyng demand for power to whch t reacts by atttudes regardng the buldng of new or mantanng exstng power statons. On the ground, nvestors are faced wth a choce between several avenues, and they perform a comparson of the Levelsed Cost of Electrcty (LCOE, defned n secton 2.3). Ths cost, however, s not strctly well defned, snce a large number of aspects of very local nature affect the cost of specfc projects, such as the length of power lnes or the cost of land, and decsonal ssues whch may or may not lead to a ratonal choce of the opton wth the lowest LCOE. We therefore take costs as beng dstrbuted when examned from an aggregated vewpont at the regonal or natonal level, n order to avod lookng nto the detals of these local aspects. Therefore, even f one partcular technology s on average less expensve than another, there are almost always ndvdual cases where t s the reverse that s true, and therefore a probablstc treatment s approprate. Thus we assume that t s the trend rather than ndvdual decsons whch leans towards technologes wth lower LCOEs. We defne f (C, C ) as a probablty dstrbuton of the cost of technology based on real recent nstances of the cost of new unts (taken from IEA (2010a)), whch has a medan value C. The probablty, when buldng a new unt, that t comes at a cost lower than the value C s the cumulatve probablty dstrbuton F (C, C ). For two technologes and j, the fracton of new unts 3 of j whose cost turn out lower than the medan value of s F j (C, C j ), whle the fracton of unts of whose cost s lower than the medan of j s F (C j, C ). These can also be expressed n terms of the cost dfference C j. Ths s shown n fgure 1. Usng these arguments, the dfferental shares equaton can be constructed. We model technology swtchng usng a parwse comparson of costs, and flows of unts of market shares S of dfferent technologes, that s, the shares of electrcty generaton capacty. We wll demonstrate that ths s equvalent to a smultaneous comparson of the costs of all optons. We frst determne the lkely number of energy producton unts flowng from category j towards category, denoted S j. The number of unts of market share that can flow away from technology j s a fracton of the exstng fleet, and therefore s proportonal to the market share tself S j. The rate, therefore, at whch power plants of category j wll be closed down, f none are scrapped before the end of ther lfetmeτ j, s S j /τ j. Meanwhle, the rate at whch technology can grow s also proportonal to ts share of the market. Ths s smlar to populaton growth: the technology s assumed to breed tself, an assumpton whch reflects the growng ablty of an ndustry to expand as ts sales ncrease, and forms the bass of the logstc based market competton lterature (Bass, 1969, Sharf and Kabr, 1976, C., 1989, Morrs and Pratt, 2003). The rate at whch technology can grow s nversely proportonal to the technology specfc buldng capacty expanson tme constant t, and s therefore S /t. In a unt of tme dt, where t s measured n years throughout, we assume that the number of unts flowng from j to s proportonal to the probablty that the cost of s less than that of j, F (C j, C ), and conversely, leadng to S j S j S j τ j S t F ( C j ) t S τ S j t j F j ( C j ) t. These aggregate movements are completely ndependent of each other: they stem from every ndvdual substtuton event where nvestors make a specfc choce and where costs dffer for local reasons whch we gnore here. There results a net flow, whch s the dfference between these two terms, denoted S j, and represents every substtuton event that occurs durng the tme nterval t. In order to fnd the total number of unts flowng towards from all other technologes, we perform the sum of S j over j, whch yelds the shares equaton central to ths model: ( S = S S j A j F ( C j ) A j F j ( C j ) ) t, (1) j where A j = K/τ t j s a matrx of substtuton frequences, whch s not symmetrcal, wth K a tme scalng constant. 1 Ths matrx contans all the natural tme constants assocated wth nduced technologcal change. Its detaled dervaton nvolves a lengthy calculaton that cannot be ncluded here but wll be the subject of a forthcomng publcaton. Ths theory yelds a changeover (transton) tme, whch depends on the strength

4 of nvestor preferences F ( C j ), and ranges between the shortest possble decommsson tmeτ(f ( C j )=1) and nfnty (F ( C j )=0). Note that A j s not symmetrcal, and ths stems from the fact that rates of technology uptake and rates of decommsson are not dentcal across technologes. For example, technology beng decommssoned at rateτ 1 and replaced by technology j at a rate proportonal to t 1 j does not occur at the same rate as the reverse process where technology j s decommssoned, at rateτ 1 j, and replaced by at rate proportonal to t 1. Addtonally, t can be seen that the equaton s nvarant over permutatons of and j, requred n order to have j ds j = 0 and S = 1. It can be shown numercally that the result does not depend on the order over whch the calculaton s performed. Ths means that, f the same result s obtaned rrespectve of the order over whch the parwse comparsons are done, that the total operaton corresponds to a smultaneous comparson of all technologes. Ths s ensured by the dfferental form of eq. (1), where exchanges are performed over nfntesmal tme ntervals, and s further demonstrated wth examples n secton Induced Technologcal Change One of the most mportant drvers of change n the electrcty sector s the cost reducton that stems from technologcal learnng, and have been ncluded n several energy models, notably n the precedng energy submodel of E3MG (Koehler et al., 2006). It s well establshed that the repettve producton of goods gves rse to mprovements n producton methods and economes of scale that lead to cost reductons, and these depend drectly on the number of unts produced (Arrow, 1962, Grubler et al., 1999). Moreover, the functonal form of the experence curve s very smple, t s that of a power law of whch the negatve exponent, denoted b here, s related to the progress rate (IEA, 2000, Berglund and Soderholm, 2006). Thus, the logarthm of the cost of producton C relatve to an ntal value depends lnearly on the logarthm of the total number of unts W sold snce the very frst one came out of the factory. In lnear terms, the experence curve s expressed as ( ) b W (t) C (t)=c 0,, (2) W 0, where C 0, s the cost assocated wth the cumulatve number of unts W 0, produced up to an arbtrary startng tme t=0. Values for W 0 are not trval to fnd n the cases of old technologes snce they nvolve the number of unts currently operatng but also the total number of unts whch have been abandoned or demolshed; learnng s however mnmal n such cases. Where learnng s mportant, wth newer technologes, these correspond approxmately to the current number of nstalled unts snce most of them are stll n operaton. The learnng process can also be seen n an ncremental way, through the accumulaton of knowledge, expressed as C = b C W W, (3) 4 whch nvolves mplctly startng values C 0, and W 0,, and s equvalent to the prevous equaton. Technology categores for learnng and those represented n a model may not necessarly concde, and thus a certan amount of mxng, or knowledge spllover, may have to be ncluded. In other words, partcular sets of categores may be closely related technologcally and a learnng spllover matrx B j should be defned n order to calculate an effectve W from ncremental postve capacty addtons: W (t)= B j j ( du j (τ) dτ +δ j U j (τ) ) dτ, δ 0 ju j (τ)dτ, t 0 t du j (τ) dτ > 0 du j (τ) dτ 0, where U j s the capacty of technology j (defned below n eq. 6). Ths equaton thus nsures that knowledge s shared between related technologes. 2 There exsts an extensve lterature on learnng and progress rates for all sorts of goods beyond the electrcty sector (see for nstance Koehler et al. (2006), Pan and Koehler (2007), Grubler et al. (1999)). Wthn the power sector, learnng rates have been compled n both IEA (2000) and McDonald and Schrattenholzer (2001). The consequences of ncludng experence curves are very mportant but lead to strong non-lneartes n any model due to ther strong assocated concavty (Berglund and Soderholm, 2006), and to systems wth multple optmum values or none altogether (Grubb et al., 2002), makng optmsaton models dffcult to solve, a problem that does not arse here. They moreover produce a strong path dependence, an effect we consder fundamental to the process of technologcal change. Snce cost reductons depend on the rato of W /W 0, accordng to an nverse power law, new technologes see hgher cost changes than mature ones for dentcal numbers of addtonal sales. New expensve technologes have the potental to become less expensve than mature ones gven that enough unts are produced. Therefore, n a market where nvestors choose accordng to costs, new technologes requre government R&D nvestment and subsdes n order to become compettve, so-called market push polces (Nemet, 2009). However, wth learnng, such technologes may n tme pass a compettveness threshold, after whch an avalanche effect can happen, where learnng results n more unts sold, whch results n more learnng and so on, and produces an nduced technology transton. No sngle optmum soluton exsts snce costs depend on paths, and an unlmted number of paths can be followed, whch depend prmarly on the assortment of possble polcy decsons. These effects are nherently dynamc n nature, and are thus approprately represented by a dynamc set of logstc dfferental equatons connected to experence curves, whch lead naturally to a behavour very close to the observed S -shaped technologcal transtons, as explored for nstance by Grubler et al. (1999) Levelsed cost of electrcty Investors n the power sector face complex decsons nvolvng a very large number of parameters. However, a standard framework exsts whch s used by most of the ndustry worldwde, that of the LCOE. It s not the focus n ths work to revew ths framework, and we refer the reader to the recent report by

5 the IEA (2010a) on the costs of varous electrcty generaton technologes. We thus use the LCOE n the followng form, LCOE (t)= τ T I (t)+om (t)+fc (t)+cc (t) t=0 (1+r) t τ, (4) EP (t) t=0 (1+r) t where T I denotes the specfc technology nvestment cost, OM the operaton and mantenance costs, FC the fuel costs, CC the carbon cost component assocated wth emssons allowances or taxes where applcable, r the technology or regon specfc dscount rate and EP s the energy that the power staton s expected to produce. The LCOE s the cost whch s compared n the shares equaton. Snce all components of the LCOE can be taken as dstrbutons, the LCOE s also tself a dstrbuton. Defnngα as the emsson factors, GHG emssons and cumulatve emssons are E(t) = α G (t), (9) E tot (t) = t 0 α G (t )dt. (10) The carbon cost component of the LCOE s proportonal toα tmes the prce of carbon, where we neglect effcency mprovements wthn technology categores Techncal Constrants 2.4. Complete set of varables The evoluton of the shares S, as gven above, correspond to changes n the composton of the electrcty generaton G, n GWh, G (t) = U (t)cf (t)=u tot S (t)cf (t), (5) U tot = U, where, U s the capacty n unts of GW, CF s the capacty (or load) factor, defned as the average fracton of tme that a unt runs at full output, tmes a constant of converson between G and U of 8760 h/y. The capacty U can therefore be expressed as U (t) = S (t)d(t), (6) CF(t) CF(t) = S (t)cf (t), where D(t) s the electrcty demand at tme t and CF s the weghed average of the capacty factor over the whole electrcty sector. Changes n capacty stem from ndependent changes n three varables, S, D and CF, emphassng three possble processes. Ths can be expressed n dfferental form, du = S CF dd+ D CF ds S D dcf, (7) 2 CF where the frst term expresses changes n capacty due to changes n electrcty demand, the second due to changes n the composton of the electrcty sector, and the thrd due to changes n the effcency at whch the electrcty sector s used. Investment I(t) n new generaton capacty correspond to postve changes n U tmes ther cost C : I (t)= C (t) ( du (t) dt +δ U (t) ), C (t)δ U (t), du (t) dt > 0 du (t) dt < 0 whereδ s the rate of power plant decommssonng equal to the nverse of the lfetmeτ. Fnally, from G (t) and emssons factors obtaned from the IPCC Gudelnes (IPCC, 2006), t s straghtforward to calculate the GHG emssons that stem from electrcty producton. (8) 5 Fgure 2: Sketch of a hypothetcal profle of daly power demand U D (t D ) as a functon of the tme of day t D, expressed n capacty unts (GW), and how t mght be met by varous types of supply. Areas represent amounts of energy. The green area s the supply of varable renewable energy, whle the red rectangle s the base load supply. The dfference between the sum of these and the demand must be met by flexble systems that can vary ther output from zero to a large fracton of the total demand, shown as a blue area. D(t)/T D s the average daly power demand and T D s the length of a day. Double arrows on the rght correspond to the shares of capacty. The shares equaton treats every technology equally. However, real world electrcty grds face complex optmsaton problems due to fluctuatons n tme of both demand and supply. Whle the demand vares accordng to consumer habts, the supply fluctuates accordng to the varyng nature of some natural resources such as wnd and solar rradaton. For grd stablty, demand and supply must be met at every second, and therefore some flexblty must exst to make both meet. Ths s done wth power plants whch have the ablty to vary ther output rapdly, such as combned cycle gas turbnes (CCGT), desel generators or hydroelectrc dams. However, a large fracton of the electrcty generaton comes from power statons whch do not have the ablty to vary ther output, but have constant power producton (.e. base load systems), such as nuclear reactors or coal power statons, or have an uncontrollable varable output, such as solar panels or wnd turbnes. Fgure 2 sketches a typcal daly demand and supply stuaton, where n green s represented the

6 varable renewables contrbuton, whle n red s shown the base load generaton. U D s the power demand durng a day, and therefore energy unts correspond to areas. We observe that the area n blue, whch durng a day may vary between zero and a large fracton of the total output, must be covered by flexble energy systems wth fast slew rates. These flexble systems can, however, also contrbute to base load producton. We defne here very smple rules that approprately captures the relevant dynamcs wthout requrng complex hourly based smulatons of the power system mert order. We frst classfy generaton technology nto three categores, S Base, S Var and S Flex, whch correspond n order to the base load technologes, the varable renewables and the flexble systems. In real natonal power systems, the herarchcal choce of power statons follows the so-called mert order prncple (see for nstance IEA (2010a)), where power unts are nvoked n order of cost as the demand grows, startng wth base load power statons, costly desel unts beng reserved for relatvely rare demand spkes. We requre here, however, to avod smulatng such a system and replace t wth a set of smple rules. Flexble systems are necessary n order to mantan grd stablty, but may n some cases ncur hgher operatng costs than base load systems. They are n general however rewarded by hgher electrcty prces. The total share of flexble systems can be constraned to reman above a certan necessary value that depends on the amount of avalable storage, the sze of the peak load demand and the amount of renewables present n the power system. Wthn ths necessary share, a choce of flexble technologes exst between whch exchanges can occur, wthn a restrcted electrcty market that meets peak load demand. Smlarly, the total share of varable renewable systems can be made to reman below a certan lmt that s gven by the amount of avalable storage, the sze of the peak load demand and the total capacty of flexble systems. Changes n the capacty of flexble systems can allow expansons of renewables and, conversely, reductons n the capacty of varable renewables enable reductons n the capacty of flexble systems. In order to generate these propertes, we defne as U D the heght of the daly peak demand functon U D (t D ) averaged yearly, where t D s the tme of day, relatve to the lowest daly demand value. We assume that the demand D(t)/T D s the average daly power demand at year t. Therefore, f enough storage was avalable to cover any area above D(t)/T D, denoted D, no flexble systems would be needed, but t s unlkely n most regons of the world gven the large scale of ths amount of energy. 3 Note that changes n avalable storage, dstrbuted generaton, demand response effects and smart grds, whch can nfluence the ssues treated n ths secton through effcency savngs, are consdered outsde of the scope of the present work, but ther effect can be ncluded exogenously or dynamcally from another submodel through a changes n the demand D or n the parameters U D, D, U s and E s. We thus assgn to the storage generaton capacty the varable U s (n GW) and to the energy storage the varable E s (n GWh). 6 The constrant s then expressed as three nequaltes: ( D S Flex CF Flex + S Var CF Var CF D + U VarT D E ) s, D D ( UD S Flex S Var S Base + S Var U tot U ) s, U tot ( CF 1 U D + U s 2 U tot U tot ), (11) where D s the yearly average demand, the same value as n eq. 6, and CF Flex (and smlarly for CF Var ) s calculated as an average of the capacty factor over all flexble systems: CF Flex = 1 S Flex =Flex S CF. The frst these nequaltes expresses the constrant on amounts of energy, where flexble systems are requred to produce the blue area n fgure 2, whle the second expresses the requrement n terms of mnmum capacty of flexble systems, whch must cover at least U D plus the complete capacty of varable systems for when the natural resource s not avalable (.e. for nstance when the wnd does not blow). The requrement on shares of flexble capacty gven by the second nequalty s ndcated wth the double arrow on the rght of fgure 2. The thrd nequalty refers to the maxmum allowed capacty for base load and varable systems, whch cannot exceed the area below the lowest daly demand value. 4 These nequaltes can be transformed nto lmts for the value of shares. They can be nterpreted n several ways. For nstance, we may ask what s the mnmum share of capacty allowed of flexble systems S Flex gven the amount of storage U s and the share of capacty of varable systems S Var. But we may also ask what s the maxmum for S Var gven S Flex and the other parameters. Ths s summarsed by the followng defntons for the lmts on share values (denoted wth a hat): Ŝ = [( UD Ŝ =± U tot U s U tot )+S Var S Flex ] + S, [( CF 1 U D + U ] s ) S Base S Var + S, (12) 2 U tot U tot where n the frst equaton the frst term s postve f refers to flexble output and negatve f t refers to varable renewables 5. The most constranng of both lmts s taken, whch can change accordng to the stuaton. The constrants summarsed by eq. 11 and 12 lead to specfc capacty factors for flexble sources, descrbed n Mercure (2011). These constrants must be respected n order to mantan grd stablty, and thus must be appled to the shares equaton (eq. 1). It s done wth the nterpretaton that nvestment nto specfc technologes slows down when nvestors perceve a rsk of seeng ther capacty unused for grd stablty ssues, and restrct ther technology choces. It s thus through the nvestor fear of stranded assets that the lmtng naturally occurs, an ssue of current sgnfcant mportance for renewable systems (see for

7 nstance the report by Poyry (2009)). Thus we can defne a probablty of reducng nvestment g(s, Ŝ ) near the lmt gven by eq. 12, and assocated cumulatve probablty dstrbuton of nvestng gven whether the system has past that lmt or not, G(S, Ŝ ), whch goes to zero past the lmt. Systems can have ether an upper lmt or a lower lmt dependng on ther nature. The central term of the shares equaton can thus be rewrtten as A j F ( C j )G max G mn j A j F j ( C j )G max j G mn, where G max represents the probablty of nvestng gven S and an upper lmt Ŝ for technology, f t has one (.e. f s a varable or base load system), whle G mn represents a lower lmt, f applcable (.e. f s a flexble system). Near an upper lmt, unts of shares cannot flow nto a category but can flow away from t. Conversely, near a lower lmt share unts cannot flow away from a category but can flow n. These constrants, when ntroduced nto the shares equaton, lead naturally to a splttng of the electrcty generaton technology market when the system s close to the lmt, for nstance constraned by the amount of flexblty. In such a splt market, submarkets forms n whch competton occurs at dfferent cost levels Smple numercal examples wth four categores Share S 1 1 a) b) S 1 S 2 S 3 S 4 c) d) 1 Share S 1 are gven wthout share lmts, n a) and b), and two wth some lmtng, n c) and d). In a), most of the market s gven to S 1, the most expensve technology, wth 97% and 1% for the other three. We observe that wth tme, S 1 gradually loses the market whle the three compettors ncrease ther share exponentally, consstent wth the small share lmt of a logstc equaton. However, competton later arses between these three, where only one, the less expensve S 4, wns the whole market. Therefore, wth an unconstraned shares equaton wth constant costs, the technology wth lowest cost always eventually wns. Ths corresponds to an equlbrum state whch occurs at a tme longer than several tmes the natural tme constant assocated wth these technology changes, of order of several decades. In b), a larger share s gven to S 2. We observe that S 2 ntally ncreases ts share despte the fact that two other technologes have lower costs. Ths s consstent wth a populaton growth phenomenon, where a larger populaton, albet havng a lower brth rate, s able to grow faster than a small populaton. Gven enough tme, however, the technology wth lowest cost stll wns n the equlbrum lmt. These examples demonstrate the equvalence between a parwse comparson of technologes performed wth a dfferental equaton and a smultaneous comparson of all technologes by exhbtng true multple nteractons, provded by the coupled dfferental form of the shares set of equatons. We now ntroduce lmtng, frst n c), where upper lmts of 20% and 30% were gven to S 3 and S 4. We observe that these saturate at ther lmt, and the remanng market s taken by S 2, the second most expensve technology. Therefore, wth lmts, t s not the least expensve that wns the whole market, and we obtan a stable heterogenous system. Addng a lower lmt of 40% for S 1, we see that all lmted technologes become locked to ther lmt and that the unlmted technology S 2 becomes stable. Share S 0 Tme Tme Fgure 3: Examples of calculatons usng the shares equaton n dfferent stuatons, where constant costs are assumed, wth C 1 = 1.3C 4, C 2 = 1.1C 4 and C 3 = 1.05C 4, assocated wth S 1, S 2, S 3 and S 4. a) Startng shares of 97% for S 1, 1% for the others, no lmts are mposed. b) Startng shares of 69% for S 1, 29% for S 2 and 1% for the others, no lmts are mposed. c) Shares as n a) usng upper lmts of 20% and 30% for S 3 and S 4. d) Shares as n a) usng upper lmts as n c) wth the addtonal lower lmt of 40% for S 1. We present a bref exploraton of market competton as modeled by the shares equaton, shown n fgure 3. The system s defned wth four technology categores, the shares of whch sum to one, denoted S 1 to S 4, usng dentcal cost dstrbutons profles but wth dfferent medan values, and A j = 1. The constant medan cost values are gven n terms of the lowest value C 4 : C 1 = 1.3C 4, C 2 = 1.1C 4 and C 3 = 1.05C 4. Two stuatons Share S Natural resource use and depleton 3.1. The cost-supply probablty dstrbuton framework Regonal energy landscapes are defned by the local avalablty of natural resources. Local resources can be assessed and ranked n order of cost of extracton. Naturally, a regonal electrcty sector wll develop accordng to these, unless tradng wth neghbourng regons provdes a more cost-effectve soluton. The feasblty of a global transton from fossl fuels to sustanable energy sources strongly depends on assessments of global natural resources. These, as gven n for nstance the World Energy Assessment (UNDP, 2000, 2004) are ncomplete when they do not provde assocated cost dstrbutons. When natural resources occur wth low energy denstes, n some cases, the underlyng assumpton to the dea of harnessng such sources mples unrealstcally large captal nvestments, such as, for nstance, nstallng large numbers of wnd turbnes on many low wnd speed stes. The resultng expected aggregated energy producton may be very large, but n realty no nvestor would be wllng to undertake such an unproftable venture that nvolves large producton costs per unt energy.

8 a) b) Techncal Potental n ( C ) N ( C ) c) C d) C C ( N ) Techncal Potental C ( N ) N N Fgure 4: a) Sketch of a hypothetcal densty n(c) of cost ranked energy or energy flow unts avalable wthn cost range C up to C + C, wth assocated uncertanty ndcated usng error bars. b) Number of energy or energy flow unts N(C) avalable for a cost of extracton below the arbtrary value C. c) The cost C(N) of extractng an addtonal energy or energy flow unt gven that N were already exploted, commonly called the cost-supply curve d) Cost-supply curve defned as a probablty dstrbuton, where the red curves ndcate the 95% confdence level range and the blue curve the most probable cost-supply curve. The framework of cost-supply curves, wdely used after the work of Rogner (1997) and taken agan by Hoogwjk et al. (2004, 2005, 2009), Hoogwjk (2004) for ther assessments of wnd, solar and bomass global energy potentals and used for nstance n the TIMER model (Bouwman et al., 2006) addresses ths queston. We extend here ths approach to nclude the mssng treatment of uncertanty and the avalablty and evoluton of knowledge wth respect to natural resources. We use these defntons as one of the central aspects drvng the shares equaton n FTT:Power, n order to produce meanngful regonal energy landscapes. Ths framework s depcted n fgure 4. 8 We defne the densty n(c) of unts of energy, for nonrenewable sources, or unts of energy flows, for renewable sources, as a functon of ther cost of extracton C (see fgure 4 panel a). Ths corresponds to a hstogram of the number of energy unts that can be extracted for a cost wthn a certan range C. Data ponts n ths densty functon possess an uncertanty value (schematcally depcted usng error bars), correspondng to the lack of exact knowledge over the number of energy unts that can be extracted wthn a partcular cost range. The number of unts that can be extracted for a cost lower than the arbtrary value C s the ntegral of the densty n(c) up to C, that we denote N(C ) (panel b). Ths functon converges towards a unque value at hgh C, whch corresponds to the total techncal potental, and s the total area under the densty n(c). Note that the assocated uncertanty grows cumulatvely. Ths relatonshp can be nverted n order to express the cost of extractng an addtonal unt of energy or energy flow after N have already been exploted, denoted C(N) (panel c). Ths cost-supply relatonshp dverges at the total techncal potental. The dvergence corresponds to, for example, the nstallaton of a dvergng number of wnd turbnes on a very large number of stes wth vanshngly small average wnd speeds, but nevertheless gvng a fnte amount of energy flows. Usng uncertanty values, one can n prncple determne three cost-supply curves n order to defne a confdence range for where the real actual costsupply curve may le. Ths replaces the tradtonal cost-supply curve by a probablty dstrbuton that reflects mperfect knowledge. Panel d depcts ths, where the red curves determne the 95% confdence level regon n cost-quantty space, and the blue curve corresponds to the most probable set of values. Uncertanty n the determnaton of natural resource avalablty s notable n the case of fossl fuel reserves and resources. Rogner (1997) pants a clear pcture of the process of fossl reserve expanson. In ths vew, dscoveres of deposts of hydrocarbons expand known resources, but ther costs of extracton are hghly uncertan. 7 However, wth technologcal learnng and addtonal exploraton, knowledge over cost values mproves but costs also gradually decrease, resultng n a flow from uncertan and costly resources to known and economc reserves. Thus, as tme progresses nto the future, the cost-supply curves of fossl resources become better and better defned, an effect that can strctly only be descrbed by a probablty dstrbuton defned now based on current knowledge, whch determnes the set of probable future paths.

9 3.2. Cost-supply curves n FTT:Power Cost-supply curves are necessary to defne the level of use of each natural resource n a model for electrcty producton such as FTT:Power. The underlyng assumpton s that wthn a specfc world regon, unts of a partcular resource wll be used n order of cost to the best of current knowledge. In the case of wnd power for nstance, wnd stes wth partcularly good wnd speed dstrbutons and low turbulence propertes are lkely to be developed frst, furthermore n areas where they are lkely to be socally accepted and where land s affordable. All these aspects should therefore be ncluded n such assessments, an effort whch has been done by Hoogwjk et al. n an extensve assessments of global wnd potentals (Hoogwjk et al., 2004). However, the level up to whch resources are used depends entrely on the avalablty and cost of all alternatves whch can be used to produce electrcty n that regon n order to meet the demand. Ths s precsely what s done n FTT:Power, where cost supply curves are used for every calculaton of the LCOE as defned by eq. 4. Ths s an aspect of central mportance, snce the complete set of cost-supply curves mantan the calculaton wthn boundares defned accordng to a complete set of total techncal potentals and thus yelds reasonable results of natural resource use. Addtonally, a margnal cost of electrcty producton can smply be derved usng the LCOE averaged over all technologes, S LCOE (S beng the shares of capacty), and used wth addtonal nformaton regardng taxes and proft margns to determne the prce of electrcty. The prce of energy carrers beng mportant for the global economy, as represented for nstance n E3MG, ths aspect of the model has a deep nfluence over future projectons of energy demand and supply and assocated effects on economc actvty. Moreover, cost-supply curves nherently track closely resource depleton. The set of cost-supply curves n a specfc world regon defne a margnal cost of electrcty producton, and therefore results n all world regons havng dfferent prces of electrcty. In a real world, ths nevtably results n electrcty tradng between regons whch have ther grds connected, from resource poor to resource rch regons. Examples of ths are the sale of mportant amounts of Canadan or Scandnavan hydroelectrcty n Amercan or European markets respectvely. Thus except for closed energy markets, the level of use of natural resources for electrcty producton determned through cost-supply curves of world regons should therefore addtonally depend on each other. The detals and effects of energy tradng between the 20 E3MG regons wthn FTT:Power are complex and wll be explored elsewhere. Moreover, wth the ntroducton of other sectoral models of technology of the FTT famly nto E3MG, competton for resources wll arse, an aspect that wll be dscussed n later work. FTT:Power uses determnstc cost-supply curves, or n other words, a sngle curve amongst all possble paths defned by the probablty dstrbuton of fgure 4 d). Although the decson makng underlyng the shares equaton (eq. 1) should nclude an assessment of rsk and uncertanty n nvestment, a subject whch s not the scope of the present work but wll be the theme 9 of a subsequent paper, the uncertanty over natural resources does not nfluence the choces of nvestors drectly. Ths s due to the fact that the cost-supply curve always has a vanshng uncertanty at the current level of natural resource use, exemplfed by a well defned prce. It s the level of future use whch s uncertan, generally beyond the typcal tme horzon for capacty nvestment. Thus, at every tme step of the smulaton, resource avalablty and costs as seen by nvestors should be well defned quanttes, wthout uncertanty. It s the number of possble futures whch s large. The determnaton of uncertanty over predctons therefore requres runs of the model usng large numbers of possble determnstc cost-supply curves derved from probablty dstrbutons akn to that n fgure 4 d), for all sectors. Monte-Carlo technques are thus the approprate way forward n order to defne confdence levels attrbuted to future predctons, an aspect we treat as essental n the formulaton of ths model, and wll be explored n detal elsewhere. 4. Results 4.1. Model parameters and assumptons The mplementaton of FTT:Power whch was used to produce the results presented here was made wth the defnton of 24 energy technologes and one sngle world regon, along wth global cost-supply curves for 13 types of resources, some of whch are used by several competng energy technologes. A table summarsng these s provded n the appendx. As a standalone model, t essentally requres a curve for energy demand and startng values for the capactes U, and calculates energy capacty and generaton from eq. 5 and 6, wth eq. 1 drvng the changes n the energy mx. For the sake of explorng the propertes of FTT:Power and presentng results whch stem from the smplest assumptons, we have excluded any feedback wth the global economy, from electrcty demand and through electrcty prces, by runnng the model by tself wthout E3MG. Therefore, we used an exogenous electrcty demand curve. Components of the LCOE, gven n the appendx, were derved from both recent IEA data (IEA, 2010a) and from a set of global cost-supply curves (a detaled analyss s forthcomng, n Mercure and Salas, 2012). These curves were constructed from an ensemble of sources n order to produce a complete set (IAEA, 2009, Rogner, 1997, Hoogwjk et al., 2004, 2005, 2009, Hoogwjk, 2004, Themels and Ulloa, 2007, Mock et al., 1997, IJHD, 2011, Lako et al., 2003). In the case of Bomass, the global potental depends hghly onto assumptons regardng world populaton and global food demand, and the results of Hoogwjk et al. (2009) for the B1 SRES scenaro were used. For ths paper, as a sngle regon global model, the assumpton s made of perfect electrcty tradng worldwde, such that every type of resource s avalable everywhere, an assumpton whch s done here for the convenence of usng a sngle world regon for the purpose of the demonstraton of the valdty of the model equatons, but whch wll be dropped after the subsequent defnton of cost-supply curves for 20 world regons that concde wth those of E3MG, an aspect of FTT:Power whch s not the focus of the present work. The results that we present,

10 Nuclear Ol Coal Coal+CCS IGCC IGCC+CCS CCGT CCGT+CCS Sold Bomass S Bomass+CCS BIGCC BIGCC+CCS Bogas Bogas+CCS S Hydro L Hydro Onshore Offshore Solar PV CSP Geothermal Wave Fuel Cells CHP Capacty ( TW ) Electrcty Generaton ( PWh ) Capacty ( TW ) Electrcty Generaton ( PWh ) Nuclear Ol Coal Coal+CCS Gas Gas+CCS Bomass Bomass+CCS Hydro Wnd Solar Geothermal Wave Fgure 5: Model results for two sets of assumptons, a baselne (left panels) and a mtgaton (rght panels) scenaros. The top panels show the capactes U for 24 energy technologes, whle the bottom two show electrcty generaton G as areas that sum up to the total electrcty demand. For G, related categores were aggregated nto a smaller set for clarty. 10 although hghly aggregated, are of nterest n themselves and n how they depct the technology substtuton processes through whch a transton towards a low carbon electrcty sector may occur n the future. We present here results for two smple scenaros, a baselne where no mtgaton effort s made but CO 2 prcng exsts, and a scenaro where emssons are decreased through reducng the number of allowances such that the prce for allowances ncreases exponentally (n real 2008 dollars). We thus assumed a constant prce of carbon of 22$/t CO 2 n the baselne scenaro, whle n the mtgaton scenaro the prce of carbon s the same ntally but ncreases by 1% per year. Ths was done n order to have an dentcal startng pont, such that a dvergng evoluton of the nature of new energy capacty can be observed, entrely attrbutable to the dfference n the prce of carbon. 8 Wth FTT:Power as part E3MG, the prce of carbon wll be endogenous. The prce of carbon s fed back nto the system through the carbon component of the LCOE (eq. 4), where a hgh dscount rate of 10% was used. 9 Therefore, an economc ncentve to depart from heavly emttng technologes and edge towards more effcent systems s created. For smplcty, however, no feed-n tarffs or subsdes were assumed for any of the renewable categores, and therefore wnd and solar systems do not appear n the results. Both scenaros use an dentcal demand functon, whch was derved from that of the New Polces Scenaro of the World Energy Outlook 2010 (IEA, 2010b). It was extended beyond 2035 assumng a slghtly decreasng growth rate, chosen as a good example rather than for accuracy. In the real world, the demand s expected to be reduced n mtgaton scenaros where ncreases n producton costs due to an ncreasng prce of carbon or technology swtchng s lkely to be passed on to the consumer through ncreasng electrcty prces. Such an effect would consttute a drect result of feedback between FTT:Power and E3MG Model Outputs Fgure 5 presents the results for the capacty U (t) (top panels) and the power generaton G (t) (bottom panels), the latter beng represented as coloured areas whch sum up to the total electrcty demand. The baselne s shown n the left panels, whle the mtgaton scenaro s shown n the rght panels. The baselne scenaro s domnated by coal fueled power statons, except for a szable amount of hydroelectrcty, whch s seen to converge towards a capacty of about 0.7 TW worldwde. Ths s due to the fact that a large fracton of the world s total hydroelectrc potental s already developed, and that costs rse very rapdly wth addtonal dams beng bult. Thus, ths partcular sngle world regon calculaton predcts a scarcty of sutable economc stes for addtonal hydroelectrc projects. 10 The remanng demand s covered mostly by coal power statons and combned cycle gas turbnes (CCGT). The mtgaton scenaro exhbts more rchness wth a larger varety of substtuton events trggered by the ncreasng prce of carbon. Hydroelectrcty appears hardly affected by the dfference n assumptons. However, as the prce of carbon ncreases, coal s gradually phased out, peakng n around 2040, where CCGT takes the lead. Ths s due to the much ncreased effcency of CCGT, whch emts less GHGs per unt of energy compared to conventonal coal. Smple CCGT systems are however n turn gradually phased out from 2070 onwards