Određivanje Vetroenergetskog Potencijala Vojvodine Linearnom Interpolacijom Podataka Sa Meteoroloških Stanica

ODREĐIVANJE VETROENERGETSKOG POTENCIJALA VOJVODINE LINEARNOM INTERPOLACIJOM PODATAKA SA METEOROLOŠKIH STANICA ASSESSMENT OF WIND POTENTIAL IN VOJVODINA BY LINEAR INTERPOLATION OF METSTATIONS DATA N. Vasiljević* and M. Zlatanović ** PDRB Kolubara D.O.O. Lazarevac* Elektrotehnički fakultet, 11120 Beograd, Bulevar Kralja Aleksandra 73** Abstract: By linear interpolation of data from meteorological stations the wind potential map of Vojvodina, north region of Serbia, was generated. As a reference source of long term wind data eight met stations were selected. From this data and local terrain roughness, orography and obstacles data in surrounding area, the regional wind climate tables were calculated using WAsP methodology and licensed software. Digital topography maps were used to generate local wind climate atlases. Validity of each local wind atlas in 20km*20km area around the met station was assumed based on mostly plain and homogeneous orography of Vojvodina. For the regions outside 20km*20km areas covered by met station data, linear interpolation of wind data was used between minimum two, and maximum three nearest wind atlases to complete Vojvodina wind atlas. As the result, 66 single wind atlases were derived, 58 of them by linear interpolation method and 8 from original met station data. By assembling regional wind atlases the wind atlas of Vojvodina was generated. Key words: Observed Wind Climate, Renewable Energy, WAsP, Wind Atlases, Wind Energy 1. INTRODUCTION In April 2008 the total world installed capacity of wind turbines reached the figure of 100GW. Today wind energy supplies 5% of the European Union electricity while 12% is planned for 2020, which assumes installation of 180GW wind turbines power. The benefits of using wind energy include economy, employment, security of supply, electricity market and environment benefits. As the world leader in wind energy sector EU has today 150 000 jobs and 380 000 are expected in 2030. General opinion based on existing data is that all European countries have good wind energy potential. European wind atlas prepared by Risø Institute in 1989 was completed by atlases of new EU and candidate countries [1]. The wind atlas of Serbia prepared by the same methodology still is not published. Some reports contain preliminary data based on meteorological stations observation with indicated areas of good wind potential [2,3]. Two projects based on specialized wind measurement stations data financed by Serbian Ministry of Science and Spanish government have objective to produce wind and solar energy potential atlas of Serbia, but the final results are still not presented. Wind potential map of South East Banat prepared by the team of professor Zlatanovic confirms good wind energy potential of selected region according to measurements performed on 50m height wind mast and on two stations at 33m [4]. Questions related to sensors properties and validation of meteorological data collected from different internet web sites were discussed in references [5,6]. Wind map of former Yugoslavia region using model of the most probable atmospheric state carried out by the ARMINES Institute does not predict very good wind potential of Serbia. 1 Most existing data predict significant wind potential of Vojvodina, the region in which the wind energy was used back in history. The terrain of Vojvodina is very open and flat or gently undulating and corresponds to roughness class 1 (z0=0.03 m) according European wind atlas, except Vršac’s mountain and Fruška Gora that together make no more than 10% of whole area. We present wind atlas of Vojvodina generated by licensed WAsP software tool and linear interpolation method of metstations data [7,8]. The application of WAsP is fully justified since the terrain is mostly homogeneous and accuracy of simulation of wind over this region is accurate. 2. METHOD OF WIND ATLAS GENERATION Publicly available data were used in wind map of Vojvodina preparation and processed according to European Wind Atlas methodology including verification procedure. Several databases were created such as wind characteristics, orography, roughness and obstacles. The internet is selected as the only source of data including digital orography maps, so the reported results may be reproduced by using WAsP software tools and linear interpolation method. The step by step procedure of map generation and verification is described. 2.1 Meteorological data As the member of world network of meteorological stations the Republic Hydrometeorological Service of Serbia sends minimum 8 hours values of meteorological data for international exchange. This data are publicly available for example on the internet site http://meteo.infospace.ru as integer values. The validity of such data is discussed in Ref. [5]. In our procedure we used data from 8 meteorological stations located in the region of Vojvodina. These stations are: Palić, Sombor, Kikinda, Novi Sad, Sremska Mitrovica, Banatski Karlovac, Vršac and Surčin. Data were collected from 3.3.2000 to 3.3.2008 and, in the best case, consist of data recorded every 3 hours. Using the software tool OWC Wizard, which is the part of WAsP Climate Analysis, OWC (Observed Wind Climate) files based on filtered wind speed and direction data were calculated. 2 Fig. 1. Orography map of Vojvodina with the locations of meteorological stations. The color scale represents relative terrain altitude These files represent wind climate data base for certain location, i.e. contain 16 sectors wind rose with mean speed data for each sector and overall mean speed for all sectors. Also they consist of statistical Weibull distribution parameters, like shape parameter k (non-dimensional quantity) and the scale parameter A (dimensions of speed m/s), for each sector as well as overall non directional data. 2.2 Orography For preparation of Vojvodina’s orography map, data were downloaded from NASA’s ’’Shuttle Radar Topography Mission’’ web site, actually its ftp server ftp://e0srp01u.ecs.nasa.gov. After importing these data in to the application Global Mapper, Vojvodina’s orography map was derived as shown in Fig. 1. Obtained orography map were converted in the form that can be used in WAsP software tool by generating contours in required form. Procedure of generating contours was done also in Global Mapper so that contour lines were generated with 10m spacing to make suitable vector terrain map. 2.3 Generating wind atlases For terrain roughness of Vojvodina region generalization was applied. Since this region is a plain which surface is mostly farmed for cereal cultures, roughness class 1 (z0=0.03 m) was used in general. Description of obstacle groups surrounding meteorological stations was made on the basis of satellite pictures using Google Earth. After all of these steps databases were generated and prepared to be imported in WAsP for the purpose of wind atlas generation. Each wind atlas project file in WAsP contains the following layers: OWC file, contour vector map, rose of terrain roughness and obstacle group surrounding the meteorological station. The last step prior to generation of wind atlas in WAsP was locating meteorological station on contour vector map from which data were derived in the form of OWC file. After this last step WAsP was fully prepared for selected region wind atlas calculation. Wind atlases contain the site-independent or regional wind climate, derived from the wind measurements at a meteorological station through the wind atlas analysis. The observed wind climate (OWC) has thereby been reduced to certain standard conditions, i.e. wind roses and wind speed distributions for five standard heights and four roughness classes in a number of sectors (usually 12). Since we had wind data from 8 meteorological stations together with terrain description we derived 8 so called primary or reference wind atlases. 2.4 Interpolation Area of Vojvodina amounts approximately 24 500 km2. In the best case if we have terrain of the lowest roughness, which is equivalent to the calm sea surface, and using the assumption that the terrain is absolutely flat, single wind atlas can be applied on to the area no more than 50km*50km=2500 km2. In our case if we presume that meteorological stations are located ideally, we need at least data from 10 meteorological stations. Even in this idealized case we still have lack of data. Because of this reason, linear interpolation independent of direction was applied between existing wind atlases for the purpose of filling the lack of data for specific areas of Vojvodina. The whole region of Vojvodina was divided into 66 sectors of 20km*20km area covering complete Vojvodina (Fig. 2). In Figure 2 the primary atlas regions are indicated by pink squares. They represent areas in centre of which the meteorological stations are located and wind atlases for these areas were derived from local meteorological stations data. For each sector map beyond the extent of primary atlases, wind atlases were derived using linear interpolation between minimum two, and 3 maximum three nearest wind atlases, which were derived from meteorological stations data and terrain description. The squares of different colors are indicated in Fig. 2 which center represents virtual locations of wind atlases derived using linear interpolation. Fig. 2. Map of Vojvodina divided into 66 sectors of 20km*20km areas. Reference sectors are pink colored (see the text) For the illustration on how interpolated atlases were generated the next example is given. Wind atlas derived from Sombor meteorological station data (UTM coordinates: 356149.9, 5069788) fully describes regional wind climate for area of 20km*20km, which center is in fact location of this meteorological station. However, wind climate of areas surrounding Sombor wind atlas cannot be describe only by this atlas data due to the influence of the other wind atlases. Location SNP1 with UTM coordinates 377835.398, 5070068.243 is also center of an 20km*20km area positioned east of meteorological station Sombor and has the regional wind climate that can be described taking into consideration two more wind atlases beside Sombor atlas (Fig.

Load more