1. Introduction
The combined heat and power (CHP) system produces heat and electricity with high efficiency by consuming oil, natural gas, and biomass,
etc. [
1]. The well-known cogeneration system has been considered to be the most promising alternative to traditional energy supplying systems. Compared to conventional generation of heat and electricity in a decoupled way, the overall efficiency of the co-generating technique can reach up to 70%–80% [
2,
3]. As a low-carbon, cost-effective, and high energy conversion efficiency technology, the total capacity of CHP is expected to reach 483.7 GW in 2023 worldwide [
4].
Generally, CHPs’ operating strategies are categorized into two groups: heat lead (HL) and electricity lead (EL) [
5]. For HL CHPs, they are operated mainly to satisfy heat load and the heat deficit is imported from auxiliary boilers, district heating or heating grid [
6]. While for EL CHPs, electric loads are satisfied first and the deficient demand is provided from the power grid. The main drawback for the two types of CHPs is the inflexibility between electric and heat output. Heat to power ratio (HTPR) is introduced to characterize the heat and electric output proportion, which is also closely linked with overall efficiency and loading level of CHPs [
7,
8]. Thus, CHPs can be operated at different loading levels to serve various heat and electric loads.
The existing literatures on CHP operation mainly concentrate on achieving economic goals (e.g., low operation cost, reduction of fuel consumption) and environmental goals (e.g., low carbon dioxide emission). The output of CHP is optimized to minimize the annual operating and maintaining costs of the whole system in [
9]. The heat and electric loads are supposed to be known, which might lead to inaccuracy if actual demand diverts greatly from the forecasted demand. The least-cost operating strategy for CHPs is studied in [
10], and results show that the strategy does not mean the lowest carbon dioxide emission. Thus in [
11], the optimal strategy for CHPs is studied by taking several indicators into account, including the least energy consumption, the lowest operation cost and the lowest carbon dioxide emission. The main disadvantage of the mentioned literatures is that energy prices are supposed to be constant for a long period of time. As a result, the CHPs are not optimized in real-time, failing to reflect the actual conditions.
In order to address the problem, various optimization algorithms have been proposed to achieve real-time energy management. A comprehensive review of modeling methods, planning approaches (various indicators are introduced) and energy management algorithms for a combined cooling, heat, and power (CCHP) microgrid are presented in [
12]. An online algorithm is devised in [
13] for real-time energy management. The algorithm is superior to conventional dynamic programming approaches in realization of arbitrary energy prediction errors. An optimal power management strategy for hybrid energy systems is presented in [
14], aiming to minimize the total cost and fuel emissions. In [
15], the optimal size of a CHP-plant with thermal store is analyzed in the Germany spot market to maximize the profits, where only the electricity spot price is considered and the other prices (those of heat and gas) are supposed to be constant. In [
16], a CHP-based district heating system with renewable energy source and energy storage system is optimized by using LP (linear programming) method, where the overall costs of net heat and power acquisition is minimized effectively. A model predictive control strategy is proposed for the CHP to achieve effective demand response in [
17]. The results show the cost reductions for the households reach their highest when the real-time energy prices fluctuate strongly. Reference [
18] proposes a smart hybrid renewable energy for communities (SHREC) system, where both thermal and electricity markets are considered. A planning model is developed for the SHREC system and optimized by using the LP algorithm, whose effectiveness and flexibility are verified through specific calculation results.
All the aforementioned literatures proposed effective approaches to design optimal operating strategy for the CHP from different aspects. While during the research, the real-time energy prices are assumed to have been already known. However, in real-time markets, the prices of electricity, heat, and natural gas are not known in advance. The operation schedules of CHPs are pre-developed based on forecasted prices or load conditions. On the other hand, the influence of loading level on overall efficiency and HTPR of CHPs is not considered either. The overall efficiency is generally in positive correlation with the loading level of CHPs, i.e., a high loading level means high overall efficiency. Another problem for existing researches is that the profit optimization of CHP considering the forecast prices is neglected, where the profits of pre-developed operating strategy differ greatly from those of adjusted strategy in response to actual prices. In deregulation energy markets, the profits are key incentive for CHP owners. Thus, to optimize the profits of CHPs in real-time, there are several challenges: (i) to determine the heat and electricity output proportion, i.e., HTPR; (ii) to forecast energy prices; (iii) to adjust pre-developed operating strategies to comply with actual conditions.
A discrete operation optimization model is developed in this paper to devise the optimal operation strategy for CHPs in real-time, where the profits reach the maximum. First, the SIMO (single input and multi-output) CHP model is established, where both the overall efficiency and HTPR are in variation with the loading level. Then, energy prices, including that of heat, electricity and natural gas, are forecasted by the gray forecasting model GM(1,n) (i.e., first-order, n variables grey model) and revised with actual prices by using the least square method. Finally, a discrete optimization model is developed to obtain the optimal operating points, which may produce maximum profits. The profits are optimized every 30 min, which determine the heat and electricity output in real-time. Then the operating strategy is pre-developed based on the forecasted prices. A dynamic programming algorithm is proposed to adjust the pre-developed strategy after actual prices are known, where the maximum profits are achieved through reducing the profit loss during the modification to the least. The proposed model and optimization approach are demonstrated with a 1MW CHP under different loading levels.
The novelty of this paper is that it: (i) introduces varying HTPR with respect to the loading level of CHP, which is fixed in most previous research; (ii) proposes a discrete optimization model to project the indication optimal operating points which may produce potential maximum profits; (iii) proposes a dynamic programming algorithm to obtain the maximum profits based on actual and forecasted energy prices.
The remaining parts of this paper are organized as follows. In
Section 2, the SIMO model for CHP is established, whose overall efficiency and HTPR both vary with the loading level. The price forecasting method is presented in
Section 3. Then the discrete optimization model for CHP is developed in
Section 4, followed by the dynamic programming algorithm. In
Section 5, a 1 MW CHP is optimized with the proposed approach. Finally, conclusions are drawn in
Section 6.
3. Price Forecasting
When optimizing the operation strategy of CHPs to achieve maximum profits, energy prices should also be determined. In the wholesale market, electric price varies at 30-min resolution. The real-time price of natural gas and heat are supposed to be varying every 30 min as well in the study. Each price interval is also called a dispatching step and therefore there are 48 dispatching steps in one day. Usually, the real-time prices of the kth dispatching step are not known until beforehand. Thus, the operation schedules of CHPs should be pre-determined based on forecasted prices. The real-time energy prices fluctuate with the load of the energy system, the climate and so on. However, there is still strong regularity compared to the history prices for the same time in one day.
For the energy prices are relevant to several factors, the gray forecasting model GM(1,
n) [
27] is adopted in this paper to forecast the energy prices. In the model, “1” means the first-order differential equations are used and “
n” is the number of relevant variables. Suppose the price to be forecasted at
is
, and the historical price series is
, where
is the historical price for different days at
. The relevant variable (including load, climate and
etc.) sequences are:
where
is the historical value series for the
ith variable in
M days at the fixed time
t;
is historical value of the
ith variable in the
mth day.
Then calculate the 1-AGO (accumulated generating operation) sequence for each relevant variable to reduce the randomness of the data:
where
and
,
.
The adjacent mean value sequence for
is calculated as follows:
where
and
,
k = 1, 2, …,
M − 1.
Based on and , is obtained with the model. The forecasting method mentioned above is applied in each dispatching step to get the forecasted prices for heat, electricity, and natural gas in one day.
After the actual prices for heat, electricity, and natural gas are known at the start of
kth dispatching step, the forecasted prices for the (
k + 1)th dispatching step should be revised with certain method. Suppose the forecast prices for
k + 1 dispatching steps are
, which are obtained from the gray forecasting model, and the actual prices for
k dispatching steps are
. The fitting price curves for
and
can then be obtained by using the least squares method, as shown in
Figure 4.
Figure 4.
Fitting price curves.
Figure 4.
Fitting price curves.
In
Figure 4,
is the fitting curve for
and a new forecasted price
(denoted by a small circle) for the (
k + 1)th dispatching step is obtained.
in red is the fitting curve for
. The blue curve
is obtained by using the least square method, whose square of distance to
and
is the minimum. At last, the revised forecasted price
(denoted by a small solid square) is obtained, which is in the curve
.
4. Operation Optimization of Combined Heat and Power
4.1. Discrete Optimization Model for Combined Heat and Power
In the wholesale market, electricity, heat and natural gas prices all vary at 30-min resolution, denoted by , and , respectively. Thus, CHPs could be operated according to the combination of the three prices to maximize benefits in the 48 dispatching steps.
In the
dispatching step, the profit is calculated by:
where PRO is abbreviation of profit and
denotes the profit in the
kth step;
and
are the income from selling heat and electricity, respectively;
is the cost of buying natural gas.
Equation (12) can be further written as follows:
where
is the volume of natural gas consumed in the
dispatching step;
,
and
are the prices of heat, electricity and natural gas in the
dispatching step respectively.
Usually, the energy contained in a cubic meter of natural gas is a constant, denoted by
, in kJ/m
3. Thus the total energy injected into the CHP in the
step is expressed as:
Through substituting Equation (6) into Equation (7),
and
are obtained and shown as follows:
By substituting Equations (5), (14) and (15) into Equation (13), the objective function of CHP is obtained, given by:
To solve the objective function Equation (16), we define a new variable
, which means the profit obtained through consuming a cubic meter of natural gas by CHP:
Obviously, as long as reaches its maximum, the maximum value of will be obtained. In Equation (17), after , , , and are all determined in the dispatching step, is just a function of η, which is quite easy to solve.
There are also a set of constraints which must be satisfied during optimization.
4.1.1. Output Constraints
The heat and electricity outputs of CHP must vary between their maximum and minimum capacities, respectively:
where
and
are allowable minimum and maximum heat outputs of CHP;
and
are the allowable maximum electricity outputs of CHP.
4.1.2. Input Constraints
The input of CHP is the natural gas only and the input should be limited between its minimum and nominal value due to capacity limits of CHP:
where
is the rated capacity of CHP;
and
are the allowable minimum and nominal input volume of natural gas respectively.
4.1.3. Ramp Constraints
The output power of CHP between two consecutive dispatching steps cannot be modified randomly, which is constrained by a certain ramp rate of CHP. The whole process is modeled by:
where
and
are the useful energy converted from natural gas by CHP in the
and
dispatching step, respectively;
is the ramp capacity of CHP, in kW/min;
is the discrete time step length of 30 min.
By submitting Equations (4) and (15) into Equation (20), the heat and electricity outputs ramp constraints of CHP are:
4.2. Solution for the Model
In actual application of CHPs, η is described with a series of data from field test. A discrete method is proposed to solve the optimization model. In this paper, the typical loading levels in
Figure 2 and
Figure 3 are shown in
Table 1.
Table 1.
Typical loading levels in Figure 2 and Figure 3.
Typical Loading Levels in Figure 2 | Typical Loading Levels in Figure 3 |
---|
| 40 | | 40 |
| 65 | | 60 |
| 78 | | 80 |
| 90 | | 90 |
| 100 | | 100 |
By combining
Figure 2 and
Figure 3, the relation between
and η is rather clear, as shown in
Figure 5. In
Figure 5, there are mainly seven intervals where CHP could be operated at corresponding efficiency and HTPR. The seven intervals are shown in
Table 2.
Figure 5.
Optimization intervals division.
Figure 5.
Optimization intervals division.
Table 2.
Operation intervals for CHP. Heat to power ratios: HTPR.
Table 2.
Operation intervals for CHP. Heat to power ratios: HTPR.
Intervals | Overall Efficiency | HTPR | Interval type |
---|
I | f1(L) | ζ1 | Zero |
II | f2(L) | ζ2 | Non-zero |
III | f2(L) | ζ3 | Non-zero |
IV | f3(L) | ζ3 | Non-zero |
V | f4(L) | ζ3 | Non-zero |
VI | f4(L) | ζ4 | Non-zero |
VII | f5(L) | ζ5 | Non-zero |
In
Table 2, there are one zero interval (Interval I) and six non-zero intervals (from Interval II to Interval VII). In the
dispatching step, the forecasted
,
and
are all constant. Through solving Equation (17) in the six non-zero intervals, respectively, the maximum profit for each interval will be obtained.
Take the case in Interval II as an example. In this interval, a fixed step of
L is adopted, which is denoted by
. Then a series profit can be obtained, shown as:
The results obtained from Equation (22) failing to satisfy Equations (18)–(21) are deleted first, and then the maximum profit for Interval II is obtained and denoted by .
The steps of L in each interval are usually different, which is determined by the loading level resolution of the field test results. The maximum profits for the other intervals are calculated similar to that of Interval II.
Then the maximum profit of the
dispatching step can be expressed as:
where
is the maximum value of
in the
dispatching step;
is the maximum value of
in the
interval.
4.3. Dynamic Programming
Based on the forecasted prices of heat, electricity, and natural gas, the indication optimal operation points for CHP in one day can be obtained and corresponding operating schedules are pre-developed. Since the actual prices for each dispatching step are known at its beginning, certain modification has to be made to the operating schedule to obtain the actual maximum profits for CHP.
From Δ
T to 2Δ
T, the optimal operation points with maximum profits are obtained first based on forecasted prices. Around Δ
T, the loading level is to be modified to the optimal loading level of the following dispatching step. The whole process is presented in
Figure 6.
is the optimal loading level obtained based on the actual prices of 0–Δ
T, where the profit of CHP is the maximum. The loading level is modified from 70% to
at the ramp rate of CHP. At Δ
T, the actual prices for Δ
T–2Δ
T are known immediately and the corresponding optimal loading level is
. The load has the characteristic of retaining unchanged during a short period, which is reflected as the retaining characteristic of price. Thus, the CHP should be modified to the target optimal operation point when it is different from the current one.
Figure 6.
Dynamic modification of CHP based on actual and forecasted prices.
Figure 6.
Dynamic modification of CHP based on actual and forecasted prices.
For different relationship between the actual prices and forecasted prices, the actual maximum loading level may be higher or lower than the forecasted maximum loading level, named as Case 1 and Case 2, respectively.
4.3.1. Case 1
In Case 1, the forecasted optimal loading level
, obtained from the revised forecasted prices, is lower than
. There are two potential operation routes:
- ■
In the first route, is modified to in advance, and then modified to at ΔT directly.
- ■
In the second route, CHP is operated at until ΔT, and then modified to directly.
To decide which route is better, the profit during the two modification processes should be calculated first. The two routes are denoted by and , respectively.
(1) Route
In 0–Δ
T, the pre-modification time Δ
T’ can be calculated as follows:
where
is the ramp rate and
.
At Δ
T, the loading level is modified to
directly and ends at
F. The overall profit during the modification process from
A’ to
F is:
where
is the profit-loading curve in 0–Δ
T,
is the profit-loading curve in Δ
T–2Δ
T.
(2) Route
The loading level keeps unmodified until
D and modified to
. During the whole process, the forecast optimal loading level
does not affect the operating schedule. The overall profit is:
Based on and , the optimal modification route with the maximum profit can then be determined.
If , the CHP is modified according to .
If , the CHP is modified according to .
4.3.2. Case 2
In Case 2, the forecast optimal loading level is
, which is also obtained from the revised forecasted prices, is higher than
. There are also two potential operation routes.
- ■
In the first route, is pre-modified to , and then decreased to at ΔT.
- ■
In the second route, CHP is operated at until ΔT, and then modified to . The second route is almost the same with that in Case 1.
The two routes are denoted by and , respectively.
(1) Route
In 0–Δ
T, the pre-modification time Δ
T’’ is determined as follows:
The actual prices for Δ
T–2Δ
T are known at
B. Then the loading level of CHP is modified to
G and ended at
F. The overall profit for the modification process is expressed as:
(2) Route
The route is almost same with that in Case 1, with a different beginning point
A. The overall profit for the modification process is calculated by:
Based on and , the optimal modification route can be obtained.
If , the CHP is modified according to .
If , the CHP is modified according to .
The same way is applied to each 2ΔT time window, which starts at and ends at , to obtain the optimal modification route for one day. During the whole process, the price forecasting is considered, which is more practical in real-time markets.
5. Demonstration Examples
In this paper, a 1 MW CHP is optimized based on the proposed method. The other technical parameters of the CHP are presented in
Table A1 in the appendix.
The prices of electricity, heat, and natural gas are forecasted by using the gray forecasting model and revised in real-time with the least square method, where the actual prices are from [
28]. The results are shown in
Figure 7.
As shown in
Figure 7, the energy prices are forecasted accurately by using the proposed method. The variances between the actual and forecasted price of electricity, heat and gas are 0.023, 0.00011, and 0.00018, respectively. In
Figure 7a, the forecasted price is almost the same with its actual value when the price curve is flat, between 0:00 and 7:00, 16:00 and 19:30, 21:00 and 24:00 When the price changes suddenly, which reflects the sudden change of the electric loads, there will exist certain error between the actual and forecasted prices. The case happens at 10:00 and 13:00. The similar phenomenon also occurs in
Figure 7b,c.
Figure 7.
Actual and forecasted prices of electricity, heat, and natural gas. (a) Electricity price; (b) heat price; and (c) natural gas price.
Figure 7.
Actual and forecasted prices of electricity, heat, and natural gas. (a) Electricity price; (b) heat price; and (c) natural gas price.
Based on the actual and forecasted prices, the optimal operation points, where the profits of CHP reach their maximum, are obtained and shown in
Figure 8. At 0:00, The CHP is operated at the loading level of 70%, which is obtained from the day-ahead operating data. In
Figure 8, the profits of CHP operated at different loading levels at different time in one day are presented. The white dotted line is the indication operation points for CHP in one day, which is obtained based on actual prices.
From
Figure 8, it can be found that the maximum profit occurs at 10:00, about 12 Euro shown in deep red, when the CHP is in full load. The relatively high profit for CHP is between 8:00 and 16:00, 20:00 and 21:00. From
Figure 7, it is found the electric price is much higher than that of heat in between 8:00 and 16:00 and, thus, the profit of CHP is relatively higher than other time. Between 0:00 and 8:00, 16:00 and 20:00, when the electric price is the same with or lower than heat price, the profits for CHP are relatively low. So the profit of CHP is mainly determined by the electric price, for its usually high value than heat.
Figure 8.
Profit distribution of the CHP.
Figure 8.
Profit distribution of the CHP.
Based on the indication operation points, specific operation schedule can be developed after considering the ramp constraints. The green dotted line in
Figure 8 is the forecasted operating route obtained from the forecasted prices. The red line in
Figure 8 in is the optimal operating route obtained from the proposed dynamic programming method. It can be found that the red line and the white dotted line almost coincide with each other. Namely, the optimized profits comply quite well with those derived from actual conditions. At 1:30, the loading level for maximum profit point got from the forecasted prices is 90%, and it is modified to 90% before 1:30. However, the loading level for the maximum profit point got from the actual price is 80% at 1:30, thus the loading level has to be modified back to 80%. During the process, a lot of profits are lost for CHP. The optimal operating route is staying unchanged, as shown in the red line. At 2:30, the operating route is modifying CHP before 2:30 in the green dotted line, while the optimal route is modifying CHP after 2:30, shown in the red line. Sometimes, the operating route obtained from the forecasted prices happen to be the optimal one, such as at 6:30, 16:30, 19:00, 21:30, and 22:30. Take the case at 6:00 as an example, the loading level of CHP should be modified to 100% before 6:30, and the results of the green dotted line are consistent with that of the red line. The profits of the routes in green and red are shown in
Table 3, respectively.
From
Table 3, it can be found that the profits for CHP of the optimal route reach their maximum. The maximum profit in one day is about 14 Euro, between 10:00 and 11:00. The maximum profit is just about 1 Euro in the deep night. The high loading level for CHP does not mean high profits, either. Between 0:00 and 2:30, 2:30 and 6:00, the CHP is operated at the loading level of 80% and 90% to obtain the maximum profits. The similar case occurs between 16:00 and 18:00, 21:00 and 24:00. The reason is that the price of electricity is low and the CHP is operated at relatively low loading level to gain a high value of HTPR.
Table 3.
Maximum profits of the forecasted route and optimal route.
Table 3.
Maximum profits of the forecasted route and optimal route.
Time | Maximum Profit from the Forecasted Route (Euro) | Maximum Profit from the Optimal Route (Euro) | Time | Maximum Profit from the Forecasted Route (Euro) | Maximum Profit from the Optimal Route (Euro) |
---|
0:30 | 1.01 | 1.01 | 12:30 | 6.11 | 6.11 |
1:00 | 1.05 | 1.07 | 13:00 | 12.37 | 12.37 |
1:30 | 1.00 | 1.03 | 13:30 | 10.13 | 10.13 |
2:00 | 0.97 | 0.97 | 14:00 | 7.30 | 7.30 |
2:30 | 0.95 | 0.95 | 14:30 | 7.24 | 7.24 |
3:00 | 0.95 | 0.95 | 15:00 | 5.93 | 5.93 |
3:30 | 0.95 | 0.95 | 15:30 | 4.52 | 4.52 |
4:00 | 1.01 | 1.01 | 16:00 | 2.38 | 2.38 |
4:30 | 1.04 | 1.04 | 16:30 | 2.05 | 2.05 |
5:00 | 1.01 | 1.03 | 17:00 | 0.98 | 0.99 |
5:30 | 0.97 | 1.00 | 17:30 | 1.05 | 1.08 |
6:00 | 0.92 | 0.95 | 18:00 | 1.12 | 1.15 |
6:30 | 0.81 | 0.81 | 18:30 | 1.54 | 1.54 |
7:00 | 0.81 | 0.81 | 19:00 | 1.91 | 1.91 |
7:30 | 1.98 | 1.98 | 19:30 | 1.74 | 1.74 |
8:00 | 2.67 | 2.67 | 20:00 | 3.88 | 3.88 |
8:30 | 6.71 | 6.71 | 20:30 | 4.64 | 4.64 |
9:00 | 8.15 | 8.15 | 21:00 | 1.94 | 1.94 |
9:30 | 12.09 | 12.09 | 21:30 | 1.27 | 1.27 |
10:00 | 14.05 | 14.05 | 22:00 | 1.18 | 1.18 |
10:30 | 14.03 | 14.03 | 22:30 | 0.99 | 0.99 |
11:00 | 10.70 | 10.70 | 23:00 | 0.88 | 0.90 |
11:30 | 8.29 | 8.29 | 23:30 | 0.91 | 0.93 |
12:00 | 5.74 | 5.74 | 24:00 | 0.38 | 0.38 |