#### Abstract

With the increase in the proportion of clean energy connected to the grid, the effective coordination of the operation of various energy power has become a new challenge facing the current power system scheduling. The coordinated operation of the clean energy power generation system can alleviate the contradiction between power generation and output power fluctuations and overcome the bottleneck of new energy development. Considering the natural characteristics of clean energy, this paper aims to make full use of clean energy, reduce system operating costs, increase system power generation, and reduce output fluctuations; we establish a multiobjective optimization model for coordinated scheduling of clean energy power systems. The model seeks to maximize power generation and minimize output fluctuations, power purchase costs, and maintenance costs under the constraints of the grid structure. In this paper, the joint algorithm has an accelerated effect on the target optimization calculation, and then the superiority of the algorithm is verified by the standard system. The standard IEEE39 node test system is used to verify the rationality and feasibility of the model built and provides a reference strategy for the coordinated operation mechanism of the clean energy system. According to the model, in the example in this paper, the maximum value of photovoltaic power prediction is 1290 MW, and the minimum value is 210 MW; the maximum value of wind power prediction is 780 MW, and the minimum value is 28 MW; the minimum cost of power purchase and maintenance is 56,950.395; the maximum generating capacity is 5.045 GW; the minimum output fluctuation is 0.120 GW.

#### 1. Introduction

Large-scale wind and photovoltaic grid integration is of great significance to power energy saving and emission reduction, but the uncertainty and limited storability of wind-solar-water energy power generation have led to the complexity of power system scheduling and output fluctuations and power reduction related issues [1], which severely hinder the use of renewable energy and slow down the speed of the reduction of carbon emissions. The coordinated operation of hydropower and intermittent renewable energy can effectively solve these problems. Through modeling and optimizing the coordinated operation of the wind-solar-water system [2], a multisource coordination mechanism was discussed and studied to improve the utilization rate of clean energy and achieve a large-scale energy optimization configuration.

Currently, the research on wind-photoelectric fields in China and abroad is mainly based on economic dispatch [3, 4]. In recent years, attention has been focused on the theoretical characterization of the stochastic characteristics of wind power. Focusing on water-wind, water-photovoltaic, or water-wind-photovoltaic system, the coordinated operation was modeled and optimized [5]. The goals of the current model usually include power generation, economic benefits, emission reductions, water supply, power generation costs, power reduction, production fluctuations, and flow patterns of hydropower generation. At present, the goal of the model is usually to reflect the effect of power system coordinated operation in many aspects, such as power generation, economic benefits, emission reduction, water supply, power generation cost, power reduction, output fluctuation, and flow pattern of hydropower flow [6, 7]. Since the cost of hydropower generation system is low and almost constant in the short term, most of the optimization objectives of the model are focused on economic dispatch [8]. The main purpose of multienergy coordination is to ensure the stability of the power grid system under the premise of regulating the generation of various energy sources and reducing output fluctuations and to provide a theoretical basis for the construction and development of integrated energy systems. With the rapid development of Chinese wind-solar-water-heat clean energy and the lag of power grid construction, power grid constraints, including the power framework and power transmission capacity, have become the main constraints for coordinated operation.

Most of the earlier studies focused on increasing the profit of power generation and reducing output fluctuations but did not consider the relationship between power generation and output fluctuations [9]. Actually, due to the randomness of wind and solar radiation, generating more power led to increasing output fluctuations. In order to meet the needs of users for various energy sources, in literature [10], the authors established a joint scheduling model of regional integrated energy system ( for short) with renewable energy and combined cooling heating and power ( for short) and energy storage system, which mainly focused on the energy conversion efficiency and economic benefits of energy conversion equipment and energy storage equipment, but they did not consider the impact of the volatility of different output objects on the whole energy system. In [11], taking the constraints of multiple types of power generation and energy storage into consideration, a model with renewable energy and cooling, heating, and power was established to minimize the overall operating cost, but the randomness and output fluctuation of clean energy were not considered, which increased the difficulty of scheduling, and seriously affected the balance between power generation and output fluctuation. The coordinated operation of clean energy followed the correlation to effectively balance the relationship between multiple targets and the degree of dependence [12], but the study did not take into account different resource changes in different seasons. Moreover, when the transmission capacity of the power grid reached the limit, there was no reduction in power distribution among the various power stations, which was unavoidable for the lag of power grid construction; at the same time, the network architecture of the power system studied was mostly single-layer, which could not reflect the embedding structures of the power grid.

Based on the above literatures, the current research objectives are mainly concentrated in a single dimension. In the current situation of large-scale integrated energy access to the grid, this paper will consider multidimensional research from the energy system to further improve the whole energy system architecture, with multiple objectives of energy system stability, economic benefits, and system power generation as the research object. This paper will mainly study from the following three dimensions: (1) In the case of a variety of energy into the grid, we further discuss the relationship between the energy generation and output fluctuation in order to minimize the output fluctuation. (2) To optimize the economic benefits of the whole system, the main objective is to minimize the power purchase cost and equipment maintenance cost of the main power grid so as to improve the economic benefits of system. (3) This paper explores the detailed coordinated operation mode of the multienergy system, comprehensively considers the typical resource conditions in different seasons, extracts the multilayer power grid structure, and establishes the multiobjective model of the coordinated operation of the integrated energy system. In view of the above problems, considering the advantages of large-scale integrated energy access to the grid, IEEE39 will be adopted to test the rationality and feasibility of the model and algorithm. In this paper, we hold a view to providing data support and application guidance for the coordinated operation of the future energy interconnection system.

The rest of this paper is arranged as follows: in Section 2, the description and establishment of the basic principles of the model are involved; in Section 3, the solution steps of the model are demonstrated in detail; in Section 4, taking the IEEE39 node test system as an example, we comprehensively discuss the research results of the energy system coordination mechanism; finally, we come to the conclusions in this paper.

#### 2. Multiobjective Model of the Clean Energy System

The goal of the water-wind-solar clean energy system coordination model is to maximize the generation of clean energy systems and minimize output fluctuations to maximize economic benefits as well [13, 14]. Constraint types include resource types (water flow, wind speed, and solar radiation), grid types (power balance, transmission capacity, and power reduction allocation), and power plant types (water level, flow, capacity, and energy utilization factor).

##### 2.1. Objective Function

In this section, we study the objective functions that are crucial for optimal control from three aspects, respectively, that is, power generation, power generation cost, and power system volatility.

First of all, we put forward objective function as shown in equation (1) corresponding to mathematical model (2); the goal in the objective function is to maximize the power generation of the clean energy power generation system. In order to overcome the random characteristics of wind-solar energy, we actively coordinate the power generation of the entire system by reducing the discharge water flow and using the water saved in the reservoir as the residual income of the current optimization range. Because the saved water can be used for power generation in the future and is the result of coordination within the current range, it is converted into power generation and included in the total power generation.satisfyingwhere represents the number of periods: , and is the biggest number of periods; represents the time length of a single period; denotes total output of hydropower, wind power, and photovoltaic power generation in the *i*-th period; denotes hydropower in the *i*-th period, is short for hydropower; denotes wind power in the *i*-th period, is short for wind power; denotes photovoltaic power in the *i*-th period, is short for photovoltaic power; denotes the remaining water volume in the reservoir; is conversion coefficient; denotes water density; is the amount of water remaining in the reservoir at the end of the optimization; denotes the acceleration of gravity; denotes the average head of remaining water.

Furthermore, in order to realize the operation of the system economically, efficiently, and environmentally friendly, we can follow the policy of giving priority to wind power and full access to the grid in accordance with the requirements of the Renewable Energy Law. Therefore, we set the main power grid purchase cost and maintenance cost as the second optimization goals as follows [15]:satisfyingwhere represents the cost of purchasing electricity and represents the cost of equipment maintenance; , , and denote the power purchased, the price of electricity purchased, and the price of electricity sold, respectively, during the period; denotes the unit maintenance cost of device ; represents the output of device in the period.

Finally, in the current energy system, it is not enough to consider only the output fluctuations of clean energy, since this will bring about huge errors and irrationality to the subsequent research. Therefore, we should take the output fluctuations of thermal power into consideration to derive the minimum output fluctuations of the clean energy system. Thus, we further put forward the objective function as follows:satisfyingwhere denotes quantitative fluctuations and contour fluctuations of , respectively, is the mean of , the definitions of parameters such as can be found in [12], is the minimum output fluctuation of the thermal power unit in the adjacent period, is the output power of the thermal power unit in the period, and represents the total number of thermal power units.

##### 2.2. Constraints

As mentioned before, constraints are divided into three types [16]. Large and complex grids are composed of simple regional grids with hierarchical interconnections, resulting in coordination between water-wind-solar energy sources that may occur at different grid levels. Therefore, there are mainly the following aspects when considering the constraints of multilayer power balance.(i)Power balance constraints: where represents the load of the layer in the period; is the power exchange between layer and layer in the period. While , it indicates downloading power from the layer; while , it indicates uploading power from the layer; while , it indicates no power exchange. denote the output of electricity from a hydraulic, wind, photovoltaic, thermal power generation, and other power sources at a node of the *m*-th layer in the period, respectively.(ii)System backup constraints: where is the on-off state of the thermal power unit in the period, 1 stands for booting and 0 stands for out of service; represents the maximum output power of thermal power unit ; represents the maximum output power of hydropower power unit ; denote the total numbers of hydropower units and wind power units, respectively. represent load fluctuation coefficient and wind power fluctuation coefficient, respectively, and we value them as and . The reserve capacity is required to meet the random fluctuations of wind power and load.(iii)The upper and lower limits of thermal power output: where and represent the minimum and maximum output power of thermal power unit , respectively.(iv)Hydropower generation capacity and conversion efficiency constraints: where and represent the minimum and maximum output power of the hydropower unit , respectively; represents hydroelectric conversion constant; denotes the water consumption of reservoir in the period; represents the head height of the reservoir in the period.(v)Wind power capacity constraints: where and denote the upper and lower limits of the wind turbine output , respectively; and denote the cut-in and cut-out wind speeds of wind turbine generator, respectively; describes the relationship between wind energy output and wind speed.(vi)Photovoltaic power generation capacity constraints:where and are the upper and lower limits for photovoltaic output, respectively; represent photovoltaic array coefficient, inverter coefficient, and connection power coefficient, respectively; denotes solar radiation intensity.

##### 2.3. Research Objectives

In order to discuss the mechanism and superiority of coordinated operation of multitarget clean energy power generation system, the analysis based on the above model was studied. In this model, the goal is to maximize the power generation and minimize output fluctuations through the coordinated operation of clean energy power generation systems, while minimizing the cost of power purchase and maintenance. The objective functions are as follows:

#### 3. Model Solving Algorithm

The clean energy power generation system is a large-scale system with multiobjective, nonlinear, and strong coupling characteristics. The coordinated operation of the system includes decoupling the network architecture from the lower layer to the upper layer and independent configuration according to the characteristics of the power station, which leads to the problem of multiobjective collaborative optimization which is difficult to converge. Using the elite strategy of simultaneous cloning and evolution and the idea of similar mating [17], the genetic algorithm () and particle swarm optimization () algorithm are combined to further improve its convergence. Due to the low complexity, fast convergence, and uniform distribution of the Pareto boundary, this paper applies the joint algorithm to solve multiobjective optimization problems. In this study, the MATLAB multiobjective optimization toolbox was used to determine the best coordinated operation of multiobjective optimization.

joint algorithm implementation steps are as follows: *Step 1 (Initialization).* For an -dimensional problem, a group with a size of is generated within the feasible domain. In this paper, the clean energy output ratio is used as the initial group. When the grid automatically stratified number , some points are selected from the random samples as the initial point, and the individuals constitute the evolutionary operation. *Step 2 (Sorting)*. Calculate the fitness values of individuals and rank them from good to bad. *Step 3 (GA Method)*. The individuals above are cross-mutated with the real-coded GA operator to generate new individuals. *Step 4.* Calculate the fitness value of the new chromosomes generated by GA and replace the worst individual with the elite parent individual; that is, the fitness value is inferior to that of the parent elite individual, and place it in the corresponding position. *Step 5*. The above chromosomes are sorted according to the fitness value and paired with each other. Each team forms a neighborhood and determines the optimal individual. *Step 6 (PSO Method).* The best fitness individual in GA algorithm is regarded as the global extremum of the PSO algorithm, and the better fitness value in each neighborhood is regarded as the individual extremum of the neighborhood. Compared with the previous neighborhood optimal particle, the optimal particle is updated in real time. *Step 7*. Determine whether the convergence conditions are met. If the convergence condition is satisfied or the number of iterations is reached, the algorithm ends; otherwise, return to step 2 and repeat the above process until the convergence condition is satisfied.

The logic diagram of the joint algorithm is shown in Figure 1.

#### 4. Example Analysis

We adopt the IEEE14 node test system appearing in [18] as an example and consider the lowest purchase cost and maintenance cost as the goal. Suppose that the population size is 14, the maximum number of iterations is 100, and the calculation time is 25. GA, PSO, and GA_PSO algorithms are used, respectively, for optimization and comparison.

It can be seen from Figure 2 that the algorithm has much better robustness, faster convergence, and more reasonable calculation results; calculation functions often have premature phenomena, and the system convergence is rather weak; optimization algorithm has the best iterative convergence, and the system tends to be more stable, so the optimization algorithm has an accelerated effect on the optimization calculation of the function. From the convergence of the above three algorithms, it can be seen that when the lowest power purchase cost and maintenance cost are selected as the objective functions, the proposed algorithm is the most feasible to optimize this model.

#### 5. Applications on Models and Algorithms

We adopt IEEE10 machine 39-node standard system model [19] and standard data of the system and use part of the existing data [20] to solve the objective functions proposed in this paper.

Figure 3 shows the predicted power value and the load value in each cycle of the wind farm.

The predicted maximum value of wind power is 780 MW, accounting for 78% of the total installed capacity of the wind farm; the minimum value is 28 MW, accounting for 2.4% of the total installed capacity of the wind farm. It can be seen that the peak-valley difference of the random fluctuation of wind power output is very large, and it can also be seen from Figure 1 that wind power has a significant inverse peak-shaving characteristic relative to the load.

Through the statistics of photovoltaic output, it can be seen from Figure 4 that the photovoltaic power distribution is close to a normal distribution, and, at the same time, the photovoltaic output has strong trend characteristics and obvious random fluctuations. The above analysis provides basic data and reliable theoretical basis. The predicted maximum value of photovoltaic power is 1290 MW, and the minimum value is 210 MW.

By comparing the data in Table 1, it can be seen that the three optimization algorithms have obtained better results, which fully shows that the joint algorithm increases the particle diversity and effectively improves the defect that the algorithm is easy to fall into the local optimum. Among them, the final optimization result of power purchase cost and maintenance cost based on joint algorithm is dollars.

As seen from Figure 5, in view of the surplus water resources in the flood season, the dispatch center considers that, from the safety point of view, the hydropower plant must be full, so the hydropower in the flood season has no regulation capacity. Clean photovoltaic and wind power are given priority, so the optimized variable wind-light output has the smallest fluctuations and ensures the maximum power generation. The output of photovoltaics fluctuates significantly, with a large peak-to-valley difference of , which accounts for about of the maximum output. This shows that, in the case of photovoltaics that do not have the ability to regulate, wind power must undertake the dual task of making up for the random fluctuations of photovoltaics and the difference between peaks and valleys in load, which leads to frequent adjustments of wind power output.

In Table 2, the symbol “COV” means “cost optimization value”; the symbol “OVPG” means “optimal value of power generation”; the symbol “OVOF” means “optimal value of output fluctuation.”

For the multiobjective optimization of the clean energy system, the relationship between the multiobjective optimizations can be obtained through the literature [21] to satisfy the Pareto relationship. Therefore, this paper uses the joint algorithm to optimize power purchase costs and maintenance costs, power generation, and output fluctuations for multiobjective optimization. From Table 2, we can see that the minimum power purchase cost and maintenance cost is ; the maximum power generation is ; the output fluctuation is minimized to .

#### 6. Conclusions

In view of the limitations of the traditional power system clean energy joint scheduling in large-scale clean energy centralized access system and the contradiction between multienergy power generation and output fluctuations, in this paper, we establish a clean energy power generation system model in order to study the coordinated operation mechanism of different energy. Our objectives are to make full use of clean energy, reduce system operating costs, and make use of mutual compensation between different clean energy sources to to avoid output fluctuations and power reduction. The simulation results of the example show the following:(1)The algorithm uses selection, crossover, and mutation operator operations for function optimization and directly uses the target function of the model as search information in a probabilistic manner, thus enhancing the global optimization ability of the optimization algorithm, speeding up the evolution of the algorithm, and improving the convergence accuracy.(2)The model built in this paper can provide some data supporting for the scheduling strategy of the clean energy system and a strong guarantee for the full coordination of clean energy power generation systems, which can further improve the economic benefits of the power grid and ease the pressure of peak and frequency regulation of the power grid and ensure the maximum utilization of the energy system as well.(3)Adopting the coordinated operation mechanism of the clean energy power generation system can improve the grid’s ability to accept clean energy under the premise of ensuring the operation efficiency and stability of the power grid system and can provide reference strategies for the coordinated operation of other power sources. Additionally, coordinated operation mechanism of new energy power generation also can provide data support and application guidance for the construction and coordinated operation of the energy Internet.

#### Data Availability

Part of the data is acquired from other published papers (cited), and some data are estimated in a reasonable range.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was supported by the National Natural Science Foundation of China (NSFC) under Grant no. 11661050.