Multi-Objective Sizing Optimization of Hybrid Renewable Energy Microgrid in a Stand-Alone Marine Context

: Islands are the main platforms for exploration and utilization of marine resources. In this paper, an island hybrid renewable energy microgrid devoted to a stand-alone marine application is established. The speciﬁc microgrid is composed of wind turbines, tidal current turbines, and battery storage systems considering the climate resources and precious land resources. A multi-objective sizing optimization method is proposed comprehensively considering the economy, reliability and energy utilization indexes. Three optimization objectives are presented: minimizing the Loss of Power Supply Probability, the Cost of Energy and the Dump Energy Probability. An improved multi-objective grey wolf optimizer based on Halton sequence and social motivation strategy (HSMGWO) is proposed to solve the proposed sizing optimization problem. MATLAB software is utilized to program and simulate the optimization problem of the hybrid energy system. Optimization results conﬁrm that the proposed method and improved algorithm are feasible to optimally size the system, and the energy management strategy effectively matches the requirements of system operation. The proposed HSMGWO shows better convergence and coverage than standard multi-objective grey wolf optimizer (MOGWO) and multi-objective particle swarm optimization (MOPSO) in solving multi-objective sizing problems. Furthermore, the annual operation of the system is simulated, the power generation and economic beneﬁts of each component are analyzed, as well as the sensitivity. and 13a,b. The proposed HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal front of HSMGWO is more distributed and broader than that of MOGWO and MOPSO. HSMGWO has obvious performance advantages in optimizing system size. The results indicate that the proposed HSMGWO has good performance in solving such kind of multi-objective problems.


Introduction
Under the multiple pressures of economic growth, energy shortage and environmental protection, marine resources will be an important guarantee for human survival and sustainable development in the coming century. The 13th Five-Year Plan for China's Marine Economic Development [1] points out that the ocean is an important strategic space for China's economic and social development. It is necessary to seize the great opportunity of the "21st Century Maritime Silk Road" [2], adhere to the basic principle of "green development and ecological priority", strengthen the development and construction of islands, accelerate the development and utilization of marine renewable energy. China has more than 3.4 million km 2 of sea area and 7300 islands larger than 0.5 km 2 , including more than 400 inhabited islands, facing a serious shortage of water and electricity supply [3,4]. Islands are indispensable platforms for shipping industry, global trade, national security, scientific research and marine exploration. It has great significance to study and build a clean and reliable island microgrid. Nevertheless, the proper sizing of system components is a key factor for system's techno-economic feasibility.
Governments and research institutes around the world have carried out a series of studies, demonstrations and experiments on sizing optimization of island hybrid renewable energy systems [5]. Eco-grid EU project had been built and studied on the Danish island of Bornholm for smart grid demonstration, including a hybrid wind, biogas, photovoltaic, cogeneration, diesel, fuel cell and storage bank energy system [6]. Literature [7] introduced the energy system integration of Maui Island, Hawaii, which consists of oil-fired power plants, large diesel-fueled units, biomass-fueled power plant, wind plants, grid-connected PV systems, and battery energy storage system. A large amount of wind power produced on Maui is curtailed because the power cannot be accommodated by the existing power grid. Literature [8] shared the real case of the Caribbean island of St. Eustatius, presented the design and performance results of a MW-sized photovoltaic (PV), battery storage system (BSS), and diesel genset hybrid system, which has been now in operation for over 2 years. Literature [9] analyzed the economic benefits of replacing pure diesel generators with a hybrid photovoltaic diesel system in a remote island of the Philippines and showed that the energy cost could be reduced by 70%. Research on sizing optimization and system construction of island hybrid energy microgrids had also been carried out by many other countries and regions according to the literature [10][11][12].
In terms of sizing optimization objectives and methods, Abo-Elyousr et al. presented a method for the optimal configuration of a hybrid microgrid system, which investigated a bi-objective based optimization algorithm to reduce cost, emissions, and a combined solution between cost and emissions [13]. Tawil et al. presented a case study of a hybrid microgrid on the island of Ouessant in France, including wind turbines, tidal turbines, diesel generators and pumped storage power stations. The main sizing constraint is to achieve a compromise between the hybrid system overall cost, for a period of 15 years, and the CO 2 emission minimization [14]. Mahesh et al. proposed a photovoltaic, wind, storage grid connected system, the genetic algorithm was used to solve the sizing optimization problem with the objective of minimizing the average energy cost [15]. Zhu et al. proposed a sizing optimization method based on improved grey wolf optimization algorithm for island hybrid energy microgrid. The optimization objective was to compromise the minimum annualized cost of system and the minimum deficiency of power supply probability [16]. Mahesh et al. proposed an advanced Pareto-front non-dominated sorting multi-objective particle swarm optimization method for optimal configuration (placement and sizing) of distributed generation (DG) in a radial distribution system [17]. Yilmaz et al. optimized a hybrid energy system composed of photovoltaic, diesel generator and battery bank by using the Homer software, and the optimal scale of the system was solved with the objective of minimum investment cost [18].
The existing research on sizing optimization and calculation method of island hybrid renewable energy systems has made significant progress but shows some limitations in the following aspects. Firstly, the present power supply of offshore islands mostly depends on diesel generators, which causes high cost of fuel transportation and serious environmental pollution. This is contrary to the principles of low-carbon development and ecological protection advocated worldwide [19]. Secondly, current sizing optimization methods of hybrid energy system mainly use single objective optimization, some use bi-objective optimization, few studies consider tri-objective optimization. Island energy system sizing is an extremely complex problem, which requires comprehensive consideration of various factors in order to adapt to local conditions and obtain the most appropriate optimization results. Thirdly, in marine conditions, in view of the special geographical environment of islands, few studies have considered the contradiction between the lack of land resources and the land demand for generation plants construction. Most of the current studies only consider solar energy, onshore wind energy, hydro energy, pumped storage energy, solar thermal energy, biomass energy and other conventional land-based energy, while ignoring the abundant marine energies, such as tidal current energy, which does not suffer the scarcity of land needed for building the generation plants [20]. Taking the South China Sea as an example, the area of most major islands in the South China Sea is no more than 0.5 km 2 . Many islands come from reclamation, having no possibility to build large-scale land-based power stations [21]. The strategic status of islands highlights its important roles in marine economy, marine exploration, scientific research and militarization. Precious land resources should be the key consideration in sizing optimization of island microgrid. With the improvement of marine energy technology and equipment manufacturing technology [22], marine energy holds the highest opportunity given the wide abundance of the energy source, and has become the best choice on the table for building island renewable energy microgrid.
In this paper, based on the current research results, land demand is taken as one of the key considerations of component selection in the energy system. An island hybrid renewable energy microgrid composed of wind turbines, tidal current turbines, and battery storage systems is constructed based on the abundant wind and tidal energy resources in the ocean. A tri-objective optimization model including minimizing the Loss of Power Supply Probability (LPSP), the Cost of Energy (COE) and the Dump Energy Probability (DEP) is established. An improved Multi-Objective Grey Wolf Optimizer based on Halton sequence and social motivation strategy (HSMGWO) is proposed to solve the proposed sizing optimization problem. MATLAB software is utilized to program and simulate the optimization problem of the hybrid energy system. A case study is carried out and analyzed to verify the feasibility of the sizing optimization method.
This paper is organized as follows. Section 2 describes the mathematical model of each component of the hybrid energy system. The energy management strategy, description of objective functions and constraints, and modeling of optimization problems are presented in Section 3. Section 4 is dedicated to the improvement of multi-objective grey wolf optimizer to solve the tri-objective optimization problem, and its performance verification. The simulation results are discussed in Section 5 to illustrate the performance of the proposed method, the benefit analysis and sensitivity analysis of each component are carried out. The conclusion is presented in the last section. The workflow of this paper is presented in Figure 1.

Selection and Modeling of System Components
In this paper, an island in Zhoushan Archipelago of China is taken as a case study to build an island hybrid renewable energy microgrid. The primary task is to select the

Selection and Modeling of System Components
In this paper, an island in Zhoushan Archipelago of China is taken as a case study to build an island hybrid renewable energy microgrid. The primary task is to select the appropriate energy combination scheme according to the natural renewable resources. Figure 2 shows the map of the studied site [23]. Considering the limitation of the precious island land resources, the offshore generation components should be selected in the first instance. Therefore, land-based photovoltaic, solar thermal, hydropower, pumped storage, tidal energy and other conventional renewable energies that need a lot of land resources will not be given priority. The available marine renewable energy resources of the studied site mainly include solar energy, wind energy, tidal current energy and wave energy. According to the evaluation of references [24,25], the island is rich in wind energy resources with great wind energy quality. A large tidal energy reserve exists in the sea area, which has great exploration and utilization value. Figure 3 presents the annual wind speed curve of a typical year at the studied site. The historical data are obtained from the database of NOAA's National Centers for Environmental Information [26]. Figure 4 presents the annual tidal current speed curve, which is simulated according to the historical data in reference [25]. It is worth mentioning that the solar energy and wave energy resources of the studied site have certain exploration and utilization value. However, offshore floating photovoltaic and wave energy generation would not be taken into consideration because of the low technology maturity and high energy cost [27].

Resources Characterization and Components Selection
In this study, offshore wind energy and tidal current energy are selected as distributed generation sources. Meanwhile, in order to avoid the continuous utilization of diesel generators, the mature battery storage system is selected to stabilize the fluctuation of renewable energy and improve the system reliability and resources utilization rate. The proposed island renewable energy microgrid is presented in Figure 2, which consists of offshore wind turbines (WT), tidal current turbines (TCT), battery storage systems (BSS) and other indispensable components such as inverters and transformers. The detailed modeling of the system major components is presented hereafter.
Electronics 2021, 10, x FOR PEER REVIEW 5 of 24 diesel generators, the mature battery storage system is selected to stabilize the fluctuation of renewable energy and improve the system reliability and resources utilization rate. The proposed island renewable energy microgrid is presented in Figure 2, which consists of offshore wind turbines (WT), tidal current turbines (TCT), battery storage systems (BSS) and other indispensable components such as inverters and transformers. The detailed modeling of the system major components is presented hereafter.

Wind Turbine Model
The model of power output generated by the wind turbine can be expressed as the following equation.
where C pw is the power coefficient of the wind turbine, ρ w is the air density, r is the radius of wind turbine blades, v w is the wind speed. The relationship between the power output of wind turbine and wind speed can be expressed as a piecewise function, which is determined by different wind speed intervals. When the wind speed is in the range of cut-in speed and rated speed, the output power is the product of rated power and efficiency coefficient. When the wind speed is in the range of rated speed and cut-out speed, the wind turbine outputs rated power [28].
where v is the wind speed, v ci , v co and v r are the cut-in wind speed, rated wind speed and cut-out wind speed, P r is the rated electrical power.

Tidal Current Turbine Model
The astronomical nature of tides is determined by the gravitational interaction between the earth, the moon and the sun. There are two flood tides, and two ebb tides every day. The spring tide and neap tide appear periodically. The strength of tides varies according to the distance, relative position and angle between the moon, the sun and the earth. The power generated by a tidal current turbine can be estimated by the following equation [29].
where C p is the power coefficient of TCT, which is essentially the ratio of the mechanical power obtained by TCT, taking into account the conversion efficiency. ρ is the density of seawater, A is the cross-sectional area of TCT, v tid is the tidal current velocity. The periodic variation of tidal current velocity v tid can be simplified as the following equation.
where T is the period between spring tide and neap tide, T is the period of tidal current, V m is the average velocity of tidal current calculated by (V s + V n )/2, V s and V n are the maximum velocity of spring tide and neap tide, K is the wave coefficient of tidal current, calculated by (V s − V n )/(V s + V n ).

Battery Storage System Model
Energy storage system can smooth the fluctuation of distributed generation, enhance the system schedulability and improve the reliability of system operation. BSS is the most flexible, reliable and mature system among various storage technologies. The capacity of the battery storage system can be expressed by the following equation [30].
where A d is the desired autonomy day, which means the maximum number of days that the battery can provide continuous energy without charging, P l is the daily power demand, η bat is the battery efficiency, η inv is the conversion efficiency of the inverter, DOD is the maximum allowable depth of discharge. SOC is the battery state of charge, which is usually used to indicate the remaining capacity of the battery. The SOC(t) at time t is given by [31]: where σ and η bat are the self-discharge rate and efficiency, P b (t) is the difference between total generation power and the demand power P l considering the efficiency factor, given by:

Optimization Objectives
The purpose of establishing renewable energy microgrid with WTs, TCTs and BSSs is to provide electricity for isolated islands and ensure that the energy system has high reliability level, low energy cost and high energy utilization rate. Review literatures [32] summarized the sizing optimization objectives and methods of hybrid energy microgrid, among which economy, reliability, energy utilization and pollution indexes are most commonly utilized. In this paper, the Loss of Power Supply Probability (LPSP), the Cost of Energy (COE) and the Dump Energy Probability (DEP) are selected as optimization objectives.

Loss of Power Supply Probability
Reliability is the basis of system operation. In this paper, the Loss of Power Supply Probability (LPSP) is used as the evaluation index of system reliability. LPSP indicates the Electronics 2021, 10, 174 7 of 24 extent to which power generation fails to match the load demand. LPSP can be defined as total power supply failure time divided by the total reporting time [33]: where P loss (t) is the probability that the load exceeds the power supply at time t, ∆t is the step size that is usually taken as 1h; T is the total duration that is generally taken as 8760h.

Cost of Energy
The Cost of Energy (COE) is the ratio of the annual cost of the system to the power generation. COE comprehensively considers the system investment, operation and maintenance costs, annual power generation, and economic factors such as annual interest rate. COE can reasonably reflect the economic benefits of the system. The mathematical model can be expressed as follows [34]: where x i is the number of component i, n is the number of component types, I C,i is the annualized capital cost, R C,i is the annualized replacement cost, O MC,i is the annualized operation and maintenance cost, P O,i is the annual power generation. CRF is the capital recovery factor and SFF is the sinking fund factor calculated by the following equations.
where r is the discount rate, which is calculated as 6.7% according to the interest rate of the PBOC, y is the lifespan of the project, which is 20 years in this paper, y i is the lifespan of the component.

Dump Energy Probability
Due to the inherent intermittently of renewable energy, excessive pursuit of reliability will lead to large-scale abandonment of wind energy and tide current energy. The more balanced way of optimization is to take the renewable energy utilization rate into consideration [35]. In this paper, it is measured by Dump Energy Probability (DEP), which can be calculated by the following equation.
where P dp (t) is the dumped power of energy system at time t, P re (t) is the total power generated by energy system at time t.

Objective Function and Constraints
Based on the above discussion, the optimization objective function of island renewable energy system based on minimum LPSP, minimum COE and minimum DEP is as follows: The multiple objectives of sizing optimization are subjected to a number of constraints including supply-demand balance, generation unit boundaries and the BSS constraints. These constraints are explained in detail below.

Variable Constraints
The system variables, i.e., the number of system components (WT, TCT, BSS), should be defined with reasonable maximum value, expressed as:

Generation Unit Boundaries
The output of each generation unit at any time shall be within the range of minimum and maximum values, expressed as:

Supply-Demand Balance
The total power provided by the microgrid system should always match the load demand, expressed as:

BSS Constraints
The SOC of the BBS remains between the minimum and maximum values, which can be expressed as:

Energy Management Strategy
The criteria for the establishment of island renewable energy microgrids include maximum utilization of renewable energy, reliable power supply, reduction of energy cost, etc. The inherent intermittently of renewable energy will lead to the mismatch between the total power generated by the system and the load demand. In some cases, the load demand may be exceeded, and in other cases, the load demand may not be matched. Energy management of hybrid energy system is the basic guarantee for economical operation. BSS can store the surplus energy when the power supply is excessive, and timely supplement the power shortage when the power supply is insufficient. The energy management strategy proposed in this paper is as follows.

1.
Renewable energy power generation is enough to match the load demand.
• If the BSS is not fully charged at this time, the excess power generation will charge the BSS.

•
If the BSS is fully charged at this time, the renewable energy output will be reduced to match the load. Record the dump energy.

2.
Renewable energy power generation is not enough to match the load demand.
• If the BSS is available at this time, the BSS will discharge to match the load demand.

•
If the BSS is not available at this time, the unsatisfiable load will be cut off. Record the failure time of power supply.

Improved Multi-Objective Grey Wolf Optimizer
The sizing optimization of a hybrid energy system is a typical multi-objective optimization problem that contains mutually exclusive optimization objectives, multiple variables and constraints. Multi-objective optimization can be formulated as a maximization or minimization problem. The solution set of multi-objective optimization problem can be represented by Pareto optimal front, which is a set containing the corresponding objective values of Pareto optimal solutions. In order to obtain the Pareto optimal front of the multi-objective sizing optimization problems, an improved Multi-Objective Grey Wolf Optimizer based on Halton sequence and social motivation strategy (HSMGWO) is proposed.

Multi-Objective Grey Wolf Optimizer
Grey wolf optimizer (GWO) is a powerful meta-heuristic algorithm [36], which has competitiveness in terms of solution accuracy, minimum computational effort, and aversion of premature convergence. According to the social grade of grey wolf colony, GWO divides the wolf pack into four pyramid population hierarchies, namely, α, β, δ, and ω. The global optimization is realized by imitating the predatory behaviors of grey wolves such as searching, tracking, encircling, chasing and attacking. In the GWO, the position of a wolf is the vector of system variables, and the prey corresponds to the objective function. The algorithm simulates the hunting process by constantly updating the location of wolf pack (set of the variables) in each iteration to obtain the Pareto optimal front (non-dominated solutions set). The optimization process is mainly guided by the best three solutions (α, β, and δ). The mathematical model of predatory behaviors is expressed as follows: where t is number of iterations, A and C are the coefficient vectors; X p is the position vector of the prey, X(t) is the position vector of the wolf, Max iter is the total number of iterations allowed for the optimization, r 1 and r 2 are the random vector in [0, 1]. Figure 5 illustrates the location updating mechanism of grey wolf. The mathematical model can be expressed as follows: where X α , X β , X δ are the positions of α, β, δ, X represents the position of the wolf, D α , D β , D δ , respectively represent the distance between the current candidate and the three optimal wolves. In order to solve the multi-objective optimization problem, Mirjalili et al. developed a multi-objective grey wolf optimizer (MOGWO) by integrating two new components into GWO algorithm [37]. The first component is the external archive, which stores the nondominated Pareto optimal solutions obtained so far. The second component is the leader selection mechanism, which selects the appropriate leader from all the non-dominated optimal solutions obtained in the external archive by roulette. The probability of being selected can be calculated by the following equation.
where C is a constant greater than 1, and can be configured according to requirements, N j is the total number of individuals in this segment.
are the positions of , , , represents the position of the wo , , respectively represent the distance between the current candidate and the optimal wolves. In order to solve the multi-objective optimization problem, Mirjalili et al. deve a multi-objective grey wolf optimizer (MOGWO) by integrating two new componen GWO algorithm [37]. The first component is the external archive, which stores the dominated Pareto optimal solutions obtained so far. The second component is the selection mechanism, which selects the appropriate leader from all the non-dom optimal solutions obtained in the external archive by roulette. The probability of selected can be calculated by the following equation.

Improvement of MOGWO
MOGWO algorithm is essentially an intelligent algorithm based on population. The quality of initial population has a great impact on the optimization effect and solution set. However, population diversity cannot be guaranteed by random initialization. Meanwhile, the wolf pack location is updated strictly by the guidance of α, β, δ wolves, which makes the algorithm easy to fall into local optimum. Aiming at these weaknesses, this paper introduces Halton sequence initialization and social motivation strategy to propose an improved multi-objective grey wolf optimization algorithm (HSMGWO).

Halton Sequence Initialization
In this paper, Halton sequence is used to improve population initialization instead of the original random initialization. Halton sequence is a well-known multi-dimensional low-discrepant sequence which can generate evenly distributed set of points in the search space [38]. An M-dimensional Halton sequence can be expressed as follows: where p i is a prime number, φ p i (n) is the base inverse function of p i , which can be expressed as follows: For a grey wolf population with dimension M, the location range is [x min , x max ], h m is a random number generated by Halton sequence, and the individual initial position can be expressed as follows: x The two-dimensional point sets generated by the random method and Halton sequence are shown in Figure 6. Compared with the random method, the individuals generated by Halton sequence are more evenly distributed without aggregation, and the numbers of individuals in different intervals are more average. The results show that Halton sequence can generate individuals with better distribution uniformity.
The two-dimensional point sets generated by the random method and Halton sequence are shown in Figure 6. Compared with the random method, the individuals generated by Halton sequence are more evenly distributed without aggregation, and the numbers of individuals in different intervals are more average. The results show that Halton sequence can generate individuals with better distribution uniformity.

Social Motivation Strategy
In MOGWO, wolves are only guided by three leaders ( , , ). This monotonous location update strategy will lead to the loss of population diversity, poor global

Social Motivation Strategy
In MOGWO, wolves are only guided by three leaders (α, β, δ). This monotonous location update strategy will lead to the loss of population diversity, poor global exploration ability and slow convergence speed. In this paper, social motivation strategy is introduced to update the position of the wolf pack in each iteration referring to the idea of Chimp optimization algorithm (ChOA) [39]. Social motivation includes social status competition and sexual motivation, which leads to the chaos of gray wolves. The wolves may not accept the guidance of the leaders and release their responsibilities during the hunt. This chaotic behavior helps the algorithm to improve the ability to jump out of local optimum and deal with high-dimensional problems. According to the evaluation of reference [40], the Gauss/mouse chaotic map with an initial value of 0.7 is selected to improve the algorithm, which can be expressed as follows: In order to simulate this social motivation strategy, it is assumed that there is a 50% probability for wolves to choose between the normal update model and the chaotic update model in the optimization process. The mathematical model is represented by the following equation.
where µ is a random number between [0, 1], m is the chaotic vector calculated based on the Gauss/mouse chaotic map.

Verification of HSMGWO
Two standard multi-objective test functions are selected to verify the performance of HSMGWO. The benchmark problems are provided as follows. DTLZ1: Electronics 2021, 10, 174 12 of 24 DTLZ2: In order to evaluate the performance of the improved algorithm, the HSMGWO is compared to the MOGWO and the well-known algorithm MOPSO with the performance index of Inverted Generational Distance (IGD). IGD can reflect the convergence and coverage of the algorithm. The smaller value of IGD indicates that the calculation result of the algorithm is closer to the real Pareto front. The mathematical formula of this index is as follows.
where n is the number of true Pareto optimal solutions and d i indicates the Euclidean distance between the ith true Pareto optimal solution and the closest obtained Pareto optimal solution in the reference set.
where is the number of true Pareto optimal solutions and indicates the Euclidean distance between the th true Pareto optimal solution and the closest obtained Pareto optimal solution in the reference set.
The  Obviously, the improved HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal solutions of HSMGWO are closer to the true Pareto optimal front and highly distributed along all objectives. Although there where is the number of true Pareto optimal solutions and indicates the Euclidean distance between the th true Pareto optimal solution and the closest obtained Pareto optimal solution in the reference set.
The  Obviously, the improved HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal solutions of HSMGWO are closer to the true Pareto optimal front and highly distributed along all objectives. Although there are discontinuities on the Pareto optimal front of HSMGWO, the coverage of the whole Obviously, the improved HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal solutions of HSMGWO are closer to the true Pareto optimal front and highly distributed along all objectives. Although there are discontinuities on the Pareto optimal front of HSMGWO, the coverage of the whole front is broader on all the test functions. The solution sets provide superior convergence and uniform distribution, due to the improvement of the Halton Sequence Initialization and the social motivation strategy. The IGD results are listed in Table 1. In the initialization stage of MOGWO, the position of wolf pack is generated randomly, which can easily lead to uneven population distribution, and seriously affect the solution quality and optimization efficiency. The population initialization method improved by the low-differential Halton sequence can obtain a more uniform initial population, significantly improve the population diversity and the optimization efficiency of the algorithm, effectively avoid falling into local optimization, and ensure the basic solution quality of the algorithm.
Improved by social motivation strategy, the position updating with a certain probability of chaotic state as the external performance prevents wolves from quickly becoming trapped in local optimum, improves the global exploration ability and convergence efficiency of the algorithm, and significantly improves the optimization effect of the algorithm for multi-dimensional complex optimization problems. In essence, social motivation strategy provides soft transition between global and local search ability. In this way, wolves tend to broadly discover promising regions of search space and exploit the best one. These superior results are based on the special form of the Gauss/mouse map. In early stage it has large and extremely variable amplitude, while its amplitude and variableness decrease severely in the final stages. Wolves change abruptly in the early stages of the hunting process and then gradually converge. However, there is no additional computational cost for the proposed algorithm.

9.
Add the non-dominated solutions to the archive according to the archiving rules and remove the dominated solutions. Remove the excess solutions when the population of external archive is full. 10. Use Max iter to judge whether the algorithm should be terminated. The algorithm ends and all the non-dominated solutions are output if Max iter is reached. Conversely, return to Step 6. social motivation strategy provides soft transition between global and loc In this way, wolves tend to broadly discover promising regions of search s the best one. These superior results are based on the special form of the Ga In early stage it has large and extremely variable amplitude, while its variableness decrease severely in the final stages. Wolves change abru stages of the hunting process and then gradually converge. Howe additional computational cost for the proposed algorithm.

Flowchart of Optimization
The flowchart of sizing optimization using HSMGWO is shown in F

Case Studies
The proposed sizing optimization method has been applied to analyze a hybrid renewable energy system of an island in Zhoushan Archipelago, China. The average load of the island microgrid is 200 kW, and the peak load is about 500 kW. The load curve is shown in Figure 10. At present, the traditional diesel generators are used to supply electricity to the island residents, causing serious environmental pollution. The stability of island energy system is seriously affected by diesel transportation.
The proposed sizing optimization method has been applied to analyze a hybrid renewable energy system of an island in Zhoushan Archipelago, China. The average load of the island microgrid is 200 kW, and the peak load is about 500 kW. The load curve is shown in Figure 10. At present, the traditional diesel generators are used to supply electricity to the island residents, causing serious environmental pollution. The stability of island energy system is seriously affected by diesel transportation. The distributions of annual wind speed and tidal current speed are shown in Figure  11 according to the historical data. The available hours of wind and tidal energy at the study site are considerable, which account for 80.42 and 77.73% of the annual 8760 h, respectively. Wind and tidal current energy can be used to generate electricity for most of the year. Due to the high availability of resources, conventional offshore wind turbine and tidal current turbine can be selected as the main components of the energy system. As a The distributions of annual wind speed and tidal current speed are shown in Figure 11 according to the historical data. The available hours of wind and tidal energy at the study site are considerable, which account for 80.42% and 77.73% of the annual 8760 h, respectively. Wind and tidal current energy can be used to generate electricity for most of the year. Due to the high availability of resources, conventional offshore wind turbine and tidal current turbine can be selected as the main components of the energy system. As a special research object with low load level, island microgrid is vulnerable to marine natural disasters. In order to ensure the reliability of the system, components with appropriate rated capacity should be considered to avoid large-scale power outage caused by equipment failure. The selection of system components refers to the project data of State Grid Corporation of China, and is adjusted according to references [41,42]. The economic and technical parameters of the major components are shown in Table 2.

R PEER REVIEW
16 of 24 special research object with low load level, island microgrid is vulnerable to marine natural disasters. In order to ensure the reliability of the system, components with appropriate rated capacity should be considered to avoid large-scale power outage caused by equipment failure. The selection of system components refers to the project data of State Grid Corporation of China, and is adjusted according to references [41,42]. The economic and technical parameters of the major components are shown in Table 2.

Results and Discussions
The sizing optimization of island renewable energy system based on HSMGWO has been implemented by the MATLAB (R2018b, MathWorks, Natick, MA, USA) simulation platform on an 8-core Lenovo computer (T490, Lenovo, Beijing, China) with 16 GB memory and 3.4 GHz CPU clock speed, running Windows 10 enterprise operating system (64-bit). The parameters of HSMGWO are shown in Table 3. The optimization process is carried out in hourly steps to simulate the annual operation of the energy system. Different sizing configurations correspond to different variable inputs.

Comparison of Pareto Fronts
In order to further test the performance of the algorithm, HSMGWO, MOPSO and MOGWO algorithms are utilized to solve the optimization objectives. The Pareto fronts obtained by HSMGWO, MOPSO and MOGWO are shown in Figures 12 and 13a,b. The proposed HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal front of HSMGWO is more distributed and broader than that of MOGWO and MOPSO. HSMGWO has obvious performance advantages in optimizing system size. The results indicate that the proposed HSMGWO has good performance in solving such kind of multi-objective problems.
In order to further test the performance of the algorithm, HSMGWO, MOPSO and MOGWO algorithms are utilized to solve the optimization objectives. The Pareto fronts obtained by HSMGWO, MOPSO and MOGWO are shown in Figures 12 and 13a,b. The proposed HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal front of HSMGWO is more distributed and broader than that of MOGWO and MOPSO. HSMGWO has obvious performance advantages in optimizing system size. The results indicate that the proposed HSMGWO has good performance in solving such kind of multi-objective problems.  MOGWO algorithms are utilized to solve the optimization objectives. The Pareto fronts obtained by HSMGWO, MOPSO and MOGWO are shown in Figures 12 and 13a,b. The proposed HSMGWO shows better convergence and coverage than standard MOGWO and MOPSO. The Pareto optimal front of HSMGWO is more distributed and broader than that of MOGWO and MOPSO. HSMGWO has obvious performance advantages in optimizing system size. The results indicate that the proposed HSMGWO has good performance in solving such kind of multi-objective problems.   Figure 12 shows the relationship between LPSP, COE and DEP of the studied hybrid renewable energy system. Table 4 lists the optimal solutions of each objective and the optimal equilibrium solution of the system. The optimal equilibrium solution is the solution with the minimum value after normalizing all the non-dominated solutions, which is calculated by the following equation.

Optimization Results
where f i is the normalized value of solution i, f i (n) is the function value of objective n of solution i, f min (n) denotes the minimum function value of objective n, f max (n) denotes the maximum function value of objective n.  19.6%, the corresponding DEP is 13.5%. Solution 3: The minimum DEP is 4.4%, the corresponding LPSP is 26.3%, the corresponding COE is 0.19 $/kWh. Solution 4 (the optimal equilibrium solution): The LPSP is 9.8%, the corresponding COE is 0.18 $/kWh, the corresponding DEP is 16.6%. Obviously, LPSP, COE and DEP are three contradictory objectives of energy system sizing. Reliability is guaranteed at the expense of high energy cost and resource dumping. In addition, this study simulated a microgrid of 100% renewable energy penetration, without diesel generators and conventional landbased power plants considering the scarcity of island land resources. Although it is in line with the principle of low-carbon and sustainable development, the reliability of the renewable energy system cannot reach the level of conventional power grid, which can be seen from the higher average value of LPSP in the results. Therefore, the accurate decision is based on the fully evaluation of the actual demand and the detailed analysis of Pareto solution set. From the perspective of decision makers, it is necessary to balance the weight of reliability, economy and resource utilization, further determine the appropriate type and size of system components to avoid excessive investment, substandard reliability, or energy dumping. Further analysis is based on Solution 4. In this case, the components of hybrid renewable energy system include four wind turbines, three tidal current turbines and five battery storage systems. Figure 14 demonstrates the annual operation simulation results of hybrid renewable energy system. Figure 15 demonstrate the zoomed version for one typical week. The result shows that TCTs provide stable power output, while WTs provide intermittent power output, which is determined by the characteristics of different resources. WTs provide higher peak of power output than TCTs. The BSS and distributed generation form an effective complementarity. When the distributed generation is sufficient, the BSS is charged to improve the utilization of renewable energy. When the distributed generation is insufficient, the BSS discharges to provide power output. The system realizes energy translation in time dimension. The energy management strategy fully matches the operation requirement of island microgrid and ensures reliable and sustainable power supply.      Figure 16 shows the total electricity generation of WTs and TCTs in each month, with the monthly load as reference. The results show that WTs provide more electricity than TCTs. Although the power curve demonstrates the volatility of wind energy, the monthly power generation of WTs is relatively average due to the abundant wind energy resources in the studied site. The TCTs provide stable power output, and the monthly power generation is uniformly distributed with its resource curve. The total power output of WTs and TCTs exceeds the load demand in each month. However, the renewable energy and load both have inherent randomness, the load cannot always be matched with the  Figure 16 shows the total electricity generation of WTs and TCTs in each month, with the monthly load as reference. The results show that WTs provide more electricity than TCTs. Although the power curve demonstrates the volatility of wind energy, the monthly power generation of WTs is relatively average due to the abundant wind energy resources in the studied site. The TCTs provide stable power output, and the monthly power generation is uniformly distributed with its resource curve. The total power output of WTs and TCTs exceeds the load demand in each month. However, the renewable energy and load both have inherent randomness, the load cannot always be matched with the power output of energy system, which leads to lower reliability compared with grid-connected microgrid. power output of energy system, which leads to lower reliability compared with gridconnected microgrid.  Figure 17 shows the generation and investment percentages of different components in the system. WTs contributed 49.83% power output with 48.47% investment, TCTs contributed 29.21% power output with 35.07% investment and BSSs contributed 20.96% power output with 16.46% investment. Obviously, WTs and BSS have the best investment benefit, the investment benefit of TCTs is ordinary. The utilization of BSSs significantly  Figure 17 shows the generation and investment percentages of different components in the system. WTs contributed 49.83% power output with 48.47% investment, TCTs contributed 29.21% power output with 35.07% investment and BSSs contributed 20.96% power output with 16.46% investment. Obviously, WTs and BSS have the best investment benefit, the investment benefit of TCTs is ordinary. The utilization of BSSs significantly improves the economy of the system.  Figure 17 shows the generation and investment percentages of different components in the system. WTs contributed 49.83% power output with 48.47% investment, TCTs contributed 29.21% power output with 35.07% investment and BSSs contributed 20.96% power output with 16.46% investment. Obviously, WTs and BSS have the best investment benefit, the investment benefit of TCTs is ordinary. The utilization of BSSs significantly improves the economy of the system. In order to further study the impact of different components on different optimization objectives in isolated island microgrid, sensitivity analysis is carried out by changing the size of various components. The quantitative sensitivity analysis results of WTs, TCTs and BSSs are shown in Figures 18, 19 and 20  In order to further study the impact of different components on different optimization objectives in isolated island microgrid, sensitivity analysis is carried out by changing the size of various components. The quantitative sensitivity analysis results of WTs, TCTs and BSSs are shown in Figures 18-20, respectively.

Components Sensitivity Analysis
The sensitivity analysis results show that LPSP decreases with the increase in all components sizes. For WTs, LPSP decreases from the maximum 52.56% to the minimum 3.91% as the size increases. For TCTs, LPSP decreases from the maximum 29.57% to the minimum 3.81% as the size increases. For BSSs, LPSP decreases from the maximum 29.59% to the minimum 8.24% as the size increases. The size of WTs has a great impact on LPSP, while the sizes of TCTs and BSSs have the relatively smaller impact on LPSP. This further reflects that WTs and TCTs provided the most basic and indispensable power output and have an incredible contribution to the reliability of renewable energy system. At the same time, BSSs can well suppress the fluctuation and randomness of renewable energy and reduce the risk of power shortage.         The size of WTs affects COE seriously. For renewable energy system without WTs, the COE can reach as high as 1.72 $/kWh. Once WTs is incorporated into the system, the COE will reduce greatly to 0.39 $/kWh and then decreases steadily to the minimum 0.16 $/kWh with the increase in size. On the contrary, the COE increases steadily from the minimum 0.12 $/kWh to the maximum 0.24 $/kWh as the size of TCTs increases, which is related to the high cost of tidal current power generation. With the development of technology and the policy guidance of new energy, the cost of tidal current power generation is expected to be reduced to the same level as offshore wind power generation. The size of BSSs has a tiny impact on COE with the minimum 0.17 $/kWh and the maximum 0.18 $/kWh. Taking 2 BSSs as the watershed, the COE of the system decreases slightly first and then increases slightly with the increase in the size of BSSs.
The results show that DEP increases with the increase in the sizes of WTs and TCTs. For WTs, DEP increases from the minimum 0.41% to the maximum 27.03% as the size increases. For TCTs, DEP increases from the minimum 2.11% to the maximum 42.28% as the size increases. Different characteristics of renewable energies determine the different impact on DEP. The size of TCTs contributes more to the increase in DEP than that of WTs due to the stability of tidal current energy. DEP decreases from the maximum 32.75% to the minimum 15.56% as the size of BSSs increases, which reflects that BSSs can greatly reduce energy dumping and improve the utilization of renewable energy.

Conclusions
In this paper, taking land demand as one of the key considerations of component selection, an island hybrid renewable energy microgrid, composed of wind turbines, tidal current turbines, battery storage systems, is established in a stand-alone marine context. A tri-objective sizing optimization method is proposed, comprehensively considering the economy, reliability and energy utilization indexes. Three optimization objectives are presented: minimizing the Loss of Power Supply Probability, the Cost of Energy and the Dump Energy Probability. An improved Multi-Objective Grey Wolf Optimizer based on Halton sequence and social motivation strategy (HSMGWO) is proposed to solve the sizing optimization problem. The system is modeled and simulated by MATLAB software. The major conclusions in this paper are listed as follows: 1.
The simulation results verify the feasibility of the proposed optimization method. The optimal size of the island renewable energy system is obtained. WT, TCT, and BBS are closely complementary base on an effective and reliable energy management strategy.

2.
The simulation results verify the feasibility and advancement of the proposed HSMGWO in solving the multi-objective sizing optimization problem of renewable energy system. The HSMGWO shows better convergence and coverage of Pareto front compared with MOGWO and MOPSO. 3.
The case study proves that it is feasible to construct an island micro grid composed of offshore renewable energy components on the basis of specific resource evaluation, taking valuable land resources as one of the primary considerations.

4.
The results of operation and sensitivity analysis indicate that in the studied area, the investment and power generation of various components of the microgrid are relatively balanced. WTs provide the most power and have the greatest impact on the economy and reliability indexes. TCTs and BSSs have a great impact on the reliability index. BSSs can well shift the system's generation and steadily improve the reliability and resource utilization.
The multi-objective sizing optimization method of island renewable system proposed in this paper is also applied to other similar optimization problems of isolated microgrid. This paper has presented a methodological study, further investigation of operation strategy should be carried out to improve the system reliability and reduce energy cost.
Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the confidentiality of involved data.

Conflicts of Interest:
The authors declare no conflict of interest.