Next Article in Journal
Global Health Perspectives on Cigarette Butts and the Environment
Next Article in Special Issue
Development of a WebGIS-Based Analysis Tool for Human Health Protection from the Impacts of Prescribed Fire Smoke in Southeastern USA
Previous Article in Journal
The Implication of Land-Use/Land-Cover Change for the Declining Soil Erosion Risk in the Three Gorges Reservoir Region, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Metamodeling Framework for Quantifying Health Damages of Power Grid Expansion Plans

1
Department of Supply Chain Management, Rutgers Business School, Newark, NJ 07102, USA
2
Department of Industrial & Systems Engineering, Rutgers University, Piscataway, NJ 07102, USA
3
Department of Industrial Engineering, Tsinghua University, 30 Shuangqing Rd, Haidian District, Beijing 10084, China
4
Center for Energy, Economic & Environmental Policy, Rutgers University, New Brunswick, NJ 07102, USA
5
Department of Chemistry, University of California-Irvine, Irvine, CA 92697, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(10), 1857; https://doi.org/10.3390/ijerph16101857
Submission received: 2 April 2019 / Revised: 19 May 2019 / Accepted: 21 May 2019 / Published: 26 May 2019
(This article belongs to the Special Issue Air Quality and Health Predictions)

Abstract

:
In this paper, we present an analytical framework to establish a closed-form relationship between electricity generation expansion planning decisions and the resulting negative health externalities. Typical electricity generation expansion planning models determine the optimal technology–capacity–investment strategy that minimizes total investment costs as well as fixed and variable operation and maintenance costs. However, the relationship between these long-term planning decisions and the associated health externalities is highly stochastic and nonlinear, and it is computationally expensive to evaluate. Thus, we developed a closed-form metamodel by executing computer-based experiments of a generation expansion planning model, and we analyzed the resulting model outputs in a United States Environmental Protection Agency (EPA) screening tool that approximates the associated human health externalities. Procedural guidance to verify the accuracy and to select key metamodel parameters to enhance its prediction capability is presented. Specifically, the metamodel presented in this paper can predict the resulting health damages of long-term power grid expansion decisions, thus, enabling researchers and policy makers to quickly assess the health implications of power grid expansion decisions with a high degree of certainty.

1. Introduction

Since 2001, electricity generation in the United States has followed an increasing trajectory as there has been increasing demand for uninterrupted access to electricity for various global industries that are vital to economic growth [1]. In order to satisfy this increasing demand reliably and economically, fossil fuels, such as coal and natural gas, are used as the primary sources of electricity, accounting for nearly two-thirds of the total energy supply [1]. However, since fossil fuels are burned to produce steam that is then used to power turbines that drive electrical generators, this produces air emissions such as CO2, SOX, NOX, volatile organic compounds (VOCs), and particulate matter [2]. As a result, ground-level ozone concentrations in the tropospheric region of the atmosphere increase as a result of photochemical reactions with emissions from fossil fuel plants. These by-products are health hazards, as they contribute to smog, which may ultimately lead to heart and chronic lung conditions like asthma, bronchitis, and emphysema [3]. Furthermore, air emissions from fossil fuels react with atmospheric water molecules to produce various forms of acidic precipitation and dry acidic deposition, known as acid rain. This occurs when sulfur dioxide and nitric oxides react with atmospheric water to generate sulfuric acid and nitric acid, which creates acid rain [4]. While acid rain does not pose any direct health hazards, fine acidic particulates in air contribute to heart and lung conditions, which result in health externalities or costs absorbed by consumers [4].
Based on the literature, one of the earliest approaches to quantifying health damages from the electricity sector was developed by Rowe et al. In this research, critical factors in computing externalities for electricity generation sources were assessed by utilizing the New York State Environmental Externalities Cost Study and computerized externality model (EXMOD). Using response surfaces in combination with EXMOD model outputs, the researchers were able to develop a model to predict health externalities from electricity generation by varying 15 different factors, including the selection of generation type, location, and operating characteristics [5]. Their research suggested that the most critical factors in developing a model framework were the selection and application of air dispersion models, selection of air pollution thresholds for health impacts, reduced life span risks associated with ozone exposure and long-term exposure to particulate matter, values for CO2 damages, and the value to be applied to increased risks of reduced life span for individuals age 65 or older. This research was extended by Thanh and Lefevre, where they applied an impact pathway approach (IPA) to relate the origin of environmental burdens to human health consequences. This relationship enabled the estimation of health damages of sulfur dioxide and particulate matter emissions from the electricity generation output of four power units using different fuels (lignite, oil, natural gas, and coal) at four locations in Thailand. Using response surfaces, their results suggested that these externalities were relatively small, but not negligible, as they ranged from 0.006 to 0.05 U.S. cents per kilowatt-hour (in 1995 dollars) [6].
Conversely, researchers have also studied the health benefits associated with emissions reductions from the electricity sector. For example, Burtaw et al. quantified the economic benefits of reduced air pollution in the U.S. related to greenhouse gas emissions mitigation policies in the electricity sector. In this work, researchers leveraged historical emissions data and a United States Environmental Protection Agency (EPA) screening tool to generate regression models of human health benefits. By conducting hourly, day-ahead electricity market simulations, and linking this output to the EPA screening tool, their research suggested that a tax of $25 per metric ton of carbon emissions would yield NOX-related health benefits of about $8 per metric ton of carbon reduced in the year 2010 (1997 dollars) [7].
While the electricity sector is the prime culprit of air emissions in the United States, researchers have also worked to study the broader impact of air quality on human health across multiple sectors including transportation and agriculture. As part of these efforts, Voorhees et al. conducted a sensitivity analysis using an industry standard air quality model, BenMAP, and analyzed the associated human health externalities. BenMAP, which stands for the Environmental Benefits Mapping and Analysis Program, is an EPA economic model used for estimating health effects and the associated externalities associated with changes in air quality. BenMAP achieves this by leveraging a geographic information system-based program to estimate population-level exposure rates and changes in incidences of health outcomes [8]. In the work done by Voorhees et al., they leveraged the BenMAP tool to model various climate change scenarios, and they quantitatively assessed the output [9]. Furthermore, the management consulting firm ICF International also used BenMAP as a resource to create a damage function, which was a response surface used to evaluate changes in health indices as a result of the emissions reduction goals proposed in the Clean Power Plan [10].
State-of-the-art literature related to this particular domain primarily focuses on quantifying the impact of power generation, primarily from fossil fuel sources, on air quality as well as the aggregate societal benefits of avoided emissions resulting from the integration of renewables into the generation portfolio. For instance, in China, Qin et al. found that replacing coal with synthetic natural gas coupled with carbon capture and storage technology in the residential sector would avoid approximately 32,000 pollution-related premature deaths by 2020 [11]. Huang et al. qualitatively studied the relationship between fuel consumption for electricity and the resulting contributions to ambient concentrations of major air pollutants from the power sector [12]. They extended their work quantitatively by studying five development scenarios in an air quality model [13]. Additionally, Peng et al. evaluated various decarbonization scenarios for the electric power grid in order to develop pathways to mitigate peaking emissions by 2030 in China while simultaneously maximizing air quality and health benefits [14]. More aligned to the research in this paper, Millstein et al. estimated the benefits associated with reduced air emissions resulting from increased wind and solar generation using various EPA screening tools [15]. Additionally, Kerl et al. introduced the air pollutant optimization model (APOM), which solved the hourly unit commitment model in the state of Georgia with the assistance of a reduced form air quality model [16]. To estimate their health impact costs, a baseline hourly unit commitment model, along with the inputs and outputs of two EPA models, were synchronized linearly to develop response surfaces, or regression curves, that predicted health impact costs as a function of decision variables in the unit commitment model. In the long-term planning space, Rodgers et al. studied the health benefits of adding emissions limits to multi-period, multi-region, generation expansion planning (GEP) problems using an EPA screening tool [3]. Furthermore, Rodgers et al. extended this research by applying a simulation-based optimization framework to solve the GEP problem that simultaneously minimized market costs and health damages from air emissions in the objective function. To include health damages in their model, a metamodel was developed to approximate health damages from NOX and SO2 emissions associated with power grid capacity expansion decisions [17]. Specifically, health damages from simulated power grid capacity expansion plans were evaluated in the co-benefits risk assessment (COBRA) model. COBRA is an EPA screening model that helps state and local governments assess human health damages from air emissions at the county, state, regional, or national level [18]. After executing these computer-based experiments, a closed-form kriging metamodel that approximated health damages as a function of capacity expansion decisions was developed, which was used as a surrogate for health damages. While this research provides guidance on how to obtain an expansion plan under these circumstances, there is need for systematic framework to evaluate a metamodel and simultaneously select its parameters.
Drawing from the aforementioned state-of-the-art research, and because of the highly stochastic, nonlinear relationship between electricity generation and the resulting health damages, researchers often apply black box modeling techniques in these instances with little understanding of how decision variables impact the response in question [17]. This lack of system knowledge is intensified as the magnitude and scope of power systems planning decisions grows from hourly dispatch decisions to long-term capacity expansion decisions across multiple regions. In this paper, our primary objective is to close this gap by developing a step-by-step, analytical framework that accomplishes the following:
  • Leverages computer-based simulations of a generation expansion planning model, where the model outputs are used as inputs into the COBRA tool to quantify the resulting health damages;
  • Establishes a closed-form expression to predict the aforementioned health damages as a function of power grid expansion decisions using a kriging metamodel; and
  • Provides a procedure to select metamodel parameters that maximize prediction accuracy.

2. Materials and Methods

To address the aforementioned challenges associated with approximating health damages from generation expansion planning decisions, we extended the work done by Rodgers et al. in [3,17]. Specifically, we proposed a detailed, systematic framework to establish a mathematical relationship that quantified health damages as a function of power grid expansion decisions while simultaneously considering prediction accuracy. Following the approach developed by Rodgers et al., we first formulated and solved a GEP model to determine the optimal technology capacity investment strategy that minimized market costs, including capital costs, fixed operation and maintenance costs, and variable costs of generation. Using this GEP model, we simulated key model parameters from a normal distribution, and we executed multiple computer-based experiments from the optimization model. We then evaluated the outputs of our simulations in the EPA’s COBRA tool.. Next, we applied a kriging interpolation model to our dataset to express the relationship between the covariates—electricity expansion decisions—and the response—human health damages. This kriging interpolation model, which is a geospatial estimation method that predicts the values of a random field at unobserved locations based on an interpolated function of observed samples, served as a metamodel or surrogate function that approximated the output of COBRA using input data generated from GEP simulations, thus, it simplified the ability to quantify human health damages resulting from power grid expansion plans. Since our metamodel was a function of computer-based experiments, there was a quantifiable amount of error in prediction. Lastly, to account for this, we proposed a cross-validation procedure to select the metamodel with the least prediction error.
In order to accomplish our research objectives, our first task was to establish a mathematical relationship to predict human health damages associated with electric power grid capacity expansion decisions, per the guidance of Rodgers et al., as given in Figure 1 [17].
To summarize, we generated expansion planning scenarios to be evaluated in an EPA screening tool, COBRA, to quantify human health externalities. A kriging metamodel was then applied to this data set, which predicted human health damages as a function of electricity dispatch decisions.

2.1. Simulation and Experimentation Procedure

To create a metamodel of human health damages in the context of GEP, we needed to generate a diverse sample space of data points from our expansion planning model to be evaluated in the COBRA tool. We referred to this process as seed generation.
To initiate the seed generation process, as shown in Figure 2, first we solved a baseline GEP model. We then solved several experimental GEP models where key parameters were simulated from a normal distribution. For each experimental trial, we then compared the expansion decisions and obtained a change in MWh for each unit available within the system. Simultaneously, for each trial we obtained the percent change in NOX and SO2 emissions from the baseline, which were used as inputs to be evaluated in the COBRA model.

Generation Expansion Planning (GEP) Model Formulation

We applied the GEP model used by Rodgers et al. to conduct our GEP simulation trials. The objective function of this GEP model was the sum of total market costs, which included investment costs, fixed operation and maintenance costs, and variable operation and maintenance costs, which were discounted over a predetermined time horizon [17]. The elements of the objective function are given below.
Investment Cost:
O I = y Y 1 ( 1 + r ) y r 1 R i I α y i y y r 1 i
Investment costs are the capital costs associated with capacity expansion. y y r 1 i is the capacity expansion decision, expressed as a continuous decision variable and given in MW, by generation unit type i in year y in region r 1 . α y i is the corresponding investment cost ($/MW). Additionally, we defined r as the interest or discount rate, Y as the set of all years in the planning horizon, R as the set of all regions in the network, and I as the set of all generating unit options.
Fixed operation and maintenance costs:
O M = y Y 1 ( 1 + r ) y [ r 1 R i I ( u ( u y ) Y β u i y u r 1 i ) + r 1 R i I * β y i q r 1 i ]
Fixed operation and maintenance costs for unit i in year y are given by β y i in units of $/MW. In Equation (2), I* is the set of all existing generating unit options, I is the set of all available generating unit options, and q r 1 i is the initial capacity of unit i in region r1. Additionally, we used the alias u as a surrogate notation for the year, y , to calculate the fixed operation and maintenance costs for all new capacity investments up to year y .
Variable operation and maintenance costs:
O G = y Y 1 ( 1 + r ) y t T r 1 R i I v y t i x y t r 1 i
Variable operation and maintenance costs (including fuel costs) of generation from unit i during period t of year y are given by v y t i in units of $/MWh. Each year is subdivided into six distinct time periods consisting of: (i) spring/fall-offpeak; (ii) spring/fall-peak; (iii) summer-offpeak; (iv) summer-peak; (v) winter-offpeak; and (vi) winter-peak, which span the set of all time intervals, T. x y t r 1 i is the amount of aggregate electricity dispatched (MWh) from unit i in region r 1 in period t of year y .
Steam generation revenue:
O R = y Y 1 ( 1 + r ) y t T r 1 R i I R E V ϕ y i x y t r 1 i
Annual steam generation (or cogeneration) revenue from combined heat and power (CHP) units is given by ϕ y i in terms of $/MWh, and x y t r 1 i is the amount of electricity dispatched from cogeneration units (in MWh) from in region r 1 period t of year y .
The total cost objective function, z , is given by Equation (5).
z = O G + O I + O M O R .
In our models, we minimized z across the network, subjected to system constraints. Definitions for GEP model sets, parameters, and decision variables are given in Table 1, Table 2 and Table 3.
The generation expansion planning formulation is given as follows:
min   z = O G + O I + O M O R ,
and is subject to:
r 1 Ω ( y , r 1 , r 2 ) l t w y t r 1 r 2 r 1 Ω ( y , r 1 , r 2 ) w y t r 2 r 1 + i I x y t r 1 i = d y t r 1         y Y , t T , r 2 R
x y t r 1 i ( q r 1 i + u ( u y ) Y y u r 1 i ) δ t , i h t           y Y , t T , r 1 R , i I N D
x y t r 1 i ( q r 1 i + u ( u y ) Y y u r 1 i ) δ t i h t η t r 1 i         y Y , t T , r 1 R , i I N D
w y t r 1 r 2 w y t r 2 r 1 χ y r 1 r 2 h t         ( y , r 1 , r 2 ) Ω ( y , r 1 , r 2 ) , Y , t T , i I
w y t r 2 r 1 w y t r 1 r 2 χ y r 2 r 1 h t   ( y , r 1 , r 2 ) Ω ( y , r 1 , r 2 ) , Y , t T , i I
i I N D q r 1 i + u ( u y ) Y i I N D y u r 1 i + i I N D q r 1 i θ i + u ( u y ) Y i I N D y u r 1 i θ i d y i * m r 1       y Y , r 1 R
r 1 R y y r 1 i κ y , i         y Y , i I M A X
y Y r 1 R y y r 1 i κ i *         i I M A X
x y t r 1 i 0 , w y t r 1 r 2 0 , y y r 1 i 0           ( y , r 1 , r 2 ) Ω ( y , r 1 , r 2 ) , Y , t T , i I
Equation (6) is the demand (or energy balance) constraint. For each of the periods in a given year, the total generation and transmission from both existing and new generating units should be at least as much as the corresponding demand in that region. Equations (7) and (8) are related to generating unit capacities for dispatchable and nondispatchable units. Equations (9) and (10) are the transmission network constraints. Equation (11) is the reserve margin constraint, and Equations (12) and (13) are investment limit constraints on an annual basis and throughout the time horizon of the model, respectively. Equation (14) specifies non-negativity for generation, transmission, and investment decision variables.
For each execution of the GEP model, the following key outputs were obtained and were used downstream as inputs into the COBRA tool:
  • Generation (in MWh) in period t of year y of unit i in region r1;
  • NOX and SO2 emissions in region r1 in year y.
As previously stated, to create a metamodel of health damages as a function of power grid expansion decisions, a diverse sample space of data points (or outputs), via simulation, from our GEP model were evaluated in the COBRA tool. Following the procedure outlined in Figure 1, first we solved a baseline GEP model, which minimized total market costs. Concurrently, we solved additional experimental GEP models with the same model parameters as the baseline model but with simulated numerical values. Specifically, to obtain a diverse sample space, we varied several key parameters in the experimental model by simulating from a normal distribution with a coefficient of variation following a uniform random variable on the interval (0, 1). The parameters we varied in the GEP model were as follows:
  • Demand;
  • Peak demand;
  • Reserve margin;
  • Minimum and maximum limits;
  • Minimum renewable generation levels;
  • Emissions rates
  • Emissions limits;
  • Costs (investment, fixed, and variable);
  • Steam revenue;
  • Transmission capacity;
  • Capacity factors for nondispatchable units; and
  • Derating values.
For each GEP experimental trial, we computed the following values for every year in the planning horizon:
Δ x y r 1 i = t T ( x y t r 1 i ,     e x p e r i m e n t a l * x y t r 1 i ,     b a s e l i n e * )
Δ N O X y ,   r 1 ( % ) = N O X y ,     e x p e r i m e n t a l N O X y ,     b a s e l i n e N O X y ,     b a s e l i n e
Δ S O 2 y ,   r 1 ( % ) = S O 2 y ,     e x p e r i m e n t a l S O 2 y ,     b a s e l i n e S O 2 y ,     b a s e l i n e
Equation (15) defines the annual deviation in electricity dispatch of a given unit in a specified region of the experimental GEP model from the baseline model. Equations (16) and (17) quantify the annual percent deviation of NOX and SO2 emissions, respectively, in a given region from the experimental trial to the baseline trial. We used the metrics from Equations (16) and (17) as inputs to the COBRA model, which allowed users to generate scenarios by specifying the changes in NOX and SO2 emissions on an annual basis for the appropriate region. COBRA then estimated the corresponding changes in equivalent particulate matter (PM2.5) concentrations, using a source–receptor matrix, and calculated the estimated health damages based on embedded epidemiological and economic functions. Computer-based experiments in COBRA can be executed fairly quickly; however, the relationship between GEP decisions and the resulting health damages is highly stochastic and nonlinear, and thus is difficult to express in closed-form. As previously mentioned, black box methods are often utilized in these instances with little understanding of how decision variables impact the response in question. Using this computer-based experimentation approach, this resulted in a data set that was used for metamodeling, which established a closed-form relationship and closed this gap of information regarding how key variables impacted the response variable.

2.2. Kriging Metamodel Details

Kriging is an interpolation method that uses the observed data at all sample points to provide a statistical prediction of an unknown function by minimizing its mean squared error (MSE) [19]. To develop a kriging metamodel, the output of a deterministic computer experiment is treated as a realization from a stochastic process, which is then defined as the sum of a global trend function and a Gaussian stochastic process.
In a general sense, various metamodeling techniques have been used to solve complex optimization and decision-making problems. The simplest of these metamodel techniques leverages regression theory. For instance, by applying regression theory, experimental design, and feasible region partitioning approaches, Roux et al. investigated methods for obtaining more metamodels to approximate the value of an objective function [20]. Kalil et al. applied regression techniques along with factorial design in an industrial bioprocess optimization problem to maximize yield and productivity [21]. Quanhong and Calil applied regression methods to quantify the effect of liquid:solid ratio, NaCl concentration, and reaction time in order to maximize protein production from germinant pumpkin seeds [22]. They then used this regression model as a surrogate objective function in an optimization problem.
When the cost of evaluating the response of a system by experiments is computationally expensive, it is critical to use the most accurate metamodel for prediction purposes. In comparison to other metamodeling methods used in optimization, kriging provides the best prediction accuracy and serves as a more accurate and reliable surrogate objective function [23]. Furthermore, these models can be applied to surrogate systems to reduce the total cost of response evaluation [24]. Such metamodels have been deployed in black-box systems [24], metal-forming processes [25], aerodynamic design applications [26], and even multiobjective optimization models [27].
In our work, the kriging metamodel we utilized to quantify human health damages as a function of GEP expansion decisions was given by Equation (18) [17].
γ ^ y r 1 i ( Δ x y r 1 i ) = μ ( Δ x y r 1 i ) + k = 1 n ( Δ x y r 1 i ) ω k [ γ y i ( Δ x y r 1 i k ) μ ( Δ x y r 1 i k ) ]
γ ^ y r 1 i ( Δ x y r 1 i ) is the predicted value of human health damages in year y for unit i in region r1, which is a function of a vector of unobserved dispatch values, Δ x y r 1 i . Δ x y r 1 i k is a vector of observed dispatch deviations obtained from the experimental procedure, indexed by k. The observed value of Δ x y r 1 i k is given by γ ^ y r 1 i ( Δ x y r 1 i ) , and n ( Δ x y r 1 i ) is the number of nearest neighbors to consider in the model. Kriging weight, ω k , is derived from the residuals, and the means of the predicted response and observed response are given by μ ( Δ x y r 1 i ) and μ ( Δ x y r 1 i k ) , respectively.
Residuals, R ( Δ x y r 1 i ) , in a kriging model have a stationary mean and covariance, where R ( Δ x y r 1 i ) = γ y i ( Δ x y r 1 i k ) μ ( Δ x y r 1 i ) , with E [ R ( Δ x y r 1 i ) ] = 0 , and C o v { R ( Δ x y r 1 i ) , R ( Δ x y r 1 i + h ) } = C ( h ) for some lag h. In this case, C ( h ) = C ( 0 ) + S V ( h ) , where we defined C ( 0 ) as the sill and S V ( h ) as the semivariogram. For this particular application, we employed the Gaussian semivariogram in Equation (19).
S V ( h ) = ( s i l l n u g g e t ) × ( 1 e x p ( 3 h 2 r a n g e 2 ) ) + n u g g e t .
The Gaussian semivariogram function has two parameters that the user can specify—the sill and nugget. The sill is a stationary value of the semivariogram, where the slopes of the tangent lines become zero. The nugget is the value of the semivariogram function with zero lag. Additionally, the range is the lag value when the slope of the tangent lines becomes zero.
To select the parameters to be used in our metamodel, and to assess prediction accuracy, we applied a modified cross-validation procedure as follows:
  • Step 1: We randomly partitioned the data set into eleven (11) equal subsets, designating one subset as a testing set and the remaining ten as validation sets to fit metamodels.
  • Step 2: Using the validation sets identified in Step 1, kriging metamodels were fitted with a Gaussian variogram.
  • Step 3: For each validation set, the mean absolute prediction error (MAPE) was computed against the test set.
  • Step 4: All validation sets were aggregated into a single data set. This set was then used to build a full metamodel of the test set data and calculate the model error, as described in Step 3.
We repeated this procedure for using sill-to-nugget (STN) ratios ranging from 1 to 1.9 to fit our metamodels. We then plotted the MAPE as a function of the STN ratio, and we selected the STN value that minimized the MAPE for the full metamodel.

3. Results and Discussion

For the numerical example presented in this paper, we considered the northeastern United States, as shown in the representation in Figure 3 [3,17].
Our test network had six regions, given as follows:
  • NE: New England (Maine, Vermont, New Hampshire, Massachusetts, Connecticut, and Rhode Island);
  • NY: New York State (excluding New York City);
  • NYC: New York City;
  • NJ: New Jersey;
  • MD & DE: Maryland, Delaware, and the District of Columbia;
  • Rest of PJM: Illinois, Indiana, Kentucky, Michigan, North Carolina, Ohio, Tennessee, Virginia, and West Virginia.
Each region in the network had demand for electricity and also was capable of generating electricity. The available pool of generating unit technologies considered in our model comprised cycle gas turbines, coal, natural gas turbines, hydro, nuclear, petroleum, solar, biomass, on-shore wind, and off-shore wind. Additionally, connections between regions defined the transmission network. We considered a 25 y planning horizon (2015 through 2040). GEP model cost inputs are given in Table 4 (refer to Appendix A for a full list of GEP model assumptions and data sources).
Using the simulation procedure and the GEP model presented in Section 2, we executed ten simulation trials of our experimental GEP model, and we compared the results to a baseline GEP model to obtain changes in investment and dispatching decisions by unit and changes in NOX and SO2 emissions. All GEP models were formulated as linear programs and solved using a CPLEX solver (an optimization package that leverages the simplex algorithm via the C programming language) in the generalized algebraic modeling system (GAMS) [30]. We then input these emissions changes into the COBRA model to quantify the health damages for each trial. Once we generated our data set, we then fit our kriging metamodel that predicted health damages. Our kriging metamodel was created using XonGrid, which is an open source interpolation tool [31].
Applying the cross-validation procedure presented in Section 2, we obtained the following prediction errors as a function of STN ratios displayed in Table 5.
As evidenced by Table 5 and Figure 4, we executed the cross-validation procedure for various sill and nugget parameters in the Gaussian semivariogram. Based on the individual validation sets, on average, as we increased the STN-ratio, the prediction error decreased. Variability of the prediction error calculations was approximately 2%, which was relatively minor in the context of this problem. Upon applying the full data set to fit our metamodel, however, we noticed that the prediction error was significantly reduced and remained relatively stable for all values of k. Based on the full model, since the smallest error value occurred at STN = 1.5, this was the parameter we applied to the semivariogram and, thus, was used as an input parameter in our metamodel.
For comparative purposes, we displayed the results of the baseline model. Figure 5 displays the annual aggregate dispatch strategy associated with our baseline GEP model. In our network, base load units were nuclear, coal, hydro, and natural gas. Additional amounts of combined heat and power (or steam), solar, and wind (on land) were dispatched to satisfy the remainder of the demand. Table 6 gives the capacity investment strategy associated with our baseline GEP model, which suggested that we added over 17,000 MW of capacity to our network. The majority of these investments were in fossil fuel sources, as they were able to satisfy demand economically and reliably.
As given in the cost summary presented in Figure 6, market costs (specifically investment costs, fixed operation and maintenance costs, variable operation and maintenance costs, and combined heat and power (CHP) revenue) summed to roughly $465B, most of which were variable operation and maintenance costs. These costs were included in the objective function of our GEP and, thus, were optimized via the CPLEX tool. However, health damages of the resulting expansion plan were nearly three times the value of the optimized market costs. Additionally, the health damages predicted from our kriging metamodel were, globally, within 2.5% of the true values from COBRA.
Figure 7 and Table 7 further detailed the breakdown of health damage values by region from the output of our baseline GEP model. Overall, there were areas where the model either over-predicted or under-predicted, depending on the region and generating unit. Regionally, the model over-predicted by approximately 11% in the NYC region and 4% in the MD & DE region, which was driven by the over-prediction of coal and nuclear damages, respectively. Additionally, the model slightly under-predicted health damages in the NY region by approximately 2%, which was driven by the under-prediction of natural gas damages.

4. Conclusions

This paper presents an extension of the research done by Rodgers et al. in [3,17], where health damages from generation expansion planning decisions are approximated. Specifically, an analytical framework to establish a mathematical relationship that quantifies health damages as a function of power grid expansion decisions, while simultaneously considering prediction accuracy, is presented. As demonstrated by the results, our metamodel successfully demonstrates the ability to approximate health damages as a function of GEP decisions with approximately 2.5% error, globally. This allows researchers and policy makers to quantify health damages as a function of power grid expansion decisions using a closed-form function, and it enables them to make more informed decisions on expanding power grid capacity.
Considering the output of the baseline model presented in Section 3, if not included in the objective function, health damages are nearly triple the optimized market costs. With the metamodeling framework presented in this paper, researchers can incorporate health damages in the objective function of a cost minimization GEP with the certainty of a high degree of prediction accuracy. Ultimately, this may yield significant investments in cleaner sources of energy, such as wind and solar technology, thus reducing health damages.
Another potential research extension would allow for a more detailed assessment of air quality and health damages by systematically linking more sophisticated air quality and economic models used by the EPA to our generation expansion planning model. One of the primary drawbacks of the reduced form air quality model COBRA is that it is based on simplified functions that translate air emissions, such as NOX and SO2, into equivalent PM2.5 concentrations. Thus, in order to better estimate health damages, a more rigorous suite of screening tools can be used in place of the COBRA model to address this gap. Specifically, in place COBRA, we may substitute the EPA’s sparse matrix operator kernel emissions (SMOKE) model, which is a tool that allocates emissions both spatially and temporally. Using the SMOKE model output, we may use this data as inputs into the community multiscale air quality (CMAQ) model, which computes pollutant concentrations by using continuity equations [32]. Economic health implications of the associated pollutant concentrations can be assessed in the EPA’s BenMAP (Environmental Benefits Mapping and Analysis Program) tool, which approximates the health damages as a function of air quality effects [8]. Systematically linking these EPA models to our generation expansion planning model via the framework presented in this paper would not only yield more accurate emissions predictions, but it would also provide decision makers with a more detailed perspective of the health implications of expansion plans.

Author Contributions

Conceptualization, M.D.R., D.W.C., F.A.F. and A.G.C.; methodology, M.D.R., D.W.C., F.A.F. and A.G.C.; software, M.D.R.; validation, M.D.R., D.W.C. and F.A.F.; formal analysis, M.D.R.; investigation, M.D.R.; resources, M.D.R., D.W.C., F.A.F. and A.G.C.; data curation, M.D.R., D.W.C. and F.A.F.; writing-original draft preparation, M.D.R.; writing-review and editing, M.D.R., D.W.C. and F.A.F.; visualization, M.D.R.; supervision, D.W.C., F.A.F. and A.G.C.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Additional Model Assumptions

Figure A1. Load Profile by Region (U.S. Energy Information Administration, 2015).
Figure A1. Load Profile by Region (U.S. Energy Information Administration, 2015).
Ijerph 16 01857 g0a1
Table A1. Initial Capacity by Region (in MW) (Cambridge Energy Solutions, 2011).
Table A1. Initial Capacity by Region (in MW) (Cambridge Energy Solutions, 2011).
RegionNuclearCombined Cycle Gas TurbinesCoal (Low)Coal (High)Gas TurbinesWind (On Land)PetroleumHydroSolar
NE2873203241412,44721100365610,342303
NY52825313446327633051274995705-
NYC-22167765-5426-107--
NJ4119605512412063481684740613
MD & DE5943-34191340--5905348128
RoPJM95838041491117,143351570420422812
Totals27,79921,82824,21235,75616,1721985470324,082446
Table A2. Transmission Limits by Region (in MWh) (Cambridge Energy Solutions, 2011).
Table A2. Transmission Limits by Region (in MWh) (Cambridge Energy Solutions, 2011).
RegionReceiving Region
NENYNYCNJMD & DERoPJM
Transmitting RegionNE-1420----
NY---950--
NYC430--1685--
NJ------
MD & DE-----5150
RoPJM-2000 9268--
Table A3. Unit Derating and Capacity Value Percentages (Cambridge Energy Solutions, 2011) (U.S. Energy Information Administration, 2015).
Table A3. Unit Derating and Capacity Value Percentages (Cambridge Energy Solutions, 2011) (U.S. Energy Information Administration, 2015).
Unit TypeDerating Factor (%)Capacity Value (%)
Nuclear91%100%
Combined Cycle Gas Turbines84%100%
Gas Turbines86%100%
Wind (On Land)97%32%
Gas Turbines97%100%
Combined Heat & Power86%100%
Biomass89%100%
Coal (High)86%100%
Coal (Low)86%100%
Solar95%14%
Petroleum86%100%
Hydro95%100%
Table A4. Unit Capacity Factors by Region (Cambridge Energy Solutions, 2011) (U.S. Energy Information Administration, 2015).
Table A4. Unit Capacity Factors by Region (Cambridge Energy Solutions, 2011) (U.S. Energy Information Administration, 2015).
RegionWind (On Land)Wind (Off Shore)Solar
NE31%59%13%
NY31%57%12%
NYC31%57%14%
NJ21%39%14%
MD & DE25%39%14%
RoPJM27%43%13%
Table A5. Reserve Margin Percentages.
Table A5. Reserve Margin Percentages.
RegionReserve Margin Capacity (% of Peak Demand)
NE116%
NY117%
NYC117%
NJ115%
MD & DE115%
RoPJM115%
Table A6. Regional Greenhouse Gas Initiative (RGGI) CO2 Emissions Limits.
Table A6. Regional Greenhouse Gas Initiative (RGGI) CO2 Emissions Limits.
YearTotal Annual CO2 Emissions Limits in RGGI Regions (in lbs, Billions)
2015366.6
2016357.4
2017348.5
2018 through 2040 annually339.8
Table A7. Regional Emissions Limits (U.S. Environmental Protection Agency, 2013).
Table A7. Regional Emissions Limits (U.S. Environmental Protection Agency, 2013).
RegionNOX Emissions Limits (in lbs, Millions)SO2 Emissions Limits (in lbs, Millions)
NE120.000.54
NY120.000.26
NYC120.000.26
NJ120.000.13
MD & DE120.000.36
RoPJM120.001.07
Table A8. Annual Emissions Limits (U.S. Environmental Protection Agency, 2013).
Table A8. Annual Emissions Limits (U.S. Environmental Protection Agency, 2013).
YearTotal Annual SO2 Emissions Limits (in lbs, Billions)
2015 through 2040 annually17.9
Table A9. Transmission Losses.
Table A9. Transmission Losses.
PeriodTransmission Losses by Period (%)
Summer-Peak91%
Summer-Off-Peak91%
Winter-Peak95%
Winter-Off-Peak95%
Spring/Fall-Peak93%
Spring/Fall-Off-Peak93%
Table A10. Renewables Trading Network (General Electric International, 2014).
Table A10. Renewables Trading Network (General Electric International, 2014).
RegionReceiving Region
Transmitting RegionNENE
NYNY
NJ
MD & DE
RoPJM
NYCNYC
NE
NJNJ
RoPJM
MD & DEMD & DE
NJ
NY
RoPJM
RoPJMRoPJM
Table A11. Available Renewables by Region (General Electric International, 2014).
Table A11. Available Renewables by Region (General Electric International, 2014).
RegionRenewables Available by Region
MD & DESolar
Biomass
Wind (On Land)
Wind (Off Shore)
NEBiomass
Wind (On Land)
Wind (Off Shore)
Solar
NJBiomass
Wind (On Land)
Wind (Off Shore)
Solar
NYWind (On Land)
Wind (Off Shore)
Biomass
Solar
NYCBiomass
Wind (Off Shore)
Solar
RoPJMBiomass
Wind (On Land)
Solar
Table A12. Minimum Percentage of Total Annual Dispatch from Renewables by Region (General Electric International, 2014).
Table A12. Minimum Percentage of Total Annual Dispatch from Renewables by Region (General Electric International, 2014).
YearNENYNYCNJMD & DERoPJM
201511%6%6%12%12%11%
201612%6%6%13%14%14%
201713%6%6%14%15%14%
201814%6%6%16%17%15%
201915%6%6%18%18%15%
202016%6%6%20%19%16%
202116%6%6%23%19%18%
202216%6%6%23%20%18%
202316%6%6%23%20%18%
202417%6%6%23%20%18%
202517%6%6%23%20%18%
202617%6%6%23%20%18%
202717%6%6%23%20%18%
202817%6%6%23%20%18%
202917%6%6%23%20%18%
203017%6%6%23%20%18%
203117%6%6%23%20%18%
203217%6%6%23%20%18%
203317%6%6%23%20%18%
203417%6%6%23%20%18%
203517%6%6%23%20%18%
203617%6%6%23%20%18%
203717%6%6%23%20%18%
203817%6%6%23%20%18%
203917%6%6%23%20%18%
204017%6%6%23%20%18%
Table A13. Minimum Percentage of Total Regional Dispatch from Renewable Energy Sources (General Electric International, 2014).
Table A13. Minimum Percentage of Total Regional Dispatch from Renewable Energy Sources (General Electric International, 2014).
RegionBiomassWind (On Land)Solar
MD & DE--2%
NE---
NJ--2%
NY---
NYC---
RoPJM1%8%10%
Table A14. Annual Construction Limits by Unit (in MW per Region).
Table A14. Annual Construction Limits by Unit (in MW per Region).
YearWind (On Land)Combined Heat & PowerNuclearWind (Off Shore)SolarGas TurbinesCombined Cycle Gas TurbinesPetroleum
201510,00010,0000010,00010,00010,00010,000
201610,00010,0000010,00010,00010,00010,000
201710,00010,0000010,00010,00010,00010,000
201810,00010,0000010,00010,00010,00010,000
201910,00010,0000010,00010,00010,00010,000
202010,00010,0000010,00010,00010,00010,000
202110,00010,000010,00010,00010,00010,00010,000
202210,00010,000010,00010,00010,00010,00010,000
202310,00010,000010,00010,00010,00010,00010,000
202410,00010,000010,00010,00010,00010,00010,000
202510,00010,000010,00010,00010,00010,00010,000
202610,00010,000010,00010,00010,00010,00010,000
202710,00010,00010,00010,00010,00010,00010,00010,000
202810,00010,00010,00010,00010,00010,00010,00010,000
202910,00010,00010,00010,00010,00010,00010,00010,000
203010,00010,00010,00010,00010,00010,00010,00010,000
203110,00010,00010,00010,00010,00010,00010,00010,000
203210,00010,00010,00010,00010,00010,00010,00010,000
203310,00010,00010,00010,00010,00010,00010,00010,000
203410,00010,00010,00010,00010,00010,00010,00010,000
203510,00010,00010,00010,00010,00010,00010,00010,000
203610,00010,00010,00010,00010,00010,00010,00010,000
203710,00010,00010,00010,00010,00010,00010,00010,000
203810,00010,00010,00010,00010,00010,00010,00010,000
203910,00010,00010,00010,00010,00010,00010,00010,000
204010,00010,00010,00010,00010,00010,00010,00010,000
Table A15. Overall Construction Limits by Unit.
Table A15. Overall Construction Limits by Unit.
Unit TypeMaximum Construction Limit (MW)
Wind (On Land)150,000
Combined Heat & Power150,000
Nuclear150,000
Wind (Off Shore)150,000
Solar150,000
Petroleum150,000

References

  1. U.S. Energy Information Administration. Electricity Data Browser; U.S. Energy Information Administration: Washington, DC, USA, 2015. Available online: http://www.eia.gov/electricity/data/browser/ (accessed on 6 January 2016).
  2. Farkas, C.M.; Moeller, M.D.; Felder, F.A.; Baker, K.R.; Rodgers, M.D.; Carlton, A.G. Temporalization of Peak Electric Generation Particulate Matter Emissions during High Energy Demand Days. Environ. Sci. Technol. 2015, 49, 4696–4704. [Google Scholar] [CrossRef] [PubMed]
  3. Rodgers, M.D.; Coit, D.W.; Felder, F.; Carlton, A.G. Assessing the Effects of Power Grid Expansion on Human Health Externalities. Socio-Econ. Plan. Sci. 2019, 66, 92–104. [Google Scholar] [CrossRef]
  4. Union of Concerned Scientists. The Hidden Cost of Fossil Fuels; Union of Concerned Scientists: Cambridge, MA, USA, 2013; Available online: http://www.ucsusa.org/clean_energy/our-energy-choices/coal-and-other-fossil-fuels/the-hidden-cost-of-fossil.html#.VzkEo5GDGko (accessed on 1 March 2013).
  5. Rowe, R.D.; Lang, C.M.; Chestnut, L.G. Critical factors in computing externalities for electricity resources. Resour. Energy Econ. 1996, 18, 363–394. [Google Scholar] [CrossRef]
  6. Thanh, B.D.; Lefevre, T. Assessing health impacts of air pollution from electricity generation: The case of Thailand. Environ. Impact Assess. Rev. 2000, 20, 137–158. [Google Scholar] [CrossRef]
  7. Burtaw, D.; Krupnick, A.; Palmer, K.; Paul, A.; Toman, M.; Bloyd, C. Ancillary benefits of reduced air pollution in the US from moderate greenhouse gas mitigation policies in the electricity sector. J. Environ. Econ. Manag. 2003, 45, 650–673. [Google Scholar] [CrossRef]
  8. Abt Associates Inc. Environmental Benefits Mapping and Analysis Program User’s Manual; Abt Associates Inc.: Rockville, MD, USA, 2012. Available online: http://www.epa.gov/airquality/benmap/models/BenMAPManualOct2012.pdf (accessed on 1 March 2013).
  9. Voorhees, A.S.; Fann, N.; Fulcher, C.; Dolwick, P.; Hubbell, B.; Bierwagen, B.; Morefield, P. Climate Change-Related Temperature Impacts on Warm Season Heat Mortality: A Proof-of-Concept Methodology Using Ben MAP. Environ. Sci. Technol. 2011, 45, 1450–1457. [Google Scholar] [CrossRef] [PubMed]
  10. U.S. Environmental Protection Agency. Regulatory Impact Analysis for the Clean Power Plan Final Rule; U.S. Environmental Protection Agency: Washington, DC, USA, 2015.
  11. Qin, Y.; Wagner, F.; Scovronick, N.; Peng, W.; Yang, J.; Zhu, T.; Smith, K.R.; Mauzerall, D.L. Air quality, health, and climate implications of China’s synthetic natural gas development. Proc. Natl. Acad. Sci. USA 2017, 114, 4887–4892. [Google Scholar] [CrossRef] [PubMed]
  12. Huang, L.; Hu, J.; Chen, M.; Zhang, H. Impacts of power generation on air quality in China—part I: An overview. Resour. Conserv. Recycl. 2017, 121, 103–114. [Google Scholar] [CrossRef]
  13. Hu, J.; Huang, L.; Chen, M.; He, G.; Zhang, H. Impacts of power generation on air quality in China—Part II: Future scenarios. Resour. Conserv. Recycl. 2017, 121, 115–127. [Google Scholar] [CrossRef]
  14. Peng, W.; Yang, J.; Lu, X.; Mauzerall, D.L. Potential co-benefits of electrification for air quality, health, and CO2 mitigation in 2030 China. Appl. Energy 2018, 218, 511–519. [Google Scholar] [CrossRef]
  15. Millstein, D.; Wiser, R.; Bolinger, M.; Barbose, G. The climate and air-quality benefits of wind and solar power in the United States. Nat. Energy 2017, 2, 1–10. [Google Scholar] [CrossRef]
  16. Kerl, P.Y.; Zhang, W.; Moreno-Cruz, J.B.; Nenes, A.; Realff, M.J.; Russell, A.G.; Sokol, J.; Thomas, V.M. New approach for optimal electricity planning and dispatching with hourly time-scale air quality and health considerations. Proc. Natl. Acad. Sci. USA 2015, 112, 10884–10889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Rodgers, M.D.; Coit, D.W.; Felder, F.A.; Carlton, A. Generation expansion planning considering health and societal damages—A simulation-based optimization approach. Energy 2018, 164, 951–963. [Google Scholar] [CrossRef]
  18. United States Environmental Protection Agency. User’s Manual for the Co-Benefits Risk Assessment (COBRA) Screening Model—Version 2.7; United States Environmental Protection Agency: Washington, DC, USA, 2015.
  19. Bohling, G. Kriging. 19 October 2015. Available online: http://people.ku.edu/~gbohling/cpe940/Kriging.pdf (accessed on 22 June 2012).
  20. Roux, W.J.; Stander, N.; Haftka, R.T. Response Surface Approximations for Structural Optimization. Int. J. Numer. Methods Eng. 1998, 42, 517–534. [Google Scholar] [CrossRef]
  21. Kalil, S.J.; Maugeri, F.; Rodrigues, M.I. Response surface analysis and simulation as a tool for bioprocess design and optimization. Process Biochem. 2000, 35, 539–550. [Google Scholar] [CrossRef]
  22. Quanhong, L.; Caili, F. Application of response surface methodology for extraction optimization of germinant pumpkin seeds protein. Food Chem. 2005, 92, 701–706. [Google Scholar] [CrossRef]
  23. Han, Z.-H.; Zhang, K.-S. Surrogate-Based Optimization; Roeva, D.O., Ed.; IntechOpen: London, UK, 2012. [Google Scholar]
  24. Huang, D.; Allen, T.T.; Notz, W.I.; Zeng, N. Global Optimization of Stochastic Black-Box Systems via Sequential Kriging Meta-Models. J. Glob. Optim. 2006, 34, 441–466. [Google Scholar] [CrossRef]
  25. Jakumeit, J.; Herdy, M.; Nitsche, M. Parameter optimization of the sheet metal forming process using an iterative parallel Kriging algorithm. Struct. Multidiscip. Optim. 2005, 29, 498–507. [Google Scholar] [CrossRef]
  26. Jeong, S.; Murayama, M. Efficient Optimization Design Method Using Kriging Model. J. Aircr. 2005, 42, 413–420. [Google Scholar] [CrossRef]
  27. Li, M.; Li, G.; Azarm, S. A Kriging Metamodel Assisted Multi-Objective Genetic Algorithm for Design Optimization. J. Mech. Des. 2008, 130, 031401. [Google Scholar] [CrossRef]
  28. Nuclear Energy Institute. U.S. Capacity Factors by Fuel Type; Nuclear Energy Institute: Washington, DC, USA, 2012. [Google Scholar]
  29. U.S. Environmental Protection Agency. Energy and the Environment; U.S. Environmental Protection Agency: Washington, DC, USA, 2013.
  30. GAMS. CLPEX 12. GAMS. Available online: https://www.gams.com/latest/docs/S_CPLEX.html (accessed on 7 May 2019).
  31. XonGrid. XonGrid Interpolation Add-In. XonGrid. Available online: http://xongrid.sourceforge.net/ (accessed on 7 May 2019).
  32. The Institute for the Environment—The University of North Carolina at Chapel Hill. Smoke v3.1 User’s Manual; The Institute for the Environment: Chapel Hill, NC, USA, 2012; Available online: http://www.smoke-model.org/version3.1/html/ (accessed on 1 March 2013).
Figure 1. Metamodeling Procedure Flowchart. GEP: generation expansion planning; COBRA: co-benefits risk assessment.
Figure 1. Metamodeling Procedure Flowchart. GEP: generation expansion planning; COBRA: co-benefits risk assessment.
Ijerph 16 01857 g001
Figure 2. Simulation and Experimentation Procedure.
Figure 2. Simulation and Experimentation Procedure.
Ijerph 16 01857 g002
Figure 3. Transmission Network of the northeastern United States.
Figure 3. Transmission Network of the northeastern United States.
Ijerph 16 01857 g003
Figure 4. Mean Absolute Prediction Error Values from the Cross-Validation Procedure.
Figure 4. Mean Absolute Prediction Error Values from the Cross-Validation Procedure.
Ijerph 16 01857 g004
Figure 5. Baseline Model Dispatch Summary.
Figure 5. Baseline Model Dispatch Summary.
Ijerph 16 01857 g005
Figure 6. Total Cost Output Summary.
Figure 6. Total Cost Output Summary.
Ijerph 16 01857 g006
Figure 7. Health Damages by Region.
Figure 7. Health Damages by Region.
Ijerph 16 01857 g007
Table 1. GEP Model Set Definitions.
Table 1. GEP Model Set Definitions.
YYears in the planning horizon (2015 through 2040)—indexed by y (or u as an alias)
TPeriods in the planning horizon (summer—peak/offpeak, winter—peak/offpeak, and spring/fall—peak/offpeak)—indexed by t
IAll generating units (nuclear, combined-cycle gas turbine, natural gas turbine, wind (land and offshore), biomass, coal, combined heat and power, solar, petroleum, and hydro)
I*All existing units within the network—indexed by i
IMAXUnits with construction limits—indexed by i
IREVUnits with steam revenue (combined heat and power only)—indexed by i
INDNondispatchable units (wind and solar)—indexed by i
RRegions within the northeastern US network (New England, NY, NYC, NJ, MD & DE, and Rest of PJM)—indexed by r 1 (or r 2 )
Ω ( y , r 1 , r 2 ) Set of transmission lines by year—indexed by ( y , r 1 , r 2 )
Φ ( r 1 , r 2 ) Network of regions for renewable trading—indexed by ( r 1 , r 2 )
Table 2. GEP Model Parameter Definitions
Table 2. GEP Model Parameter Definitions
η t r 1 i Capacity factor for nondispatchable units l t Transmission losses by season
θ i Capacity value by unit m r 1 Reserve margin by region
κ y , i Construction limits for each unit by yearrInterest rate (3%)
d y t r 1 Demand by period and region d y i * Peak demand by region and year
δ t i Derating value by season and unit ϕ y i Steam revenue ($/MWh)
β y i Fixed cost ($/MW) κ i * Total construction limits by unit
h t Hours in each period χ y r 1 r 2 Transmission capacity by year
q r 1 i Initial capacity by unit and region v y t i Variable costs ($/MWh)
α y i Investment costs ($/MW)
Table 3. GEP Model Decision Variables.
Table 3. GEP Model Decision Variables.
x y t r 1 i Generation in period t of year y of unit i in region r 1 (MWh)
y y r 1 i Capacity investment in year y in region r 1 of unit i (MW)
w y t r 1 r 2 Transmission of electricity in period t of year y from region r 1 to region r 2 (MWh)
Table 4. Generating Units and Their Costs, Capacities, and Emissions Rates.
Table 4. Generating Units and Their Costs, Capacities, and Emissions Rates.
Unit Type1 Aggregate Capacity (MW)2 Investment Costs ($/MW)2 Maintenance Costs ($/MW)2 Variable (Including Fuel) Costs ($/MWh)2 Fuel Costs ($/MWh)2 Steam Revenue ($/MWh)3 NOX (lbs/MWh)3 SO2 (lbs/MWh)
Biomass0$4,359,382$74,627$30.29$0.00$0.00410
Combined Cycle Gas Turbine (CCGT)21,828$1,096,358$13,546$69.37$15.69$0.0020
Combined Heat & Power (CHP)0$1,744,772$18,563$75.86$14.30$56.6320
Coal59,968$2,713,713$38,327$25.75$12.91$0.00613
Natural Gas Turbine (GT)16,172$733,926$12,198$50.07$15.69$0.0020
Hydro24,082$2,909,653$12,198$0.00$0.00$0.0000
Nuclear27,799$4,326,537$104,245$6.35$5.21$0.0000
Petroleum4703$1,114,480$14,452$99.31$64.87$0.00412
Solar446$6,989,283$13,523$0.00$0.00$0.0000
Off-Shore Wind0$4,459,051$98,446$0.00$0.00$0.0000
Wind on Land1985$2,226,694$35,088$0.00$0.00$0.0000
Notes: 1 [28]; 2 [1]; 3 [29].
Table 5. Mean Absolute Prediction Error (MAPE) Values.
Table 5. Mean Absolute Prediction Error (MAPE) Values.
SetSTN Values = Sill/Nugget
11.11.21.31.41.51.61.71.81.9
MAPE ValuesSet 115.8%15.7%14.8%13.6%14.5%14.3%14.6%15.5%12.2%12.4%
Set 217.1%16.1%14.3%15.8%12.9%12.7%12.8%13.3%14.0%15.4%
Set 313.7%14.1%15.9%11.7%12.0%13.6%14.0%12.8%12.8%13.6%
Set 414.1%12.5%14.0%11.9%11.4%12.5%12.2%14.5%12.7%12.4%
Set 519.1%19.2%21.0%20.0%22.6%21.8%20.8%19.8%17.4%16.9%
Set 618.1%19.6%19.3%15.9%16.8%16.7%15.2%13.6%15.2%14.1%
Set 720.0%18.1%16.7%16.0%18.9%15.0%15.8%16.8%15.5%17.2%
Set 814.6%14.7%12.7%12.0%12.9%13.3%13.3%12.3%12.2%11.2%
Set 917.0%14.3%15.2%15.7%16.9%16.1%15.2%15.7%16.3%14.7%
Set 1015.8%14.8%14.2%15.2%13.6%13.0%14.9%13.3%14.2%13.5%
10 Set Avg16.5%15.9%15.8%14.8%15.3%14.9%14.9%14.8%14.2%14.1%
Full Model2.7%1.8%1.8%1.9%1.8%1.7%1.9%1.9%1.8%2.0%
Table 6. Baseline Model Investment Strategy.
Table 6. Baseline Model Investment Strategy.
RegionInvestment Summary by Region (MW)
GTCoalCHPPetroleumCCGTHydroTotals
MD & DE23484046487097177385564
NJ2379144689269864906065
NYC21691420103379559206009
Total6897327025742202195773817,637
Table 7. Comparison of Health Damage Values by Unit and Region.
Table 7. Comparison of Health Damage Values by Unit and Region.
RegionNuclearGTCoal
Health Damages (COBRA) (in 2015 $, Billions)Health Damages (Metamodel) (in 2015 $, Billions)Error (%)Health Damages (COBRA) (in 2015 $, Billions)Health Damages (Metamodel) (in 2015 $, Billions)Error (%)Health Damages (COBRA) (in 2015 $, Billions)Health Damages (Metamodel) (in 2015 $, Billions)Error (%)
NE$22.2$23.13.9%$18.5$16.6−11.0%$163.3$162.3−0.6%
NY$40.8$36.2−12.6%$1.3$1.1−11.1%$117.2$118.20.9%
NYC$0.0$0.0$25.0$26.03.8%$203.7$228.610.9%
NJ$32.0$31.9−0.4%$45.8$42.6−7.4%$52.6$57.07.6%
MD & DE$46.1$49.26.2%$6.6$6.3−5.7%$89.0$92.23.4%
RoPJM$72.7$73.30.8%$30.5$28.8−5.7%$206.0$208.51.2%
Totals$213.8$213.7−0.1%$127.6$121.4−5.1%$831.9$866.74.0%

Share and Cite

MDPI and ACS Style

Rodgers, M.D.; Coit, D.W.; Felder, F.A.; Carlton, A.G. A Metamodeling Framework for Quantifying Health Damages of Power Grid Expansion Plans. Int. J. Environ. Res. Public Health 2019, 16, 1857. https://doi.org/10.3390/ijerph16101857

AMA Style

Rodgers MD, Coit DW, Felder FA, Carlton AG. A Metamodeling Framework for Quantifying Health Damages of Power Grid Expansion Plans. International Journal of Environmental Research and Public Health. 2019; 16(10):1857. https://doi.org/10.3390/ijerph16101857

Chicago/Turabian Style

Rodgers, Mark D., David W. Coit, Frank A. Felder, and Annmarie G. Carlton. 2019. "A Metamodeling Framework for Quantifying Health Damages of Power Grid Expansion Plans" International Journal of Environmental Research and Public Health 16, no. 10: 1857. https://doi.org/10.3390/ijerph16101857

APA Style

Rodgers, M. D., Coit, D. W., Felder, F. A., & Carlton, A. G. (2019). A Metamodeling Framework for Quantifying Health Damages of Power Grid Expansion Plans. International Journal of Environmental Research and Public Health, 16(10), 1857. https://doi.org/10.3390/ijerph16101857

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