Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6287291 |

Xian-hua Gao, Zhi-gang Su, "Artificial Bee Colony Optimization of NOx Emission and Reheat Steam Temperature in a 1000 MW Boiler", Mathematical Problems in Engineering, vol. 2019, Article ID 6287291, 13 pages, 2019.

Artificial Bee Colony Optimization of NOx Emission and Reheat Steam Temperature in a 1000 MW Boiler

Academic Editor: Kauko Leiviskä
Received03 Jul 2019
Revised10 Sep 2019
Accepted30 Sep 2019
Published11 Nov 2019


This paper puts forward a new viewpoint on optimization of boiler combustion, namely, reducing NOx emission while maintaining higher reheat steam temperature rather than reducing NOx emission while improving boiler efficiency like traditional practices. Firstly, a set of multioutputs nonlinear partial least squares (MO-NPLS) models are established as predictors to predict these two indicators. To guarantee better predictive performance, repeated double cross-validation (rdCV) strategy is proposed to identify the structure as well as parameters of the predictors. Afterward, some controllable process variables, taken as inputs of the predictors, are then optimized by minimizing NOx emission and maximizing reheat steam temperature via multiobjective artificial bee colony (MO-ABC). Results show that our rdCV-MO-NPLS model with MO-ABC optimization methods can reduce NOx emission synchronously and improve reheat steam temperature effectively compared with nondominated sorting genetic algorithm II (NSGA-II) and combustion adjustment experimental data on a real 1000 MW boiler.

1. Introduction

Boiler combustion optimization of reducing nitrogen oxide (NOx, for short) has become a hot topic during the past three decades. There exist several currently typical ways in coping with the boiler combustion optimization problem. The first popular method is numerical simulation that is based on computational fluid dynamics (CFD) [13]. However, the optimal set points explored by numerical simulation cannot guarantee optimal value in practice due to inconsistencies between designed conditions in CFD and practical conditions. To consider this, it is interesting to take into account the abundant information originating from practical running data. Under such circumstance, the data-driven method is motivated.

The development of data-driven methods can be divided into two stages. At the first stage, researchers only focus on NOx emission reduction by using evolutional algorithms (EAs) to search the optimal set point of controllable process variables. The typical EAs include genetic algorithm (GA), particle swarm optimization (PSO), and ant colony optimization (ACO). Among them, GA is frequently used on boiler combustion optimization. For example, in [4, 5], artificial neural network (ANN) combined with GA is proposed to reduce NOx emission. And in [6], support vector machine (SVM) and conventional GA (CGA) are used on NOx emission reduction. Researchers can notice that only NOx emission is considered as the optimization objective in the above researches, but ignoring that the combustion system is a hybrid system which is coupled with other subsystems. The reduction of NOx usually plays great influence on other important parameters such as reheat steam temperature (Tr, °C). That is, the optimization of boiler combustion system is essentially a multiobjective optimization problem.

This issue with multiobjective optimization in boiler combustion systems has gained great attention at the second stage. Traditionally, boiler combustion optimization is implemented to reduce NOx emission and increase boiler efficiency by EAs (for example, see works [710]). ANNs with multiobjective GA (MOGA) are employed on boiler combustion optimization. Although the above methods have made significant contributions on reducing NOx emission while improving boiler efficiency, there still exist some unconsidered issues. Firstly, reheat steam temperature is always lower than its set value due to the requirement of meeting electric peak regulation. This phenomenon not only imposes a bad influence on the whole plant efficiency but also poses a threat to safe control of the boiler combustion system. Secondly, the precise calculation of boiler efficiency is hard to complete actually due to some practical issues. The third problem is the lack of effective modelling data, as historical running data are heavily coupled and incomplete. Therefore, it is urgent to improve reheat steam temperature.

Up to date, to the best of our knowledge, there is few literature on boiler combustion optimization to reduce NOx emission and improve reheat steam temperature simultaneously. Motivated by the above statements, this paper aims to reduce NOx emission while achieving higher reheat steam temperature by proposing a jointed optimization method. More precisely, a Multiple Output version of Nonlinear Partial Least Squares (MO-NPLS) regression model is first proposed to predict behaviors of NOx emission and reheat steam temperature synchronously. To achieve better performance, the structure and parameters of MO-NPLS are identified by repeated double cross validation (rdCV) rather than traditional k-fold cross-validation strategy. Then, a multiobjective artificial bee colony (MO-ABC) algorithm is adopted to search optimal inputs of predictors that can achieve lower NOx emission and higher reheat steam temperature from Pareto front. The modelling data are from a combustion adjustment experiment conducted on a real 1000 MW coal-fired boiler. It will be seen that our proposed method outperforms combustion adjustment experimental data and commonly used nondominated sorting genetic algorithm II (NSGA-II). Therefore, the main contributions of this paper can be summarized in three-fold:(1)Both NOx emission and reheat steam temperature are firstly proposed in boiler combustion optimization. This strategy can reduce NOx emission while maintaining higher reheat steam temperature to guarantee the safety in boiler combustion system control as well as the improvement of whole plant efficiency.(2)A MO-NPLS method with rdCV strategy is raised to establish the predictive model of both NOx emission and reheat steam temperature. This algorithm is powerful due to its strong generalization ability especially with small number of samples.(3)A MO-ABC-based multiobjective optimization method is proposed on boiler combustion optimization. On the basis of the rdCV-MO-NPLS modelling method, our proposed MO-ABC algorithm can reduce NOx emission synchronously and improve reheat steam temperature effectively compared with NSGA-II and the experimental data.

The rest of the paper is organized as follows. In Section 2, the 1000 MW boiler combustion system and data sources are described. Section 3 establishes predictors used to predict behaviors of NOx emission and reheat steam temperature by rdCV-MO-NPLS. Then, Section 4 presents the joint optimization method by MO-ABC and the comparison with NSGA-II and experimental data. Section 5 concludes this paper.

2. Description of Combustion System

In this section, the boiler combustion system is introduced, together with the dynamic characteristics of combustion system. And then, the data source used for modelling establishment and validation is presented.

2.1. Boiler Combustion System

Figure 1 illustrates a 1000 MW tower-type ultrasupercritical boiler adopting the reheat system and air-staged combustion technique. As can be seen from Figure 1, a selective catalytic reduction (SCR) system is installed at the exit of convection flue gas and two air preheaters are arranged beneath the SCR system. Along vertical direction, the entire furnace is divided into three zones including main combustion zone, reduction zone, and burning-out zone. Furthermore, in the main combustion zone, twelve-level pulverized coal burners are equipped and correspond to six mills indexed by letters from A to F that are with burners A1 and A2 at the bottom to F1 and F2 at the top, respectively.

Generally, different coal allocations decided by the coal feeding rate of each mill will lead to different temperature fields and flame centers in the furnace. For instance, if more coal is supplied by bottom mills, the flame center will move downward. As a result, NOx production will reduce due to the extension of reductive zone; however, the reheat steam temperature (Tr) will decrease too because of less quality of heat. In the present work, although the coal feeding rate is a directly manipulated parameter, thermal load () is monitored instead by considering the diversity of coal types as follows:where is the total coal fed by mills and indicates the total thermal load.

In contrast to coal allocations, primary air is usually not adjusted because its main task is to transport pulverized power from mills, i.e., to match the pulverized coal. Secondary air, including that in the main combustion zone and separated overfire air (SOFA) in the burning-out zone, can be adjusted to improve combustion conditions in the boiler. To enhance reducibility in the main combustion zone so as to reduce NOx production, one can decrease secondary air in the main combustion zone with increase of SOFA ratio () as the total air volume remains constant. Evidently, this parameter represents the ratio of secondary air to total air. More interestingly, higher reheat steam temperature can be achieved with higher , but excessive will affect the combustion performance in the main combustion zone. However, the volume of secondary air in the main combustion zone is rather difficult to monitor due to tens of dampers in power plant. Instead, is used and defined aswhere total air = primary air + secondary air. Considering that there is about 20% running oxygen () in total air. Therefore, the higher the running is, the greater the total air volume will be.

In addition, notice that the higher the running value is, the better combustion effect will be. Then, both NOx emission and reheat steam temperature will increase subsequently and vice versa. This means there exists an optimal α value, namely, optimal running .

In summary, to reduce NOx emission and increase reheat steam temperature synchronously heavily relies on the distribution of thermal loads (), running , and in the furnace.

Remark 1. The burners at vertical direction are tilted within a range of ±30° to tune the reheat steam temperature () [11]. At low load condition, even if burners have been already fully titled up, the reheat steam temperature will be still significantly lower than its desirable setpoint. This will impose a bad influence on whole power plant efficiency (the whole plant efficiency where is plant thermal efficiency which depends on reheat steam temperature).

2.2. Data Source

The value of reheat steam temperature from actual operational conditions is not accurate as it is controlled by a closed-loop temperature system. Therefore, we adopted the data collected from a combustion adjustment experiment shown in Supplementary file of the 1000 MW power plant described in Section 2.1. Consider the problem of excessive NOx emission and lower reheat steam temperature is more severe under lower loads. As a result, only 650 MW, 800 MW, and 850 MW conditions are discussed in the following. Corresponding data are listed in Appendix A (Tables A1–A3) in the Supplementary file.

Remark 2. Reheat steam temperature in the experiment was not controlled by a closed-loop way, and hence, its values are real. Therefore, all collected experimental data are more useful than data collected in practical running way when validating the performance of our method that will be presented later.

3. Estimations of NOx and Reheat Steam Temperature Based on MO-NPLS with rdCV

PLS/NPLS algorithms have been applied in engineering (for example, see [1214]). However, due to the limited number of samples, the prediction accuracy of PLS/NPLS based on traditional k-fold cross validation may decrease. Motivated by this and literature [15], we propose a multiple output version of NPLS (MO-NPLS) model with repeated double cross-validation (rdCV) strategy in Section 3.1. The predictive model is necessarily used to predict the behaviors of NOx emission and reheat steam temperature and then identify and validate it in different cases of loads in Section 3.2.

3.1. MO-NPLS Model with rdCV

Give inputs and outputs with . Then, let and be the number of n observations for inputs and outputs, respectively. Evidently, indicates the observation of the input, and this is held on for output Y. The superscript “T” represents transpose. To make the results independent from units, all attributes of data are normalized into range according to the z-score rule and are denoted by with and .

Define and are the principal components of X and Y, respectively. Similarly as traditional nonlinear PLS [16, 17], rdCV-MO-NPLS is in nature a recursive modelling strategy, aiming to establish relationship between and firstly and then implement the regression relationship of X and Y.

Let be any column of at the beginning. Then, one can estimate of the first components (, ) as follows.

Firstly, is calculated bywhere is the weight of matrix and obtained by .

Then, one can fit the polynomial relationship between and by the least squares method. Specifically, assume the relationship between and is quadratic, namely,where is the residual vector which is independent and identically distributed. is a vector of ones. Take and , then (4) can be rewritten by a compact form:

So, one can get the polynomial coefficient by the least squares method:where () represents the inverse matrix of .

After coefficients are fitted by formula (6), then , i.e., the estimate of , can be calculated by

Then, update (mark ) bywhere are weights of matrix and determined by . Next, needs to be updated according to [18].

At this point, we can recalculate the vector (mark ) by formula (3) with the updated .

Check the convergence of : if , refit the polynomial coefficient between and . Otherwise, calculate the regression coefficient by

Finally, the residual of and (mark and ) that will be used in the next iteration can be induced bywith regression coefficient .

Replace and with and and formulas (3)–(10) are repeated until all m principal components () with are extracted.

After all m principal components have been extracted, optimum number () can be selected. Usually, k-fold strategy is applied, whereas rdCV strategy is adopted here for achieving more robust predictive ability, which has been proven of PLS in [15]. This is the reason why we call our model rdCV-MO-NPLS. Once is determined, the prediction of outputs Y can thus be obtained according towhere means point-product operation of two vectors with equal dimensionality; with ; and with .

From the above interpretations, the rdCV-MO-NPLS is summarized in Algorithm 1.

Input: data and parameters ;
Output: optimal number of principal component and rdCV-MO-NPLS model.
for to do
 Split dataset S equally into groups, i.e., ;
for to do
  Select as Test set, and construct Calibration sets  = ;
  Split Calibration set equally into groups, ;
  for to do
   Select as Validation set and as Training set;
   Fit MO-NPLS based on with, respectively, of principal components,
   Apply the MO-NPLS models to and get predictive for ;
   Calculate mean square error , with the number of objects in the used validation set     and the output objects of validation set ;
  Estimate optimum principal components , …, according to based on standard error method [15] (Here, more than   one may be selected (i.e., ) because of different confidence interval);
  Make MO-NPLS models based on Calibration set with ;
  Test fitted models on and obtain a group of k predictions as well as k biases;
  Find the smallest bias and determine the optimal principal component .
 One can get with the number of after completing the outer loop.
Totally, after a complete rdCV run, we can get with the number of (). The final optimum of principal component is the one with highest frequency in .
Identify rdCV-MO-NPLS model with database S and to get model parameters [, p, c, q] based on (3)–(10).
3.2. Model Identification and Validation
3.2.1. Model Identification

In this section, rdCV-MO-NPLS is identified, respectively, at 650 MW, 800 MW, and 850 MW. In each case, identification as well as rdCV data are shown in Table A1 (from lines 1 to 45), Table A2 (lines 1 to 30), and Table A3 (lines 1 to 90). In cases of 650 MW and 850 MW, six inputs were considered, including C-F level thermal loads ( to , MW), running (), and (, ), whereas seven inputs were considered in 800 MW, which are B-F level thermal loads ( to , MW), running (), and (, ). In all cases, the outputs are

Figure 2 displays the results of applying rdCV. It can be seen that the highest percentage in three cases are, respectively, 0.32 (), 0.38 (), and 0.256 ( = 2). Hence, the final optimum principal component in each case is 3, 2, and 2, respectively. With the help of , the rdCV-MO-NPLS model can be identified as Case 1 to Case 3.

For better comparison, 5-fold cross validation (We call 5-fold-MO-NPLS model) is conducted by the same identification data as the rdCV method. After 5-fold cross validation process, the optimum principal component of 650 MW, 800 MW, and 850 MW cases is 4, 1, and 4, respectively. Then, the validation results of 5-fold-MO-NPLS model with their will be seen in the next section.

Case 1. (rdCV-MO-NPLS of 650 MW load condition). Since  = 3, the first three components (, ), (, ), and (, ) are selected. Normalize the training data presented in Table A1 (from lines 1 to 45) into and according to rule. Corresponding parameters of () shown in Table 1 as well as (, ), (, ), and (, ) are identified based on (3)–(11). The rdCV-MO-NPLS predictive model is in form of




Case 2. (rdCV-MO-NPLS of 800 MW load condition). The first two components (, ) and (, ) are chosen. Based on data as shown in Table A2 (lines 1 to 30), and identified parameters of rdCV-MO-NPLS shown in Table 2, the rdCV-MO-NPLS model can be established as



Case 3. (rdCV-MO-NPLS of 850 MW load condition). Similar as Case 2, the first two components (, ) and (, ) were selected on the basis of Table A3 (from lines 1–90). Additionally, identified parameters of rdCV-MO-NPLS are shown in Table 3. The rdCV-MO-NPLS model is



3.2.2. Model Validation

In general, generalization ability of a prediction model is more important than its fitting ability. Therefore, it is necessary to validate aforementioned models to get their generalization ability. Validation samples in three cases are seen in Table A1 (from lines 46 to 53), Table A2 (from lines 31 to 36), and Table A3 (from lines 91 to 96) in Supplementary file, respectively.

To state the superiority of the rdCV-MO-NPLS model, another three well-known methods: 5-fold-MO-NPLS, support vector regression (SVR), and artificial neural network (ANN) are employed to establish models of NOx emission and reheated steam temperature. All methods are on the same identification and validation data. Figures 35 show the predictive results for different methods. The prediction error sum ( with prediction sample number ) representing generalization ability and simulation error sum ( with simulation sample number ) reflecting approximate ability is presented in Table 4.


650 MWrdCV-MO-NPLS0.13460.01930.15390.84580.11380.9596

800 MWrdCV-MO-NPLS0.09850.01400.11250.53680.11390.6507

850 MWrdCV-MO-NPLS0.09310.01410.10722.44690.25172.6986

1The values in bold are the optimal in their columns. 2; .

It can be seen from Figures 35 that the prediction of rdCV-MO-NPLS approaches the real values the most while the deviation among 5-fold-MO-NPLS, SVR, or ANN and real value is larger in the three working load conditions. Furthermore, from Table 4, it can be concluded that the value of rdCV-MO-NPLS are lower than the other three methods in all the three cases, which means rdCV-MO-NPLS has stronger predictive capability. This verified that firstly, rdCV can achieve more robust ability than the k-fold strategy; secondly, SVR can only establish the SISO model, as the boiler combustion system is strongly coupled, so the accuracy of predictive model will be reduced due to decoupling. Finally, as for ANN, it can deal with the MIMO system but is prone to overfitting easily.

Remark 3. Compared with k-fold cross validation, rdCV can achieve more robust ability. This also shows that rdCV can be successfully used in NPLS. Besides, rdCV-MO-NPLS is more powerful than SVR and ANN in generalization ability on the modelling of NOx emission and reheat steam temperature.

4. Combustion Optimization

Artificial bee colony (ABC) is a well-known intelligent optimization algorithm proposed in 2005 [19]. It has been applied more and more on engineering optimization in recent years [2023] because of its simplicity to implement as it uses fewer control parameters [24, 25] and robust ability to get out of a local minimum [26]. In this section, a multiobjective version of ABC (MO-ABC) [25] is used to optimize NOx emission and reheat steam temperature.

4.1. Joint Optimization via Multiobjective Artificial Bee Colony Algorithm

The objective of combustion optimization is to minimize NOx emission and maximize reheat steam temperature (Tr) through searching optimal inputs (x) which is m-dimensional, namely, thermal loads (), running , and based on rdCV-MO-NPLS models established in Section 3. Considering that these two objectives are contradictory, we can only find a set of nondominated solutions. Objective function can be described as follows:where is the NOx emission. is the inverse of Tr by converting the problem of maximizing Tr to minimize its inverse. The values of and the inverse of of 650 MW, 800 MW, and 850 MW can be calculated from formulas (13)–(15), respectively. and are the lower and upper bounds of inputs in the rdCV-MO-NPLS model.

To solve problem 4.1, a MO-ABC is applied step by step as follows:

Step 1. (solutions initialization). New variables named and FoodNumber will be defined as the population size and the number of food sources, respectively. Since each employed bee represents a food source and the number of employed bees and onlooker bees is equal, the FoodNumber can be set to half of the population (i.e., ) [25]. Each food source will be initialized through a function randomly generated in input space according towhere i is the index of the food source with the range of FoodNumber; d is the dimension index of each food source with the range of ; rand [0,1] is a random number distributed uniformly over the interval [0,1]; and and are the lower and upper bounds of the dth dimension, respectively.
In addition, assign each food source with a trial index , FoodNumber; is initialized to 0 for each food source.

Step 2. (solutions evolve in employed bee phase). The initial solutions in (17) first evolved through the help of employed bees according towhere and and . is a random number distributed uniformly between [−1, 1]. If dominates , i.e., and (denoted as ), update to . If not, the update was unsuccessful, the trial for the old food source, , is incremented by one.

Step 3. (solutions evolve in onlooker bee phase). The quality of solution , i.e., fitness of solution , is evaluated bywhere is the number of solutions that are dominated by food source . The probability for each food source advertised by the corresponding employed bee will be calculated in a roulette way according toOnce a solution is selected, this solution will evolve to a new position according to formula (18) if this new position dominates its previous one. Otherwise, its trial value is increased by one.

Step 4. (solutions evolve in scout bee phase). There is at most one scout bee in the colony. This means if the maximal value in the colony reaches the limit, then the corresponding food source will turn to a scout bee and after doing a random search according to formula (17). At the same time, its value will be reset to zero. The limit is determined by

Step 5. (solution archive updates). A fixed-size archive (e.g., ) is conserved to hold the best solutions. Specifically, after the th cycle, we can get the nondominated solutions. Combine them with solutions which are held in archive at th cycle. Then, select the top solutions by nondomination strategy to update archive.
All the above steps are repeated until the maximum iteration, , is met. The main process of MO-ABC is summarized in Figure 6.

Remark 4. The specific values of and are shown in Appendix B (Tables B1B3) in Supplementary part, and it is determined by history running and experimental data as well as security conditions as the usual practice of data-driven methods (for example, see work [4]).

4.2. Optimizing Results and Analysis

To implement MO-ABC as shown above, some parameters should be preset including and . Here, was set to 100 by a trivial way. was set to 300 (at 650 MW and 850 MW) and 350 (at 800 MW), respectively, based on formula (21).

Furthermore, another popular multiobjective EA (MOEA), NSGA-II, is also conducted for comparison. Values of and [] in NSGA-II are set the same as MO-ABC.

In the present work, we use the number of function evaluations () as a measure criterion to determine . Set that both algorithms are terminated after 100,000 . Define and as the maximum iteration of NSGA-II and MO-ABC algorithms, respectively. Since , then we can get  = 1000. In addition, the upper bound of the fitness evaluations in basic ABC is defined and given as . In each iteration, there are employed bees, onlooker bees, and at most 1 scout bee [27]. In this paper, substitute “” into above equation, then we can get .

The optimization results of MO-ABC and NAGA-II are presented in Figures 79.

To assess the solution sets shown in Figures 79, a comprehensive quality indicator (), i.e., hypervolume (HV) [28], is applied here. The HV indicator represents the volume of the area enclosed by solution set and a reference point, which can combine the quality of convergence, spread, cardinality, and uniformity of solution set and can be defined as follows.

Given a solution set A and a reference point r, HV can be calculated aswhere λ is the Lebesgue measure. Put it simply, the HV value of a set can be seen as the volume of the union of the hypercubes determined by each of its solutions and the reference point (as the left-bottom vertex and the right-top vertex, respectively).

Specifically, for a two-objective optimization problem, HV can be calculated by the following steps. Firstly, the nondominated points in set A are sorted by descending order on their first dimensional objective values. And then, HV is calculated bywhere represents the value of the jth solution in A on their ith dimensional objective. is set to . n is the number of solutions in set A.

The r and HV are illustrated in the third and fourth column of Table 5.

LoadsMOEAsrHVTime consuming (s)

650 MWMO-ABC(1200, 0.002)0.357436.503
NSGA-II(1200, 0.002)0.289237.021

800 MWMO-ABC(500, 0.002)0.148138.105
NSGA-II(500, 0.002)0.106237.929

850 MWMO-ABC(500, 0.002)0.123930.856
NSGA-II(500, 0.002)0.059029.499

Table 5 shows that the HV indicators of MO-ABC are higher than NSGA-II in all three load conditions. That is, the solution sets of MO-ABC all outperform NSGA-II (illustrated in Figures 79) in the comprehensive performance with comparable time consumption (shown in the last column of Table 5). This means MO-ABC is more applicable to boiler combustion optimization than NSGA-II.

Furthermore, we can present our findings in our current study whether the nondominated solutions by MO-ABC shown in Figures 79 are better than real values (i.e., experimental values). Among them, six recommended nondominated solutions with “+” (the principle of selecting these recommended nondominated solutions is to keep NOx as low as possible while maintaining the reheat steam temperature above its designed value) are shown in Tables 68 (No. 1–6) as well as the experiment values of NOx emission and reheat steam temperature (No. 7).

NOxThermal load (C)Thermal load (D)Thermal load (E)Thermal load (F)Running


NOxThermal load (B)Thermal load (C)Thermal load (D)Thermal load (E)Thermal load (F)Running


NOxThermal load (C)Thermal load (D)Thermal load (E)Thermal load (F)Running


In order to quantify the degree of reduction of NOx and increase of Tr relative to the experimental value, percentage improvement ( and ) are defined as follows:where and are the NOx values obtained by MO-ABC and experiment, respectively. and are Tr values obtained by MO-ABC and experiment, respectively. Results are shown in the last two columns of Tables 68.

By comparing results in Figures 79 and Tables 68, we can conclude that optimal NOx emission and reheat steam temperature achieved by MO-ABC is better than that explored directly from experimental ones (i.e., real data) based on the percentage improvement. This verified that our new viewpoint on combustion optimization is meaningful. To this end, joint optimization of NOx emission and reheat steam temperature can bring much more benefits to the combustion system with the points of synthetic optimization view.

5. Conclusions

This paper puts a new point on optimizing NOx emission and reheat steam temperature simultaneously using a joint optimization method, in which a set of rdCV-MO-NPLS strategies were proposed as predictors in three cases of load such as 650 MW, 800 MW, and 850 MW. Compared with another three well-known methods (k-fold-MO-NPLS, SVR, and ANN), this method can achieve higher predictive accuracies. Then, MO-ABC was applied to search the optimal set point of controllable process variables to reduce NOx emission and improve reheat steam temperature. Results showed that our joint optimization of boiler combustion with MO-ABC provided a set of tradeoff solutions and outperformed that obtained by NSGA-II and experimental data. This implies that our proposed method on boiler combustion optimization can guarantee higher economy as well as safety control of combustion systems.

Data Availability

The data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

The authors declare no conflicts of interest.


This work was supported by the National Natural Science Foundation of China (Grant no. 51676034) and the Fundamental Research Funds for the Central Universities (Grant no. 2242019K40002).

Supplementary Materials

Appendix A: experimental data are provided when the boiler worked stably, respectively, at 650 MW, 800 MW, and 850 MW. More details can be viewed in Tables A1–A3. Table A1: the steady-state experimental data of 650 MW. Table A2: the steady-state experimental data of 800 MW. Table A3: the steady-state experimental data of 850 MW. Appendix B: Table B1: lower and upper bounds of inputs of the rdCV-MO-NPLS model at 650 MW. Table B2: lower and upper bounds of inputs of the rdCV-MO-NPLS model at 800 MW. Table B3: lower and upper bounds of inputs of the rdCV-MO-NPLS model at 850 MW. (Supplementary Materials)


  1. Z. Wei, X. Li, L. Xu, and C. Tan, “Optimization of operating parameters for low NOx emission in high-temperature air combustion,” Energy & Fuels, vol. 26, no. 5, pp. 2821–2829, 2012. View at: Publisher Site | Google Scholar
  2. S. Belosevic, V. Beljanski, I. Tomanovic, N. Crnomarkovic, D. Tucakovic, and T. Zivanovic, “Numerical analysis of NOx control by combustion modifications in pulverized coal utility boiler,” Energy & Fuels, vol. 26, no. 1, pp. 425–442, 2012. View at: Publisher Site | Google Scholar
  3. S. Dal Secco, O. Juan, M. Louis-Louisy, J. Y. Lucas, P. Plion, and L. Porcheron, “Using a genetic algorithm and CFD to identify low NOx configurations in an industrial boiler,” Fuel, vol. 158, pp. 672–683, 2015. View at: Publisher Site | Google Scholar
  4. H. Zhou, K. Cen, and J. Fan, “Modeling and optimization of the nox emission characteristics of a tangentially fired boiler with artificial neural networks,” Energy, vol. 29, no. 1, pp. 167–183, 2004. View at: Publisher Site | Google Scholar
  5. P. Ilamathi, V. Selladurai, K. Balamurugan, and V. T. Sathyanathan, “ANN-GA approach for predictive modeling and optimization of NOx emission in a tangentially fired boiler,” Clean Technologies and Environmental Policy, vol. 15, no. 1, pp. 125–131, 2013. View at: Publisher Site | Google Scholar
  6. Z. Wei, X. Li, L. Xu, and Y. Cheng, “Comparative study of computational intelligence approaches for NOx reduction of coal-fired boiler,” Energy, vol. 55, no. 18, pp. 683–692, 2013. View at: Publisher Site | Google Scholar
  7. P.-H. Wang, L.-L. Li, Q. Chen, and Y.-H. Dong, “Research on applications of artificial intelligence to combustion optimization in a coal-fired boiler,” Proceedings of Chinese Society for Electrical Engineering, vol. 24, no. 4, p. 5, 2004. View at: Google Scholar
  8. X. Peng and P. Wang, “An improved multiobjective genetic algorithm in optimization and its application to high efficiency and low NOx emissions combustion,” in Proceedings of the 2009 Asia-pacific Power and Energy Engineering Conference, pp. 1–4, IEEE, Wuhan, China, March 2009. View at: Publisher Site | Google Scholar
  9. H. Zhou, K. Cen, and J. Fan, “Multi-objective optimization of the coal combustion performance with artificial neural networks and genetic algorithms,” International Journal of Energy Research, vol. 29, no. 6, pp. 499–510, 2005. View at: Publisher Site | Google Scholar
  10. Y. Zhang, Y. Ding, Z. Wu, L. Kong, and T. Chou, “Modeling and coordinative optimization of NOx emission and efficiency of utility boilers with neural network,” Korean Journal of Chemical Engineering, vol. 24, no. 6, pp. 1118–1123, 2007. View at: Publisher Site | Google Scholar
  11. Y. Zheng, X. Gao, and C. Sheng, “Impact of co-firing lean coal on NOx emission of a large-scale pulverized coal-fired utility boiler during partial load operation,” Korean Journal of Chemical Engineering, vol. 34, no. 4, pp. 1273–1280, 2017. View at: Publisher Site | Google Scholar
  12. F. M. Ham and I. Kostanic, “Partial least-squares: theoretical issues and engineering applications in signal processing,” Mathematical Problems in Engineering, vol. 2, no. 1, pp. 63–93, 1996. View at: Publisher Site | Google Scholar
  13. Z. Meng, S. Zhang, Y. Yang, and M. Liu, “Nonlinear partial least squares for consistency analysis of meteorological data,” Mathematical Problems in Engineering, vol. 2015, Article ID 143965, 8 pages, 2015. View at: Publisher Site | Google Scholar
  14. X. Dong and J. Wang, “Modeling and optimization for piercing efficiency and energy consumption based on mean value substaged KELM-PLS and GA method,” Mathematical Problems in Engineering, vol. 2014, Article ID 132654, 12 pages, 2014. View at: Publisher Site | Google Scholar
  15. P. Filzmoser, B. Liebmann, and K. Varmuza, “Repeated double cross validation,” Journal of Chemometrics, vol. 23, no. 4, pp. 160–171, 2009. View at: Publisher Site | Google Scholar
  16. S. Wold, N. Kettaneh-Wold, and B. Skagerberg, “Nonlinear PLS modeling,” Chemometrics and Intelligent Laboratory Systems, vol. 7, no. 1-2, pp. 53–65, 1989. View at: Publisher Site | Google Scholar
  17. R. Rosipal, “Nonlinear Partial Least Squares: An Overview,” in Chemoinformatics and Advanced Machine Learning Perspectives: Complex Computational Methods and Collaborative Techniques, pp. 169–189, IGI Global, Hershey, PA, USA, 2011. View at: Google Scholar
  18. G. Baffi, E. B. Martin, and A. J. Morris, “Non-linear projection to latent structures revisited: the quadratic PLS algorithm,” Computers & Chemical Engineering, vol. 23, no. 3, pp. 395–411, 1999. View at: Publisher Site | Google Scholar
  19. D. Karaboga, An Idea Based on Honey Bee Swarm for Numerical Optimization (Technical Report-tr06), Erciyes University, Engineering Faculty, Computer Engineering Department, Kayseri, Turkey, 2005, Technical Report.
  20. H. Wang, H. Liang, and L. Gao, “Using an improved artificial bee colony algorithm for parameter estimation of a dynamic grain flow model,” Mathematical Problems in Engineering, vol. 2018, Article ID 2132963, 11 pages, 2018. View at: Publisher Site | Google Scholar
  21. F. Wahid and D. H. Kim, “An efficient approach for energy consumption optimization and management in residential building using artificial bee colony and fuzzy logic,” Mathematical Problems in Engineering, vol. 2016, Article ID 9104735, 13 pages, 2016. View at: Publisher Site | Google Scholar
  22. G. Deng, H. Yang, and S. Zhang, “An enhanced discrete artificial bee colony algorithm to minimize the total flow time in permutation flow shop scheduling with limited buffers,” Mathematical Problems in Engineering, vol. 2016, Article ID 7373617, 11 pages, 2016. View at: Publisher Site | Google Scholar
  23. L. Sun, J. Hu, and H. Chen, “Artificial bee colony algorithm based on K-means clustering for multiobjective optimal power flow problem,” Mathematical Problems in Engineering, vol. 2015, Article ID 762853, 18 pages, 2015. View at: Publisher Site | Google Scholar
  24. D. Karaboga and B. Akay, “A comparative study of artificial bee colony algorithm,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 108–132, 2009. View at: Publisher Site | Google Scholar
  25. R. Akbari, R. Hedayatzadeh, K. Ziarati, and B. Hassanizadeh, “A multi-objective artificial bee colony algorithm,” Swarm and Evolutionary Computation, vol. 2, no. 1, pp. 39–52, 2012. View at: Publisher Site | Google Scholar
  26. D. Karaboga and B. Basturk, “A powerful and efficient algorithm for numerical function optimization: artificial bee colony (ABC) algorithm,” Journal of Global Optimization, vol. 39, no. 3, pp. 459–471, 2007. View at: Publisher Site | Google Scholar
  27. M. Mernik, S.-H. Liu, D. Karaboga, and M. Črepinšek, “On clarifying misconceptions when comparing variants of the artificial bee colony algorithm by offering a new implementation,” Information Sciences, vol. 291, pp. 115–127, 2015. View at: Google Scholar
  28. M. Li and X. Yao, “Quality evaluation of solution sets in multiobjective optimisation: a survey,” ACM Computing Surveys (CSUR), vol. 52, no. 2, pp. 1–38, 2019. View at: Publisher Site | Google Scholar

Copyright © 2019 Xian-hua Gao and Zhi-gang Su. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.