Wave energy array development siting method

Abstract

The application discloses a wave energy array development site selection method, and belongs to the field of wave energy development and utilization, comprising the following steps: 1) processing geographic distribution data by using a spatial analysis method of GIS to obtain a restricted development area layer; 2) constructing an evaluation index system for wave energy array development site selection, and performing standardization processing to obtain a comprehensive evaluation result; 3) calculating an offshore distance to generate a refined developable range layer; 4) carrying out statistical analysis on a wave direction in the developable range to obtain an array development range; 5) determining a potential development site on the basis of the array development range of a target sea area; and 6) taking the potential development site as a center point to perform quantitative analysis and trend analysis on array development capacity, and determining a key development site. The application integrates multidimensional evaluation indexes such as environmental compatibility, resource endowment, technical feasibility and economic rationality, and provides a scientific basis for wave energy power station array development site selection.

Description

Technical Field

[0001] This invention relates to a wave energy development and utilization technology, and more particularly to a wave energy array development and site selection method. Background Technology

[0002] With economic and social development, the consumption of fossil fuels such as coal and oil has intensified, causing serious pollution to the ecological environment. Developing renewable energy sources such as ocean energy, solar energy, and wind energy, and increasing the proportion of clean energy, are important measures to solve the problems of energy shortage and environmental pollution. Among them, wave energy has advantages such as high energy density and wide distribution, and has been widely studied and applied both domestically and internationally.

[0003] However, wave energy is unstable and exhibits seasonal and regional variations. Therefore, assessing the development potential and stability of wave energy resources, and scientifically delineating the development scope of wave energy resources, can provide important data support for the site selection of wave energy power plants.

[0004] Chinese patent CN118396433A discloses a site selection method for joint development of marine energy, including: 1) determining the latitude and longitude range of the simulation area and dividing it into grids, and using ERA5 reanalysis wind field data to drive the SWAN wave numerical model to obtain wave field data; 2) determining the evaluation indicators for joint wind and wave development, mainly including four aspects: resource reserves, resource stability, resource complementarity, and natural conditions; 3) standardizing each indicator and calculating the weight of each indicator using the analytic hierarchy process (AHP); 4) calculating the joint wind and wave development coefficient and drawing a cloud map; 5) based on the drawn cloud map, delineating potential development areas, selecting representative test points in each area, and using cross-spectral analysis to compare and select the optimal site for joint wind and wave development. This patent establishes a systematic evaluation indicator system for joint wind and wave development, uses cross-spectral analysis to determine the optimal development site, and proposes a relatively comprehensive optimized site selection framework for joint wind and wave development. However, the following shortcomings exist: (1) The patent focuses on the joint development of wind and waves and does not provide a systematic method for the site selection of wave energy array development; (2) Although an evaluation index system has been established, there is a lack of relevant content on the preprocessing of the research sea area based on the exclusion index.

[0005] Chinese patent CN116341952A discloses a method for assessing the joint development of marine wind and wave resources, comprising the following steps: 1) collecting wind energy resource data and wave energy resource data of the sea area to be assessed, and determining the power data of wind energy and wave energy based on the wind energy resource data and wave energy resource data; 2) determining the inherent stability of wind energy and wave energy based on the wind energy resource data and wave energy resource data; 3) determining the correlation between wind energy and wave energy based on the wind energy resource data and wave energy resource data, and determining the complementarity of wind energy and wave energy; 4) determining the potential level of joint development of offshore wind energy and wave energy based on the inherent stability, correlation, and complementarity of wind energy and wave energy. This patent proposes a framework for assessing the joint development of marine energy based on resource stability, correlation, and complementarity, providing methodological support for the quantitative assessment of joint development potential. However, there are the following shortcomings: (1) The research mainly focuses on potential assessment, lacking specific site selection and in-depth quantitative analysis; (2) A complete process from macro-constraint screening to refined site selection has not been established; (3) The research focuses on resource assessment and does not take into account the engineering and technical requirements involved in array development. Summary of the Invention

[0006] The purpose of this invention is to overcome the shortcomings of the existing technology and provide a site selection method for the array-based development of wave energy. This method innovatively constructs a comprehensive evaluation model for the site selection of array-based wave energy development based on a Geographic Information System (GIS). This model breaks through the limitations of traditional site selection methods, systematically integrating multi-dimensional evaluation indicators such as environmental compatibility, resource endowment, technical feasibility, and economic rationality, providing a scientific basis for the site selection of array-based wave energy power plants.

[0007] To achieve the above objectives, the present invention adopts the following technical solution:

[0008] A method for site selection for array-based development of wave energy includes:

[0009] 1) Use GIS spatial analysis methods to process geographic distribution data, exclude restricted areas, areas with unsuitable water depth and distance from the shore, and areas with excessively high or low wave heights, to obtain a restricted development area layer;

[0010] 2) Construct an evaluation index system for the selection of wave energy array development sites, standardize the system, and sum the weighted values ​​of each index to obtain a comprehensive evaluation result;

[0011] 3) Calculate the offshore distance, and within the area less than the set threshold, exclude areas with low comprehensive evaluation results and restricted development areas to generate a refined development range layer;

[0012] 4) Conduct wave direction statistical analysis within the exploitable range to exclude areas with overly dispersed wave direction distribution, thereby further determining the array-based exploitable range;

[0013] 5) Based on clarifying the target sea area for arrayed development, conduct site selection studies for arrayed wave energy development and identify potential development sites;

[0014] 6) Using potential development sites as the center, rationally delineate development areas, conduct quantitative and trend analysis of array-based development capabilities, and identify key development sites.

[0015] In step 2), the evaluation indicators include wave power density, the proportion of available wave height, coefficient of variation, and monthly variation index. The weights of each indicator are calculated using the analytic hierarchy process.

[0016] In step 3), the offshore distance is calculated using the Haversine formula:

[0017] The offshore distance D is calculated using the Haversine formula:

[0018] ,

[0019] Δ φ = φ 2 - φ 1 ,

[0020] Δ λ = λ 2 - λ 1 ,

[0021] in, r For the Earth's radius, φ 1 and φ 2 These are the latitudes of the two points (in radians). λ 1 and λ 2 These are the longitudes of the two points (in radians).

[0022] In step 4), firstly, the wave directions within the target sea area are statistically analyzed, dividing the wave directions into N, NNE, NE, etc., with each wave direction interval being 22.5°. The frequency of the dominant wave direction is statistically analyzed, and wave direction areas below a set threshold are excluded from the exploitable range. Simultaneously, the frequency difference and direction difference between the dominant and secondary dominant wave directions are further calculated. Areas with a frequency difference less than the threshold and a direction difference greater than the threshold are also excluded from the exploitable range, further obtaining the array-based exploitable range. Areas with concentrated wave directions are retained, providing favorable conditions for the array-based deployment of the device. The exclusion criteria can be expressed as:

[0023] ,

[0024] in, P m The frequency of the main wave direction, P s The frequency of the secondary main wave direction, P th,m The frequency threshold of the main wave direction, P th,d The frequency difference threshold between the primary wave direction and the secondary primary wave direction. θ m The main wave is oriented towards the center angle. θ s The angle of the secondary main wave towards the center. θ th,d This is the threshold for the difference in center angles between the primary and secondary primary wave directions. The threshold is determined through statistical analysis of the wave directions in the target sea area.

[0025] Building upon this foundation, the study couples marine environmental dynamic resources with the reconstructed energy-concentrating structure, considering reflection and transmission characteristics. Small-scale simulations are conducted within the array-based development range to overcome the limitations of pure resource assessment and further evaluate the resource characteristics of the reconstructed wave field. In the SWAN model, by setting different combinations of reflection and transmission coefficients for the energy-concentrating structure, the actual gain effect of the structure on wave energy resources under various combinations is systematically analyzed. This establishes a mapping relationship from laboratory conditions to complex actual sea conditions, providing a quantitative basis for practical engineering design. Given that the research focuses on the enhancement effect on wave energy resources within the influence range of the energy-concentrating structure, the absolute size of the structure has a relatively small impact on the research conclusions; therefore, the interference caused by size factors can be ignored in the analysis.

[0026] Energy gain is quantitatively evaluated by comparing the rate of change in wave power density before and after the placement of the energy-concentrating structure, thus assessing the improvement effect of array-based development. Energy gain coefficient. β The calculation formula is:

[0027] ,

[0028] in, P wave,e To determine the wave power density after the energy-concentrating structure is installed, P wave,s The wave power density is the value of the wave without a condenser structure.

[0029] In step 5), the spatial scale of the arrayed deployment is first determined based on the number of wave energy conversion devices to be deployed. Considering the reconfiguration effect of the energy-concentrating structure on the local wave field, the actual theoretical potential of wave energy resources is calculated using the following formula: P a =(1+β ) P wave B，

[0030] in, P a For actual theoretical reserves, B For the width of the array deployment, P wave The wave power density.

[0031] By comparing the scale, energy conversion efficiency, and wave field resource reconstruction distribution characteristics of typical wave energy device arrays, key parameters and site selection criteria for array development are determined. A spatial analysis model is constructed, and an evaluation index ψ for array-based development is proposed, with the calculation formula as follows:

[0032] ,

[0033] In the formula, P r This represents the minimum required installed capacity for array deployment. η The energy conversion efficiency of the array arrangement.

[0034] Within the scope of array-based development, potential development sites that meet the minimum installed capacity requirements are selected by calculating array-based development evaluation indicators. This selection criterion not only meets the basic requirements for resource development scale but also prioritizes areas with high resource conversion efficiency and outstanding economic viability. Ultimately, the selected sites will possess significant resource advantages and development potential, providing a scientific basis and key decision support for subsequent detailed design and implementation, while ensuring the sustainability and economic rationality of the development project.

[0035] In step 6), the wave energy array development capability is quantitatively analyzed using the spectral width index to further clarify the resource potential of the selected area. ,

[0036] In the formula, ε 0 represents spectral width. m y For wave spectrum y Step moment, m 0 The zeroth moment of the wave spectrum, m -1 The -1st moment of the wave spectrum, m -2 The -2nd moment of the wave spectrum, S For the spectrum, f For frequency, Δ f For frequency intervals, N The number of discrete parts of the frequency.

[0037] Spectral width, as a key indicator, is used to characterize the diffusion characteristics of wave energy in frequency distribution. By calculating and comparing the spectral width of various potential development areas, the resource development potential of these areas can be further verified and refined.

[0038] A long-term trend analysis of the effective wave height was conducted for each potential development area. This trend analysis focuses on predicting future resource changes and helps assess the future sustainable development potential of each area. The analytical method is as follows:

[0039] (1) Theil-Sen Median trend analysis:

[0040] Theil-Sen Median method, also known as Sen slope estimation, is a robust nonparametric statistical trend calculation method. Its calculation method is as follows: ,

[0041] in, n For the sample size, H s For the effective wave height, M Median refers to the median of the slopes of n(n-1) / 2 data combinations; Median means to calculate the median.

[0042] (2) Mann-Kendall test:

[0043] The Mann-Kendall test is a non-parametric method for testing trends in time series data. It does not require that the measurements follow a normal distribution and is unaffected by missing or outlier values. The Mann-Kendall test is suitable for testing the significance of trends in long-term data series. The test process involves determining the random arrangement of data in the series and whether a significant upward or downward trend exists. The test statistic S... m The calculation formula considers the magnitude relationship between data pairs and uses a two-sided trend test to determine the significance of the trend. Its test statistic St m The calculation method is as follows:

[0044] ,

[0045] Where sgn is the sign function:

[0046] ,

[0047] When the statistic S m When the sequence follows a normal distribution, and without considering the existence of data points with equal values, its mean is 0, and the variance of the statistic is calculated by the following formula:

[0048] ,

[0049] Use test statisticK The trend test is performed using the following calculation method:

[0050] ,

[0051] By applying a two-tailed test, for a specified confidence level α The significance of the trend can be assessed; at that time This indicates that the time series is at a confidence level. α There is a significant trend; when Values ​​greater than 1.65, 1.96, and 2.58 indicate that the trend has passed the significance tests at confidence levels of 90%, 95%, and 99%, respectively. This method provides a given confidence level... α The value is 0.05, and the threshold is 1.96.

[0052] The beneficial effects of this invention are:

[0053] 1) By using GIS spatial analysis methods, areas with unsuitable constraints, water depth, and distance from the shore, as well as areas with excessively high or low wave heights, were excluded. This enabled the integrated processing of multi-dimensional constraints, laying a solid foundation for subsequent refined analysis.

[0054] 2) By introducing wave direction statistical analysis methods, quantitative analysis of the frequency distribution characteristics of the main and secondary main wave directions is conducted to systematically identify and retain the advantageous areas where wave directions are concentrated, providing more favorable marine dynamic environment conditions for the arrayed arrangement of wave energy devices.

[0055] 3) By coupling the reflection and transmission characteristics of the energy-concentrating structure with the SWAN wave numerical model, the energy gain effect prediction and evaluation from laboratory conditions to complex sea conditions were realized. This solved the technical limitation that pure resource assessment could not consider the influence of the device, and provided a reliable theoretical basis and technical support for the prediction of the engineering application and array development effect of the energy-concentrating structure.

[0056] 4) Combining the number of wave energy devices to be deployed and the spatial scale requirements, the entire process from regional selection to specific site determination has been optimized. While meeting the minimum installed capacity requirements, sites with high resource conversion efficiency and outstanding economic benefits are given priority, providing a scientific site selection method and decision support tool for large-scale wave energy development. Attached Figure Description

[0057] Figure 1 This represents the wave power density in the embodiment of the invention without a focused energy structure.

[0058] Figure 2 The wave power density in this embodiment of the invention is the area where the energy-concentrating structure is placed.

[0059] Figure 3 This refers to the energy gain ratio in the embodiments of the present invention;

[0060] Figure 4 This is the wave height variation trend of potential development area 1 in this embodiment of the invention;

[0061] Figure 5 This shows the wave height variation trend of potential development area 2 in this embodiment of the invention. Detailed Implementation

[0062] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0063] The structures, proportions, and sizes illustrated in the accompanying drawings are merely for illustrative purposes and to aid those skilled in the art in understanding and reading the invention. They are not intended to limit the scope of the invention and therefore have no substantial technical significance. Any modifications to the structure, changes in proportions, or adjustments to size, provided they do not affect the effectiveness or purpose of the invention, should still fall within the scope of the technical content disclosed herein. Furthermore, the terms "upper," "lower," "left," "right," "middle," and "one" used in this specification are merely for clarity and not intended to limit the scope of the invention. Changes or adjustments to their relative relationships, without substantially altering the technical content, should also be considered within the scope of the invention's implementation.

[0064] The simulation area in this embodiment covers latitude and longitude ranges of 117°E-127°E and 35°N-41°N, using a rectangular grid with a resolution of [resolution missing]. Wave field simulation was performed using the SWAN wave numerical model driven by ERA5 reanalysis wind field data, obtaining wave parameters for the entire year from March 1, 2023 to February 29, 2024. Furthermore, the method proposed in this invention was used to conduct research on the selection of wave energy array development sites.

[0065] Using GIS spatial analysis methods, restricted areas, areas with unsuitable water depth and offshore distance, and areas with excessively high or low wave heights were excluded to obtain a restricted development area layer. Specifically, the restricted areas must fully consider the marine spatial planning requirements of Shandong Province, Hebei Province, Tianjin Municipality, and Liaoning Province. According to relevant planning documents, important functional areas such as ecological protection zones and shipping channels were first excluded to avoid adverse impacts on the ecosystem and shipping safety, ensuring that the project layout is consistent with the regional marine functional positioning.

[0066] An evaluation index system for the selection of wave energy array development sites was constructed, mainly including wave power density, the proportion of available wave height, coefficient of variation, and monthly variation index, and standardized. The weights of each index were calculated using the analytic hierarchy process (AHP) to be 0.466, 0.254, 0.175, and 0.105, respectively, and then weighted and summed to obtain the comprehensive evaluation result.

[0067] Calculate the offshore distance using the Haversine formula:

[0068] ,

[0069] Δ φ = φ 2 - φ 1 ,

[0070] Δ λ = λ 2 - λ 1 ,

[0071] in, r For the Earth's radius, φ 1 and φ 2 These are the latitudes of the two points (in radians). λ 1 and λ 2 These are the longitudes of the two points (in radians).

[0072] In this embodiment, the offshore distance threshold is set to 50km, and within areas smaller than this threshold, areas with low comprehensive evaluation results and restricted development areas are excluded to generate a refined development range layer.

[0073] Statistical analysis of wave directions within the target sea area was conducted, dividing wave directions into N, NNE, NE, etc., with each wave direction interval being 22.5°. The frequency of the dominant wave direction was statistically analyzed, and areas with wave directions accounting for less than 15% of the exploitable range were excluded. Simultaneously, the frequency and direction differences between the dominant and secondary dominant wave directions were further calculated. Areas with a frequency difference less than 4% and a direction difference greater than 90° were also excluded from the exploitable range, further determining the array-based exploitable range. Areas with concentrated wave directions were retained to provide favorable conditions for the array-based deployment of the equipment. The exclusion criteria can be expressed as: ,

[0074] in, P m The frequency of the main wave direction, P s The frequency of the secondary main wave direction, P th,m The frequency threshold of the main wave direction, P th,d The frequency difference threshold between the primary wave direction and the secondary primary wave direction. θ m The main wave is oriented towards the center angle. θ s The angle of the secondary main wave towards the center.θ th,d This is the threshold for the difference in center angles between the primary and secondary primary wave directions. The threshold is determined through statistical analysis of the wave directions in the target sea area.

[0075] Building upon this foundation, the study couples marine environmental dynamic resources with the reconstructed energy-concentrating structure, considering reflection and transmission characteristics. Small-scale simulations are conducted within the array-based development range to overcome the limitations of pure resource assessment and further evaluate the resource characteristics of the reconstructed wave field. In the SWAN model, by setting different combinations of reflection and transmission coefficients for the energy-concentrating structure, the actual gain effect of the structure on wave energy resources under various combinations is systematically analyzed. This establishes a mapping relationship from laboratory conditions to complex actual sea conditions, providing a quantitative basis for practical engineering design. Given that the research focuses on the enhancement effect on wave energy resources within the influence range of the energy-concentrating structure, the absolute size of the structure has a relatively small impact on the research conclusions; therefore, the interference caused by size factors can be ignored in the analysis.

[0076] Energy gain is quantitatively evaluated by comparing the rate of change in wave power density before and after the placement of the energy-concentrating structure, such as... Figures 1 - 3 As shown, without the wave-focusing structure, the wave power density in the selected area is around 6 kW / m. After installing the wave-focusing structure, the maximum wave power density within its influence range reaches 10 kW / m, with the energy gain ratio generally between 20% and 60%. The formula for calculating the energy gain coefficient is: Energy Gain Coefficient β The calculation formula is:

[0077] ,

[0078] in, P wave,e To determine the wave power density after the energy-concentrating structure is installed, P wave,s The wave power density is the value of the wave without a condenser structure.

[0079] Based on the number of wave energy conversion devices to be deployed, the spatial scale of the array deployment is determined. Considering the reconfiguration effect of the energy-concentrating structure on the local wave field, the actual theoretical potential of wave energy resources is calculated using the following formula:

[0080] P a =(1+ β ) P wave B，

[0081] in, P a For actual theoretical reserves, B For the width of the array deployment,P wave The wave power density.

[0082] By comparing the scale, energy conversion efficiency, and wave field resource reconstruction distribution characteristics of typical wave energy device arrays, key parameters and site selection criteria for array development are determined. A spatial analysis model is constructed, and an evaluation index ψ for array-based development is proposed, with the calculation formula as follows:

[0083] ,

[0084] In the formula, P r This represents the minimum required installed capacity for array deployment. η The energy conversion efficiency of the array arrangement.

[0085] Within the scope of array-based development, potential development sites that meet the minimum installed capacity requirements are selected by calculating array-based development evaluation indicators. This selection criterion not only meets the basic requirements for resource development scale but also prioritizes areas with high resource conversion efficiency and outstanding economic viability. Ultimately, the selected sites will possess significant resource advantages and development potential, providing a scientific basis and key decision support for subsequent detailed design and implementation, while ensuring the sustainability and economic rationality of the development project.

[0086] Based on the above analysis, two areas with high development potential were selected from the research sea area. Table 1 lists the geographical coordinates of these two potential development areas.

[0087] Table 1 Geographic coordinate information of potential development areas

[0088] area longitude range latitude range 1 122.85°-122.90°E 37.35°-37.40°N 2 121.50°-121.55°E 36.15°-36.20°N

[0089] By using the spectral width index, a quantitative analysis of the wave energy array development capability is conducted to further clarify the resource potential of the selected area.

[0090] ,

[0091] In the formula, ε 0 represents spectral width. m y For wave spectrum y Step moment, m 0 The zeroth moment of the wave spectrum, m -1 The -1st moment of the wave spectrum, m -2 The -2nd moment of the wave spectrum, S For the spectrum, f For frequency, Δ f For frequency intervals, NThe number of discrete parts of the frequency.

[0092] Spectral width, as a key indicator, is used to characterize the diffusion characteristics of wave energy in frequency distribution. By calculating and comparing the spectral width of various potential development areas, the resource development potential of these areas can be further verified and refined.

[0093] Long-term trend analysis of effective wave height in various potential development areas is conducted. This trend analysis focuses on predicting future resource changes, helping to assess the future sustainable development potential of each region, such as... Figure 4 and Figure 5 As shown. The analysis method is as follows:

[0094] (1) Theil-Sen Median trend analysis:

[0095] Theil-Sen Median method, also known as Sen slope estimation, is a robust nonparametric statistical trend calculation method. Its calculation method is as follows: ,

[0096] in, n For the sample size, H s For the effective wave height, M Median refers to the median of the slopes of n(n-1) / 2 data combinations. Median means to calculate the median.

[0097] (2) Mann-Kendall test:

[0098] The Mann-Kendall test is a non-parametric method for testing trends in time series data. It does not require that the measurements follow a normal distribution and is unaffected by missing or outlier values. The Mann-Kendall test is suitable for testing the significance of trends in long-term data series. The test process involves determining the random arrangement of data in the series and whether a significant upward or downward trend exists. The test statistic S... m The calculation formula considers the magnitude relationship between data pairs and uses a two-sided trend test to determine the significance of the trend. Its test statistic St m The calculation method is as follows:

[0099] ,

[0100] Where sgn is the sign function:

[0101] ,

[0102] When the statistic S m When the sequence follows a normal distribution, and without considering the existence of data points with equal values, its mean is 0, and the variance of the statistic is calculated by the following formula:

[0103] ,

[0104] Use test statistic K The trend test is performed using the following calculation method:

[0105] ,

[0106] By applying a two-tailed test, for a specified confidence level α The significance of the trend can be assessed; at that time This indicates that the time series is at a confidence level. α There is a significant trend; when Values ​​greater than 1.65, 1.96, and 2.58 indicate that the trend has passed the significance tests at confidence levels of 90%, 95%, and 99%, respectively. This method provides a given confidence level... α The value is 0.05, and the threshold is 1.96.

[0107] While the specific embodiments of the present invention have been described above in conjunction with the accompanying drawings, this is not intended to limit the scope of protection of the present invention. Those skilled in the art should understand that various modifications or variations that can be made by those skilled in the art without creative effort based on the technical solutions of the present invention are still within the scope of protection of the present invention.

Claims

1. A method of site selection for arrayed exploitation of wave energy, characterized in that, Comprise: (1) The spatial analysis method of GIS is used to process the geographical distribution data, and the restricted area, water depth, and the area with inappropriate offshore distance, as well as the area with excessively high or low wave height are excluded to obtain the restricted development area layer; (2) An evaluation index system for wave energy array development site selection is constructed, and the indexes are standardized and weighted to obtain the comprehensive evaluation result; the evaluation indexes include wave power density, available wave height proportion, coefficient of variation, and monthly variation index, and the weights of the indexes are calculated by using the analytic hierarchy process; (3) The offshore distance is calculated, and in the area less than the set threshold, the area with low comprehensive evaluation result and the restricted development area are excluded to generate the fine developable range layer; The offshore distance D adopts the Haversine formula: ， wherein r R is the radius of the earth, (4) In the developable range, statistical analysis of wave direction is carried out, and the area with excessively dispersed wave direction distribution is excluded to further obtain the array development range; 1 and Firstly, the wave direction in the target sea area is statistically analyzed, and the wave direction is divided according to the direction, and each wave direction interval is 22.5 degrees; the frequency of the main wave direction is counted, and the wave direction area with a frequency less than the set threshold is excluded in the developable range; at the same time, the frequency difference and direction difference between the main wave direction and the secondary main wave direction are further calculated, and the area with a frequency difference less than the threshold and a direction difference greater than the threshold is also excluded in the developable range to further obtain the array development range; the area with concentrated wave direction is retained to provide favorable conditions for wave energy array arrangement; 2 φi and φ2are the latitudes of the two points, in radians, The exclusion condition is expressed as: 1 and (5) On the basis of the array development range of the target sea area, the selection of wave energy array development site is studied to determine the potential development site; 2 λi and λ2are the longitudes of the two points, in radians. Firstly, the spatial scale of array arrangement is determined according to the number of wave energy conversion devices to be arranged; considering the reconstruction effect of the energy concentrating structure on the local wave field, the actual theoretical reserves of wave energy resources are calculated, and the calculation formula is: By comparing the scale layout, energy conversion efficiency, and wave field resource reconstruction distribution characteristics of the typical wave energy device array, the key parameter indexes and site selection standards of array development are determined, a spatial analysis model is constructed, and the array development evaluation index ψ is calculated, and the calculation formula is: In the array development range, the array development evaluation index is calculated to select the potential development site that meets the minimum installed capacity requirement; ， wherein, P m a frequency of the primary wave direction, P s a frequency of the secondary wave direction, P th,m a frequency threshold of the primary wave direction, P th,d a frequency difference threshold of the primary wave direction and the secondary wave direction, (6) The development area is reasonably delineated with the potential development site as the center point, the quantitative analysis and trend analysis of array development capacity are carried out, and the key development site is determined; m a central angle of the primary wave direction, ​ s a central angle of the secondary wave direction, ​ th,d a central angle difference threshold of the primary wave direction and the secondary wave direction. ​ ​ ， wherein, P wave,e Pw is the wave power density for the unplaced energy focusing structure, P wave,s Pw is the wave power density for the unplaced energy focusing structure, ​ ​ ， wherein, P a is the actual theoretical in-place reserves, B is the width of the arrayed layout, P wave,s is the wave power density without the placement of the energy focusing structure; ​ ， wherein P r lower limit of the required installed capacity for the array arrangement, ​ energy conversion efficiency for the array arrangement; ​ ​ The wave energy array development capacity is quantitatively analyzed by using the spectral width index to further clarify the resource potential of the selected area. ， In the formula, ε 0 represents spectral width. m y For wave spectrum y Step moment, m 0 The zeroth moment of the wave spectrum, m -1 The -1st moment of the wave spectrum, m -2 The -2nd moment of the wave spectrum, S i For the first i One spectrum f y i For the first i Frequency y The power, Δ f i For the first i One frequency interval, N The number of discrete parts of the frequency.

Images