Optimal Scheduling of a Microgrid Using Multiobjective Biogeography-Based Optimization Model and Algorithm with Adaptive Migration

Aiming at the important research topic of optimal scheduling in the microgrid field, the general model for multiobjective dynamic optimal scheduling of a microgrid is established with the objective of minimizing economic and environmental costs. On this basis, the model is organically integrated with constraint handling technology, multiobjective optimization, and biogeographybased optimization algorithm, and then a constrained multiobjective evolutionary model suitable for biogeography-based optimization is further established. ,e corresponding constraint handling mechanism, the determination method of habitat suitability index, and migration strategy are improved, and the convergence performance and the distribution uniformity of Pareto frontier for multiobjective evolutionary algorithm are effectively enhanced. Applied to the optimal scheduling of typical microgrid systems, the effectiveness of the proposed model and method is verified.


Introduction
Due to the global energy crisis and environmental pollution, the renewable clean energy such as solar and wind power has developed rapidly [1]. However, uncertain factors such as the intermittence and randomness of wind and light power output and the load fluctuation threaten the safe and stable operation of the power system [2,3]. Meanwhile, it is difficult to meet the power supply requirements in isolated islands or other remote areas. erefore, microgrid, the small power generation and distribution system combining distributed power supply, energy storage device, load, and protection device, can effectively solve the above problems [4][5][6]. How to properly manage the distributed generation of microgrid and maximize the economic and environmental benefits of microgrid has become an important research topic [7,8].
On account of the diversity of energy structure, distributed power types, and control methods in microgrid, the energy management and optimal operation of microgrid are of high complexity. Many scholars are committed to the research of microgrid energy management system. In [9], a simulation-based optimization method is proposed for the design of energy management system in photovoltaic gridconnected microgrid with heterogeneous occupancy. In [10], a sustainable scheduling problem for household microgrids with unscheduled distributed energy and responsive load demand is studied and the energy management system optimization algorithm based on mixed-integer linear programming (MILP) is proposed. In [11], the plug-in electric vehicles (PEVs) are introduced into microgrid, and a two-stage energy management strategy is proposed for the contribution of PEVs in the demand response plan of commercial building microgrid. In [12], a distributed requirements management system (DMS) is established to enhance any rule-based DMS through feedback behavior. Meanwhile, considering the multiobjective attribute of microgrid itself, the traditional single-objective optimal scheduling is difficult to adapt to the complex microgrid environment. At present, most of the research on the optimal operation of microgrid system focuses on the optimization of the model and the improvement of scheduling strategy and evolutionary strategy. In terms of solving the optimal scheduling model, the operation optimization of microgrid system is a nonlinear optimization problem with multiobjective and multiconstraints.
At present, the common intelligent optimization algorithms, genetic algorithm (GA) [13], particle swarming algorithm (PSO) [14], bacterial foraging algorithm (BFA) [15], and so on, are favored by scholars for their good parallelism, robustness, and global optimization ability but have their own advantages and disadvantages. It is still one of the problems worth studying to explore the optimization algorithm which is more suitable for solving the problem of optimal operation of microgrid. In [16], a multiobjective optimization model for microgrid operation management is built from three aspects of operating cost, pollutant emission, and operational risk level, and a niche multiobjective particle swarm optimization algorithm is proposed to solve the problem. In [17], a multiobjective economic scheduling model of cogeneration microgrid system is proposed, which considers the active and reactive power output and heat generation income. However, the above method is mainly optimized for the single objective and does not consider multiple objectives at the same time. In [18], the concept of energy and environmental benefits is introduced to modify the optimal scheduling model of power system including wind farms, and an improved particle swarm optimization algorithm combining tabu search and fuzzy processing technology is proposed. In [19], an improved Pareto genetic multiobjective optimization algorithm with vector module fitness function is proposed, which combines the idea of nondominated sorting, elitist strategy, and improved niche technology. In [20], a multiobjective optimization model of microgrid considering both the minimum O & M cost and the lowest environmental pollution is established, the multiobjective problem is transformed into a single-objective problem by using the linear weighted sum method, and the influence of different weight coefficients on the optimization results is discussed. In [21], the mathematical model of the optimal operation of microgrid based on improved gravitational search algorithm is proposed, taking into account several operation indexes such as economic costs, environmental benefits, and network losses. In [22], the multiobjective optimization model is solved by using the intensity Pareto evolutionary algorithm with three indexes of active network loss, pollution gas emission, and system voltage stability as optimization objectives. In [23], the mathematical model of microgrid operation is established, and the problem of optimal operation of microgrid is solved using improved adaptive firefly algorithm. However, the above method is mainly optimized for the single objective and does not consider multiple objectives at the same time. Numerous multiobjective intelligent optimization algorithms, such as nondominated sorting genetic algorithm (NSGA-II) [24], multiobjective particle swarm algorithm (MOPSO) [25][26][27], nondominated sorting differential evolution (NSDE) [28], and game theory algorithms [29], are, respectively, used to solve the multiobjective optimal scheduling of microgrid. Although these methods belong to multiobjective optimization, they still suffer from the problem of falling into the local optimum frontier, which is unable to take into account both convergence and uniformity of the distribution.
In this paper, a general model for multiobjective dynamic optimal scheduling of microgrid is established with the objective of minimizing economic and environmental costs on the premise of satisfying load demand and microsource constraints. On top of this, the model is organically integrated with constraint handling technology, multiobjective optimization, and biogeography-based optimization (BBO) [30] algorithm. en a constrained multiobjective evolutionary model suitable for BBO is further established and the corresponding constraint handling mechanism, the determination method of habitat suitability index, and adaptive migration strategy are effectively improved. Taking the typical microgrid system including photovoltaic array, wind turbine, fuel cell, microgas turbine, diesel engine, and storage battery as an example, the validity and reliability of the model and the algorithm proposed are analyzed.

Optimal Scheduling Model of the
Microgrid System e microgrid system used in this paper, which includes renewable and uncontrollable clean units such as photovoltaic array (PV) and wind turbine (WT), controllable nonclean units with high efficiency and small pollution such as fuel cell (FC), microgas turbine (MT), and diesel generator set (DE), and energy storage battery, is in island operation, and the internal load of microgrid is powered by each distributed generation unit. e goal of microgrid optimal scheduling is to optimize the output of microgrid power generation units under various constraints, so as to achieve the best economy and environmental protection of the system.

Wind Power Generation Model.
e specific formula of the output power of wind turbine is shown as follows: where v is the actual wind speed; v ci and v co are the cut-in and the cut-out wind speed. When the actual wind speed is lower than v ci or higher than v co , the turbine stops; P v is the rated output power of wind turbine within the normal wind speed.

Photovoltaic Generation Model.
Considering the randomness and intermittence of light and ambient temperature, the output power is calculated as follows: In the formula, P PV is the output power of photovoltaic cell; P s is the maximum output power under standard conditions; S A is the ambient light intensity; S o is the light intensity under standard conditions; k is power temperature coefficient; T C is the working temperature of battery board; T r is the reference temperature.

Microgas Turbine Fuel Cost Model.
e fuel cost of microgas turbine (MT) is shown by where F MT is the fuel cost of microgas turbine; C fuel is the price of natural gas; L is the low calorific value of natural gas; P MT and η MT are the output power of MT and efficiency.

Fuel Cell
Model. e energy consumption characteristic of the fuel cell model is similar to that of MT, and the fuel consumption can be expressed as shown in where F FC is the fuel cost of fuel cell; P FC and η FC are the output power and efficiency of FC.

Diesel Generator Fuel Cost
Model. e fuel cost of diesel generator is shown by where F DE is the fuel cost of diesel generator; α, β, and c are the fuel cost coefficient; P DE is the output power of DE.

Storage Battery Model.
In the economic dispatching of microgrid, the storage battery plays a role of peak load shifting in the load of microgrid and ensures continuous power supply of the system. e state of charge of a battery (SOC) is the ratio of its remaining capacity to its rated capacity. SOC at time t + 1 is determined by battery capacity and SOC at time t as follows; SOC charging formula can be expressed as where σ is the self-discharge rate; P ES,t is the charge-discharge power at time t; η IN is the charge efficiency; and C is total battery capacity within unit time Δt. SOC discharging formula can be expressed as where η OUT is the discharge efficiency.

Optimization
Model of the System. e purpose of microgrid optimal scheduling is to make the operating cost minimum of system and the environmental benefit maximum, by optimizing the output of each power supply in the microgrid with satisfying the load requirements and the constraints on the output of each distributed power.

Objective Function
(1) Minimize the Operating Cost e operating cost is an important index to measure the economic operation of microgrid, which mainly includes fuel cost, maintenance cost, and low-voltage distribution and online electricity cost. It is shown in the following formula: where F 1 is the operating cost of microgrid; T is the number of scheduling periods in the scheduling period; N is the number of distributed power generation units; C F,j,t and C OM,j,t are the fuel cost and maintenance cost of the i-th power generation unit at time t; C G,t is the cost of power interaction between microgrid and large grid at time t; n is the cost factor in the system, when running in island mode, n � 0.

(2) Minimize the Pollutant Treatment Cost
With the implementation of the strategy of sustainable development, environmental benefits are increasingly valued. In this paper, the pollution degree is defined to measure the harm degree of CO 2 , SO 2 , NO x , and other pollution gases to the environment caused by the equivalent discharge from MT, DE, and absorbing power of low-voltage distribution grids. Since PV and WT, as new energy units, do not produce pollution gas during operation, the pollution degree is not considered. e object function formula is as follows: where F 2 is the treatment cost of pollutant discharge; M is the type of pollutants discharged such as CO 2 , SO 2 , and NO x ; C j is the cost of discharged pollutants; K ij is the discharge coefficient of pollutant j in the i-th power generation unit; and P i,t is the generation capacity of the i-th unit at time t.

Constraints
(1) e power balance constraint is as follows: where P i,t is the power output of each micropower source at time t; P ES,t is the output power of energy storage battery; when P ES,t > 0, it means discharging; P L,t is the microgrid load at time t.

Mathematical Problems in Engineering
(2) e output power constraint of micropower supply is as follows: where P imin and P imax are the minimum and maximum output power, respectively. (3) e climbing speed constraint is as follows.
When the load increases, When the load is reduced, Here,d R iup and R idown are the upper and lower limits of climbing rate of the i-th power generation unit. (4) e charging and discharging constraints of energy storage battery are as follows.
To ensure the battery life, charging and discharging must follow the following formula: where SOC imin and SOC imax are the minimum and maximum capacity of energy storage battery.

Constrained Multiobjective Optimization Problem.
A constrained multiobjective optimization problem (CMOP) with minimization can be defined as defines the function which is mapped from the decision space to the objective space; g(x) and h(x) define q inequality constraints and m − q equality constraints.

Modified Multiobjective Biogeography-Based Optimization Model and Algorithm.
In view of the fact that the BBO algorithm itself does not have the ability to deal with multiobjective optimization and constrained multiobjective optimization problems, a multiobjective optimization model is proposed in this paper, which is suitable for the BBO algorithm based on the new multiobjective constraint handling mechanism. In addition, considering the greater effect of the performance of the evolutionary strategy on the overall performance of the algorithm, in order to improve the constrained optimization ability of multiobjective optimization, the efficiency of BBO algorithm must be ensured.

Constrained Multiobjective Handling
Mechanism. e main idea of ε constraint handling technology is to measure the degree of the individual by comparing the level parameter ε with the individual constraint violation degree, enlarge the constraint region effectively under the allowable constraint violation degree, and promote more excellent and infeasible individuals with lower constraint violation degree to participate in the evolution, so as to guide the whole group to approach the feasible region. For the individual en, the comprehensive constraint violation degree of the individual ϕ(x) is defined as the linear addition of each constraint violation degree.
However, the ε constraint handling method based on linear addition and simple constraint violation has its drawback in solving CMOPs that ignore the difference of the degree of difficulty satisfying the different constraints. In order to reflect the real constraint violation degree more objectively, a new self-adaptive weighted ε comprehensive constraint violation degree determination method is proposed as shown in where NP is the total number of individuals in the current population; n j fea is the number of individuals satisfying the jth constraint (φ j (x) � 0); ε j (t) is the level parameter of the jth constraint at t time, as determined by formula (5) in the following chapter. In order to further measure the quality of individuals, we also need to integrate the multiobjective Pareto dominance into the comparison among individuals according to the characteristics of CMOPs. erefore, for the individual x 1 and x 2 in population H, the ε constraint dominance relationship combining the Pareto dominance with the adaptive constraint violation degree ϕ(x) is shown in In the new ε constraint multiobjective handling mechanism proposed above, each constraint condition j needs to be given a level parameter epsilon ε j . In view of the above problems of the original level parameter ε which is based on the nonlinear decrease of iterations, a new method to determine ε is shown in where φ t j (x)is the constraint violation of individual X i at time t; φ t jmax and φ t jmin are the maximum and minimum of φ t j (x). For the jth constraint, all of ε j of each generation is adaptively adjusted for the changes of feasible region according to the jth constraint violation degree of the individual x. φ j (x) participates in regulating the level parameter in the form of its proportion in the violation degree of individuals. e larger φ j (x) may affect the level parameter ε j with a smaller proportion, and vice versa.

Modified Determination Method of Habitat Suitability
Index. In the BBO algorithm, habitat suitability index (HSI) is an indicator of measuring the habitat quality and also an important factor affecting the distribution and migration of habitat populations. For CMOPs, in consideration of the characteristics of multiobjectives, it is necessary to determine HSI by combining Pareto dominance among habitat individuals. However, it is far from enough to define the HSI which is not enough to evaluate the quality of the habitat in this way and the constraints should also be taken into account fully.
us, it is ideal to design a determination mechanism of HSI, which can take into account the Pareto domination relationship and the degree of constraint satisfaction of the individual itself.
For this purpose, this paper proposes a new determination method of habitat suitability index HSI on the basis of the idea of ε constraint domination mentioned above. First of all, the individuals in habitat population H � (x i , i � 1, 2, . . . , NP) are ranked in ascending order in accordance with the adaptive ε weighting constraint violation and form a new habitat population H.
where x i , x j , x k ∈ H is the habitat individual including ndimensional suitability index variables (SIV). Finally, the HSI G i of the individual x i in the new H is determined by

Determination Method of Migration
Rate. e determination mechanism for the migration rate of the original BBO does not apply to CMOPs, which can only reflect the order of individuals but ignores the current evolutionary state.
us, it is impossible to effectively evaluate the quality of individuals. When the gap between individuals is large or the distribution is unbalanced, it is very easy to make a wrong assessment of the quality of individuals, resulting in the fact that the information of better individuals cannot be retained or worse individuals participate in evolution too early.
In order to reflect the quality of individuals in real time and ensure that the information of the better individuals can be made full use of in different evolutionary stages, a method of determination of immigration rate λ and the emigration rate μ is shown in the following formula: where G i is the HSI of x i ; α is the nondominated rate of H εfea , the proportion of nondominated ε feasible individuals in the number n εfea of H εfea ; G max and G min are the maximum and minimum HSI of G i ; η is the ε feasible rate, the proportion of n εfea in NP. e HSI of normalized individuals is used to adjust the immigration rate in combination with the proportion of ε feasible individuals and the proportion of ε feasible nondominated individuals in the current population. According to the difference of the evolutionary stages, the role of adaptive migration rate determination Mathematical Problems in Engineering 5 mechanism in the evolution of populations can be attributed to the following three scenarios: (a) At the beginning of population evolution, when the number of ε feasible individuals is small, the adaptive migration rate can give more opportunities to ε feasible individuals or ε infeasible individuals with smaller constraint violation to guide the evolution of the population according to the quality of individuals and rapidly evolve the individuals with larger constraint violation to approach the feasible region from different directions (b) As the number of ε feasible individuals gradually increases, the adaptive migration rate can guide the ε feasible individuals to move towards the direction with better nondomination degree and guide the remaining ε infeasible individuals to approach the feasible region or even enter the feasible region, so as to efficiently evolve the less nondominant individuals and the remaining individuals with ε infeasible (c) When the population is composed of feasible individuals and the number of nondominant feasible individuals gradually increases or even spreads over the population, the more concentrated nondominant feasible individuals have more opportunities to participate in evolution and thus become more evenly distributed

Modified Migration Strategy.
For CMOPs, when the search space is nonconvex or disconnected, the algorithm often falls into a certain subregion and the solution found is a local Pareto optimal solution, which makes it very difficult for the algorithm to solve CMOPs and poses higher requirements for the evolutionary mechanism of the algorithm. In order to overcome the effects of the constraints, the immigration rate with ε feasible and Pareto dominant information is used as a control parameter to adjust the current habitat individual x i and the emigrated habitat individual x k adaptively. In the meantime, the difference vector of exploiting the new solution is added and the disturbance of which is dynamically adjusted by the time factor. e modified dynamic migration strategy is shown in where x ij is the jth-SIV of x i ; r is the time factor that can control the amplitude of the difference vector perturbation; r max and r min are the maximum and the minimum of r; r 1 and r 2 are the random number between 1 and NP; and β is the specified lower normal number. e proposed dynamic migration operator can dynamically adjust the individual's evolutionary trend at different evolutionary periods. At the beginning of evolution, the first half of formula (23) is used to rapidly guide the evolution of the population towards the better direction, and the second half is used to perturb the population by the small amount of difference vector. In this way, when the number of ε feasible individuals is small, the population can fully absorb the information of the better ε infeasible individuals or even ε feasible individuals, which ensures the rapid population evolution towards the ε feasible region. With the continuous evolution, as the number of ε feasible individuals in the population increases and the presence of nondominant ε feasible individuals gradually increases with evolution, the population can continue to converge towards the nondominant frontier of the feasible region at high speed.
However, when the population has evolved to a certain stage where the population is almost entirely composed of nondominant feasible individuals, in order to make the nondominant feasible individuals closer to the Pareto frontier and obtain a better distribution, at this time, the algorithm mainly uses the time factor to dynamically adjust the difference vector in the second half of formula (23) to realize the full exploitation of new individuals. e gradual increase in the magnitude of the deviation perturbation can accelerate the spread of individuals to sparse regions, improve the population diversity, and provide more opportunities for escaping from local extremes.

e Improved Algorithm Flow.
According to the analysis of multiobjective optimization model above, the optimal scheduling of a microgrid is a complex nonlinear optimization problem with multiobjective and multiple constraints.
e specific procedure of dynamic optimal scheduling of microgrid based on improved constraint multiobjective biogeography-based optimization model and algorithm with adaptive migration (MOBBO-AM) is shown in Figure 1.
Step 1. Set the relevant data and operating parameters of the microgrid, which involves the thermoelectric load, the power generation of wind turbines and photovoltaic cells, and so on. And set the relevant constraints of each micropower source.
Step 2. Initialize the relevant parameters of the proposed algorithm, which involves m max , r max , r min , and β. Randomly generate NP individuals as the initial population H � (x i , i � 1, 2, . . . , NP). For each H, it is generated by a micropower in different periods and its dimensions are N × T. N is the number of distributed generation units to be optimized in the microgrid, and T is the total scheduling period. Step Step 4. Calculate the ε value under each constraint of x i , and determine the ε constraint violation degree ϕ(x i ) of x i . Sort x i by its ϕ(x i ) in ascending order, and reform a new population H. According to ϕ(x i ), determine the feasible population H εfea in H, calculate the nondominated degree F i 6 Mathematical Problems in Engineering of individual in H εfea , and thus obtain the HSI value G i of each habitat individual x i .
Step 5. Sort the habitat individuals of H by HSI in ascending order, and reserve the NP individuals with better distribution and nondominated as current optimal habitat population H.
Step 6. Determine the species immigration rate λ i and the emigration rate μ i of x i in H, and calculate the species probability P i and mutation rate m i . On the basis of immigration rate λ i and the emigration rate μ i , make migration and mutation, and thus obtain new population H. Judge whether it satisfies the termination condition or not. If it satisfies, then go to Step 7; otherwise, go to Step 3.
Step 7. Output the nondominated feasible population H.

Basic Data of Calculation Example and Parameter Setting.
In order to verify the effectiveness of the above-mentioned model, a typical microgrid system, including photovoltaic cells, wind turbines, fuel cells, microgas turbines, diesel generator sets, energy storage battery, and load, carried out the optimal scheduling. In this paper, the output of each distributed generation (DG) unit in 24 periods of a day is taken as the optimization object, and its parameters are shown in Table 1. e typical daily load data is selected, and the prediction curves of load, temperature, wind speed, and light intensity are shown in Figure 2. Meanwhile, on the premise of maximizing the use of renewable energy, PV and WT work in maximum power point trace mode (MPPT), and the output curve of WT and PV is shown in Figure 3. e discharge data, pollutant value standard, and penalty level of pollutants are shown in [16]. e relation curves of the environmental treatment cost and the output power, the fuel cost, and the output power of controllable microsource are shown in Figures 4 and 5, respectively. e parameters of MOBBO-AM algorithm proposed in this paper are set as follows: the initial population size NP is set as 100, m max � 0.005, r max � 0.9, r min � 0.4, and β � 0.02, and the maximum number of iterations is 500.
From Figures 4 and 5, it can be seen that FC has the lowest environmental pollution treatment cost during operation, followed by MT, and DE has a higher environmental cost. For the operating cost, the fuel cost of MT and FC is higher in operation, in which MT is slightly higher than FC, while the fuel cost of DE is lower.

Simulation Results and Analysis.
In order to realize the comprehensive optimal scheduling of economy and environmental protection of microgrid system, the above microoptimal scheduling model combined with improved algorithm is used for experiments, and the optimal Pareto frontier of operating cost and environment cost is shown in Figure 6.
As can be seen from Figure 6, the environmental cost range is about $40 to $140, and the operating cost range is about $270 to $470 on the Pareto frontier. When the operating cost is low, the environmental cost may rise, and vice versa. If the microgrid operates at the lowest operating cost, the total cost of simple summation is about $410, while the total cost of operating microgrid at the lowest environmental cost is about $510. erefore, taking into account the environmental cost and operating cost, it is likely to obtain the optimal solution with lower comprehensive operation cost or equivalent comprehensive cost but less environmental cost by choosing the compromise solution on Pareto frontier. In this example, when the operating cost is $290, the environmental cost is only $83, for a total cost of $373, which not only reduces the environmental impact but also maximizes the economy. It is not difficult to see from the figure that the relationship between operating cost and environmental cost is contradictory, and with the increase of operating cost of microgrid, the pollutant treatment cost of microgrid is relatively reduced. It is necessary to comprehensively consider various factors to choose a compromise between them. e algorithm in this paper can not only obtain the effective convergence solution set but also obtain the Pareto frontier evenly and widely distributed, which provides a more Judge whether it satisfies the termination condition or not.
Output the nondominated feasible population H.
End Y N Sort the habitats by HSI in ascending order, and reserve the individuals with better distribution and nondominated as current optimal habitat population.
Calculate the feasible population and the nondominated degree of individual to obtain the HSI value of each habitat. Mathematical Problems in Engineering comprehensive and reasonable choice for the economic and environmental operation and dispatching of microgrid. ree representative solutions are chosen on the Pareto frontier in Figure 6 and their corresponding microgrid scheduling scenarios are analyzed briefly, where Figures 7 to 9 are, respectively, the cases with lower operating cost, lower environmental cost, and minimum comprehensive cost. Afterward, the economic scheduling with the minimum operating cost, the environmental scheduling with the minimum environmental cost, and the comprehensive optimal scheduling are analyzed, respectively, and the results are shown in Table 2.
It is not difficult to see from Figure 7 that, considering the economic cost and environmental protection cost, the operation of microgrid needs to be analyzed in different periods. At 1:00 to 7:00 of electricity valley time, there is no light but there is wind at this moment; PV does not generate electricity, but WT can supply power continuously. Considering the supply and demand of electricity, FC and MT keep the low power operation, DE replenishes the remaining  electricity demand, and the battery remains charged. After 8: 00, electricity consumption continues to rise; although the light intensity gradually increases, PV and WT still cannot meet the demand for electricity. When DE power gradually increases, the battery starts to discharge, which is used to meet the load demand for electricity. After 11:00, when the power consumption continues to rise to the point where DE is close to full load, FC begins to replenish the power demand while MT still operates at a lower power.
At 15:00 to 17:00, there is a local electricity valley time. At this time, DE is still running at higher power, while the power of FC is reduced. After the battery is discharged at the peak of power consumption in the morning, it is charged at this stage to store the energy for the late peak. After 18:00, electricity consumption begins to rise again and reaches the    Mathematical Problems in Engineering peak time in the whole day at 20:00, when there is no light and no PV output, and WT remains continuously powered. At this time, FC and DE are basically close to power operation in full load. e battery begins to discharge, while MT also operates with a higher power to meet the demand for electricity load. After 23:00, the electricity consumption decreases, the remaining capacity of the battery reaches the lower limit, and the battery starts to be charged. e operating scenario in Figure 7 is the case of minimized operating cost. It can be seen that when cooperating with the PV and WT output, DE with the lowest operating cost is preferred, followed by FC, and finally MT is used as the compensation power demand in peak period. is scenario makes extensive use of DE, so its environmental cost is higher.
As can be seen from Figure 8, this operating scenario is to minimize the environmental cost. Unlike Figure 6, when FC can meet the power demand, MT and DE operate at lower power. And when FC power is insufficient, MT is called upon first to supplement power generation, and DE participates in power generation only during the peak period of the whole day to ensure the electricity demand. In this scenario, the DE with the highest environmental cost is used with the least possibility, but FC and MT are kept running at higher power, so the operating cost is higher. e operating scenario in Figure 9 is the case of optimal comprehensive cost, that is, the compromise solution in Pareto front. It can be seen that when the demand for electricity is low, the output of using FC to cooperate with PV and WT is the main output, while when the demand for electricity is higher, the output of FC is the main output, followed by DE, and MT is used as a deficiency supplement.
is scenario takes into account the low operating cost of DE and the low environmental cost of FC and MT, which is the comprehensive optimal scenario combining the operating cost and environmental cost.
It can be seen from Table 2 that, when the operating cost is minimum, the power generation cost of microgrid increases because of increasing the input cost of diesel generator sets with high pollutant emissions, which leads to the increase of environmental cost. When the environmental cost is minimum, the cost of pollutant treatment leads to the relative increase in operating cost. When the comprehensive cost objective is minimum, its environment cost and operating cost are in a limited range, and then it is necessary to balance various factors such as the distribution of wind and light resources, fuel prices, and pollutant treatment costs and choose a reasonable dispatching scheme when optimizing scheduling. It is likely to obtain the optimal solution with lower comprehensive operation cost or equivalent comprehensive cost but less environmental cost by selecting the compromise solution on Pareto frontier, which not only reduces the environmental impact but also maximizes the economy.
en, in order to further verify the performance of the proposed MOBBO-AM algorithm, it is compared with two typical constrained multiobjective optimization algorithms. And for the purpose of proving the effectiveness and advancement of the proposed algorithm, we first compare the algorithm on a standard test set, and its performance is compared with that of NSGA-II [24] and NSDE [28] on the representative standard test problems [31] of TNK, SRN, CONSTR, and OSY. e convergence and distribution performance of the algorithms are tested using two evaluation criteria [32]: generational distance metric (GD) and    to quantitatively measure the convergence performance of the solution set. e statistical box plots of MOBBO-AM and the other two algorithms are shown in Figure 10, where MOBBO-AM, NSGA-II, and NSDE are represented by numbers 1 to 3, respectively. It can be seen from Table 3 and Figure 10 that, for the 4 test functions, MOBBO-AM can not only obtain the minimum average GD value compared to NSGA-II and NSDE but also obtain the minimum GD value which is significantly smaller than that obtained by other algorithms. is shows that the results obtained by MOBBO-AM are closer to the ideal Pareto frontier and has an obvious advantage in convergence. Meanwhile, Figure 10 also shows that the statistical data obtained by MOBBO-AM also has fewer outliers than the other algorithms, which indicates that MOBBO-AM has the best convergence stability among the 3 algorithms.
For the distribution of the solution sets obtained by the three algorithms, the experimental statistical results of SP index are given in Table 4, and the statistical boxplots are shown in Figure 11.
From Table 4 and Figure 11, it is not difficult to see that MOBBO-AM obtains the smallest SP value and the least outliers in the statistical data for the 4 test functions, which indicates that MOBBO-AM has more advantages than the other two algorithms in terms of the uniformity of distribution and can make the distribution of approximate Pareto optimal feasible solution set more uniform.
Summarizing the above experiments on the standard test functions, MOBBO-AM outperforms the other 2 algorithms in terms of both convergence and uniformity of distribution. Next, we use the algorithm for solving the optimized scheduling model for microgrids. e MOBBO-AM algorithm proposed in this paper is shown in Figure 6 and the simulation results of the NSGA-II and NSDE algorithms involved in the comparison are shown in Figures 12 and 13.
As can be seen from the figure, the optimization result of the proposed algorithm is that the solution sets are very evenly distributed on Pareto optimal frontier, while the Pareto optimal frontier of microgrid optimal scheduling  based on NSGA-II and NSDE makes the assemblage and discontinuity of solution set in different degrees, and the distribution uniformity is obviously inferior to the proposed MOBBO-AM. erefore, compared with the comparative algorithm, the proposed algorithm can more effectively and comprehensively reflect the contradiction between the environmental cost and the operation cost of microgrid. As the basis for guiding the multiobjective optimization operation and research of microgrid, the optimization quality of the Pareto optimal solution set directly determines the selection of operation scheme. Hence, obtaining the Pareto optimal frontier with excellent uniformity and convergence performance can provide a more complete selection for the operation scheme, without missing the better compromise solution which is more in line with the actual demand owing to the discontinuity of the solution set.

Conclusion
In this paper, the dynamic optimal scheduling problem of microgrid in island operation is studied, and a multiobjective optimal scheduling model considering operating cost and environmental treatment cost is established. Taking   Operating cost ($) Figure 13: Pareto optimal front obtained by NSDE.
Mathematical Problems in Engineering 13 a specific microgrid system including wind, light, storage battery, microgas turbine, fuel cell, diesel engine, and load as an example, constraint multiobjectives and biogeographybased optimization algorithm are organically integrated, and the corresponding constraint handling mechanism and the determination method of habitat suitability index with adaptive migration strategy are effectively improved. e simulation results show that the optimal scheduling algorithm based on multiobjective biogeography-based optimization with adaptive migration proposed in this paper is significantly superior and it can simply and quickly achieve global constrained multiobjective optimization to obtain accurate and complete Pareto frontier. According to biobjective operation results of operating cost and environmental cost, it can effectively better guide the optimal operation of microgrid according to the actual situation and play a certain guiding role in ensuring the islanded operation of microgrid and uninterrupted power supply for sensitive loads in microgrid, which can provide a reference for the economic and environmental protection scheduling of energy conservation and emission reduction.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest.