Wave Energy Resource Assessment off the Coast of China around the Zhoushan Islands

Based on ERA-Interim reanalysis wave field data for the 36 years from 1979 to 2014, the temporal and spatial distributions and development potential of wave energy are studied in detail in the offshore and relatively nearshore waters adjacent to the Zhoushan Islands. The results show that areas of relatively high wave energy are located in the offshore waters to the east and southeast of the Zhoushan Islands. The potential wave energy in relatively nearshore waters (water depths from 10 m to 65 m) is relatively higher and especially in the nearshore waters to the southeast of the city of Taizhou, Zhejiang Province, which are suitable locations for wave energy development. The conclusions provide scientific guidance for wave energy development in the sea areas adjacent to the Zhoushan Islands.


Introduction
Ongoing social and economic development is increasing the demand for energy. The energy supply problems arising from a shortage of conventional fossil energy and the pollution problems resulting from the overuse of conventional fossil energy are intensifying. Wave energy is a type of clean and renewable energy that does not produce carbon dioxide. The development of wave energy can effectively ease the energy crisis and environmental pollution, but before the development of wave energy, the reserves and the temporal and spatial distributions of wave energy in interesting areas must be reliably assessed. Many wave energy assessment studies have been conducted in different countries with respect to global sea areas and large-scale sea areas in different regions. Cornett [1] investigated global wave energy resources based on a third-generation model, Wave Watch III (WAVEWATCH-III). Barstow et al. [2] and Mork et al. [3] assessed the wave energy potential of global oceans from deep waters to relatively nearshore waters based on the global wave dataset Worldwaves. Arinaga and Cheung [4] employed 10 years of wave field data from WAVEWATCH-III to evaluate the global wave energy. Zheng et al. [5] evaluated the wind wave energy and swell energy for global oceans in detail based on global wind wave fields and swell fields from European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 reanalysis data over the last 45 years. Pontes et al. [6] and Pontes [7] evaluated the wave energy in Europe's coastal waters in the Atlantic and Mediterranean by using the wave model Wave Modelling project (WAM). Soares et al. [8] assessed the wave energy in detail along the Atlantic European coast by using WAVEWATCH-III and Simulating Waves Nearshore (SWAN). Henfridsson et al. [9] studied the wave energy potential of the Baltic Sea based on a hindcast dataset. Zheng et al. [10] and Zheng et al. [11] comprehensively and systematically assessed the wave energy resources in the China Sea by using high-resolution wave field data that were simulated by the WAVEWATCH-III wave model and obtained the reserves and distribution of the wave energy in a large-scale area in the China Sea. Zheng et al. [12] and Zheng et al. [13] analyzed the wind energy in the global ocean and the wave energy in the China Sea in detail based on Cross Calibrated Multi-Platform (CCMP) wind data and WAVEWATCH-III wave data. Wave energy and wind energy assessments in large-scale areas are essential to understand the temporal and spatial distributions of wave energy and can provide macroscopic references for decision makers to evaluate wave energy projects. However, for practical engineering applications, wave energy assessments in large-scale areas cannot meet the requirements for the siting of wave farms and the design of wave energy converters (WECs). Some wave energy assessments in relatively nearshore, small-scale sea areas were carried out for siting wave farms. Kim et al. [14] assessed the wave energy in offshore deep waters in the sea areas adjacent to Korea based on a hindcast dataset from 1979 to 2003 and evaluated the wave energy in the relatively nearshore waters of Hongdo by using the SWAN model. From 2009 to 2011, Iglesias et al. [15][16][17][18][19][20] conducted a refined evaluation of several small-scale areas in the sea areas adjacent to Spain based on wave model data. Liang et al. [21] assessed the wave energy for two important areas in the adjacent seas to the east of China based with 22 years of wave data by using the SWAN model. Zheng et al. [22] assessed the wave energy in key locations along the relatively nearshore coastline waters of Taiwan based on buoy data from adjacent sea areas in 12 regions of Taiwan. The above works provide good references for wave energy assessment in relatively nearshore, small-scale sea areas.
The region of interest in this study is the sea area adjacent to the Zhoushan Islands, which is shown in Figure 1 (rectangular region in green). The Zhoushan Islands, located in the East China Sea, are the largest islands in China and constitute 20% of the country's islands, with 1390 islands and a population of one million. Fish production is important in the sea areas adjacent to the Zhoushan Islands. In this area, the demand for energy is enormous because of the power that is required by residents and for fish production and sea water desalination. Wave energy is exploitable in the sea areas adjacent to the Zhoushan Islands. Reasonable wave energy development could provide a sustainable energy supply for the development of the Zhoushan Islands and its adjacent regions, however, a comprehensive understanding of the development value of wave energy in the sea areas adjacent to the Zhoushan Islands is lacking. Therefore, based on high-resolution and high-accuracy ERA-Interim reanalysis wave field data from the ECMWF from 1979 to 2014, the temporal and spatial distributions and the development potential of wave energy are studied in the sea areas adjacent to the Zhoushan Islands. The conclusions can provide reliable references for practical engineering applications, such as siting wave farms and the design and operation of WECs in the sea areas adjacent to the Zhoushan Islands.
In the remainder of this paper, the data and their verification are described in Section 2. In Section 3, the distributions of the wave energy are analyzed in the offshore, large-scale areas of these sea areas to help readers understand the overall distribution of the wave energy in the sea areas adjacent to the Zhoushan Islands. In Section 4, we establish a novel criterion for the regional division of wave energy that is suitable for local sea areas and identify the best regions for wave energy in the relatively nearshore waters of the Zhoushan Islands with a focus on relatively nearshore applications. In Section 5, key single locations are selected in key regions, for which the wave energy is analyzed in detail, and important reference information for the practical development of wave energy is described. The last section provides the conclusions.

Data and Their Verification
The significant wave height (Hs, m) and the energy period (Te, s) are two basic parameters for calculating the characteristic quantity of wave energy. The energy period is a mean wave period that can be calculated by m−1/m0 which are spectral moments of the wave spectrum and is related to energy [23]. The wave data in this paper were ERA-Interim reanalysis wave field data over the last 36 years, which were collected by ECMWF with a spatial resolution of 0.125 × 0.125 and a temporal resolution of 6 h. The ECMWF is one of the main data reanalysis centers in the world. The ERA-Interim reanalysis data and its previous edition, ERA-40, have been widely used [5,24,25]. The accuracy of wave energy assessment is directly related to the accuracy of the wave parameters, so we verified the accuracy of the ERA-Interim wave field data based on data from two buoys located in the area of interest. The buoys were Directional Waverider MkIII buoys from Datawell BV (Haarlem, The Netherlands). Directional Waverider buoys are the most advanced buoys available worldwide and dominate the international market for wave buoys. The wave data that were measured by these buoys were averaged values from 10 min of spectral data and were calculated from the spectral moment.
The locations of the buoys are shown in Figure 2, and the specifications of the buoys are shown in Table 1. The locations and observation times of the buoys were different from the ERA-Interim data, so the ERA-Interim data, which were gridded data, were interpolated in space and time to collocate with the buoy data. The spatial resolution and the temporal resolution of the ERA-Interim data were 0.125 × 0.125 and 6 h, respectively, and the locations and observation times of the buoys were fixed values. First, the 4 gridded ERA-Interim point data closest to the buoy were interpolated according to the longitude and latitude of the buoy by using bilinear interpolation to obtain collocated ERA-Interim data at the buoy location, achieving interpolation in space. Then, the ERA-Interim data at the 2 closest times to the observation time of the buoy were interpolated by using linear interpolation to obtain collocated ERA-Interim data at the observation time of the buoy, achieving interpolation in time. The two parameters Hs and Te, were verified; the results are shown in Figure 3, and the error indices are shown in Table 2. The root mean square error (RMSE) and the correlation coefficient (CC) can be calculated as follows:

Data and Their Verification
The significant wave height (H s , m) and the energy period (T e , s) are two basic parameters for calculating the characteristic quantity of wave energy. The energy period is a mean wave period that can be calculated by m −1 /m 0 which are spectral moments of the wave spectrum and is related to energy [23]. The wave data in this paper were ERA-Interim reanalysis wave field data over the last 36 years, which were collected by ECMWF with a spatial resolution of 0.125 × 0.125 and a temporal resolution of 6 h. The ECMWF is one of the main data reanalysis centers in the world. The ERA-Interim reanalysis data and its previous edition, ERA-40, have been widely used [5,24,25]. The accuracy of wave energy assessment is directly related to the accuracy of the wave parameters, so we verified the accuracy of the ERA-Interim wave field data based on data from two buoys located in the area of interest. The buoys were Directional Waverider MkIII buoys from Datawell BV (Haarlem, The Netherlands). Directional Waverider buoys are the most advanced buoys available worldwide and dominate the international market for wave buoys. The wave data that were measured by these buoys were averaged values from 10 min of spectral data and were calculated from the spectral moment.
The locations of the buoys are shown in Figure 2, and the specifications of the buoys are shown in Table 1. The locations and observation times of the buoys were different from the ERA-Interim data, so the ERA-Interim data, which were gridded data, were interpolated in space and time to collocate with the buoy data. The spatial resolution and the temporal resolution of the ERA-Interim data were 0.125 × 0.125 and 6 h, respectively, and the locations and observation times of the buoys were fixed values. First, the 4 gridded ERA-Interim point data closest to the buoy were interpolated according to the longitude and latitude of the buoy by using bilinear interpolation to obtain collocated ERA-Interim data at the buoy location, achieving interpolation in space. Then, the ERA-Interim data at the 2 closest times to the observation time of the buoy were interpolated by using linear interpolation to obtain collocated ERA-Interim data at the observation time of the buoy, achieving interpolation in time. The two parameters H s and T e , were verified; the results are shown in Figure 3, and the error indices are shown in Table 2. The root mean square error (RMSE) and the correlation coefficient (CC) can be calculated as follows: Energies 2017, 10, 1320 where x i denotes the wave parameters from the buoys; y i denotes the wave parameters from the ERA-Interim; x and y denote the mean value of the wave parameters from the buoys and the ERA-Interim respectively; and n is the number of samples. The units of measurement of the RMSE were meters and seconds for H s and T e , respectively. The correlation coefficient is a statistical index to reflect the close degree of correlation between wave parameters from the buoys and those from the ERA-Interim. where denotes the wave parameters from the buoys; denotes the wave parameters from the ERA-Interim; ̅ and denote the mean value of the wave parameters from the buoys and the ERA-Interim respectively; and is the number of samples. The units of measurement of the RMSE were meters and seconds for Hs and Te, respectively. The correlation coefficient is a statistical index to reflect the close degree of correlation between wave parameters from the buoys and those from the ERA-Interim. The RMSE of Hs for buoy 006 was 0.31 m and the CC was 0.90; the RMSE for Te was 0.56 s and the CC was 0.78. The RMSE of Hs for buoy kzszs was 0.36 m and the CC was 0.78; the RMSE for Te was 0.81 s and the CC was 0.70. The results show that the wave field data from the ERA-Interim were accurate in the sea areas adjacent to the Zhoushan Islands and are suitable for wave energy assessment.

Wave Energy Analysis in Offshore Waters
Indices such as the distributions of the significant wave height, energy period, wave power density, Pw (kW/m) and stability of the wave energy (ST) were employed to analyze the temporal and spatial distributions of the wave climate and wave energy in the offshore large-scale areas in the sea areas adjacent to the Zhoushan Islands and can provide macroscopic references for wave energy development.

Temporal and Spatial Distributions of the Wave Climate
Before assessing the wave energy, researchers always analyze the wave climate distribution of the area of interest and to better understand the distribution of the wave energy resources [4,26]. Therefore, we analyze the distributions of the significant wave height and energy period in detail based on the ERA-Interim data from 1979 to 2014.
3.1.1. Temporal and Spatial Distributions of the Significant Wave Height Figure 4 shows the annual average and seasonal average distributions of Hs in the sea areas adjacent to the Zhoushan Islands, which were calculated based on the ERA-Interim data from 1979 to 2014. The results show that the annual average Hs (approximately 0.7-1.6 m) gradually decreased from the offshore deep water areas to the relatively nearshore shallow water areas, which is the direction in which the waves are travelling. Areas with large annual average Hs values (approximately 1.4-1.6 m) were primarily located in the eastern and southeastern offshore sea waters of the Zhoushan Islands, where the sea-states are higher. Areas with small annual average Hs values (below 0.7 m) were primarily located in the relatively nearshore waters of the Changjiang Estuary and Hangzhou Bay, where the sea-states were lower. In this area, the seasons are influenced by the East Asian monsoon climate and extreme sea conditions from the ocean. Winter is from December to February; spring is from March to May; summer is from June to August; and autumn is from September to November. The average Hs showed pronounced seasonal variations in the sea areas adjacent to the Zhoushan Islands. Hs was higher during autumn and winter than during spring and summer because of the monsoon. The sea-states were highest during winter, when Hs reached 2 m; this value was slightly lower during autumn. Hs had similar values during spring and summer.

Wave Energy Analysis in Offshore Waters
Indices such as the distributions of the significant wave height, energy period, wave power density, P w (kW/m) and stability of the wave energy (ST) were employed to analyze the temporal and spatial distributions of the wave climate and wave energy in the offshore large-scale areas in the sea areas adjacent to the Zhoushan Islands and can provide macroscopic references for wave energy development.

Temporal and Spatial Distributions of the Wave Climate
Before assessing the wave energy, researchers always analyze the wave climate distribution of the area of interest and to better understand the distribution of the wave energy resources [4,26]. Therefore, we analyze the distributions of the significant wave height and energy period in detail based on the ERA-Interim data from 1979 to 2014.
3.1.1. Temporal and Spatial Distributions of the Significant Wave Height Figure 4 shows the annual average and seasonal average distributions of H s in the sea areas adjacent to the Zhoushan Islands, which were calculated based on the ERA-Interim data from 1979 to 2014. The results show that the annual average H s (approximately 0.7-1.6 m) gradually decreased from the offshore deep water areas to the relatively nearshore shallow water areas, which is the direction in which the waves are travelling. Areas with large annual average H s values (approximately 1.4-1.6 m) were primarily located in the eastern and southeastern offshore sea waters of the Zhoushan Islands, where the sea-states are higher. Areas with small annual average H s values (below 0.7 m) were primarily located in the relatively nearshore waters of the Changjiang Estuary and Hangzhou Bay, where the sea-states were lower. In this area, the seasons are influenced by the East Asian monsoon climate and extreme sea conditions from the ocean. Winter is from December to February; spring is from March to May; summer is from June to August; and autumn is from September to November. The average H s showed pronounced seasonal variations in the sea areas adjacent to the Zhoushan Islands. H s was higher during autumn and winter than during spring and summer because of the monsoon. The sea-states were highest during winter, when H s reached 2 m; this value was slightly lower during autumn. H s had similar values during spring and summer.     Figure 5 shows the annual average and seasonal average distributions of Te in the sea areas adjacent to the Zhoushan Islands, which were calculated based on the ERA-Interim data from 1979 to 2014. The results show that the annual average Te values were relatively small (approximately 4.0-5.5 s) in the relatively nearshore waters and were relatively large (approximately 5.5-6.5 s) in the offshore waters. Te was steady in most of the sea areas. The results of the seasonal average show that the seasonal variation in Te was not obvious: Te was approximately 5.5-7.0 s in most of the sea areas and was slightly higher during winter and autumn. The seasonal spatial distributions of Te were similar to the annual average distributions. A relatively steady Te is advantageous for wave energy stability.  Figure 5 shows the annual average and seasonal average distributions of T e in the sea areas adjacent to the Zhoushan Islands, which were calculated based on the ERA-Interim data from 1979 to 2014. The results show that the annual average T e values were relatively small (approximately 4.0-5.5 s) in the relatively nearshore waters and were relatively large (approximately 5.5-6.5 s) in the offshore waters. T e was steady in most of the sea areas. The results of the seasonal average show that the seasonal variation in T e was not obvious: T e was approximately 5.5-7.0 s in most of the sea areas and was slightly higher during winter and autumn. The seasonal spatial distributions of T e were similar to the annual average distributions. A relatively steady T e is advantageous for wave energy stability.

Calculation Method for the Wave Power Density
The wave power density (Pw) is the most important index for wave energy assessment and is related to the accuracy of the wave energy assessment. The water depth is shallow in the sea areas adjacent to the Zhoushan Islands. The influence of the water depth must be considered when calculating Pw to improve the accuracy of Pw because of shoaling. In this paper, we considered the influence of the water depth when calculating Pw. The wave power density for all water depths can be calculated as follows [27]: where = g is the wave energy density (J/m 2 ), is the density of sea water (kg/m 3 ), g is the acceleration from gravity (m/s 2 ), = is the wavenumber (m −1 ), is the wavelength (m), d is the water depth (m), and * = (1 + ).

Calculation Method for the Wave Power Density
The wave power density (P w ) is the most important index for wave energy assessment and is related to the accuracy of the wave energy assessment. The water depth is shallow in the sea areas adjacent to the Zhoushan Islands. The influence of the water depth must be considered when calculating P w to improve the accuracy of P w because of shoaling. In this paper, we considered the influence of the water depth when calculating P w. The wave power density for all water depths can be calculated as follows [27]: where E = 1 16 ρgH 2 s is the wave energy density (J/m 2 ), ρ is the density of sea water (kg/m 3 ), g is the acceleration from gravity (m/s 2 ), k = 2π λ is the wavenumber (m −1 ), λ is the wavelength (m), d is the water depth (m), and P * = 1 2 (1 + 2kd sin h2kd ). The bathymetry data were obtained from the University of California San Diego (UCSD) website and had a spatial resolution of 1 × 1 [28]. λ refers to the average wavelength according to the wave number and is estimated by using the mean wave period; to obtain accurate results, λ was calculated according to a relationship that was presented by Wen et al. [29], which was obtained by experimentation. λ can be calculated as follows: In deep water, the wave power density can often be calculated as follows: Equation (5) is the formula for deep water. This formula has been widely adopted [5][6][7]10,11,24,30] and has high accuracy, which can reduce the computational efforts.

Temporal and Spatial Distributions of the Wave Power Density
The annual average and seasonal average distributions of P w in the sea areas adjacent to the Zhoushan Islands, which were calculated based on the ERA-Interim data from 1979 to 2014, are shown in Figure 6. According to the results of the annual average, P w (approximately 2-13 kW/m) gradually decreased from the offshore deep water areas to the relatively nearshore shallow water areas because the wave power declines as the resource moves into shallow water. Areas with relatively large annual average P w values (approximately 10-13 kW/m) were primarily located in the eastern and southeastern offshore sea waters of the Zhoushan Islands, where wave energy was relatively abundant. Areas with small annual average P w values (below 2 kW/m) were primarily located in the relatively nearshore waters of the Changjiang Estuary and Hangzhou Bay, where the development value of wave energy was low. It is very important to note that compared with global sea areas [31], the Zhoushan Islands coastal average wave energy resource is weak as it only reaches 10 kW/m nowhere along the coastline, the wave resource investigated here is weak in comparison to the global resource. So we declare the results discussed in this paper only aim at local sea areas around the Zhoushan Islands and want to provide references for wave farm construction and siting in this local areas. All conclusions in the next sections about abundant level of wave energy are relative to local sea areas around the Zhoushan Islands not relative to global scale. Pronounced seasonal variations in the wave energy were present in the sea areas adjacent to the Zhoushan Islands because of the East Asian monsoon climate and extreme sea conditions from the ocean. P w was higher during autumn and winter than during spring and summer because of this monsoon climate. Wave energy was most abundant in winter, when P w reached 18 kW/m. P w was slightly lower during autumn than during winter. Spring and summer had similar P w values, although P w was slightly higher during summer, when typhoons from the West Pacific influence this area. Wave energy was usable, and has certain development value in the sea areas adjacent to the Zhoushan Islands, with autumn and winter being the key periods for energy extraction.

Stability of Wave Energy
Several studies have shown that we also must consider the stability when assessing wave energy resources. The stability of wave energy means that the wave energy is sustainable for a long time and the variation in the wave energy is small. Abundant and stable wave energy is advantageous for exploitation and utilization [1]. Therefore, studying the stability of the wave energy in areas of interest is important. The stability index (ST) can express the stability of wave energy and can be calculated as follows [24]: where _ is the mean wave power density in the minimum energy month (kW/m) and is the annual mean wave power density over the last 36 years (kW/m). Equation (6) shows that the mean wave power density in the minimum energy month is closer to the annual mean wave power density when ST is higher, indicating less variation in the wave power density. Under this circumstance, the wave energy is more stable and sustainable. Figure 7 shows the distributions of ST in the sea areas adjacent to the Zhoushan Islands, which were calculated based on ERA-Interim data over the last 36 years. This figure shows that the stability was relatively higher along relatively nearshore coastline areas, with the ST values all above 0.45; the stability was relatively lower in offshore sea areas, with the ST values are below 0.45. The stability further decreased in the sea areas to the east of the

Stability of Wave Energy
Several studies have shown that we also must consider the stability when assessing wave energy resources. The stability of wave energy means that the wave energy is sustainable for a long time and the variation in the wave energy is small. Abundant and stable wave energy is advantageous for exploitation and utilization [1]. Therefore, studying the stability of the wave energy in areas of interest is important. The stability index (ST) can express the stability of wave energy and can be calculated as follows [24]: where P M_min is the mean wave power density in the minimum energy month (kW/m) and P year is the annual mean wave power density over the last 36 years (kW/m). Equation (6) shows that the mean wave power density in the minimum energy month is closer to the annual mean wave power density when ST is higher, indicating less variation in the wave power density. Under this circumstance, the wave energy is more stable and sustainable. Figure 7 shows the distributions of ST in the sea areas adjacent to the Zhoushan Islands, which were calculated based on ERA-Interim data over the last 36 years. This figure shows that the stability was relatively higher along relatively nearshore coastline areas, with the ST values all above 0.45; the stability was relatively lower in offshore sea areas, with the

Discussion on Wave Energy in Offshore Waters
In Section 3, we analyzed the characteristics of the wave energy in the offshore waters in the sea areas adjacent to the Zhoushan Islands in detail. The results showed that wave energy is relatively weak in the sea areas adjacent to the Zhoushan Islands compared with global ocean, for the most part less than 10 kW/m annually. Although the wave energy is not abundant, it is usable and has certain potential for development. In the sea areas around the Zhoushan Islands, the areas with relatively high annual average wave power density (approximately 10-13 kW/m) are primarily located in the eastern and southeastern offshore sea waters of the Zhoushan Islands. The wave power density is higher during autumn and winter than during spring and summer because of the monsoon climate. Autumn and winter are the key periods for energy extraction. In addition, the wave energy is relatively more stable and sustainable in the relatively nearshore waters than in the offshore waters. However, the relatively nearshore waters have less wave energy than the offshore waters. The utilization and development of wave energy are mainly located in relatively nearshore waters because of the current limits of technology and the costs of wave farms. Therefore, the relatively nearshore sea areas to the east of the Zhoushan Islands, which have an annual average wave power density of 5-8 kW/m and a stability index of 0.4-0.45, are the most suitable areas for wave energy development in this areas.

Division of Key Wave Energy Regions in Relatively Nearshore Waters
At present, the development of wave energy is mainly located in relatively nearshore local regions because of the technological limits and the cost of development. Accurately evaluating wave energy in relatively nearshore local regions is very important to support practical engineering in the development of wave energy. A single location mode is usually adopted in relatively nearshore wave energy assessment. That is, we can select locations (fixed locations) in relatively nearshore waters and study the distributions of the wave energy in each location. Before selecting key single locations in relatively nearshore waters, we must first carry out a regional division for wave energy in relatively

Discussion on Wave Energy in Offshore Waters
In Section 3, we analyzed the characteristics of the wave energy in the offshore waters in the sea areas adjacent to the Zhoushan Islands in detail. The results showed that wave energy is relatively weak in the sea areas adjacent to the Zhoushan Islands compared with global ocean, for the most part less than 10 kW/m annually. Although the wave energy is not abundant, it is usable and has certain potential for development. In the sea areas around the Zhoushan Islands, the areas with relatively high annual average wave power density (approximately 10-13 kW/m) are primarily located in the eastern and southeastern offshore sea waters of the Zhoushan Islands. The wave power density is higher during autumn and winter than during spring and summer because of the monsoon climate. Autumn and winter are the key periods for energy extraction. In addition, the wave energy is relatively more stable and sustainable in the relatively nearshore waters than in the offshore waters. However, the relatively nearshore waters have less wave energy than the offshore waters. The utilization and development of wave energy are mainly located in relatively nearshore waters because of the current limits of technology and the costs of wave farms. Therefore, the relatively nearshore sea areas to the east of the Zhoushan Islands, which have an annual average wave power density of 5-8 kW/m and a stability index of 0.4-0.45, are the most suitable areas for wave energy development in this areas.

Division of Key Wave Energy Regions in Relatively Nearshore Waters
At present, the development of wave energy is mainly located in relatively nearshore local regions because of the technological limits and the cost of development. Accurately evaluating wave energy in relatively nearshore local regions is very important to support practical engineering in the development of wave energy. A single location mode is usually adopted in relatively nearshore wave energy assessment. That is, we can select locations (fixed locations) in relatively nearshore waters and study the distributions of the wave energy in each location. Before selecting key single locations in relatively nearshore waters, we must first carry out a regional division for wave energy in relatively nearshore waters and identify key regions of wave energy. Then, we can select key locations for studies in these regions.

Division Criterion
A novel criterion for regional division that is suitable for local sea areas was established. This criterion is based on three parameters: the annual average wave power density (P w , kW/m), annual average effective time (T E , h) and total wave energy (TWE, J). The annual average effective time is the annual average duration for the wave energy, with 1 m ≤ H s ≤ 4 m. The reason is wave energy converters will possibly be damaged and have higher cost for maintenance in high sea states, for example Wave Star will operate in storm protection mode and stop generating power when significant wave height is above 3.0 m [32]. At the same time, when significant wave height is less than 1 m, it may not be possible to economically harvest energy with present state WECs. The total wave energy can be calculated by using two main methods: the synoptic method and climatology method [33,34]. The climatology method is the better choice because it uses a long-term average, which can represent a long-term trend. The regional division criterion is established as follows: (1) Nine small continuous areas that were 0.5 × 0.5 in size were selected in the relatively nearshore waters of the Zhoushan Islands, which are shown in Figure 8. First, the annual average P w and annual average T E of every grid were calculated for each small area. Then, the results of every grid in each small area were averaged, providing the annual average P w and the annual average T E of each small area. Finally, the TWE for each small area was calculated. (2) The comprehensive division coefficient (CDC) was defined as follows: When the CDC is higher, wave energy is abundant and the potential of development is higher. (3) Five grades were established according to the CDC. The general criterion of the regional division is shown in Table 3. The potential of the wave energy increases from levels 1 to 5, which indicate poor, available, good, better and best potential. The ranges of each index in Table 3 exhibit different value ranges for the different areas. The threshold value a-c can be determined by equal division as follows (with P w as an example): annual average P w interval ( P w _interval) = P wmaxaa − P wminaa 5 Then: a1 = P wminaa + P w _interval (9) a2 = P wminaa + 2 × P w _interval (10) where P wmaxaa is the maximum annual average P w and P wminaa is the minimum annual average P w .
The method for determining the other threshold values is similar to that of P w . (4) Thresholds were calculated for each index, as shown in Table 4.  (5) The criterion of the regional division was established. This criterion is suitable for the relatively nearshore waters of the Zhoushan Islands. The results are shown in Table 5. (4) Thresholds were calculated for each index, as shown in Table 4.  (5) The criterion of the regional division was established. This criterion is suitable for the relatively nearshore waters of the Zhoushan Islands. The results are shown in Table 5.

Division Results
We calculated the CDC and performed a regional division for nine small areas with a size of 0.5 × 0.5 in the relatively nearshore waters of the Zhoushan Islands according to the regional division criterion in Section 4.1; the results are shown in Figure 8. This figure shows that the grade of area A4 was one (CDC of 3518), showing poor development potential. The grades of areas A1 and A8 were two (CDCs of 12,400 and 9901, respectively), showing available development potential. The grades of areas A2, A3 and A7 were three (CDCs of 23,213, 16,324 and 18,697, respectively), showing good development potential. The grades of areas A5, A6 and A9 were five (CDCs of 54,187, 54,342 and 54,084, respectively), showing the best development potential. In this paper, we selected areas A5, A6 and A9 as the key areas for wave energy development in these relatively nearshore waters.

Division Results
We calculated the CDC and performed a regional division for nine small areas with a size of 0.5 × 0.5 in the relatively nearshore waters of the Zhoushan Islands according to the regional division criterion in Section 4.1; the results are shown in Figure 8. This figure shows that the grade of area A4 was one (CDC of 3518), showing poor development potential. The grades of areas A1 and A8 were two (CDCs of 12,400 and 9901, respectively), showing available development potential. The grades of areas A2, A3 and A7 were three (CDCs of 23,213, 16,324 and 18,697, respectively), showing good development potential. The grades of areas A5, A6 and A9 were five (CDCs of 54,187, 54,342 and 54,084, respectively), showing the best development potential. In this paper, we selected areas A5, A6 and A9 as the key areas for wave energy development in these relatively nearshore waters.

Wave Energy Analysis in the Key Locations in the Relatively Nearshore Waters
We selected key locations after determining the key areas. Indices that involved siting wave farms and the design of WECs were analyzed to provide reliable references for engineering in the development of wave energy in relatively nearshore waters.

Selecting Key Locations in the Relatively Nearshore Waters
We referred to the CDC in Section 4.1 when selecting key locations in these relatively nearshore waters and calculated the CDC for every grid in three key areas, where the TWE was the total wave energy in the range of 0.125 × 0.125 around each grid. The results are shown in Figure 9 (CDC contour plot). Figure 9 shows that the locations in each key area with the largest CDC, which denotes the highest potential for development, were selected as the key locations. The information of each location is shown in Figures 8 and 9 and Table 6. These three key locations are S1 (123.250 • E, 30.

Wave Energy Analysis in the Key Locations in the Relatively Nearshore Waters
We selected key locations after determining the key areas. Indices that involved siting wave farms and the design of WECs were analyzed to provide reliable references for engineering in the development of wave energy in relatively nearshore waters.

Selecting Key Locations in the Relatively Nearshore Waters
We referred to the CDC in Section 4.1 when selecting key locations in these relatively nearshore waters and calculated the CDC for every grid in three key areas, where the TWE was the total wave energy in the range of 0.125 × 0.125 around each grid. The results are shown in Figure 9 (CDC contour plot). Figure 9 shows that the locations in each key area with the largest CDC, which denotes the highest potential for development, were selected as the key locations. The information of each location is shown in Figures 8 and 9 and   Next, the annual average Pw was calculated for each key location based on the ERA-Interim data from the last 36 years, and the results are shown in Table 6. According to the results, the annual average Pw was 6.53 kW/m or higher in the three key locations and was higher than the criterion of the exploitable wave energy (Pw ≥ 2 kW/m) [35]. The results show that the wave energy has certain potential for development in the relatively nearshore waters around the Zhoushan Islands.

Directionality of the Wave Energy's Propagation
The directionality of the wave energy's propagation is very meaningful to improve the conversion efficiency of WECs. Generally, the incident wave, which is orthometric with WECs, has the highest absorbance efficiency [27]. The directionality of the wave energy's propagation is only of interest for certain WECs, such as oscillating wave surge converters (OWSCs), but not for axisymmetric devices. In addition, the wave energy in extreme sea-states can never be exploited.  Next, the annual average P w was calculated for each key location based on the ERA-Interim data from the last 36 years, and the results are shown in Table 6. According to the results, the annual average P w was 6.53 kW/m or higher in the three key locations and was higher than the criterion of the exploitable wave energy (P w ≥ 2 kW/m) [35]. The results show that the wave energy has certain potential for development in the relatively nearshore waters around the Zhoushan Islands.

Directionality of the Wave Energy's Propagation
The directionality of the wave energy's propagation is very meaningful to improve the conversion efficiency of WECs. Generally, the incident wave, which is orthometric with WECs, has the highest absorbance efficiency [27]. The directionality of the wave energy's propagation is only of interest for certain WECs, such as oscillating wave surge converters (OWSCs), but not for axi-symmetric devices. In addition, the wave energy in extreme sea-states can never be exploited. Therefore, we must consider the influence of extreme sea-states when analyzing the directionality of exploitable wave energy. This study examined the directionality of the wave energy's propagation under operational and extreme sea-states for key locations in the relatively nearshore waters of the Zhoushan Islands based on ERA-Interim data over the last 36 years for certain WECs. An operational sea-state is the sea-state with H s < 4 m, indicating that the wave energy is exploitable, and an extreme sea-state is a sea-state with H s ≥ 4 m, indicating that the wave energy is not exploitable. The percentage of the energy distribution in each direction was calculated separately for operational and extreme sea-states, and the wave power roses are plotted in Figures 10 and 11. The direction of the dominant wave power for the operational sea-states at Location S1 was NNW-NE (clockwise direction), and the wave power accounted for approximately 63.08% of the total exploitable wave energy. At Location S2, the direction of the dominant wave power was NNW-NE, and the wave power accounted for approximately 62.83% of the total exploitable wave energy. At Location S3, the direction of the dominant wave power was N-ENE, and the wave power accounted for approximately 71.67% of the total exploitable wave energy. The directionality of the exploitable wave energy's propagation was basically consistent in the key locations in the relatively nearshore waters of the Zhoushan Islands. The direction of the dominant wave power was mainly in the N and NNE directions because the influence of the northeastern monsoon is prominent in these areas; the waves are rough in the northeastern direction, and wave energy is relatively abundant. These results can serve as a reference to determine the operational direction of WECs. The direction of extreme waves for extreme sea-states at Location S1 was N-SE, and the wave power accounted for approximately 86.26% of the total extreme wave energy. At Location S2, the direction of extreme wave was N-SE, and the wave power accounted for approximately 87.68% of the total extreme wave energy. At Location S3, the direction of extreme wave was NNE-SSE, and the wave power accounted for approximately 91.80% of the total extreme wave energy. These results show that extreme waves mainly originated from the N-SSE directions because of cold swells from northeastern monsoons during winter and storms from the western Pacific during summer. This finding is of interest to developers. Therefore, we must consider the influence of extreme sea-states when analyzing the directionality of exploitable wave energy. This study examined the directionality of the wave energy's propagation under operational and extreme sea-states for key locations in the relatively nearshore waters of the Zhoushan Islands based on ERA-Interim data over the last 36 years for certain WECs. An operational sea-state is the sea-state with Hs < 4 m, indicating that the wave energy is exploitable, and an extreme sea-state is a sea-state with Hs ≥ 4 m, indicating that the wave energy is not exploitable. The percentage of the energy distribution in each direction was calculated separately for operational and extreme sea-states, and the wave power roses are plotted in Figures 10 and 11. The direction of the dominant wave power for the operational sea-states at Location S1 was NNW-NE (clockwise direction), and the wave power accounted for approximately 63.08% of the total exploitable wave energy. At Location S2, the direction of the dominant wave power was NNW-NE, and the wave power accounted for approximately 62.83% of the total exploitable wave energy. At Location S3, the direction of the dominant wave power was N-ENE, and the wave power accounted for approximately 71.67% of the total exploitable wave energy. The directionality of the exploitable wave energy's propagation was basically consistent in the key locations in the relatively nearshore waters of the Zhoushan Islands. The direction of the dominant wave power was mainly in the N and NNE directions because the influence of the northeastern monsoon is prominent in these areas; the waves are rough in the northeastern direction, and wave energy is relatively abundant. These results can serve as a reference to determine the operational direction of WECs. The direction of extreme waves for extreme sea-states at Location S1 was N-SE, and the wave power accounted for approximately 86.26% of the total extreme wave energy. At Location S2, the direction of extreme wave was N-SE, and the wave power accounted for approximately 87.68% of the total extreme wave energy. At Location S3, the direction of extreme wave was NNE-SSE, and the wave power accounted for approximately 91.80% of the total extreme wave energy. These results show that extreme waves mainly originated from the N-SSE directions because of cold swells from northeastern monsoons during winter and storms from the western Pacific during summer. This finding is of interest to developers.  Therefore, we must consider the influence of extreme sea-states when analyzing the directionality of exploitable wave energy. This study examined the directionality of the wave energy's propagation under operational and extreme sea-states for key locations in the relatively nearshore waters of the Zhoushan Islands based on ERA-Interim data over the last 36 years for certain WECs. An operational sea-state is the sea-state with Hs < 4 m, indicating that the wave energy is exploitable, and an extreme sea-state is a sea-state with Hs ≥ 4 m, indicating that the wave energy is not exploitable. The percentage of the energy distribution in each direction was calculated separately for operational and extreme sea-states, and the wave power roses are plotted in Figures 10 and 11. The direction of the dominant wave power for the operational sea-states at Location S1 was NNW-NE (clockwise direction), and the wave power accounted for approximately 63.08% of the total exploitable wave energy. At Location S2, the direction of the dominant wave power was NNW-NE, and the wave power accounted for approximately 62.83% of the total exploitable wave energy. At Location S3, the direction of the dominant wave power was N-ENE, and the wave power accounted for approximately 71.67% of the total exploitable wave energy. The directionality of the exploitable wave energy's propagation was basically consistent in the key locations in the relatively nearshore waters of the Zhoushan Islands. The direction of the dominant wave power was mainly in the N and NNE directions because the influence of the northeastern monsoon is prominent in these areas; the waves are rough in the northeastern direction, and wave energy is relatively abundant. These results can serve as a reference to determine the operational direction of WECs. The direction of extreme waves for extreme sea-states at Location S1 was N-SE, and the wave power accounted for approximately 86.26% of the total extreme wave energy. At Location S2, the direction of extreme wave was N-SE, and the wave power accounted for approximately 87.68% of the total extreme wave energy. At Location S3, the direction of extreme wave was NNE-SSE, and the wave power accounted for approximately 91.80% of the total extreme wave energy. These results show that extreme waves mainly originated from the N-SSE directions because of cold swells from northeastern monsoons during winter and storms from the western Pacific during summer. This finding is of interest to developers.

Distribution of the Wave Energy Density According to the Wave Condition
When designing WECs, we must consider the range of wave conditions to improve the conversion efficiency. These conditions are the dominant contribution to the total wave energy density (the mean of the sum of P w multiplied by the number of hours in a year MWh/(m·year)) and can be an important index for the design of WECs [24]. Based on the ERA-Interim data over the last 36 years by using an H s interval of 0.5 m and a T e interval of 1 s, we calculated the percentage of the total wave energy density under each wave condition to the total wave energy density under all wave conditions; the results are shown in Figure 12. In this figure, the color code denotes the percentage and the number denotes the annual average number of times for each wave condition. At Locations S1, S2 and S3, the dominant wave conditions all had H s values of 0.5-4.0 m and T e values of 4-9 s (the dominant wave condition is the range of wave conditions in which the wave energy accounts for 85% of the total wave energy density). Within this wave state, the wave energy accounted for 92.58%, 90.73% and 85.96% of the total wave energy density, respectively. The dominant wave conditions in the key relatively nearshore locations were the same and thus can act as a uniform standard for the design of WECs. WECs that are designed according to this standard can achieve higher conversion efficiency. In addition, the number of occurrences of wave condition with an H s value of 0.5-1.0 m and a T e value of 5-6 s at Locations S1 and S2 was the highest, occurring 289 and 274 times and accounting for 4.95% and 4.50% of the total wave energy density, respectively. At Location S3, the number of occurrences of the wave condition with an H s value of 1.0-1.5 m and a T e value of 5-6 s was the highest, occurring 211 times and accounting for 6.63% of the total wave energy density. However, the ratio of the wave energy under the most common wave conditions to the total wave energy density was not the highest. At Locations S1, S2 and S3, the wave condition with the highest ratio of the total wave energy density was that with an H s value of 2.0-2.5 m and a T e value of 6-7 s. At each location, the annual average occurrences of the wave condition with the highest ratio of the total wave energy density were 70, 74 and 88, accounting for 11.84%, 11.67% and 11.37% of the total wave energy density, respectively. These results show that the contribution of the wave conditions to the most frequent ratio was not the most important factor.

Distribution of the Wave Energy Density According to the Wave Condition
When designing WECs, we must consider the range of wave conditions to improve the conversion efficiency. These conditions are the dominant contribution to the total wave energy density (the mean of the sum of Pw multiplied by the number of hours in a year MWh/(m·year)) and can be an important index for the design of WECs [24]. Based on the ERA-Interim data over the last 36 years by using an Hs interval of 0.5 m and a Te interval of 1 s, we calculated the percentage of the total wave energy density under each wave condition to the total wave energy density under all wave conditions; the results are shown in Figure 12. In this figure, the color code denotes the percentage and the number denotes the annual average number of times for each wave condition. At Locations S1, S2 and S3, the dominant wave conditions all had Hs values of 0.5-4.0 m and Te values of 4-9 s (the dominant wave condition is the range of wave conditions in which the wave energy accounts for 85% of the total wave energy density). Within this wave state, the wave energy accounted for 92.58%, 90.73% and 85.96% of the total wave energy density, respectively. The dominant wave conditions in the key relatively nearshore locations were the same and thus can act as a uniform standard for the design of WECs. WECs that are designed according to this standard can achieve higher conversion efficiency. In addition, the number of occurrences of wave condition with an Hs value of 0.5-1.0 m and a Te value of 5-6 s at Locations S1 and S2 was the highest, occurring 289 and 274 times and accounting for 4.95% and 4.50% of the total wave energy density, respectively. At Location S3, the number of occurrences of the wave condition with an Hs value of 1.0-1.5 m and a Te value of 5-6 s was the highest, occurring 211 times and accounting for 6.63% of the total wave energy density. However, the ratio of the wave energy under the most common wave conditions to the total wave energy density was not the highest. At Locations S1, S2 and S3, the wave condition with the highest ratio of the total wave energy density was that with an Hs value of 2.0-2.5 m and a Te value of 6-7 s. At each location, the annual average occurrences of the wave condition with the highest ratio of the total wave energy density were 70, 74 and 88, accounting for 11.84%, 11.67% and 11.37% of the total wave energy density, respectively. These results show that the contribution of the wave conditions to the most frequent ratio was not the most important factor.

Inter-Annual Variation in the Total Wave Energy
The inter-annual variation in the total wave energy can effectively express the duration and stability of the wave energy over the long term at each location. This index is important to reasonably position WECs. Based on ERA-Interim data over the last 36 years, the total wave energy in the range of 0.125 × 0.125 was calculated for each year in the relatively nearshore key locations. The inter-annual variation in the total wave energy for each location is shown in Figure 13. This figure shows that the

Inter-Annual Variation in the Total Wave Energy
The inter-annual variation in the total wave energy can effectively express the duration and stability of the wave energy over the long term at each location. This index is important to reasonably position WECs. Based on ERA-Interim data over the last 36 years, the total wave energy in the range of 0.125 × 0.125 was calculated for each year in the relatively nearshore key locations. The inter-annual Energies 2017, 10, 1320 20 of 25 variation in the total wave energy for each location is shown in Figure 13. This figure shows that the total wave energy of Locations S1, S2 and S3 in each year is higher than 1.4 × 10 11 , 1.5 × 10 11 , and 1.7 × 10 11 , respectively. The wave energy in each location and year was sustainable, and the variation in the wave energy was small, which is advantageous to the development of wave energy projects.
Energies 2017, 10, 1320 20 of 25 total wave energy of Locations S1, S2 and S3 in each year is higher than 1.4 × 10 11 , 1.5 × 10 11 , and 1.7 × 10 11 , respectively. The wave energy in each location and year was sustainable, and the variation in the wave energy was small, which is advantageous to the development of wave energy projects.

Survivability and Use Ratio for Wave Energy Development
We must consider the survivability and practical use ratio for wave energy development in addition to the abundance level, the direction distributions, the wave energy distributions according to the wave condition, and the total wave energy.
The survivability refers to the allowable maximum sea state in which a WEC can operate normally. The maximum sea state must be defined by WEC developers so that they can design their machines to withstand, but this has a cost. When the sea state is higher, the cost of designing a WEC for survivability is also higher. The maximum sea state can be expressed by the maximum wave power density (Pmax, kW/m). When Pmax is higher, the cost of designing a WEC for survivability is also higher. Based on ERA-Interim data over the last 36 years, Pmax was calculated in the relatively nearshore key locations. The results are shown in Table 6. At each location, Pmax showed small differences, and the maximum Pmax was 385.30 kW/m, which can be used as a standard for developers to design WECs for survivability. The results in this paper can be a reference for the design of WECs in the relatively nearshore key locations of the Zhoushan Islands.
The practical use ratio is the percentage of the energy of the dominant wave condition that can be effectively converted by the WECs designed according to the standard in Section 5.3. The practical use ratio for the wave energy can be expressed as the ratio of the wave energy in the dominant wave condition to the total wave energy, where the dominant wave condition is defined as in Section 5.3. Based on the ERA-Interim data over the last 36 years, this value was calculated for each location, and the results are shown in Section 5.3. The results show that at each location, the practical use ratios were all greater than 85%. Thus, the practical use ratio for wave energy is high, and the development value of the wave energy is high in the key relatively nearshore locations of the Zhoushan Islands.

Survivability and Use Ratio for Wave Energy Development
We must consider the survivability and practical use ratio for wave energy development in addition to the abundance level, the direction distributions, the wave energy distributions according to the wave condition, and the total wave energy.
The survivability refers to the allowable maximum sea state in which a WEC can operate normally. The maximum sea state must be defined by WEC developers so that they can design their machines to withstand, but this has a cost. When the sea state is higher, the cost of designing a WEC for survivability is also higher. The maximum sea state can be expressed by the maximum wave power density (P max , kW/m). When P max is higher, the cost of designing a WEC for survivability is also higher. Based on ERA-Interim data over the last 36 years, P max was calculated in the relatively nearshore key locations. The results are shown in Table 6. At each location, P max showed small differences, and the maximum P max was 385.30 kW/m, which can be used as a standard for developers to design WECs for survivability. The results in this paper can be a reference for the design of WECs in the relatively nearshore key locations of the Zhoushan Islands.
The practical use ratio is the percentage of the energy of the dominant wave condition that can be effectively converted by the WECs designed according to the standard in Section 5.3. The practical use ratio for the wave energy can be expressed as the ratio of the wave energy in the dominant wave condition to the total wave energy, where the dominant wave condition is defined as in Section 5.3. Based on the ERA-Interim data over the last 36 years, this value was calculated for each location, and the results are shown in Section 5.3. The results show that at each location, the practical use ratios were all greater than 85%. Thus, the practical use ratio for wave energy is high, and the development value of the wave energy is high in the key relatively nearshore locations of the Zhoushan Islands.

Performance Assessment of Wave Energy Converters
Compared to theoretical wave energy, decision makers that evaluate wave energy projects care more about the conditions of exploitable wave energy, which is more valuable for wave energy development. In this field, Folley and Whittaker [36] indicated that the exploitable wave energy resources in relatively nearshore waters can be defined by the average wave power that propagates in a fixed direction and is limited to a multiple of four times the average wave power. This work can provide reference for exploitable wave energy assessment. In Sections 5.2-5.5, we evaluated some characteristics of the theoretical wave energy, which may be useful for the design, placement and operation of WECs in key locations. However Rusu and Soares [37] indicated that the actual wave power yield depends on the particular WEC device as rated by its wave power generation because each technology has different operational ranges and different efficiencies under various sea states. We must estimate the energy production of some proposed state of the art WECs to develop a pilot plant of wave energy in these key locations, which may be more useful for the practical development of wave energy. In this paper, we refer to the work of Vannucchi and Cappietti [38], who evaluated the performances of six different state of the art WECs to find which device was the best suited for Italian areas. This work is meaningful for wave energy development.
In this paper, we selected five different state of the art WECs that were used by Vannucchi and Cappietti [38], including AquaBuoy (>50 m) [39], AWS (>50 m) [40], Pelamis (>50 m) [41], Oyster (≈15 m) [42] and Wave Star (Nearshore) [43], by considering the water depths of S1, S2 and S3 (62 m, 60 m, 16 m, respectively) and the suitable water depths of these WECs. The main features of these WECs are shown in Vannucchi and Cappietti [38]. Then, we employed three indices, including the electricity production (P e , kW), capacity factor (C f , %) and capture width (C w , m), which were indicated by Vannucchi and Cappietti [38] to analyze the energy production and determine the best suited WEC in each key location off the coast of China around the Zhoushan Islands. The formulas of these three indices are shown in Equations (13)-(15) [38]: where nT is the number of T e classes and nH is the number of H s classes. P ij is the power output of the device, which originates from the power matrix for the specific WEC. The power matrices of the WECs that were selected in this paper were collected from the open literature [32,44,45]. f ij is the occurrence frequency for various sea states, which was collected from scatter diagrams [37]. The rated power (unit: kW) is the rated power of the WEC. The rated power of the WECs in this paper were 250 kW (AcquaBuoy); 2000 kW (AWS); 750 kW (Pelamis); 800 kW (Oyster) and 600 kW (Wave Star). P w (unit: kW/m) is the annual average flux of energy in each key location. It is very important to note that because different devices have vastly different size, the capture width is not a good measure for comparison. So we employed the relative capture width (R cw , %) instead of the capture width to compare the performance of WECs. The relative capture width is defined as Equation (16): where main_dimension is a main dimension of WECs. A main dimension of AcquaBuoy is about 20 m [39]; AWS is 144 m [46]; Pelamis is 180 m [47]; Oyster is 26 m [42]; WS is 70 m [32]. Based on the above methodology and the ERA-Interim reanalysis wave field data from 1979 to 2014, we calculated the mean power output, capacity factor and relative capture width in each key location. The results are shown in Table 7. The highest capacity value was obtained for Pelamis in Location S1, where the capacity factor was 5.09%, the mean power output was 38.16 kW, but the relative capture width was 3.24%. The best relative capture width in Location S1 was obtained for AB with 8.05%. The results in Location S2 were similar to those in Location S1. Locations S1 and S2 were located in the relatively nearshore waters to the east of the Zhoushan Islands, where the water depth are approximately 60 m, so Pelamis looks like the best suited WEC when considering capacity factor and AB looks like the best suited WEC when considering relative capture width for this area. For Location S3, which is located in the nearshore waters to the southeast of the city of Taizhou and where the water depth is approximately 16 m, Wave Star looks like the best suited WEC when considering capacity factor, with a capacity factor of 18.38%, mean power output of 110.31 kW; Oyster looks like the best suited WEC when considering relative capture width, with relative capture width of 21.58%. But unfortunately, the rather poor capacity factors reported in Table 7 for the different WECs support the non-suitability of this wave climate for present state of the art WECs. In order to achieve cost efficient wave energy conversion, we must develop more suitable WECs for this areas. In addition, the wave farms are mainly constructed in nearshore shallow waters because of the limitation of technology and the cost of energy transmission. When comparing the performances of the state of the art WECs in Locations S1, S2 and S3, the mean power output was higher in Location S3 than in S1 and S2, and the capacity factor and relative capture width were higher too for all WECs. Table 7. Mean power output, capacity factor and relative capture width in each key location.

Discussion for Wave Energy in Key Locations in the Relatively Nearshore Waters
The key locations in the relatively nearshore waters were in areas A5, A6 and A9, which included the relatively nearshore waters to the east of the Zhoushan Islands and to the relatively nearshore waters southeast of the city of Taizhou, Zhejiang Province. The direction of the dominant wave power under operational sea states was mainly in the N and NNE directions, which can serve as a reference to determine the operational direction of certain WECs when acquiring wave energy. Extreme waves mainly originated from the N-SSE directions because of cold swells and storms, which is of interest to developers. The dominant wave conditions all had an H s value of 0.5-4.0 m and a T e value of 4-9 s, and the wave energy accounted for 85% of the total wave energy density; which can be a uniform standard for the design of new WECs that are suitable for these local areas. The wave energy in these key locations was sustainable for each year, and the variation in the wave energy was small, which is advantageous for the development of wave energy projects. The maximum wave power density was 385.30 kW/m, which can be used as a standard for developers to design WECs for survivability. Additionally, the practical use ratios for the wave energy were all greater than 85%. Furthermore, at present, there is no suitable state of the art WECs for cost efficient wave energy conversion in 3 key locations in the relatively nearshore waters of the Zhoushan Islands.

Conclusions
Based on high-resolution and high-accuracy ERA-Interim reanalysis wave field data from 1979 to 2014, the temporal and spatial distributions and the development potential of wave energy were Energies 2017, 10, 1320 23 of 25 studied in detail in offshore and relatively nearshore (key region) waters of the sea areas adjacent to the Zhoushan Islands. The conclusions are as follows: (1) The most suitable wave energy development areas in the offshore waters of the Zhoushan Islands are the eastern relatively nearshore sea areas of the Zhoushan Islands. Wave farms are more difficult to construct in offshore deep water areas and the cost of energy transmission is higher. Therefore, the east relatively nearshore sea areas of the Zhoushan Islands are a better choice for wave energy development. Although wave energy is not abundant in this area compared to global wave energy resources. Nonetheless, the wave energy here is still usable and stable and can serve as a suitable ocean renewable energy for the energy supply of the Zhoushan Islands. The best location for wave energy development is located in the relatively nearshore waters to the southeast of the city of Taizhou (Location S3).