Next Article in Journal
R&D Intensity and Its Curvilinear Relationship with Firm Profitability: Perspective from the Alternative Energy Sector
Next Article in Special Issue
Energy Transition at Home: A Survey on the Data and Practices That Lead to a Change in Household Energy Behavior
Previous Article in Journal
Growth of Firms in a Fragmented Cultural Industry: Italian Commercial Art Galleries’ Competitive Strategies
Previous Article in Special Issue
Event Matching Classification Method for Non-Intrusive Load Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flexible Continuous-Time Modeling for Multi-Objective Day-Ahead Scheduling of CHP Units

Faculty of Electrical and Computer Engineering, University of Tabriz, Tabriz 5166616471, Iran
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(9), 5058; https://doi.org/10.3390/su13095058
Submission received: 17 March 2021 / Revised: 25 April 2021 / Accepted: 26 April 2021 / Published: 30 April 2021
(This article belongs to the Special Issue Energy Systems Integration: From Policy-Makers to Consumers)

Abstract

:
Increasing applications of CHP units have turned the problem of finding the best optimization model into a significant subject for scholars. In this respect, this paper is aimed at driving a novel formulation to the multi-objective day-ahead scheduling of CHP units using Bernstein polynomials, which more optimally schedules power and heat generations as well as ramping trajectories. This procedure includes yielding an affine function that closely approximates real-time net-load and generation trajectories, which is demonstrated to have a superior performance to the conventional hourly day-ahead scheduling of CHP units based on discrete-time approximation. The problem of how to handle various objective functions by function space method is also addressed. The simulations conducted on the sample test systems, which consist of CHP systems, thermal and heat-only units, as well as thermal and electrical loads, show that the suggested multi-objective model can perfectly cover the total heat and electrical loads in terms of economic and environmental criteria. More importantly, the results indicate that the accuracy of the proposed approach renders cost saving of 1.67% and emission saving of 1.46% in comparison with the conventional hourly-based model, apart from leading to fewer ramping scarcities in real-time operations.

1. Introduction

Nowadays, the ever-increasing need for electric energy power, especially in developing countries, has led to an increase in the consumption of fossil fuels. The growing tendency for fossil fuel together with the environmental concerns they induce by air pollution in power generation have compelled operators to use distributed energy resources [1,2], such as wind power, solar energy, etc. Among all these resources, thermal power plants have been widely regarded. Flexibility and reliable performance of combined heat and power (CHP) [3] units with great efficiency and high security are among the main reasons these units have been holding so much appeal among all available thermal power units [4,5]. The goal of the CHP system is to simultaneously provide electrical and thermal energies for loads. Hence, while thermal units operate, the CHP system captures squandered heat and utilizes this energy to satisfy thermal loads [6]. Furthermore, the extra heat from an electric power plant can be utilized for industrial aims or other goals, which boosts the efficiency of the system [7]. For these reasons and more, CHP units are widely used in the planning and operation of power systems. Given that biomass is being promoted to substitute coal for CHP plants, the importance of large-scale CHP plants in a 100% renewable energy context will be fully recognized in the near future [8].
Accordingly, in order to gain an acceptable power system operation and planning, disparate features such as minimizing fuel costs, losses, and environmental pollution, must be considered, to name but a few. In this respect, some studies have been conducted to co-optimize the CHP units. For instance, a multi-objective stochastic model is proposed for scheduling CHP units based on microgrid (MG) considering uncertainties in load and wind power in [5]. Another multi-objective electric model is developed in [9] based on the integration of thermal units, heat-only, and CHP units, where the objective functions include minimization of total cost and maximization of heat generation, and the normal boundary intersection method is applied to solve the model. Moreover, a novel optimal operation of a CHP-based MG is presented in [10] regarding energy storage system and demand response in the form of a multi-objective self-scheduling optimization problem. Both convex and nonconvex models, together with the multi-objective issues, need suitable approaches to cope with the multi-objective issue and both rely on finding a well-organized set of Pareto optimal solutions [11]. Therefore, various methods are expounded to solve multi-objective optimization problems (MOOPs), among all of which, the ε-constrained method is highly practical in handling multi-objective problems [12].
The ε-constrained method is an extensive problem-solving technique to cope with multi-objective issues in which the Pareto front can be generated exactly by changing the ε-vector, and all but one of the objective functions are transformed into constraints. By way of example, the ε-constraint method is utilized to solve the problem of self-scheduling home energy management systems, which is defined as a multi-objective optimization problem in [13].
On the other side, some studies have been conducted on the scheduling of CHP units. For example, the problem of optimal day-ahead scheduling of CHP units with electric and thermal energy storages regarding security constraints is discussed in [14], which includes the deterministic optimal day-ahead unit commitment of generation resources relying on system loads. In another study, a whale optimization algorithm is presented as a new method for solving CHP economic dispatch issues [15]. Moreover, a new perspective is proposed based on microgrid in [16] to resolve the economic dispatch of CHP units, thermal units, heat-only units, wind turbines, photovoltaic systems, and battery storage, where AC optimal power flow is utilized. In addition, the information gap theory is applied to handle the uncertainties of generating units, market price and load. In [17], optimal operation of CHP units in a microgrid is discussed based on minimizing the total costs and contenting heat demand. Likewise, it is indicated that taking the feasible operation region (FOR) of CHPs into consideration, heat and power outputs are inseparable and rest on each other in a manner that one output effects the other one. The stochastic fractal search, which is a novel optimization solver, is utilized to solve CHP economic dispatch in [18]. On the other hand, some studies have been carried out on the linearized formulation of thermal units [19] and CHP units either in the economic dispatch issue or unit commitment [20,21]. A number of state-of-the-art optimal dispatch reviews of CHP in integrated energy systems and CHP operation can be found in [22,23,24].
Continuous-time modeling, rather than discrete-time modeling, is a new subfield of power systems modeling in which the problem is solved in a continuous-time framework with linear objective functions and linear constraints. A trajectory of the system’s state over time is smoothly given everywhere in the phase space by integrating a continuous-time model over time, leading to the fact that instantaneous abrupt changes, which are inevitable in discrete-time models, would not be seen in continuous-time models.
This technique was firstly presented in [25,26]. In this way, to solve the continuous-time infinite-dimensional problem, developing a numerically stable way based on decreasing the dimensionality of the continuous-time problem is necessary. The suggested method in these references is a function space-based method with Bernstein polynomials which has been suggested for modeling the other upcoming problems. For example, an optimization scheduling of thermal generating units integrated with energy storage in a continuous-time framework is discussed in [27], where a continuous-time Gaussian technique is applied to cope with load uncertainty. The outcomes show the proficiency of their presented technique, which reduces total cost. Additionally, a novel formulation is presented in [28] for the locational marginal price of electricity, which is the precisely combined fluctuation of load demand and operation constraints of power network in evaluating market price. Their suggested method converts the problem into a continuous-time optimal control problem, where the numerical results depict the flexibility and accuracy of the proposed method. Furthermore, a new optimization model is represented in [29] for balancing energy and flexible ramp products.
Despite the differences in formulation and solution techniques, current CHP scheduling models have deployed the discrete-time load and generation trajectories which lead to an obvious contrast between real-time operation and day-ahead operation. Therefore, in order to attain any degree of flexibility or economical operation of CHP units, accurate day-ahead scheduling modeling is demanded.
To bridge the above research gaps, a novel and flexible scheduling model is presented in this study for CHP units, which contains a linear formulation of thermal heat-only units and CHP units, while a continuous-time optimization model is applied to schedule CHP units for the first time in literature. To be more specific, a novel technique has been applied as a continuous-time model for the CHP systems in day-ahead power systems operation. As a matter of fact, the multi-objective scheduling of CHP integrated with thermal and heat-only units is solved by means of the epsilon constraint method in a continuous-time framework, where the objective functions comprise cost and emission. In fact, a flexible optimal model is found to schedule the generating units in comparison with day-ahead scheduling. It is worth noting that although solution time of the represented model is a bit more compared to the hourly model, it engenders a higher accuracy, not to mention that the total cost is considerably lower compared to the hourly model. The presented co-objective optimization problem is resolved by employing the ε-constraint method, which is a powerful approach capable of solving issues with non-convex Pareto front. The presented technique produces a single-objective function by transforming all but one of the objectives into constraints. A function space method is then suggested to handle the problem.
The rest of the paper is organized as follows: The general framework of the discrete-time and continuous-time multi-objective CHP planning are presented in Section 2, while Section 3 presents the use of Bernstein polynomials to model the problem in detail. The numerical results of applying the suggested model on a test system are presented in Section 4, followed by summarizing the conclusions in Section 5.

2. Problem Modeling

2.1. Discrete-Time Multi-Objective CHP Scheduling

The goal here is to optimize the multi-objective operation problem of a CHP-based system to obtain an optimal resource scheduling of the system over the day-ahead scheduling horizon 𝓣 = [ 0 , T ] . Discrete-time scheduling subdivides the day-ahead time horizon into M  intervals, 𝓣 = m 0 M 1 𝓣 m ,   𝓣 m = [ t m , t m + 1 ) with the same length T m = t m + 1 t m where t 0 = 0 and t M = T . The objective of this model is to minimize the overall generation cost for energy as well as the atmospheric pollutants, supplying both thermal and electrical load demands with respect to several constraints. In this study, it is assumed that spare heat is wasted if the produced heat is more than the heat demand, and it is provided with auxiliary boilers if it is less than the heat demand. Accordingly, three types of units are considered: power-only units, CHP-based units, and heat-only units.

2.1.1. Multi-Objective Programming Based on ε-Constraint Method

Decision-making problems mostly rely on multiple criteria and can infrequently be defined based on a well-defined single objective. In fact, although, multiple facets of a decision process can be put into a single objective, this simplification hardly captures the complexity of real-world decision applications. Given that solving multi-criteria decision issues is not as simple as single-criterion issues, finding a single solution that satisfies all constraints and optimizes different objective functions is not simply possible. As a result, providing a set of non-dominated solutions and finding the “most preferred” alternative among the obtained solutions can be the best bet. As a matter of fact, the ε-constraint method is one of the best intuitive approaches to produce the Pareto optimal solutions, which takes into account the user’s preferences. This method keeps one of the objectives only (the most important one) and restricts the rest of the objectives within user-specific values ( f i ( x , u ) ε i ). A generalized hourly trigeneration planning framework to arrange the optimal operating strategy is formulated as follows.

2.1.2. Objective Functions

1. Minimizing the total cost of operation ( F Cos t Total ) : The primary objective function aims at minimizing the operational cost of the system including the production cost of thermal units, CHP units, and heat-only units for a definite time horizon. Such an objective function can be described as:
min t m 𝓣 [ F Cos t Total ( P ( t m ) , H ( t m ) , c , x ( t m ) ) , F Emission Total ( P ( t m ) , H ( t m ) ) ]
where
F Cos t Total ( P ( t m ) , H ( t m ) , c , x ( t m ) ) = f Cos t TH ( P TH ( t m ) ) + f Cos t CHP ( c , x ( t m ) ) + f Cos t H ( H H ( t m ) ) , t m 𝓣
and
F Emission Total ( P ( t m ) , H ( t m ) ) = f Emission TH ( P TH ( t m ) ) + f Emission CHP ( P CHP ( t m ) ) + f Emission H ( H H ( t m ) ) , t m 𝓣
The total non-linear cost of traditional thermal power units can be approximated either as a quadratic convex or as piecewise-linear functions. The quadratic form can be formulated as below [30]:
f Cos t TH ( P TH ( t m ) ) = i 𝓘 c 2 i , TH P i , TH 2 ( t m ) + c 1 i , TH P i , TH ( t m ) + c 0 i , TH , t m 𝓣
The heat and power produced by CHP plants are highly coupled and are not generated at similar efficiencies and prices. The maximum efficiency of a cogeneration unit can be reached at the full extent of power and proportional heat. The so-called characteristic points of the feasible area comprise these kinds of operating points. A representative operating area is illustrated in Figure 1 with 6 characteristic points in terms of cost ( c ) , power ( p ) , and heat ( h ) . An operating region is convex if, given any two points, it contains the entire line segment that joins them. Accordingly, the feasible area is supposed to be convex and the cost can be also considered a convex function of the generated power and heat [31]. Based on this assumption, the total cost of a CHP unit can be defined as convex combinations of the corner points as below [32].
f Cos t CHP ( c , x ( t m ) ) = j 𝓙 l 𝓛 c j l x j l ( t m ) ,    t m 𝓣
Here, the production cost of the CHP f Cos t CHP can be defined as a combination of fuel costs, service costs, taxes, etc. This single plant model can be used as a basis for both the hourly multiple plant CHP planning model and the multi-period planning model.
Generally, in order to supply high heat load demands and increase the flexibility of CHP units, heat-only units are added to the system. The quadratic convex cost function of these units can be defined as follows [9]:
f Cos t H ( H H ( t m ) ) = k 𝓚 c 2 k , H H k , H 2 ( t m ) + c 1 k , TH H k , H ( t m ) + c 0 k , H , t m 𝓣
It should be noticed that the quadratic cost functions of thermal (2) and heat-only units (4) can be approximated by a piecewise linear cost function to keep the linearity of the formulation [25].
2. Minimizing the total emission ( F Emission Total ) : The subsidiary objective function which is taken into consideration is the total pollutant emission. In general, this function is described as the second-order polynomial or linear forms. The first part of the second objective function is the minimization of the total emission of traditional thermal units described as the second-order polynomial form [33]:
f Emission TH ( P TH ( t m ) ) = i 𝓘 e 2 i , TH P i , TH 2 ( t m ) + e 1 i , T H P i , TH ( t m ) + e 0 i , T H , t m 𝓣
Minimization of the amount of pollutants released into the air by CHP plants, defined as the second part of the additional objective function, can be formulated as below [9]:
f Emission CHP ( P CHP ( t m ) ) = j 𝓙 e 1 j , CH P P j , CHP ( t m ) + e 0 j , CHP , t m 𝓣
Finally, the last part of the secondary objective function includes the mitigation of pollutants emitted by heat-only units, which can be described in a linear form as follows [9]:
f Emission H ( H H ( t m ) ) = k 𝓚 e 1 k , H H k , H ( t m ) + e 0 k , H ,    t m 𝓣
Considering the ε-constraint, Equation (1) can be reformulated as below:
min t m 𝓣 F Cos t Total ( P ( t m ) , H ( t m ) , c , x ( t m ) )
subject to
t m 𝓣 F Emission Total ( P ( t m ) , H ( t m ) ) ε 2 , n 2 ,    n 2 = 0 , 1 , , p 2
Here, p2 denotes the density of Pareto solutions associated with the minimization of pollutant emission.

2.2. Continuous-Time Multi-Objective CHP Scheduling Model

As mentioned before, the traditional CHP planning problem includes discrete-time generation schedules as a subdivision of 𝓣 to the same M intervals 𝓣 m = [ t m , t m + 1 ) ,   𝓣 = m 0 M 1 𝓣 m with length T m = t m + 1 t m , e.g., hourly. In this way, the upshot power and heat generation schedule to meet discrete-time load samples at minimum cost and emission and sequentially to solve the traditional M single-period economic dispatch problems [34]. The use of a discrete-time CHP planning model implies that units must track piecewise constant generation trajectories from one schedule to the next, where the units’ ramping are obtained from the finite difference between two successive time steps. Obviously, not only can continuous-time scheduling appropriately model faster variations of the net-load to avoid ramping scarcity events, but it can also allow a greater flexibility for generating units.
Defining the continuous-time generation trajectories P ( t ) = [ P TH ( t ) ,   P CHP ( t ) ] and H ( t ) = [ H CHP ( t ) ,   H H ( t ) ] , the continuous-time multi-objective CHP scheduling problem over the day-ahead scheduling horizon 𝓣 subject to the operating constraints can be formulated as follows:
min 𝓣 ( F Cos t Total ( P ( t ) , H ( t ) , c , x ( t ) ) ) d t
subject to
𝓣 ( F Emission Total ( P ( t ) , H ( t ) ) ) d t ε 2 , n 2 , t 𝓣 , n 2 = 0 , 1 , , p 2
d d t P TH ( t ) = P TH ( t ) , t 𝓣
d d t ( times ( p CHP , x CHP ( t ) ) ) = P CHP ( t ) ,    t 𝓣
P load ( t ) = 1 N TH T P TH ( t ) + 1 N CHP T P CHP ( t ) , t 𝓣
H heat ( t ) = 1 N H T H H ( t ) + 1 N CHP T H CHP ( t ) , t 𝓣
P _ TH P TH ( t ) P ¯ TH , t 𝓣
P _ CHP ( H CHP ( t ) ) P CHP ( t ) P ¯ j , CHP ( H CHP ( t ) ) , t 𝓣
H _ CHP ( P CHP ( t ) ) H CHP ( t ) H ¯ CHP ( P CHP ( t ) ) , t 𝓣
H _ H H H ( t ) H ¯ H , t 𝓣
P CHP ( t ) = diag ( p CHP x CHP T ( t ) ) ,    t 𝓣
H CHP ( t ) = diag ( h CHP x CHP T ( t ) ) ,    t 𝓣
x CHP ( t ) 1 L = 1 ,    t 𝓣
0 x CHP ( t ) 1 , t 𝓣
P ˙ _ TH P ˙ TH ( t ) P ˙ ¯ TH , t 𝓣
P ˙ _ CHP P ˙ CHP ( t ) P ˙ ¯ CHP , t 𝓣
Here, the ramping trajectories of thermal and CHP units are described in (14) and (15), wherein times ( p CHP , x CHP ( t ) ) is a matrix obtained by multiplying the p CHP and x CHP ( t ) arrays element by element. The total operation cost (12) is equal to (2), representing the real power generation cost of thermal units (4) plus the total cost of the CHP units (5) and heat-only units (6). The ε-constraint is applied to the secondary objective function by (13). (15) and (16) represent the continuous-time real power and heat balance constraints, respectively, wherein 1 N TH , 1 N CHP , and 1 N H signify N TH , N CHP , and N H -dimensional vectors of ones. (18)–(21) are the box constraints of the power and heat trajectories.
Moreover, (22) and (23) embody the relationship between power production and characteristic points as well as heat production and characteristic points of the CHP units, wherein diag() yields a column vector of the main diagonal elements of multiplying h CHP and x CHP T ( t ) . Convexity and non-negativity constraints for CHPs are also enforced in (24) and (25), wherein 1 L represents the L-dimensional vector of ones. (26) and (27) signify the ramping trajectories of thermal and CHP units over 𝓣 , where P ˙ TH ( t ) and P ˙ CHP ( t ) are obtained as follows,
d d t P TH ( t ) = lim T m 0 P TH ( t m + 1 ) P TH ( t m ) T m , t 𝓣
d d t P CHP ( t ) = lim T m 0 ( times ( p CHP , x CHP ( t m + 1 ) ) ) ( times ( p CHP , x CHP ( t m ) ) ) T m ,   t 𝓣

3. Cubic Spline Model for Multi-Objective Scheduling of CHPs

Although modeling in infinite-dimensional spaces is more precise, it is, in general, more computationally intractable. Under suitable assumptions, let us develop an approximation bridge from the infinite-dimensional linear programming to tractable finite convex models in which the performance of the approximation is computed clearly. In this paper, modeling is proposed by means of the cubic spline function space spanned by Bernstein polynomials, which enables making up a function and refining the approximation to any desired precision while also providing additional flexibility to fit the continuous-time load variations [25]. Bernstein polynomials are first defined to provide a foundation for the rest of the paper. The Bernstein-basis polynomials of degree Q include Q + 1 polynomials can be described as:
b q ( Q ) ( t ) = ( q Q ) t q ( 1 t ) ( Q q ) , q = 0 , , Q ,    t [ 0 , 1 ]
where ( q Q ) is the binomial coefficient.
A convex combination of Bernstein-basis polynomials is called a Bernstein polynomial. The Bernstein polynomial (degree Q) of a continuous-time function on the interval [0, 1] is described as:
B Q ( f ) ( t ) = q = 0 Q f ( q Q ) b q ( Q ) ( t )
Figure 2a simply displays Bernstein-basis polynomials of degree 3 on the unit square and Figure 2b displays an approximation of a continuous function as a linear combination of Bernstein basis polynomials of degree 3.
Notation 1. The smallest convex set l (shaded in Figure 2b) that contains the control points (marked by squares) is the so-called “convex hull”, inside which the whole curve is confined.
Notation 2. Considering M subdivisions of the scheduling horizon, a vector of Bernstein-basis functions over the entire 𝓣 can be obtained as:
v m ( Q + 1 ) + q ( Q ) ( t ) = b q ( Q ) ( t t m T m ) ,    t [ t m , t m + 1 ] ,
Notation 3. It is well worth considering that the derivatives of the Qth degree Bernstein polynomials are also polynomials of degree (Q − 1) and can be expressed as linear combination of Bernstein polynomials.
d d t b q ( Q ) ( t ) = Q ( b q 1 ( Q 1 ) ( t ) b q ( Q 1 ) ( t ) )

3.1. Modeling the Generation/Demand Trajectories and Operating Constraints

3.1.1. Continuous-Time Modeling of the Generation/Demand Trajectories

The projection of generation/demand trajectory on the Bernstein function space is defined in (34)–(36):
P TH ( t ) = P TH v ( Q ) ( t ) ,     P CHP = P CHP v ( Q ) ( t )
H H ( t ) = H H v ( Q ) ( t ) ,     H CHP = H CHP v ( Q ) ( t )
P load ( t ) = P load v ( Q ) ( t ) ,    H heat ( t ) = H heat v ( Q ) ( t )
where P TH , P CHP , and H H are the N TH × ( M ( Q + 1 ) ) , N CHP × ( M ( Q + 1 ) ) , and N H × ( M ( Q + 1 ) ) mapping matrices associated with the thermal units’ power, CHP units’ power and heat, and heat-only units’ heat generation, respectively. Bernstein coefficient row vectors of the electric and heat demand are described by the ( M ( Q + 1 ) ) -dimensional vectors, respectively.

3.1.2. Continuous-Time Modeling of the Power Balance Equations

Considering these coefficients, power balance equations can be reformulated based on Bernstein polynomials as (37) and (38) and then (39) and (40):
P load v ( Q ) ( t ) = 1 N TH T P TH v ( Q ) ( t ) + 1 N CHP T P CHP v ( Q ) ( t ) ,     t 𝓣
H heat v ( Q ) ( t ) = 1 N H T H H v ( Q ) ( t ) + 1 N CHP T H CHP v ( Q ) ( t ) ,     t 𝓣
Bernstein-basis functions can be removed from both sides of (33) and (34). Thus, the following equations are obtained:
P load = 1 N TH T P TH + 1 N CHP T P CHP ,     t 𝓣
H heat = 1 N H T H H + 1 N CHP T H CHP ,     t 𝓣

3.1.3. Continuous-Time Modeling of the Box Constraints

The inequality constraints in the Bernstein function space can be reformulated in the following. Satisfying the convex hull property is one of the strengths of Bernstein polynomials, which keeps the continuous-time trajectories inside the convex hull [35].
P _ TH P TH v ( Q ) ( t ) P ¯ TH ,    t 𝓣
P _ CHP ( H CHP v ( Q ) ( t ) ) P CHP v ( Q ) ( t ) P ¯ j , CHP ( H CHP v ( Q ) ( t ) ) , t 𝓣
H _ H H H v ( Q ) ( t ) H ¯ H ,    t 𝓣

3.1.4. Continuous-Time Modeling of the CHP Characteristic Points

The relationship between power, heat, and the coordinates of the corner points and other constraints of CHP units in the Bernstein function space can be defined as follows:
P CHP v ( Q ) ( t ) = diag ( p CHP ( x CHP v ( Q ) ( t ) ) T ) ,    t 𝓣
H CHP v ( Q ) ( t ) = diag ( h CHP ( x CHP v ( Q ) ( t ) ) T ) ,    t 𝓣
( x CHP v ( Q ) ( t ) ) 1 L = 1 ,    t 𝓣
0 x CHP v ( Q ) ( t ) 1 ,    t 𝓣

3.1.5. Modeling Ramping Constraints

The ramping power constraints of the units, which restrict the value power production of the thermal or CHP units, in terms of the Bernstein representation over 𝓣 , are equal to:
P ˙ _ TH P ˙ TH v ( Q ) ( t ) P ˙ ¯ TH ,    t 𝓣
P ˙ _ CHP P ˙ CHP v ( Q ) ( t ) P ˙ ¯ CHP ,    t 𝓣
According to Notation 2, the derivative of power trajectories of the Qth degree Bernstein polynomials is polynomials of degree (Q − 1). This significant property helps to define the continuous-time ramping trajectory of the power and generating units in the Bernstein function space as:
P ˙ TH ( t ) = P TH v ˙ ( Q ) ( t ) = P TH P TH v ( Q 1 ) ( t ) = P ˙ TH v ( Q 1 ) ( t )
P ˙ CHP ( t ) = P CHP v ˙ ( Q ) ( t ) = P CHP P CHP v ( Q 1 ) ( t ) = P ˙ CHP v ( Q 1 ) ( t )
where P TH and P CHP are the M ( Q + 1 ) × M Q mapping matrices, and P ˙ TH and P ˙ CHP are matrices of order N TH × M Q and N CHP × M Q , respectively.

3.2. Modeling the Objective Functions

3.2.1. Continuous-time modeling of the cost function as the main objective

In order to benefit from the Bernstein polynomials, linear approximation of the non-linear functions is essential. To accomplish of which, the piecewise linear version of (4) and (6) are first described as (51) and (52):
f Cos t TH ( P TH ( t ) ) = c 2 i , TH P _ i , TH 2 + c 1 i , TH P _ i , TH + c 0 i , TH + s = 0 N S 1 s i , TH s p i , TH s ( t ) ,    i 𝓘 , t 𝓣
where
P i , TH ( t ) = P _ i , TH + s = 0 N S 1 p i , TH s ( t ) ,   i 𝓘 , t 𝓣
f Cos t H ( H H ( t ) ) = c 2 k , H H _ k , H 2 + c 1 k , H H _ k , H + c 0 k , H + s = 0 N S 1 s k , H s h k , H s ( t ) ,    k 𝓚 , t 𝓣
where H k , H ( t ) = H _ k , H + s = 0 N S 1 h k , H s ( t ) ,    k 𝓚 , t 𝓣 and h k , H s ( t ) in the Bernstein function-space can be expressed as:
p i , TH s ( t ) = p i , TH s T v ( Q ) ( t )
h k , H s ( t ) = h k , H s T v ( Q ) ( t )
where p i , TH s and h k , H s are the M ( Q + 1 ) -dimensional vectors of Bernstein coefficients of the thermal and CHP units’ cost function coefficients, respectively. Accordingly, the cost functions obtained from the proposed function space-based method can be reformulated as (57), where 1 M ( Q + 1 ) is the M ( Q + 1 ) -dimensional all-ones vector.
𝓣 ( F Cos t Total ( P ( t ) , H ( t ) , c , x ( t ) ) ) d t = m = 0 M 1 [ i = 1 N TH T m ( c 2 i , TH P _ i , TH 2 + c 1 i , TH P _ i , TH + c 0 i , TH ) + k = 1 N H T m ( c 2 k , H H _ k , H 2 + c 1 k , H H _ k , H + c 0 k , H ) + T m s = 0 N S 1 [ i = 1 N TH s i , TH s ( p i , TH s T 1 M ( Q + 1 ) Q + 1 ) + k = 1 N H s k , H s ( h k , H s T 1 M ( Q + 1 ) Q + 1 ) ] + T m ( diag ( p CHP x CHP T ) Q + 1 ) ]

3.2.2. Continuous-Time Modeling of the Emission Function Based on ε-Constraint

Let us project the total pollutant emission function converted to an inequality constraint by virtue of the ε-constraint in the space spanned by v ( Q ) ( t ) as (58),
𝓣 ( F Emission Total ( P ( t ) , H ( t ) ) ) d t = [ m = 0 M 1 [ T m i = 1 N TH ( e 2 i , TH P _ i , TH 2 + e 1 i , TH P _ i , TH + e 0 i , TH ) + T m j = 1 N CHP c 0 j , CHP + T m k = 1 N H ( c 1 k , H H _ k , H + c 0 k , H )      + T m s = 0 N S 1 [ i = 1 N TH s i , TH s ( p i , TH s T 1 M ( Q + 1 ) Q + 1 ) + k = 1 N H s k , H s ( h k , H s T 1 M ( Q + 1 ) Q + 1 ) ] + T m c 1 j , CHP ( diag ( p CHP x CHP T ) Q + 1 ) ] ] ε 2 , n 2 , n 2 = 0 , 1 , , p 2

3.2.3. Enforcing the Continuity Constraints

Let us confirm that the generating and ramping trajectories projected on the function space are continuous. Hence, C0 continuity is required to be imposed to check the continuity of the generation trajectories at the edge points of the intervals, while C1 continuity is crucial to be enforced to guarantee at interval connection points.
P i , TH , m ( Q + 1 ) + Q = P i , TH , ( m + 1 ) ( Q + 1 ) ,    i 𝓘 , m = 0 , , M 1 ,
P j , CHP , m ( Q + 1 ) + Q = P j , CHP , ( m + 1 ) ( Q + 1 ) ,    j 𝓙 , m = 0 , , M 1 ,
1 T m ( P i , TH , m ( Q + 1 ) + Q P i , TH , m ( Q + 1 ) + Q 1 ) = 1 T m + 1 ( P i , TH , ( m + 1 ) ( Q + 1 ) + 1 P i , TH , ( m + 1 ) ( Q + 1 ) ) ,    i 𝓘 , m = 0 , , M 1 ,
1 T m ( P j , CHP , m ( Q + 1 ) + Q P j , CHP , m ( Q + 1 ) + Q 1 ) = 1 T m + 1 ( P j , CHP , ( m + 1 ) ( Q + 1 ) + 1 P j , CHP , ( m + 1 ) ( Q + 1 ) ) ,     j 𝓙 , m = 0 , , M 1 ,

4. Numerical Results and Discussion

To assess the proposed multi-objective model, a system [36] including four thermal units, two CHP units, and two heat-only units was considered over a 24 h time interval. The limits of thermal units and heat-only units are given in the Appendix A, Table A1 and Table A2, respectively. Further, the characteristic points of the CHP units depending on power and heat generation are described in Table A3 in the Appendix A. Netload data of the California ISO (CAISO) for 22 June 2018 [37] were used for continuous-time power and heat demand. The corresponding data are scaled back to the peak loads of 2048.84 MW and 1024.42 MW for the electric and heat demands, respectively. Bernstein polynomials were used to construct the proposed model, and all simulations were performed in a General Algebraic Modeling System (GAMS) [38] environment with a CPLEX solver on a system with a Core i7-7500U processor, 64 × 2.90-GHz, 12 GB of RAM.
The resultant electrical load profile over the Bernstein function space and hourly-based model are captured in Figure 3. The resultant heat load profile is also similar to the electrical load profile with different peak loads. To compare the results more clearly, simulations were conducted in two different cases, based on the load profiles of Figure 3, categorized as studies I and II:
Study I: In this case, the problem is defined as the minimization of total cost and pollutant emission as two competing objectives. Discrete-time modeling is the basis of this study, which aims to evaluate the leading properties of continuous-time modeling. Generally, to solve MO optimization problems, the optimization process is required to [39]:
  • Keep the solution points, the decision space, and the non-dominated solutions in the objective space;
  • Preserve the algorithmic progress toward the Pareto front;
  • Keep the diversity of Pareto front solutions;
  • Provide a large enough but limited number of solutions for the Pareto front;
  • Search the best compromise solutions (BCSs) among the solution points of the Pareto front.
In this work, the ε-constraint method is suggested for providing the Pareto front, and Fuzzy multi-criteria decision-makers are utilized for selecting the optimal compromise solutions of the two case studies in the presence of two criteria. Detailed information regarding this procedure can be found in [40]. Using a fuzzy decision maker, hourly scheduling obtained the best trade-off between objectives at a point with 4,048,725 (USD) and 385,565 (kg). The generation trajectories of the best compromise solution associated with the thermal, CHP, and heat-only units based on the traditional hourly model are displayed in Figure 4 and Figure 5, and the obtained results are listed in Table 1.
Study II: The next set of numerical analyses includes applying the bi-objective Bernstein-polynomial-based modeling. More detailed results are given in Table 1, in which the difference between the real-time operations is clearly seen. A demonstrative set of Pareto fronts for the total day-ahead emission versus the total day-ahead cost provided by hourly-based and continuous-time scheduling is demonstrated in Figure 6.
The horizontal axis represents the cost minimization as the main objective, while the minimization of emission pollutants has been shown on the vertical axis. To better clarify, the location of the best compromises of both models is shown clearly and redrawn inside the figure. Further, according to Figure 7, the presented CHP and thermal units model furnishes a continuous-time schedule for generating units in the electrical part that ably utilizes their capacity to chase the continuous-time net-load, whenever energy loss is reduced. In Figure 8, the continuous-time schedule from CHP and heat-only units in the heat part is captured. Evaluating the above results and the interpretation of related Pareto front approximations, the following significant points are revealed:
(a)
Even though the solution points in both Pareto curves in Figure 6 are well distributed, it is clear that the solutions of the continuous-one outperform the non-dominated solutions of the hourly model. However, according to the results of Table 1, the BCS of continuous-time modeling is obtained with a trade-off between 3269 (USD) increase in total day-ahead cost and 434 (kg) increase in total emission with regards to the hourly scheduling. Meanwhile, the real-time operation cost and emission are 113,180 (USD) and 10,655 (kg) for Case I as well as 40,227 (USD) and 4443 (kg) for Case II, respectively. This result highlights cost-saving of 69,684 (USD) (1.67%) and emission saving of 57,778 (kg) (1.46%) in day-ahead scheduling, which is substantial.
(b)
Ramping scarcities can increase overall system costs because of the penalty prices due to the ramping capability shortage. Prevention of more ramp scarcities is another fundamental characteristic of the Bernstein expansion, which delivers accuracy and tighter function. For evaluating this attribute, the ramping values of the thermal units are determined by scaling down the maximum value of the units by a ratio of fifteen. In this case, the results represent 6 ramping scarcity events for the hourly based model in comparison with no scarcity events for the proposed model.
(c)
The last row in Table 1 shows the computing time of the case studies in a minute, while the upper bound on the duality gap is set to zero. The execution time of the hourly day-ahead operation is 0.3648 min, which is lower than that of the proposed function space-based model, 1.9433 min. It is clear that the proposed method is approximately slower by 5 due to having more variables and constraints in the Bernstein polynomial-based modeling. This shows the main drawback of this expansion, especially when the accuracy of polynomials is increased. However, reducing the accuracy of the polynomials can bring about a significant reduction in the computational burden. Hence, relying on the problem, a tradeoff between accuracy and time should be determined.

5. Conclusions

This paper represents a different way of optimization for continuous-time scheduling of CHP regarding thermal and heat-only units in day-ahead power systems. A flexible optimal model is proposed to optimize cost and emission, simultaneously, in the function space considering the ramping requirements of thermal and CHP units. The paper starts with reformulating the problem as a linear scheduling programming, followed by subsequently evolving a polynomial-based method to solve the continuous-time problem, which includes modeling the parameters in a finite order function space formed by Bernstein polynomials. Numerical results using real load data obtained from the ISO California show that the presented model decreases the total operation cost and emission as well as the number of ramping scarcities in the real-time operations compared with the hourly model. Extending the proposed model to one that includes uncertain parameters is the scope of our future study.

Author Contributions

Investigation, M.A.; Methodology, S.B.-S. and B.M.-I.; Project administration, B.M.-I. and M.A.; Resources, M.A.; Software, E.D.; Writing—original draft, E.D. and S.B.-S.; Writing—review & editing, B.M.-I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Research Affairs Office of the University of Tabriz, Tabriz, Iran, under grant number S-27-8 as part of a postdoctoral research project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors appreciate the support of Research Affairs Office of the University of Tabriz, Tabriz, Iran (grant number S-27-8).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A. 
Objective Functions
F Cos t Total Total generation cost of units (USD)
F Emission Total Total emission of units (kg)
f Cos t TH , f Cos t CHP and f Cos t H Generation cost of thermal, CHP and heat-only units
f Emission TH , f Emission CHP and f Emission H CO2 emission of thermal, CHP and heat-only units
B. 
Vectors/Matrices
P ( t ) Power vector including power generation of units P ( t ) = [ P TH ( t ) , P CHP ( t ) ] (MW)
P load ( t ) / P load ( t n ) Vector of electric load demand (MW)
P TH ( t ) / P TH ( t n ) Vector of thermal units’ power (MW)
P CHP ( t ) Vector of CHP units’ power (MW)
H ( t ) / H ( t n ) Heat vector including heat generation of units H ( t ) = [ H CHP ( t ) , H H ( t ) ] (MW)
H heat ( t ) / H heat ( t n ) Vector of heat demand (MW)
H CHP ( t ) / H CHP ( t n ) Vector of CHP units’ heat (MW)
H H ( t ) / H H ( t n ) Vector of heat-only units’ heat (MW)
P ˙ TH ( t ) , P ˙ CHP ( t ) Time derivatives of the thermal and CHP units’ power generation
P ˙ ¯ TH , P ˙ _ TH Vector of thermal units’ ramp up and ramp down, respectively (MW/h)
P ˙ ¯ CHP , P ˙ _ CHP Vector of CHP units’ ramp up and ramp down, respectively (MW/h)
c CHP , p CHP , h CHP , x CHP ( t / t m ) Cost, power, heat, and coefficient matrices ( N CHP -by- L ) of characteristic points for CHP units, respectively
1 R R-dimensional vector of ones
I R Identity matrix of order R
e R Standard vector (with a 1 in the Rth coordinate and zeros elsewhere)
C. 
Parameters
N TH The number of thermal units
N CHP The number of CHP units
N H The number of heat-only units
N S The number of piecewise segments
L The number of characteristic points
Q Degree of Bernstein polynomials
e m , q Q ( t ) Bernstein-basis functions of degree Q
D. 
Variables
P i , TH ( t ) Power generation output of thermal unit i
P j , CHP ( t ) Power generation output of CHP unit j
H j , CHP ( t ) Heat production of CHP unit j
H k , H ( t ) Heat production of heat-only unit k
P ˙ ¯ j , CHP , P ˙ _ j , CHP CHP units’ ramp up and ramp down, respectively (MW/h)
x j , CHP l Decision variable encoding the convex combination of the operating region of CHP unit j
c j , CHP l , p j , CHP l , h j , CHP l Characteristic points of CHP unit j
p i , TH ( t ) , h k , H ( t ) Auxiliary power and heat generation variables of thermal and heat-only units, respectively
E. 
Symbol, Sets and Indices
i Index of thermal units
j Index of CHP units
l Index of characteristic points of CHP units
k Index of heat-only units
f Index of linearization segments
m Index of the hourly time interval
𝓘 Index set of thermal units
𝓙 Index set of CHP units
𝓛 Index set of characteristic points of CHP units
𝓚 Index set of heat-only units
𝓕 Index set of linearization segments
𝓜 Index set of hourly time interval
𝓣 Day-ahead scheduling horizon
F. 
Superscripts and subscripts
TH or TH Thermal
CHP or CHP Combined heat and power
H or H Heat-only
_ / ¯ Minimum/maximum magnitude operator
G. 
Constants
c 2 i , TH , c 1 i , TH and c 0 i , TH Cost coefficients of thermal units i
c 2 k , H , c 1 k , H and c 0 k , H Cost coefficients of heat-only units k
e 2 i , TH , e 1 i , TH and e 0 i , TH Emission coefficients of thermal units i
e 1 k , H and e 0 k , H Emission coefficients of heat-only units k
e 1 j , CHP and e 0 j , CHP Emission coefficients of CHP units j
H. 
Operators
T Transpose operator
diag ( ) diag’s operator, which returns a column vector of the main diagonal elements of
times ( , ) Element-wise multiplication operator which multiplies arrays and , element by element

Appendix A

Table A1. Required data of thermal units.
Table A1. Required data of thermal units.
Thermal Unit c 2 i , TH c 1 i , TH c 0 i , TH e 2 i , TH e 1 i , TH e 0 i , TH P ¯ TH P _ TH
TH 10.0206914.8357.110.00683−0.5455140.26698000
TH 20.0323218.5457.110.00683−0.5455140.26695000
TH 30.0106560.26126.80.00461−0.5111642.89555000
TH 40.0422221.1957.110.00461−0.5111642.89554500
Table A2. Required data of heat-only units.
Table A2. Required data of heat-only units.
Heat-Only Unit c 2 i , H c 1 i , H c 0 i , H e 1 i , H e 0 i , H P ¯ H P _ H
H 10.010510.5523.4260.3276713.859323000
H 20.02999.2110.7210.3276713.859322700
Table A3. Required data of CHP units.
Table A3. Required data of CHP units.
CHP Unit p CHP A , B , C , D p CHP A , B , C , D c CHP A , B , C , D e 1 i , CHP e 0 i , CHP
CHP 1247, 215, 118, 890, 240, 104.8, 01306.85, 1954.2646, 705.6, 982−0.5734311.57
CHP 2125, 110, 114, 750, 235, 35, 01030, 2123, 1120, 1419−1.7669821.65

References

  1. Zappa, W.; Junginger, M.; van den Broek, M. Is a 100% renewable European power system feasible by 2050? Appl. Energy 2019, 233, 1027–1050. [Google Scholar] [CrossRef]
  2. Renau, J.; García, V.; Domenech, L.; Verdejo, P.; Real, A.; Giménez, A.; Sánchez, F.; Lozano, A.; Barreras, F. Novel Use of Green Hydrogen Fuel Cell-Based Combined Heat and Power Systems to Reduce Primary Energy Intake and Greenhouse Emissions in the Building Sector. Sustainability 2021, 13, 1776. [Google Scholar] [CrossRef]
  3. Kim, D.; Kim, K.-T.; Park, Y.-K. A Comparative Study on the Reduction Effect in Greenhouse Gas Emissions between the Combined Heat and Power Plant and Boiler. Sustainability 2020, 12, 5144. [Google Scholar] [CrossRef]
  4. Mohammadi, A.; Venkatesh, G.; Sandberg, M.; Eskandari, S.; Joseph, S.; Granström, K. A Comprehensive Environmental Life Cycle Assessment of the Use of Hydrochar Pellets in Combined Heat and Power Plants. Sustainability 2020, 12, 9026. [Google Scholar] [CrossRef]
  5. Pourghasem, P.; Sohrabi, F.; Abapour, M.; Mohammadi-Ivatloo, B. Stochastic multi-objective dynamic dispatch of renewable and CHP-based islanded microgrids. Electr. Power Syst. Res. 2019, 173, 193–201. [Google Scholar] [CrossRef]
  6. Gbadamosi, S.L.; Nwulu, N.I. Optimal power dispatch and reliability analysis of hybrid CHP-PV-wind systems in farming applications. Sustainability 2020, 12, 8199. [Google Scholar] [CrossRef]
  7. Mohammadi-Ivatloo, B.; Moradi-Dalvand, M.; Rabiee, A. Combined heat and power economic dispatch problem solution using particle swarm optimization with time varying acceleration coefficients. Electr. Power Syst. Res. 2013, 95, 9–18. [Google Scholar] [CrossRef]
  8. Lund, R.; Mathiesen, B.V. Large combined heat and power plants in sustainable energy systems. Appl. Energy 2015, 142, 389–395. [Google Scholar] [CrossRef]
  9. Ahmadi, A.; Moghimi, H.; Nezhad, A.E.; Agelidis, V.G.; Sharaf, A.M. Multi-objective economic emission dispatch considering combined heat and power by normal boundary intersection method. Electr. Power Systems Res. 2015, 129, 32–43. [Google Scholar] [CrossRef]
  10. Aghaei, J.; Alizadeh, M.-I. Multi-objective self-scheduling of CHP (combined heat and power)-based microgrids considering demand response programs and ESSs (energy storage systems). Energy 2013, 55, 1044–1054. [Google Scholar] [CrossRef]
  11. Soroudi, A.; Mohammadi-Ivatloo, B.; Rabiee, A. Energy hub management with intermittent wind power. In Large Scale Renewable Power Generation; Springer: Berlin/Heidelberg, Germany, 2014; Chapter 16; pp. 413–438. [Google Scholar]
  12. Mavrotas, G. Effective implementation of the ε-constraint method in multi-objective mathematical programming problems. Appl. Math. Comput. 2009, 213, 455–465. [Google Scholar] [CrossRef]
  13. Javadi, M.; Lotfi, M.; Osório, G.J.; Ashraf, A.; Nezhad, A.E.; Gough, M.; Catalão, J.P. A multi-objective model for home energy management system self-scheduling using the epsilon-constraint method. In Proceedings of the 2020 IEEE 14th International Conference on Compatibility, Power Electronics and Power Engineering (CPE-POWERENG), Setubal, Portugal, 8–10 July 2020; pp. 175–180. [Google Scholar]
  14. Kia, M.; Nazar, M.S.; Sepasian, M.S.; Heidari, A.; Siano, P. Optimal day ahead scheduling of combined heat and power units with electrical and thermal storage considering security constraint of power system. Energy 2017, 120, 241–252. [Google Scholar] [CrossRef]
  15. Nazari-Heris, M.; Mehdinejad, M.; Mohammadi-Ivatloo, B.; Babamalek-Gharehpetian, G. Combined heat and power economic dispatch problem solution by implementation of whale optimization method. Neural Comput. Appl. 2019, 31, 421–436. [Google Scholar] [CrossRef]
  16. Nazari-Heris, F.; Mohammadi-ivatloo, B.; Nazarpour, D. Network constrained economic dispatch of renewable energy and CHP based microgrids. Int. J. Electr. Power Energy Syst. 2019, 110, 144–160. [Google Scholar] [CrossRef]
  17. Kopanos, G.M.; Georgiadis, M.C.; Pistikopoulos, E.N. Energy production planning of a network of micro combined heat and power generators. Appl. Energy 2013, 102, 1522–1534. [Google Scholar] [CrossRef]
  18. Alomoush, M.I. Optimal Combined Heat and Power Economic Dispatch Using Stochastic Fractal Search Algorithm. J. Mod. Power Syst. Clean Energy 2020, 8, 276–286. [Google Scholar] [CrossRef]
  19. Carrión, M.; Arroyo, J.M. A computationally efficient mixed-integer linear formulation for the thermal unit commitment problem. IEEE Trans. Power Syst. 2006, 21, 1371–1378. [Google Scholar] [CrossRef]
  20. Alipour, M.; Zare, K.; Seyedi, H. A multi-follower bilevel stochastic programming approach for energy management of combined heat and power micro-grids. Energy 2018, 149, 135–146. [Google Scholar] [CrossRef]
  21. Alipour, M.; Zare, K.; Zareipour, H.; Seyedi, H. Hedging strategies for heat and electricity consumers in the presence of real-time demand response programs. IEEE Trans. Sustain. Energy 2018, 10, 1262–1270. [Google Scholar] [CrossRef]
  22. Wang, J.; You, S.; Zong, Y.; Træholt, C.; Zhou, Y.; Mu, S. Optimal dispatch of combined heat and power plant in integrated energy system: A state of the art review and case study of Copenhagen. Energy Procedia 2019, 158, 2794–2799. [Google Scholar] [CrossRef]
  23. Wang, J.; You, S.; Zong, Y.; Træholt, C.; Dong, Z.Y.; Zhou, Y. Flexibility of combined heat and power plants: A review of technologies and operation strategies. Appl. Energy 2019, 252, 113445. [Google Scholar] [CrossRef]
  24. Bagherian, M.A.; Mehranzamir, K. A comprehensive review on renewable energy integration for combined heat and power production. Energy Convers. Manag. 2020, 224, 113454. [Google Scholar] [CrossRef]
  25. Parvania, M.; Scaglione, A. Unit commitment with continuous-time generation and ramping trajectory models. IEEE Trans. Power Syst. 2015, 31, 3169–3178. [Google Scholar] [CrossRef]
  26. Parvania, M.; Scaglione, A. Generation ramping valuation in day-ahead electricity markets. In Proceedings of the 2016 49th Hawaii International Conference on System Sciences (HICSS), Koloa, HI, USA, 5–8 January 2016; pp. 2335–2344. [Google Scholar]
  27. Li, B.; Bagherinezhad, A.; Khatami, R.; Parvania, M. Continuous-Time Look-Ahead Optimization of Energy Storage in Real-Time Balancing and Regulation Markets. IEEE Syst. J. 2020. [Google Scholar] [CrossRef]
  28. Khatami, R.; Parvania, M. Continuous-Time Locational Marginal Price of Electricity. IEEE Access 2019, 7, 129480–129493. [Google Scholar] [CrossRef]
  29. Bagherinezhad, A.; Khatami, R.; Parvania, M. Continuous-time look-ahead flexible ramp scheduling in real-time operation. Int. J. Electr. Power Energy Syst. 2020, 119, 105895. [Google Scholar] [CrossRef]
  30. Davoodi, E.; Babaei, E.; Mohammadi-Ivatloo, B.; Rasouli, M. A novel fast semidefinite programming-based approach for optimal reactive power dispatch. IEEE Trans. Ind. Inform. 2019, 16, 288–298. [Google Scholar] [CrossRef]
  31. Rong, A.; Hakonen, H.; Lahdelma, R. An efficient linear model and optimisation algorithm for multi-site combined heat and power production. Eur. J. Oper. Res. 2006, 168, 612–632. [Google Scholar] [CrossRef]
  32. Wang, H.; Yin, W.; Abdollahi, E.; Lahdelma, R.; Jiao, W. Modelling and optimization of CHP based district heating system with renewable energy production and energy storage. Appl. Energy 2015, 159, 401–421. [Google Scholar] [CrossRef]
  33. Davoodi, E.; Babaei, E.; Mohammadi-Ivatloo, B.; Shafie-Khah, M.; Catalão, J.P. Multiobjective Optimal Power Flow Using a Semidefinite Programming-Based Model. IEEE Syst. J. 2020, 15, 158–169. [Google Scholar] [CrossRef]
  34. Parvania, M.; Khatami, R. Continuous-time marginal pricing of electricity. IEEE Trans. Power Syst. 2016, 32, 1960–1969. [Google Scholar] [CrossRef]
  35. Dierckx, P. Curve and Surface Fitting with Splines; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  36. Soroudi, A. Power System Optimization Modeling in GAMS; Springer: Berlin/Heidelberg, Germany, 2017; Volume 78. [Google Scholar]
  37. California ISO Open Access Same-Time Information System (OASIS). June 2018. Available online: http://oasis.caiso.com (accessed on 29 April 2021).
  38. Brook, A.; Kendrick, D.; Meeraus, A. GAMS, a user’s guide. ACM Signum Newsl. 1988, 23, 10–11. [Google Scholar] [CrossRef]
  39. Chiandussi, G.; Codegone, M.; Ferrero, S.; Varesio, F.E. Comparison of multi-objective optimization methodologies for engineering applications. Comput. Math. Appl. 2012, 63, 912–942. [Google Scholar] [CrossRef] [Green Version]
  40. Aghaei, J.; Amjady, N.; Shayanfar, H.A. Multi-objective electricity market clearing considering dynamic security by lexicographic optimization and augmented epsilon constraint method. Appl. Soft Comput. 2011, 11, 3846–3858. [Google Scholar] [CrossRef]
Figure 1. Characteristic points (A–E) and approximate heat-power operating region for a CHP unit.
Figure 1. Characteristic points (A–E) and approximate heat-power operating region for a CHP unit.
Sustainability 13 05058 g001
Figure 2. (a) Bernstein-basis polynomials of degree 3; (b) approximation of a continuous function using Bernstein basis polynomial of degree 3.
Figure 2. (a) Bernstein-basis polynomials of degree 3; (b) approximation of a continuous function using Bernstein basis polynomial of degree 3.
Sustainability 13 05058 g002
Figure 3. Real-time, continuous-time, and hourly electrical load curves.
Figure 3. Real-time, continuous-time, and hourly electrical load curves.
Sustainability 13 05058 g003
Figure 4. Day-ahead hourly power schedule for the best compromise solution.
Figure 4. Day-ahead hourly power schedule for the best compromise solution.
Sustainability 13 05058 g004
Figure 5. Day-ahead hourly heat schedule for the best compromise solution.
Figure 5. Day-ahead hourly heat schedule for the best compromise solution.
Sustainability 13 05058 g005
Figure 6. Pareto optimal solutions of the total day-ahead emission versus total day-ahead cost.
Figure 6. Pareto optimal solutions of the total day-ahead emission versus total day-ahead cost.
Sustainability 13 05058 g006
Figure 7. Day-ahead continuous-time power schedule for the best compromise solution.
Figure 7. Day-ahead continuous-time power schedule for the best compromise solution.
Sustainability 13 05058 g007
Figure 8. Day-ahead continuous-time heat schedule for the best compromise solution.
Figure 8. Day-ahead continuous-time heat schedule for the best compromise solution.
Sustainability 13 05058 g008
Table 1. Day-ahead multi-objective scheduling.
Table 1. Day-ahead multi-objective scheduling.
Hourly Model (Study I)Continuous-Time Model (Study II)
Total Day-ahead Cost (USD)4,048,7254,051,994
Real-time Operation Cost (USD)113,18040,227
Total Operation Cost (USD)4,161,9054,092,221
Cost Saving (USD)-69,684 (1.67%)
Total Day-ahead Emission (kg)385,565385,999
Real-time Emission (kg)10,6554443
Total Emission (kg)396,221390,443
Emission Reduction (kg)-57,778 (1.46%)
Computed time (min)0.36481.9433
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Davoodi, E.; Balaei-Sani, S.; Mohammadi-Ivatloo, B.; Abapour, M. Flexible Continuous-Time Modeling for Multi-Objective Day-Ahead Scheduling of CHP Units. Sustainability 2021, 13, 5058. https://doi.org/10.3390/su13095058

AMA Style

Davoodi E, Balaei-Sani S, Mohammadi-Ivatloo B, Abapour M. Flexible Continuous-Time Modeling for Multi-Objective Day-Ahead Scheduling of CHP Units. Sustainability. 2021; 13(9):5058. https://doi.org/10.3390/su13095058

Chicago/Turabian Style

Davoodi, Elnaz, Salar Balaei-Sani, Behnam Mohammadi-Ivatloo, and Mehdi Abapour. 2021. "Flexible Continuous-Time Modeling for Multi-Objective Day-Ahead Scheduling of CHP Units" Sustainability 13, no. 9: 5058. https://doi.org/10.3390/su13095058

APA Style

Davoodi, E., Balaei-Sani, S., Mohammadi-Ivatloo, B., & Abapour, M. (2021). Flexible Continuous-Time Modeling for Multi-Objective Day-Ahead Scheduling of CHP Units. Sustainability, 13(9), 5058. https://doi.org/10.3390/su13095058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop