The energy generation in Brazil is dominated by hydropower plants with large reservoirs which are arranged in cascades. It is forecasted to represent 11.6% of the Brazilian installed capacity by 2024  . According to  , the increase of wind power capacity worldwide in 2017 was 535 GW, reaching a total of 539 GW. China alone installed 19.7 GW in 2017, reaching 188 GW in total. Brazil was sixth in percentage increase in 2017 and is in the top ten countries in terms of total installed capacity.
In September 2016, the capacity factor registered by the state of Rio Grande do Norte was almost 57% and all states of the Northeast had capacity factor higher than 50%  . Wind power in the Northeast was able to meet 39% of the demand in this same month, when the reservoir levels were experiencing a shortage (Figure 1).
The Brazilian System Operator (ONS) dispatches the power plants to meet the demand. Cepel, the Brazilian Electric Energy Research Center, implemented the computational program Newave with the stochastic dual dynamic programming (SDDP) technique  to solve the large-scale long-term hydrothermal problems for the Brazilian system. The model aims to determine the optimal amount of hydro and thermal resources for operation planning and minimize the expected value of the operation cost  . In this study, we use the Newave model for the simulations.
The inflows are represented stochastically in the models as the future is unknown. Uncertainty is represented by a tree of scenarios where each branch of the tree is a hydrologic scenario. The Newave considers the variables’ uncertainty through synthetic inflow scenarios to all subsystems  . GEVAZP, the synthetic inflow scenarios generation model, is connected to the Newave. It selects a stochastic time series PAR (p) algorithm to guarantee similarity between historical and synthetic series  . The PAR (p) is an auto regressive periodic function where p can vary from 1 - 6 (months), so each stochastic inflow can be dependent on the inflow that occurred in the same places up to 6 months before  .
 shows wind and hydrologic generation scenarios with the same stochastic model PAR (p) for the operation planning and recommends the use of the wind series for the Brazilian planning hydrothermal operation.  developed a new methodology to build the scenario trees, removing non-linearity from the equations, which could make the future cost calculation in Newave an unfeasible task. Even though some studies recommend the stochastic model PAR (p), there are some limitations and questions remain whether this model is ideal for wind power series generation.
In  wind power generation is represented as a small hydro power generation and shows that the complementarity between wind power and hydropower results. Even though this study shows the same focus of our work, it does not detail the methodology for representing the stochastic series, how to consider the wind power plants as small hydro or even use the official model.
Figure 1. Reservoirs and wind power generation (Source: ONS, 2016-Juliana Mummey et al.).
 analyses the impact of considering the uncertainty of wind power generation scenarios using the model SDDP. The results show an increase of the costs and point out that by changing the future cost methodology the impact could be reduced.
Given the wind power generation variability and considering the increase in the share of wind power in the Brazilian electricity matrix, this paper aims to present a study about the stochastic representation of wind power generation in the Brazilian official models. The impacts on the other source’s generation and when wind power plants are represented as run-of-river power plants, will be discussed ahead.
The current study differs from the existing literature because it uses an optimization simulation model (Newave) to represent the wind and hydro generation stochasticity, while previous work doesn’t represent explicitly the wind and hydro production variability. In fact, the correlation between the power generation of these two sources is obtained from the historical registered behavior of the water flows at the main Brazilian river basins and series of wind speed, which were reconstructed based on data obtained from meso-scale meteorological models  . In addition to that, this study shows the production complementarity between wind and hydro generation as a result of the simulation model when using the stochastic formulation for both kind of sources. It is worthwhile to point out that currently the Brazilian System Operator model considers the wind power deterministically and therefore it does not represent the variability of the wind power source.
This introduction describes the problem focused in this work, a brief survey of the state of the art of wind power in Brazil, expansion perspectives and studies about wind power series. Section 2 shows the methodology involved such as: 1) wind speed database and wind power generation, and 2) wind power representation in the Brazilian system. Section 3 shows the results of the simulations: 1) the rebuilt historic of wind power generation (considering the occurrence of hydraulic and wind power production varying simultaneously) and comparison to the official historical simulation (average monthly historical values, each month is the same every year); 2) based on “1” and a hypothetic scenario of an increase of wind power capacity and an increase of the Northeast and South demands; and 3) same as “1” but considering synthetic series. Section 4 concludes this article.
This work uses two sources of wind speed data representing 16 sites in the Northeast and the South of Brazil. We plotted the Weibull distribution to fit the wind speed distribution. Then, a power curve of each site is used to calculate the wind power generation. All this data represents an input to the optimization model Newave in order to run the program. Figure 2 represents the schematic diagram of the methodology.
2.1. Wind Speed and Wind Power Generation
The methodology of wind speed historical data reconstruction and its power generation was based on  . Two databases were considered: Vortex (Mesoscale atmospheric model on-line that estimates wind speed for places without measurements) and NOAA (numerical model from National Oceanic and Atmospheric Administration).
Figure 2. Schematic diagram of the methodology.
The Vortex database provides the hourly wind speed from 1982 to 2014 for every 10 meters of height, starting at 50 meters up to 150 meters. The NOAA database presents daily values of wind speed at 42 meters since 1947. There were 16 sites considered and Figure 3 shows the seasonality of hourly wind speed data of the site at Paracuru, as an example. It can be observed that there is a trend of higher wind speeds and shorter variability in the second six months of the year.
We transformed both series (NOAA and Vortex) into the same period to be comparable, using the same years from 1982 to 2010 and daily wind speeds. We also considered all series with cross-correlation greater than 0.8 and then a vertical extrapolation to transform the NOAA data from 10 meters to 120 meters high with the logarithmic curve with the Vortex data available. We calculated monthly standard deviation with the Vortex data and based on the daily average wind speed we plotted the Weibull distribution (the Weibull probability density function was used to fit the wind speed distribution) for each month. Then we calculated the energy generated with the power curve of each site and the Weibull distribution fit to the historical data. Based on the daily wind power generation, we added all days of the month to transform into monthly wind power generation from 1948 to 2014.
2.2. Wind Power Representation
The model Newave has four main subsystems and nine regions: Paraná, Itaipu, Madeira, Teles Pires and Southeast (five regions in the Southeast subsystem); North and Belo Monte (two regions in the North subsystem); Northeast subsystem (one region) and South subsystem (one region). The maximum capacity of regions for the model is 15, so we included four regions in the Northeast: CE (Ceará), BA (Bahia), PI + PE (Piauí + Pernambuco) and RN (Rio Grande do Norte), and two regions in the South region: SC (Santa Catarina) and RS (Rio Grande do Sul), as Figure 4. These new regions are wind basins and represent areas of the system with high concentration of sites and a correlated wind regime.
Each region has the wind power installed capacity of the state and the respective expansion according to the Monitoring Department of Electric System (DMSE). Each wind power site was included in the program as if it was a run-of-river power plant. As there are many sites, with different seasonalities inside the same state, we put the values of the closest sites and their expansion according to where they are in the Abeeólica (Wind Power Association) data  .
There are three simulations with the Newave model for this study: 1) considering the re-constructed wind power historical generation and a comparison with the official historical case of August 2016; 2) based on “1” but considering an increase of the wind power capacity and also an increase of the Northeast and South demands; and 3) based on “1” but considering synthetic series. With these
Figure 3. Wind speed seasonality.
Figure 4. Representation of regions.
cases, we analyze the complementarity of wind power generation with hydropower, verifying in dry seasons a) if the wind power can generate more; b) if with more wind power capacity, what is the role of wind power generation to meet the demand; and c) if the behavior with synthetic series is compatible with the observed one when using historical series.
We chose two historical critical (dry years) periods to analyze the data. The first one refers to the years 1951 to 1955 (the horizon of Newave is 5 years). These years presented low inflow in the Northeast and re-constructed wind power generation with some months below the average. The second period refers to 2010 to 2014. These years presented low inflow in the Northeast and re-constructed wind power generation above average. Considering the synthetic series, the data are analyzed through the average values.
3.1. Simulation Considering the Wind Power in Regions
The marginal costs of the simulation, when representing variability with wind power in regions and with scenarios coming from the historical series, tend to be lower than the official case, see Table 1. Considering the period of 1951 to 1955, as it is a dry period, there is a need to dispatch more thermal power plants to save more water in the reservoirs for the future, which is why there is an increase of the costs in both the official case and in the simulation. However, the costs of
Table 1. Marginal costs and generation.
the simulation are lower than in the official case. This happens because, as in the official case, the wind power generation is deterministic and, hence, there is no change of this generation if it is a dry season or not. Considering the simulation, in February of 1955, for example, even though there is a low hydraulic inflow in the Northeast and the re-constructed wind power is also low, there is a higher wind power generation, with less need of thermal power generation, reducing the marginal costs of operation, see Figure 5. It shows the complementarity of the hydro with the wind power in the Northeast in the simulation. Considering the period of 2010 to 2014, as the re-constructed wind power was higher than in 1951 to 1955, wind power generation could complement even more with the hydropower generation, see Table 1.
As an example, the last February of the period (Table 2), we can see that despite the wind power generation in the average simulation being lower than the official deterministic case as February is a month with low wind power generation, the marginal costs of the simulation are still lower than the official case. This happens because there was more hydropower generation with the increase of wind power generation in previous months, allowing more water in the reservoirs. Considering the hydrologic scenarios corresponding to the historical years of 1951-1955 and 2010-2014, which represent the dry seasons, there is an increase in the wind power generation in the simulation, complementing the hydropower generation and reducing the marginal costs, even in a month with low wind power generation.
Figure 5. Marginal costs of the simulation vs. official case.
Table 2. Marginal costs and generation.
3.2. Simulation Considering Wind Power in Regions and an Increase of the Wind Power Capacity
This hypothetical case adds 10 times the wind power capacity to the system proportionally to both the Northeast and South subsystems. We used this increase as the wind power capacity today represents a tenth of the hydroelectric plants in Brazil and we wanted to understand the impact that the same proportion of both sources would represent. It also increases the demand of both subsystems, otherwise the model would be optimistic. Considering the second six months of the years, the wind power generation can be the main source to attend the demand because in this part of the year there are higher wind speeds and greater wind power generation. However, considering the first six months, wind power generation could not be enough to meet the demand and there could be a need to increase hydro and thermal power generation, increasing the marginal cost of the Northeast, once it becomes the region with the most increase of the demand. In the example of February, in Table 3, even though it is February and the wind speeds are not high, considering the dry season of 1951-1955, the wind power generation was greater than the average. Considering the years of 2010-2014, the wind power generation is much greater, showing the complementarity of the case.
Table 3. Marginal costs and generation.
3.3. Simulation Considering Wind Power in Regions and Synthetic Series
In the case of synthetic series, the marginal costs of the simulation also present reductions in almost all months for all subsystems compared to the official case, see Figure 6. However, considering February of the last year, which in this case is February of 2020, wind power generation was lower than in the official case with synthetic series in the Northeast and higher in the South, see Table 4.
Even with the lower wind power generation in this specific month, the marginal cost was lower than the official case with synthetic series. This is because as there is more wind power generation in the previous months, the reservoir storage levels are higher and consequently there is more availability of hydropower to generate energy and reduce the costs.
Figure 7 shows that there is more wind power generation in almost all months in the simulation compared to the official case in the Northeast and in the South. However, in the wet season (period between December and April) of the years, wind power in the Northeast showed lower values, emphasizing the seasonality. There is a decrease of marginal costs in the simulation because of the higher hydropower generation. This happens because the storage level of the reservoirs is higher than the official case, with the higher wind power in previous months.
Figure 8 shows the complementarity of the wind power generation with the hydropower generation in the Northeast comparing the official case and the simulations and considering synthetic series. As the official case is deterministic, wind power generation considers the same seasonality during those years.
Table 4. Marginal costs and generation.
Figure 6. Marginal costs of the simulation vs. Official case with synthetic series.
Figure 7. Wind power generation vs. official case with synthetic series.
Figure 8. Complementarity of hydropower and wind power generation.
4. Final Remarks and Conclusions
The aim of this work was to present a study to stochastically represent wind power generation in the optimization model (Newave) of the Brazilian electricity system, also, to evaluate how the complementarity of wind power with hydropower can influence the hydro and thermal power dispatches considering variability of the wind power source, as well as the consequences for the short-term marginal prices and spot market prices.
There were three simulations comparing the official case with the wind power represented in regions with historical and synthetic series. The results presented were similar from the conceptual and qualitative point of view. In most of these cases, there is a reduction in the marginal costs of operation with the simulation and less variability in the costs. It happened because there is more complementarity in these cases, reducing the need of thermal dispatch to meet the demand. In a deterministic simulation, as nowadays state of the art, the huge variability of the wind power production is not considered and, as a consequence, the complementarity effect between hydro and wind generation is not taken into account as an important benefit to the System.
When wind power capacity was increased, it was used as the main source to meet the demand. In the second six months, there were no problems for supplying the demand, because the wind power generation is high in these periods. However, considering the first six months, the simulation showed that the wind power was not enough to attend the demand, lessening the system reliability. These scenarios of different seasons are not represented in the deterministic models, which show the relevance of considering such scenarios.
Considering synthetic series there is also an increase of wind power generation in most of the scenarios, reducing the marginal costs of operation. However, we observe that in some wet periods of the years, wind power generation was lower than the official case, showing the complementarity and seasonality of the source, highlighted even more with the stochastic representation of wind power.
This is a case study specific for the Brazilian system but presents methodology that can be replicated to problems of the same nature in any part of the world. To conduct a more detailed analysis we recommend 1) the use of a database with wind speed measurements with anemometers together with a re-analysis database; and 2) the use of a higher number of sites, covering greater parts of the regions.
Finally, it is important to emphasize that this work does not intend to forecast wind power generation, but rather to show that the deterministic methodology considered in the optimization models of the Brazilian system needs to be reevaluated.
 Maceira, M.E.P., Duarte, V.S., Penna, D.D.J., Moraes, L.A.M. and Melo, A.C.G. (2008) Ten Years of Application of Stochastic Dual Dynamic Programming in Official and Agent Studies in Brazil—Description of the NEWAVE Program. 16th Power Systems Computation Conference, Glasgow, July 2008.
 Maceira, M.E.P. and Mercio, C.M.V.B. (1997) Stochastic Streamflow Model for Hydroelectric Systems. 5th International Conference PMAPS—Probabilistic Methods Applied to Power Systems, Vancouver, September 1997.
 Castro, C.M.B., Marcato, A.L.M., Souza, R.C., Silva Junior, I.C., Oliveira, F.L.C. and Pulinho, T. (2015) The Generation of Synthetic Inflows via Bootstrap to Increase the Energy Efficiency of Long-Term Hydrothermal Dispatches. Electric Power Systems Research, 124, 33-46.
 Cesar, T.C., David, P.A.M.-S., Pereira, A.O., Souza, R.A. and Carvalho, R.N.F. (2011) Regularization of Energy Supply—The Role of Complementarity. Electrical Systems Planning Study Group XXI SNPTEE.
 El-Heri, Y.S., Borba, B.S., Bezerra, B., Carvalho, M.R.M. and Dall’orto, C.E.R.C. (2016) Analysis of the Energy Impact of the Variability of Wind Energy Production in the Brazilian Electric System. Brazil Wind Power.
 Witzler, L.T., Ramos, D.S., Camargo, L.A.S. and Guarnier, E. (2016) Reconstruction of Wind Generation Historical Series Aiming at the Analysis of Energy Complementarity: Methodology and Applications. 13th International Conference on the European Energy Market, Porto, 6-9 June 2016.