1,253
Views
1
CrossRef citations to date
0
Altmetric
Article

Role of Variable Renewable Energy Penetration on Electricity Price and its Volatility across Independent System Operators in the United States

, , , , , , ORCID Icon, & show all
Article: 2158145 | Received 22 Dec 2021, Accepted 29 Nov 2022, Published online: 22 Feb 2023

Figures & data

Table 1. Temporal range, resolution, and VRE composition for ISO data.

Table 2. Literature summary for US studies (domestic) on VRE impact on price. The specified VRE, solar or wind, is given in the rightmost column.

Figure 1. Map of independent system operators (ISOs) in the United States. ISOs are responsible for coordinating, controlling, and monitoring the electricity grid within their region of control. There are seven ISOs in the United States, and the regions covered by these ISOs account for 67% of the total electricity demand in the country. Adapted from Climate-XChange (Citation2021).

Figure 1. Map of independent system operators (ISOs) in the United States. ISOs are responsible for coordinating, controlling, and monitoring the electricity grid within their region of control. There are seven ISOs in the United States, and the regions covered by these ISOs account for 67% of the total electricity demand in the country. Adapted from Climate-XChange (Citation2021).

Figure 2. The energy price data for each ISO is detrended by the average price for hour, season, and weekend trends. The observed price data, shown with black dots, is plotted by day of the year. The red line indicates the estimated trend. The resulting detrended price is the observed value minus the average. Note there are repeated observations for multiple years and the y-axis is clipped to the range $10–150.

Figure 2. The energy price data for each ISO is detrended by the average price for hour, season, and weekend trends. The observed price data, shown with black dots, is plotted by day of the year. The red line indicates the estimated trend. The resulting detrended price is the observed value minus the average. Note there are repeated observations for multiple years and the y-axis is clipped to the range $10–150.

Figure 3. CAISO raw system prices by seasons and hour of day. The plot shows the estimated diurnal pattern (red line) by season. For plotting purposes, prices greater than 150 are excluded. Orange dots depict negative prices and teal dots depict positive prices (all in USD).

Figure 3. CAISO raw system prices by seasons and hour of day. The plot shows the estimated diurnal pattern (red line) by season. For plotting purposes, prices greater than 150 are excluded. Orange dots depict negative prices and teal dots depict positive prices (all in USD).

Figure 4. Quantiles of skew t-distribution based on the estimated pointwise means of the parameters listed in . Points correspond to seasonally detrended system price and the y-axis is trimmed to the inner 99% of detrended price; red points correspond to a negative price before seasonal detrending.

Figure 4. Quantiles of skew t-distribution based on the estimated pointwise means of the parameters listed in Table A1. Points correspond to seasonally detrended system price and the y-axis is trimmed to the inner 99% of detrended price; red points correspond to a negative price before seasonal detrending.

Figure 5. VRE (solar and wind) penetration. Shown here is the percentage of VRE generation to total generation in each of the six ISOs considered. Southwest Power Pool (SPP) achieved the highest VRE penetration based on the data used and for the study period considered, followed by CAISO, MISO, ISONE, NYISO, and PJM.

Figure 5. VRE (solar and wind) penetration. Shown here is the percentage of VRE generation to total generation in each of the six ISOs considered. Southwest Power Pool (SPP) achieved the highest VRE penetration based on the data used and for the study period considered, followed by CAISO, MISO, ISONE, NYISO, and PJM.

Table 3. Descriptive Statistics Summary. This shows the minimum, median, maximum, and mean statistics for the detrended system price (in USD) and temporal price volatility (in USD).

Figure 6. VRE penetration (in %) and detrended electricity price (in $). The six panels show the relationship between VRE penetration and detrended electricity price for NYISO, MISO, PJM, CAISO, ISONE, and SPP, respectively. The continuous lines show the quantiles with the colors yellow, green, blue, and purple, representing the 0.25th, 0.50th (median), 0.75th, and the 0.90th quantiles, respectively. Each panel’s axis limits are adjusted separately in order to clearly show the quantiles’ trend direction.

Figure 6. VRE penetration (in %) and detrended electricity price (in $). The six panels show the relationship between VRE penetration and detrended electricity price for NYISO, MISO, PJM, CAISO, ISONE, and SPP, respectively. The continuous lines show the quantiles with the colors yellow, green, blue, and purple, representing the 0.25th, 0.50th (median), 0.75th, and the 0.90th quantiles, respectively. Each panel’s axis limits are adjusted separately in order to clearly show the quantiles’ trend direction.

Figure 7. The derivative of the detrended system price (in USD) and detrended system price volatility (temporal volatility in USD) with respect to the percentage of VRE for the PJM independent system operator.

Figure 7. The derivative of the detrended system price (in USD) and detrended system price volatility (temporal volatility in USD) with respect to the percentage of VRE for the PJM independent system operator.

Figure 8. Estimated effects of VRE percentage on the skewness (ν) parameter for a skew t-distribution for detrended price using a log link. Gray bands indicate the 95% confidence intervals. The red reference line at 1 denotes symmetry.

Figure 8. Estimated effects of VRE percentage on the skewness (ν) parameter for a skew t-distribution for detrended price using a log link. Gray bands indicate the 95% confidence intervals. The red reference line at 1 denotes symmetry.

Figure 9. VRE penetration (in %) and temporal price volatility (in USD). The panels show the relationship between VRE penetration and temporal detrended system price volatility for NYISO, MISO, PJM, CAISO, ISONE, and SPP, respectively. The continuous lines show the quantiles with the colors yellow, green, and purple, representing the 0.50th (median), 0.75th and the 0.90th quantiles.

Figure 9. VRE penetration (in %) and temporal price volatility (in USD). The panels show the relationship between VRE penetration and temporal detrended system price volatility for NYISO, MISO, PJM, CAISO, ISONE, and SPP, respectively. The continuous lines show the quantiles with the colors yellow, green, and purple, representing the 0.50th (median), 0.75th and the 0.90th quantiles.

Figure 10. Difference between 0.75th and 0.25th quantile (yellow), 0.75th and 0.50th (teal), and 0.90th and 0.50th (purple) as plotted with approximate 95% credible intervals. It is important to note that for several of the ISOs, a decrease in the spread of the distribution (volatility with respect to VRE) is observed in the quantiles as measured by all three differences, but SPP and CAISO have curvilinear trends in the difference between 0.75th and 0.25th, which is caused by the skewness changing from positive to negative with increasing VRE penetration.

Figure 10. Difference between 0.75th and 0.25th quantile (yellow), 0.75th and 0.50th (teal), and 0.90th and 0.50th (purple) as plotted with approximate 95% credible intervals. It is important to note that for several of the ISOs, a decrease in the spread of the distribution (volatility with respect to VRE) is observed in the quantiles as measured by all three differences, but SPP and CAISO have curvilinear trends in the difference between 0.75th and 0.25th, which is caused by the skewness changing from positive to negative with increasing VRE penetration.

Figure A1. The three panels show an illustration of the detrending procedure for a snapshot of data from CAISO in the second half of July 2019. In (A) the raw system price is shown (dots) along with the estimated seasonal trend (red line). We detrend the raw system price by subtracting the trend to obtain the detrended system price (B). The detrended system price represents how much the system price differs from the seasonal average; note that most of the time it is near 0. From (B), we calculated the exponentially weighted moving average standard deviation (C) that is a temporal measure of the variability of the detrended system price- in other words, How much does the detrended price vary from the detrended prices in the recent past?

Figure A1. The three panels show an illustration of the detrending procedure for a snapshot of data from CAISO in the second half of July 2019. In (A) the raw system price is shown (dots) along with the estimated seasonal trend (red line). We detrend the raw system price by subtracting the trend to obtain the detrended system price (B). The detrended system price represents how much the system price differs from the seasonal average; note that most of the time it is near 0. From (B), we calculated the exponentially weighted moving average standard deviation (C) that is a temporal measure of the variability of the detrended system price- in other words, How much does the detrended price vary from the detrended prices in the recent past?

Table A1. Parameter values for the generalized skew t-distribution hourly detrended price. An example of interpretation of the effect of VRE on the centrality and skewness of the distribution of hourly price is an increase in VRE decreases both centrality and skewness of hourly price in CAISO (β1μ<0 and β1ν<0). The interpretation also holds for ISONE, NYISO and SPP.

Data availability statement

Raw data and code can be found here: https://github.com/Qunlexie/VRE-Impact-on-Price-and-Volatility.