Hierarchical Multi-Objective Fuzzy Collaborative Optimization of Integrated Energy System under Off-Design Performance

: In order to solve the frequently occurred nonlinear working characteristics problem of the devices in integrated energy system (IES), where this nonlinear problem causes the offset of operation characteristics and design points of those devices, in this paper, the hierarchical multi-objective fuzzy collaborative optimization model of IES under off-design performance is proposed. Firstly, the high-order nonlinear models of devices under off-design performance, including fuel cell (FC), gas turbine (GT), absorption lithium bromide chiller (ABS), etc., are presented considering renewable energy and diversiﬁed storage devices, which can more accurately describe the actual working condition of devices. Secondly, according to the needs of different engineering application modes, the lower level collaborative optimization model of energy hub (EH) with economic, eco-friendly, primary energy saving, and renewable energy accommodation rate as optimization objectives is proposed, which can well adapt to and reﬂect real world energy system. Thirdly, based on the relationship between master and slave dispatch centers, the upper level modeling and optimal scheduling model of IES consisting of coupling power system, natural gas system, and heat system is proposed. Then, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, where continuous differentiable Sigmoid function is taken as the membership function. Finally, simulation results show that the proposed models and optimal dispatch method can effectively solve the frequently occurred non-linear working characteristics problem of the devices in IES. The optimal results are suitable for planning, calculation, operation, dispatch of IES, which can reasonably reﬂect the operation characteristics of IES.


Introduction
The integrated energy system (IES) has become a vital strategic research direction in international energy field, due to its great significance in improving the comprehensive utilization efficiency of energy, the utilization efficiency of social infrastructure, and the energy supply security [1][2][3].In addition, the IES also can promote the large-scale development of renewable energy and achieving energy conservation and emissions reduction in China [2,3].At the same time, the remarkable change of energy production and consumption mode requires changing the construction path and development mode of traditional energy system and constructing IES [4][5][6].Constructing IES from economic, management and planning perspective is an important part of "building a low-carbon, clean, safe and efficient modern energy system" [6,7].
Energies 2019, 12, 830 3 of 27 (1) High-order nonlinear models of devices under off-design performance are developed considering renewable energy and diversified storage devices.(2) EH under off-design performance is introduced based on geographical resource endowment, then, a lower level multi-objective collaborative optimization model for such case is proposed.
(3) Upper level modeling and multiple optimal scheduling modes of IES with the coupling and interactions between heat system, power system, and natural gas system, are proposed, fully considering the key characteristic variables of the energy subsystems.(4) A hierarchical multi-objective collaborative optimization model for EH and IES is developed, facing engineering application modes and based on the relationship between master and slave dispatch centers.( 5) An improved multi-objective fuzzy optimization method is presented, where continuous differentiable Sigmoid function is introduced as an improved membership function.

Lower Level Modeling and Optimization
EH is defined as an input-output port model that describes the exchange and coupling relationship between source, load, storage, and network in IES.The EHs are a sub-part of the IES, where the EH is an equivalent multi-energy demand load terminal of IES.EH can act as a mutual conversion, distribution, consumption, and storage of all kinds of energy, with a strong flexibility [6,8,10].

Energy Device Unit Model of Energy Hub under Off-Design Performance
The main frame model of EH is relatively mature, and the detailed understanding can be referred to [10,[14][15][16].Next, this paper mainly introduces the key devices model under off-design performance in EH.

Fuel Cell (FC)
FC is a kind of power generation and heat generation device that directly converts the chemical energy of fuel into electric energy.That chemical energy of hydrocarbon fuel such as natural gas and hydrogen is directly converted into electric energy through electrochemical reaction, where its load response is fast and the operation quality is high.FC can also achieve good coupling and complementarity with GT and power generation system, and form a multi-energy flow supply system.It follows from [12,17] that the FC model under off-design performance is shown as: where η t FC and γ t FC denote the generation efficiency and heat-to-electric ratio of FC at time t, respectively.ϕ η (•) is the mathematical functional relation.a n and b n are the constant coefficients of the nth power item in different functions.F t FC , P t FC and P t FC represent the input fuel amount, output electric power and part load rate of FC at time t, respectively.X t FC is the binary-variable of FC operation state, with 1 for run and 0 for shutdown.
The established FC model above has extensive universality.A typical FC case under off-design performance is shown in Figure 1, in which the data is from the reference [17].As can be seen in Figure 1, the generation efficiency and heat efficiency of FC is changing under different part load rates.Obviously, a large error will be produced by using the modeling method of the constant generation efficiency without considering the off-design performance characteristics.

Gas Turbine (GT)
Regarding GT as the main energy supply equipment of the multi-energy flow system with coupling natural gas, electricity and heat, its operation efficiency varies greatly under part load rates [10,14].The third-order generation efficiency function is used to fit the model [18], the GT model under variable operating conditions is given by: ( ) ( ) ( )  REC centralizes recovery of waste heat generated by GT and FC to achieve further EH energy utilization efficiency [18,19].Common REC is mainly composed of regenerators, or absorption refrigeration equipment, or heat exchangers, etc.The heat energy output by REC is related to the heat recovery efficiency, which is also affected by the part load rate, which can be expressed by: ( ) ( )

Gas Turbine (GT)
Regarding GT as the main energy supply equipment of the multi-energy flow system with coupling natural gas, electricity and heat, its operation efficiency varies greatly under part load rates [10,14].The third-order generation efficiency function is used to fit the model [18], the GT model under variable operating conditions is given by: where F t GT , P t GT , P t GT and H t GT are the input fuel amount, output electric power, part load rate and output thermal power of GT at time t, respectively.η t GT and η loss represent the generation efficiency and heat loss rate of GT at time t, respectively.c n is the constant coefficient of the nth power item for GT model.X t GT denotes the binary-variable of GT operation state, with 1 for run and 0 for shutdown.

Waste Heat Recovery System (REC)
REC centralizes recovery of waste heat generated by GT and FC to achieve further EH energy utilization efficiency [18,19].Common REC is mainly composed of regenerators, or absorption refrigeration equipment, or heat exchangers, etc.The heat energy output by REC is related to the heat recovery efficiency, which is also affected by the part load rate, which can be expressed by: where F t REC , H t REC and H t REC are the input waste heat power, output heat power, and part load rate of REC at time t, respectively.η t REC is the heat recovery efficiency of REC at time t.d n represents the constant coefficient of the nth power item for REC model.X t REC denotes the binary-variable of REC operation state, with 1 for run and 0 for shutdown.

Gas Boiler (GB)
GB is usually used as the heat source of EH to provide or supplement heat energy, and its working characteristics are relatively stable.The heat generated by GB is related to the heating efficiency of Energies 2019, 12, 830 5 of 27 boiler, it follows from [18][19][20] that the second-order variable working condition characteristic model of GB is formulated as: where F t GB , H t GB and H t FC are the input waste heat power or fuel amount, output heat power, and part load rate of GB at time t, respectively.η t GB represents heating efficiency of GB at time t. e n is the constant coefficient of the nth power item for GB model.X t GB denotes the binary-variable of GB operation state, with 1 for run and 0 for shutdown.

Ground Source Heat Pump (GSHP)
GSHP usually takes shallow geothermal resources as the heat source, which has the advantages of energy saving, environmental protection and convenient regulation.It can be used for heating in winter and cooling in summer [10,18].Compared to traditional heating methods that can only convert a small part of the electrical energy into heat energy, which is not economic, GSHP's output (thermoelectric or cold energy), by absorbing low temperature heat from the environment or releasing heat energy into the environment, can be relatively high with only consuming small amount of electric power as an input, thereby significantly reducing fossil energy consumption [18,21].However, its working time pattern is usually limited for reducing damage to the energy balance of the environment, thus GSHP works as follows: where P t GSHP , H t GSHP , and H t GSHP are the input electrical power, output heat power, and part load rate of GSHP at time t, respectively.COP t GSHP denotes coefficient of performance of GSHP at time t.T tem ev is the average water temperature at the inlet and outlet of evaporator.T tem cd denotes the average water temperature at the inlet and outlet of the condenser.T tem ev and T tem cd are affected by seasonal climate and environmental temperature, therefore, for the sake of simplicity, and this paper does not model temperature changes and sets them as given the variables.f 0 ∼ f 2 are the temperature fitting constant coefficients.g n denotes the constant coefficient of the n-th power item for GSHP model.X t GSHP is the binary-variable of GSHP operation state, with 1 for run and 0 for shutdown.T start and T end are the starting and ending time of GSHP's working period, respectively.ε acc denotes the ratio of the maximum working time of GSHP in the optimized scheduling cycle T.

Micro-Gas Turbine (MT)
MT can use natural gas, biogas and other fuels for power generation and heating.Generally, the generating efficiency of MT is improved with the increase of power generation, and reaches the maximum efficiency at the rated load [17].The load characteristics of MT describe the relationship between the equipment efficiency, smoke discharge temperature, smoke discharge flow rate and other parameters when the output power changes at a constant speed.It follows from [18,22] that the relationship can be expressed by the following equations in the case of capacity determination: where F t MT , P t MT and P t MT denote the input fuel amount, output electrical power, and part electric load rate of MT at time t, respectively.η t MT represents electrical efficiency of MT at time t, respectively.T tem,t MT denotes the reduced smoke exhaust temperature, that is, the ratio between the actual smoke exhaust temperature of MT and its calibration value.m tem,t MT is the reduced flue gas flow rate, that is, the ratio between the actual flue gas flow rate of MT and its calibration value.h 0 ∼ h 3 are the efficiency fitting constant coefficients.j n and k n stand for the constant coefficient of the n-th power item for smoke exhaust temperature, and flue gas flow rate model.

Absorption Lithium Bromide Chiller (ABS)
The high temperature flue gas discharged from MT is available to be used for refrigeration, heating and supply of domestic hot water through the flue gas ABS, so as to realize the cascade utilization of different quality energy [12,18].When the output thermal power of MT changes, the smoke exhaust temperature and gas flow rate also change, which lead to the change of the cooling or heat production for ABS [23].The functional relationship between the output thermal power of ABS and flue gas flow rate, and flue gas temperature can be described by: where H t ABS denotes part heat load rate of ABS at time t.T tem,t ABS denotes the reduced temperature, that is, the ratio of ABS smoke inlet temperature to its calibration value.m tem,t ABS stands for the reduced flue gas flow rate, that is, the ratio of ABS gas flow rate to its calibration value.m 0 ∼ m 5 are the efficiency fitting constant coefficients.
2.1.8.Electric Boiler/Power to Gas/Transformer/Heat Exchanger (EB/P2G/TRANS/HE) EB, P2G, TRANS, and HE are less affected by changes in the external environment and part load rate, and their off-design performance are ignored in this paper [8,10].It follows from [24,25] that the relationship between their input and output can be uniformly expressed as: where F in,t type , H out,t type , and η type represent the input power, output power and energy conversion efficiency of type, respectively.2.1.9.Electricity/Heat/Natural Gas Storage (ES/HS/GS) The energy storage unit realizes energy transfer across time period and coordinates the imbalance between "source and load" in EH, which has the function of cutting peak and filling valley, and restraining the fluctuation of renewable energy, hence, promotes IES to operate economically and efficiently [19,22].In EH, electricity storage (ES), heat storage (HS) and natural gas storage (GS), three different types of energy storage equipment, are used to ensure the flexibility of EH's optimized operation strategy.Meanwhile, the time mismatch between heat load and power load, and wind power and photovoltaic is shifted [24,25].ES, HS and GS models are similar, and their general models can be summarized as: where P cha,t type , and P dis,t type are the charging, and discharging power values of type at time t, respectively.W t type denotes the energy storage capacity of type at time t.η cha type and η dis type represent the working efficiency of charging and discharging of type, respectively.σ type is the energy loss rate or self-released energy rate of type during the course of charging and discharging energy.∆t stands for the simulation step size of EH optimal scheduling.2.1.10.Wind Turbine/Photovoltaic Cell/Solar Collector (WT/PV/SC) Due to the randomness and volatility of wind and solar energy, they are often not fully absorbed.In order to fully exploit renewable energy, assuming that both WT and PV work in the maximum power point tracking mode, a prediction curve of renewable energy can be used to represent the maximum output in a day.In order to make EH optimal scheduling policy more flexible, the optimization strategy is allowed to determine whether WT and PV are put into operation or not, as well as the degree of absorption.In this paper, the WT model [24] and PV model [19] are utilized, which can be expressed by: where P t wt denotes the WT predicted maximum output power at time t.P r is the rated output power of WT. v(t) represents the actual predicted wind speed at time t.v in , v r and v out are the cut in wind speed, rated wind speed and cut out wind speed, respectively.x, y, and z stand for a polynomial constant coefficients: where P t pv and G t pv denote the PV predicted output maximum power, and the actual solar radiation intensity at time t, respectively.P stc and G stc are the output power, and the solar radiation intensity under standard test conditions, respectively.δ pv stands for the power drop factor of PV, which is used to characterize the PV surface grime and aging as well as other factors that cause the decline in output power, and it is generally valued at 0.9.k pv is the power-temperature coefficient.T tem,t pv denotes actual surface temperature of PV at time t.T re f represents the temperature for PV under standard test conditions.
Solar collector (SC) is a kind of green device which can convert solar energy into heat energy by using only a small amount of electric energy, and it can effectively use solar energy to make heat, and has the advantages of cleanliness, environmental protection and pollution-free [5,6,26].The model of SC can be described by: where P t sc denotes the predicted maximum output thermal power of SC at time t.η t sc is the photo-thermal conversion efficiency of SC at time t.A sc denotes the effective heat collecting area of SC.

I t
sc is the solar radiation value at time t.η 0 sc denotes the standard photo-thermal conversion efficiency or heat collection efficiency of SC. k sc is the temperature difference coefficient.T tem,t sc,a stands for the temperature of SC under actual operation at time t.T tem,t sc,o denotes the operating ambient temperature of SC at time t.
Remark 1.In this paper, considering nonlinear working characteristics problem of the devices, the high order off-design performance model of key devices is proposed, and the utilization of waste heat and environmental factor are fully considered.Compared with the results in [4,8,10,22], the newly developed model can better reflect the actual operation characteristics of device in engineering applications.At the same time, the interaction and coupling between output and input for devices are also highlighted in these models.

Multi-Objective Collaborative Optimization Model of Energy Hub
As is well known that optimization problem is one of the main problems in engineering applications and scientific research, in which the multi-objective optimization problem is very commonly issue, and the objective functions involved are usually not comparable (conflicting with each other in most cases) [6,22].Especially in EH, economic dispatch, energy saving and emission reduction, and improvement of renewable energy permeability are challenging multi-objective optimization problems [11,20].Therefore, it is necessary to develop a multi-objective operation optimization model considering the operating cost, consumption of renewable energy, pollution emission and consumption of fossil energy in EH [12,16,27].In this paper, four optimization objective functions are presented, respectively.

Economic Optimization Objective
The operation cost of EH in the dispatch cycle is usually concerned to achieve the optimal operation cost.This optimization function aims to minimize the operation cost of EH: where λ t g , λ t e , and λ t h represent the unit power price of natural gas, electric energy and thermal energy purchased by EH at time t, respectively.F t g , P t e , and H t h are the amount of natural gas injected into EH from the natural gas network, the electric energy injected into EH from the electric network, and the heat energy injected into EH from the thermal network at time t, respectively.T stands for optimized scheduling cycle of EH.

Environmental Protection Optimization Objective
It is an immediate priority for many entrepreneurs to effectively improve the efficiency of energy utilization, as well as environmental protection.Furthermore, the following optimization objective function is mainly to minimize EH emissions: where λ i type denotes pollution gas emission factor i = CO 2 , SO 2 , NO x of gas equipment ype = FC, GT, GB, MT at time t.κ i e , and κ i h represent the equivalent pollutant gas emission factor i that consumes unit electric energy, and heat energy, respectively.F t type represents the amount of fuel consumed by gas equipment type at time t.

Primary Energy Optimization Objective
The minimum primary energy (fossil energy) consumption is selected as the energy objective, which can be expressed as: where η g,e , and η t,e represent the average generation efficiency of thermal power generating unit, and the transmission efficiency of power grid, respectively.
where P a,t wt , P a,t pv , and P a,t sc represent the actual power injected into EH by WT, PV and SC at time t, respectively.

Optimization Constraints for Energy Hub
In order to ensure the safety and stability of EH optimal operation, relevant constraints, such as equipment unit and energy balance, need to be comprehensively considered [7,10].

Energy Bus Equilibrium Constraints
For energy bus type including electric bus, heat bus and natural gas bus, the sum of energy ∑ P in,t b,type flowing into the bus at time t is equal to the sum of energy P out,t b,type flowing out of the bus, and it can be summarized as follows:

Device Unit Operation Constraints
For device type, including FC, GT, GB, REC, P2G, E2H, HX, MT, ABS, and GSHP, we mainly consider their climbing rate constraints, minimum and maximum load rate constraint and maximum startup times N max type constraint within the optimization scheduling cycle, which can be expressed by: where δ o f f type denotes the operation excision coefficient of device type.For energy storage unit, in addition to physical constraints on charging power, discharging power and energy storage level, in order to ensure that energy storage unit can flexibly participate in the optimal operation scheduling of the next cycle, it is necessary to restrict the state of energy storage level at the beginning and at the end, which is given as: where ζ cha E/H , and ζ dis E/H are the maximum charging and discharging power factors, respectively.

Tie-Line Power Constraints
EH buys electricity from power system and natural gas from the natural gas system, which is constrained by the maximum input power rate in unit time.In this paper, WT, PV, and SC adopt the maximum power point tracking mode, and their actual power inputs of EH cannot exceed the maximum predicted power value.Those constraints can be summarized as follows: In general, the tie-line power is a variable to be optimized, and its minimum value can be equal to zero.According to the actual different engineering application requirements, different EHs characteristics and other factors, the minimum value of tie-line power may be set to a non-zero constant, but the set values in this paper are 0, that is, F min g = 0, P min e = 0 and H min h =0.

Upper Level Modeling and Optimization
Different from the modeling and optimization of a single energy supply system, IES includes many different forms of energy systems, such as electricity, natural gas, heat, and so on [4,10].The characteristics of various forms of energy system are also different, and there are complex coupling relationships among various forms of energy system at the same time, thus, modeling and optimization of IES is more complex [6,21].

Subsystems Model of Integrated Energy System
Next, power flow calculation model of multi-energy flow for IES will be developed, which mainly includes power system, natural gas system, heat system and their key equipment components.

Power System
Similar to the traditional power flow model of power system, power system in IES is described by the classical AC power flow model [3,10], and the power balance equation of its nodes is expressed as follows: where P and Q denote the active power and reactive power of nodes, respectively.Y is the node admittance matrix.
. U represents the vector of nodes' voltage.

Natural Gas System
Natural gas system delivers energy to the load through natural gas pipelines, and its source is usually natural gas well or natural gas storage device, while the load is usually a variety of gas equipment [1,10].The gas flow in pipeline is mainly controlled by adjusting the pressure of pressure regulating valves at different nodes.Next, it follows from [2,22,24] that, the steady state model of natural gas system is established.
For the ideal adiabatic natural gas pipeline, its start and end nodes are i and j, respectively, and its steady state pipeline flow q ij is described by: where p i and p j represent the pressure of node i and node j, respectively.k ij denotes a pipeline characteristic constant related to pipeline diameter, pipeline length, pipeline transmission efficiency, natural gas temperature, gas compression factor, etc., and the specific value can be found in [10,24,28].
π ij is the variable reflecting the direction of natural gas flow in the pipeline.
Due to the existence of friction resistance in the pipeline, transmission loss will be generated, which will lead to a certain pressure loss [16,29].In order to ensure the reliability of long-distance energy transmission for natural gas pipeline network, a certain number of compressors should be installed in the natural gas system to compensate the pressure loss caused by friction.
Figure 2 shows schematic diagram of the pipeline with compressor driven by gas turbine.The mathematical model of natural gas pipeline with compressor can be expressed as follows: where f jk denotes the natural gas flow consumed by the compressor.k com is the compressor characteristic constant.q jk denotes the natural gas flow through the compressor whose inlet node is j and outlet node is k.r jk = p k /p j represents compression ratio of compressor.α denotes compressor variable index.According to the law of flow conservation [22,30], the flow balance equation of each node in the natural gas network can be expressed in matrix form: where c A represents the node-branch complete incidence matrix of natural gas network.s A is the starting nodal association matrix of the natural gas network.q  denotes the flow vector of natural gas pipeline or branch.0 q  represents the node net load flow vector in natural gas system.f  denotes the flow vector consumed by compressor in natural gas pipeline.
After developing the natural gas pipeline model and pipeline with compressor driven by gas turbine model, the directed graph is adopted to represent the topology structure of natural gas pipeline network based on the idea of graph theory, and the branch reference direction is set in advance.Assuming that the compressors are located at the first end of the pipeline and that the compressors do not affect the pipeline length, the matrix operation model of the natural gas system is formulated as: where e A is the end nodal association matrix.
denotes the vector associated with the node pressure.
represents the vector associated with branch gas flow.
( ) , , , b diag r r r = Diag  denotes the diagonal matrix for the compression ratio.
The subscript b is the branch or pipeline number.The subscript node is the node number.′ ⋅ According to the law of flow conservation [22,30], the flow balance equation of each node in the natural gas network can be expressed in matrix form: where A c represents the node-branch complete incidence matrix of natural gas network.A s is the starting nodal association matrix of the natural gas network.
. q denotes the flow vector of natural gas pipeline or branch.
. q 0 represents the node net load flow vector in natural gas system.
. f denotes the flow vector consumed by compressor in natural gas pipeline.
After developing the natural gas pipeline model and pipeline with compressor driven by gas turbine model, the directed graph is adopted to represent the topology structure of natural gas pipeline network based on the idea of graph theory, and the branch reference direction is set in advance.Assuming that the compressors are located at the first end of the pipeline and that the compressors do not affect the pipeline length, the matrix operation model of the natural gas system is formulated as: Energies 2019, 12, 830 where A e is the end nodal association matrix.
node denotes the vector associated with the node pressure. .
b represents the vector associated with branch gas flow.Diag = diag r 2 1 , r 2 2 , • • • , r 2 b denotes the diagonal matrix for the compression ratio.The subscript b is the branch or pipeline number.The subscript node is the node number.(•) and [•] denote the transpose of the matrix.

Heat System
Heat system is usually composed of heat source, heat network and heat load, which assume the functions of heat energy production, transmission and exchange.The heat network includes a supply water pipeline network and a return water pipeline network with water or steam as the heat transfer medium, which connects the heat source with the heat load [22].Taking the supply water pipeline network as an example, the topology of heat network is described in this paper based on graph theory and the mature network topology description method in power system analysis [23,31].The established heat system model is suitable for IES analysis, by combining the analysis method of hydraulic balance and thermal balance process in traditional thermal network.
Heat transfer medium flow in the network should meet the basic law of network: each nodal flow should meet the flow continuity equation [10,32,33].In a closed loop of pipelines, the sum of the pressure head losses of heat transfer medium flow in each pipe is 0, which can be described by: where A h is the node-branch complete incidence matrix of heat system.B denotes the basic loop vector of heat system.
. m denotes the pipeline or branch flow vector.
. m 0 represents the net load flow of each node in heat system.∆ .H b is the pressure head losses vector of pipeline or branch in heat system, and the calculation method can be found in [10,28,33].
. H p represents the head vector of circulating pump in pipeline or branch in heat system.
The variables involved in the thermal model are nodal supply water temperature T tem s , nodal return water temperature T tem r and nodal heat power [10,22], and heat power includes heat output ∆h l of the heat source and heat demand ∆h s of the load, which can be described by: where C p is the specific heat capacity of heat transfer medium.m s 0 and m l 0 are the flow through the heat source and the heat load, respectively.
Considering the heat energy loss of the pipeline, the temperature drop of the heat transfer medium in its transmission is shown as: where T tem start,p and T tem end,p are the heat transfer medium temperature at the starting point of the pipeline, and the heat transfer medium temperature at the ending point of the pipeline, respectively.T tem a is the ambient temperature.m p denotes the flow rate of the pipeline, where the pipeline length is L p and the heat transfer coefficient is λ p .
For any node in the heat network, the sum of the incoming thermal power is equal to the sum of the outgoing thermal power [3,10], so the temperature after convergence can be calculated by: Energies 2019, 12, 830 13 of 27 where m out , and T out are the flow rate and temperature at the node outlet, respectively.m in and T in are the flow rate and temperature at the node entrance, respectively.
The heat network usually transmits heat energy to the heat load through the heat exchanger which is the equivalent heat load [33,34].Another type of heat load is heating buildings, such as residential areas, huge shopping malls, large libraries, etc., which are here called 'buildings'.This kind of heat load has strict requirements on temperature and needs to be modeled separately [21,25].The following equivalent thermal parameter model can be employed to describe the characteristics of buildings, where the equation of state is discretized during programming: where T in , T m,in and T out are indoor temperature, indoor air temperature, and outdoor environment temperature, respectively.R 1 and R 2 are the gain or heat loss coefficient.C a and C m are the heat capacity of air, and the mass heat capacity of buildings and their content, respectively.H l denotes the heating power of heating ventilation air conditioning system.u is the first order square matrix, and the element value is 1.

Remark 2.
Compared with [6,13,20], the developed IES model simultaneously coupling electric system, heat system and natural gas system can fully consider the key characteristic variables of each system, which really reflects "system level" operational characteristics.Furthermore, the matrix operation models about systems are also deduced firstly, as shown in ( 22)- (24), which can greatly facilitate the engineering calculation and application of large systems.
Remark 3. On the basis of heat exchanger as heat load, an equivalent thermal parameter model of heat load is proposed, compared with the models developed in [10,23,31], the proposed model can better reflect the continuity and correlation of heat load in time, and it has more depth and accuracy to measure the characteristics of heat load by temperature, thus which can improve the flexible adjustment ability of heat load, and better reflects the heat load characteristic.

Multi-Objective Collaborative Optimization Model of Integrated Energy System
Through collaborative optimization of various energy and utilizing the coupling mechanism between multi-energy systems in time and space, IES reduces the utilization of fossil energy, and improves the utilization rate of renewable energy, thus reducing operating cost and improving the comprehensive utilization level of energy [26,30,35].Therefore, an optimization model will be developed considering economic cost and renewable energy accommodation in the next part.

Economic Optimization Objective
As a complex energy system, IES has a variety of evaluation criteria.Under the background of current energy shortage, most analysis is based on reducing operating costs, in order to improve user economy.The specific optimization objective is given by: Energies 2019, 12, 830 14 of 27 where F t n,s,g denotes the total amount of natural gas purchased by IES from natural gas source n at time t, and the price of natural gas is λ t n,s,g .P t s,e denotes the electricity purchased by IES from the large power grid at time t, and time-of-use electricity price (TOU) is expressed in λ t s,e .

Renewable Energy Accommodation Rate Optimization Objective
Improving the contribution of renewable energy and promoting the clean energy of using are also of vital importance to the rapid development of IES.Therefore, wind power absorption is taken as the optimization objective of IES.The objective function is given by: where P s,a,t wt and P s,t wt denote the actual power injected into IES by WT, and the WT predicted output maximum power in the period t, respectively.

Optimization Constraints for Integrated Energy System
Power system constraints include tie line power exchange P t s,e constraints, CHP unit active power P t CHP , reactive power Q t CHP and heat power H t CHP upper and lower bound constraints, node voltage V t node upper and lower bound constraints and branch capacity S t b constraints, which are formulated as: where α k , β k and γ k are the kth inequality constraint coefficients of CHP operating region.
Natural gas system constraints mainly include gas source injection power F t n,s,g constraint, upper and lower limits of node pressure p t i , compressor compression ratio r t b constraints, which can be summarized as: In heating system, in order to ensure the high heating quality of the heat source and an efficient energy experience for terminal heat users, supply water temperature T tem,t s,l of the heat source, return water temperature T tem,t s,l of the heat load and flow rate m t p in pipeline need to be limited.In addition, it is necessary to ensure that the node pressure head h t b in heat network changes within the set range.The pressure head h t p lifted by the circulation pump should be also limited to the maximum head range.It can be described by: Terminal heat users or heat loads have a certain degree of fuzziness to the perception of comfortable temperature.In particular, for the buildings, in order to ensure the convenience of living, it is assumed that the flexible adjustment ability of buildings meets the following requirements:

Solution Approach
The lower level multi-objective optimization model for EH and the upper level multi-objective optimization model for IES under off-design performance established above are typical 0-1 mixed integer non-linear optimization problems [10,30], which can be solved by using the method of 0-1 mixed integer non-linear programming (MINLP).In this paper, the solution approach to solve problems of lower and upper level, single and multi-objective optimization, and optimization process, can be clearly demonstrated and illustrated in Figure 3.

Solution Approach
The lower level multi-objective optimization model for EH and the upper level multi-objective optimization model for IES under off-design performance established above are typical 0-1 mixed integer non-linear optimization problems [10,30], which can be solved by using the method of 0-1 mixed integer non-linear programming (MINLP).In this paper, the solution approach to solve problems of lower and upper level, single and multi-objective optimization, and optimization process, can be clearly demonstrated and illustrated in Figure 3. From Figure 3, optimization of each level includes single objective optimization and multiobjective optimization in this paper.Next, the related processes and concepts as shown in Figure 3 will be further explained.

Fuzzification of Objective Function
It is generally known that there are many methods to solve multi-objective optimization problem, and fuzzy optimization is one of the common methods.It is also universal to use fuzzy optimization to solve multi-objective optimization problems in academic papers, and this method is also employed in this paper.
Traditional fuzzy algorithm for solving multi-objective problems generally adopts increasing linear function or decreasing linear function as shown in Figure 4 ( n κ and n Δ denote the characteristic parameter of increasing linear function and decreasing linear function).In this paper, the following continuous differentiable Sigmoid function or reverse Sigmoid function are used as membership function, the parameters of Sigmoid membership function are also determined.Because of applying improved membership functions, nonlinear programming methods are suitable to solve above multi-objective optimization model: ( ) From Figure 3, optimization of each level includes single objective optimization and multi-objective optimization in this paper.Next, the related processes and concepts as shown in Figure 3 will be further explained.

Fuzzification of Objective Function
It is generally known that there are many methods to solve multi-objective optimization problem, and fuzzy optimization is one of the common methods.It is also universal to use fuzzy optimization to solve multi-objective optimization problems in academic papers, and this method is also employed in this paper.
Traditional fuzzy algorithm for solving multi-objective problems generally adopts increasing linear function or decreasing linear function as shown in Figure 4 (κ n and ∆ n denote the characteristic parameter of increasing linear function and decreasing linear function).In this paper, the following continuous differentiable Sigmoid function or reverse Sigmoid function are used as membership function, the parameters of Sigmoid membership function are also determined.Because of applying improved membership functions, nonlinear programming methods are suitable to solve above multi-objective optimization model: where ω 1 and ω 2 denote characteristic parameter of Sigmoid and reverse Sigmoid function.The Sigmoid function curve has the same semi-open shape as increasing linear function curve, so it is suitable for the membership function described as the "maximum" issue.In order to make the shape of Sigmoid function as close as possible to the original increasing linear function, in particular, to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multiobjective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree  is the minimum value in the membership function of all objective functions, that is, According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by equations ( 1)-( 20) and (37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about EH is optimized including single objective optimization and multi-objective optimization, where the various detailed data are obtained from the lower level optimization results.Secondly, facing different engineering application modes, the lower level dispatching center uploads the data of the lower level optimization results to the upper level dispatching center, which prepare for upper level optimization.Then, the upper level optimization also includes single objective optimization and The Sigmoid function curve has the same semi-open shape as increasing linear function curve, so it is suitable for the membership function described as the "maximum" issue.In order to make the shape of Sigmoid function as close as possible to the original increasing linear function, in particular, to keep the coincidence of the points (waist values) of the two curves at µ(F n ) = 0.5, the points at µ(F n ) = 0.5 and 0.9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.
The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multi-objective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree The Sigmoid function curve has the same semi-open shape as increasing linear function curve, so it is suitable for the membership function described as the "maximum" issue.In order to make the shape of Sigmoid function as close as possible to the original increasing linear function, in particular, to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multiobjective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree  is the minimum value in the membership function of all objective functions, that is, According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by equations ( 1)-( 20) and (37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about EH is optimized including single objective optimization and multi-objective optimization, where the various detailed data are obtained from the lower level optimization results.Secondly, facing different engineering application modes, the lower level dispatching center uploads the data of the lower level optimization results to the upper level dispatching center, which prepare for upper level optimization.Then, the upper level optimization also includes single objective optimization and The Sigmoid function curve has the same semi-open shape as increasing linear function curve, so it is suitable for the membership function described as the "maximum" issue.In order to make the shape of Sigmoid function as close as possible to the original increasing linear function, in particular, to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multiobjective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree  is the minimum value in the membership function of all objective functions, that is, According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by equations ( 1)-( 20) and (37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about EH is optimized including single objective optimization and multi-objective optimization, where the various detailed data are obtained from the lower level optimization results.Secondly, facing different engineering application modes, the lower level dispatching center uploads the data of the lower level optimization results to the upper level dispatching center, which prepare for upper level = min {µ(F n )}.According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of The optimization problem with fuzzy objective is called fuzzy optimization problem.When there is objective fuzzy optimization problem.According t method, the satisfaction degree  is the minimum v functions, that is, According to the theory, the multi-objective optimization problem maximization with meeting all constraints, that is, th Based on fuzzy theory, a multi-objective and IES is developed, the continuous differentiable membership function.Therefore, the nonlinear pro suitable to solve above multi-objective optimization functions.

Optimization Framework
The framework of hierarchical multi-objectiv provided in Figure 5. From Figure 5, the high-orde performance and optimization models are presented by equations ( 1)-( 20) and (37)-(39).And subsyste models are proposed in upper level for IES, which c on the relationship between master and slave dispatc EH is optimized including single objective optimizat various detailed data are obtained from the lowe different engineering application modes, the lower maximization with meeting all constraints, that is, the max  The Sigmoid function curve has the same semi-open shape as increasing linear functio so it is suitable for the membership function described as the "maximum" issue.In order to m shape of Sigmoid function as close as possible to the original increasing linear function, in p to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one is called fuzzy optimization problem.When there is more than one objective, it is then calle objective fuzzy optimization problem.According to the maximum minimum satisfactio method, the satisfaction degree  is the minimum value in the membership function of all functions, that is, According to the maximum and minimum principle of theory, the multi-objective optimization problem can be transformed into the problem maximization with meeting all constraints, that is, the max  problem, which is described Based on fuzzy theory, a multi-objective fuzzy collaborative optimization mode and IES is developed, the continuous differentiable Sigmoid function is introduced as an im membership function.Therefore, the nonlinear programming methods are better reason suitable to solve above multi-objective optimization model for applying the improved mem functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization m provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under o performance and optimization models are presented in lower level for EH, which can be rep by equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and opti models are proposed in upper level for IES, which can be represented by equations ( 21)-(3 on the relationship between master and slave dispatching center, Firstly, the lower level mod EH is optimized including single objective optimization and multi-objective optimization, w various detailed data are obtained from the lower level optimization results.Secondl problem, which is described by: max The Sigmoid function curve has the same semi-open shape as increasing linear function curve, so it is suitable for the membership function described as the "maximum" issue.In order to make the shape of Sigmoid function as close as possible to the original increasing linear function, in particular, to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multiobjective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree  is the minimum value in the membership function of all objective functions, that is, According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about The optimization problem with fuzzy objective function and constraint condition or one of is called fuzzy optimization problem.When there is more than one objective, it is then called m objective fuzzy optimization problem.According to the maximum minimum satisfaction d method, the satisfaction degree  is the minimum value in the membership function of all obj functions, that is, According to the maximum and minimum principle of fuz theory, the multi-objective optimization problem can be transformed into the problem maximization with meeting all constraints, that is, the max  problem, which is described by Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model fo and IES is developed, the continuous differentiable Sigmoid function is introduced as an imp membership function.Therefore, the nonlinear programming methods are better reasonabl suitable to solve above multi-objective optimization model for applying the improved membe functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization meth provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-d performance and optimization models are presented in lower level for EH, which can be repres by equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and optimiz models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).The Sigmoid function curve has the same semi-open shape as increasing linear function curve so it is suitable for the membership function described as the "maximum" issue.In order to make th shape of Sigmoid function as close as possible to the original increasing linear function, in particula to keep the coincidence of the points (waist values) of the two curves at ( )=0.The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called mult objective fuzzy optimization problem.According to the maximum minimum satisfaction degre method, the satisfaction degree  is the minimum value in the membership function of all objectiv functions, that is, According to the maximum and minimum principle of fuzzy se theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improve membership function.Therefore, the nonlinear programming methods are better reasonable an suitable to solve above multi-objective optimization model for applying the improved membershi functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method i provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-desig performance and optimization models are presented in lower level for EH, which can be represente by equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and optimizatio Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by Equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by Equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about EH is optimized including single objective optimization and multi-objective optimization, where the various detailed data are obtained from the lower level optimization results.Secondly, facing different engineering application modes, the lower level dispatching center uploads the data of the lower level optimization results to the upper level dispatching center, which prepare for upper level optimization.Then, the upper level optimization also includes single objective optimization and multi-objective optimization, after optimizing the upper level model about IES, and the upper level dispatching center feeds back the updated scheduling plan to the lower level dispatching center.This "hierarchical collaborative optimization" strategy can well adapt to and reflect real world energy system where engineering application modes are diversified.
Energies 2019, 12, x FOR PEER REVIEW 18 of 28 multi-objective optimization, after optimizing the upper level model about IES, and the upper level dispatching center feeds back the updated scheduling plan to the lower level dispatching center.This "hierarchical collaborative optimization" strategy can well adapt to and reflect real world energy system where engineering application modes are diversified.In view of the hierarchical multi-objective fuzzy collaborative optimization model developed in this paper, where it's essentially MINLP, general commercial optimization software, such as LINGO, GAMS, CPLEX and so on, can be used to solve it.Programs are written based on the LINGO18.0 optimization software platform (LINDO & LINGO (TIANJIN) TECH.CO., LTD, TIANJIN, CHINA) and the global solver is utilized to solve the above problems.

System Description
The typical days in winter are taken as the research object, which is divided into 24 periods (that is, an optimal dispatch cycle).The specific structure of the model is shown in Figure 6.
IES consists mainly of 13 power system nodes, eight natural gas system nodes, 12 heat system nodes and three EHs (recorded as EH1, EH2 and EH3, respectively).The three systems are directly coupled together through CHP unit, in which the fuel required by CHP is mainly provided by biomass and garbage, and the natural gas imported to CHP unit plays a role of afterburning.Natural gas system and heat system are directly coupled through GB.To take full advantage of the geographical resource endowment, wind farm (WT3) provides heat energy to heat system through EB during the day.In order to ensure safe and reliable power supply of power system, the power system interconnects with the large power grid through node 13.Natural gas system mainly relies on node 1 gas source and node 2 gas source to obtain natural gas.Three systems are indirectly coupled with EH1, EH2 and EH3, and the electricity, heat and gas input into EH work together.Three EH in lower layer and three systems in the upper layer all work together to maximize the comprehensive benefit.In view of the hierarchical multi-objective fuzzy collaborative optimization model developed in this paper, where it's essentially MINLP, general commercial optimization software, such as LINGO, GAMS, CPLEX and so on, can be used to solve it.Programs are written based on the LINGO18.0 optimization software platform (LINDO & LINGO (TIANJIN) TECH.CO., LTD, TIANJIN, CHINA) and the global solver is utilized to solve the above problems.

System Description
The typical days in winter are taken as the research object, which is divided into 24 periods (that is, an optimal dispatch cycle).The specific structure of the model is shown in Figure 6.
IES consists mainly of 13 power system nodes, eight natural gas system nodes, 12 heat system nodes and three EHs (recorded as EH1, EH2 and EH3, respectively).The three systems are directly coupled together through CHP unit, in which the fuel required by CHP is mainly provided by biomass and garbage, and the natural gas imported to CHP unit plays a role of afterburning.Natural gas system and heat system are directly coupled through GB.To take full advantage of the geographical resource endowment, wind farm (WT3) provides heat energy to heat system through EB during the day.In order to ensure safe and reliable power supply of power system, the power system interconnects with the large power grid through node 13.Natural gas system mainly relies on node 1 gas source and node 2 gas source to obtain natural gas.Three systems are indirectly coupled with EH1, EH2 and EH3, and the electricity, heat and gas input into EH work together.Three EH in lower layer and three systems in the upper layer all work together to maximize the comprehensive benefit.
The maximum startup times of FC, GT, GB, GSHP, MT, and so on in a scheduling cycle are all 4. The parameter settings of each device under off-design performance are all from the corresponding references.All the factors of equations ( 1)-( 7) are reported in Table A1 of Appendix A, and IES subsystem topology information are showed in Table A2 of Appendix B. GSHP is allowed to work during periods from 7 to 24.TRAN's efficiency and TOU references [10,12,23].The units of power in this paper are MW unless otherwise specified.The objectives selected for EH1 are F 1 , F 2 and F 3 .For EH2, F 1 and F 4 are selected as its objectives.For EH3, F 1 and F 4 are selected for its objectives.f 1 and f 2 are selected as IES's objectives.The maximum startup times of FC, GT, GB, GSHP, MT, and so on in a scheduling cycle are all 4. The parameter settings of each device under off-design performance are all from the corresponding references.All the factors of equations ( 1)-( 7) are reported in Table A1 of Appendix A, and IES subsystem topology information are showed in Table B1 of Appendix B. GSHP is allowed to work during periods from 7 to 24.TRAN's efficiency and TOU references [10,12,23].The units of power in this paper are MW unless otherwise specified.The objectives selected for EH1 are 1 F , 2 F and 3 F .
For EH2, 1 F and 4 F are selected as its objectives.For EH3, 1 F and 4 F are selected for its objectives. 1 f and 2 f are selected as IES's objectives.

Results and Analyses
The results of single objective and multi-objective fuzzy optimization are shown in Table 1.

Results and Analyses
The results of single objective and multi-objective fuzzy optimization are shown in Table 1., the points at ( )=0.5 n F μ and 0.9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.
The optimization problem with fuzzy objective function and constraint condition or one of them is called fuzzy optimization problem.When there is more than one objective, it is then called multiobjective fuzzy optimization problem.According to the maximum minimum satisfaction degree method, the satisfaction degree  is the minimum value in the membership function of all objective functions, that is, According to the maximum and minimum principle of fuzzy set theory, the multi-objective optimization problem can be transformed into the problem of  maximization with meeting all constraints, that is, the max  problem, which is described by: ( . .0 1 Remark 4. Based on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH and IES is developed, the continuous differentiable Sigmoid function is introduced as an improved membership function.Therefore, the nonlinear programming methods are better reasonable and suitable to solve above multi-objective optimization model for applying the improved membership functions.

Optimization Framework
The framework of hierarchical multi-objective fuzzy collaborative optimization method is provided in Figure 5. From Figure 5, the high-order nonlinear models of devices under off-design performance and optimization models are presented in lower level for EH, which can be represented by equations ( 1)-( 20) and ( 37)-(39).And subsystem "energy system" models and optimization models are proposed in upper level for IES, which can be represented by equations ( 21)-(39).Based on the relationship between master and slave dispatching center, Firstly, the lower level model about EH is optimized including single objective optimization and multi-objective optimization, where the various detailed data are obtained from the lower level optimization results.Secondly, facing different engineering application modes, the lower level dispatching center uploads the data of the lower level optimization results to the upper level dispatching center, which prepare for upper level optimization.Then, the upper level optimization also includes single objective optimization and of multi-objective fuzzy optimization is 0.705874.At this time, operating fuel cost of multi-objective fuzzy optimization increases by 3.84% compared with F 1 , pollutant gas emission increases by 4.85% compared with F 2 , and primary fossil energy consumed increases by 1.6% compared with F 3 .Multi-objective fuzzy optimization balances the contradictions and conflicts among each single objective, and takes into account the comprehensive requirements of each objective, and obtains an optimal scheduling scheme which is satisfactory to all parties.Multi-objective fuzzy optimization can more accurately reflect the actual operation of EHs and IES, and has obvious advantages over single objective optimal scheduling model.
In this paper, only FC is taken as an example to analyze the variation of device energy efficiency under off-design performance and the difference of device operation mode under different optimization objectives.As shown in Figure 7, when the generation power of FC is high, its power generation efficiency will be reduced accordingly.During the time period from 1 to 10, the power generation efficiency of FC under different optimization objectives is relatively stable, but there is a significant difference between power generation efficiency under F 1 and F 2 , and the cost of operating fuel under F 1 is the least and the fuel can be saved under the running mode of high generation efficiency.However, the minimum emission is sought under F 2 , and the CO 2 emission factor of FC is obviously smaller than that of GT.It is obvious that the full power operation of FC under F 2 is much helpful to reduce emissions.the contradictions and conflicts among each single objective, and takes into account the comprehensive requirements of each objective, and obtains an optimal scheduling scheme which is satisfactory to all parties.Multi-objective fuzzy optimization can more accurately reflect the actual operation of EHs and IES, and has obvious advantages over single objective optimal scheduling model.
In this paper, only FC is taken as an example to analyze the variation of device energy efficiency under off-design performance and the difference of device operation mode under different optimization objectives.As shown in Figure 7, when the generation power of FC is high, its power generation efficiency will be reduced accordingly.During the time period from 1 to 10, the power generation efficiency of FC under different optimization objectives is relatively stable, but there is a significant difference between power generation efficiency under 1 F and 2 F , and the cost of operating fuel under 1 F is the least and the fuel can be saved under the running mode of high generation efficiency.However, the minimum emission is sought under 2 F , and the 2 CO emission factor of FC is obviously smaller than that of GT.It is obvious that the full power operation of FC under 2 F is much helpful to reduce emissions.Under the multi-objective fuzzy collaborative optimization mode, change of energy level including electrical and heat energy storage equipment and renewable energy accommodation for EH1 during each period are shown in Figure 8. Results from Figure 8 show that renewable energy is completely absorbed, although multi-objective optimization of EH1 does not take into account the optimal problem about renewable energy accommodation.Increasing the penetration of renewable energy by EH1 can reduce the cost of purchasing electricity and natural gas while EH1 runs under the actual optimized scheduling and reduce fossil fuel combustion emissions, at the same time, increasing the given renewable energy can be equivalent to reduce the consumption of fossil fuels.As shown in Figure 8, the time periods from 1 to 9 are the periods of valley-flat electricity price, and ES is always in the charging state.The periods from 10 to16 are the periods of flat-peak electricity price, ES is always in the state of discharge.PV power can be fully absorbed during time periods from 17 to 19.ES can effectively reduce operating costs, pollution emissions and fossil energy consumption in optimal scheduling cycle.HS is mainly in the state of heat storage from time period 10 to time period 17, at the same time, and this period is the low period of heat load.It is in the exothermic state Under the multi-objective fuzzy collaborative optimization mode, change of energy level including electrical and heat energy storage equipment and renewable energy accommodation for EH1 during each period are shown in Figure 8. Results from Figure 8 show that renewable energy is completely absorbed, although multi-objective optimization of EH1 does not take into account the optimal problem about renewable energy accommodation.Increasing the penetration of renewable energy by EH1 can reduce the cost of purchasing electricity and natural gas while EH1 runs under the actual optimized scheduling and reduce fossil fuel combustion emissions, at the same time, increasing the given renewable energy can be equivalent to reduce the consumption of fossil fuels.As shown in Figure 8, the time periods from 1 to 9 are the periods of valley-flat electricity price, and ES is always in the charging state.The periods from 10 to16 are the periods of flat-peak electricity price, ES is always in the state of discharge.PV power can be fully absorbed during time periods from 17 to 19.ES can effectively reduce operating costs, pollution emissions and fossil energy consumption in optimal scheduling cycle.HS is mainly in the state of heat storage from time period 10 to time period 17, at the same time, and this period is the low period of heat load.It is in the exothermic state during the time periods from 21 to 6, where this period is the peak and valley period of heat load.It can be seen that HS can effectively deflate and bridge the peak and valley difference of heat load.Energy flows mode of EH are optimized by the coordination of various energies and equipment, which highlights the benefits of coordinated optimal operation.In EH2, GB is the direct coupling device of natural gas and heat energy, while MT is the direct key coupling device of electricity, heat and natural gas.GB and MT work together, in which the amount of natural gas entering MT from natural gas node 8 can be reflected by the allocation proportion factor ( t MT v ). Figure 9 shows change of power generation efficiency and allocation proportion factor of MT, and power of SC for EH2 under different optimization objectives in different time periods.As can be seen from Figure 9, as long as MT works, no matter in which optimization mode, t MT v 's value is greater than 0.8, which indicates that, compared with GB, MT's operation plays a major role in consuming natural gas, while GB actually plays a role in afterburning.In addition, t MT v of MT is relatively stable under 1 F , that is, natural gas node 8 supply of natural gas will be relatively stable.In general, power generation efficiency of MT is relatively stable under different optimization scenarios, all of which are around 0.38.The stable power generation efficiency of MT is also good for ABS stabilization work, which is beneficial to the coordination of MT and ABS. Figure 9 shows that under different optimization modes, the actual output of SC is significantly different.There are no energy storages equipment in EH2, compared with EH1, and it is clear that there is no free and flexible absorption of solar energy.But under 4 F , the actual power accepted output of SC has reached 78.58%, which is relatively acceptable in actual operation, and multi-objective fuzzy optimization has achieved the desired results.In EH2, GB is the direct coupling device of natural gas and heat energy, while MT is the direct key coupling device of electricity, heat and natural gas.GB and MT work together, in which the amount of natural gas entering MT from natural gas node 8 can be reflected by the allocation proportion factor (v t MT ). Figure 9 shows change of power generation efficiency and allocation proportion factor of MT, and power of SC for EH2 under different optimization objectives in different time periods.As can be seen from Figure 9, as long as MT works, no matter in which optimization mode, v t MT 's value is greater than 0.8, which indicates that, compared with GB, MT's operation plays a major role in consuming natural gas, while GB actually plays a role in afterburning.In addition, v t MT of MT is relatively stable under F 1 , that is, natural gas node 8 supply of natural gas will be relatively stable.In general, power generation efficiency of MT is relatively stable under different optimization scenarios, all of which are around 0.38.The stable power generation efficiency of MT is also good for ABS stabilization work, which is beneficial to the coordination of MT and ABS. Figure 9 shows that under different optimization modes, the actual output of SC is significantly different.There are no energy storages equipment in EH2, compared with EH1, and it is clear that there is no free and flexible absorption of solar energy.But under F 4 , the actual power accepted output of SC has reached 78.58%, which is relatively acceptable in actual operation, and multi-objective fuzzy optimization has achieved the desired results.
Node 1 in the natural gas system is the pressure reference point, which is also the gas source point.The results of Figure 10 show that the pressure values of each node are (0.8p.u., 1.2p.u.) in the whole IES optimal scheduling cycle as opposed to the pressure of the reference point, which fully meet the pressure setting requirements.The nodal pressure of the compressor downstream is mainly adjusted by the compressor to maintain its safety and reasonable range.The nodal pressure of compressor upstream is mainly regulated by natural gas source 2, that is, node 2. From Figure 10, it can be found that the pressure fluctuation of node 2 is relatively large, and the natural gas source 2 plays a role in balancing the pressure stability of the whole natural gas system's node to a certain extent.Node 1 in the natural gas system is the pressure reference point, which is also the gas source point.The results of Figure 10 show that the pressure values of each node are (0.8p.u., 1.2p.u.) in the whole IES optimal scheduling cycle as opposed to the pressure of the reference point, which fully meet the pressure setting requirements.The nodal pressure of the compressor downstream is mainly adjusted by the compressor to maintain its safety and reasonable range.The nodal pressure of compressor upstream is mainly regulated by natural gas source 2, that is, node 2. From Figure 10, it can be found that the pressure fluctuation of node 2 is relatively large, and the natural gas source 2 plays a role in balancing the pressure stability of the whole natural gas system's node to a certain extent.Node 1 in the natural gas system is the pressure reference point, which is also the gas source point.The results of Figure 10 show that the pressure values of each node are (0.8p.u., 1.2p.u.) in the whole IES optimal scheduling cycle as opposed to the pressure of the reference point, which fully meet the pressure setting requirements.The nodal pressure of the compressor downstream is mainly adjusted by the compressor to maintain its safety and reasonable range.The nodal pressure of compressor upstream is mainly regulated by natural gas source 2, that is, node 2. From Figure 10, it can be found that the pressure fluctuation of node 2 is relatively large, and the natural gas source 2 plays a role in balancing the pressure stability of the whole natural gas system's node to a certain extent.Figure 11 shows that the actual temperature in the room well tracks the expected temperature well during IES optimization scheduling cycle, which indicates that the thermal parameter equivalent model developed for buildings is effective.The change of the actual temperature in the room not only meets the requirements of thermal users, but also increases the flexibility of IES scheduling and endows heat load with certain flexible characteristics.
In Figure 6, the heat system is a central heating system with double heat sources, multiple heat load types and loops.The loop can guarantee the reliability of heating, and some pipelines on the loop are often in a standby state in practical operation.From Figure 12, it can be found that heat source #1 mainly provides heat energy to nearby heat load #3 and heat load #4, while heat source #2 provides heat energy to nearby heat load #5, heat load #6 and heat load #7.This kind of nearby heating characteristic of heat source is very common under the practical engineering application, which gives priority to nearby heating, avoids the long distance transmission heat transfer medium flow, and reduces the heat energy loss.In Figure 6, the heat system is a central heating system with double heat sources, multiple heat load types and loops.The loop can guarantee the reliability of heating, and some pipelines on the loop are often in a standby state in practical operation.From Figure 12, it can be found that heat source #1 mainly provides heat energy to nearby heat load #3 and heat load #4, while heat source #2 provides heat energy to nearby heat load #5, heat load #6 and heat load #7.This kind of nearby heating characteristic of heat source is very common under the practical engineering application, which gives priority to nearby heating, avoids the long distance transmission heat transfer medium flow, and reduces the heat energy loss.

Conclusions
According to the requirements of practical engineering applications in real world situation, this paper proposed a hierarchical multi-objective fuzzy collaborative optimization model for IES under off-design performance.Lower level multi-objective fuzzy collaborative optimization model mainly centered on EH, where models of various devices under off-design performance were established in detail.Lower level model is mainly oriented to regional or industrial park IES, which can adapt to the optimization and regulation analysis of IES in the real world.Upper level multi-objective fuzzy collaborative optimization model of IES including power system, natural gas system, and heat system was proposed based on the relationship between master and slave dispatch centers.Upper level model has large potential and wide applications in analyzing optimal operation of regional multienergy coupling systems in the real world.
The models were applied and verified by a chosen case study of IES, where the case and optimization model in this paper can well adapt to and reflect real world energy system.The results  In Figure 6, the heat system is a central heating system with double heat sources, multiple heat load types and loops.The loop can guarantee the reliability of heating, and some pipelines on the loop are often in a standby state in practical operation.From Figure 12, it can be found that heat source #1 mainly provides heat energy to nearby heat load #3 and heat load #4, while heat source #2 provides heat energy to nearby heat load #5, heat load #6 and heat load #7.This kind of nearby heating characteristic of heat source is very common under the practical engineering application, which gives priority to nearby heating, avoids the long distance transmission heat transfer medium flow, and reduces the heat energy loss.

Conclusions
According to the requirements of practical engineering applications in real world situation, this paper proposed a hierarchical multi-objective fuzzy collaborative optimization model for IES under off-design performance.Lower level multi-objective fuzzy collaborative optimization model mainly centered on EH, where models of various devices under off-design performance were established in detail.Lower level model is mainly oriented to regional or industrial park IES, which can adapt to the optimization and regulation analysis of IES in the real world.Upper level multi-objective fuzzy collaborative optimization model of IES including power system, natural gas system, and heat system was proposed based on the relationship between master and slave dispatch centers.Upper level model has large potential and wide applications in analyzing optimal operation of regional multienergy coupling systems in the real world.
The models were applied and verified by a chosen case study of IES, where the case and optimization model in this paper can well adapt to and reflect real world energy system.The results

Conclusions
According to the requirements of practical engineering applications in real world situation, this paper proposed a hierarchical multi-objective fuzzy collaborative optimization model for IES under off-design performance.Lower level multi-objective fuzzy collaborative optimization model mainly centered on EH, where models of various devices under off-design performance were established in detail.Lower level model is mainly oriented to regional or industrial park IES, which can adapt to the optimization and regulation analysis of IES in the real world.Upper level multi-objective fuzzy collaborative optimization model of IES including power system, natural gas system, and heat system was proposed based on the relationship between master and slave dispatch centers.Upper level model has large potential and wide applications in analyzing optimal operation of regional multi-energy coupling systems in the real world.
The models were applied and verified by a chosen case study of IES, where the case and optimization model in this paper can well adapt to and reflect real world energy system.The results demonstrated the availability of regarding the proposed multi-objective fuzzy collaborative optimization under off-design performance conditions as the basis of performance analysis of devices, calculation and operation analysis of EH, optimal scheduling, collaborative operation of hybrid multi-energy system, the planning of different operation modes, and engineering calculation analysis of IES.The i-th optimization objective function The unit power price of energy i (CNY/MW) κ i j Equivalent pollutant gas emission factor i that consumes unit energy j (t/MW) The maximum charging power factor for device i (%) The maximum discharging power factor for device i (%) p i The pressure of node i in the natural gas system (bar) f The natural gas flow consumed by the compressor in the natural gas system (SCF/h) The variable reflecting the direction of natural gas flow in the pipeline q Steady state pipeline flow in the natural gas system (SCF/h) k com The compressor characteristic constant on fuzzy theory, a multi-objective fuzzy collaborative optimization model for EH loped, the continuous differentiable Sigmoid function is introduced as an improved nction.Therefore, the nonlinear programming methods are better reasonable and e above multi-objective optimization model for applying the improved membership n Framework work of hierarchical multi-objective fuzzy collaborative optimization method is ure 5. From Figure 5, the high-order nonlinear models of devices under off-design d optimization models are presented in lower level for EH, which can be represented 1)-( 20) and (37)-(39).And subsystem "energy system" models and optimization posed in upper level for IES, which can be represented by equations ( 21)-(39).Based hip between master and slave dispatching center, Firstly, the lower level model about including single objective optimization and multi-objective optimization, where the d data are obtained from the lower level optimization results.Secondly, facing eering application modes, the lower level dispatching center uploads the data of the imization results to the upper level dispatching center, which prepare for upper level hen, the upper level optimization also includes single objective optimization and The satisfaction degree Appendix A

Figure 1 .
Figure 1.A typical FC case under off-design performance: (a) Schematic diagram of FC working principle; (b) Off-design performance of FC.

Figure 1 .
Figure 1.A typical FC case under off-design performance: (a) Schematic diagram of FC working principle; (b) Off-design performance of FC.

Figure 2 .
Figure 2. Pipeline with compressor driven by gas turbine.

Figure 2 .
Figure 2. Pipeline with compressor driven by gas turbine.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.
9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.
9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.

5
9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.

Energies 2019 ,Figure 4 . 9 (
Figure 4. Fuzzy membership function: (a) Increas Decreasing linear function and reverse Sigmoid func The Sigmoid function curve has the same semi so it is suitable for the membership function describe shape of Sigmoid function as close as possible to the to keep the coincidence of the points (waist values ( )=0.5 n F μ

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function Decreasing linear function and reverse Sigmoid function.

5
9 (shoulder value) in this paper are taken as coincident points, which can to obtain the Sigmoid characteristic parameters, and vice versa.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.
9 (shoulder value) in this paper are taken as coincident points, which can be used to obtain the Sigmoid characteristic parameters, and vice versa.

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b Decreasing linear function and reverse Sigmoid function.The Sigmoid function curve has the same semi-open shape as increasing linear function c so it is suitable for the membership function described as the "maximum" issue.In order to mak shape of Sigmoid function as close as possible to the original increasing linear function, in parti to keep the coincidence of the points (waist values) of the two curves at ( )=0.5 n F μ , the poi ( )=0.5 n F μ

Figure 4 .
Figure 4. Fuzzy membership function: (a) Increasing linear function and Sigmoid function; (b) Decreasing linear function and reverse Sigmoid function.
9 (shoulder value) in this paper are taken as coincident points, which can be use to obtain the Sigmoid characteristic parameters, and vice versa.

Figure 5 .
Figure 5. Framework of proposed fuzzy collaborative optimization method for IES.

Figure 5 .
Figure 5. Framework of proposed fuzzy collaborative optimization method for IES.

Figure 6 .
Figure 6.Schematic diagram of the integrated energy system.

Figure 6 .
Figure 6.Schematic diagram of the integrated energy system.

Figure 7 .
Figure 7. Change of power generation efficiency and electrical output for FC under different optimization objectives in different time periods.

Figure 7 .
Figure 7. Change of power generation efficiency and electrical output for FC under different optimization objectives in different time periods.

Energies 2019 ,
12, x FOR PEER REVIEW 21 of 28during the time periods from 21 to 6, where this period is the peak and valley period of heat load.It can be seen that HS can effectively deflate and bridge the peak and valley difference of heat load.Energy flows mode of EH are optimized by the coordination of various energies and equipment, which highlights the benefits of coordinated optimal operation.

Figure 8 .
Figure 8. Change of energy storages and renewable energy accommodation for EH1 under multiobjective fuzzy optimization mode in different time periods.

Figure 8 .
Figure 8. Change of energy storages and renewable energy accommodation for EH1 under multi-objective fuzzy optimization mode in different time periods.

Figure 9 .
Figure 9. Change of power generation efficiency and allocation proportion factor of MT, and power of SC for EH2 under different optimization objectives in different time periods.

Figure 10 .
Figure 10.Change of nodal pressure in natural gas system under 12 f .

Figure 11
Figure11shows that the actual temperature in the room well tracks the expected temperature well during IES optimization scheduling cycle, which indicates that the thermal parameter equivalent model developed for buildings is effective.The change of the actual temperature in the room not only meets the requirements of thermal users, but also increases the flexibility of IES scheduling and endows heat load with certain flexible characteristics.

Figure 9 .
Figure 9. Change of power generation efficiency and allocation proportion factor of MT, and power of SC for EH2 under different optimization objectives in different time periods.

Figure 9 .
Figure 9. Change of power generation efficiency and allocation proportion factor of MT, and power of SC for EH2 under different optimization objectives in different time periods.

Figure 10 .
Figure 10.Change of nodal pressure in natural gas system under 12 f .

Figure 11
Figure11shows that the actual temperature in the room well tracks the expected temperature well during IES optimization scheduling cycle, which indicates that the thermal parameter equivalent model developed for buildings is effective.The change of the actual temperature in the room not only meets the requirements of thermal users, but also increases the flexibility of IES scheduling and endows heat load with certain flexible characteristics.

Figure 10 .
Figure 10.Change of nodal pressure in natural gas system under f 12 .

Figure 12 .
Figure 12.Changes in the mass flow rate of the pipelines connecting sources and loads under 12f .

Figure 11 .
Figure 11.Indoor temperature changes of buildings under f 12 .

Figure 11 .
Figure 11.Indoor temperature changes of buildings under 12 f .

Figure 12 .
Figure 12.Changes in the mass flow rate of the pipelines connecting sources and loads under 12f .

Figure 12 .
Figure 12.Changes in the mass flow rate of the pipelines connecting sources and loads under f 12 .

.fUmb
The flow vector consumed by compressor in natural gas pipeline YThe node admittance matrix.The vector of nodes' voltage.The pipeline or branch flow vector in heat system m The flow rate in heat system (kg/s)λ pThe heat transfer coefficient in heat system (J/(kg• • C) Φ l ,Φ s Heat output of the heat source and heat demand of the load (MW) α k ,β k ,γ k The kth inequality constraint coefficients of CHP operating region F t g The amount of natural gas (MW) P t e The amount of electric energy (MW) H t h The amount of heat energy (MW) V t node Node voltage at time t (p.u.) S t b Branch capacity at time t (MW•h) r t Compressor compression ratio at time t κ n Characteristic parameter of increasing linear function ∆ n Characteristic parameter of decreasing linear function function: (a) Increasing linear function and Sigmoid function; (b) linear function and reverse Sigmoid function.idfunction curve has the same semi-open shape as increasing linear function curve, for the membership function described as the "maximum" issue.In order to make the id function as close as possible to the original increasing linear function, in particular, ncidence of the points (waist values) of the two curves at ( )=0.5 n F μ , the points at 0.9 (shoulder value) in this paper are taken as coincident points, which can be used gmoid characteristic parameters, and vice versa.ization problem with fuzzy objective function and constraint condition or one of them optimization problem.When there is more than one objective, it is then called multioptimization problem.According to the maximum minimum satisfaction degree isfaction degree  is the minimum value in the membership function of all objective is, the maximum and minimum principle of fuzzy set lti-objective optimization problem can be transformed into the problem of  ith meeting all constraints, that is, the max  problem, which is described by: Difference of indoor room actual temperature T a,t in and the expected temperature T e,t in is in a certain range [−T min in ,T max in ].(2)The total actual temperature of indoor rooms ∑ T a,t in is equal to the total expected temperature of indoor rooms ∑ T e,t in within T period.It is given as:

Single Objective Optimal Results of Single Objective Optimal Results of Multi-Objective Fuzzy Optimization Result of 1
. Single-objective and multi-objective fuzzy collaborative optimization results.EH1

Single Objective Optimal Results of Single Objective Optimal Results of Multi-Objective Fuzzy Optimization Result of F 1 Result of F 2 Result of F 3
. Single-objective and multi-objective fuzzy collaborative optimization results.EH1

Single Objective Optimal Results of Single Objective Optimal Results of Multi-Objective Fuzzy Optimization Result of f 1 Result of f 2
n F μ

i
Electric part load rate of device i at time t (MW/MW) X t Binary-variable of device i operation state at time t Constant coefficient of the n-th power item for device I model H t

i
Output thermal power of device i at time t (MW) H t

i
Thermal part load rate of device i at time t (MW/MW) W t type Energy storage capacity of device type at time t (MW•h) P cha,t type Charging power values of device type at time t (MW) Discharging power values of device type at time t (MW) Certain characteristic temperature of device or system i ( • C) F i