# A Complex Evaluation and Optimization Approach for Oxygen-Enriched Combustion Characteristics of Blended Fuels Based on Response Surface Methodology

Authors: Y. C. Liu; H. Zhang; Z. M. Luo; S. Qing; A. M. Zhang; S. P. Yang
Journal: Complexity
 (2020)

---

## Abstract

With the energy consumption increasing, the coal supply in China has been becoming tight, which has made it difficult for thermal power generation in Yunnan Province. Making full use of local inferior coal and biomass resources in Yunnan can remedy the lack of fuel in power plants. In this paper, an oxygen-rich atmosphere thermogravimetric experiment was performed for a blended sample of Xiaolongtan lignite, Yiliang tobacco rod, and Fuyuan bituminous coal. The combustion characteristics of the mixed fuel under several key operating parameters (i.e., mass ratios, oxygen concentration, and heating rates) were studied. The response surface methodology was used to determine the optimal blending ratio of the three fuels. The results show that the ignition and burnout temperature of coal decrease and the combustion time diminishes with the enrichment of oxygen. The optimal oxygen concentration in the practical application is around 30%. The activation energy and preexponential factor increase with the enlargement of oxygen concentration. Such complex evaluation and optimization approach ensure the stable operation of thermal power plant production.


---

## Body

## 1. Introduction

Combustion of solid fuel plays an essential role for the social development as it provides the energy needed in the practical operation [1–3]. With the development of the electric power industry, more and more power plants have been developed by using the blended fuels combustion technology in recent years. Hence, both the academic and industrial communities have paid increasing attention to the study of blended fuels. Meanwhile, the optimization of coal blending technology can provide additional benefits such as environmental protection and effective utilization of solid fuel resources. According to the experimental observation with regard to the blended combustion, different types of coal will intensely interact [4]. The ignition temperature of blended coal varies with different amounts of blended coal at low temperatures. Increasing the proportion of high volatile coal improves the ignition, stability, and burnout characteristics of blended coal. Chen et al. [5, 6] studied the mixed combustion characteristics of two kinds of pulverized coals with different ignition characteristics and found that the ignition characteristics of the blended coal are better than those of the single coal. The ignition temperature of blended coal tends to be higher as the proportion of pulverized coal with good ignition characteristics increases. Gao et al. [7] studied the combustion characteristics of blended lignite, bituminous coal, and anthracite by thermogravimetry. It was found that when the low activation energy bituminous coal or lignite is mixed with the anthracite with high activation energy, the ignition characteristic of blended coal is the best. However, when the activation energy of two kinds of coals is similar, the activation energy of blended coal tends to be the coal with relatively high activation energy. From the combustion experiments of blended coal via the simulating combustion conditions, similar conclusions can be drawn [8–10]. The characteristics of blended coal are not a simple linear superposition of single coal. The bigger the difference of coal quality is, the more noticeable this situation will be. Peratla et al. [11] studied the interaction of blended coal and found that the coal quality determines the combustion characteristics of blended coal, which depends on the proportion of combustible coal. Some researchers [12, 13] pointed out that the blended coal is the most critical factor in the coal combustion, and the burnout characteristics rely on the type of coal that burns slowly. At present, the studies for coal blending mainly focus on the blending of different coals while few efforts have been made to study the proportion of biomass smoke rods and coals. In addition, the performance of biomass is different from that of coal because the biomass mainly consists of high volatile substances, which have a low ignition point, rapid burning, and short burnout time. After the biomass is mixed with coal, the blended fuels have complex combustion characteristics. For the superposition of fuels, it needs to consider the coupling and synergy of the mixed fuels, which all depend on the proportion of different combustibles. In other words, the relationship between the coal blending and the main indexes of each coal composition is complexly nonlinear, which can be identified by nonlinear optimization models via response surface method (RSM) and fuzzy logic systems. Moreover, the blending ratio of blended fuels in most of the actual process is ambiguous so that the combustion performance of the blended fuels cannot be guaranteed.Motivated by the above requirements, this article will study the optimization problem of the fuels blending via the response surface method, aiming at finding an optimal mix scheme for different fuels with the best combustion performance achieved. The flammability index, comprehensive combustion characteristic index, burnout temperature, and economic costs are chosen to characterize the advantages and disadvantages of different mixing schemes with different parameters. After validated by extensive experiments, the proposed optimization technique is utilized to find the best ratio of different fuels. The experiments are divided into 10 groups, and each group has three kinds of fuels (i.e., lignite of Xiaolongtan, tobacco rod of Yiliang, and bituminous coal of Fuyuan) with different proportions. The advantages and disadvantages of using the response surface method are analyzed, and the different proportional reference assignment schemes are discussed. The analysis and experimental results will shed light on the blending of three kinds of fuels to achieve an excellent combustion efficiency.This paper is organized as follows: Section2 gives the problem formulation. Section 3 gives the blended coal studies, which includes raw materials and instruments illustration, thermogravimetric analysis, and chemical structure analysis. Section 4 gives the decision analysis of the response surface method, which includes the optimization of oxygen concentration and heating rate and optimizing the proportion of blended fuel. Section 5 gives the conclusions.


## 2. Problem Formulation

Combustion of blended fuels should consider its ignition characteristics, volatile precipitation characteristics, slagging, burn-up, and pollutant emissions. However, the blending ratio of blended fuels in most of the actual process is blind and arbitrary so that the combustion performance of the blended fuels cannot be guaranteed. This article aims at discussing the pros and cons of different proportions of blended fuels by using the response surface method, including its flammability index, comprehensive combustion characteristics index, burnout temperature, and economic cost, and finally getting the best ratio of three kinds of fuels blended coal combustion.Under different O2/CO2 atmospheres, thermogravimetric experiments were performed on blending fuels for three samples. The response surface method was used to analyze and optimize the pros and cons of different proportioning schemes to find the optimal ratio. First, the three kinds of fuels are divided into several groups according to different proportions. The different groups include one kind of fuel and a mixture of three kinds of fuel. The air atmosphere is set as 30% O2 atmosphere, 50% O2 atmosphere, and the heating rate of 5 K/min, 10 K/min, and 20 K/min, respectively.


## 3. Blended Fuel Combustion Studies

### 3.1. Raw Materials and Instruments

Xiaolongtan lignite (X), Yiliang tobacco rod (Y), and Fuyuan bituminous coal (F) are ground by ball mill to obtain the blended fuel with a diameter of 53∼75 μm. The atmosphere in the furnace is air and O2/CO2 mixture containing 30% O2, and the protective gas is high purity argon (i.e., 99.9%).The experimental instrument is STA 449 F3 Jupiter synchronous thermal analyzer (TGDSC, Germany, NETZSCH) and corundum crucible.


### 3.2. Thermogravimetric Analysis

Three sample ratios were selected randomly to carry out the experiments. The results showed that the ratio of mixed coal has no influence on the pyrolysis property, while the concentration of O2 affects this aspect. Taking the X : Y : F mass ratio of 1 : 1:1 and 2 : 2 : 1 in the experiment, the thermogravimetric curves of blended coal samples under different atmospheres with a heating rate of 20 K/min are shown in Figure 1. Higher oxygen concentration in the atmosphere gives rise to a more inclined thermogravimetric curve of blended coal and increases the maximum weight loss rate, indicating that the blended coal segregates in the oxygen-rich atmosphere. With the increase of oxygen concentration, the ignition temperature and burnout temperature of the blended coal decrease dramatically, which indicates that the combustion reaction of blended coal can be completed at a lower temperature and the combustion time can be reduced, which is consistent with the results of references [14, 15].Figure 1
Effect of oxygen concentration on pyrolysis of blended coal with different mass ratios.Figure2 shows the thermogravimetric curves of blended coal with different mass ratios at the heating rate of 20 K/min. It is noted that the curves have a similar trend. The curve is the most inclined to the high temperature zone with the X : Y : F mass ratio of 0 : 2 : 1 while the curve is the most inclined to the low temperature zone with the X : Y : F mass ratio of 2 : 1 : 1, indicating that the burning speed of blended coal is the fastest. Moreover, the better rate is conducive to its complete combustion [16, 17].Figure 2
TG curves of mixing coals with different mass ratios.


### 3.3. Chemical Structure Analysis

Fuel is mainly composed of the carbon skeleton, which is surrounded by various alkyl side chains and oxygen-containing functional group branched chains. The chemical properties of the carbon skeleton are stable at low temperature, while the properties of side chains and branched chains are more active. During the combustion, the groups of side chains and branched chains are easy to participate in the chemical reactions, mainly the oxidation of active functional groups and the precipitation of volatile compounds in coal. The content of aliphatic side chains and oxygen functional groups are higher than that of lignite, which is also an essential factor leading to a better combustion performance of bituminous coal than lignite. Therefore, with the increase of the bituminous coal ratio, the combustion performance of blended fuel is gradually improved [18].


## 3.1. Raw Materials and Instruments

Xiaolongtan lignite (X), Yiliang tobacco rod (Y), and Fuyuan bituminous coal (F) are ground by ball mill to obtain the blended fuel with a diameter of 53∼75 μm. The atmosphere in the furnace is air and O2/CO2 mixture containing 30% O2, and the protective gas is high purity argon (i.e., 99.9%).The experimental instrument is STA 449 F3 Jupiter synchronous thermal analyzer (TGDSC, Germany, NETZSCH) and corundum crucible.


## 3.2. Thermogravimetric Analysis

Three sample ratios were selected randomly to carry out the experiments. The results showed that the ratio of mixed coal has no influence on the pyrolysis property, while the concentration of O2 affects this aspect. Taking the X : Y : F mass ratio of 1 : 1:1 and 2 : 2 : 1 in the experiment, the thermogravimetric curves of blended coal samples under different atmospheres with a heating rate of 20 K/min are shown in Figure 1. Higher oxygen concentration in the atmosphere gives rise to a more inclined thermogravimetric curve of blended coal and increases the maximum weight loss rate, indicating that the blended coal segregates in the oxygen-rich atmosphere. With the increase of oxygen concentration, the ignition temperature and burnout temperature of the blended coal decrease dramatically, which indicates that the combustion reaction of blended coal can be completed at a lower temperature and the combustion time can be reduced, which is consistent with the results of references [14, 15].Figure 1
Effect of oxygen concentration on pyrolysis of blended coal with different mass ratios.Figure2 shows the thermogravimetric curves of blended coal with different mass ratios at the heating rate of 20 K/min. It is noted that the curves have a similar trend. The curve is the most inclined to the high temperature zone with the X : Y : F mass ratio of 0 : 2 : 1 while the curve is the most inclined to the low temperature zone with the X : Y : F mass ratio of 2 : 1 : 1, indicating that the burning speed of blended coal is the fastest. Moreover, the better rate is conducive to its complete combustion [16, 17].Figure 2
TG curves of mixing coals with different mass ratios.


## 3.3. Chemical Structure Analysis

Fuel is mainly composed of the carbon skeleton, which is surrounded by various alkyl side chains and oxygen-containing functional group branched chains. The chemical properties of the carbon skeleton are stable at low temperature, while the properties of side chains and branched chains are more active. During the combustion, the groups of side chains and branched chains are easy to participate in the chemical reactions, mainly the oxidation of active functional groups and the precipitation of volatile compounds in coal. The content of aliphatic side chains and oxygen functional groups are higher than that of lignite, which is also an essential factor leading to a better combustion performance of bituminous coal than lignite. Therefore, with the increase of the bituminous coal ratio, the combustion performance of blended fuel is gradually improved [18].


## 4. Decision Analysis of the Response Surface Method

### 4.1. Optimization of Oxygen Concentration and Heating Rate

Miscellaneous hybrid design of response surface optimization method in Design-Expert 8.0.6 software was used to design and optimize oxygen concentration and heating rate. The optimal mixing scheme was obtained by taking the combustion characteristic index and burnout temperature as the evaluation indexes. The corresponding response was filled with the calculated results based on experimental data. The design scheme and experimental evaluation indexes are shown in Table1.Table 1
Miscellaneous design and evaluating indicator.
Sequence no.Oxygen concentration,A (%, φ)Heating rate,B (K/min)Comprehensive combustion characteristic index,S (×10−8)Burnout temperature (°C)150.520.020.97516.3250.512.59.26472.9380.012.513.12419.8480.05.05.91406.2550.55.03.36422.9650.512.59.81483.7780.020.039.02427.821.020.013.84616.6950.512.58.87503.51021.05.01.92453.71121.012.57.33526.8


#### 4.1.1. Combustion Characteristic Index and Variable Relation Model

Tables2 and 3 are the variance and evaluation of the combustion characteristic indices fitted by various models, respectively. Table 4 is the confidence degree of the quadratic model. From Table 2, it is noted that the quadratic equation fits best among six models, while the cubic equation fits the worst. From Table 3, the quadratic equation has the best fitting result. Table 4 shows that the quadratic equation has a significant fitting effect and the heating rate has a significant effect on the comprehensive combustion.Table 2
The variance of comprehensive combustion characteristic index of different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value1816.1711816.17Linear model763.112381.5621.270.00042FI10011008.120.0261Quadratic equation83.31241.669.480.0183RecommendedCubic equation30.1725.0923.360.0041PoorResidual deflection3.6250.72Grand total2796.3813215.11Table 3
Comprehensive evaluation ofR2 of comprehensive combustion characteristic index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model4.210.79840.76020.5012541.182FI3.180.86170.85790.4631606.15Quadratic equation1.960.93660.95170.6982327.08RecommendedCubic equation0.770.99700.99140.6771392.12PoorTable 4
The confidence of the two equations for comprehensive combustion characteristic index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue8.971.000.836.7310.76Oxygen concentration,A5.311.000.773.377.321.12Heating rate,B9.951.000.918.2611.881.12AB4.831.001.022.267.271.12A21.471.001.03−1.284.311.03B23.941.001.151.087.061.03Figure3 shows the distribution of studentized residuals for the comprehensive combustion characteristic index fitting model. The dots represent comprehensive combustion characteristic indices. The red values stand for high probability while the blue values represent low probability. It is noted that the residual points are basically distributed near the straight line, indicating that the model has a good applicability. The role of the oxygen concentration and heating rate in the comprehensive combustion characteristic index is shown in Figure 4. As the oxygen concentration increases, the role of the heating rate decreases. As the heating rate increases, the role of the oxygen concentration decreases.Figure 3
Distribution of the residuals of the combustible index.Figure 4
Response surface of comprehensive combustion characteristic index and oxygen.


#### 4.1.2. Relation Model between Burnout Temperature and Variables

Tables5 and 6 present the variance comparisons of burnout temperatures fitted by various models. From the tables, the fitting effect of the quadratic equation is the best. Table 7 shows the confidence of the quadratic equation of the burnout temperature, which shows that the fitting effect of the quadratic model is remarkable, and the heating rate has the greatest influence on burnout temperature.Table 5
The variance of burnout temperature with different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value3.173 × 10613.173 × 106Linear model33677.22216838.6132.26<0.00012FI3987.6513987.6526.820.0004Quadratic equation838.422419.216.330.02271RecommendedCubic equation128.92264.461.080.37528PoorResidual deflection2.9650.59Grand total2786.2513214.33Table 6
Comprehensive evaluation ofR2 of burnout temperature.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model22.880.83490.81340.652512953.472FI11.920.94370.93910.91423388.57Quadratic equation8.020.98170.97990.95473446.29RecommendedCubic equation7.340.99170.98070.0783134337.29PoorTable 7
The confidence of the two equations for burnout temperature.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue492.7713.55483.48502.87Oxygen concentration,A−60.1813.46−64.35−49.381.00Heating rate,B50.8313.4641.0752.771.00AB−29.9814.26−42.66−23.531.00A2−13.7315.18−27.81−2.641.12B2−4.7215.18−15.866.8321.12Figure5 shows the distribution of residuals of the fitting model in Design-Expert 8.0.6. The ignition represents the burnout temperature. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the fitting effect of the model is satisfactory.Figure 5
Distribution of the residuals of burnout temperature.The control power of oxygen concentration and heating rate on burnout temperature is shown in Figure6. With the increase of oxygen concentration, the control power of the heating rate increases. With the increase of heating rate, the control power of oxygen concentration increases. Therefore, in order to get the lowest burnout temperature, the heating rate and oxygen concentration can be reduced simultaneously.Figure 6
Response surface of the relationship between burnout temperature and oxygen concentration and heating rate.


#### 4.1.3. Optimizing Process Control Parameters

To make the comprehensive combustion characteristic index the largest and the burnout temperature the smallest, the experiment under certain oxygen concentration and heating rate should be optimized. The range of response index is designed in Numerical Criteria. The analysis via software results of the optimum prediction scheme is shown in Table8. It is noted that when the oxygen concentration is 80%(φ) and the heating rate is 20 K/min, the expectation is the highest. In this case, the best results can be obtained by combining the combustion characteristic index and burnout temperature.Table 8
Comprehensive optimization scheme of various indexes of coal under different oxygen concentration and heating rate.
Sequence numberOxygen concentration (%,φ)Heating rate (K/min)Comprehensive combustion characteristic index,S (×10−8)Burnout temperature (°C)Degree of expectation180.0020.0035.9821438.8480.63274.8220.0032.6452457.6460.61370.2120.0033.0137470.3250.53421.005.003.59823469.0180.21Figures7 and 8 are the optimized contours and their three-dimensional response surface diagrams according to the relationship between the oxygen concentration and heating rate using the optimized scheme 1. Contours represent the specific values of the significant factors in different intervals of the target value, and their shapes reflect the interaction between the independent variables to control the effect of variables. Figure 7 shows the contours change from elliptic arc to monotonic shape with oxygen concentration when the expectation is 1. The role of heating rate increases first and then decreases when the degree of expectation is less than 1. Meanwhile, the role of heating rate increases gradually with the increase of oxygen concentration. The response surface graph of Figure 8 shows the change of expectation degree with oxygen concentration and heating rate.Figure 7
The optimized contour map of oxygen concentration and heating rate.Figure 8
Expected response surface map between desirability of oxygen concentration and heating rate using optimization scheme 1.


### 4.2. Optimizing the Proportion of Blended Fuel

Design-Expert 8.0.6 software was used to design the scheme and analyze the results of the thermogravimetric experiment of blended fuel. The four evaluation indices, that is, combustibility index, comprehensive combustion characteristic index, burnout temperature, and economic cost, were optimized and the optimal mixing scheme was obtained.For blended coal of Xiaolongtan lignite, Yiliang tobacco rod, and Fuyuan bituminous, a total of 27 groups of orthogonal experiments with 3 factors and 3 levels are needed. The Design-Expert 8.0.6 Box-Behnken is used to carry out the partial orthogonal experiments. The experimental points are the skeleton points of the experimental model according to the software and model analysis. As shown in Figure9, there are four repetitive points on each surface of the model skeleton. The designed repetitive experimental group also needs repetition. If the same set of experimental data is used for analysis, the results of the model are inaccurate.Figure 9
Test point distribution of three-factor BBD experimental design.According to the design scheme, the corresponding response column is filled according to the calculation results of the experimental data. The designed scheme and the experimental evaluation index are shown in Table9. It can be seen that the trend of the change of the combustibility index is generally the same as that of the comprehensive combustion characteristic index, and the combustibility index is taken for analysis.Table 9
Design and results of Box-Behnken orthogonal experiment.
Experiment no.Mass ratio ofX : Y : FCombustible index (×10−7)Comprehensive combustion characteristic index,S (×10−10)Burnout temperature (°C)Cost (¥/t)11 : 1 : 135.9313.2567370138321 : 1 : 134.6273.1052870038331 : 1 : 134.1533.2652769538341 : 1 : 133.6473.3248970138351 : 1 : 135.7433.2543270338362 : 1 : 042.5844.3939569129671 : 2 : 227.7092.4858969539880 : 1 : 011.1471.2316270728692 : 0 : 130.3282.98336679363102 : 1 : 250.6325.46327678402111 : 0 : 011.4751.18341668276121 : 2 : 037.8713.42526676265131 : 0 : 233.0133.68168709488140 : 2 : 127.8711.99386716378152 : 2 : 131.9872.68282718351160 : 1 : 230.1662.25216709497170 : 0 : 14.7390.41678742606


#### 4.2.1. Relation Model between Combustibility Index and Variables

Table10 shows the flammability index variance of blended coal of different proportions, which can be fitted by various models. From the table, the quadratic model has the best fitting effect, which is recommended for the system. The fitting effect of the cubic equation is the poorest. Table 11 compares the standard deviation, R2, and its correction value, predicted value, and squared sum of predicted residual of four fitting polynomial models. Apparently, the results show that the quadratic equation has the best fitting results.Table 10
The variance of combustibility index using different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionQuadratic sum15973.686115973.686Linear model1298.8973432.96576.37690.00592FI368.1963122.7322.43670.1213Quadratic equation356.9623118.98738.8960.0078RecommendedCubic equation96.981332.327137.2340.0002PoorResidual deflection0.92440.231Grand total18095.65171064.45Table 11
Comprehensive evaluation ofR2 of combustibility index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model7.9860.5990.5380.2191706.4862FI7.0380.7360.6670.0812003.612Quadratic equation3.4760.9460.8780.3021613.548RecommendedCubic equation0.5131.0000.998PoorTable12 shows the results of software analysis on the quadratic equation and the confidence of the influencing factors. It shows that the quadratic equation fitting effect is remarkable, and the interaction between rich source bituminous coal and rich source bituminous coal has the greatest influence on the flammability index.Table 12
The confidence of the two equations for combustibility index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue3.2510.552.264.58Xiaolongtan lignite,A−0.5110.44−1.320.631.00Yiliang tobacco stem,B−1.8610.44−2.58−0.591.00Fuyuan bituminous coal,C−2.3310.44−3.02−1.321.00AB0.6110.63−0.871.791.00AC2.5210.630.993.581.00BC2.3110.630.713.271.00A20.0210.60−1.371.541.01B2−0.4110.60−1.651.221.01C22.6710.601.243.761.01


#### 4.2.2. Optimizing the Proportion of Blended Fuel

Figure10 is the distribution of the residual of the flammability index fitting model. The point represents the flammability index. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the simulation effect of the model is very good.Figure 10
Distribution of the residuals of the combustible index.Figure11 shows the relationship between the ratio of the other two samples and the flammability index when the ratio of Fuyuan bituminous coal is 50%. When the flammability index is less than 20, the control power of Yiliang tobacco rod decreases with the increase of Xiaolongtan lignite. When the flammability index is 30, the content of Xiaolongtan lignite decreases first and then increases with the increase of Yiliang smoke stem.Figure 11
Response surface of the relationship between combustible index of mixing fuels and mass ratios of Xiaolongtan lignite and Yiliang tobacco rod.Figure12 shows the relationship between the proportion of the other two kinds of coals and the combustibility index when the proportion of Yiliang tobacco rod is 50%. The combustibility index increases with the increase of Xiaolongtan lignite and the combustibility index increases with the increase of rich source bituminous coal. Based on the fitting of experiment and model, the experiment is optimized by using Design-Expert 8.0.6 software, and the combustibility index, burnout temperature, and economic cost are optimized. In the best case, we can find the optimal mixing ratio scheme. Design Criteria of Design-Expert 8.0.6 software are also used to design the range of values of each factor and set the flammability index and combustion characteristic index to be the largest and the burnout temperature and investment cost to be the smallest.Figure 12
Response surface of the relationship between combustible index of mixing coals and mass ratios of Xiaolongtan lignite and Fuyuan bituminous coal.Through the analysis, the optimal prediction scheme is shown in Table13. From the optimization results, scheme 1 is the optimal mixing ratio, which is also suggested by the software. Xiaolongtan lignite: Yiliang tobacco rod: Fuyuan bituminous coal mass ratio is 100 : 55.12 : 1. Under this condition, the model predicts that the combustibility index of blended fuel is 43.3574 × 10−7, the combustion characteristic index is 4.27 × 10−10, the burnout temperature is 678.112°C, and the cost is 251.688 ¥/t, demonstrating that the expectation is the highest.Table 13
Comprehensive optimization scheme of each index.
Sequence no.Mass ratio (%)Combustible index × 10−7Combustion performance index × 10−10Burnout temperature (°C)Investment cost (¥/t)Degree of expectationXiaolongtan ligniteYiliang tobacco stemFuyuan bituminous coal1100.0055.121.0043.35744.19875678.112251.6880.8852100.0055.641.0043.48664.14176678.332251.5340.885399.5855.121.0043.33054.32126678.862251.3780.884499.3756.371.0043.49274.38726678.375250.0230.885598.2256.811.0043.39934.32462678.562251.6380.884697.6256.131.0043.27844.14324679.012251.8680.883796.8857.661.0043.26424.14389679.326250.6580.881895.7967.861.0043.98734.19868685.465241.6780.8729100.0035.871.3239.98553.96782670.483266.7140.84310100.0046.0236.8343.01174.15865692.821322.9330.80411100.0036.33100.0046.98655.19546682.757438.2680.79312100.0035.05100.0048.09635.29232682.962439.3930.7931399.2736.78100.0047.67525.18383682.878439.8870.79414100.0033.8987.27945.88924.98712690.009420.6530.78815100.0041.0271.4344.56714.53527692.656386.8480.782


## 4.1. Optimization of Oxygen Concentration and Heating Rate

Miscellaneous hybrid design of response surface optimization method in Design-Expert 8.0.6 software was used to design and optimize oxygen concentration and heating rate. The optimal mixing scheme was obtained by taking the combustion characteristic index and burnout temperature as the evaluation indexes. The corresponding response was filled with the calculated results based on experimental data. The design scheme and experimental evaluation indexes are shown in Table1.Table 1
Miscellaneous design and evaluating indicator.
Sequence no.Oxygen concentration,A (%, φ)Heating rate,B (K/min)Comprehensive combustion characteristic index,S (×10−8)Burnout temperature (°C)150.520.020.97516.3250.512.59.26472.9380.012.513.12419.8480.05.05.91406.2550.55.03.36422.9650.512.59.81483.7780.020.039.02427.821.020.013.84616.6950.512.58.87503.51021.05.01.92453.71121.012.57.33526.8


### 4.1.1. Combustion Characteristic Index and Variable Relation Model

Tables2 and 3 are the variance and evaluation of the combustion characteristic indices fitted by various models, respectively. Table 4 is the confidence degree of the quadratic model. From Table 2, it is noted that the quadratic equation fits best among six models, while the cubic equation fits the worst. From Table 3, the quadratic equation has the best fitting result. Table 4 shows that the quadratic equation has a significant fitting effect and the heating rate has a significant effect on the comprehensive combustion.Table 2
The variance of comprehensive combustion characteristic index of different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value1816.1711816.17Linear model763.112381.5621.270.00042FI10011008.120.0261Quadratic equation83.31241.669.480.0183RecommendedCubic equation30.1725.0923.360.0041PoorResidual deflection3.6250.72Grand total2796.3813215.11Table 3
Comprehensive evaluation ofR2 of comprehensive combustion characteristic index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model4.210.79840.76020.5012541.182FI3.180.86170.85790.4631606.15Quadratic equation1.960.93660.95170.6982327.08RecommendedCubic equation0.770.99700.99140.6771392.12PoorTable 4
The confidence of the two equations for comprehensive combustion characteristic index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue8.971.000.836.7310.76Oxygen concentration,A5.311.000.773.377.321.12Heating rate,B9.951.000.918.2611.881.12AB4.831.001.022.267.271.12A21.471.001.03−1.284.311.03B23.941.001.151.087.061.03Figure3 shows the distribution of studentized residuals for the comprehensive combustion characteristic index fitting model. The dots represent comprehensive combustion characteristic indices. The red values stand for high probability while the blue values represent low probability. It is noted that the residual points are basically distributed near the straight line, indicating that the model has a good applicability. The role of the oxygen concentration and heating rate in the comprehensive combustion characteristic index is shown in Figure 4. As the oxygen concentration increases, the role of the heating rate decreases. As the heating rate increases, the role of the oxygen concentration decreases.Figure 3
Distribution of the residuals of the combustible index.Figure 4
Response surface of comprehensive combustion characteristic index and oxygen.


### 4.1.2. Relation Model between Burnout Temperature and Variables

Tables5 and 6 present the variance comparisons of burnout temperatures fitted by various models. From the tables, the fitting effect of the quadratic equation is the best. Table 7 shows the confidence of the quadratic equation of the burnout temperature, which shows that the fitting effect of the quadratic model is remarkable, and the heating rate has the greatest influence on burnout temperature.Table 5
The variance of burnout temperature with different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value3.173 × 10613.173 × 106Linear model33677.22216838.6132.26<0.00012FI3987.6513987.6526.820.0004Quadratic equation838.422419.216.330.02271RecommendedCubic equation128.92264.461.080.37528PoorResidual deflection2.9650.59Grand total2786.2513214.33Table 6
Comprehensive evaluation ofR2 of burnout temperature.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model22.880.83490.81340.652512953.472FI11.920.94370.93910.91423388.57Quadratic equation8.020.98170.97990.95473446.29RecommendedCubic equation7.340.99170.98070.0783134337.29PoorTable 7
The confidence of the two equations for burnout temperature.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue492.7713.55483.48502.87Oxygen concentration,A−60.1813.46−64.35−49.381.00Heating rate,B50.8313.4641.0752.771.00AB−29.9814.26−42.66−23.531.00A2−13.7315.18−27.81−2.641.12B2−4.7215.18−15.866.8321.12Figure5 shows the distribution of residuals of the fitting model in Design-Expert 8.0.6. The ignition represents the burnout temperature. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the fitting effect of the model is satisfactory.Figure 5
Distribution of the residuals of burnout temperature.The control power of oxygen concentration and heating rate on burnout temperature is shown in Figure6. With the increase of oxygen concentration, the control power of the heating rate increases. With the increase of heating rate, the control power of oxygen concentration increases. Therefore, in order to get the lowest burnout temperature, the heating rate and oxygen concentration can be reduced simultaneously.Figure 6
Response surface of the relationship between burnout temperature and oxygen concentration and heating rate.


### 4.1.3. Optimizing Process Control Parameters

To make the comprehensive combustion characteristic index the largest and the burnout temperature the smallest, the experiment under certain oxygen concentration and heating rate should be optimized. The range of response index is designed in Numerical Criteria. The analysis via software results of the optimum prediction scheme is shown in Table8. It is noted that when the oxygen concentration is 80%(φ) and the heating rate is 20 K/min, the expectation is the highest. In this case, the best results can be obtained by combining the combustion characteristic index and burnout temperature.Table 8
Comprehensive optimization scheme of various indexes of coal under different oxygen concentration and heating rate.
Sequence numberOxygen concentration (%,φ)Heating rate (K/min)Comprehensive combustion characteristic index,S (×10−8)Burnout temperature (°C)Degree of expectation180.0020.0035.9821438.8480.63274.8220.0032.6452457.6460.61370.2120.0033.0137470.3250.53421.005.003.59823469.0180.21Figures7 and 8 are the optimized contours and their three-dimensional response surface diagrams according to the relationship between the oxygen concentration and heating rate using the optimized scheme 1. Contours represent the specific values of the significant factors in different intervals of the target value, and their shapes reflect the interaction between the independent variables to control the effect of variables. Figure 7 shows the contours change from elliptic arc to monotonic shape with oxygen concentration when the expectation is 1. The role of heating rate increases first and then decreases when the degree of expectation is less than 1. Meanwhile, the role of heating rate increases gradually with the increase of oxygen concentration. The response surface graph of Figure 8 shows the change of expectation degree with oxygen concentration and heating rate.Figure 7
The optimized contour map of oxygen concentration and heating rate.Figure 8
Expected response surface map between desirability of oxygen concentration and heating rate using optimization scheme 1.


## 4.1.1. Combustion Characteristic Index and Variable Relation Model

Tables2 and 3 are the variance and evaluation of the combustion characteristic indices fitted by various models, respectively. Table 4 is the confidence degree of the quadratic model. From Table 2, it is noted that the quadratic equation fits best among six models, while the cubic equation fits the worst. From Table 3, the quadratic equation has the best fitting result. Table 4 shows that the quadratic equation has a significant fitting effect and the heating rate has a significant effect on the comprehensive combustion.Table 2
The variance of comprehensive combustion characteristic index of different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value1816.1711816.17Linear model763.112381.5621.270.00042FI10011008.120.0261Quadratic equation83.31241.669.480.0183RecommendedCubic equation30.1725.0923.360.0041PoorResidual deflection3.6250.72Grand total2796.3813215.11Table 3
Comprehensive evaluation ofR2 of comprehensive combustion characteristic index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model4.210.79840.76020.5012541.182FI3.180.86170.85790.4631606.15Quadratic equation1.960.93660.95170.6982327.08RecommendedCubic equation0.770.99700.99140.6771392.12PoorTable 4
The confidence of the two equations for comprehensive combustion characteristic index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue8.971.000.836.7310.76Oxygen concentration,A5.311.000.773.377.321.12Heating rate,B9.951.000.918.2611.881.12AB4.831.001.022.267.271.12A21.471.001.03−1.284.311.03B23.941.001.151.087.061.03Figure3 shows the distribution of studentized residuals for the comprehensive combustion characteristic index fitting model. The dots represent comprehensive combustion characteristic indices. The red values stand for high probability while the blue values represent low probability. It is noted that the residual points are basically distributed near the straight line, indicating that the model has a good applicability. The role of the oxygen concentration and heating rate in the comprehensive combustion characteristic index is shown in Figure 4. As the oxygen concentration increases, the role of the heating rate decreases. As the heating rate increases, the role of the oxygen concentration decreases.Figure 3
Distribution of the residuals of the combustible index.Figure 4
Response surface of comprehensive combustion characteristic index and oxygen.


## 4.1.2. Relation Model between Burnout Temperature and Variables

Tables5 and 6 present the variance comparisons of burnout temperatures fitted by various models. From the tables, the fitting effect of the quadratic equation is the best. Table 7 shows the confidence of the quadratic equation of the burnout temperature, which shows that the fitting effect of the quadratic model is remarkable, and the heating rate has the greatest influence on burnout temperature.Table 5
The variance of burnout temperature with different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionAverage value3.173 × 10613.173 × 106Linear model33677.22216838.6132.26<0.00012FI3987.6513987.6526.820.0004Quadratic equation838.422419.216.330.02271RecommendedCubic equation128.92264.461.080.37528PoorResidual deflection2.9650.59Grand total2786.2513214.33Table 6
Comprehensive evaluation ofR2 of burnout temperature.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model22.880.83490.81340.652512953.472FI11.920.94370.93910.91423388.57Quadratic equation8.020.98170.97990.95473446.29RecommendedCubic equation7.340.99170.98070.0783134337.29PoorTable 7
The confidence of the two equations for burnout temperature.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue492.7713.55483.48502.87Oxygen concentration,A−60.1813.46−64.35−49.381.00Heating rate,B50.8313.4641.0752.771.00AB−29.9814.26−42.66−23.531.00A2−13.7315.18−27.81−2.641.12B2−4.7215.18−15.866.8321.12Figure5 shows the distribution of residuals of the fitting model in Design-Expert 8.0.6. The ignition represents the burnout temperature. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the fitting effect of the model is satisfactory.Figure 5
Distribution of the residuals of burnout temperature.The control power of oxygen concentration and heating rate on burnout temperature is shown in Figure6. With the increase of oxygen concentration, the control power of the heating rate increases. With the increase of heating rate, the control power of oxygen concentration increases. Therefore, in order to get the lowest burnout temperature, the heating rate and oxygen concentration can be reduced simultaneously.Figure 6
Response surface of the relationship between burnout temperature and oxygen concentration and heating rate.


## 4.1.3. Optimizing Process Control Parameters

To make the comprehensive combustion characteristic index the largest and the burnout temperature the smallest, the experiment under certain oxygen concentration and heating rate should be optimized. The range of response index is designed in Numerical Criteria. The analysis via software results of the optimum prediction scheme is shown in Table8. It is noted that when the oxygen concentration is 80%(φ) and the heating rate is 20 K/min, the expectation is the highest. In this case, the best results can be obtained by combining the combustion characteristic index and burnout temperature.Table 8
Comprehensive optimization scheme of various indexes of coal under different oxygen concentration and heating rate.
Sequence numberOxygen concentration (%,φ)Heating rate (K/min)Comprehensive combustion characteristic index,S (×10−8)Burnout temperature (°C)Degree of expectation180.0020.0035.9821438.8480.63274.8220.0032.6452457.6460.61370.2120.0033.0137470.3250.53421.005.003.59823469.0180.21Figures7 and 8 are the optimized contours and their three-dimensional response surface diagrams according to the relationship between the oxygen concentration and heating rate using the optimized scheme 1. Contours represent the specific values of the significant factors in different intervals of the target value, and their shapes reflect the interaction between the independent variables to control the effect of variables. Figure 7 shows the contours change from elliptic arc to monotonic shape with oxygen concentration when the expectation is 1. The role of heating rate increases first and then decreases when the degree of expectation is less than 1. Meanwhile, the role of heating rate increases gradually with the increase of oxygen concentration. The response surface graph of Figure 8 shows the change of expectation degree with oxygen concentration and heating rate.Figure 7
The optimized contour map of oxygen concentration and heating rate.Figure 8
Expected response surface map between desirability of oxygen concentration and heating rate using optimization scheme 1.


## 4.2. Optimizing the Proportion of Blended Fuel

Design-Expert 8.0.6 software was used to design the scheme and analyze the results of the thermogravimetric experiment of blended fuel. The four evaluation indices, that is, combustibility index, comprehensive combustion characteristic index, burnout temperature, and economic cost, were optimized and the optimal mixing scheme was obtained.For blended coal of Xiaolongtan lignite, Yiliang tobacco rod, and Fuyuan bituminous, a total of 27 groups of orthogonal experiments with 3 factors and 3 levels are needed. The Design-Expert 8.0.6 Box-Behnken is used to carry out the partial orthogonal experiments. The experimental points are the skeleton points of the experimental model according to the software and model analysis. As shown in Figure9, there are four repetitive points on each surface of the model skeleton. The designed repetitive experimental group also needs repetition. If the same set of experimental data is used for analysis, the results of the model are inaccurate.Figure 9
Test point distribution of three-factor BBD experimental design.According to the design scheme, the corresponding response column is filled according to the calculation results of the experimental data. The designed scheme and the experimental evaluation index are shown in Table9. It can be seen that the trend of the change of the combustibility index is generally the same as that of the comprehensive combustion characteristic index, and the combustibility index is taken for analysis.Table 9
Design and results of Box-Behnken orthogonal experiment.
Experiment no.Mass ratio ofX : Y : FCombustible index (×10−7)Comprehensive combustion characteristic index,S (×10−10)Burnout temperature (°C)Cost (¥/t)11 : 1 : 135.9313.2567370138321 : 1 : 134.6273.1052870038331 : 1 : 134.1533.2652769538341 : 1 : 133.6473.3248970138351 : 1 : 135.7433.2543270338362 : 1 : 042.5844.3939569129671 : 2 : 227.7092.4858969539880 : 1 : 011.1471.2316270728692 : 0 : 130.3282.98336679363102 : 1 : 250.6325.46327678402111 : 0 : 011.4751.18341668276121 : 2 : 037.8713.42526676265131 : 0 : 233.0133.68168709488140 : 2 : 127.8711.99386716378152 : 2 : 131.9872.68282718351160 : 1 : 230.1662.25216709497170 : 0 : 14.7390.41678742606


### 4.2.1. Relation Model between Combustibility Index and Variables

Table10 shows the flammability index variance of blended coal of different proportions, which can be fitted by various models. From the table, the quadratic model has the best fitting effect, which is recommended for the system. The fitting effect of the cubic equation is the poorest. Table 11 compares the standard deviation, R2, and its correction value, predicted value, and squared sum of predicted residual of four fitting polynomial models. Apparently, the results show that the quadratic equation has the best fitting results.Table 10
The variance of combustibility index using different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionQuadratic sum15973.686115973.686Linear model1298.8973432.96576.37690.00592FI368.1963122.7322.43670.1213Quadratic equation356.9623118.98738.8960.0078RecommendedCubic equation96.981332.327137.2340.0002PoorResidual deflection0.92440.231Grand total18095.65171064.45Table 11
Comprehensive evaluation ofR2 of combustibility index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model7.9860.5990.5380.2191706.4862FI7.0380.7360.6670.0812003.612Quadratic equation3.4760.9460.8780.3021613.548RecommendedCubic equation0.5131.0000.998PoorTable12 shows the results of software analysis on the quadratic equation and the confidence of the influencing factors. It shows that the quadratic equation fitting effect is remarkable, and the interaction between rich source bituminous coal and rich source bituminous coal has the greatest influence on the flammability index.Table 12
The confidence of the two equations for combustibility index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue3.2510.552.264.58Xiaolongtan lignite,A−0.5110.44−1.320.631.00Yiliang tobacco stem,B−1.8610.44−2.58−0.591.00Fuyuan bituminous coal,C−2.3310.44−3.02−1.321.00AB0.6110.63−0.871.791.00AC2.5210.630.993.581.00BC2.3110.630.713.271.00A20.0210.60−1.371.541.01B2−0.4110.60−1.651.221.01C22.6710.601.243.761.01


### 4.2.2. Optimizing the Proportion of Blended Fuel

Figure10 is the distribution of the residual of the flammability index fitting model. The point represents the flammability index. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the simulation effect of the model is very good.Figure 10
Distribution of the residuals of the combustible index.Figure11 shows the relationship between the ratio of the other two samples and the flammability index when the ratio of Fuyuan bituminous coal is 50%. When the flammability index is less than 20, the control power of Yiliang tobacco rod decreases with the increase of Xiaolongtan lignite. When the flammability index is 30, the content of Xiaolongtan lignite decreases first and then increases with the increase of Yiliang smoke stem.Figure 11
Response surface of the relationship between combustible index of mixing fuels and mass ratios of Xiaolongtan lignite and Yiliang tobacco rod.Figure12 shows the relationship between the proportion of the other two kinds of coals and the combustibility index when the proportion of Yiliang tobacco rod is 50%. The combustibility index increases with the increase of Xiaolongtan lignite and the combustibility index increases with the increase of rich source bituminous coal. Based on the fitting of experiment and model, the experiment is optimized by using Design-Expert 8.0.6 software, and the combustibility index, burnout temperature, and economic cost are optimized. In the best case, we can find the optimal mixing ratio scheme. Design Criteria of Design-Expert 8.0.6 software are also used to design the range of values of each factor and set the flammability index and combustion characteristic index to be the largest and the burnout temperature and investment cost to be the smallest.Figure 12
Response surface of the relationship between combustible index of mixing coals and mass ratios of Xiaolongtan lignite and Fuyuan bituminous coal.Through the analysis, the optimal prediction scheme is shown in Table13. From the optimization results, scheme 1 is the optimal mixing ratio, which is also suggested by the software. Xiaolongtan lignite: Yiliang tobacco rod: Fuyuan bituminous coal mass ratio is 100 : 55.12 : 1. Under this condition, the model predicts that the combustibility index of blended fuel is 43.3574 × 10−7, the combustion characteristic index is 4.27 × 10−10, the burnout temperature is 678.112°C, and the cost is 251.688 ¥/t, demonstrating that the expectation is the highest.Table 13
Comprehensive optimization scheme of each index.
Sequence no.Mass ratio (%)Combustible index × 10−7Combustion performance index × 10−10Burnout temperature (°C)Investment cost (¥/t)Degree of expectationXiaolongtan ligniteYiliang tobacco stemFuyuan bituminous coal1100.0055.121.0043.35744.19875678.112251.6880.8852100.0055.641.0043.48664.14176678.332251.5340.885399.5855.121.0043.33054.32126678.862251.3780.884499.3756.371.0043.49274.38726678.375250.0230.885598.2256.811.0043.39934.32462678.562251.6380.884697.6256.131.0043.27844.14324679.012251.8680.883796.8857.661.0043.26424.14389679.326250.6580.881895.7967.861.0043.98734.19868685.465241.6780.8729100.0035.871.3239.98553.96782670.483266.7140.84310100.0046.0236.8343.01174.15865692.821322.9330.80411100.0036.33100.0046.98655.19546682.757438.2680.79312100.0035.05100.0048.09635.29232682.962439.3930.7931399.2736.78100.0047.67525.18383682.878439.8870.79414100.0033.8987.27945.88924.98712690.009420.6530.78815100.0041.0271.4344.56714.53527692.656386.8480.782


## 4.2.1. Relation Model between Combustibility Index and Variables

Table10 shows the flammability index variance of blended coal of different proportions, which can be fitted by various models. From the table, the quadratic model has the best fitting effect, which is recommended for the system. The fitting effect of the cubic equation is the poorest. Table 11 compares the standard deviation, R2, and its correction value, predicted value, and squared sum of predicted residual of four fitting polynomial models. Apparently, the results show that the quadratic equation has the best fitting results.Table 10
The variance of combustibility index using different models.
Source of varianceQuadratic sumDegree of freedomMean squareF-valueProbability >FScheme selectionQuadratic sum15973.686115973.686Linear model1298.8973432.96576.37690.00592FI368.1963122.7322.43670.1213Quadratic equation356.9623118.98738.8960.0078RecommendedCubic equation96.981332.327137.2340.0002PoorResidual deflection0.92440.231Grand total18095.65171064.45Table 11
Comprehensive evaluation ofR2 of combustibility index.
ModelStandard deviationR2R2 corrected valueR2 predicted valuePredicted residual sum of squaresScheme selectionLinear model7.9860.5990.5380.2191706.4862FI7.0380.7360.6670.0812003.612Quadratic equation3.4760.9460.8780.3021613.548RecommendedCubic equation0.5131.0000.998PoorTable12 shows the results of software analysis on the quadratic equation and the confidence of the influencing factors. It shows that the quadratic equation fitting effect is remarkable, and the interaction between rich source bituminous coal and rich source bituminous coal has the greatest influence on the flammability index.Table 12
The confidence of the two equations for combustibility index.
FactorParameter estimationDegree of freedomStandard deviation95% confidence interval lower control limit95% confidence interval upper control limitSignificant factorValue3.2510.552.264.58Xiaolongtan lignite,A−0.5110.44−1.320.631.00Yiliang tobacco stem,B−1.8610.44−2.58−0.591.00Fuyuan bituminous coal,C−2.3310.44−3.02−1.321.00AB0.6110.63−0.871.791.00AC2.5210.630.993.581.00BC2.3110.630.713.271.00A20.0210.60−1.371.541.01B2−0.4110.60−1.651.221.01C22.6710.601.243.761.01


## 4.2.2. Optimizing the Proportion of Blended Fuel

Figure10 is the distribution of the residual of the flammability index fitting model. The point represents the flammability index. The red values stand for high probability while the blue values represent low probability. It is noted that the residuals are basically distributed near the straight line, which shows that the simulation effect of the model is very good.Figure 10
Distribution of the residuals of the combustible index.Figure11 shows the relationship between the ratio of the other two samples and the flammability index when the ratio of Fuyuan bituminous coal is 50%. When the flammability index is less than 20, the control power of Yiliang tobacco rod decreases with the increase of Xiaolongtan lignite. When the flammability index is 30, the content of Xiaolongtan lignite decreases first and then increases with the increase of Yiliang smoke stem.Figure 11
Response surface of the relationship between combustible index of mixing fuels and mass ratios of Xiaolongtan lignite and Yiliang tobacco rod.Figure12 shows the relationship between the proportion of the other two kinds of coals and the combustibility index when the proportion of Yiliang tobacco rod is 50%. The combustibility index increases with the increase of Xiaolongtan lignite and the combustibility index increases with the increase of rich source bituminous coal. Based on the fitting of experiment and model, the experiment is optimized by using Design-Expert 8.0.6 software, and the combustibility index, burnout temperature, and economic cost are optimized. In the best case, we can find the optimal mixing ratio scheme. Design Criteria of Design-Expert 8.0.6 software are also used to design the range of values of each factor and set the flammability index and combustion characteristic index to be the largest and the burnout temperature and investment cost to be the smallest.Figure 12
Response surface of the relationship between combustible index of mixing coals and mass ratios of Xiaolongtan lignite and Fuyuan bituminous coal.Through the analysis, the optimal prediction scheme is shown in Table13. From the optimization results, scheme 1 is the optimal mixing ratio, which is also suggested by the software. Xiaolongtan lignite: Yiliang tobacco rod: Fuyuan bituminous coal mass ratio is 100 : 55.12 : 1. Under this condition, the model predicts that the combustibility index of blended fuel is 43.3574 × 10−7, the combustion characteristic index is 4.27 × 10−10, the burnout temperature is 678.112°C, and the cost is 251.688 ¥/t, demonstrating that the expectation is the highest.Table 13
Comprehensive optimization scheme of each index.
Sequence no.Mass ratio (%)Combustible index × 10−7Combustion performance index × 10−10Burnout temperature (°C)Investment cost (¥/t)Degree of expectationXiaolongtan ligniteYiliang tobacco stemFuyuan bituminous coal1100.0055.121.0043.35744.19875678.112251.6880.8852100.0055.641.0043.48664.14176678.332251.5340.885399.5855.121.0043.33054.32126678.862251.3780.884499.3756.371.0043.49274.38726678.375250.0230.885598.2256.811.0043.39934.32462678.562251.6380.884697.6256.131.0043.27844.14324679.012251.8680.883796.8857.661.0043.26424.14389679.326250.6580.881895.7967.861.0043.98734.19868685.465241.6780.8729100.0035.871.3239.98553.96782670.483266.7140.84310100.0046.0236.8343.01174.15865692.821322.9330.80411100.0036.33100.0046.98655.19546682.757438.2680.79312100.0035.05100.0048.09635.29232682.962439.3930.7931399.2736.78100.0047.67525.18383682.878439.8870.79414100.0033.8987.27945.88924.98712690.009420.6530.78815100.0041.0271.4344.56714.53527692.656386.8480.782


## 5. Conclusion

In this work, the response surface methodology was used to optimize the oxygen concentration, heating rate, and pulverized coal ratio. The effects of coal ratio on the combustibility index, combustion characteristic index, burnout temperature, and the investment cost of blended coal were analyzed. The following conclusions can be drawn:(1)
In the oxygen-enriched atmosphere, the overall precipitation of mixed coal becomes better. With the increase of oxygen concentration, the combustion rate increases, and the influence of oxygen concentration is greater than that of mixing proportion. Meanwhile, the ignition temperature and burnout temperature of pulverized coal decrease significantly, indicating that the pulverized coal can be burned at the lower temperature and combustion time can be shortened under the oxygen-enriched condition.(2)
When the expectation is 1, the role of heating rate increases first and then decreases with the increase of oxygen concentration. When the expectation is less than 1, the control force of heating rate increases gradually with the increase of oxygen concentration.(3)
Setting the flammability index and combustion characteristic index to be the largest, the burnout temperature and investment cost will be the smallest. The optimal prediction scheme 1 is the case that the mass ratio of Xiaolongtan lignite: Yiliang tobacco rod: Fuyuan bituminous coal is 100 : 55.12 : 1. Under this condition, the model predicts that the combustibility index of blended fuel is 43.3574 × 10−7, the combustion characteristic index is 4.27 × 10−10, the burnout temperature is 678.112°C, and the cost is 251.688 ¥/t.



---

*Source: 2487310-2020-07-15.xml*