Johns Hopkins Bloomberg School of Public Health
Quantify the Impacts of NOx Cap on Generators’ Emission Rate using Multiple Level Models
Abstract
In this paper, we applied the multiple-level model (MLM) to quantify the degree of stringency to which the US Environmental Protection Agency (EPA) NOx cap-and-trade imposes on the electricity generators. The US EPA NOx is implemented throughout the summer season from May to September every year. The design of the program is resemble to the crossover design of the case-control studies in epidemiology field. The strength of this design is that the non-summer season serves as the self-control groups where the time-invariant latent variables at generators’ level are automatically controlled. The hierarchical structure of our model includes generators at the lowest level, ownership at the second level and finally plant location at the upper level. The covariates that we control for are capacity, vintage year, output level and the fuel type at generators’ level, and price at location level. Using the MLM provides a more flexible modeling framework and allows us to decompose the residual variance which otherwise would not be explained by non-MLM approaches into variance at the various levels.
The application is to Pennsylvania – New Jersey – Maryland Interconnection market. The results show that significant variance reduction is achieved by including random-effects at three levels. Together with information about permits prices, we are able to estimate the effects of the regulatory programs on power industries in the monetary terms.
Introduction
Various studies have shown that the increases of mortality rate and other health-related endpoints are associated with the elevation levels of ozone concentration during summer time [1-5]. For example, Bell et al. [5] found that a 10-ppb increase in the previous week’s ozone was associated with a 0.52% increase in daily mortality and a 0.64% increase in cardiovascular and respiratory mortality, respectively. The formation of ozone is a complicated photochemical process, which involves the reactions of two precursors – NOx and volatile organics chemicals (VOCs) – in the presence of sunlight. The peak ozone concentration usually occurs during the early afternoon, i.e., 14:00 – 16:00, in summer time when adequate amount of biological VOCs reacts with vehicles-exhausted NOx in the morning traffic. The well-known nonlinearity in the isopleths diagram indicates whether the reduction of NOx or VOCs is an effective means to reduce ozone depends on the ratio of the two precursors [6]. This poses substantial challenges for the regulatory agencies to implement control strategies. The Mid-Atlantic region is one of the areas that are frequently in violation of ozone National Ambient Air Quality Standards (NAAQS) according to the US Environmental Protection Agency (EPA) [7]. Interestingly, unlike regulating SO2, which is led by the federal government, the establishment of ozone transportation commission (OTC) in 1997 is a joint effort by industries, environmental groups and the authorities at the state level [8]. Because the contribution of ozone precursors is mainly originated from the upwind states, notably power plants located in Ohio and Kentucky, where there is less incentive for them to reduce emissions since they are constantly in compliance with the NAAQS. The federal government seemingly should be in a great position to take the leadership.
The stewardship of the OTC led to the creation of regional cap-and-trade program for NOx in 1999 (OTC NOx Budget Program). The theoretical rationale behind the cap-and-trade programs is that the successful implementations will result in a least-cost solution for a given emission reduction target. While the geographical scope of the OTC program has evolved over time to encompass more states, the central piece of program remains intact.[1] First, the program is only enforced in the summer period from May 1st to September 30th. The program first gives out tradable permits to the affected facilities based on the historical heat input (three-year moving average) multiplied with a predetermined input emission rate.[2] One unit of tradable permit entitles the owners to emit one ton of NOx that year. Besides, the flexibility of the program allows the owners of permits freely trade them in the secondary markets; and the excess permits can be banked for future use but are subject to certain restrictions, which limit the upper bond of banking amount. Overall, the program has demonstrated great success over years.
One challenge faced by the policymakers is to know how effective the program is and to quantify the degree to which the introduced program forces the polluters to reduce their emission. Apparently, the aggregate emission cannot exceed the emission cap; but under what circumstances will a firm reduce their emission? In general, the operation of generators depends on several factors: wholesale price level, market structure, engineering constraints and etc. In particular, if marginal cost of a generator is below predicted wholesale price, the generator will be brought online given that transmission network is feasible. A generator may be operated differently if market is less competitive when demand is inelastic. Moreover, operating generators is subject to engineering constraints such as ramp rate, which limits how fast a generator can increase its output of given time. Therefore, the operating of generators is a function of where it is located, what kind of the technology it uses and what the objective the firm has.
To answer the aforementioned question turns out to be a daunting task, in part, because there are many ways an industry can response to the policy changes. In the short run, firms can modify their production factors such as labors, materials or the production processes; in the long run, it can engage in the research for advanced control technologies. In addition, the market structure and the nature of regulatory regime also have significant influence on it. For example, a monopoly under cost-of-service pricing scheme has little incentive to improve its performance since the incurred costs can be fully passed through on consumers [10][11]. Therefore, a full analysis requires considering the factors that affect generators’ operations at the market’s, the ownership’s and the generators’ level.
In compared with structural analysis, we take an empirical approach and model the resulting changes of generators’ performance as a consequence of the interplay of various factors at the different levels. The report is organized as following. In the section two, we will introduce the analytical framework of multi-level models. This includes defining the dependent variables and the covariates at each level. The data sources that we use for this analysis will be briefly discussed in the third section. We present our preliminary results in the section four, and enclose this report with a discussion.
Analytical Framework
To simplify the analysis, we consider covariates at the three levels[3] including a generator’s location or regional market, type of the ownership and the generation technology, which are indicated by l, f and g, respectively. Following the convention, the alphabetical letters are used for donating variables while the Greek variables are used for presenting coefficients. The capital letters indicate the time-invariant variables. The input NOx emission rates are assumed to have a normal distribution with mean [pic] and variance [pic]. We begin with the simplest multi-level model in which the only source of random effects is from the intercepts related to a generator’s location, type of the ownership and the generation technology. (The fixed-effect model can be obtained by ignoring the random-effect coefficients: [pic], [pic] and [pic].)
[pic] (1)
[pic] (2)
[pic] (3)
[pic] (4)
[pic] (5)
[pic] (6)
The dependent variable [pic] is the hourly input NOx emission rate in lbs/mmBTU of g-type generator owned by the f-type of ownership located at the location l during hour t. The nonlinear relationship between heat rate and output level is well known by power engineers. This is in particular the case for the peaking units such as combined-cycle and combustion turbines. However, it is unclear whether the nonlinearity is apparently for the relationship between the emission rate and output. As a beginning point, the model is analyzed with the linearity assumption[4]. To illustrate it, in Figure 1 we plot the NOx emission rate vs. output level of two generators, representing steam and combined-cycled technology. The first plot shows that steam unit’s input emission rate is scattered between 0.1 and 0.4 lbs/mmBTU, while in the second plot, the emission rate of combined-cycled unit is centered around 0.1 lbs/mmBTU. In both cases, there is no apparently nonlinearity existed. The distribution of the input NOx emission rate and output is displayed in Figures 2-3. A rather bi-mode distribution is shown in the steam unit and a single mode is with the combined-cycled unit. The possible explanation is the ramping-up process for steam unit has different effects on NOx emission rate. Figure 3 indicates that a significant amount of time steam unit is online at its capacity, while combined-cycled is operating only half of the time. The variable [pic] is an indicator variable, in which its value is equal to 1 when t is under ozone period, and 0, otherwise. The variable [pic] is the power price in $/MWh at the location l during t period; the variable [pic] indicates the output level in MW by the generators located at l with f-type of ownership using g-type technology in the period [pic]; the variable [pic] donates the rated capacity in MW of the generator located at l with f-type ownership using g-type technology; the variable [pic] is the vintage year of the generator located at l with f-type of ownership using g-type technology. The variable [pic]is an indicator with value equals to 1 for utility and 0 otherwise. The location variables [pic] are also indicator variables that are used to represent the ten distinct locations in the network. Finally, the variable [pic] is the error term and we assume it follows the AR(1) process (equation (3)). The random-effect coefficients ([pic], [pic] and [pic]) which are used to capture the variation of the latent variables are assumed to have normal distributions with mean 0 and variance [pic], [pic] and [pic], respectively. The coefficient [pic] is of our interest and it summaries the extra reduction in hourly NOx emission rate (if negative) in the logarithm scale.
Data Sources
There are two primary data sources supporting this analysis. The generation characteristics are obtained from Energy Information Administration (EIA). This includes an inventory of generating units and their NOx input emission rate, vintage year, generation capacity, generation technology and type of ownership. The information is then cross-validated with the information from the EPA’s Emissions and Generation Resource Integrated Database (eGRID)[5]. The data retrieved from the US EPA Continuous Emission Monitoring system (CEMS) provides the hourly measurements of NOx emission rate and generating output level. Of the total 265 emission sources listed in the CEMS in 2003 for Pennsylvania – New Jersey – Maryland (PJM) Interconnection, we identified 134 generating units, which have one year worth data. (Some emission sources are industrial boilers or oil refineries, which are not considered in the analysis.) The set of generating units accounts for 43% of total generation capacity but more than 70% of fossil-fueled capacity.
The PJM is spatially comprised with the ten zones/nodes: METED, PPL, BGE, PEPCO, JCPL, DPL, PECO, PSEG, AECO, and PENELEC. The difference in nodal power prices reflects the network topology and the locations of the generation capacity. In general, the power prices during the ozone period are roughly 30% to 40% higher than that during non-ozone period. This indicates fewer peaking units are operated during non-ozone period; and this could be problematic for our analysis because only a smaller portion of the generating units operating frequently in both periods. To account for this, we selected two types of generators: steam turbines and combined-cycle generators. The steam-turbine units are baseload units, which are operated through out the year; and combined-cycle generators are intermediate units, which are operated more frequently than other non-steam units. This left us a total of 88 generating units. Due to the large dataset, we decided to select the hour of 14:00 as the sample of our daily measurement and further sampled a subset of dataset as our sample.[6] This reduces our sample size to a manageable scale. The other approaches such as taking daily average will lead to the substantial loss of information and statistical relevance. The final dataset is in a relative large scale with roughly 603 observations for each of 28 variables.
Results
Table 1 provides the summary statistics of our database grouped by location, technology, fuel and ownership. The dataset accounts for a total of 24,000 MW of generating capacity, of which 23,227 MW (96.7%) is owned by non-utilities and 768 MW (3.3%) is owned by utilities.[7] Moreover, a total of 85 generating units are steam-turbine units and their capacity accounts for 92.6% of the total capacity. The number of units fueled by coal, oil and gas is 55, 18 and 36, respectively. The capacity of them is 64%, 16% and 20% for three fuel sources, respectively.
We conducted the MLM analysis using both STATA(8) and SAS (v8) statistical packages. The initial data cleaning and examination was done in SAS, which handled large datasets more efficiently. In addition, Proc Mixed with random statement was attempted to illuminate the MLM structures of the specified relationship (Equation (2)). However, due to limited proficiency with the software, we could only specify generators as the grouping variable (i.e., data within each generator is correlated assuming AR(1)) but not able to include higher level indicators (e.g., location, ownership, and technology). We resolved to include these indicators as level-2 covariates in the regression model and presented the SAS results in Table 2.
Table 3 summarized our modeling results using STATA. The fixed- and random-effects models give not compatible coefficient estimations for fixed effects. For instance, the coefficient of load estimated by fixed – and random-effects models is -.0000344 and .0001144, respectively. The inconsistency could be due to different modeling specification. In general, we pay with variance explained by the model (larger SE’s) to gain generalizability in random effects models.
The coefficient of ozone is the most policy relevant. Based on random-effects model, it suggests that in average, the implement of NOx results in a reduction of 0.105316 lb/mmBTU in hourly NOx emission rate. If we assume the heat rate is 10000 BUT/kWh, the estimated incurred cost is roughly 1.05316 $/MWh. For a coal unit with capacity equal to 1 MW and a capacity factor of 0.8, the annual cost becomes 7,381 $/year. Given its 95% CI, the estimated cost will be between 9461.1 and 5300 $/year.
Discussion
We presented an application of using MLM to quantify the policy impacts on energy industries. To solving models using STATA is subject to curse of dimension. As number of observation and variables increases, the integration procedure embedded in the package becomes extremely slow. For instance, it takes more than three hours to finish the run for our simplest dataset. However, our sample procedure should provide unbiased estimation of coefficients since the sample scheme is independent from other factors.
The comparison between fixed – and random-effects clearly shows that the estimations from fixed-effects models could sometimes be misleading. This can be attributed to the fact that fixed-effects model fail to take into accounts the underlying hierarchical nature of the problems. The results also show that the latent variables effectively explain a significant amount of variance, which is not explained in the fixed-effect model.
[pic]
Figure 1. Input NOx emission rate vs. output level (ST: steam and CC: combined cycle)
[pic]
Figure 2. Histograms of input NOx emission rate (ST: steam and CC: combined cycle)
[pic]
Figure 3. Distribution of Output (ST: steam and CC: combined cycle)
Table 1. Generator characteristics (n=134), EPA emission data of PJM Region, 2003
| |Total (n=134) |Study (n=88) |
|Characteristics | | |
| |Frequency |% |Frequency |% |
|Location | | | | |
|METED |6 |4.48 |5 |5.68 |
|PPL |14 |10.45 |10 |11.36 |
|BGE |12 |8.96 |11 |12.50 |
|PEPCO |15 |11.19 |8 |9.09 |
|JCPL |17 |12.69 |9 |10.23 |
|DPL |14 |10.45 |11 |12.50 |
|PECO |10 |7.46 |5 |5.68 |
|PSEG |29 |21.64 |18 |20.45 |
|AECO |10 |7.46 |6 |6.82 |
|PENELEC |7 |5.22 |5 |5.68 |
|Technology | | | | |
|ST |85 |63.43 |67 |76.14 |
|CC |24 |17.91 |21 |23.86 |
|CT |5 |3.73 | | |
|GT |20 |14.93 | | |
|Fuel type | | | | |
|Coal |55 |41.04 |44 |50.00 |
|Oil |25 |18.66 |13 |14.77 |
|Gas |54 |40.30 |31 |35.23 |
|Utility | | | | |
|Yes |123 |91.79 |81 |92.05 |
|No |11 |8.21 |7 |7.95 |
Table 2. SAS Proc Mixed output of the specified model (Equation (2)) with random intercept and assuming autoregression correlation between observations
Covariance Parameter Estimates
Cov Parm Subject Estimate
Variance gid 0.01339
AR(1) gid 0
Residual 0.01269
Solution for Fixed Effects
Standard
Effect LOCATION Estimate Error DF t Value Pr > |t|
Intercept 10.0557 3.0227 72 3.33 0.0014
price -0.00003 0.000037 17E3 -0.76 0.4483
load 0.000257 0.000012 17E3 21.06 |z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
load | .0001144 .0000819 1.40 0.162 -.000046 .0002749
capacity | .0001008 .0000992 1.02 0.309 -.0000936 .0002952
nyear | -.005696 .0017011 -3.35 0.001 -.00903 -.0023619
tech_st | .0269664 .0945173 0.29 0.775 -.1582842 .2122169
price | .0001833 .0003079 0.60 0.552 -.0004201 .0007868
oil | -.0128868 .0517603 -0.25 0.803 -.1143351 .0885614
gas | -.1504444 .0626391 -2.40 0.016 -.2732147 -.0276741
ozone | -.105316 .015147 -6.95 0.000 -.1350036 -.0756284
_cons | 11.53186 3.332747 3.46 0.001 4.999797 18.06392
------------------------------------------------------------------------------
Variance at level 1
------------------------------------------------------------------------------
.01434908 (.00127843)
Variances and covariances of random effects
------------------------------------------------------------------------------
***level 2 (gid)
var(1): .01376243 (.00307607)
***level 3 (location)
var(1): .00295807 (.00270522)
***level 4 (utility)
var(1): 1.011e-18 (5.230e-11)
------------------------------------------------------------------------------
References
[1] Samet, J., F. Dominici, et al. (2000). "Fine Particulate Air Pollution and Mortality in 20 U.S. Cities, 1987-1994," N Engl J Med 343(24): 1742-1749.
[2] Katsouyanni, K., G. Touloumi, C. Spix et al. (1997). "Short-term effect of ambient suphur dioxide and particulate matter on mortality in 12 European cities: results from times series data from the APHEA project: Air Pollution and Health: a European Approach." British Medical Journal, (314)1658-63.
[3] Thurston, G.D., K. Ito (1999). "Epidemiological studies of ozone exposure effect." In: Holgate ST., Samet JM, Koren HS, Maynard RL, eds. Air pollution and health. San Diego, Calif.: Academic Press, 177-85.
[4] Touloumi, G., K. Katsouyanni, D. Zmirou, et al. (1997). "Short-term effects of ambient oxidant exposure on mortality: a combined analysis within the APHEA Project," American Journal of Epidemiology, 146: 177-85.
[5] Bell, LM, McDermott, A, Zeger, SL, Samet, JM, and Dominci, F (2004). “Ozone and Short-term Mortality in 95 US Urban Communities, 1987-2000,” JAMA, 292(19) 2372-2378.
[6] Seinfeld, J.H. and S.N. Pandis (1997). "Atmospheric Chemistry and Physics." Willey Interscience.
[7] US Environmental Protection Agency (EPA),
[8] Ferrel A., Carter R., and Raufer R. (1999). “ The NOx budget: market-based control of tropospheric ozone in the northeastern United States,“ Resources and Energy Economics, (21) pp 103-124
[9] Stavins, RN (1998) "What Can We Learn from the Grand Policy Experiment? Positive and Normative Lessons from SO2 Allowance Trading." Journal of Economic Perspectives, Volume 12, Number 3, pages 69-88, Summer 1998.
[10] Joskow, P. and Rose, N (1990). "The Diffusion of New Technology: Evidence From the Electric Utility Industry", Rand Journal of Economics, 21:3, 1990, pp 354-373.
[11] Kira M. Fabrizio, Nancy L. Rose, Catherine Wolfram, "Does Competition Reduce Costs? Assessing the Impact of Regulatory Restructuring on U.S. Electric Generation Efficiency" (December 2004), available at
[12] James B. Bushnell and Catherine Wolfram, "Ownership Change, Incentives and Plant Efficiency: The Divestiture of U.S. Electric Generation Plants" (March 2005), available at
Appendix 1
The analysis of this project is also implemented in the STATA package. The full code includes the fixed- and random-effects on intercepts on generators, ownership and location is displayed as below.
*** 5/19/2005
clear
set memory 100m
capture log close
set more off
log using "C:\My Stuffs\My Research\course\2005 spring\bio140.656 mlm\project\code\model1.log", replace
use "C:\My Stuffs\My Research\course\2005 spring\bio140.656 mlm\project\data\all\long1400.dta", clear
*** fixed effect*
reg nox load utility nyear tech_st capacity price oil gas ozone loc1 loc2 loc3 loc4 loc5 loc6 loc7 loc8 loc9 loc10
*** random intercept model ***
gllamm nox load capacity nyear tech_st price oil gas ozone, fam(gaussian) link(id) i(gid location utility) adapt
-------------------------------------------------------------------------------------------------
Appendix 2
Excerpt of SAS program files and output of the multilevel analysis carried out in this project.
-----------------------------------------------------------------------------------------------
*** Reverse reference group ***;
DATA long;
SET mlm.long;
ozone1=1-ozone;
oil1=1-oil;
gas1=1-gas;
tech_st1=1-tech_st;
utility1=1-utility;
DROP ozone oil gas tech_st utility;
run;
PROC MIXED data=long METHOD=REML;
CLASS location gid;
MODEL nox = price load ozone1 capacity nyear oil1 gas1 tech_st1 utility1/solution;
random int / SUBJECT=gid TYPE=AR(1) g;
TITLE 'Model 1: Random intercept, AR(1), SUBJECT=gid';
run;
PROC MIXED data=long METHOD=REML;
CLASS location gid;
MODEL nox = price load ozone1 capacity nyear oil1 gas1 tech_st1 utility1 location/solution;
random int / SUBJECT=gid TYPE=AR(1) g;
TITLE 'Model 2: Random intercept, AR(1), SUBJECT=gid, location included';
run;
-------------------------------------------------------------------------------------------------
Model 1: Random intercept, AR(1), SUBJECT=gid
The Mixed Procedure
Model Information
Data Set WORK.LONG
Dependent Variable nox
Covariance Structure Autoregressive
Subject Effect gid
Estimation Method REML
Residual Variance Method Profile
Fixed Effects SE Method Model-Based
Degrees of Freedom Method Containment
Class Level Information
Class Levels Values
LOCATION 10 1 2 3 4 5 6 7 8 9 10
gid 88 1 2 3 4 5 6 7 8 9 10 11 12 13
14 15 16 17 18 19 21 22 23 24
31 32 33 34 35 36 44 45 53 54
55 56 57 58 59 60 61 62 64 65
66 67 68 69 70 71 73 74 75 78
79 80 81 82 83 84 85 88 89 96
97 98 100 103 104 106 109 111
112 113 114 115 116 118 119
120 125 127 128 129 130 131
132 133 134
Dimensions
Covariance Parameters 3
Columns in X 10
Columns in Z Per Subject 1
Subjects 88
Max Obs Per Subject 365
Observations Used 17308
Observations Not Used 14812
Total Observations 32120
Iteration History
Iteration Evaluations -2 Res Log Like Criterion
0 1 -11985.25046117
1 3 -25913.47094925 0.00346465
2 1 -25914.73778762 0.00017066
3 1 -25914.83206194 0.00000118
4 1 -25914.83279734 0.00000000
Convergence criteria met but final hessian is not positive
definite.
The Mixed Procedure
Estimated G Matrix
Row Effect gid Col1
1 Intercept 1 0.01672
Covariance Parameter Estimates
Cov Parm Subject Estimate
Variance gid 0.01672
AR(1) gid 0
Residual 0.01269
Fit Statistics
Res Log Likelihood 12957.4
Akaike's Information Criterion 12954.4
Schwarz's Bayesian Criterion 12950.7
-2 Res Log Likelihood -25914.8
Null Model Likelihood Ratio Test
DF Chi-Square Pr > ChiSq
2 13929.58 |t|
Intercept 8.6521 2.9789 81 2.90 0.0047
price -0.00003 0.000037 17E3 -0.77 0.4397
load 0.000253 0.000012 17E3 20.87 ................
................
In order to avoid copyright disputes, this page is only a partial summary.
To fulfill the demand for quickly locating and searching documents.
It is intelligent file search solution for home and business.
Related download
- lecture 1 random vectors and multivariate normal distribution
- computational modeling of residual stress formation during
- 9 0 generator exciter and voltage
- residual scatter plot generator goa
- unit 1 exploring and understanding data 25 days
- the practice of statistics
- johns hopkins bloomberg school of public health
- the practice of statistics ap statistics ap statistics
Related searches
- illinois department of public health worker registry
- list of public health issues
- illinois department of public health cna
- department of public health complaints
- list of public health laws
- examples of public health issues
- master of public health jobs
- department of public health vaccines
- illinois department of public health license
- department of public health texas
- illinois department of public health forms
- department of public health and environment