A Multi-Objective Evaluation Method for Distributed ...

0 downloads 0 Views 388KB Size Report
Initial investment/Operation cost. inj. M. Cost of investment for equipment j in the system. R. Capital recovery factor. unitj. M. Equipment's initial investment for unit ...
$0XOWL2EMHFWLYH(YDOXDWLRQ0HWKRGIRU'LVWULEXWHG ,QWHJUDWHG(QHUJ\6\VWHP Cong Chen

Yixun Xue

Xinwei Shen

Qinglai Guo

Tsinghua-Berkeley Shenzhen Institute

Hongbin Sun

Shenzhen, China

Dept. of Electrical Engineering, Tsinghua University

[email protected]

Beijing, China

Abstract—In this paper, a multi-objective mathematical

a / b / c1 / c2 Coefficient of the linear model in the CCHP.

model concerning the optimization of ecology, economy and

Q11

environment benefits in a distributed integrated energy system

Q6 / Q7 / Q8 Cooling/heating/power load of the system.

(IES) is proposed. Compared with the traditional evaluation methods, e.g. Analytic Hierarchy Process (AHP), this method is based on the Pareto Front (PF) of the multi-objective optimization model, in which the multi-objective potential variant is defined to evaluate the potential of the system making the evaluation more convincing and indisputable. Moreover, to avoid the phenomenon

quantity and ignoring the quality, we utilize the Exergy Efficiency to

amend

the

evaluation

based

on

theories

D 6 /D 7 'T Ti M in / M run

Conversion factor.

M inj

Cost of investment for equipment j in the system.

R M unitj

that we generally tend to overweight the significance of energy’s

(EE)

Power output of the grid.

Absolute value for temperature change. Temperature for the output water or vapor. Initial investment/Operation cost. Capital recovery factor. Equipment’s initial investment for unit capacity in unit time.

in

Pj

Rated capacity for equipment j.

Particle Swarm Optimization (MOPSO) algorithm is applied to

r

Annual interest rate.

the distributed IES. By this method, the PF of the multi-objective

n

Life-cycle of the facilities.

thermodynamics. During operation simulation, a Multi-Objective

M runj /M wj Equipment’s fixed operation costs/ variable

model can be easily found.

operation costs for unit capacity in unit time.

Index Terms—multi-objective evaluation; Pareto Front (PF);

M gas , elec ,i

Exergy Efficiency (EE); distributed integrated energy system (IES)

the subsidy for unit capacity of the devices in unit

NOMENCLATURE

time.

Indices:

i

P k  Ck

Index for the evaluating time stage. Index for serial number

by the sequence of N O , SO2 ,

i

i

influence/ the greenhouse effect/ the dust

2

NOx , CH 4 , CO2 , PM 2.5 . j

emission effect of the gas k .

Index for serial number by the sequence of photovoltaic

K / M /W

power generations, wind power generations, the small

Efficiency/economy/environment objective.

Variables:

hydropower stations, electric storages, CCHPs, gas boilers,

Pl (i )

electric storages, heat pumps, heat storages, electric

Power output of the CCHP in unit time.

H l (i ) / Cl (i ) Heating/cooling output of the CCHP in unit time.

coolers, the distributed electric grids, the cooling grids and

Fl (i ) Qg i Fg (i)

the heating grids. Input Parameters and Functions:

Pl min / Pl max Minimum/maximum power output of the CCHP. Hl min / Hl max Minimum/maximum power output of the CCHP. 1

‹,(((

Emission factor/ amount of pollutant k .

w1 / w2 / w3 Weight coefficients. I / U / T i Factor denoting the atmospheric acidification

T Index for the evaluating time. k

Price of gas, the price of electricity in the grid or

Consume of gas in the CCHP in unit time. Heat output of the boiler in unit time. Consume of gas in the boiler in unit time.

Pdiselec,char/ heat i CEC (i )

Pwind (i)

Ppv (i)

Power input/output of the electric/heat storage in

Strength Pareto Evolutionary Algorithm (SPEA) algorithm [15]

unit time.

and some other Sequence Quadratic Programming (SQP) methods [16-17] are used to solve the multi-objective

Cooling output of electric cooler in unit time. Power output of wind power generation in unit

optimization problem by getting the Pareto Front (PF), and then

time.

analyze the results. That way, the PF can act as a perfect reference to evaluate the planning and operation of a system.

Power output of photovoltaic power generation in

However, it’s rather difficult to get an accurate PF when the

unit time.

Phd (i)

problems are sophisticated especially in distributed IES.

Power output of hydropower station in unit time.

Moreover, how to utilize the PF to conduct a multi-objective

I. INTRODUCTION

evaluation is seldom researched.

Under the circumstances of globalized energy crisis and the

With considerations above, in this paper, we propose a

environmental problem, many countries have decided to

method to realize the multi-objective evaluation for a

upgrade the energy system and to include the clean energy in

distributed IES, in which the mathematic model of distributed

the way of distributed Integrated Energy System (IES), which

IES is simplified, a Multi-Objective Particle Swarm

is also known as Energy Hub and Multiple energy carriers [1-

Optimization (MOPSO) algorithm is utilized to get the PF and

2].To establish a distributed IES which is optimal in terms of

the multi-objective potential variant is brought about to

economic, environment friendliness and renewable energy

evaluate the system based on the PF.

integration is one of the major concerns in many researches.

The remaining content is organized as below: In Section II,

However, an accurate way to evaluate the distributed IES is still

the model formulation of distributed IES is introduce. Then, a

under discussion because of the complexity of the distributed

mathematical model for multi-objective evaluation method is

IES and the coupling between different energy flows [3-4].

presented in Section III, as well as the multi-objective potential

Taking all these factors into consideration, this study put

variant, which is defined to illustrate how to evaluate the

forward a multi-objective evaluation plan for the disposition of

system based on the PF. In Section IV, a MOPSO algorithm is

distributed IES based on the analysis of the distributed IES in

applied to search for the PF of the model. Case studies are

Guizhou Province, China.

carried out in Section V to demonstrate the feasibility and

Many scholars have done a lot of research on the single

improvement of the multi-objective evaluation method.

objective evaluation of the CCHP [5-13] and the chosen

II. MODELING OF DISTRIBUTED IES.

objective are paramount to the evaluation. The objectives like

A. Combined Cooling Heating and Power (CCHP) System

the energy efficiency, the exergy efficiency [5] and the

This paper mainly focused on the optimization of ecology,

Coefficient of Performance (COP) [6] are utilized to evaluate

economy and environment in a distributed IES, thus the output

the interest of ecology of a distributed IES system. In [7-9], the

and the electric-thermal coupling relationship are mainly

interest of economy is mainly concerned and the emission of

concerned in the mathematic model.

the CO2, the influence of the greenhouse effect, the pollution

 Pl min d Pl i d Pl max 

of the emission, etc. [10] are utilized to evaluate the interest of

(1)

H l min d H l i d H l max

environment. However, concerning the complication of the

 aPl i  bH l i

distributed IES single objective analysis is far from enough.

(2)

Fl i 

Hl i

There are also some literatures on multi-objective

(3) c1 Pl i  c2 

evaluation for distributed IES [11-24]. Generally, there are two

(4)

ways. Firstly, the integrated evaluation [11], AHP-entropy

B. Gas Boiler, Storage, Electric cooler and Other Equipment

weight method [12], the hierarchal analysis [13], the fuzzy

When the power output of the CCHP is limited, the system

method [14], etc. are applied to transform multi-objective

may need the gas boiler to help fulfill the heat load. And in a

problem into single objective problem making the problem

distributed IES, the electric storage devices will charge and

much easier to solve. Nevertheless, it’s controversial to add

discharge for at least one time during one day. It can be charged

objectives with different properties together. Secondly, the

when the electricity price is low and discharged to meet the 2

need of the power load in system. The electric cooler’s cooling

The energy efficiency can be described as follow.

power is related to the power flow and the efficiency. The

K=

model for gas boiler, storage, electric cooler and other

Q6  Q7  Q8 Q2  Q3  Q4  Q11  Q5

(5)

equipment like photovoltaic power generation, wind power

In the calculation of the renewable energy, the formula leaves

generation and small hydropower station are illustrated in [11].

out the transformations from wind energy, solar energy and

III. MODELING OF MULTI-OBJECTIVE EVALUATION

hydro energy to power energy, for the transformation efficiency

This paper give a multi-objective model concerning the

is so low that the value will be too small to analysis. Besides,

interest of ecology, economy and environment. And the

Q1 , Q9 and Q10 , denoting power supply for electric storage,

evaluation are given based on the multi-objective potential

heat supply for the heat pump and heat supply for the heating

variant which is defined in this paper to illustrate to potential of

storage, are not included in the formula for their operational

the designing plan and operational plan in the system according

cycle time are far shorter than the evaluating time.

to the PF. This evaluation will be conducted every half year and

The exergy efficiency can be described as follow.

can bring about a revisionary plan.

K=

A. Model of ecology objective The distributed IES combines the power, heating and cooling

Q6 ˜ D 6  Q7 ˜ D 7  Q8  Q2  Q3  Q4  Q11  Q5 'T , i 6,7  Ti

Di

together, as is shown in Fig.1. It consists of several facilities,

(6) (7)

the model of which are illustrated in Section II. Q2 is equal to Fl . Q3 , Q4 and Q5 is respectively equal to Pwind , Ppv and Phd .

The electricity energy is of high quality, for which the

The green lines show the power flow of renewable energy, the

amend the plan and make high quality use of the energy.

orange lines show the heat flow, the blue lines show the cooling

B. Model of economy objective

conversion factor is 1. The exergy efficiency can be used to

flow, the grey lines show the other power flow and the dotted

The costs and subsidies of a distributed IES can be used for

lines show the energy supply inside the system, which are not

the economy quantitative analysis of a system. In this paper, the

included in Q6 , Q7 and Q8 .

model of economy objective consists of two parts, the investment cost and the operation cost.

M in  M run

M 1.

(8)

The initial investments The initial investment can be calculated by: 10

¦M

M in

inj



(9)

j 1

M inj

R

2.

R ˜ M unitj ˜ Pj

r ˜ (1  r ) n  (1  r ) n  1

(10) (11)

The operation cost The operation cost of a system consists of three parts:1) the

fixed operation costs include the fee to hire employee; 2) the variable operation costs such as the maintenance costs; 3) other

Fig.1 The energy flow in distributed IES

costs or subsidies like the fuel cost, the fee to buy electricity

Both the energy efficiency and the exergy efficiency can be

from the grid and the subsidies for renewable energy.

used to evaluate the interest of ecology in a system. Whereas,

10

the former focuses on evaluating the efficiency of energy

M run

10

runj

j 1

quantity and the later focuses on evaluating the efficiency of

˜ Pj  ¦ T ˜ M wj ˜ Pj + j 1

5

¦¦ M

energy quality. 1.

¦T ˜ M

i 1 T

The Energy Efficiency 3

gas ,elec ,i

˜ Fj ˜ Pj

(12)

And the constraints of load can be illustrated below including the load balance of power, heating and cooling.

C. objective function about environment impact

¦ T

This paper is mainly concerned about the environment pollution of several gas like N O , SO2 , NOx etc., which can bring about atmospheric acidification and CH 4 , CO2 , PM 2.5 etc.,

l

l

EC

6

T

( Fl (i )  Fg (i )) 

(13)



{

TABLE I. THE VALUE OF EMISSION FACTOR (g/kWh)

Pk

7

T

T

The multi-objective potential variant can be calculated by:

T

Pollutant

heat dis , char

g

T

amount of the pollutants can be calculated by: k

(17)

¦ ( H (i)  Q (i)  P (i)) ¦ Q (i) (18) ¦ (C (i)  C (i)) ¦ Q (i)  (19)

which can bring about the greenhouse effect. The emission

¦P

Q8 (i) ¦ T

elec  Phd (i)  Pdis , char (i ))

2

Ck

( Pl (i)  Q11 (i)  Pwind (i)  Ppv (i)

'x

,

'y

,

'z

xmax  xmin ymax  ymin z max  z min

}  (20)

CO2

CH 4

N 2O

SO2

NOx

PM 2.5

where x, y and z denotes three objectives. And 'x , 'y and

203.74

0.015

0.0004

0.011

0.202

0.0012

'z is the distance between this solution and the PF of whole

The influence of the pollutant can be divided by the influence

solution set in the direction of x, y and z, the formulas for which

of atmospheric acidification, the greenhouse effect and the dust

can be found.

emission effect [11], which can be illustrated by this formula:

W

6

6

6

i 1

i 1

i 1

w1 ¦ Ii ˜ Ci  w2 ¦ Ui ˜ Ci  w3 ¦Ti ˜ Ci

(14)

TABLE II. THE VALUE OF POLLUTION FACTOR Pollutant

U (g CO -

CO2 CH 4 N 2O SO2 NOx PM 2.5

equiv./g pollutant) 1 21 310 -

i

2

I

(g SO2 i equiv./g pollutant) 0.7 1 0.7 -

T i (g PM 2.5

equiv./g pollutant) 1.9 0.3 1

Fig.2 The multi-objective potential variant

D. Model of Multi-Objective and Multi-Objective Potential

The multi-objective potential variant can be utilized to

Variant 1.

evaluate the potential of the plan and find out the best way to

The Multi-Objective Model

improve the plan, which means if the potential variant on one

The multi-objective model including all the objectives above,

direction is larger, the plan should be adjusted on that interest.

which is a vector:

IV. SOLUTION STRATEGY

F {K , M ,W }

The algorithm to multi-objective evaluation includes the

(15)

following steps: Step1) Multi-Objective optimization:

This can be used to compare different plans of planning and

Calculate the optimal solution sets of the multi-objective

operation of a system. 2.

problem with (15). Step2) Pareto Front: Calculate the PF of

The Multi-Objective Potential Variant

the system by algorithms like MOPSO if the model is non-

The multi-objective potential variant is defined in this paper

convex; and many more simple algorithms can be used if the

to efficiently make use of the Pareto Front (PF) in the procedure

model is convex. Step3) Multi-Objective Potential Variant:

of evaluation, which can evaluate the potential waiting to be

Calculate the potential variant with (20). Step4) Exergy

improved in the plan.

Efficiency (EE): Calculate the EE to amend the plan if the plan

The PF can be found by solving this multi-objective

need to make high quality use of the energy. Step5) Evaluation

optimization problem. The objective function is: min{K ( x ), M ( x ), W ( x )}

Outcome: Evaluate and improve the plan based on F, σ and EE (16)

respectively gotten in step1), 3) and 4). The improvement can

The constraints consist of two parts, the constraints of load

be conducted by analysis the relationships between different

and the constraints of the equipment, illustrated in Section II.

decision variables and different objective functions. 4

can be thoroughly evaluated and the improvement direction of the plan are clearer.

V. CASE STUDIES

B. Evaluate and Improve the Energy Quality Efficiency

This paper use the distributed IES being constructed in

EE

Guizhou Province, China as a case.

In some circumstances, the system need to be planned and

Firstly, this paper evaluates the plan and brings up an

operated with high energy quality efficiency making it

improved plan without the amend of EE like in most real case

necessary to supervised the improved plan with EE. After

the energy efficiency is enough for evaluate ecology property,

amended, in the supervised plan E=0.9915 and F={0.8001,

thus making the evaluation method feasible and illustrating

3.900×107, 2.472×109}, while in the improved plan EE=0.9716

how to improve the plan based on the multi-objective potential

and F={0.883, 3.560×107, 2.094×109}. Apparently, to improve

variant. Then, another plan amended by EE is proposed to show

EE, the other objectives will be slightly lowered. Nevertheless,

how the EE criteria can be effectively utilized to evaluated and

the EE can be improved to almost 100%, which is of great

improve the energy quality efficiency. The initial plan is not

significance to energy conservation and emission reduction.

located on PF and the evaluation time ranges from May to Nov.

TABLE III. THE CONSTRUCTION OF EQUIPMENT AND LOADS IN

A. Evaluate and Improve the Plan without EE

THE SYSTEM

The construction of facilities and loads in the distributed IES is give in table III. The evaluation can be conducted by

Equipment &Loads

the algorithm in Section IV and as the model is non-convex,

CCHP1

the MOPSO is utilized to find PF efficiently. The evaluation

CCHP 2

outcome of the initial plan is F={0.8474, 7.330×108, 2.472×109} and σ={0.231,0.563,0.245}, which shows the initial plan can be improved by 23.1% in ecology objective, 56.3% in economy objective and 24.5% in environment objective. Based on the evaluation above, the economy interest of the initial has a huge potential to be improved. After improvement the evaluation outcome is F={0.883, 3.560×107, 2.094×109}

with

CCHP 3

Capacity

Remarks

1.2MW /0.6MW 1.6MW /0.8MW 2MW /1MW 4MW

The max power output is 1.2MW; the max heating/cooling output is 0.6MW The max power output is 1.6MW; the max heating/cooling output is 0.8MW The max power output is 2MW; the max heating/cooling output is 1MW Being Under planning

Wind Farm Small Hydropower 12MW Station Electric Storage 2MWh

Fixed capacity in the system

PV

4MW

Power Load Heating Load Cooling Load

8395MW 1080MW 55500MWh

and σ={0.012, 0.185,0.04}, which shows the improved plan is

The max storage capacity is 2MWh The max power output of the distributed PVs is 4MW The max power load is reached in summer The max heating load is reached in winter The max cooling load is reached in summer

VI. CONCLUSION

nearly on PF, which can be seen on Fig.3.

A multi-objective model and a multi-objective potential variant is brought up in this paper to evaluate the interest of ecology, economy and environment of a distributed IES. As is shown in case studies, this method based on the PF of the multievaluation model are more convincing and indisputable, which can evaluated the initial plan thoroughly and point out a clearer improvement direction of the plan. Moreover, to fulfill the need to make high quality use of the energy, this paper utilizes the EE to amend the evaluation, which is of great significance to energy conservation and emission reduction. Based on the proposed multi-objective model, several

Fig.3 The values of the decision variables And by comparing the values of σ, the potential of the system

interesting directions are open for future study. How to find the

has be dug out by 21.9% in ecology objective, 37.8% in

PF of highly complicated distributed IES model more

economy objective and 20.5% in environment objective.

effectively with guaranteed convergence is still a note-worthy

Apparently, with the multi-objective model and the multi-

topic. The analysis of the accuracy of the PF and the

objective potential variant model in this paper, the initial plan

quantitative analysis of the relationships between different

5

decision variables and different objectives are worth further

[11] LI Miao, Integrated Evaluation of Distributed Energy System Based on 3E

investigating.

Benefit, Dalian: Dalian University of Technology, 2015. [12]

ACKNODEGMENT

DONG

Fugui,

ZHANG

Ye,

SHANG

Meimei.

Multi-criteria

Comprehensive Evaluation of Distributed Energy System [J]. Proceedings of the CSEE,2016,36(12):3214-3222.

The authors gratefully acknowledge the National Natural Science Foundation of China (51537006) and the Shenzhen

[13] HUANG-FU Yi, WU Jing-Yi, WANG Ru-Zhu, HUANG Xing-Hua. Study

Municipal Development and Reform Commission, Shenzhen

on Comprehensive Valuation Model for Combined Cooling Heating and Power

Environmental Science and New Energy Technology

System(CCHP)[J].Journal of Engineering Thermophysics, 2005,26(9):13-16.

Engineering Laboratory, Grant Number: SDRC [2016]172

[14] Youyin Jiang, He Bai, Jiangjiang Wang. A fuzzy multi-criteria decision

that provided the funding.

making model for CCHP systems driven by different energy sources[J]. Energy policy, 2012(42): 286-296.

REFERENCES

[15] WANG Yi, CHENG Hao-zhong. Pareto Optimality Based Multi-objective

[1] M. Geidl, G. Koeppel, P. Favre-Perrod, B. Kliockl, G. Andersson and K.

Transmission Planning Considering Transmission Congestion [J]. Proceedings

Frohlich,“Energy hubs for the future,” IEEE power and energy magazine, vol.

of the CSEEˈ2008ˈ28(13):132-138.

5, no. 1, pp. 24–30, Jan. 2007.

[16] ZHENG Zhang-hua, AI Qian, GU Cheng-hong, JIANG Chuan-wen.

[2] M. Geidl and G. Andersson,“Optimal power flow of multiple energy

Multi-objective

carriers,” IEEE Transactions on Power Systems, vol. 22, no. 1, pp. 145–155,

Allocation

of

Distributed

Generation

Considering

Environmental Factor [J]. Proceedings of the CSEE,2009,29(13):23-28.

Feb. 2007.

[17] M. H. Nehrir,C. Wang,K. Strunz, et al. A review of hybrid

[3] H. Sun, Q. Guo and Z. Pan,“Energy Internet: Concept, Architecture and

renewable/alternative energy systems for electric power generation:

Frontier Outlook,” Automation of Electric Power Systems, vol. 39, no. 10, pp.

configurations, control, and applications [J]. IEEE Transactions on Sustainable

1–8, Oct. 2015.

Energy, 2011, 2(4): 392-403.

[4] J. S, Z. Ma, Y. Shang and D. Liu,“Technical Review of Round Table Session of CICED 2016,” Power System Technology, vol. 40, no. 11, pp. 3368-3374, 2016. [5] Y. Xue, Q. Guo, H. Sun, X. Shen and L. Tang, “Comprehensive energy utilization rate for park-level integrated energy system”, Electric Power Automation Equipment, vol.37, no. 6, pp.117-123, 2017,. [6] De Angelis Francesco and Grasselli Umberto, “The next generation green data center: A multi-objective energetic analysis for a traditional and CCHP cooling system assessment”, IEEE Conference Proceedings, pp.1-6, Jun. 2016. [7] X. Shen, M. Shahidehpour, Y. Han, S. Zhu and J. Zheng, “Expansion planning of active distribution networks with centralized and distributed energy storage systems,” IEEE Transactions on Sustainable Energy, vol. 8, no. 1, pp. 126-134, Jan. 2017. [8] X. Shen, S. Zhu, J. Zheng, et al, “Applying phase-change energy storage in active distribution system planning,” CIRED Workshop 2016, pp.186, Jun. 2016. [9] X. Shen, Y. Han, S. Zhu, et al, “Comprehensive power-supply planning for active distribution system considering cooling, heating and power load balance,” Journal of Modern Power Systems and Clean Energy, vol. 3, no. 4, pp. 485493, 2015. [10] Y. Jing, H. Bai and J. Wang, “Multi-objective optimization design and operation strategy analysis of BCHP system based on life cycle assessment,” Energy, vol.37, pp. 405-416, Dec. 2012.

6