Research Article  Open Access
Generation Expansion Models including Technical Constraints and Demand Uncertainty
Abstract
This article presents a Generation Expansion Model of the power system taking into account the operational constraints and the uncertainty of longterm electricity demand projections. The model is based on a discretization of the load duration curve and explicitly considers that power plant ramping capabilities must meet demand variations. A model predictive control method is used to improve the longterm planning decisions while considering the uncertainty of demand projections. The model presented in this paper allows integrating technical constraints and uncertainty in the simulations, improving the accuracy of the results, while maintaining feasible computational time. Results are tested over three scenarios based on load data of an energy retailer in Colombia.
1. Introduction
Decisions related to investments in power plants, known as Generation Expansion Planning (GEP), are often based on the load duration curve (LDC). Linear optimization models are already used for a long time to find the optimal mix of generation technologies supplying this demand. The LDCbased solutions describe the required installed generation capacity of each technology and their operating hours, in order to minimize the total cost, that is, the investment and operational cost. Since traditional models do not include the variations of the demand over time, they do not take into account the ramping and cycling capabilities of power plants. However, not including these technical constraints of power systems underestimates the flexibility needs of the system. This fact becomes increasingly relevant facing the rising trend of participation of renewable generation characterized by a variable profile of supply. Traditional approaches based on the LDC are not considering the potential value of flexibility provided by different generation technologies, and this results in a suboptimal solution. On the other hand, solutions that include these technical elements require high computational efforts. The fact that traditional solutions do not take into account the uncertainty of the demand forecast results in overcapacity or reliability issues regarding the installed capacity. As a consequence, suboptimal decisions of investments are associated with the longterm use of power plants.
The need to include technical constraints and uncertainty in the GEP requires the integration of new elements. However, these new elements increase the complexity of the problem. The same difficulty is faced while using scenarios which consider uncertainty but improve the ability to deal with longterm planning decisions and make the results robust in front of future trends in the power system. The state of the art in the generation investment problem shows different approaches to take the integration of new variables and uncertainties into account. A first set of approaches deals with the LDC model itself. In [1], a linear piecewise function is used to represent the German and UK LDC; a residual LDC is used to integrate renewable sources in the part of the curve where most changes of demand take place. Technical improvements such as ramp rates and the maximum on time of a generating units are suggested to be included in the LDC methodology [2]. In [3], a robust optimization that includes uncertainty of demand, the electricity price, and generation cost is proposed. A second set of approaches deals with the problem of technical constraints integration into the classical economic generation investment model, that is, based on the time series representation of the demand. Authors of [4] include several additional variables and constraints related to the power system to solve a single period problem. In [5], a relevant discussion about how dynamic constraints such as ramping rates and startup cost are changing the optimal point of the classical generation investment solution is presented. The work stresses the need of methodologies exploring the solution in a prediction time horizon also called multiperiod solutions. A mixed solution proposed in [6] considers the integration of renewable generation in the GEP by means of net demand curve; to deal with prediction errors, an operating reserve constraint is imposed. The solution includes the ramp rate of the conventional technologies to cover the uncertainty of the wind generation power.
Finally, a set of detailed solutions makes use of stochastic and time series based models to solve the GEP. A generation plan expansion methodology using a multiperiod horizon solved by a bilevel problem is proposed in [7]. The problem maximizes the investment profits in the top level and the second level deals with technical constraints such as generation, transmission constraints, and load flows. The impact of renewable sources in the GEP is discussed in [8] with a model optimizing all the yearly investments, fixed and variable operational cost, reserve costs, and ramping and startup costs constraints. The problem deals with a new optimal point to balance the increasing amount of renewable with reserves to meet the system reliability requirements. Last, in book [9], the longterm investments considering the GEP problem integrated with transmission investments are discussed. This book presents a full methodology based on stochastic processes. The solutions presented in this book take into account the integration of several planning horizons by means of multistage models with detailed technical constraints. Nevertheless, obtained results are based on decision trees and exhaustive solutions and, due to the very detailed models, results obtained in this book require a very high computational effort. Now, as explored in this review, time series and detailed solutions lead to use of complex methodologies that imply a high computational cost. The discussion of stochastic methods used to forecast power variables is not presented in this review. The scope of this state of art is to tackle the investment problem from the generation mix design. Concluding, new methodologies must deal with a tradeoff between the detailed level and the required time to solve the problem in a large time horizon. Details are not paramount, planning exercises and scenarios analysis can use more efficiently solutions with less technical and forecast details. This work focuses on an optimal, flexible, and efficient multiperiod methodology to solve the GEP problem. The solution is inspired by the MPC (model predictive control) methodology. The solution includes ramping constraints creating a compromise between the power system flexibility and the installed capacity. This leads to a new optimal point which guarantees a reliable operation minimizing the fixed and variable investments cost. Four scenarios are simulated with uncertainties in the variables showing the advantages as a planning tool. The paper is organized as follows: Section 2 presents a brief description of the methodology. Sections 2.1 and 2.2 present the used LDC discretization and the dynamic analysis made to include the ramping constraints. Section 3 presents the comparison between the proposed model that includes the ramping constraints and the classic economic solution. Section 4 presents a brief description of the MPC and the approach used for the multiperiod GEP methodology. Lastly, Section 5 presents the used data, the test scenarios, and the model performance comparison of Section 3 finishing with the results of the MPC planning.
2. Methodology
As mentioned, the classic solution of the GEP based on the LDC only considers economic variables. Additionally, the method does not consider the time series of the demand to create the LDC information. New generation expansion plans should consider a wide selection of technologies and uncertainties; each technology has their own technical constraints and economic cost. An optimal temporal LDC partition allows distributing the GEP problem in several blocks or demand levels. Each level has his own features: required power, uncertainties, and power changes. Including the demand levels in the problem can be used to include dynamic analysis in the solution as technical constraints for the yearly operation. The constraint guarantees a minimum power change response for each block. The proposed model provides the optimal generation expansion plan cost with a reliable operation design considering the power changes in each block. Notice that, considering a maximum power change for each block, all the feasible changes are covered, and as a consequence, the LDC solution is robust for all possible system changes in the planning period for each block. Therefore, the optimal economic power mix problem is solved integrating the ramp velocities in the economic model. A tradeoff between the economic and the technical part is created. Finally, the proposed model only solves the oneyear problem and it is extended by means of the model predictive control (MPC). The solution of the problem in a prediction horizon provides a new optimal solution that considers the predicted changes in the variables. This methodology is easily improved by including additional constraints and models into the problem.
2.1. Load Duration Curve Discretization
The LDC curve is a continuously decreasing function that represents the amount of time that a certain electricity demand is requested to the power system (see Figure 1). In this section, a LDC discretization is proposed, based on [10, 11]. The proposed LDC discretization minimizes the difference between the area and the area composed by several blocks based on temporal decomposition. As shown in Figure 1, using the time axis with values as optimization variables, the width of each discretized block corresponds to the time duration of the block along one year. Analogously, the height of each block represents the maximum required power (or installed capacity) and is given by the function. Now, with and already defined, the approximation of , as a discrete function considering that is a finite number, and , will have an error associated that can be expressed asIf function is approximated piecewise with trapezoids, the accuracy is improved and the result preserves the LDC decreasing characteristic. Hence, it is possible to say that
With the described trapezoidal approximation, the named LDC block discretization will always be greater than the LDC function for a finite number of blocks. This provides a confidence interval to generate a maximum amount of energy. The next step minimizes the approximation error between the trapezoid LDC approximation and the energy blocks. If are the mesh points for the time discretization and the hours of the year, then, it is possible to define the intervals in hours [h] as optimization variable. Now, with the LDC function, the power value for the point is obtained. In consequence, the product between and will produce block that represents the energy required by the country along the year. The optimization problem is formulated as follows:
Using the optimal values that minimizes the LDC discretization error, for simplicity, the optimal trapezoids will be named as energy blocks and represented as follows:
This work considers three blocks to discretize the LDC function . The optimal discrete LDC approach presents the following characteristics for each block: block has the lowest variation and the minor approximation error. From the financial and operative point of view, this block represents the load related to forward energy agreements. This energy has a small volatility and is dispatched with base load generation technologies. Block represents an increasing variable energy consumption; this block includes load composed by energy agreements, residential customers, and companies. The volatility of this block is bigger than block and is correlated with the consumption patterns of regulated users. Also block energy demand is usually generated with technologies which had affordable prices for big amounts of energy but also should be flexible and with a fast response to the load variations. Finally, the demand block represents the higher demand points. In one year, these peak demand points have a duration of no more than one hour, and their appearances during the day are few. These low frequency peaks create a flexibility market opportunity; peak demand points imply the need to have a corresponding generation installed capacity available, even if it is only used few hours in the year. It is a duty from the generator and the energy retailers to have the required energy available at any time. In this sense, demand response and load shift among other smart grid technologies can be used to supply the required energy at peak hours; for now, this work wants to deal only with installed capacity to supply the energy.
2.2. Ramping Constraint Model
The discretization proposed in Section 2.1 allows including particular constraints for each block. Particularly, it is possible to find the optimal economic mix of generation technologies including technical constraints of the generation portfolio. The proposed management strategy optimizes the generation technologies mix according to the ramp velocities criteria to meet the power changes values on each block. Considering the notion that every change in the load implies a change in the generation system is true that results in a change in the operational point of one or several power plants of the system to meet the demand denoted as . Then, it is possible to establish that for each . The generation technology or mix must be able to reach the new operation point in a proper time. Considering a power generation matrix with generation technologies available, it is possible to meet the operation changes in the system with a mix of the capacities of each technology represented as a weighted sum, where is the participation of the th technology in the power change and is the ramp rate of each technology: Defining as the produced energy for a generation technology, the variable is correlated with the generation cost and installed capacity. The weighted portfolio approach formulated can be included according to a particular design parameter in the GEP changing the obtained solution. Now, this work will use the load time series analysis described in Section 2.3 to find the values for demand levels. These values will be used as the minimum ramping value that the solution must guarantee as follows:
2.3. Methodology Used for Time Series Analysis and Description of the Data Used
In order to identify the value in each demand level as measurement of the higher energy changes in each block, the methodology calculates the time ahead demand time series derivative. To do this, the time series is matched with each demand level obtained in Section 2.1 (see Figure 2). The derivatives of each block are shown in Figure 3. Finally, the values are obtained from the maximum derivative value for each block and presented in Table 1. The comparison between the identified and the maximum values for each demand level is shown in Table 1. The maximum values represent the biggest change defined by the block size. The data used is described in Section 5.

3. Ramping Constrained Model Approach
With the values of each demand block, the next step is to find the optimal combination of generation technologies to supply the demand at the lowest cost taking into account the ramping constraints. The formulation presented in (6) allows combining the technologies ramp velocities by means of the values to guarantee a desired performance for each demand block. The integration is made balancing the generation per block of each generation technology with the values. The equality guarantees the installed capacity of each block. The dynamic constraint puts the desired operative level for the ramping combination in each block. Finally, considering the operation of each plant and the installed capacity are also restricted by the variables; this constraint joint with extends the impact to the installed capacity. The reference for the used variables is presented in Table 2 and the classic and dynamic formulations are presented in (9). To solve the optimization problems, letwhere represents the system variables: installed capacity, use of each technology, and the weights of the dynamic solution for each technology.

where is the vector with the generation technologies fixed and variable cost per year. The values represent the annualized installation capacity and fixed operation cost for each generation technology; represent the sum of the operation and maintenance (O&M) cost and the fuel cost for each technology. The use of the annualized costs shares the installed capacity price over the lifespan of each technology. Notice that the cost of the values is zero. In order to compare the classic and the ramping constrained models, the same cost function is used.
4. Model Predictive Control
The solution of the generation investment planning in a prediction time window is relevant as mentioned in [5]. Here, the merit order dispatch is discussed in front of the emergent dynamic constrains of the power system. The inclusion of variables such as green generation technologies, demand response plans, fuel prices, and storage technologies among others variables becomes necessary to deal with the future scenarios of the power system. The integration of these new variables increases the amount and complexity of the models used in the power system planning and operation. Control applications and control based management strategies had proven their efficiency to achieve good and desired performance in several dynamic and complex systems. The dynamic constrained generation expansion problem formulated in Section 2 found an optimal economic point that meets the imposed LDC dynamics for a single year. The new challenge is to propose an analog dynamic problem formulation that considers the future scenarios and the previous installed capacities in a set of optimal solutions. A wellknown tool used to solve constrained dynamic linear and nonlinear problems providing the optimal trajectory for based on the future problem states is the model predictive control (MPC); using the receding horizon principle after the computation of the optimal control sequence applies only the first control actions; then, the next iteration solution provides a new control sequence including the new information available and the previous system states. The loop continues until a desired horizon is reached. This kind of solution provides a planning tool able to analyze several scenarios at the same time, providing a set of optimal trajectories for the variables and control actions that allow analyzing and determining, for each time, the generation system evolution and his related returns including fixed and variable cost.
4.1. Model Predictive Controller Formulation
To create an optimal generation planning tool, this work implements an MPC controller. MPC makes use of the receding horizon strategy that allow solving the ramping constrained model in a prediction horizon. The solution calculates the optimal decisions for the GEP considering future LDC values. Using the model proposed in (9) as an economical objective function and including the proper constraints, the optimal solution formulation is presented considering a time discrete state space model:where are the system states, are the system outputs, and is the current control vector. Using (6), the state space representation of the problem for generation technologies and a LDC discretized in blocks can be written as follows:
As mentioned, infrastructure investments have a natural condition for the variables: they cannot decrease. Negative changes in the installed capacity would mean a “destruction” or dismantling of the plant; this in the real operation implies to shut down the plant or decrements in the plant efficiency. In this work, negative values are not considered, plant aging is reduced by considering the maintenance cost in the fixed operation cost. However, the aging parameter can be easily included in the space state model. where vector represents the control actions taken to modify the variables; represent the decisions taken to increase the installed capacity; represent the expected use of the generations technologies; and represent the weights of the generation technologies according to the ramp velocities. The interactions between the and the variables are described by and matrices as follows:where represents the states evolution. Thus, the installed capacities are the only variables considered persistent in the system (integration variables). The installed capacity at step depends only on the previous state and the control actions explicit in the matrix. is an identity matrix. The matrix indicates the relationship between the decision variables and the state variables. In this case, there is only a direct influence of each control variable over its corresponding state.The vector is the same as described in (8); the use of the same cost vector in all optimizations cost functions allows directly comparing the results of the presented methodologies. The relationships related to the constraints are explicit in the MPC as hard constraints. Changes in the model output as a consequence of the control actions are not considered; therefore . Now, with the dynamic state space system written, the traditional economic MPC problem based in [12] is defined as follows:
At time step , let and be the state trajectory and the control sequences, with being the prediction horizon and being the economic stage cost. The system is subject to hard constraints on state , output , and input for all , where , , and are closed sets. The objective functions of (9) are clearly a linear combination of the states. Using the MPC methodology, a great flexibility in the objective function is available. In order to calculate the optimal control solution , the cost function will be based on the state space output described in (10) as and the economic MPC formulation is given by
5. Case Study
As an application case, the 2013 year demand data from EPM (Empresas Publicas de Medellin) that is one of the biggest energy retailers from the center region of Colombia sharing the 23.9% of the maximum 25% allowed by the country regulation is used. The demand is composed by regulated and nonregulated users; in Colombia, there are not differential energy prices, and even regulated costumers ignore the realtime energy price. Therefore, the regulated load does not have an energy saving culture making use of the energy at any time along the day increasing the load volatility. Renewable generation technologies are not mature in the market, neither are demand response or similar technologies. The generation technologies used in the problem are selected considering competitive technologies in terms of cost and ramping rates into the problem; the technical values are obtained from [8, 13] and presented in Table 3. The results compare the traditional economical LDC solution without ramping constraints and the ramping constrained LDC methodology described in Section 2.1. Three additional artificial scenarios shown in Figure 4 are created based on the data set. The first scenario is taken from the demand time series presented in Figure 3 and the demand levels obtained from the discretization presented in Section 2.1. The values of the first scenario correspond to the values presented in Table 1. The remaining scenarios are hypothetical and they are designed to stress the methodology in front of positive and negative LDC and dynamics changes as follows: second scenario presents a 10% growth in all the blocks for the and values assuming a positive development in the country. The third scenario assumes an additional 10% increase in block 1 and a 30% reduction in block 2. Finally, the fourth scenario simulates a 10% growth in block 3 load.

5.1. Economic and Ramping Constrained Comparison
The installed capacities for each solution are presented in Figure 5 and the generation mix use for each scenario is presented in Figures 6 and 7; the traditional model without ramping constraints and the proposed ramping constrained model provides the same solution for the first scenario: a base load technology operating in blocks and with higher energy demand. An OCGT (opencycle gas turbine) generation plant is the choice in the block providing the energy with a higher cost but the lowest install capacity cost (see Figures 5 and 6). In this solution, the same technology mix meets the ramping requirements. Then, in scenario number two, the absence of ramping constrains increases the installed capacity mix of scenario 1 proportionally to the increase of each demand level. The lack of ramping constraints makes the solution blind to the dynamic changes of the demand. In contrast, a noticeable change in the generation mix is made by the ramping constrained solution compensating the demand increment using CCGT (combinedcycle gas turbine) in the solution. This change increases the cost but represents the lowest cost possible to meet the ramping requirements. The tradeoff between capacity and the response flexibility of the system has been shown. Scenario three presents a particular case where the classic and the ramping constraint solutions are the same. The derivative value in this scenario is included in the feasible region where the classic solution has his optimal value. Scenario three shows that, in relaxed cases where the dynamics are not so severe (low load values), the classic and the ramping constrained solutions could be the same. Finally, scenario four shows how a increment in the base load can move the optimal point increasing the coal installed capacity and generating the energy in all the blocks only with the cheapest technology. The result generates an interesting solution where a big installed capacity of a slow generation technology could meet the imposed load dynamics.
5.2. Model Predictive Control Results
The use of MPC as generation investment planning tool results are presented. As designed in matrix, the MPC should keep the installed capacity between scenarios and increases the values if it is needed. Comparing the installed capacity of the ramping constrained solution with the MPC results in Figure 8, the previous capacity values are preserved. When a previous installed capacity is already built, the optimization has the choice to turn off a particular technology and to pay the fixed O&M cost and supply the energy with another technology. The second scenario as seen in Section 5.1 pushes the system to another technology mix including CCGT as base generation plant increasing the energy production cost. The MPC keeps the previous coal capacity of scenario one and installs a small CCGT installed capacity amount to fulfill the problem requirements. This behavior creates the opportunity to include and manage demand response technologies designed to supply small amounts of energy with low installation cost. Scenarios three and four present smaller load values. The MPC keeps the previously installed capacities avoiding the negative changes presented in Figure 5.
With the installed capacities described in Figure 8, the MPC proposes a different energy generation mix for each scenario presented in Figure 9. Now, having the installed capacities and fixed O&M cost, stabilized generation costs are the key to achieve the economic optimal. The dynamic imposed to the installed capacity allows deciding, in the planning horizon, which technology is the most appropriated to work in a specific scenario selecting which plants should be turned on or off. The strategy offers a planning tool making use of the cheapest technology efficiently and covers high load dynamics with fastest but expensive technologies.
Finally, the economic cost presented solutions are compared in Table 4. The presented values correspond to the objective function of each solution. As expected, the traditional solution is always the cheapest solution. The ramping constraint increased the solution cost in scenarios two and four where the load and ramping condition increased. Finally, the MPC faced the scenarios keeping the solution cost close to the traditional solution but at the same time meeting the solution cost. Finally, the MPC solution minimizes the economic overall cost meeting the ramping constrains and staying close to the traditional methodology cost.

6. Conclusions
The above presented work proposed an optimal generation planning tool based on a modified GEP. The GEP is improved including a technical constraint that considers the generation technologies ramp velocities and the load dynamics in the optimization. The constraint includes the ramping values of the generation technologies and mix them to satisfy a design operating condition optimizing to the lowest cost of the installed capacity and fixed and variable values cost of the generation portfolio. Then, using MPC theory, the optimal longterm planning problem is solved providing a new generation mix that is reliable in front of load dynamics and keeps a minimal economic cost. The LDC discretization used in the GEP is improved minimizing the difference between the LDC and three demand levels. As a result, generated levels are used to classify the load time series dynamics and also represent the expected energy required by each level per year. This representation resumes the load and the demand peaks per year with a few values allowing the easy integration of several data sources as stochastic models, expert knowledge, or even fictional scenarios. This provides a planning tool useful to analyze several cases with low computational time. Presented discretization is not suitable for realtime operation due to the assumption of the ramp values in MW/h, when they usually are given in orders of MW/min. Last, solving the proposed ramping constrained model and comparing the solution with the exclusive economic model applied to four scenarios, a comparison of the performance is made. The comparison shows that considering the generation technology ramp velocities in the technology mix problem defines a new optimal solution cost. The new solution provides a reliable generation mix design able to fulfill the system power changes imposed by the load dynamics. Then, using the MPC methodology, the model is used to create an optimal generation investment plan. The MPC formulation provides a vector with the optimal solutions for all the scenarios in the order of seconds. The methodology deals with the uncertainties of the LDC forecast considering past and future values of the variables to make an technical and economical GEP solution. The MPC strategy found an optimal economic solution including the ramping constraints into the problem. Mixing the operation of slow and fast generation technologies, the strategy minimized the installed capacity changes, meeting the load requirements by keeping the economic cost close to the traditional solution.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This article is sponsored by COLCIENCIASConv. 528; KU Leuven, ELECTA Group; and Universidad Nacional de Colombia, GAUNAL Group.
References
 F. Ueckerdt, R. Brecha, G. Luderer et al., “Representing power sector variability and the integration of variable renewables in longterm energyeconomy models using residual load duration curves,” Energy, vol. 90, pp. 1799–1814, 2015. View at: Publisher Site  Google Scholar
 T. Boßmann and I. Staffell, “The shape of future electricity demand: exploring load curves in 2050s Germany and Britain,” Energy, vol. 90, pp. 1317–1333, 2015. View at: Publisher Site  Google Scholar
 K. Akbari, M. M. Nasiri, F. Jolai, and S. F. Ghaderi, “Optimal investment and unit sizing of distributed energy systems under uncertainty: a robust optimization approach,” Energy and Buildings, vol. 85, pp. 275–286, 2014. View at: Publisher Site  Google Scholar
 P. N. Georgiou, “A bottomup optimization model for the longterm energy planning of the Greek power supply sector integrating mainland and insular electric systems,” Computers and Operations Research, vol. 66, pp. 292–312, 2016. View at: Publisher Site  Google Scholar
 I. Staffell and R. Green, “Is there still merit in the merit order stack? the impact of dynamic constraints on optimal plant mix,” IEEE Transactions on Power Systems, vol. 31, no. 1, pp. 43–53, 2016. View at: Publisher Site  Google Scholar
 C. De Jonghe, E. Delarue, R. Belmans, and W. D'haeseleer, “Determining optimal electricity technology mix with high level of wind power penetration,” Applied Energy, vol. 88, no. 6, pp. 2231–2238, 2011. View at: Publisher Site  Google Scholar
 J. Valinejad and T. Barforoushi, “Generation expansion planning in electricity markets: a novel framework based on dynamic stochastic MPEC,” International Journal of Electrical Power and Energy Systems, vol. 70, pp. 108–117, 2015. View at: Publisher Site  Google Scholar
 A. van Stiphout, K. Poncelet, G. Deconinck, and K. De Vos, “The impact of operating reserves in generation expansion planning with high shares of renewable energy sources,” in Proceedings of the IAEE European Energy Conference, pp. 1–6, 2014. View at: Google Scholar
 A. J. Conejo, L. B. Morales, S. J. Kazempour, and A. S. Siddiqui, Investment in Electricity Generation and Transmission, Springer, New York, NY, USA, 2016. View at: Publisher Site
 H. D. Sherali, A. L. Soyster, F. H. Murphy, and S. Sen, “Linear programming based analysis of marginal cost pricing in electric utility capacity expansion,” European Journal of Operational Research, vol. 11, no. 4, pp. 349–360, 1982. View at: Publisher Site  Google Scholar  MathSciNet
 P. Tatjewski, “Advanced control and online process optimization in multilayer structures,” Annual Reviews in Control, vol. 32, no. 1, pp. 71–85, 2008. View at: Publisher Site  Google Scholar
 J. B. Rawlings, “Tutorial overview of model predictive control,” IEEE Control Systems Magazine, vol. 20, no. 3, pp. 38–52, 2000. View at: Publisher Site  Google Scholar
 K. De Vos, Sizing and allocation of operating reserves following wind power integration [Ph.D. thesis], KU Leuven, Departement Elektrotechniek (ESAT), 2013.
Copyright
Copyright © 2017 P. Deossa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.