Next Article in Journal
Analyzing the Impacts of Financial Services Regulation to Make the Case That Buy-Now-Pay-Later Regulation Is Failing
Previous Article in Journal
Effect of Soil Management on Erosion in Mountain Vineyards (N-W Italy)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Annualization of Renewable Investment Costs for Finite Horizon Electricity Pricing and Cost Recovery

1
Institute for Research in Technology, Technical School of Engineering, Comillas Pontifical University, 28015 Madrid, Spain
2
Institute for Systems and Computer Engineering, Technology and Science, 465 Porto, Portugal
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(4), 1993; https://doi.org/10.3390/su13041993
Submission received: 17 December 2020 / Revised: 3 February 2021 / Accepted: 5 February 2021 / Published: 12 February 2021

Abstract

:
The increasing penetration of renewable electricity generation is complicating the bidding and estimating processes of electricity prices, partly due to the shift of the overall cost sensitivity from operation (fuel) costs to investment costs. However, cost minimization models for capacity expansion are frequently based on the principle that, for a perfectly adapted system allowing non-served energy, marginal remuneration allows overall operation and investments costs recovery. In addition, these models are usually formulated as finite-horizon problems when they should be theoretically solved for infinite horizons under the assumption of companies’ infinite lifespan, but infinite horizon cannot be dealt with mathematical programming since it requires finite sets. Previous approaches have tried to overcome this drawback with finite horizon models that tend asymptotically to the original infinite ones and, in many cases, the investment costs are annualized based on the plants’ lifespan, sometimes including a cost residual value. This paper proposes a novel approach with a finite horizon that guarantees the investment costs’ recovery. It is also able to obtain the marginal electricity costs of the original infinite horizon model, without the need for residual values or non-served energy. This new approach is especially suited for long-term electricity pricing with investments in renewable assets when non-served demand is banned or when no explicit capacity remuneration mechanisms are considered.
JEL Classification:
G12; O21; Q55; L94; C61

1. Introduction

1.1. Problem Description

In 2030, half of the EU’s electricity demand will most likely be met with renewable energy sources, and by 2050 the power sector will have to be almost fully renewable. Classical marginalist theory [1] states that assuming: (1) marginal cost pricing based on short-term marginal production costs; (2) a perfectly adapted system; (3) a set of technologies available for investments; and (4) a cost for the non-supplied energy, the overall cost (investment plus operational costs) is recovered for all technologies. However, this theory requires the existence of unsupplied demand at certain peak hours with very high prices, in order to allow marginal plants to recover their overall costs, which is unrealistic from a practical point of view in many power systems in developed countries. Moreover, current market mechanisms based on the marginalist theory could not provide adequate economic signals for investment, especially in systems with a high share of renewables. Therefore, remuneration mechanisms for generation, based only on short-term variable costs could imperil the 2050 EU objectives. Indeed, in practice, null variable costs of renewables lead to close-to-zero prices, making it more difficult to recover their full costs, thus compromising not only their profitability but also the economic viability of other backup technologies needed to meet demand in the absence of wind and sun and to provide regulation reserves. Additionally, the integration of renewable installations might lead to systems that suffer from an effect of profitability-loss cannibalization, since the higher the renewable production, the lower the market incomes received by all generation assets [2].
In this context, well-designed capacity remuneration mechanisms may overcome this problem, as capacity pricing can establish the amount of firm capacity that can be provided by each technology, being an essential tool for cost-effective capacity investments. The absence of these mechanisms may discourage capacity providers from investing in new generation capacity, putting at risk supply security.
This paper analyzes the problem of investment cost recovery in a power system with a large penetration of renewables without the existence of an explicit capacity remuneration mechanism. This is a critical problem for future fully renewable systems even with a capacity remuneration mechanism since, as discussed in [3], renewable generation has so far had very limited participation in capacity mechanisms.

1.2. State of the Art

Generation expansion studies used to analyze the problem of cost recovery in the absence of explicit capacity mechanisms are frequently based on the formulation of finite-horizon mathematical programming (MP) problems, even knowing that finite horizons are not suitable in the long-term planning with infinite-life generation companies, [4]. Two main approaches to overcome this problem have been proposed in the literature on this topic, the static and the dynamic approaches [5].

1.2.1. The Static Approach

This approach determines the optimal investment plant (OIP) for a single representative year in the future, dividing the investment costs by the lifespan of the generation assets. Classic works that use this approach are [6,7,8,9,10]. For example, static MP models including the network constraints are presented in [9,10]. A more recent static cost minimization MP model including the network is described in [11], while [12,13] present static MP models for the analysis of the operating reserves and the carbon dioxide emissions, respectively. While it is useful in some situations for simplified studies, this single-year static approach may be misleading in the current context of the rapid development of renewable energy technologies, as their investment costs decrease on a yearly basis. Static analyses provide very vague results regarding the year in which investments should be made since they do not provide insights on how electricity prices could evolve annually.

1.2.2. The Dynamic Approach

Unlike the static approach, the dynamic approach computes the OIP for a multi-year finite horizon. A residual value modeling the costs for the years beyond the finite horizon is usually included indefinitely replicating the cost of the last year of the horizon [14]. A residual value is alternatively considered in [15], which uses the resale value of the generation assets at the end of the study year. However, residual values may be volatile and complex to compute when considering a company with a big market share. Moreover, the replication may lead to inaccurate results, for example, if investment costs decrease significantly over time, as commonly happens with renewable technologies. Some other works such as [16,17,18,19], do not represent residual values, assuming that they can be disregarded under some circumstances such as, for example, for long finite horizons and/or high discount rate values.
Most of the above works distribute overall costs over the years by applying a particular annualization of these costs and allocating them based on the generation assets’ lifespan. The most classical annualization of costs is based on the capital recovery factor. For example, in [20] a capital recovery annualization is applied to reduce renewable energy curtailment by improving the accommodation of this generation taking into account short-term operation constraints and uncertainties from the load and renewable energy sources.

1.3. Contributions of the Proposed Dynamic Approach

This paper overcomes some of the limitations of the abovementioned approaches by analyzing the problem of determining the long-term marginal costs (hereinafter termed electricity prices, for simplicity’s sake) in order to guarantee the recovering of the investment costs of renewable technologies. It applies a dynamic approach, more suitable than the described static ones, in the absence of an explicit capacity mechanism, and when non-served demand is not allowed. Unlike the dynamic approaches presented in the literature, this research:
  • Focuses on the long-term recovering electricity prices (understood as those that, under some assumptions, guarantee the recovery of the total costs of producing electricity), computed by the dual variables of the generation-demand balance constraints of the MP problems, instead of aiming at the primal variables to estimate the OIP as other approaches do (for instance [16,17,18,19]).
  • Assuming infinite-life generation companies, proposes a finite MP problem with a novel annualization of the investment costs which ensures that the electricity prices obtained from the finite problem are the same as those obtained using the corresponding infinite model. As far as the authors know, previous works (for instance [15]) have not been able to ensure the fulfilment of this interesting and relevant property. Hereinafter, when this happens, it will be said that the annualization method verifies the finite dual pricing (FDP) property.
    More specifically, if C is the overall cost of the model for an OIP under an infinite horizon, this paper formulates a finite truncation C of C using only p years such that the derivatives of C and C with respect to the demand dp at any year pp (the dual variables of the generation-demand balance constraint) coincide. This paper also proves the well-known result that states that the straightforward truncation Cp of C (i.e., solving the problem for a finite number of years without annualization) does not verify the FDP property and leads to large runtimes when demanding a reasonable approximation of C by Cp.
  • Formulates, in order to deduce the annualization method, a stylized cost-minimization MP problem with a finite horizon without the need of considering the abovementioned residual value that current dynamic approaches need (for instance [14,15]), which makes it easier both from a theoretical and practical point of view.
  • Significantly reduces the estimation time needed to reach reasonable approximations of the electricity prices relative to other dynamic approaches based on finite MP problems (for instance [21]), due to the contribution referred in the previous points.
  • Even if the annualization proposed guarantees the total cost recovery only under some assumptions (see Section 2), in view of the lack of methods that guarantee the FDP property, reasonable price estimations could be obtained even under more relaxed hypothesis. However, additional research is needed in this point.
In short, this paper presents a relevant theoretical advance that allows simplifying an infinite OIP problem into a finite one to estimate long term electricity prices. This approach can be directly used from a practical point of view to estimate the total electricity price that guarantees investments’ recovery for the economic sustainability of the power sector. This result is useful both for the investment planning of the electricity companies and for the regulators assessment of current markets designs.
This paper is structured as follows. The main assumptions that define the presented framework are described in the following section. Section 3 is devoted to analyzing the OIP and its sensitivity with respect to the demand. From this result, the marginal prices that guarantee the investment costs recovery are computed in Section 4, including and excluding corporate taxes. Section 5 presents the annualized investment costs proposed in this paper and proves how it fulfils the FDP property. Conclusions and future developments are discussed in the final section.

2. Main Assumptions

The main assumptions of this paper are: (1) for simplicity’s sake, investments are made in renewable generation with null operational costs; (2) no capacity remuneration mechanisms are explicitly modeled nor needed; (3) no unsupplied demand is allowed, as occurs in many real systems; and (4) other remuneration mechanisms might not exist.
The following modeling hypotheses have been also assumed (parameters are written in lowercase, decision variables in uppercase, and sets with capital Greek letters):
  • There is a centralized decision-maker with infinite life that minimizes the total investment costs of generation assets with null operation costs.
  • A finite temporal horizon Γ⊆   , which is discretized into periods p whose duration is typically one year for long-term planning is considered. Periods p∈ℕ are those where investments Kp (in MW) must be decided, while Kp for p∈Γ\ℕ are investments that have already been decided (and, therefore, they are inputs, i.e., Kp = kp).
  • The lifespan of the investment Kp at period p is l > 1, i.e., the lifespan extends from p to p + l − 1.
  • The rate rp > 0 to discount a future amount of money to reflect its present value does not depend on p, i.e., rp = r.
  • Only one technology available for investment and, as other authors assume (see for example [22]), its investment cost fp (€/MW/period) increases according to inflation i:
    f p + p   =   f p ( 1 + i ) p   =   f 0 ( 1 + i ) p + p
  • The utilization ep of capacity Kp does not depend on p, i.e., ep = e. Parameter e is well-known as the utilization of the asset Kp, i.e., the ratio of the energy produced in p divided by the installed power. Therefore, a mature technology is assumed.
  • Inflation i is less than r and, therefore, the present value of the investment cost fp is decreasing:
    f p + 1 ( 1 + r )   =   ( 1 + i 1 + r ) f p   <   f p
    Hereinafter the quotient (1 + i)/(1 + r) will be called real discount factor (RDF).
  • Demand dp (MWh demanded in p) is inelastic and must always be satisfied. It is represented as a single block per year that, in capacity expansion analyses, may be inputted as the yearly peak value multiplied by the number of hours in p. Note that dp must be satisfied with the capacities Kp at periods p′ ∈ Γ with pl + 1 ≤ p′ ≤ p.
Although the hypotheses described here may appear simplistic (for example, the consideration of only one technology for investment), from a practical point of view, the annualization of the resulting investment costs, based on a thorough and theoretical rigorous mathematical analysis, must be understood as a mathematical financial tool for investment decisions that may be applied in more complex analyses (like those with more than one technology for investment). Indeed, several annualizations may be chosen but the proposed one guarantees the FDP property under the abovementioned hypotheses, which could be relaxed in the absence of more rigorous studies. Understood in this way, the annualization we propose is a powerful mathematical tool to estimate long term electricity prices, assuming that, whatever the market mechanism is, total electricity prices should guarantee the investment costs recovery, ensuring the economic sustainability of the power sector even when non-served demand is not allowed. Moreover, the proposed annualization can help not only in the investment planning of electricity companies, but also in the reforming processes of current electricity markets. These processes must be set by regulators, which need tools to assess, for instance, the need of capacity remuneration mechanisms in power systems with a large amount of renewable generation, in which energy only markets tend to fail due to the existence of very low variable costs.

3. Optimal Investment Plan (OIP)

Considering the above-mentioned assumptions, this section describes some useful properties of the OIP and how to obtain it.

3.1. OIP Computation

Appendix A proves that the OIP is unique and can be directly computed as follows:
K p   =   k p m p l + d p ( m p 1 ) l d p ( m p 1 ) l 1 e + + d p l d p l 1 e + d p d p 1 e
where:
m p   =   min { m / p m l   <   0 }
Therefore, the new capacity investments Kp at period p depend only on inputs of the model, namely, the demand in periods before or equal to p, the hours of utilization e, and the capacity kp already installed for the period p′∈Γ\ℕ closest to the first period (p = 0) of the temporal horizon. Note that each quotient in (3) represents the capacity needed to satisfy the variation of the demand in two consecutive years p′ and p′ − 1, p′ being a year before p a number of times multiple of the lifespan l, when certain quantity of the installed capacity expires.

3.2. OIP Sensitivity with Respect to the Demand

Appendix B proves that the variation ΔKp of the OIP with respect to the variation ΔDp in the demand dp at period p is:
Δ K p Δ D p   =   { 1 e   i f   p   =   p + n l   1 e   i f   p   =   p + n l + 1 0         o t h e r w i s e   for   n
This result will be used in Section 5 to analyze how ΔDp influences the total cost of the investments and means that the variation ΔKp in a future period p′ respect to the variation ΔDp in a past period p takes only values if the difference between both periods p′ and p is a multiple of the lifespan l or a multiple of l plus 1 (since, as mentioned before, the capacity Kp must satisfy the variation of the demand in two previous consecutive years p and p − 1). Moreover, when positive, the variation ΔDp relative to ΔKp is the utilization e or −e since dpdp1 is not satisfied with the OIP but with the production generated from the OIP.

4. Prices to Recover Investment Costs

Appendix C proves that, for the previously computed OIP, the energy prices Pp (€/MWh) that allow the exact recovery of the investment costs fp satisfy:
P p   =   f p e s   =   P 0 ( 1 + i ) p
where:
s   =   p   =   0 l 1 ( 1 + i 1 + r ) p
Therefore, on the one hand, price Pp is equal to fp divided by the utilization e of the investment, taking into account s, which is the summation of RDF over the generation assets’ lifespan. On the other hand, Equation (6) identifies the well-known result that states that the prices to recover investment costs growing up with the inflation i should also increase with the inflation i [23]. Taxes can be taken into account in (6) by substituting fp with fp as shown in Appendix D.

5. Finite Dual Pricing Property (FDP)

This section presents the annualized investment costs proposed in this paper and proves how the resulting electricity marginal prices fulfils the FDP property. It also presents how these prices can be computed with MP, avoiding complexities like dealing with infinite horizons.

5.1. The Infinite Horizon Problem with Cost C: Prices to Recover Investment Costs as Dual Variables

If the original infinite horizon problem is considered, the minimization of the present value of the total investment cost is:
M i n K p     0 { C   =   p   =   0 f p K p ( 1 + r ) p }
subject to the generation-demand balance constrain (dp must be satisfied with a utilization e of the non-expired investments Kp installed at periods before p, i.e., pl + 1 ≤ p′ ≤ p):
d p   =   e p   =   p l + 1 p K p
Appendix E proves that the dual variables Lp of (9) when minimizing the objective function in (8) are:
L p   =   { n { : p + n l } f 0 e ( 1 + i 1 + r ) p + n l n { : p + n l + 1 } f 0 e ( 1 + i 1 + r ) p + n l + 1 }
It also proves that Lp are the discounted prices Pp of (6) (taxes can be considered as indicated in Appendix D), i.e.,:
L p   =   P p ( 1 + r ) p
This result can be also drawn from [23]. However, C has an infinite summation that makes impossible the exact computation of primal investment decisions with MP. To overcome this drawback, the next subsection presents a finite truncation C of C that guarantees the FDP property without a residual value in its formulation.

5.2. The Straightforward Truncation Cp of C

It is well known that the straightforward truncation Cp of C up to p periods (C with ∞ = p in (8)), without including a residual value in Cp, is not valid for computing the discounted prices Pp as the dual variables Lpp of (9) minimizing Cp (note that Cp and Lpp are formulated as C and Lp substituting ∞ by p in (8) and (10), respectively). As an example, let us consider the data in Table 1.
Apart from the chosen lifespan l = 5 which has been kept small for the sake of simplicity (even though this will yield too high energy prices in this example), the rest of the data may correspond to onshore wind farm installations (according to the International Renewable Energy Agency, see [24]). With these data, (7) leads to s = 4.3237.
Using this dataset, Figure 1 depicts Lp7 (green discontinuous line named by LS7) and the discounted price P7/(1 + r)7 = 87.03 €/MWh (brown line named by P7), as a function of p. Each grey column in Figure 1 represents one of the two summations in (10) (making ∞ = p), positive columns (grey columns named by P) corresponding to the first summation (adding n ∈ ℕ such that p+n⋅5 ≤ p, with p = 7) while negative ones (grey columns named by N) to the second (n ∈ ℕ: 7 + n⋅5 + 1 ≤ p). According to (10), Lp7 is the summation of all the columns (whether positive or negative ones) before p.
Note that for any p, Lp7 does not exactly coincide with P7/(1 + r)7. However, a good approximation can be computed with a value of p sufficiently large. Figure 1 also presents the approximation error of Lp7 with respect to P7/(1 + r)7 (in %, blue line labelled “Error”). As can be seen, to estimate P7 with less than 10% error, 30 additional years are needed, leading to a total horizon length of 7 + 30 years. According to an investment model based on [21,25] with only one week per year, an hourly basis, and a time horizon of 37 years, the execution time for such an error would be approximately 45 min. Larger horizons lead easily to undesirable execution times or “out of memory” errors on a 64-bit Intel-Core CPU at 3.4 GHz. However, as the next subsection proves, the proposed truncation C of C provides the exact price P7 with a total horizon of 7+0 years in only 5 min, considering the abovementioned computer.

5.3. The valid truncation C of C for FDP

Let us consider an annualization fp,p of fp in Γp= {p, …, min(p, p + l − 1)}, p ∈ {0, 1, …, p} and p′ ∈ Γp (thus, annualization over the generation assets’ lifespan), i.e.,:
f p   =   p Γ p f p , p ( 1 + r ) p p
and such that fp,p is uniform and proportional respect to fp, i.e.,:
f p , p   =   f p s ( 1 + i ) p p p { 0 , , p _ } , p Γ p
The objective function (8) truncated to {0, 1, …, p} using fp,p becomes (taxes can be easily included by substituting fp by fp, as it is described in Appendix D):
C   =   p   =   0 p _ { K p p Γ p f p , p ( 1 + r ) p }
Appendix F proves that the dual variables Lp of (9) when minimizing the objective function (14) are:
L p   =   { f 0 s e u   =   p p _ ( 1 + i 1 + r ) u f 0 s e u   =   p + 1 p _ ( 1 + i 1 + r ) u }
Using the same inputs of Table 1, in Figure 2, grey columns represent each summand in (15), positive columns (grey columns named P) corresponding to the first summation while negatives (grey columns named N) to the second one. It also represents L7 (green discontinuous line named L7) and the discounted price P7/(1 + r)7 (brown line named P7) of (6), both as a function of p.
As can be seen only the first term in the first summand of (15) is not compensated with a term in the second summand. Therefore:
L p   =   f 0 s e ( 1 + i 1 + r ) p   =   f p s e ( 1 + r ) p   =   P p ( 1 + r ) p
Thus, L p are the discounted prices Pp of (6), and no approximation is needed, i.e.:
L p p _ L p   =   L p   =   P p ( 1 + r ) p
For generalization purposes, proof of Appendix F considers the next generalization f p , p a of fp,p (generalization if a = i) as an alternative non-uniform annualization of fp in Γp for FDP, which depends on a parameter a:
f p , p a   =   f p s p a s ( 1 + a ) p p   p { 0 , , p _ } , p Γ p
being:
s p a   =   p Γ p ( 1 + i 1 + r ) p p p Γ p ( 1 + a 1 + r ) p p

6. Conclusions

The current context of decarbonization of the energy system, with the closing of fossil fuel plants and the electrification of energy end usages, is entailing a significant increase in generation capacity based on renewable sources of energy, most of them with negligible operation costs. In this context, there is a need to assess how the electricity prices that would allow the recovery of the total generation cost (operation and investment costs) could evolve in the future years. This paper addresses this challenge when no scarcity prices (that appear in hours with non-supplied demand) are allowed, which is the practical approach of regulators and system operators of many power systems.
Based on a stylized cost-minimization model that only considers one renewable technology for investment, this paper proposes a novel annualization of the investment costs (summarized in Equations (12) and (13)) that ensures the recovery of the investment costs. Although mathematically based on a marginalist approach with no additional capacity mechanisms, the proposed model is intended to provide the total electricity price estimation independently of the market mechanism in place. The model can be solved with a finite MP problem despite the hypothesis of infinite life for companies, and the electricity prices are obtained from the dual variables of the demand balance constraints.
With regard to other similar approaches in the literature, this approach:
  • Focuses on the electricity prices instead of on the estimation of the primal investment variables, proving they are the same as those that had been obtained from the original infinite model. In addition, no residual values to represent costs in periods beyond the finite horizon are needed.
  • Allows computing the exact m marginal prices of the infinite model using the proposed model with a finite horizon of only m periods. The straightforward truncation of the infinite problems would need a much longer time horizon to reach a reasonable approximation of these same prices.
Although some of the hypotheses of the model might seem simplistic (especially the assumption of a single renewable technology for investment), the proposed annualization of the investment costs must be understood as a sound financial mechanism to estimate long-term future electricity prices according to the awaited expansion of the renewable generation capacity. These types of studies are especially essential for generation companies when planning future investment strategies, but also for regulators to design suitable subsidies, capacity mechanisms, and energy policies, especially under the expected penetration of renewable plants.
The proposed annualization is a first step for long-term prices estimation to mitigate the approximation issues of other similar approaches that truncate the original infinite horizon. Current authors’ work aims to extend this approach to represent several simultaneous technologies and generation companies, include storage facilities, and include non-constant utilization factors and inflation rates.

Author Contributions

Conceptualization, F.A.C.; Formal analysis, F.A.C.; Funding acquisition, J.V.; Investigation, F.A.C. and J.V.; Methodology, F.A.C. and J.V.; Project administration, F.A.C. and J.V.; Supervision, F.A.C. and E.C.; Writing—original draft, F.A.C.; Writing—review & editing, F.A.C., J.V. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially financed by National Funds through the Portuguese funding agency, FCT—Fundação para a Ciência e a Tecnologia—within project UIDB/50014/2020.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Indexes and sets
,   integer/natural numbers
Γfinite set of periods
p, pperiods and last period in G
Γπset of periods with capacity decided at p
Scalars
iinflation (%)
llife span of investments (periods)
Parameters
fpinvestment cost (€/MW/period)
taxcorporate taxes (%)
f p investment cost including taxes (€/MW/period)
dpdemand (MWh)
rpdiscount rate (%)
ephours of utilization (hours)
kpcapacity investments already decided (MW)
fp,p′annualized investment costs (€/MW/period)
f p , p a annualized investment costs with non-uniform annualization (€/MW/period)
apamortization of investment costs (%)
Variables (capital letters)
C¥infinite horizon cost (€)
Cpstraightforward truncation of C¥ up to p periods (€)
Cproposed truncation of C¥ up to p periods (€)
L¥p,Lpp,Lpbalance constraint dual variable minimizing C¥, Cp, and C (€/MWh)
Kpnew capacity investments (MW)
ΔΔπdemand variation between p and p + 1 (MWh)
Ppenergy prices (€/MWh)

Appendix A. OIP Computation

This subsection proves that, with the above assumptions, the OIP is unique and can be computed from inputs by directly solving algebraic equations. Since dp must be satisfied with a utilization e of the non-expired investments Kp installed at periods p′ previous to p, i.e., pl + 1 ≤ p′ ≤ p, it is:
d p   =   e p   =   p l + 1 p K p
Applying (A1) to p − 1 and subtracting it from (A1):
d p d p 1   =   e ( K p K p l )
And therefore:
K p   =   K p l + d p d p 1 e
Kp is thus given by the reposition of the capacity Kp−l just expired (no longer available at period p) plus the amount needed to satisfy the demand variation dpdp−1 (scaled by the utilization e). Since we are assuming that Kp ≥ 0, clearing from (A3), the demand dp must satisfy:
( d p d p 1 ) e K p l
Therefore, we add the hypothesis that −(dp − dp−1) must not be larger than the used capacity e⋅Kp−l just expired, even if the demand dp between two consecutive periods decreases (i.e., dpdp−1).
Applying (A5) to p = 0, it can be seen that K0 depends on d0 and d1, but also on kl which is the capacity already installed l periods before the first one p = 0. Therefore, K0 depends only on the inputs d0, d−1, and kl. Analogously, for a general p ∈ ℕ, applying (A6) recursively until a period p′ ∈ Γ\ℕ is reached, it can be seen that Kp is a function of inputs at periods lesser than or equal to p. Indeed:
K p   =   K p l + d p d p 1 e   =   K p 2 l + d p l d p l 1 e + d p d p 1 e   =   k p m p l + d p ( m p 1 ) l d p ( m p 1 ) l 1 e + + d p l d p l 1 e + d p d p 1 e
where:
m p   =   min { m / p m l   <   0 }

Appendix B. OIP Sensitivity with Respect to the Demand

This subsection obtains some mathematical relationships between the variation ΔDp in the demand dp at period p, and the variation ΔKp of the OIP. Any variation ΔDp not only affects Kp but also Kp+1 since, applying (3) for p + 1:
K p + 1   =   k p m p + 1 l + d ( p + 1 ) ( m p + 1 1 ) l d ( p + 1 ) ( m p + 1 1 ) l 1 e + + d ( p + 1 ) l d ( p + 1 ) l 1 e + d p + 1 d p e
Note that since l > 1 only one term in (A7) includes dp (the last one). Furthermore, applying (3) to p+nl, n∈ℕ, ΔDp also affects future investments Kp+n·l:
K p + n l   =   k p m p + n l l + d ( p + n l ) ( m p + n l 1 ) l d ( p + n l ) ( m p + n l 1 ) l 1 e + + d p d p 1 e + + d ( p + n l ) l d ( p + n l ) l 1 e + d p + n l d ( p + n l ) 1 e
Note that (A8) is a generalization of (3) (since (3) is (A8) for n = 0). Analogously, generalizing (A7), ΔDp also affects Kp+n·l+1, n ∈ ℕ:
K p + n l + 1   =   k p m p + n l + 1 l + d ( p + n l + 1 ) ( m p + n l + 1 1 ) l d ( p + n l + 1 ) ( m p + n l + 1 1 ) l 1 e + + d p + 1 d p e + + d ( p + n l + 1 ) l d ( p + n l ) l e + d p + n l + 1 d p + n l e
From (A8) and (A9), it is proven that for n∈ℕ:
Δ K p Δ D p   =   { 1 e   i f   p   =   p + n l   1 e   i f   p   =   p + n l + 1 0         o t h e r w i s e

Appendix C. Prices to Recover Investment Costs

Prices that recover investment costs fp over the generation assets’ lifespan should satisfy:
f p   =   τ   =   0 l 1 e P p + τ ( 1 + r ) τ
Note that fp is assumed to be satisfied with prices Pp in years p in the lifespan of the capacity installed at p′. Applying (A11) to p= p + nl and p= p + nl + 1, n∈ℕ, dividing the later by (1 + r), and subtracting both:
f p + n l f p + n l + 1 ( 1 + r )   =   e P p + n l e P p + ( n + 1 ) l ( 1 + r ) l
Applying (1) to fp+nl+1 in (A12) and then to fp+nl:
f p + n l ( 1 1 + i 1 + r )   =   e ( P p + n l P p + ( n + 1 ) l ( 1 + r ) l ) f p e ( 1 + i ) n l ( r i 1 + r )   =   P p + n l P p + ( n + 1 ) l ( 1 + r ) l
Dividing (A13) by (1+r)nl:
f p e ( 1 + i 1 + r ) n l ( r i 1 + r )   =   P p + n l ( 1 + r ) n l P p + ( n + 1 ) l ( 1 + r ) ( n + 1 ) l
Note that according to(A14), for any k ∈ ℕ:
n   =   0 k { f p e ( 1 + i 1 + r ) n l ( r i 1 + r ) }   =   P p P p + ( k + 1 ) l ( 1 + r ) ( k + 1 ) l
Therefore, taking limits in k:
f p e ( r i 1 + r ) n   =   0 ( 1 + i 1 + r ) n l   =   P p l i m k P p + ( k + 1 ) l ( 1 + r ) ( k + 1 ) l
Denoting Limp as the above limit and clearing Pp, (A16) is:
P p   =   f p e s + L i m p
where:
s   =   1 ( r i 1 + r ) n   =   0 ( 1 + i 1 + r ) n l   =   1 ( 1 + i 1 + d ) l ( r i 1 + r )   =   p   =   0 l 1 ( 1 + i 1 + r ) p
which converges if i < r (see Section 2). Note that parameter s is a function of RDF.
Substituting (A17) in (A11), and taking into account (1) and (A18):
f p   =   τ   =   0 l 1 [ e ( 1 + r ) τ ( f p + τ e s + L i m p + τ ) ]   =   f p 1 s τ   =   0 l 1 ( 1 + i 1 + r ) τ + e τ   =   0 l 1 L i m p + τ ( 1 + r ) τ   =   f p + e τ   =   0 l 1 L i m p + τ ( 1 + r ) τ
Therefore:
τ   =   0 l 1 L i m p + τ ( 1 + r ) τ   =   0
Since Limp ≥ 0, (A20) holds if and only if Limp is 0 for every period p. Therefore, from (1) and (A17), prices to recover investment costs satisfy (6).

Appendix D. Considering Corporate Taxes

This annex adapts (6) when there are corporate taxes charged on the profits generated from the exploitation of the investments. Profits are computed by subtracting book amortizations from incomes, that is, from the allocation of the investment cost fp over the periods p′ in which the assets are used.
Let us consider that the amortization period coincides with the lifespan of an asset and that the amortization is represented through a percentage aτ+p of fp for the l periods τ starting from p, τ∈{0, …, l − 1}. Then:
τ   =   0 l 1 a p + τ   =   1
Similarly to (A11), to recover investment costs, prices Pp must satisfy:
f p + t a x ( τ   =   0 l 1 e P p + τ ( 1 + r ) τ τ   =   0 l 1 f p a p + τ ( 1 + r ) τ )   =   τ   =   0 l 1 e P p + τ ( 1 + r ) τ
being tax the percentage of taxes applied to the profit, assuming that incomes are greater than or equal to the amortizations. Extracting common factor to the right-hand side of A22) and clearing, Equation (A22) is equivalent to:
f p ( 1 t a x c a p 1 t a x )   =   τ   =   0 l 1 e P p + τ ( 1 + r ) τ
where cap is the present value, at period p, of the summation of the percentages aτ+p for τ = 0, …, l − 1, i.e.,:
c a p   =   τ   =   0 l 1 a p + τ ( 1 + r ) τ
Note that, according to(A21), cap < 1 and therefore the left-hand side of (A23) is higher than fp, which can be interpreted as the investment cost fp including the effect of taxes, i.e.,:
f p   =   f p ( 1 t a x c a p 1 t a x )
Following the same reasoning as in the previous subsection, but starting from (A23) and (A25) instead of (A11), the prices Pp taking taxes into account can now be computed as in (6) but using fp instead of fp.

Appendix E. Prices to Recover Investment Costs as Dual Variables of C

By definition, Lp is:
L p   =   l i m Δ D p 0 Δ C Δ D p   =   p   =   0 f p l i m Δ D p 0 ( Δ K p / Δ D p ) ( 1 + r ) p
According to (5) (note that ΔKpDp is independent of ΔDp):
L p   =   { n   =   0 f p + n l l i m Δ D p 0 ( Δ K p + n l Δ D p ) ( 1 + r ) p + n l + n   =   0 f p + n l + 1 l i m Δ D p 0 ( Δ K p + n l + 1 Δ D p ) ( 1 + r ) p + n l + 1 }   =   { n   =   0 f p + n l ( 1 / e ) ( 1 + r ) p + n l n   =   0 f p + n l + 1 ( 1 / e ) ( 1 + r ) p + n l + 1 }
and, therefore, (1), (A18), and (6) leads to:
L p   =   { n { : p + n l } f 0 e ( 1 + i 1 + r ) p + n l n { : p + n l + 1 } f 0 e ( 1 + i 1 + r ) p + n l + 1 }
and finally:
L p   =   f 0 e ( 1 + i 1 + r ) p ( r i 1 + r ) n   =   0 ( 1 + i 1 + r ) n l   =   f 0 e ( 1 + i 1 + r ) p 1 s   =   f p s e 1 ( 1 + r ) p   =   P p ( 1 + r ) p
This last equation proves that Lp is indeed the discounted price Pp of (6).

Appendix F. Prices to Recover Investment Costs as Dual Variables of C with FDP

The dual variable Lp of (9) when minimizing C using f p , p a in (14) with Γp = {p, …, min(p, p + l − 1)}:
L p   =   l i m Δ D p 0 Δ C Δ D p   =   p   =   0 p _ { p   =   p min ( p _ , p + l 1 ) f p , p a l i m Δ D p 0 ( Δ K p / Δ D p ) ( 1 + r ) p }
Applying now (1) and (19) (note that ΔKpDp is independent of ΔDp):
L p   =   p   =   0 p _ [ f p ( 1 + r ) p s p a s ( Δ K p / Δ D p ) { p   =   p min ( p _ , p + l 1 ) ( 1 + a 1 + r ) p p } ]   =   f 0 p   =   0 p _ [ ( 1 + i 1 + r ) p s p a s ( Δ K p / Δ D p ) { p   =   p min ( p _ , p + l 1 ) ( 1 + a 1 + r ) p p } ]   =   f 0 p   =   0 p _ [ ( 1 + i 1 + r ) p 1 s ( Δ K p / Δ D p ) { p   =   p min ( p _ , p + l 1 ) ( 1 + i 1 + r ) p p } ]   =   f 0 p   =   0 p _ [ 1 s ( Δ K p / Δ D p ) { p   =   p min ( p _ , p + l 1 ) ( 1 + i 1 + r ) p } ]
From (5):
L p   =   { n { : p + n l p _ } f 0 s e p   =   p + n l min ( p _ , p + n l + l 1 ) ( 1 + i 1 + r ) p n { : p + n l + 1 p _ } f 0 s e p   =   p + n l + 1 min ( p _ , p + n l + l ) ( 1 + i 1 + r ) p }
In this expression, each term except the first one (for n = 0) in the first summand is the opposite of a term in the second summand. Effectively, if np = max{n ∈ ℕ: p + nl + 1 ≤ p}, the first summand can be reduced to:
n { : p + n l p _ } p   =   p + n l min ( p _ , p + n l + l 1 ) ( 1 + i 1 + r ) p   =   n   =   0 n p [ p   =   p + n l min ( p _ , p + n l + l 1 ) ( 1 + i 1 + r ) p ]   =   [ ( 1 + i 1 + r ) p + ( 1 + i 1 + r ) p + 1 + + ( 1 + i 1 + r ) p + l 1 ] + [ ( 1 + i 1 + r ) p + l + ( 1 + i 1 + r ) p + l + 1 + + ( 1 + i 1 + r ) p + l + l 1 ] + + [ ( 1 + i 1 + r ) p + n p l + ( 1 + i 1 + r ) p + n p l + 1 + + ( 1 + i 1 + r ) p _ ]   =   u   =   p p _ ( 1 + i 1 + r ) u
and, analogously, the second summand becomes:
n   =   0 n p [ p   =   p + n l + 1 min ( p _ , p + n l + l ) ( 1 + i 1 + r ) p ]   =   u   =   p + 1 p _ ( 1 + i 1 + r ) u
Therefore:
L p   =   { f 0 s e u   =   p p _ ( 1 + i 1 + r ) u f 0 s e u   =   p + 1 p _ ( 1 + i 1 + r ) u }
Therefore, the dual variable Lp coincides exactly with Pp/(1 + r)p, the first term in the first summand (see (6)):
L p   =   f 0 s e ( 1 + i 1 + r ) p   =   f p s e ( 1 + r ) p   =   P p ( 1 + r ) p
Note that no approximation is needed.

References

  1. Pérez-Arriaga, I.J. Regulation of the Power Sector; Springer Science & Business Media: Berlin, Germany, 2014; ISBN 978-1-4471-5034-3. [Google Scholar]
  2. Genoese, F.; Drabik, E.; Egenhofer, C. The EU Power Sector Needs Long-Term Price Signals. Social Science Research Network, SSRN Scholarly Paper ID 2782355. 2016. Available online: https://papers.ssrn.com/abstract=2782355 (accessed on 30 May 2019).
  3. Mastropietro, P. Regulatory Design of Capacity Remuneration Mechanisms in Regional and Low-Carbon Electric Power Markets. Ph.D. Thesis, KTH Royal Institute of Technology, Stockholm, Sweden, 2016. [Google Scholar]
  4. Eckaus, R.S.; Parikh, K.S. Planning for Growth; MIT Press: Cambridge, MA, USA, 1968. [Google Scholar]
  5. Conejo, A.J.; Baringo, L.; Kazempour, S.J.; Siddiqui, A.S. Investment in Electricity Generation and Transmission: Decision Making under Uncertainty; Springer: Berlin/Heidelberg, Germany, 2016; ISBN 978-3-319-29501-5. [Google Scholar]
  6. Ramos, A.; Perez-Arriaga, I.J.; Bogas, J. A nonlinear programming approach to optimal static generation expansion planning. IEEE Trans. Power Syst. 1989, 4, 1140–1146. [Google Scholar] [CrossRef]
  7. Ventosa, M.; Denis, R.; Redondo, C. Expansion planning in electricity markets. Two different approaches. In Proceedings of the 14th Power Systems Computation Conference 2002 (PSCC 2002 SEVILLE), Sevilla, Spain, 24–28 June 2002; pp. 1–7. [Google Scholar] [CrossRef]
  8. Murphy, F.H.; Smeers, Y. Generation Capacity Expansion in Imperfectly Competitive Restructured Electricity Markets. Oper. Res. 2005, 53, 646–661. [Google Scholar] [CrossRef]
  9. Wang, J.; Shahidehpour, M.; Li, Z.; Botterud, A. Strategic Generation Capacity Expansion Planning With Incomplete Information. IEEE Trans. Power Syst. 2009, 24, 1002–1010. [Google Scholar] [CrossRef]
  10. Kazempour, S.J.; Conejo, A.J.; Ruiz, C. Strategic Generation Investment Using a Complementarity Approach. IEEE Trans. Power Syst. 2011, 26, 940–948. [Google Scholar] [CrossRef]
  11. Hinojosa-Mateus, V.H. Static generation capacity expansion planning using linear transmission distribution factors. In Proceedings of the 2016 IEEE PES Transmission Distribution Conference and Exposition-Latin America (PES T&D-LA), Morelia, Mexico, 20–24 September 2016; pp. 1–6. [Google Scholar]
  12. van Stiphout, A.; Vos, K.D.; Deconinck, G. The Impact of Operating Reserves on Investment Planning of Renewable Power Systems. IEEE Trans. Power Syst. 2017, 32, 378–388. [Google Scholar] [CrossRef]
  13. Park, H.; Baldick, R. Stochastic Generation Capacity Expansion Planning Reducing Greenhouse Gas Emissions. IEEE Trans. Power Syst. 2019, 30, 1002–1010. [Google Scholar] [CrossRef]
  14. Gil, E.; Aravena, I.; Cárdenas, R. Generation Capacity Expansion Planning Under Hydro Uncertainty Using Stochastic Mixed Integer Programming and Scenario Reduction. IEEE Trans. Power Syst. 2015, 30, 1838–1847. [Google Scholar] [CrossRef]
  15. Shengyu, W.; Lu, C.; Xiaoqing, Y.; Bo, Y. Long-term generation expansion planning under uncertainties and fluctuations of multi-type renewables. In Proceedings of the 2015 IEEE 5th International Conference on Power Engineering, Energy and Electrical Drives (POWERENG), Riga, Latvia, 11–13 May 2015; pp. 612–616. [Google Scholar]
  16. Gorenstin, B.G.; Campodonico, N.M.; Costa, J.P.; Pereira, M.V.F. Power system expansion planning under uncertainty. IEEE Trans. Power Syst. 1993, 8, 129–136. [Google Scholar] [CrossRef] [Green Version]
  17. Domínguez, R.; Conejo, A.J.; Carrión, M. Investing in Generation Capacity: A Multi-Stage Linear-Decision-Rule Approach. IEEE Trans. Power Syst. 2016, 31, 4784–4794. [Google Scholar] [CrossRef]
  18. Wogrin, S.; Barquín, J.; Centeno, E. Capacity Expansion Equilibria in Liberalized Electricity Markets: An EPEC Approach. IEEE Trans. Power Syst. 2013, 28, 1531–1539. [Google Scholar] [CrossRef]
  19. Wogrin, S.; Centeno, E.; Barquin, J. Generation Capacity Expansion in Liberalized Electricity Markets: A Stochastic MPEC Approach. IEEE Trans. Power Syst. 2011, 26, 2526–2532. [Google Scholar] [CrossRef]
  20. Li, Q.; Wang, J.; Zhang, Y.; Fan, Y.; Bao, G.; Wang, X. Multi-Period Generation Expansion Planning for Sustainable Power Systems to Maximize the Utilization of Renewable Energy Sources. Sustainability 2020, 12, 1083. [Google Scholar] [CrossRef] [Green Version]
  21. González, P.; Villar, J.; Díaz, C.A.; Campos, F.A. Joint energy and reserve markets: Current implementations and modeling trends. Electr. Power Syst. Res. 2014, 109, 101–111. [Google Scholar] [CrossRef] [Green Version]
  22. Manne, A.S. Sufficient Conditions for Optimality in an Infinite Horizon Development Plan. Econometrica 1970, 38, 18. [Google Scholar] [CrossRef]
  23. Williams, J.B. The Theory of Investment Value; Fraser Publishing Co.: Burlington, VT, USA, 1997; ISBN 978-0-87034-126-7. [Google Scholar]
  24. Renewable Power Generation Costs in 2014. Available online: https://irena.org/-/media/Files/IRENA/Agency/Publication/2015/IRENA_RE_Power_Costs_2014_report.pdf (accessed on 1 January 2015).
  25. de Sisternes, M.V.F. Investment Model for Renewable Electricity Systems (IMRES): An Electricity Generation Capacity Expansion Formulation with Unit Commitment Constraints; MIT Press: Cambridge, MA, USA, 2013. [Google Scholar]
Figure 1. Approximation error of C for finite dual pricing (FDP).
Figure 1. Approximation error of C for finite dual pricing (FDP).
Sustainability 13 01993 g001
Figure 2. C is a valid truncation of C for FDP.
Figure 2. C is a valid truncation of C for FDP.
Sustainability 13 01993 g002
Table 1. Inputs of the case example.
Table 1. Inputs of the case example.
f0 (€/MW/Year)l (Years)e (Hours)i (%)r (%)
1.4⋅10652193210
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campos, F.A.; Villar, J.; Centeno, E. Annualization of Renewable Investment Costs for Finite Horizon Electricity Pricing and Cost Recovery. Sustainability 2021, 13, 1993. https://doi.org/10.3390/su13041993

AMA Style

Campos FA, Villar J, Centeno E. Annualization of Renewable Investment Costs for Finite Horizon Electricity Pricing and Cost Recovery. Sustainability. 2021; 13(4):1993. https://doi.org/10.3390/su13041993

Chicago/Turabian Style

Campos, Fco. Alberto, José Villar, and Efraim Centeno. 2021. "Annualization of Renewable Investment Costs for Finite Horizon Electricity Pricing and Cost Recovery" Sustainability 13, no. 4: 1993. https://doi.org/10.3390/su13041993

APA Style

Campos, F. A., Villar, J., & Centeno, E. (2021). Annualization of Renewable Investment Costs for Finite Horizon Electricity Pricing and Cost Recovery. Sustainability, 13(4), 1993. https://doi.org/10.3390/su13041993

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop