Wind Power Capacity Assessment

Transcription

1 Wind Power Capacity Assessment Mary Johannis, BPA, representing Northwest Resource Adequacy Forum Northwest Wind Integration Forum Technical Working Group October 29,2009

2 March 2007 NW Wind Integration Action Plan ACTION 1: By July 2007, the Northwest Resource Adequacy Forum (NWRA Forum) should reassess its 15 percent pilot sustained wind capacity value using currently available data on wind plant operation during periods of peak load. In 2008, the NWRA Forum should further refine the sustained peaking capacity value of wind power using the improved wind resource data set of Action 3 and other available data. October 29, 2009 NW Wind Integration Forum 2

3 Phase I: Reassess15% Wind Capacity Value July 25, 2007 Forum Technical Committee Meeting Wind Capacity Subgroup calculate wind capacity value based on contribution to meeting load during sustained peak period of cold snap/heat wave events Contract with BorisMetrics to translate wind speed data into simulated data, to perform quality control of existing wind generation data and to evaluate the wind capacity value January 17 & February 28, 2008 Tech Meetings BorisMetrics developed a 4 th degree polynomial constrained econometric model to backcast hourly project output as a dependent variable of Pendleton wind speeds (E. & W. Gorge areas) Concern that statistical attributes of backcast generation do not match actual wind generation attributes, i.e. many more instances of zero generation in actual records than in backcast simulation October 29, 2009 NW Wind Integration Forum 3

4 Phase I: Placeholder Wind Capacity Value of 5% Selected Probability of Occurance BPA Analysis: 1/ Median = 7.3 % 50% of observations are below this value 50% of observations are above this value Mean or Average = 17.2 % Distribution of Wind Fleet's Capacity Factors In BPA's Control Area Newer Technology Historical Data Averaged over 6 Peak Hours (n = 60) During Winter 3 Day Cold Spells Nov, Dec and Jan from Capacity Midpoint (%) 1/ 7/8/08 Forum Tech Committee Meeting Historical record is insufficient to calculate statistically significant wind capacity factor over 18 hour sustained peak period during cold snaps Median capacity factor over 6 peak hours during cold snaps is 7.3% Adverse wind capacity factor 5% October 29, 2009 NW Wind Integration Forum 4

5 PHASE II: Long-term Plan to develop Wind Capacity Value Need sufficient years of hourly wind generation by wind site for GENESYS to perform Monte Carlo picks Options: Backcast Wind Generation using historical Anemometer records Develop Temperature-Correlated Synthetic Wind Generation Records October 29, 2009 NW Wind Integration Forum 5

6 BorisMetrics Contract Identified Issues Can wind speed be used to backcast wind generation? Example: East Gorge Generation Dec 2006 Why is there so little generation when the wind is blowing? This example points outs problem with using off-site anemometer to backcast wind generation October 29, 2009 NW Wind Integration Forum 6

7 BorisMetrics Contract Identified Issues Can a unique function calculating generation based on wind speed be determined? October 29, 2009 NW Wind Integration Forum 7

8 BorisMetrics Contract Identified Issues Pendleton Anemometer Data not Clean Pendleton Ave Annual Anemometer Data Height, reading method, and meter location are noted ' manual T-107 bldg 53' man terminal bldg 20' manual terminal bldg 30' automatic terminal bldg Ave Annual Wind Speed (mph) ' automatic terminal bldg Year October 29, 2009 NW Wind Integration Forum 8

9 Vansycle Backcast Case Study Vansycle has an anemometer on-site ½ mile from the nearest generator 6 miles from the furthest generator Wind speed data is available in 10 minute intervals for period Scada data is available in 5 minute intervals for period Vansycle Backcast should be doable Relatively long-term Generation Record Relatively clean Anemometer Record Wind Turbine Power Characteristics: Cut-in wind speed 4 m/s (8.9 mph) Nominal wind speed 15 m/s (33.6 mph) Stop wind speed 25 m/s (55.9 mph) October 29, 2009 NW Wind Integration Forum 9

10 Vansycle Study Backcasting not feasible Lessons learned: High R 2 of multivariate regression (without zeros) and residual analysis indicates that Persistence is an important feature in regression Other regressions have artificially high R 2 by including zeros Prediction interval of.3 is not sufficiently tight to backcast Backcasting Wind Generation for NW is NOT feasible Even on-site wind anemometers can be miles from some wind turbines resulting in the LACK of a unique correlation Due to the persistence feature of the regression cannot use other means to reflect randomness in the correlation Insufficient on-site anemometer data to backcast the entire NW wind generation fleet Conclusion: Develop Temperature-Correlated Synthetic Wind Generation Records October 29, 2009 NW Wind Integration Forum 10

11 Synthetic Wind Generation using K th Nearest Neighbor Method What is Kth Nearest Neighbor Method? For a time series of size N, randomly select a single or two consecutive of the N observations then select the third based on how close the lag(s) for the selected observation are to the randomly selected observation(s) For example, select two hours where the capacity is 0.3 & 0.4, respectively, then pull from observations that have capacities that are close to.3 for the observation 2 hours prior and.4 for the hour prior. Creating a subset of the K closest observations to draw from maintains the structure that is expected in the time series. Methodology is undergoing peer review A cross-correlated time series synthetic study presented to joint conference of Western North American Region of the Biometric Society and the Institute of Mathematical Statistics Paper using method for wind fleet capacity factor data submitted to IEEE Kth Nearest Neighbor Method presented to NERC RIS-IVGTF team October 29, 2009 NW Wind Integration Forum 11

12 Goal of Historic Temperature Correlated Synthetic Wind Gen It has been well established that temperatures affect load where extreme high or low temperatures translate into high loads. The synthetic wind power generation data recreates certain statistical characteristics of the original or observed wind power generation data set. The characteristics to focus on are: Distribution/Density Lag Structure or Persistence Cross-Correlation The long-term temperature-correlated wind generation records will be incorporated into the existing resource adequacy studies using the GENESYS model, which will perform Monte Carlo picks on temperature-years, thus pointing to synthetic wind generation and loads October 29, 2009 NW Wind Integration Forum 12

13 19.54% Average Observed Capacity Factor 33.1% 20.7% Based on 2006 through 2008 observationsof the BPA integrated Wind Fleet during heavy load hours as def i nedby NERCand maximum regional average temperature from Seatac, Portland, and Spokane airports. October 29, 2009 NW Wind Integration Forum 13

14 Wind Generation vs. Temperature October 29, 2009 NW Wind Integration Forum 14

15 Forum Wind Methods consistent with other Wind Forums NERC: Joint Integration of Variable Generation Task Force (IVGTF) Resources Issues Subcommittee (RIS) Task 1.2 (Capacity Value) and Task 1.4 (Flexible Resources to integrate Variable Generation) Teams IVGTF Report: WECC: Variable Generation Subcommittee (VGS) Planning Work Group Northwest: PNW Resource Adequacy Forum/NW Wind Integration Forum October 29, 2009 NW Wind Integration Forum 15

16 Wind Capacity Value Methods Effective Load Carrying Capability (ELCC) approach Evaluate effective wind capacity contribution based on LOLP studies with and without wind generation & same target (GENESYS Approach) Need sufficient wind generation data to simulate full range of generation under various conditions, especially if wind and loads correlated at times Need realistic depiction of combined uncertainties Contribution of variable generation to system capacity during high-risk hours using historical data Investigate contribution of wind capacity during heat wave and cold snap events in PNW because of evidence of statistical relationship between lack of wind generation when it gets very hot or very cold (Forum Wind Capacity Subgroup Approach) Correlation between resource contribution and the resource mix by system (e.g. what is appropriate for a hydro based system) Wind may contribute more in energy-limited system if certain amount of wind generation can be counted upon during drought October 29, 2009 NW Wind Integration Forum 16

17 Counting Wind toward Capacity Adequacy in the NW January 2009 Cold Snap BPA Balancing Authority Area Load & Total Wind Generation Jan. 5-25, MW BPA TOTAL WIND GENERATION 2000 BPA BALANCING AUTHORITY AREA LOAD /5/09 1/6/09 1/7/09 1/8/09 1/9/09 1/10/09 1/11/09 1/12/09 1/13/09 1/14/09 1/15/09 1/16/09 1/17/09 1/18/09 1/19/09 1/20/09 1/21/09 1/22/09 1/23/09 1/24/09 1/25/09 Date/Time (5-min increments) October 29, 2009 NW Wind Integration Forum 17

18 Synthetic Wind Generation: Historic Cold Snaps October 29, 2009 NW Wind Integration Forum 18

19 Synthetic Wind Generation: Historic Cold Snaps October 29, 2009 NW Wind Integration Forum 19

20 Simulated Wind Generation: Historic Heat Waves October 29, 2009 NW Wind Integration Forum 20

21 Simulated Wind Generation: Historic Heat Waves October 29, 2009 NW Wind Integration Forum 21

22 Observed Wind over all hours: Regional Load Duration Treating the wind as negative load changes the duration curve. Minimum distance between the two curves is about 1.6% of the nameplate. 99.5% of the hours have a 6.8% of the nameplate or greater contribution of wind toward reducing the load durations. Regional Load (MW) Duration Curves with and without Wind Adjustment 15,000 20,000 25,000 30,000 Minimum Difference = 35.7 MW at Observed Load of 31,962.4 MW Load Less Wind Duration Observed Load Duration Duration curve based on all observed hours from 2006 through 2008, regional load provided by Council, wind fleet adjustment based on observed capacity factors multiplied by current wind fleet nameplate. 0% 20% 40% 60% 80% 100% Percent of Hours Greater Than Load October 29, 2009 NW Wind Integration Forum 22

23 NW Wind Capacity Value using quasi-elcc Approach Difference between the percentiles of the load durations show us: Between the 10 th and 90 th percentiles the contribution of the wind fleet was fairly flat with a slight trend of more energy during the lower loads. During the highest observed loads the difference is minimal. Difference Between Duration Curves (MW) Difference of Regional Load Duration Compared to Load Less Wind Duration Maximum Difference at the Lowest Observed Loads ELCC Focus Minumum Difference at Highest Observed Loads 0% 20% 40% 60% 80% 100% Duration Percentile Corresponding to Difference October 29, 2009 NW Wind Integration Forum 23

24 NW Wind Capacity Value over 6 peak hours (high risk hrs) Alternately looking at the differences between the six peak hours with and without wind yields: A minimum difference of zero. 97% of the time, contribution is only.03% of nameplate in amw. 91% of the time, contribution of 1% of the nameplate in amw. Difference (amw) Difference Between Six Peak Load Hour ObservedLoadandObservedLoadLessWind 0% 20% 40% 60% 80% 100% Percentage of Observations Less Than Difference October 29, 2009 NW Wind Integration Forum 24

25 18-hour wind capacity 2006 First Day 18hr Mean Load (amw) 18hr Mean Load Less Wind (amw) Difference (amw) % of Integrated Fleet Nameplate Jan /16/ % Feb /16/ % Mar /8/ % Apr /17/ % May /16/ % Jun /26/ % Jul /23/ % Aug /7/ % Sep /5/ % Oct /30/ % Nov /27/ % Dec /18/ % October 29, 2009 NW Wind Integration Forum 25

26 18-hour wind capacity 2007 First Day 18hr Mean Load (amw) 18hr Mean Load Less Wind (amw) Difference (amw) % of Integrated Fleet Nameplate Jan /15/ % Feb /31/ % Mar /27/ % Apr /2/ % May /30/ % Jun /19/ % Jul /10/ % Aug /13/ % Sep /10/ % Oct /31/ % Nov /26/ % Dec /10/ % October 29, 2009 NW Wind Integration Forum 26

27 18-hour wind capacity 2008 First Day 18hr Mean Load (amw) 18hr Mean Load Less Wind (amw) Difference (amw) % of Integrated Fleet Nameplate Jan /22/ % Feb /4/ % Mar /26/ % Apr /31/ % May /17/ % Jun /30/ % Jul /7/ % Aug /13/ % Sep /15/ % Oct /22/ % Nov /24/ % Dec /15/ % October 29, 2009 NW Wind Integration Forum 27

28 Status of Wind Discussions 10/16/09 Resource Adequacy Forum Technical Committee Meeting Evidence suggests that 5% Placeholder Value for Wind Capacity is too high Use different WINTER & SUMMER values for Wind Capacity in regional Planning Reserve Margin (PRM) Calculation: PRM = ( 1 18 hr regional resources regional 1 in 2 load)/ 1 18 hr regional 1 in 2 load Current physical resource adequacy thresholds are: PRM winter 23% PRM summer 24% October 29, 2009 NW Wind Integration Forum 28

29 Next Steps In the Short-term, refine Wind Capacity Value in PRM Equation SOME OPTIONS: Select 95 th percentile wind from summer & winter all hour wind vs. load duration curves ~ Slide 23 Select 95 th percentile wind from summer and winter 6 peak or 18 peak load hour duration curves ~ Slide 24 Using actual (and possibly synthetic) wind generation data over historical heat wave and cold snap events, calculate average wind capacity contribution over 18 hour sustained peak or 6 peak hour period In the Long-term, perform true ELCC evaluation using Monte Carlo picks of temperature-years to point to loads and wind generation Create additional long-term temperature-correlated synthetic wind generation records for use in GENSYS October 29, 2009 NW Wind Integration Forum 29