Next Article in Journal
Conversion of Biomass to Chemicals via Electrofermentation of Lactic Acid Bacteria
Previous Article in Journal
Development and Validation of a Nuclear Power Plant Fault Diagnosis System Based on Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nearshore Wave Energy Resource Assessment for Off-Grid Islands: A Case Study in Cuyo Island, Palawan, Philippines

by
Jonathan C. Pacaldo
1,2,3,*,
Princess Hope T. Bilgera
4,5 and
Michael Lochinvar S. Abundo
1,2,5,*
1
Center for Research in Energy Systems and Technology (CREST), University of San Carlos, Cebu City 6000, Philippines
2
Engineering Graduate Program, School of Engineering, University of San Carlos, Cebu City 6000, Philippines
3
Department of Electrical Engineering, Palawan State University, Puerto Princesa City 5300, Philippines
4
Magwayen Enterprise for Management of Environmental Systems Incorporated, Quezon City 1105, Philippines
5
OceanPixel Pte Ltd., 39 Pandan Road, Jurong, Singapore 609281, Singapore
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(22), 8637; https://doi.org/10.3390/en15228637
Submission received: 12 October 2022 / Revised: 10 November 2022 / Accepted: 11 November 2022 / Published: 17 November 2022
(This article belongs to the Section A3: Wind, Wave and Tidal Energy)

Abstract

:
Electrifying off-grid and isolated islands in the Philippines remains one of the challenges that hinders community development, and one of the solutions seen to ensure energy security, energy access and promote a low-carbon future is the use of renewable energy sources. This study determines the nearshore wave energy resource during monsoon seasons in Cuyo Island using a 40-year wave hindcast and 8-year on-site wind speed data as inputs to develop a high-resolution wave energy model using SWAN and assesses its annual energy production through matching with wave energy devices. The results show that the average significant wave height (Hs), peak period (Tp) and wave power density (Pd) during a northeast monsoon are Hs = 1.35 m, Tp = 4.79 s and Pd = 4.05 kW/m, respectively, while a southwest monsoon, which is sheltered by the mainland, results in Hs = 0.52 m, Tp = 3.37 s and Pd = 0.34 kW/m. While the simulated model was observed to overestimate the significant wave height (bias = 0.398, RMSE = 0.54 and SI = 1.34), it has a strong relationship with the “observed values” (average r = 0.9). The annual energy production for Wave Dragon, Archimedes Wave Swing and Seawave Slot-Cone Generator are highest at 1970.6 MWh, 2462.04 MWh, 62.424 MWh and 4099.23 MWh, respectively.

Graphical Abstract

1. Introduction

Resource assessment is an important tool for verifying and quantifying energy resources, and it serves as an initial step in the development of power supply operation. It is also essential in the characterization of an energy resource to support its development. In the Philippines, most off-grid island communities rely heavily on imported oils for their power generation needs [1], and off-grid communities are normally isolated island communities where it is impossible to be connected to the main grid. In 2018, 55.16% of installed capacity was coal and oil based, where coal alone shared 37.14% of these energy needs [2]. Although the Philippines is an archipelagic country, it is up for the task of 100% electrification by 2040 for off-grid areas [3], most of these being isolated small island communities.
Wave energy development can be considered as one of the options in electrifying unviable island communities which cannot easily be reached by government programs because of their geographical constraints [4,5,6,7]. Quantifying wave energy resources in these areas will be the basis of further developing and promoting renewable energy use and will also be the answer to the first three strategic directions of the energy sector in the country, which are to ensure energy security, expand energy access and promote a low-carbon future [3].
Several studies have already been conducted in assessing the Philippines’ wave energy resource. A recent study was conducted by Aminudin, Teh and Pacaldo (2021) [8] in Dumaran Island, Palawan, which assessed the offshore wave energy resources of the island using 40-year hindcast data from MetOceanView (MOV). Quitoras, Abundo and Danao (2018) [9] assessed the energy flux of forty-seven (47) coastal areas in the country, and the results show an energy flux of approximately 10–20 kW/m; this result is within the estimated global wave energy resource assessment as reported in [10,11,12]. Although the study covered a very large area, it does not include Palawan or any part of it or, in particular, the Island of Cuyo. Another study conducted by the Mindanao State University showed that ocean energy in the country can provide an estimated 17,000 megawatts of electricity, and if we can tap this energy, it would be of great help to mitigate the country’s dependency in coal and imported oils as sources of energy [13]. Moreover, a research group from the Marine Science Institute and the College of Engineering of the University of the Philippines started working together for the uptake of ocean renewable energy in the country by identifying potential sites for wave energy resources. Several hot spots had been identified in the northern part of Palawan (the Calamian group of Islands, Dumaran Island, Cuyo Island and El Nido) and in the southern part (Balabac Island and the Cagayan Islands in the Municipality of Cagayancillo) (Figure 1) [14].
However, there are significant knowledge gaps pertaining to quantification of the nearshore wave energy climate and high-resolution wave energy resource model on small islands in semi-enclosed areas that can be used to develop a wave energy project. Having sufficient information regarding wave energy resource potential for this specific type of islands in the Philippines will paved the way for an in-depth development or device solution for small-scale wave energy production to support the island’s power requirements.
Here, a high-resolution nearshore wave model was developed through SWAN (simulating wave nearshore), a third-generation numerical wave model. SWAN was used in different wave resource assessment projects: it was coupled with WAVEWATCH III to determine the wave energy resource along the northern Spanish coast, the model was validated with buoy data to evaluate its accuracy and presented statistical analysis of wave parameters and wave power results [15], the same method was used to determine the nearshore wave energy resource in Canary islands [6] and in Sicily, Italy [16], in Puerto Rico and the United States of Virgin Islands, SWAN was also used to simulate the nearshore wave energy resource for a possible wave power generation in the US Caribbean [17]. Some studies using SWAN for island nearshore wave energy assessment are in the Madiera Islands in Portugal [18], Long Island in New York [19], Hawaiian Islands [20], Azores Islands [21], Cape Verde Islands [22], Sardinia Island [23], Tenerife Island in Spain [4] and in the Baltic Sea [24].
This study used a 40-year wave hindcast from MOV (1978–2018) [25] and 8-year (2010–2018) on-site measurement of wind speed and wind direction at Philippine Atmospheric, Geophysical, and Astronomical Services Administration (PAGASA)—Cuyo Station [26]. SWAN was used to model the wave climate during monsoon seasons and simulate a one (1) year wave hindcast (2018) at six locations (stations A–F, ~12–35 km from the island); for this study, the one year hindcast will be referred to as “hindcast parameters”. The wind data from PAGASA were used as wind input to SWAN to simulate a one (1) year wave hindcast (2018) at the same locations (stations A–F) and were used for validation; for this study, they will be referred to as “observed values”. To determine the annual energy production, a calculation was made by matching a suitable type of wave energy converter (WEC) that will optimize the wave energy resource at selected stations.

2. Materials and Methods

2.1. Study Area

Cuyo Island is the largest island among the 45 islands under Cuyo Archipelago, about 278.37 km northeast of the city of Puerto Princesa, Palawan, and has a land area of 57 km2 (22 mi.2) (Figure 2). It has an estimated population of 34,556 (2015 CENSUS), which is about 4.04% of the total population of Palawan province [27,28]. The island was identified to have a good to excellent wind resource, and the average wind speed measured for 30 years was 5 m/s at 4 m of elevation [29]; given this, a potential wave energy resource as suggested in [14] is highly probable. There are three cases in which waves are propagated: (1), a large storm generates deep water waves that propagate across shallower water while the waves continue to grow due to wind, (2) in the same scenario, a storm generated winds in an area remote from the site of interest and as waves cross shallower water with negligible wind, they propagate to the site as swell, and lastly, (3) wind blows over an area of shallow water, generating waves that grow so large as to interact with the bottom [30]. This indicates that the island’s characteristic satisfies case number three (3), where a good wind resource plays a significant role in wave transformation.

2.2. 40-Year Wave Hindcast Dataset (1978–2018)

To describe the wave climate in Cuyo Island, a wave model was developed using SWAN and will be using MOV’s 40-year, 3-hourly interval wave hindcast dataset (1978–2018, with spatial resolution of 0.5° both for wave and wind) as initial condition and to describe the wave climate surrounding the island. MOV is high-resolution web-based weather forecasting developed by MetOcean Solutions in New Zealand using Ltd. WW3 Tolman Chalicov (MSL WW3 TC) wave model and NOAA Climate Forecast System Reanalysis (CFSR) for the wind model [9]. Several studies using MOV have been published in different fields of study, such as techno-economic assessment of wave energy [9], wave energy resource assessment [30], optimizing hybrid diesel–wave electrical system for an off-grid island [31], weather forecasting for marine operations [32] and marine weather monitoring [33]. For this study, nine (9) stations surrounding Cuyo Island are selected (Figure 3) to describe the offshore wave climate at ~65.7 km from the island by computing its average annual power density (Pd) and significant wave height (Hs), and among these stations, four (4) are selected (stations 4, 8, 12 and 14) that will serve as an initial condition to simulate the numerical wave model during northeast and southwest monsoons and further determine the wave energy resource of stations A–F, which are closer to the island (distance is ~23 km from the island) (Figure 4). The wave climate in these areas (stations A–F) will be used to determine the potential wave energy resource that can be harvested annually by matching it with a suitable nearshore wave energy converter.

2.3. Determination of Hs and Pd

The Hs introduces a well-defined and standardized statistic to denote the characteristic height of the random waves in a sea state [34]. In a time-domain analysis, Hs is defined as the mean height of the highest one-third of all waves, and H1/3 is denoted in
H 1 3 = 1 N 3   i = 1 N 3 H i
where N is the number of individual wave heights, and Hi is a series of wave heights ranked from highest to lowest [35].
To estimate the wave power, the most commonly used equation regardless of water depth is given by
P = 0.577 H s 2 T
where Hs is the wave height, and T is the given period [9].

2.4. Validation

To validate the accuracy of the resulting wave model, statistical analysis between the model result and the observed values will be calculated using the following statistical metrics or error statistics [15]:
Mean   of   measured   parameters ;   X ¯ = 1 n   x i , ,
Mean   of   hindcast   parameters ;   Y ¯ = 1 n   y i ,
Bias ;   b = 1 n   ( y i x i ) ,
Root   Mean   Square   Error ; = 1 n ( x i y i ) 2 ,
Scatter   Index ;   SI = R M S E X ¯ ,
Pearson s   Correlation   Coefficient ;   r = ( x i X ¯ ) ( y i Y ¯ ) ( x i X ¯ ) 2 ( y i Y ¯ ) 2 ,
Here, x i   is the Hs of the “observed values”, y i is the Hs of the “hindcast parameters” and n is total values for both parameters.

2.5. SWAN Wave Model

In this study, a nested model was used to provide the necessary boundary conditions for the Cuyo wave model. The coarse grid is a rectangular 110 × 120 grid with ~1.5 km resolution, rotated 45° to align the grid with the dominant wind and wave directions due to the northeast (NE, Amihan) and southwest (SW, Habagat) monsoons (Figure 5, white grid). On the other hand, the nested, high-resolution grid is a 123 × 93 rectangular grid with ~500 m resolution and focused on the east side of the Cuyo Archipelago (Figure 5, blue grid). Subsequently, downloaded bathymetric data from GEBCO2021 [36] with 450 m resolution were also interpolated onto the model grid (Figure 6).
Four simulations were performed to represent the NE and SW monsoon wave conditions. The boundary and wave conditions (Hs, and peak wave direction (Dp) and period (Tp) were assumed to be uniform along the specific boundary orientation) inputted into the model runs were based on the computed mean of 3-hourly MOV 2008–2018 data (Table 1) taken from the following stations (Figure 7):
  • Station 4 (southwest boundary orientation) at 10.5° N and 120.5° E.
  • Station 8 (southeast boundary orientation) at 10.5° N and 121.5° E.
  • Station 12 (northeast boundary orientation) at 11.5° N and 121.5° E.
  • Station 14 (northwest boundary orientation) at 11.5° N and 120.5° E.
Additionally, the following parameter settings were applied in the wave model:
  • Wave spectrum.
    At the wave model boundary, a JONSWAP spectrum with a peak enhancement factor of 3.3 was assumed.
    Similarly, a directional spreading of approximately 25° (power function, with power = 4) was assumed.
  • Physical parameters.
    Third-generation mode for wind growth, quadruplet interactions and whitecapping (based on Komen, Hasselmann and Hasselmann, 1984 [37]) were considered.
    Constant depth induced breaking (alpha = 1 (the coefficient for determining the rate of dissipation), and gamma = 0.73 (the value of the breaker parameter defined as Hm0/d)).
    Constant JONSWAP bottom friction (friction coefficient = 0.067 m2/s3). (The bottom friction is computed based on the empirical model of JONSWAP (Hasselmann et al., 1973) [38]. The coefficient of the JONSWAP formulation is set at 0.067 m2/s3, which is a typical default value for wind sea conditions.)
    Non-linear wave–wave interactions due to the triads were not considered.
    No diffraction.
  • Numerical parameters.
    The amount of diffusion of the implicit scheme in the directional space (through directional discretization parameter) and frequency space (through the frequency discretization) was set to the default value 0.5.
    Accuracy:
    Relative change HsTm01: 0.02.
    Relative change with respect to the mean value: 0.02 for both Hs and Tm01.
    Convergence percentage of wet grid points: 98%.
    Maximum number of iterations: 15.

2.6. Annual Energy Production (AEP)

AEP was computed using the stations’ wave scatter diagram and wave energy devices power matrix using the formula given in Equation (9) [9,39]. Wave scatter diagram is the condition of sea state at a particular location in a year, which is generated from the historical data (see Appendix A, Figure A1), while the wave energy converters’ (WEC) power matrix is the actual amount of available energy the device can capture (see Appendix B, Figure A2). This study uses Wave Dragon 5900 kW (WD5.9), Wave Dragon 7000 kW (WD7), Archimedes Wave Swing (AWS) (2470 kW) and Seawave Slot-Cone Generator (SSG) (20,000 kW) to determine the potential annual energy that can be utilized surrounding Cuyo Island. Wave Dragon is classified as a floating overtopping wave energy device that can be considered both as a nearshore (working at depth higher than 6 m) or as an offshore device (depth higher than 25 m) [39]. AWS is a heaving point absorber classified as a nearshore device; however, unlike other devices, it does not have a fixed-rated power level, but the power output continues to rise along with both wave height and period [40]. SSG is an overtopping wave energy which works under a wide spectrum of different wave conditions, giving a high overall efficiency, and it is considered a shoreline type WEC [41]. This device is considered in this study because the shoreline on the southern part of the island probably has the same intensity as stations D and E and can be considered to a much smaller island community near the stations A–F. These devices are all considered applicable to remote islands [42,43], which makes them suitable in the study area for evaluating the annual energy that can potentially be harvested.
AEP in MWh = Wave Scatter Diagram (in hours) × Power matrix (MW)

3. Results

3.1. Dominant Wind Direction (PAGASA—Cuyo Station) vs. Dominant Wave Direction (Station 10)

This study focused on the nearshore wave energy resource of the island, and the nearest dataset that will help describe the wave behavior nearshore is the 8-year (2010–2017) PAGASA—Cuyo Station wind speed and wind direction dataset and the 40-year hindcast wave height and wave direction at station 10 (see Figure 3 for the location of station 10), which is 15 km away from the island and has an annual significant wave height (Hs) of 1.2 m and annual wave power density (Pd) of 3.13 kW/m [31]. The wind and wave directions are presented using wind rose and wave rose diagram; a rose diagram represents the two-dimensional orientation of the wind and wave climate that represents the relative frequencies of different wind and wave directions and the wind and wave heights over a period of time. It displays the distribution of data in a way that can be easily understood and evaluated [44]. Figure 8 shows the wave rose diagram of the 40-year hindcast at station 10, taken every five (5) years starting from 1978 to 2018. It can be observed that the dominant wave directions are consistent and are coming from the northeastern and southwestern sides of the island. This is mainly due to the northeast monsoon, which typically occurs from the months of December–February, and the southwest monsoon, which occurs from the months of June–August. The dominant wave direction is consistent with the eight (8)-year on-site wind measurement from 2010–2018 (Figure 9). It can be seen that the dominant wind directions are also coming from the northeastern and southwestern parst of the island. These observations strengthen the use of hindcast data as initial conditions to simulate the nearshore wave climate of the island. Additionally, a high correlation (r = 0.75) between the hindcast wave data in station 10 and PAGASA on-site wind measurements was computed and presented in Table 2 [31], and the table also shows the distance of the nine (9) stations from PAGASA station, as well as the average Hs and Pd at these stations. Similarly, the coefficient of correlation® presented in the analysis coincides with other studies, which validates the quality of the MOV hindcast data (Table 3).

3.2. Wind and Wave Climatology from MOV and PAGASA—Cuyo Station Data

Using the MOV Station 10 data for 2008–2018, the monthly wave climatology in Cuyo was generated (Figure 10). The monthly variability of the significant wave heights (median) ranges from 0.1–1.1 m with extreme significant wave heights reaching almost 2.7 m (January) and outliers reaching as high as 4.1 m. The outliers are usually found all throughout the year, with the highest outliers in November and December. The seasonal signal of the Hs related to the monsoons is also observed, wherein higher wave heights are recorded during the monsoon peaks and lower wave heights during the monsoon transition period. The months with the highest extreme Hs and outliers coincide with the northeast monsoon months (December, January and February).
Likewise, the compass roses of Hs (Figure 11) and Tp (Figure 12) with propagation direction show the strong monsoonal influence with most waves clustered along the N–NNE and WSW–SW–SSW directions. It is also notable in Figure 11 that Hs is higher during the NE monsoon (~25% Hs is ≥1.5 m) compared to the SW monsoon (~5% Hs is ≥1.5 m). This observation is also consistent with Tp (Figure 12), wherein ~80% (majority of the waves) is ≥4 s during the NE monsoon, while only ~40% is ≥4 s during the SW monsoon.
In terms of wind velocities, the 2010–2018 data from MOV station 10 (3 hourly data) and PAGASA wind station at Cuyo (daily data) were used to determine the wind conditions in the area. It is to be noted that the MOV station is further offshore and located ~16 km north of the PAGASA—Cuyo Station (more sheltered since it is located on land and only at 4 m elevation), thus showing some differences in wind velocities (Figure 13). Noticeably, wind speeds from PAGASA station (majority below 8 m/s) are weaker compared to MOV data (~30% winds are ≥8 m/s) (Figure 13). In terms of direction, the northeasterly wind produces stronger winds (~25% winds are ≥8 m/s) compared to the southwesterly wind (~5% winds are ≥8 m/s) (Figure 13). Moreover, the prominent wind direction, especially from the PAGASA station (SSW and NE wind), agrees with the wave direction.

3.3. Monsoon Wave Model for Cuyo Island

MOV’s wave data and PAGASA’s wind data both agree that the observed Hs is related to the monsoons, wherein higher wave heights are recorded, and lower wave heights take place during monsoon transition periods. The wave model developed describes the wave climate surrounding the island and specifically determines the wave parameters in the six (6) points of interest (stations A–F), which are closer to the island (Figure 14). Stations A–F have an average distance of 24.3 km from the island: the closest is 12.35 km (station A), and the farthest is 35.35 km (station C). During the northeast monsoon season, Hs and Tp are the highest at station F (1.49 m and 4.87 s, respectively), followed by stations A and E, (1.43 m and 4.7 s, respectively); other stations are much lower due to the sheltering effects of nearby islets and the main island (Table 4) (Figure 15a,b). Results of these points are expected to be a little lower than the MOV stations (Table 5) because these points are shallower [49] and more sheltered from the northeast monsoon and southwest monsoon winds [50].
Figure 15, Figure 16 and Figure 17 show the wave model results for Hs, Tp and energy transfer, respectively, during the northeast monsoon, and it can be observed that in the northeastern and southeastern sides of the island, specifically on stations A, E and F, that the results are the highest, having an average Hs, Tp and Pd of 1.35 m, 4.79 s and 4.05 kW/m, respectively, the highest being at station F with Pd = 4.25 kW/m and Tp = 4.87 s. These results can be confirmed in Figure 18, which shows that the peak wind directions converged on these areas due to less obstructions from other islets, resulting in a minimal sheltering effect and its exposure to the much stronger northeast monsoon winds. Moreover, the water depth at those stations is much deeper compared to other stations, as seen in the bathymetry data, giving more active wave behavior.
Figure 19, Figure 20 and Figure 21 show the wave model results for the Hs, Tp and Pd model for southwest monsoon season, respectively, and stations A–F have an average Hs = 0.52 m, Tp = 3.37 s and Pd = 0.34 kW/m. The model results from the southwest monsoon model have lower values compared to the northeast monsoon model due to the sheltering effect of the mainland of Palawan, which makes the wind passing from the southwest decrease its magnitude before reaching Cuyo Island.

3.4. Model Validation

Table 6 shows the statistical metrics result between the “observed values” and “hindcast parameters”, and it is shown that the “hindcast parameters” were overestimated by the wave model compared to the “observed values”, having an average bias of 0.398 m, where the lowest is at station F at 0.38 m and the highest at station E at 0.42 m. A good model has an RMSE value close to zero, and the closer it is the better ability of the model to accurately predict the data. Here, the RMSE has an average value of 0.54, the lowest at 0.51 (station F) and the highest at 0.56 (station E), which in this case is an acceptable value considering the complex topography on land where the “observed values” were taken. In general, overestimation and the degree of variance on all stations may be explained by the following factors: the distance, location and topographical differences of the two parameters. PAGASA—Cuyo Station is located on land and has an average distance from stations A to F of about 24.3 km, the nearest at 12.35 km (station A) and the farthest at 35.35 km (station C), and the hub height of the equipment used to measure wind speed and wind direction in PAGASA—Cuyo Station is only at 4 m of elevation, where it is more prone to disturbances caused by nearby buildings and trees. Moreover, the locations of stations A, B and C are in between small islands, which may affect or distort the wind and wave flow on the area, whereas stations D, E and F are more exposed with less sheltering from nearby islands (see Figure 14).
Scatter index (SI) values also tell us a constant overestimation of the model, with all the results being higher than one, the lowest at station F being 1.08 and the highest being 1.58 (station C). These values may also be due to the three factors mentioned, which tend to overestimate the simulated wave height. Figure 22 shows the scatter plot graph of the two parameters; in general, the data are clustered linearly and tend to increase exponentially at higher values starting between 0.8 and 1.0 m of the observed values. Few outliers can be observed on all stations but may not affect the overall result of the wave parameters. Finally, the data show a strong positive relationship with an average r = 0.90, and the highest is at station F (r = 0.94) followed by station A (r = 0.92), and the lowest being station C (r = 0.85).

3.5. AEP

Table 7 shows the AEP of the three WECs computed using Equation (9). Stations E and F draw the highest energy production for all the WECs tested. At station E, WD5.9, WD7 and AWS draw the highest energy production with 1970.61 MWh, 2462.04 MWh and 62.42 MWh, respectively, while at station F, SSG has the highest AEP with 4099.23 MWh. The capacity factors for three WECs are <3.5%, the highest is WD7 at 3.4%, and the lowest is AWS at 0.2%. This is because most of the data are at lower values and are not within the devices’ energy production capability. Cuyo Island has an average power demand of 0.984 MW and a peak power demand of 1.4 MW [31]. Although the highest capacity factor is only at 3.4%, the AEP of WD7 can supply ~28.53% of the islands average power demand and ~20.0% of the islands peak power demand. Given this scenario and the resulting wave model, when applied to a much smaller island community in the Cuyo Archipelago, a single WD7 can perhaps cover their energy needs. A sample computation of the AEP is presented in Table 8 and the rest are in Appendix C (Figure A3, Figure A4, Figure A5 and Figure A6).

4. Conclusions

Wave energy resource assessment in the Philippines ranges from 10–20 kW/m as reported in Quitoras et al. [9], in the 47 sites under the coastal regions of Catanduanes, Samar, Siargao Island, Surigao del Sur and Western Luzon. In Wan et. al, [51] along Luzon strait, the exploitable wave energy resource is at 10–15 kW/m, which agrees with the findings in [9]. In Dumaran’s study [8], the wave energy resource surrounding the island within a 100 km radius is less than 4.5 kW/m, and this is in agreement with Mirzae et al. [52], where the semi-enclosed sea or sheltered areas have a lower probability of harnessing wave energy resource that will exceed 5 kW/m at any season. With the same topographical characteristic as Dumaran Island, the resulting nearshore Pd in Cuyo Island during the monsoon seasons is also less than 5 kW/m, the highest being 4.25 kW/m (station F) and the lowest being 2.36 kW/m (station C). During the northeast monsoon (Amihan Season) and during the southeast monsoon (Habagat Season), the highest and lowest Pd are at 0.42 kW/m and 0.19 kW/m, respectively, which is in agreement with [8] and [48]. Although the “hindcast parameters” overestimate the “observed values” with average bias, RSME and SI values of 0.398, 0.54 and 1.34, respectively, the “hindcast parameters” have a strong positive relationship with the “observed values” having an average correlation coefficient (r) of 0.90, the highest being 0.94 (station F) and the lowest being 0.85 (station C). This signifies that the measured wind data are in agreement with the simulated wave data and therefore can be used as a reference in analyzing the nearshore wave energy resource surrounding the island of Cuyo either for exploitation or testing of nearshore and shoreline WEC [53,54]. The highest AEP is at station F, SSG with 4099.23 MWh, followed by WD7 with 2462.04 MWh, but in terms of capacity factor, WD7 is the highest at 3.4%; low capacity factors happen because the majority of the data falls on lower values of Hs and Tp where the wave energy device is not capable of producing energy or the energy produced is low. Although the highest is only at a 3.4% capacity factor, it can supply 28.53% of the islands’ average power demand or 20.0% of their peak demand. If these scenarios can be realized, they will be able to lessen the country’s dependency on fossil fuels for electrifying off-grid and isolated islands. Results show that downscaling of WEC to increase the capacity factor of the device in milder resources can answer the challenges in bringing sustainable renewable energy resources to unviable islands in the Philippines or in any areas of the same climate characteristics.
To predict more accurately and enhance the high-resolution wave model of the wave energy resource in an isolated island of the same characteristics, it is recommended to have an on-site measurement of the wave parameters for a minimum of 1 year or maybe longer.

Author Contributions

Conceptualization, J.C.P. and M.L.S.A.; methodology, J.C.P. and M.L.S.A.; software, J.C.P. and P.H.T.B.; validation, J.C.P. and P.H.T.B.; formal analysis, J.C.P., P.H.T.B. and M.L.S.A.; investigation, J.C.P.; resources, J.C.P.; data curation, J.C.P.; writing—original draft preparation, J.C.P.; writing—review and editing, J.C.P., P.H.T.B. and M.L.S.A.; visualization, J.C.P. and M.L.S.A.; supervision, revision and verifying the results, J.C.P., P.H.T.B. and M.L.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

The scientific work was supported by the Department of Science and Technology—Engineering Research and Development for Technology (DOST–ERDT), Republic of the Philippines.

Data Availability Statement

The authors and MetService make no representation or warranties regarding the accuracy or completeness of the data, the use to which data may be put or the results which may be obtained from using the data, and the authors and MetService accept no liability for any loss or damage (whether direct or indirect) incurred by any person through the use of or reliance on the data.

Acknowledgments

The authors would like to thank the University of San Carlos DOST—Engineering Research and Development for Technology (ERDT) Scholarship, Center for Research and Energy Systems Technology (CREST), and Palawan State University for their support and motivation in conducting this research and also to MetOcean Solutions, for the 40-year (1978–2018) wave hindcast dataset used (produced and distributed by MetOcean Solutions, New Zealand (www.metocean.co.nz) accessed on 7 October 2019). Thanks also for the support from the following individuals: to Shalou Maratas of PAGASA—Climate and Agrometeorological Data section (CADS) for the 8-year on-site wind parameters data; and to Weslie Capute of Oceanpixel, Pte. Ltd., Singapore, for supporting this research work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Wave scatter diagram of stations A–F.
Figure A1. Wave scatter diagram of stations A–F.
Energies 15 08637 g0a1

Appendix B

Figure A2. Power matrix of WD7, WD5.9, AWS and SSG.
Figure A2. Power matrix of WD7, WD5.9, AWS and SSG.
Energies 15 08637 g0a2

Appendix C

Figure A3. AEP computations of WD7 for stations A–F.
Figure A3. AEP computations of WD7 for stations A–F.
Energies 15 08637 g0a3
Figure A4. AEP computations of WD5.9 at stations A–F.
Figure A4. AEP computations of WD5.9 at stations A–F.
Energies 15 08637 g0a4
Figure A5. AEP computations of AWS at stations A–F.
Figure A5. AEP computations of AWS at stations A–F.
Energies 15 08637 g0a5
Figure A6. AEP computations of SSG at stations A–F.
Figure A6. AEP computations of SSG at stations A–F.
Energies 15 08637 g0a6

References

  1. IRENA. Renewables Readiness Assessment; International Renewable Energy Agency: Abu Dhabi, United Arab Emirates, 2017; Available online: https://www.irena.org/publications/2017/Mar/Renewables-Readiness-Assessment-The-Philippines (accessed on 20 May 2022).
  2. Republic of the Philippines, Department of Energy, EPIMB. 2018 Power Supply and Demand Highlights. Available online: http://161.49.106.166/electric-power/2018-power-supply-and-demand-highlights (accessed on 20 May 2022).
  3. Department of Energy. Philippine Energy Plan 2017–2040: Sectoral Plans and Roadmaps. Available online: https://www.doe.gov.ph/pep/renewable-energy-roadmap-2017-2040?withshield=1 (accessed on 20 May 2022).
  4. Veigas, M.; Ramos, V.; Iglesias, G. A wave farm for an island: Detailed effects on the nearshore wave climate. Energy 2014, 69, 801–812. [Google Scholar] [CrossRef]
  5. Rusu, E.; Onea, F. An assessment of the wind and wave power potential in the island environment. Energy 2019, 175, 830–846. [Google Scholar] [CrossRef]
  6. Gonçalves, M.; Martinho, P.; Soares, C.G. Assessment of wave energy in the Canary Islands. Renew. Energy 2014, 68, 774–784. [Google Scholar] [CrossRef]
  7. Moschos, E.; Manou, G.; Dimitriadis, P.; Afentoulis, V.; Koutsoyiannis, D.; Tsoukala, V.K. Harnessing wind and wave resources for a Hybrid Renewable Energy System in remote islands: A combined stochastic and deterministic approach. Energy Procedia 2017, 125, 415–424. [Google Scholar] [CrossRef]
  8. Aminudin, A.; Teh, H.M.; Pacaldo, J. Wave Energy Assessment in Dumaran Island, Palawan, Philippines. Int. J. Coast. Offshore Eng. 2021, 6, 51–63. Available online: https://www.ijcoe.org/article_152614_2defe337bb73a06acf5d3e6a8b118777.pdf (accessed on 5 December 2021).
  9. Quitoras, M.R.; Abundo, M.L.S.; Danao, L.A. A techno-economic assessment of wave energy resources in the Philippines. Renew. Sustain. Energy Rev. 2018, 88, 68–81. [Google Scholar] [CrossRef]
  10. Asian Development Bank. Wave Energy Conversion and Ocean Thermal Energy Conversion Potential in Developing Member Countries. Available online: https://www.adb.org/publications/wave-energy-conversion-and-ocean-thermal-energy-conversion-potential-developing-member (accessed on 20 May 2022).
  11. Cornett, A.M. A global wave energy resource assessment. In Proceedings of the International Offshore and Polar Engineering Conference, Vancouver, BC, Canada, 6–11 July 2008; pp. 1–9. [Google Scholar]
  12. Mork, G.; Barstow, S.; Kabuth, A.; Meresa, P. Assessing the global wave energy potential. In Proceedings of the 29th International Conference on Ocean, Offshore Mechanics and Artic Engineering, Shanghai, China, 6–11 June 2010. [Google Scholar]
  13. Philstar Global. Available online: https://www.philstar.com/business/science-and-environment/2016/08/18/1614467/marine-scientist-eyes-rd-ocean-energy (accessed on 20 May 2022).
  14. Quirapas, M.A.J.R.; Lin, H.; Abundo, M.L.S.; Brahim, S.; Santod, D. Ocean renewable energy in Southeast Asia: A review. Renew. Sustain. Energy Dev. 2015, 41, 799–817. [Google Scholar] [CrossRef]
  15. Bryant, M.A.; Hesser, T.J.; Jensen, R.E. Evaluation Statistics Computed for the Wave Information Studies (WIS). Available online: https://apps.dtic.mil/sti/pdfs/AD1013235.pdf (accessed on 23 May 2022).
  16. Akpınar, A.; Kömürcü, M.İ. Assessment of wave energy resource of the Black Sea based on 15-year numerical hindcast data. Appl. Energy 2013, 101, 502–512. [Google Scholar] [CrossRef]
  17. Booij, N.; Ris, R.; Holthuijsen, L. A Third-Generation Wave Model for Coastal Regions. I—Model Description and Validation. J. Geophys. Res. 1999, 104, 7649–7666. [Google Scholar] [CrossRef] [Green Version]
  18. Bento, A.R.; Martinho, P.; Soares, C.G. Wave energy assessment for Northern Spain from a 33-year hindcast. Renew. Energy 2018, 127, 322–333. [Google Scholar] [CrossRef]
  19. Silander, M.F.C.; Moreno, C.G.G. On the spatial distribution of the wave energy resource in Puerto Rico and the United States Virgin Islands. Renew. Energy 2019, 136, 442–451. [Google Scholar] [CrossRef]
  20. Rusu, E.; Soares, C.G. Wave energy pattern around the Madeira Islands. Energy 2012, 45, 771–785. [Google Scholar] [CrossRef]
  21. Buonaiuto, F.S., Jr.; Slattery, M.; Bokuniewicz, H.J. Wave modeling of Long Island coastal waters. J. Coast. Res. 2011, 27, 470–477. [Google Scholar] [CrossRef]
  22. Stopa, J.E.; Filipot, J.F.; Li, N.; Cheung, K.F.; Chen, Y.L.; Vega, L. Wave energy resources along the Hawaiian Island chain. Renew. Energy 2013, 55, 305–321. [Google Scholar] [CrossRef]
  23. Rusu, L.; Soares, C.G. Wave energy assessments in the Azores islands. Renew. Energy 2012, 45, 183–196. [Google Scholar] [CrossRef]
  24. Myslenkov, S.A.; Medvedeva, A.Y. Wave energy resources of the Baltic Sea and coastal zone of the Kaliningrad Region. Фундаментальная И Прикладная Гидрoфизика 2019, 12, 34–42. [Google Scholar] [CrossRef]
  25. MetOceanView. Available online: https://www.metoceanview.com/hindcast (accessed on 7 October 2019).
  26. PAGASA National Meteorological, Hydrological and Astronomical Service Agency of the Philippines. Available online: https://bagong.pagasa.dost.gov.ph/index.php (accessed on 14 March 2020).
  27. Cuyo, Palawan Profile—PhilAtlas. Available online: https://www.philatlas.com/luzon/mimaropa/palawan/cuyo.html (accessed on 18 April 2020).
  28. Magsaysay, Palawan Profile—PhilAtlas. Available online: https://www.philatlas.com/luzon/mimaropa/palawan/magsaysay.html (accessed on 18 April 2020).
  29. Blechenger, P.; Ermino, R.; Lopez, A.; Serafica, E.; Terrado, E.; Solar Hybridization of Large-Scale NPC-SPUG Diesel Power Plants: A Programmatic Approach. European Union–Philippines Access to Sustainable Energy Program. 22 June 2018. Available online: https://www.eu-asep.ph/wp-content/uploads/2018/06/Solar-hybridization-of-large-scale-NPC-SPUG-diesel-power-plants-Programmatic-Approach-by-Philipp-Blechinger-PhD.pdf (accessed on 20 May 2022).
  30. USACE Coastal Engineering Manual. 2002. Available online: https://www.publications.usace.army.mil/USACE-Publications/Engineer-Manuals/u43544q/636F617374616C20656E67696E656572696E67206D616E75616C/ (accessed on 18 April 2020).
  31. Pacaldo, J.C.; Abundo, M.L.S.; Billotindos, L.M.; Baco, C.S., Jr. Performance Analysis of a Hybrid Diesel–Renewable Energy (RE) Electrical System in Cuyo Island, Palawan, Philippines. Int. J. Adv. Res. Eng. Innov. 2022, 4, 1–15. [Google Scholar] [CrossRef]
  32. McComb, P.J.; Johnson, D.L. A high-Resolution Weather Forecasting Tool for Marine Operations Management in Ports and Harbours. Available online: https://www.readkong.com/page/a-high-resolution-weather-forecasting-tool-for-marine-6291286 (accessed on 23 May 2022).
  33. Boulay, S.O.; Batt, L. Marine weather forecasting and monitoring at the Port of Sydney and Botany Bay, NSW, Australia. In Proceedings of the Australasian Coasts & Ports Conference, Auckland, New Zealand, 15–18 September 2015; p. 102. Available online: https://www.researchgate.net/profile/Sebastien-Boulay/publication/283824032_Marine_weather_forecasting_and_monitoring_at_the_Port_of_Sydney_and_Botany_Bay_NSW_Australia/links/564812b208ae9f9c13e97f3b/Marine-weather-forecasting-and-monitoring-at-the-Port-of-Sydney-and-Botany-Bay-NSW-Australia.pdf (accessed on 5 May 2020).
  34. Bai, Y.; Jin, W.L. Chapter 5-Wave Loads for Ship Design and Classification. In Marine Structural Design, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2016; pp. 73–93. [Google Scholar] [CrossRef]
  35. Bai, Y. Marine Structural Design; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar] [CrossRef]
  36. GEBCO. 2021. Available online: https://www.gebco.net/data_and_products/gridded_bathymetry_data/gebco_2021.html (accessed on 6 February 2022).
  37. Komen, G.; Hasselmann, S.; Hasselmann, K. On the existence of a fully developed wind–sea spectrum. J. Phys. Oceanogr. 1984, 14, 1271–1285. [Google Scholar] [CrossRef]
  38. Hasselmann, K.; Barnett, T.; Bouws, E.; Carlson, H.; Cartwright, D.; Enke, K.; Walden, H. Measerements of wind wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Dtsch. Hydrogr. Z. 1973, 8, 12. [Google Scholar]
  39. Nezhad, M.M.; Groppi, D.; Rosa, F.; Piras, G.; Cumo, F.; Garcia, D.A. Nearshore wave energy converters comparison and Mediterranean small island grid integration. Sustain. Energy Technol. Assess. 2018, 30, 68–76. [Google Scholar] [CrossRef]
  40. Bozzi, S.; Besio, G.; Passoni, G. Wave power technologies for the Mediterranean offshore: Scaling and performance analysis. Coast. Eng. 2018, 136, 130–146. [Google Scholar] [CrossRef]
  41. Vicinanza, D.; Margheritini, L.; Kofoed, J.P.; Buccino, M. The SSG wave energy converter: Performance, status and recent developments. Energies 2012, 5, 193–226. [Google Scholar] [CrossRef] [Green Version]
  42. Fadaeenejad, M.; Shamsipour, R.; Rokni, S.D.; Gomes, C. New approaches in harnessing wave energy: With special attention to small islands. Renew. Sustain. Energy Rev. 2014, 29, 345–354. [Google Scholar] [CrossRef]
  43. Silva, D.; Rusu, E.; Soares, C.G. Evaluation of various technologies for wave energy conversion in the Portuguese nearshore. Energies 2013, 6, 1344–1364. [Google Scholar] [CrossRef] [Green Version]
  44. Sanderson, D.J.; Peacock, D.C. Making rose diagrams fit-for-purpose. Earth-Sci. Rev. 2020, 201, 103055. [Google Scholar] [CrossRef]
  45. Rapizo, H.; Godoi, V.A.; Perez, J.; Durrant, T.; Guedes, R. Analysis of the New Zealand’s Taranaki regional wave climate using high-resolution modelling. Reg. Stud. Mar. Sci. 2021, 45, 101806. [Google Scholar] [CrossRef]
  46. Belibassakis, K.; Bleuanus, S.; Vermeiden, J.; Townsend, N. Combined performance of innovative biomimetic ship propulsion system in waves with Dual Fuel ship engine and application to short-sea shipping. In Proceedings of the 31st International Ocean and Polar Engineering Conference, OnePetro, Rhodes, Greece, 20–25 June 2021. [Google Scholar]
  47. Echevarria, E.R.; Hemer, M.A.; Holbrook, N.J. Seasonal variability of the global spectral wind wave climate. J. Geophys. Res. Ocean. 2019, 124, 2924–2939. [Google Scholar] [CrossRef]
  48. Rickard, G. Review of the Report. Astrolabe Reef Metocean Conditions: Wave, Ocean Current and Wind Statistics; NIWA Taihoro Nukurangi: Auckland, New Zealand, 2014. [Google Scholar]
  49. Johnson, H.K.; Kofoed-Hansen, H. Influence of bottom friction on sea surface roughness and its impact on shallow water wind wave modeling. J. Phys. Oceanogr. 2000, 30, 1743–1756. [Google Scholar] [CrossRef]
  50. De León, S.P.; Soares, C.G. On the sheltering effect of islands in ocean wave models. J. Geophys. Res. Ocean. 2005, 110. [Google Scholar] [CrossRef]
  51. Wan, Y.; Zhang, J.; Meng, J.; Wang, J. Exploitable wave energy assessment based on ERA-Interim reanalysis data—A case study in the East China Sea and the South China Sea. Acta Oceanol. Sin. 2015, 34, 143–155. [Google Scholar] [CrossRef]
  52. Mirzaei, A.; Tangang, F.; Juneng, L. Wave energy potential assessment in the central and southern regions of the South China Sea. Renew. Energy 2015, 80, 454–470. [Google Scholar] [CrossRef]
  53. Xia, G.; Draxl, C.; Raghavendra, A.; Lundquist, J.K. Validating simulated mountain wave impacts on hub-height wind speed using SoDAR observations. Renew. Energy 2021, 163, 2220–2230. [Google Scholar] [CrossRef]
  54. Wang, Z.; Dong, S.; Dong, X.; Zhang, X. Assessment of wind energy and wave energy resources in Weifang sea area. Int. J. Hydrog. Energy 2016, 41, 15805–15811. [Google Scholar] [CrossRef]
Figure 1. Probable sites for wave energy development (yellow boxes) in the province of Palawan identified by University of the Philippines’ Marine Science Institute.
Figure 1. Probable sites for wave energy development (yellow boxes) in the province of Palawan identified by University of the Philippines’ Marine Science Institute.
Energies 15 08637 g001
Figure 2. Map showing the study site, Cuyo Island at the northeastern part of Palawan.
Figure 2. Map showing the study site, Cuyo Island at the northeastern part of Palawan.
Energies 15 08637 g002
Figure 3. Cuyo Island showing the 9 stations of MOV’s 40-year hindcast wave data used to determine the average Pd and Hs surrounding Cuyo Archipelago (yellow dots) and used as initial conditions to SWAN (stations 4, 8, 12 and 14).
Figure 3. Cuyo Island showing the 9 stations of MOV’s 40-year hindcast wave data used to determine the average Pd and Hs surrounding Cuyo Archipelago (yellow dots) and used as initial conditions to SWAN (stations 4, 8, 12 and 14).
Energies 15 08637 g003
Figure 4. Chosen sites of interest (stations A–F) about 23 km from PAGASA station used to match with WEC to determine the potential annual energy production surrounding the island.
Figure 4. Chosen sites of interest (stations A–F) about 23 km from PAGASA station used to match with WEC to determine the potential annual energy production surrounding the island.
Energies 15 08637 g004
Figure 5. Nested grids (white for coarse grid, 110 × 120 with ~150 km resolution; blue for fine grid, 123 × 93 with ~500 m resolution) used for the SWAN modelling.
Figure 5. Nested grids (white for coarse grid, 110 × 120 with ~150 km resolution; blue for fine grid, 123 × 93 with ~500 m resolution) used for the SWAN modelling.
Energies 15 08637 g005
Figure 6. Downloaded GEBCO bathymetric data (adjusted to positive values in meters for Delft3d compatibility) interpolated onto the coarse (left panel) and fine (right panel) grid of the SWAN model. Red is deepest and blue is shallowest. Contour intervals for coarse grid (left panel) are every 100 m depth, while it is every 10 m of depth for fine grid (right panel).
Figure 6. Downloaded GEBCO bathymetric data (adjusted to positive values in meters for Delft3d compatibility) interpolated onto the coarse (left panel) and fine (right panel) grid of the SWAN model. Red is deepest and blue is shallowest. Contour intervals for coarse grid (left panel) are every 100 m depth, while it is every 10 m of depth for fine grid (right panel).
Energies 15 08637 g006
Figure 7. Map of MOV stations (red squares) where 3-hourly data on wind velocity, Hs, Tp and Dp were extracted. On the other hand, the yellow circle is the PAGASA Cuyo station where daily wind data were recorded.
Figure 7. Map of MOV stations (red squares) where 3-hourly data on wind velocity, Hs, Tp and Dp were extracted. On the other hand, the yellow circle is the PAGASA Cuyo station where daily wind data were recorded.
Energies 15 08637 g007
Figure 8. Wave rose diagram at station 10 (11.0° N Lat., 121.0° E Long.) having a 5-year interval (data source: MOV).
Figure 8. Wave rose diagram at station 10 (11.0° N Lat., 121.0° E Long.) having a 5-year interval (data source: MOV).
Energies 15 08637 g008
Figure 9. Wind Rose diagram of PAGASA–Cuyo Station (10.85° N Lat., 121.04° E Long.) from 2010–2017 (data source: PAGASA–Cuyo Station).
Figure 9. Wind Rose diagram of PAGASA–Cuyo Station (10.85° N Lat., 121.04° E Long.) from 2010–2017 (data source: PAGASA–Cuyo Station).
Energies 15 08637 g009
Figure 10. Box-and-whisker plot (upper panel) and time-series plot (lower panel) of Hs from MOV. Red line inside the box of the box whisker plot shows the monthly median, box edges represent the 25th and 75th percentiles, the whiskers extend to the most extreme data points not considered outliers, while the red “+” marker symbol represents the outliers.
Figure 10. Box-and-whisker plot (upper panel) and time-series plot (lower panel) of Hs from MOV. Red line inside the box of the box whisker plot shows the monthly median, box edges represent the 25th and 75th percentiles, the whiskers extend to the most extreme data points not considered outliers, while the red “+” marker symbol represents the outliers.
Energies 15 08637 g010
Figure 11. Compass rose of Hs and Dp during the northeast and southwest monsoons.
Figure 11. Compass rose of Hs and Dp during the northeast and southwest monsoons.
Energies 15 08637 g011
Figure 12. Compass rose of Tp and Dp during the northeast and southwest monsoon.
Figure 12. Compass rose of Tp and Dp during the northeast and southwest monsoon.
Energies 15 08637 g012
Figure 13. Wind rose diagram of wind velocities using MOV data (upper right panel) and PAGASA data (lower right panel) with the location map of both stations.
Figure 13. Wind rose diagram of wind velocities using MOV data (upper right panel) and PAGASA data (lower right panel) with the location map of both stations.
Energies 15 08637 g013
Figure 14. Location of the six (6) points of interest within the model domain (stations A–F).
Figure 14. Location of the six (6) points of interest within the model domain (stations A–F).
Energies 15 08637 g014
Figure 15. (a) Hs (m) model during northeast monsoon season and (b) Hs model indicating the average Hs at stations A–F.
Figure 15. (a) Hs (m) model during northeast monsoon season and (b) Hs model indicating the average Hs at stations A–F.
Energies 15 08637 g015
Figure 16. (a) Tp (s) model during northeast monsoon season and (b) Tp model indicating the average Tp at stations A–F.
Figure 16. (a) Tp (s) model during northeast monsoon season and (b) Tp model indicating the average Tp at stations A–F.
Energies 15 08637 g016
Figure 17. (a) Wave energy propagation (J/m2) model during northeast monsoon season and (b) wave energy propagation model, indicating the total wave energy propagated and indicating the average period at stations A–F.
Figure 17. (a) Wave energy propagation (J/m2) model during northeast monsoon season and (b) wave energy propagation model, indicating the total wave energy propagated and indicating the average period at stations A–F.
Energies 15 08637 g017
Figure 18. Hs (m) and peak wind directions during northeast monsoon season.
Figure 18. Hs (m) and peak wind directions during northeast monsoon season.
Energies 15 08637 g018
Figure 19. (a) Hs (m) model during southwest monsoon season and (b) Hs model indicating the average Hs at stations A–F.
Figure 19. (a) Hs (m) model during southwest monsoon season and (b) Hs model indicating the average Hs at stations A–F.
Energies 15 08637 g019
Figure 20. (a) Tp (s) model during southwest monsoon season and (b) Tp model indicating the average Tp at stations A–F.
Figure 20. (a) Tp (s) model during southwest monsoon season and (b) Tp model indicating the average Tp at stations A–F.
Energies 15 08637 g020
Figure 21. (a) Wave energy (J/m2) model during southwest monsoon season and (b) wave energy model indicating the total wave energy at stations A–F during the season.
Figure 21. (a) Wave energy (J/m2) model during southwest monsoon season and (b) wave energy model indicating the total wave energy at stations A–F during the season.
Energies 15 08637 g021
Figure 22. Scatter plot graph of stations A–F graphically showing the correlation of the two parameters.
Figure 22. Scatter plot graph of stations A–F graphically showing the correlation of the two parameters.
Energies 15 08637 g022
Table 1. Mean Hs, Tp and Dp from MOV’s 2008–2018 data for the northeast (December–January–February (DJF)/Amihan) and southwest (June–July–August (JJA)/Habagat) monsoons.
Table 1. Mean Hs, Tp and Dp from MOV’s 2008–2018 data for the northeast (December–January–February (DJF)/Amihan) and southwest (June–July–August (JJA)/Habagat) monsoons.
StationDJF Mean (Amihan)JJA Mean (Habagat)
HsTpDpHsTpDp
41.11315.259246.62630.50664.1881−167.453
80.86014.899128.29880.50444.2172−150.142
121.1064.846138.09070.49724.3267−153.027
141.11395.313562.7730.48294.4801−164.637
Table 2. Correlation between hindcast Hs data from MOV and PAGASA measured wind speed.
Table 2. Correlation between hindcast Hs data from MOV and PAGASA measured wind speed.
StationDistance from PAGASA—Cuyo Station (km)Correlation Factor (r) Average Annual Hs (m)Average Annual Pd (kW/m)
4680.711.344.28
7400.61.112.66
8660.511.173.06
9560.661.163.05
10150.751.203.13
11600.761.445.00
12920.751.384.25
13720.731.444.88
14920.621.404.88
Hs and Pd were computed using the hindcast data from MOV.
Table 3. Studies conducted that validate MOV’s data.
Table 3. Studies conducted that validate MOV’s data.
MOV’s Wave Hindcast ValidationResultsYearSource
Hindcast significant wave height was validated against (1) altimeters in the entire Tanaki grid area and (2) wave rider buoy moored ~17 m deep near the port of Tanaki.1. Model Hs does not reveal significant or systematic errors and shows satisfactory overall performance of the wave model.2021[45]
2. Results of the model data validation against the buoy were also in good agreement, especially when considering the proximity of the buoy to shore and the mooring’s relatively shallow depth (Hs, r = 0.933, Tp, r = 0.624).
The mean monthly distribution of significant wave height from MOV was validated against the wave data from FugroGeos for the Central–Southern North Sea region.1. The coefficient of correlation is satisfactory at r = 0.55.2021[46]
Assessment of the suitability of the hindcast directional wave spectra data in the northern and southern hemisphere was compared against buoy observations from National Oceanic and Atmospheric Observation’s National Data Buoy Center.1. The Hs time series obtained from the hindcast and the observations are well correlated (rHs > 0.8 for all the selected buoys).2019[47]
2. The RMSE values are within 10% of the maximum value attained by the model for all the selected buoys
3. The monthly mean spectra from model and observations are well correlated, with r ~0.9.
A report on meteorological and oceanographic conditions at Astrolabe Reef prepared by MetOcean Solutions Limited has been reviewed by NIWA. These statistics are derived from numerical model simulations, which were first validated against available data for winds, waves and currents1. The verification of significant wave heights against measurements is satisfactory. It gives some confidence in the accuracy of the hindcast winds, in the absence of direct verification at a nearby location in the Bay of Plenty.2013[48]
Table 4. Summary results of wave parameters for the NE and SW monsoons at the six (6) points of interest near the island (simulated using MOV data).
Table 4. Summary results of wave parameters for the NE and SW monsoons at the six (6) points of interest near the island (simulated using MOV data).
StaLonLatDJF-NE MonsoonJJA-SW Monsoon
Hs (m)Tp (s)Total Wave Energy (J/m2)Pd (W/m)Hs (m)Tp (s)Total Wave Energy (J/m2)Pd (W/m)
A121111.4284.77811249.043978.250.42582.5658111.0763193.8773
B120.8510.91.25474.361964.33822789.880.54493.0525181.8621379.0649
C12110.671.19624.2348876.51342361.110.54373.0582181.0445366.3609
D121.1810.781.30774.53991047.543250.230.56573.1477195.9928420.4111
E121.27111.4274.7461247.263910.130.54583.0908182.4629388.3967
F120.8210.81.48874.8711357.534246.800.49092.834147.609280.2246
Table 5. Summary results of wave parameters for the NE and SW monsoons at the nine (9) MOV stations.
Table 5. Summary results of wave parameters for the NE and SW monsoons at the nine (9) MOV stations.
StaLonLatDJF-NE MonsoonJJA-SW Monsoon
Hs (m)Tp (s)Total Wave Energy (J/m2)Pd (W/m)Hs (m)Tp (s)Total Wave Energy (J/m2)Pd (W/m)
14120.511.51.40424.70921207.843968.370.56543.0135195.8246419.0981
1312111.51.44084.66421271.413891.820.51542.7885162.6904299.367
12121.511.51.34.36251035.142968.290.55522.9305188.8014366.2484
11120.5111.48854.77511357.164280.960.58493.0785209.5412431.2502
10121111.51354.83531403.074513.250.52742.8493170.3466331.3767
9121.5111.45054.64811288.673923.060.58353.0309208.5396420.899
4120.510.51.54874.87731469.174690.350.60253.1567222.3052484.0884
712110.51.50584.81471388.734501.620.58893.0644212.4356436.9373
8121.510.51.43174.65551255.573879.120.57593.0252203.1603411.3036
Table 6. Summary of statistical metrics between “observed values” and “hindcast parameters”.
Table 6. Summary of statistical metrics between “observed values” and “hindcast parameters”.
Station x ¯ y ¯ BiasRMSESIr
A0.400.800.400.551.380.92
B0.370.790.410.541.450.88
C0.350.770.370.541.580.85
D0.410.820.410.541.320.90
E0.450.860.420.561.240.91
F0.470.860.380.511.080.94
Table 7. AEP results of Wave Dragon, AWS and SSG, in MWh.
Table 7. AEP results of Wave Dragon, AWS and SSG, in MWh.
Wave Energy ConverterAnnual Energy Production (AEP), in MWhCapacity Factor (%)
ABCDEF
Wave Dragon, 5900 kW1626.651399.191280.771705.381970.611937.013.1
Wave Dragon, 7000 kW2046.151805.071657.622174.612462.042427.663.4
Archimedes Wave Swing, 2470 kW50.75423.11818.83444.44562.42461.5840.2
Seawave Slot-Cone Generator, 20,000 kW3322.452618.6432268.443466.913170.8924099.231.8
Table 8. Sample computation of AEP for Wave Dragon (7000 kW) at Station F.
Table 8. Sample computation of AEP for Wave Dragon (7000 kW) at Station F.
Annual Energy Production
Wave Dragon, 7000 kWTp, in seconds
456789
5678910
Hs,
in meters
0121.121.51.081.081.081.08
121094396.9
23 365.4236.7
34 170.467.62
45 69.15
AEP = 2427.48 MWh
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pacaldo, J.C.; Bilgera, P.H.T.; Abundo, M.L.S. Nearshore Wave Energy Resource Assessment for Off-Grid Islands: A Case Study in Cuyo Island, Palawan, Philippines. Energies 2022, 15, 8637. https://doi.org/10.3390/en15228637

AMA Style

Pacaldo JC, Bilgera PHT, Abundo MLS. Nearshore Wave Energy Resource Assessment for Off-Grid Islands: A Case Study in Cuyo Island, Palawan, Philippines. Energies. 2022; 15(22):8637. https://doi.org/10.3390/en15228637

Chicago/Turabian Style

Pacaldo, Jonathan C., Princess Hope T. Bilgera, and Michael Lochinvar S. Abundo. 2022. "Nearshore Wave Energy Resource Assessment for Off-Grid Islands: A Case Study in Cuyo Island, Palawan, Philippines" Energies 15, no. 22: 8637. https://doi.org/10.3390/en15228637

APA Style

Pacaldo, J. C., Bilgera, P. H. T., & Abundo, M. L. S. (2022). Nearshore Wave Energy Resource Assessment for Off-Grid Islands: A Case Study in Cuyo Island, Palawan, Philippines. Energies, 15(22), 8637. https://doi.org/10.3390/en15228637

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop