# An Optimal Allocation Strategy for Multienergy Networks Based on Double-Layer Nondominated Sorting Genetic Algorithms

Authors: Min Mou; Da Lin; Yuhao Zhou; Wenguang Zheng; Jiongming Ruan; Dongdong Ke
Journal: Complexity
 (2019)

---

## Abstract

Aiming at the problems of complex structures, variable loads, and fluctuation of power outputs of multienergy networks, this paper proposes an optimal allocation strategy of multienergy networks based on the double-layer nondominated sorting genetic algorithm, which can optimize the allocation of distributed generation (DG) and then improve the system economy. In this strategy, the multiobjective nondominated sorting genetic algorithm is adopted in both layers, and the second-layer optimization is based on the optimization result of the first layer. The first layer is based on the structure and load of the multienergy network. With the purpose of minimizing the active power loss and the node voltage offset, an optimization model of the multienergy network is established, which uses the multiobjective nondominated sorting genetic algorithm to solve the installation location and the capacity of DGs in multienergy networks. In the second layer, according to the optimization results of the first layer and the characteristics of different DGs (wind power generator, photovoltaic panel, microturbine, and storage battery), the optimization model of the multienergy network is established to improve the economy, reliability, and environmental benefits of multienergy networks. It uses the multiobjective nondominated sorting genetic algorithm to solve the allocation capacity of different DGs so as to solve the optimal allocation problem of node capacity in multienergy networks. The double-layer optimization strategy proposed in this paper greatly promotes the development of multienergy networks and provides effective guidance for the optimal allocation and reliable operation of multienergy networks.


---

## Body

## 1. Introduction

Energy is the driving force and foundation for social and economic development. As traditional fossil energy sources are increasingly exhausted, improving energy efficiency, developing new energy sources, and strengthening multienergy complementarity have become inevitable choices to solve the contradiction between energy demand growth and energy shortage in the process of social and economic development. The multienergy network is a hybrid energy supply system that integrates multiple DGs, energy storage units, loads, and monitoring and protection devices. It can be flexibly connected to the grid or operated on isolated islands, which can effectively improve the flexibility, economy, and cleanliness of power system operation, and meet the requirements of users for power supply reliability, safety, and power quality. Since wind resources and solar resources are clean energy, the power generation of the micro gas turbine (MT) is stable and controllable, and the battery can suppress system power fluctuations. The node diagram of a typical multienergy network is shown in Figure1. Therefore, the multienergy network studied in this paper includes a wind turbine generator (WTG), a photovoltaic (PV) panel, an MT, and a storage battery (SB).Figure 1
Multienergy network node architecture.However, the installation location and capacity of DG and the allocation capacity of various DGs in multienergy networks will have a great impact on the operating technical indicators and economic indicators of the system. Therefore, it is of great importance and necessity to study the optimal allocation of DG in multienergy networks. At present, a lot of researches have been carried out on the optimal allocation of DGs for multienergy networks at home and abroad. For example, literature study [1] adopts a modified particle swarm optimization algorithm to optimize the cost, reliability, and pollutant emissions of hybrid power generation systems. Literature study [2] uses the harmony search algorithm to model and analyze the cost of hybrid systems which include a PV panel, a wind generator (WG), and diesel and then compares them with conventional PV panel/WG/diesel/battery systems. And literature study [3] applies a novel improved fruit fly optimization algorithm-based multiobjective optimization method to optimize the annual total cost and the pollutant emission of the stand-alone hybrid photovoltaic panel-wind generator-diesel-battery system. In addition, literature studies [4–7] have introduced multiobjective optimization algorithms for multiple distributed energy systems. With economic and environmental benefits as optimization objectives, different optimization algorithms are used to optimize the DG allocation capacity of multienergy networks. Besides, literature studies [8–10] introduce a specific allocation method of DG, which is aimed at low system loss or high voltage stability and solves the allocation capacity of DG. However, all of the above studies use the single-level multiobjective optimization function to solve the optimal allocation capacity of DG in multienergy networks when the location of DG is determined, while the installation location of DG needs to be determined by other algorithms.In this paper, firstly, the first-layer optimization objective function of the minimum active power loss and node voltage offset in the multienergy network system is constructed. The multiobjective optimization strategy based on the nondominated sorting genetic algorithm is used to solve the DG location and sizing. Secondly, in view of the power characteristics of various DGs and the optimization results of the first-layer algorithm, the second-layer algorithm takes into account the economic cost, power supply reliability, and environmental benefits and combines the multiobjective nondominated sorting genetic algorithm to optimize the allocation of various DG capacities. Finally, the proposed algorithm is applied to a multienergy network in a certain area. The results show that the operation performance, economy, reliability, and environmental protection of the multienergy network optimized by the double-layer nondominated sorting genetic algorithm are greatly improved.


## 2. Multiobjective Nondominated Sorting Genetic Algorithms

### 2.1. Mathematical Description of the Multiobjective Optimization Algorithm

Taking minimizing multiobjective function with constraints as an example, multiobjective optimization problems can be described as follows:(1)minfX=f1X,f2X,…,fnX,giX≤0,i=1,2,3,…,hjX=0,j=1,2,3,…,where fX is a vector space with n objective functions, giX is the ith inequality constraint function, hjX is the jth equation constraint, and X=x1,x2,…,xm is a vector space with m decision variables.


### 2.2. Decision-Making in Multiobjective Optimization

The multiobjective nondominated sorting genetic algorithm is used to obtain a series of discrete solutions [11–13]. How to choose an optimal solution among these discrete points is the key of the multiobjective optimization. In this paper, an optimal scheme selection method based on the geometric distance method is proposed. The central idea is that the geometric distance between the optimal solution and the positive ideal solution is the shortest, while the geometric distance between the optimal solution and the negative ideal solution is the longest. The optimal solution comes from the optimized Pareto frontier. The positive ideal solution refers to the solution that can satisfy the optimality of each objective function. And the negative ideal solution refers to the solution that does not satisfy any objective functions.The geometric distancedi+ of any solution i on the Pareto frontier to the positive ideal solution is calculated as follows:(2)di+=∑j=1mSij−Sjp−ideal2,where m is the number of objective functions, j is the dimension in which it is located, and Sij−Sjp−ideal is the linear distance of the jth dimension from the solution i to the positive ideal solution.The geometric distancedi− of any solution on the Pareto frontier to the negative ideal solution is calculated as follows:(3)di−=∑j=1mSij−Sjn−ideal2,where Sij−Sjn−ideal refers to the linear distance of the jth dimension between the solution i and the negative ideal solution.The coefficientωi is defined as follows:(4)ωi=di+di++di−.Theωi corresponding to any solutions in the Pareto frontier can be calculated by (1) to (3), and then the solution corresponding to the smallest ωi is selected, which is the final optimal solution.


## 2.1. Mathematical Description of the Multiobjective Optimization Algorithm

Taking minimizing multiobjective function with constraints as an example, multiobjective optimization problems can be described as follows:(1)minfX=f1X,f2X,…,fnX,giX≤0,i=1,2,3,…,hjX=0,j=1,2,3,…,where fX is a vector space with n objective functions, giX is the ith inequality constraint function, hjX is the jth equation constraint, and X=x1,x2,…,xm is a vector space with m decision variables.


## 2.2. Decision-Making in Multiobjective Optimization

The multiobjective nondominated sorting genetic algorithm is used to obtain a series of discrete solutions [11–13]. How to choose an optimal solution among these discrete points is the key of the multiobjective optimization. In this paper, an optimal scheme selection method based on the geometric distance method is proposed. The central idea is that the geometric distance between the optimal solution and the positive ideal solution is the shortest, while the geometric distance between the optimal solution and the negative ideal solution is the longest. The optimal solution comes from the optimized Pareto frontier. The positive ideal solution refers to the solution that can satisfy the optimality of each objective function. And the negative ideal solution refers to the solution that does not satisfy any objective functions.The geometric distancedi+ of any solution i on the Pareto frontier to the positive ideal solution is calculated as follows:(2)di+=∑j=1mSij−Sjp−ideal2,where m is the number of objective functions, j is the dimension in which it is located, and Sij−Sjp−ideal is the linear distance of the jth dimension from the solution i to the positive ideal solution.The geometric distancedi− of any solution on the Pareto frontier to the negative ideal solution is calculated as follows:(3)di−=∑j=1mSij−Sjn−ideal2,where Sij−Sjn−ideal refers to the linear distance of the jth dimension between the solution i and the negative ideal solution.The coefficientωi is defined as follows:(4)ωi=di+di++di−.Theωi corresponding to any solutions in the Pareto frontier can be calculated by (1) to (3), and then the solution corresponding to the smallest ωi is selected, which is the final optimal solution.


## 3. First layer: Location and Capacity Optimization of DG

### 3.1. Decision Variables of the First Layer

The rational planning of the DG in the multienergy network can effectively improve the operating performance of the system, so the decision variables of the first layer are the installation location and capacity of the DG.This paper uses a set of two-dimensional vectors to represent the decision variables of the first layer, namely, the optimal position and optimal capacity of DG, as shown in Table1.Table 1
Decision variables.
Unit 1
Unit 2
…
Unitn
Position
Capacity
Position
Capacity
…
Position
Capacity


### 3.2. Objective Function of the First-Layer Optimization Algorithm

Reasonable DG installation location and capacity configuration can improve power quality and optimize power flow while taking into account economic benefits. Moreover, some studies have testified that unreasonable locations or sizes of DG may result in greater system losses than the ones in the existing network [14]. The multiobjective optimization function should include economic indicators and technical indicators. The economic index of this layer mainly refers to the minimum active power loss of the multienergy network, and the technical index mainly refers to the minimum voltage offset of load nodes. The objective function is as follows (the function that has to be minimized consists of two objectives: one is economical and the other is technical):(5)minPloss=∑k=1N1Gki,jUi2+Uj2−2UiUjcosδij,minΔU=∑i=1NdUi−UispecΔUimax2,where Ploss is the active power loss of the multienergy network; ΔU is the offset of the load node voltage in the multienergy network; N1 represents the number of branches of the multienergy network; Nd represents the number of nodes; Gki,j represents the conductance of the branch k (i and j are the node numbers at both ends of the branch k); Ui and Uj represent the voltage amplitudes of the nodes i and j, respectively; δij refers to the phase angle difference of the nodes i and j; Uispec is the desired voltage value of the node i; and ΔUimax is the maximum allowable voltage deviation of the node i, i.e., ΔUimax=Uimax−Uimin.


### 3.3. Constraints of the First-Layer Optimization Algorithm

In multienergy networks, constraints are divided into equality constraints and inequality constraints.Equality constraints are mainly network power flow balancing. Equality constraints of the nodei are as follows:(6)PGi+PDGi−PLi−Ui∑j=1NdUjGijcosδij+Bijsinδij=0,QGi−QLi−Ui∑j=1NdUjGijsinδij+Bijcosδij=0,where PGi, PDGi, and PLi refer to the active power of the generator, DG, and load of the node i, respectively; Gij and Bij represent the conductance and admittance between nodes i and j, respectively; and QGi and QLi represent the reactive power of the generator and load of the node i, respectively.Inequality constraints mainly include the following aspects:(1)
Node voltage constraints:(7)Vimin≤Vi≤Vimax,i=1,2,…,Nd,whereViminandVimax are the upper and lower limits of the voltage of the node i, respectively.(2)
DG active power upper-limit constraints:(8)0≤PDGi≤PDGmax,i=1,2,…,NDG,wherePDGmax is the upper limit of the DG active power and NDG is the number of DGs.(3)
Branch transmission power constraints:(9)Pk≤Pkmax,k=1,2,…,N1,wherePkmax is the upper limit of the transmission power of the branch k.


## 3.1. Decision Variables of the First Layer

The rational planning of the DG in the multienergy network can effectively improve the operating performance of the system, so the decision variables of the first layer are the installation location and capacity of the DG.This paper uses a set of two-dimensional vectors to represent the decision variables of the first layer, namely, the optimal position and optimal capacity of DG, as shown in Table1.Table 1
Decision variables.
Unit 1
Unit 2
…
Unitn
Position
Capacity
Position
Capacity
…
Position
Capacity


## 3.2. Objective Function of the First-Layer Optimization Algorithm

Reasonable DG installation location and capacity configuration can improve power quality and optimize power flow while taking into account economic benefits. Moreover, some studies have testified that unreasonable locations or sizes of DG may result in greater system losses than the ones in the existing network [14]. The multiobjective optimization function should include economic indicators and technical indicators. The economic index of this layer mainly refers to the minimum active power loss of the multienergy network, and the technical index mainly refers to the minimum voltage offset of load nodes. The objective function is as follows (the function that has to be minimized consists of two objectives: one is economical and the other is technical):(5)minPloss=∑k=1N1Gki,jUi2+Uj2−2UiUjcosδij,minΔU=∑i=1NdUi−UispecΔUimax2,where Ploss is the active power loss of the multienergy network; ΔU is the offset of the load node voltage in the multienergy network; N1 represents the number of branches of the multienergy network; Nd represents the number of nodes; Gki,j represents the conductance of the branch k (i and j are the node numbers at both ends of the branch k); Ui and Uj represent the voltage amplitudes of the nodes i and j, respectively; δij refers to the phase angle difference of the nodes i and j; Uispec is the desired voltage value of the node i; and ΔUimax is the maximum allowable voltage deviation of the node i, i.e., ΔUimax=Uimax−Uimin.


## 3.3. Constraints of the First-Layer Optimization Algorithm

In multienergy networks, constraints are divided into equality constraints and inequality constraints.Equality constraints are mainly network power flow balancing. Equality constraints of the nodei are as follows:(6)PGi+PDGi−PLi−Ui∑j=1NdUjGijcosδij+Bijsinδij=0,QGi−QLi−Ui∑j=1NdUjGijsinδij+Bijcosδij=0,where PGi, PDGi, and PLi refer to the active power of the generator, DG, and load of the node i, respectively; Gij and Bij represent the conductance and admittance between nodes i and j, respectively; and QGi and QLi represent the reactive power of the generator and load of the node i, respectively.Inequality constraints mainly include the following aspects:(1)
Node voltage constraints:(7)Vimin≤Vi≤Vimax,i=1,2,…,Nd,whereViminandVimax are the upper and lower limits of the voltage of the node i, respectively.(2)
DG active power upper-limit constraints:(8)0≤PDGi≤PDGmax,i=1,2,…,NDG,wherePDGmax is the upper limit of the DG active power and NDG is the number of DGs.(3)
Branch transmission power constraints:(9)Pk≤Pkmax,k=1,2,…,N1,wherePkmax is the upper limit of the transmission power of the branch k.


## 4. Second Layer: Optimal Capacity Allocation for Various Types of DGs

### 4.1. DG Modeling

Multienergy networks generally include a WTG, a PV panel, an MT, and an SB [15]. This section establishes power models for various distributed power sources at first.


#### 4.1.1. WTGs

Because of the volatility of wind resources, the power output from the WTG is also unstable. According to the working principle of the WTG, when the wind speed is less than the cut-in wind speed, the WTG does not generate electricity; when the wind speed is greater than the cut-in wind speed and less than the rated wind speed, the WTG generates electricity, and the output electric power varies with the wind speed; when the wind speed is greater than the rated wind speed and less than the cut-out wind speed, the WTG outputs rated power; when the wind speed is greater than the cut-out wind speed, the WTG stops working and does not generate electricity [16].The outputPwindt of the WTG at time t is as follows [17]:(10)Pwindt=0,vt≤vinorvt≥vout,vt3−vin3vr3−vt3∗Pr,vin≤vt≤vr,Pr,vr≤vt≤vout,where vt represents the wind speed at time t; vin, vr, and vout represent the cut-in wind speed, the rated wind speed, and the cut-out wind speed, respectively; and Pr refers to the rated power of the WTG.


#### 4.1.2. PV Panels

The output voltage and current of the PV panel vary with the change of the illumination intensity and the junction temperature of the battery. They have strong nonlinear characteristics, and there is one maximum power output point under certain working conditions. In order to make full use of solar energy for PV cells, it is necessary to make the working point of PV cells fall at the maximum power point when the illumination intensity and temperature change. Under the maximum power mode, the expressions of powerPPVt, voltage VPVt, and current IPVt, that is, the optimum working point t of the PV cell, are as follows [18]:(11)PPVt=VPVt∗IPVt,VPVt=Vmp∗1+0.0539∗logHTt1000+β∗ΔTt,IPVt=ISC∗1−a∗expVPVt−ΔVtb∗VOC−1+ΔIt,a=1−ImpISC∗exp−Vmpb∗VOC,b=Vmp/VOC−1ln1−Imp/ISC,ΔTt=TAt+0.02∗HTt−25,ΔVt=VPVt−Vmp,ΔIt=α∗HTt1000∗ΔTt+HTt1000−1∗ISC,where VmpandImp represent the maximum power point voltage and current of the PV cell, respectively; VOC represents the open circuit voltage; ISC represents the short circuit current; HTt represents the illumination radiance at time t; α represents the current temperature coefficient of the PV cell module; β represents the voltage temperature coefficient of the PV cell module; and TAt represents the ambient temperature at time t.


#### 4.1.3. MTs

MTs are widely used in distributed systems because of their high reliability, long service life, low environmental pollution, and flexible control. Their working principle is that natural gas and high-pressure gas are mixed and burned in the combustion chamber to produce high-quality gas that drives compressors and generators to generate heat and electricity. Suppose that a certain type of MT is operating under full operating conditions and that the electric load is supplied while satisfying the user’s cooling and heating load demands. The output power of the MT at timet is as follows [19]:(12)PMTt=ρ∗Vint∗LH+CP∗ρ∗Vint∗T2t−T1t,where PMTt represents the power output of the MT at time t, Vint represents the fuel intake amount at time t, ρ is the gas density, LH is the low heating value, CP represents the specific heat capacity, and T1tandT2t indicate the intake air temperature and exhaust gas temperature of the fuel at time t, respectively.


#### 4.1.4. SBs

The SB can suppress the fluctuation of distributed power supply, has the performance of peak shaving and valley filling, and maintains power stability. The output power of the SB is closely related to the operation status of the system [20].Assuming that lead-acid batteries are used, the state of timet is related to the state of time t − 1 and the supply and demand of electricity from time t − 1 to time t. When the total output power of the DG is greater than the load power consumption, the battery pack is in a charging state; otherwise, the battery pack is in a discharging state. The charge of the battery pack at time t can be expressed as follows [19]:(13)PSBt=PSBt−1+Ptotalt−Ploadtηinv∗ηcharge;charge,PSBt−1−Ploadt/ηinv−Ptotaltηdischarge;discharge,where Ptotalt=Pwindt+PPVt+PMTt, PSBt represents the amount of electricity stored in the SB at time t, Ptotalt is the amount of electricity supplied by the DG in the multienergy network at time t, Ploadt is the load of the system at time t, ηinv represents the conversion efficiency of the inverter, ηcharge is the charging efficiency of the SB, and ηdischarge is the discharging efficiency of the SB.


### 4.2. Decision Variables of the Second Layer

The multiobjective optimization algorithm of the first-layer multienergy network has obtained the distributed power supply installation location and total capacity. But the construction cost, operation cost, output stability, and environmental impact of different types of distributed power supplies are different, so it is necessary to carry out the second-layer multiobjective optimization to obtain the optimal capacity allocation of various distributed power sources.The main decision variables at the second layer of the network are as follows:(14)X=Nwind,NPV,NMT,NSB,where Nwind,NPV,NMT,andNSB represent the number of WTGs, PV cells, MTs, and SBs, respectively.


### 4.3. Objective Function of the Second-Layer Optimization Algorithm

The objectives of the second-layer optimization algorithm are mainly to consider economy, power supply reliability, and environmental protection. The objective functions are as follows.


#### 4.3.1. Economy

The economy of the multienergy network is mainly reflected by the construction cost and operation cost of the DG system [21]. The objective functions of the economy are as follows:(15)minCtotal=minCCpt+CMtn,CCpt=CCptwind+CCptPV+CCptMT+CCptSB,CCptwind=Nwind∗ewind∗Pwind∗r∗1+rmwind1+rmwind−1,CCptPV=NPV∗ePV∗PPV∗r∗1+rmPV1+rmPV−1,CCptMT=NMT∗eMT∗PMT∗r∗1+rmMT1+rmMT−1,CCptSB=∗eSB∗PSB∗r∗1+rmSB1+rmSB−1,CMtn=Nwind∗μPwind+NPV∗μPPV+NMT∗μPMT+NSB∗μPSB,where Ctotal is the total cost of power generation; CCpt and CMtn represent the cost of construction and operation and maintenance of DG in the multienergy network, respectively; Pwind,PPV,PMT,andPSB represent the output power of each WTG, PV cell, MT, and SB, respectively; ewind,ePV,eMT,andeSB represent the unit cost of the WTG, PV cell, MT, and SB, respectively; μPwind,μPPV,μPMT,andμPSB represent the operation and maintenance costs of each WTG, PV cell, MT, and SB, respectively; mwind,mPV,mMT,andmSB represent the depreciation life of the WTG, PV cell, MT, and SB, respectively; and r is the discount rate of the equipment, generally 8%.


#### 4.3.2. Power Supply Reliability

Multienergy networks contain fluctuating DGs, and the load of multienergy networks also changes with time [22]. Therefore, it is necessary to evaluate the power supply reliability of the system; that is, the rate of load shortage is required. The smaller the rate of power shortage, the higher the reliability of power supply, and vice versa [23]. This paper takes one day (T = 24 h) as the period for evaluating the power supply reliability of the system and divides this period into 24 sections. Assuming that the wind speed, light, temperature, and load are constant in each section (1 h), the load shortage rate in one evaluation period of the system is as follows:(16)minLSRP=∑t=1TPloadt−PDGt∗ηinvPloadt∗Ut,PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt−PSBmin,Ut=1,PDGt<Ploadt,0,PDGt≥Ploadt,where LSRP represents the one-day load shortage rate of the system; Ploadt represents the load at time t; Pwindt,PPVt,PMTt,andPSBt refer to the output power of the WTG, PV cell, MT, and SB at time t, respectively; PSBmin represents the minimum residual power allowed by the SB; and ηinv represents the efficiency of the inverter.


#### 4.3.3. Environmental Protection

Environmental protection is mainly reflected by the emission of pollutants (NOx, CO, CO2, and SO2) from DGs. WTG power generation and PV power generation belong to clean energy generation, and they do not emit pollutants. Although the output power of the MT can be flexibly controlled, it will emit pollutants. The higher the emission of pollutants, the worse the environmental protection [24]. The environmental cost of power generation EMplt in the evaluation cycle is as follows:(17)minEMplt=∑t=1T∑i=1nNMT∗PMTt∗αi∗Qi,where αi refers to the ith emission of pollutants of unit power of the MT, Qi is the treatment fee of ith pollutants of unit power of emission, PMTt denotes the output power of the MT at time t, and NMT is the input quantity of the MT in the system.


### 4.4. Constraints of the Second-Layer Optimization Algorithm

#### 4.4.1. Power Balance Constraint

In this paper, the DG capacity of the multienergy network solved by the first-layer optimization algorithm is the sum of all kinds of DGs in the second layer so as to ensure the minimum active power loss and voltage offset of the multienergy network. The power constraints are as follows:(18)PDGt−PDG∗tPDG∗t≤δt,where PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt, PDG∗t denotes the installed capacity of DGs obtained from the first-layer optimization algorithm at time t, and δt denotes the maximum unbalanced rate of power that the system can withstand at time t.


#### 4.4.2. Bounds of Design Variables

The bounds of design variables are as follows:(19)0≤Nwind≤Nwindmax,0≤NPV≤NPVmax,NMTmin≤NMT≤NMTmax,0≤NSB≤NSBmax,where Nwindmax,NPVmax,NMTmax,andNSBmax, respectively, represent the maximum number of WTGs, PV cells, MTs, and SBs that meets the system load demand and NMTmin is the minimum number of MTs that meets the condition of “fixing electricity by heat.”


#### 4.4.3. Battery Constrains

The battery constraints are as follows:(20)PSBmin≤PSBt≤PSBmax,where PSBmax represents the maximum storage capacity allowed by the SB.


## 4.1. DG Modeling

Multienergy networks generally include a WTG, a PV panel, an MT, and an SB [15]. This section establishes power models for various distributed power sources at first.


### 4.1.1. WTGs

Because of the volatility of wind resources, the power output from the WTG is also unstable. According to the working principle of the WTG, when the wind speed is less than the cut-in wind speed, the WTG does not generate electricity; when the wind speed is greater than the cut-in wind speed and less than the rated wind speed, the WTG generates electricity, and the output electric power varies with the wind speed; when the wind speed is greater than the rated wind speed and less than the cut-out wind speed, the WTG outputs rated power; when the wind speed is greater than the cut-out wind speed, the WTG stops working and does not generate electricity [16].The outputPwindt of the WTG at time t is as follows [17]:(10)Pwindt=0,vt≤vinorvt≥vout,vt3−vin3vr3−vt3∗Pr,vin≤vt≤vr,Pr,vr≤vt≤vout,where vt represents the wind speed at time t; vin, vr, and vout represent the cut-in wind speed, the rated wind speed, and the cut-out wind speed, respectively; and Pr refers to the rated power of the WTG.


### 4.1.2. PV Panels

The output voltage and current of the PV panel vary with the change of the illumination intensity and the junction temperature of the battery. They have strong nonlinear characteristics, and there is one maximum power output point under certain working conditions. In order to make full use of solar energy for PV cells, it is necessary to make the working point of PV cells fall at the maximum power point when the illumination intensity and temperature change. Under the maximum power mode, the expressions of powerPPVt, voltage VPVt, and current IPVt, that is, the optimum working point t of the PV cell, are as follows [18]:(11)PPVt=VPVt∗IPVt,VPVt=Vmp∗1+0.0539∗logHTt1000+β∗ΔTt,IPVt=ISC∗1−a∗expVPVt−ΔVtb∗VOC−1+ΔIt,a=1−ImpISC∗exp−Vmpb∗VOC,b=Vmp/VOC−1ln1−Imp/ISC,ΔTt=TAt+0.02∗HTt−25,ΔVt=VPVt−Vmp,ΔIt=α∗HTt1000∗ΔTt+HTt1000−1∗ISC,where VmpandImp represent the maximum power point voltage and current of the PV cell, respectively; VOC represents the open circuit voltage; ISC represents the short circuit current; HTt represents the illumination radiance at time t; α represents the current temperature coefficient of the PV cell module; β represents the voltage temperature coefficient of the PV cell module; and TAt represents the ambient temperature at time t.


### 4.1.3. MTs

MTs are widely used in distributed systems because of their high reliability, long service life, low environmental pollution, and flexible control. Their working principle is that natural gas and high-pressure gas are mixed and burned in the combustion chamber to produce high-quality gas that drives compressors and generators to generate heat and electricity. Suppose that a certain type of MT is operating under full operating conditions and that the electric load is supplied while satisfying the user’s cooling and heating load demands. The output power of the MT at timet is as follows [19]:(12)PMTt=ρ∗Vint∗LH+CP∗ρ∗Vint∗T2t−T1t,where PMTt represents the power output of the MT at time t, Vint represents the fuel intake amount at time t, ρ is the gas density, LH is the low heating value, CP represents the specific heat capacity, and T1tandT2t indicate the intake air temperature and exhaust gas temperature of the fuel at time t, respectively.


### 4.1.4. SBs

The SB can suppress the fluctuation of distributed power supply, has the performance of peak shaving and valley filling, and maintains power stability. The output power of the SB is closely related to the operation status of the system [20].Assuming that lead-acid batteries are used, the state of timet is related to the state of time t − 1 and the supply and demand of electricity from time t − 1 to time t. When the total output power of the DG is greater than the load power consumption, the battery pack is in a charging state; otherwise, the battery pack is in a discharging state. The charge of the battery pack at time t can be expressed as follows [19]:(13)PSBt=PSBt−1+Ptotalt−Ploadtηinv∗ηcharge;charge,PSBt−1−Ploadt/ηinv−Ptotaltηdischarge;discharge,where Ptotalt=Pwindt+PPVt+PMTt, PSBt represents the amount of electricity stored in the SB at time t, Ptotalt is the amount of electricity supplied by the DG in the multienergy network at time t, Ploadt is the load of the system at time t, ηinv represents the conversion efficiency of the inverter, ηcharge is the charging efficiency of the SB, and ηdischarge is the discharging efficiency of the SB.


## 4.1.1. WTGs

Because of the volatility of wind resources, the power output from the WTG is also unstable. According to the working principle of the WTG, when the wind speed is less than the cut-in wind speed, the WTG does not generate electricity; when the wind speed is greater than the cut-in wind speed and less than the rated wind speed, the WTG generates electricity, and the output electric power varies with the wind speed; when the wind speed is greater than the rated wind speed and less than the cut-out wind speed, the WTG outputs rated power; when the wind speed is greater than the cut-out wind speed, the WTG stops working and does not generate electricity [16].The outputPwindt of the WTG at time t is as follows [17]:(10)Pwindt=0,vt≤vinorvt≥vout,vt3−vin3vr3−vt3∗Pr,vin≤vt≤vr,Pr,vr≤vt≤vout,where vt represents the wind speed at time t; vin, vr, and vout represent the cut-in wind speed, the rated wind speed, and the cut-out wind speed, respectively; and Pr refers to the rated power of the WTG.


## 4.1.2. PV Panels

The output voltage and current of the PV panel vary with the change of the illumination intensity and the junction temperature of the battery. They have strong nonlinear characteristics, and there is one maximum power output point under certain working conditions. In order to make full use of solar energy for PV cells, it is necessary to make the working point of PV cells fall at the maximum power point when the illumination intensity and temperature change. Under the maximum power mode, the expressions of powerPPVt, voltage VPVt, and current IPVt, that is, the optimum working point t of the PV cell, are as follows [18]:(11)PPVt=VPVt∗IPVt,VPVt=Vmp∗1+0.0539∗logHTt1000+β∗ΔTt,IPVt=ISC∗1−a∗expVPVt−ΔVtb∗VOC−1+ΔIt,a=1−ImpISC∗exp−Vmpb∗VOC,b=Vmp/VOC−1ln1−Imp/ISC,ΔTt=TAt+0.02∗HTt−25,ΔVt=VPVt−Vmp,ΔIt=α∗HTt1000∗ΔTt+HTt1000−1∗ISC,where VmpandImp represent the maximum power point voltage and current of the PV cell, respectively; VOC represents the open circuit voltage; ISC represents the short circuit current; HTt represents the illumination radiance at time t; α represents the current temperature coefficient of the PV cell module; β represents the voltage temperature coefficient of the PV cell module; and TAt represents the ambient temperature at time t.


## 4.1.3. MTs

MTs are widely used in distributed systems because of their high reliability, long service life, low environmental pollution, and flexible control. Their working principle is that natural gas and high-pressure gas are mixed and burned in the combustion chamber to produce high-quality gas that drives compressors and generators to generate heat and electricity. Suppose that a certain type of MT is operating under full operating conditions and that the electric load is supplied while satisfying the user’s cooling and heating load demands. The output power of the MT at timet is as follows [19]:(12)PMTt=ρ∗Vint∗LH+CP∗ρ∗Vint∗T2t−T1t,where PMTt represents the power output of the MT at time t, Vint represents the fuel intake amount at time t, ρ is the gas density, LH is the low heating value, CP represents the specific heat capacity, and T1tandT2t indicate the intake air temperature and exhaust gas temperature of the fuel at time t, respectively.


## 4.1.4. SBs

The SB can suppress the fluctuation of distributed power supply, has the performance of peak shaving and valley filling, and maintains power stability. The output power of the SB is closely related to the operation status of the system [20].Assuming that lead-acid batteries are used, the state of timet is related to the state of time t − 1 and the supply and demand of electricity from time t − 1 to time t. When the total output power of the DG is greater than the load power consumption, the battery pack is in a charging state; otherwise, the battery pack is in a discharging state. The charge of the battery pack at time t can be expressed as follows [19]:(13)PSBt=PSBt−1+Ptotalt−Ploadtηinv∗ηcharge;charge,PSBt−1−Ploadt/ηinv−Ptotaltηdischarge;discharge,where Ptotalt=Pwindt+PPVt+PMTt, PSBt represents the amount of electricity stored in the SB at time t, Ptotalt is the amount of electricity supplied by the DG in the multienergy network at time t, Ploadt is the load of the system at time t, ηinv represents the conversion efficiency of the inverter, ηcharge is the charging efficiency of the SB, and ηdischarge is the discharging efficiency of the SB.


## 4.2. Decision Variables of the Second Layer

The multiobjective optimization algorithm of the first-layer multienergy network has obtained the distributed power supply installation location and total capacity. But the construction cost, operation cost, output stability, and environmental impact of different types of distributed power supplies are different, so it is necessary to carry out the second-layer multiobjective optimization to obtain the optimal capacity allocation of various distributed power sources.The main decision variables at the second layer of the network are as follows:(14)X=Nwind,NPV,NMT,NSB,where Nwind,NPV,NMT,andNSB represent the number of WTGs, PV cells, MTs, and SBs, respectively.


## 4.3. Objective Function of the Second-Layer Optimization Algorithm

The objectives of the second-layer optimization algorithm are mainly to consider economy, power supply reliability, and environmental protection. The objective functions are as follows.


### 4.3.1. Economy

The economy of the multienergy network is mainly reflected by the construction cost and operation cost of the DG system [21]. The objective functions of the economy are as follows:(15)minCtotal=minCCpt+CMtn,CCpt=CCptwind+CCptPV+CCptMT+CCptSB,CCptwind=Nwind∗ewind∗Pwind∗r∗1+rmwind1+rmwind−1,CCptPV=NPV∗ePV∗PPV∗r∗1+rmPV1+rmPV−1,CCptMT=NMT∗eMT∗PMT∗r∗1+rmMT1+rmMT−1,CCptSB=∗eSB∗PSB∗r∗1+rmSB1+rmSB−1,CMtn=Nwind∗μPwind+NPV∗μPPV+NMT∗μPMT+NSB∗μPSB,where Ctotal is the total cost of power generation; CCpt and CMtn represent the cost of construction and operation and maintenance of DG in the multienergy network, respectively; Pwind,PPV,PMT,andPSB represent the output power of each WTG, PV cell, MT, and SB, respectively; ewind,ePV,eMT,andeSB represent the unit cost of the WTG, PV cell, MT, and SB, respectively; μPwind,μPPV,μPMT,andμPSB represent the operation and maintenance costs of each WTG, PV cell, MT, and SB, respectively; mwind,mPV,mMT,andmSB represent the depreciation life of the WTG, PV cell, MT, and SB, respectively; and r is the discount rate of the equipment, generally 8%.


### 4.3.2. Power Supply Reliability

Multienergy networks contain fluctuating DGs, and the load of multienergy networks also changes with time [22]. Therefore, it is necessary to evaluate the power supply reliability of the system; that is, the rate of load shortage is required. The smaller the rate of power shortage, the higher the reliability of power supply, and vice versa [23]. This paper takes one day (T = 24 h) as the period for evaluating the power supply reliability of the system and divides this period into 24 sections. Assuming that the wind speed, light, temperature, and load are constant in each section (1 h), the load shortage rate in one evaluation period of the system is as follows:(16)minLSRP=∑t=1TPloadt−PDGt∗ηinvPloadt∗Ut,PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt−PSBmin,Ut=1,PDGt<Ploadt,0,PDGt≥Ploadt,where LSRP represents the one-day load shortage rate of the system; Ploadt represents the load at time t; Pwindt,PPVt,PMTt,andPSBt refer to the output power of the WTG, PV cell, MT, and SB at time t, respectively; PSBmin represents the minimum residual power allowed by the SB; and ηinv represents the efficiency of the inverter.


### 4.3.3. Environmental Protection

Environmental protection is mainly reflected by the emission of pollutants (NOx, CO, CO2, and SO2) from DGs. WTG power generation and PV power generation belong to clean energy generation, and they do not emit pollutants. Although the output power of the MT can be flexibly controlled, it will emit pollutants. The higher the emission of pollutants, the worse the environmental protection [24]. The environmental cost of power generation EMplt in the evaluation cycle is as follows:(17)minEMplt=∑t=1T∑i=1nNMT∗PMTt∗αi∗Qi,where αi refers to the ith emission of pollutants of unit power of the MT, Qi is the treatment fee of ith pollutants of unit power of emission, PMTt denotes the output power of the MT at time t, and NMT is the input quantity of the MT in the system.


## 4.3.1. Economy

The economy of the multienergy network is mainly reflected by the construction cost and operation cost of the DG system [21]. The objective functions of the economy are as follows:(15)minCtotal=minCCpt+CMtn,CCpt=CCptwind+CCptPV+CCptMT+CCptSB,CCptwind=Nwind∗ewind∗Pwind∗r∗1+rmwind1+rmwind−1,CCptPV=NPV∗ePV∗PPV∗r∗1+rmPV1+rmPV−1,CCptMT=NMT∗eMT∗PMT∗r∗1+rmMT1+rmMT−1,CCptSB=∗eSB∗PSB∗r∗1+rmSB1+rmSB−1,CMtn=Nwind∗μPwind+NPV∗μPPV+NMT∗μPMT+NSB∗μPSB,where Ctotal is the total cost of power generation; CCpt and CMtn represent the cost of construction and operation and maintenance of DG in the multienergy network, respectively; Pwind,PPV,PMT,andPSB represent the output power of each WTG, PV cell, MT, and SB, respectively; ewind,ePV,eMT,andeSB represent the unit cost of the WTG, PV cell, MT, and SB, respectively; μPwind,μPPV,μPMT,andμPSB represent the operation and maintenance costs of each WTG, PV cell, MT, and SB, respectively; mwind,mPV,mMT,andmSB represent the depreciation life of the WTG, PV cell, MT, and SB, respectively; and r is the discount rate of the equipment, generally 8%.


## 4.3.2. Power Supply Reliability

Multienergy networks contain fluctuating DGs, and the load of multienergy networks also changes with time [22]. Therefore, it is necessary to evaluate the power supply reliability of the system; that is, the rate of load shortage is required. The smaller the rate of power shortage, the higher the reliability of power supply, and vice versa [23]. This paper takes one day (T = 24 h) as the period for evaluating the power supply reliability of the system and divides this period into 24 sections. Assuming that the wind speed, light, temperature, and load are constant in each section (1 h), the load shortage rate in one evaluation period of the system is as follows:(16)minLSRP=∑t=1TPloadt−PDGt∗ηinvPloadt∗Ut,PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt−PSBmin,Ut=1,PDGt<Ploadt,0,PDGt≥Ploadt,where LSRP represents the one-day load shortage rate of the system; Ploadt represents the load at time t; Pwindt,PPVt,PMTt,andPSBt refer to the output power of the WTG, PV cell, MT, and SB at time t, respectively; PSBmin represents the minimum residual power allowed by the SB; and ηinv represents the efficiency of the inverter.


## 4.3.3. Environmental Protection

Environmental protection is mainly reflected by the emission of pollutants (NOx, CO, CO2, and SO2) from DGs. WTG power generation and PV power generation belong to clean energy generation, and they do not emit pollutants. Although the output power of the MT can be flexibly controlled, it will emit pollutants. The higher the emission of pollutants, the worse the environmental protection [24]. The environmental cost of power generation EMplt in the evaluation cycle is as follows:(17)minEMplt=∑t=1T∑i=1nNMT∗PMTt∗αi∗Qi,where αi refers to the ith emission of pollutants of unit power of the MT, Qi is the treatment fee of ith pollutants of unit power of emission, PMTt denotes the output power of the MT at time t, and NMT is the input quantity of the MT in the system.


## 4.4. Constraints of the Second-Layer Optimization Algorithm

### 4.4.1. Power Balance Constraint

In this paper, the DG capacity of the multienergy network solved by the first-layer optimization algorithm is the sum of all kinds of DGs in the second layer so as to ensure the minimum active power loss and voltage offset of the multienergy network. The power constraints are as follows:(18)PDGt−PDG∗tPDG∗t≤δt,where PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt, PDG∗t denotes the installed capacity of DGs obtained from the first-layer optimization algorithm at time t, and δt denotes the maximum unbalanced rate of power that the system can withstand at time t.


### 4.4.2. Bounds of Design Variables

The bounds of design variables are as follows:(19)0≤Nwind≤Nwindmax,0≤NPV≤NPVmax,NMTmin≤NMT≤NMTmax,0≤NSB≤NSBmax,where Nwindmax,NPVmax,NMTmax,andNSBmax, respectively, represent the maximum number of WTGs, PV cells, MTs, and SBs that meets the system load demand and NMTmin is the minimum number of MTs that meets the condition of “fixing electricity by heat.”


### 4.4.3. Battery Constrains

The battery constraints are as follows:(20)PSBmin≤PSBt≤PSBmax,where PSBmax represents the maximum storage capacity allowed by the SB.


## 4.4.1. Power Balance Constraint

In this paper, the DG capacity of the multienergy network solved by the first-layer optimization algorithm is the sum of all kinds of DGs in the second layer so as to ensure the minimum active power loss and voltage offset of the multienergy network. The power constraints are as follows:(18)PDGt−PDG∗tPDG∗t≤δt,where PDGt=Nwind∗Pwindt+NPV∗PPVt+NMT∗PMTt+NSB∗PSBt, PDG∗t denotes the installed capacity of DGs obtained from the first-layer optimization algorithm at time t, and δt denotes the maximum unbalanced rate of power that the system can withstand at time t.


## 4.4.2. Bounds of Design Variables

The bounds of design variables are as follows:(19)0≤Nwind≤Nwindmax,0≤NPV≤NPVmax,NMTmin≤NMT≤NMTmax,0≤NSB≤NSBmax,where Nwindmax,NPVmax,NMTmax,andNSBmax, respectively, represent the maximum number of WTGs, PV cells, MTs, and SBs that meets the system load demand and NMTmin is the minimum number of MTs that meets the condition of “fixing electricity by heat.”


## 4.4.3. Battery Constrains

The battery constraints are as follows:(20)PSBmin≤PSBt≤PSBmax,where PSBmax represents the maximum storage capacity allowed by the SB.


## 5. Algorithmic Flow

The flow chart of the optimal allocation strategy for the multienergy network based on the double-layer nondominated sorting genetic algorithm proposed in this paper is shown in Figure2.Figure 2
Flow chart of the algorithm.


## 6. Case Study

In this paper, a 9-node multienergy network system is taken as an example to verify the optimization effect of the double-layer nondominated sorting genetic algorithm. The topological structure of this multienergy network system is shown in Figure3. Among them, node 1 is a balanced node, node 3 is a PV node, and the rest nodes are all PQ nodes. The grid parameters (unit values) are shown in Table 2. The evaluation period T of this example is chosen as one day (24 h), and the load data (per unit) and meteorological data in the evaluation period are shown in Figures 4–7 separately. The parameters of lines, WTG, PV panel, MT, and SB and pollutant treatment fee are shown in Tables 2–8, respectively [25].Figure 3
Multienergy network architecture.Table 2
Line parameters.
Initial points of the line
End of the line
Resistance
Inductance
Capacitance
1
4
0.000
0.025
0.000
4
5
0.007
0.040
0.017
5
6
0.016
0.073
0.039
3
6
0.000
0.025
0.000
6
7
0.005
0.043
0.023
7
8
0.004
0.031
0.016
8
2
0.000
0.027
0.000
8
9
0.014
0.070
0.033
9
4
0.004
0.037
0.019Figure 4
Active power data of the load at each node.Figure 5
Reactive power data of the load at each node.Figure 6
Wind resource data.Figure 7
Solar resource data. (a) Radiation. (b) Temperature.
(a)
(b)Table 3
Wind turbine generator parameters.
P
r (kW)
V
in (m/s)
V
r (m/s)
V
out (m/s)
N
wind
max
500
3
12
25
8Table 4
PV panel parameters.
P
r (kW)
V
OC (V)
I
sc (A)
V
mp (V)
I
mp (A)
Series number
Parallel number
N
PV
max
100
64.6
6.14
54.7
5.76
5
64
30Table 5
MT parameters.
P
r (kW)
C
P (kJ/kg·K)
LH (MJ/m³)
ρ (kg/m³)
T
1 (°C)
T
2 (°C)
N
MT
max
N
MT
min
300
2.16
32.6
0.75
25
280
20
1Table 6
Battery parameters.
P
r (kW)
P
SB ini (kW)
P
SB
min (kW)
P
SB
max (kW)
η
charge
η
discharge
η
inv
N
SB
max
24
9.6
4.8
24
0.85
0.98
0.87
40Table 7
Cost values.
Type
Life span (years)
Unit cost ($/kW)
Maintenance cost ($/kW)
WTG
20
523
38
PV panel
20
508
29
MT
20
1160
92
SB
10
156
22Table 8
Pollutant emission of the MT.
Type
NOx
CO2
CO
SO2
Value (g/kWh)
0.6188
184.0829
0.1702
0.000928
Maintenance ($/kg)
0.250
0.00125
0.020
0.125The optimal allocation strategy proposed in this paper is used to optimize the 9-node multienergy network. The population size of the double-layer optimization algorithm is 200, the number of iterations is 100, and there are two locations in the network where DG can be installed. The first-layer multiobjective optimization algorithm optimizes the location capacity of 24-hour DG in real time based on the 24-hour real-time load data of the network. The second-layer multiobjective optimization algorithm solves the capacity optimization of various DGs with the constraints of the first-layer optimization results, equipment characteristics, and design.The results of the first-layer optimization algorithm are shown in Tables9 and 10, and the Pareto solution of 100 generations of evolution at a certain time (such as at 13:00) is shown in Figure 8.Table 9
The first-layer decision variables.
Time
First location
First sizing
Second location
Second sizing
1:00
8
0.7587
9
0.2803
2:00
8
0.2184
9
0.8181
3:00
8
0.4325
9
0.5976
4:00
8
0.7628
9
0.2839
5:00
8
0.7074
9
0.3300
6:00
8
0.7541
9
0.2678
7:00
8
0.2695
9
0.7638
8:00
8
0.7190
9
0.7229
9:00
8
0.8485
9
0.7799
10:00
8
0.8922
9
0.7764
11:00
8
0.8608
9
0.8028
12:00
8
0.8887
9
0.8094
13:00
8
0.3797
9
0.4001
14:00
8
0.7833
9
0.9528
15:00
8
0.8063
9
0.9492
16:00
8
0.9274
9
0.7632
17:00
8
0.8883
9
0.7713
18:00
8
0.7370
9
0.8519
19:00
8
0.5159
9
0.7931
20:00
8
0.4600
9
0.7055
21:00
8
0.8561
9
0.3772
22:00
8
0.7965
9
0.3319
23:00
8
0.6882
9
0.3499
24:00
8
0.5463
9
0.4537Table 10
Comparison of the optimization index.
System loss (×10−3)
Voltage deviation
Before optimization
8.435
0.416
After optimization
0.413∼5.603
0.015∼0.351
Final optimal solution
3.815
0.215Figure 8
The Pareto front result at 13:00.The above operation results show that, after adopting the optimal allocation strategy of the multienergy network based on the double-layer nondominated sorting genetic algorithm, the active power loss of the multienergy network is reduced by 54.77% and the system voltage deviation is reduced by 48.32%. This proves that this strategy can effectively improve the operating conditions of the system, guarantee the stability and economic benefits, and improve the power quality.The results of the second-layer optimization algorithm are shown in Table11, and the Pareto solution of the 100-generation evolution is shown in Figure 9. Comparing the optimization solution of this paper with the nonoptimization solution, the results are shown in Table 12. Nonoptimization solution only consists of a single type of DG such as WTG, PV, MT, and SB, and its capacity is capped.Table 11
The second-layer decision variables.
DG location
N
wind
N
PV
N
MT
N
SB
First
2
7
17
7
Second
1
8
13
23Figure 9
The second Pareto front result.Table 12
Comparison of the optimization index.
N
wind
N
PV
N
MT
N
SB
C
total ($)
LSRP
EM
plt ($)
Without optimization
16
0
0
0
949.63
6.64
0.0
0
60
0
0
561.77
10.24
0.0
0
0
40
0
2656.32
3.49
111.84
0
0
0
80
127.96
19.57
0.0
16
60
40
80
4295.67
0.0
111.84
With optimization
3
15
30
30
3448.83
0.26
81.08When there are no optimization and only a single type of DG, although the costs of construction, operation, and maintenance are relatively low, the power shortage rate and the cost of pollutant treatment are relatively high. When there are no optimization and all kinds of DGs, although the load power shortage rate is 0, the cost of construction, operation, maintenance, and pollutant treatment is relatively high. Considering the economy, reliability, and environmental protection of the system at the same time, it is necessary to optimize DG allocation of the multienergy network by using the double-layer nondominated sorting genetic algorithm. After optimization, the relative environmental cost and load power shortage rate of the multienergy network are maintained at a relatively low level, which improves the economic effect, environmental protection benefits, and power supply reliability.


## 7. Conclusion

Aiming at the complex and changeable multienergy network, this paper proposes a multienergy network optimal allocation strategy based on the double-layer nondominated sorting genetic algorithm. Firstly, according to the network load and grid structure, the first-layer optimization algorithm calculates the optimal location and capacity of DG in the multienergy network with the purpose of reducing the network loss and the system voltage offset. Subsequently, the second-layer optimization algorithm takes the first-layer calculation results as constraints and calculates the optimal allocation of different types of DGs, whose purpose is to improve economy, power supply reliability, and environmental protection. The double-layer optimization algorithm uses the optimal scheme selection method based on the geometric distance method to solve the decision variables. Finally, the case study in this paper shows that the double-layer optimal allocation strategy can effectively reduce network loss, improve system voltage, reduce system cost, improve system power supply reliability, increase system environmental benefits, and provide guidance for the design and optimal operation of multienergy networks.



---

*Source: 5367403-2019-10-29.xml*