Flame Failures and Recovery in Industrial Furnaces : A Neural Network Steady-State Model for the Firing Rate Setpoint Rearrangement

Burner failures are common abnormal conditions associated with industrial fired heaters. Preventing from economic loss and major equipment damages can be attained by compensating the lost heat due to burners’ failures, which can be possible by defining appropriate setpoints to rearrange the firing rates for healthy burners. In this study, artificial neural network models were developed for estimating the appropriate setpoints for the combustion control system to recover an industrial fired-heater furnace from abnormal conditions. For this purpose, based on an accurate high-order mathematical model, constrained nonlinear optimization problems were solved using the genetic algorithm. For different failure scenarios, the best possible excess firing rates for healthy burners to recover the furnace from abnormal conditions were obtained and data were recorded for training and testing stages. (e performances of the developed neural steady-state models were evaluated through simulation experiments. (e obtained results indicated the feasibility of the proposed technique to deal with the failures in the combustion system.


Introduction
Fired heaters are critical units in refineries and petrochemical industries, used for direct heating of a fluid stream in hydrocarbon processing.Hot gases from combustion are directly supplied to the multiple stream tubes by means of convective and radiative heat transfer mechanisms. is considerably increases the temperature of the crude process, which can be fed into distillation columns to be separated into a variety of components [1,2].In order to maintain the quality of products, the stream temperature should be controlled.Regulating the firing rate is normally the main tool for achieving target temperature [3].
Nonuniform temperature distribution inside furnaces due to the burners' failures and combustion malfunctions may cause major quality defects in distillation products [4].
e effects of burners' failures may appear in different forms such as flame impingement, flame instability, hot spots, and so on [5].Ignoring this situation could increase the danger of major equipment damages.Reducing the faulty burner's fuel flow is the main correcting action [6].In case of failure of a single burner or losing a part of burner firing capacity, the furnace can still continue to function; however, compensating the lost heat to achieve thermal equilibrium would be the main issue.In such circumstances, which are called abnormal conditions, corrective actions should be immediately taken in order to avoid coke formation inside the tubes and consequent possible damages [6,7].
e strategy to recover the system from abnormal conditions can be characterized by defining appropriate setpoints for the controllers at the lower layer.Integrating expert systems into the monitoring and supervisory control systems has been regarded as a solution to improving the capability of industrial control systems in reacting to unexpected events [8][9][10].Such systems should be able to acquire data from instruments, detect and diagnose the faults, and optimize the process conditions by adjusting new setpoints for control loops to move away from abnormal conditions [11,12].Designing an effective expert system usually requires adequate information from the plant variables, real data from critical components conditions, and the knowledge of expert operators about certain abnormal situations [13].However, the required information from real system performances during the failures of burners is barely in hand.Furthermore, applying artificial failures to inservice fired-heaters to acquire the necessary data is practically impossible [14,15].In this case, mathematical models could be employed to investigate the effects of many different failures and firing rate scenarios on heat distribution inside the furnace and optimize the process conditions to cope with abnormal situations [16][17][18].
Many different model-based optimization techniques have been reported for improving the operation of industrial processes in the literature [19][20][21].First-principle models (FPMs) are the common type of steady-state models used for plant monitoring, process optimization, and open-loop decision supports [22].Due to model uncertainties and disturbances, online steady-state data are required to identify the models.Employing steady-state models could considerably reduce the computational efforts, which is so valuable in real-time optimization (RTO) [23].Nonlinear mathematical models have been also used for optimizing the operation of industrial processes, which are generally used as the core of nonlinear model predictive control (NMPC) [24,25].
In recent years, data-based modeling approaches such as artificial neural networks have been widely employed for optimization of industrial processes [18,26,27].
e capability of characterizing the behaviors of complicated systems without prior knowledge, low computational costs, and good interpolating performances is the main advantage of neural network models [28].Different structures of artificial neural networks have been used for optimizing the operation of industrial processes; however, in general, two distinct categories can be identified in their applications: static neural network models are more common in real-time optimization (RTO) applications, whereas steady-state models are required for optimizing the operational conditions at the upper control layers; (2) dynamic neural models are used as the core of model predictive control at the lower control levels [29][30][31].
In this study, a model-based optimization approach is proposed to find the optimal operational conditions for an industrially fired-heater furnace to be recovered from possible burner failures.A nonlinear mathematical model developed by Chaibakhsh et al. was extended and used for optimization process [32].Different failure scenarios in the combustion system were applied to the model, and their effects on the tubes and furnace outputs were investigated.
en, constrained nonlinear optimization problems were solved to plan the required excess firing rate for healthy burners to recover from abnormal conditions.e obtained results were employed for training the neural network (NN) estimators.e developed models were used as the core of an abnormal management system (AMS), which was responsible for generating the reference setpoints for the burners' controllers according to the failure situations.e obtained results from the optimization process and the assessment of the effectiveness of the proposed approach were presented for different failure situations.

System Description
A cabinet-type, single-firebox furnace shown in Figure 1 is considered for investigation in this study.Five floormounted burners are arranged in a row to supply the required heat for convection and radiation sections.e crude process is fed from the top of the furnace and first passes through the convection sections 2 and 1, where its temperature increases from 504.15 K to about 561.65 K, and the fluid often remains as a one-phase liquid.In the radiation section, tubes are exposed to high heat fluxes.e crude process begins with boiling and vaporizing and the temperature increases to about 637 K at the end of this section.
At each wall, the process fluid is divided into two counterflow streams, which are heading to the middle side of the furnace.Each stream also comprises 24 rows of tubes.
e required firing rate is about 0.06842 kg/s for each burner at the nominal furnace load of 11.26 kg/s.In Table 1, the nominal parameters of the furnace are presented [32,33].
An accurate mathematical model was developed by Chaibakhsh et al., which was employed for the optimization process [32].In order to consider the effects of flame height on the heat absorbed by each tube, the firebox of the fired heater has been divided into five thermal zones, which are shown in Figure 1. e view factor for each subsection of tubes was calculated with respect to their location and the average flame height.

Optimization Problem: Objective Function and Genetic Algorithm
In this section, different burners' firing capacity shortage/failure scenarios are applied to the developed model of the furnace.en, an optimization approach based on the genetic algorithm (GA) was applied to the model in order to find the appropriate firing rate for each burner.e maximum firing rate capacity for each burner, the overall fuel consumption, the tubes' wall temperature, and the furnace outlet temperature are considered as constrains.e procedure of model-based optimization is shown in Figure 2.

Optimization Problem.
e average temperature of two streams (on both sides) is the most important variable that should remain around 637 K (Table 1); this is necessary to achieve the desired product quality and should not be exceeded significantly.us, by defining error as the difference between the average outlet temperature of the crude process and the desired temperature (637 K), the optimum excess firing rate for each burner could be calculated to reduce the outlet temperature errors, according to the constraints on variables.

2
International Journal of Chemical Engineering e objective function for the optimization process can be defined as follows: where T oil and T tube denote the state variables oil and tube surface temperature, respectively.In addition, _ m fuel and _ m oil are fuel and crude process flow rates fed to the furnace, respectively.is objective function is used for each zone of each tube, and the optimization process is performed for all abnormal conditions, which may occur during systems operation.
e considered functions for the optimization problem definition are as follows: e function of the tube's temperature in the radiation section is presented in (2), whereas the following equation shows the relation between crude process's temperature and the heat absorbed from the tubes: Cov. sec.
( International Journal of Chemical Engineering (3) In the convection section, (4) is used for estimating the tube's temperature, and (5) is the crude process's temperature function:

International Journal of Chemical Engineering
In addition, a limitation is considered on the excess firing rate for each burner as follows: e maximum tube design temperature is restricted as e inputs and outputs of all passes are located in Zone 1; accordingly, the process output temperature T k would be calculated as In case of failure in a burner, heat loss should be compensated by increasing the firing rate in the other burners.It should be noted that the furnace's thermal zones would change due to any variation in flame height.In order to determine the heat absorbed by each section of the tubes, the view factor for each subsection of tubes should be calculated with respect to their location and the average flame height.
e geometric description for the tube and flame and the effects of flame's height on the view factors are presented in Appendix.Increasing the temperature of tube surfaces is the most significant effects of overfiring, which should be limited due to its critical effects on the safety of the furnace.e maximum tube design temperature is about 900 K; however, increasing the surface temperature of the tube could increase the danger of coke formation and internal tube fouling.
erefore, the maximum target surface temperature is considered to be 800 K as follows: e appropriate firing rate for each burner can be measured by solving the nonlinear optimization problem in the form of (2) to (5) subject to inequality constrains of ( 6) and ( 7) using the genetic algorithm.In order to keep solutions in the feasible region, it is required to convert a constrained optimization problem to an unconstrained optimization [34].Employing penalty functions are a common approach to dealing with constrains in solving nonlinear optimization problems by GA.It is applied as follows: e penalty function P S k is simplified for the stream k as where T maxS k is the maximum surface temperature of tubes (24 passes for each stream pass) and is given as e penalty functions could prevent the tubes' surface temperature from reaching the maximum values.ey have a zero value if the surface temperature is less than the maximum value.By excessing the considered limit, the penalty function would have a nonzero (positive) value, which could be magnified by considering a large value for penalty constant α.
In order to compensate for the heat loss due to burner failures, the sum of excess mass flow rates in healthy burners should be equal to the missed flow rates in faulty burners.However, by increasing the burners' firing rate, the heat transfer performance will increase.
erefore, this value could be a little (about 3 percent) less than the required amount.e maximum value of excess firing rate is also limited to 30 percent.It is applied by e penalty function for limiting the excess fuel flow rate can be defined through the following equation: e penalty function, P fuel , is considered for economical operation, and the optimization process is performed based on reducing fuel consumption.However, in case of the failure/firing capacity shortage of two burners, it is not economically justifiable.In these cases, the burner's firing rate may be close to their maximum capacity, meaning that the total fuel flow rate might be more than its amount in normal operating conditions.
As a result, the fitness function is defined as where P(.) is a penalty function and α and β are penalty constants.is is a constrained optimization problem, which was transformed to an unconstrained one using an additive penalty function P(.) [34].

Genetic Algorithm (GA).
Genetic algorithm (GA) methodologies have been investigated as efficient approaches for obtaining solutions very close to optimum points [35,36].
GA is a nongradient-based stochastic technique, which is capable of exploring and exploiting the search space for the optimal solution.e GA operators including crossover and mutation allow the searching algorithm to escape from possible local minima trapping points.In addition, it is very robust to dynamic changes in the environment.is advantage has made it an appropriate tool for the optimization of systems working under different operating conditions with different behaviors [37].Different applications of GA were reported for designing decision-making systems based on nonlinear optimization problems [38].
e optimization processes were performed by using MATLAB ® Optimization Toolbox ( e MathWorks, Inc., Natick, MA, USA, ver.2014).An appropriate selection of GA parameters could improve the performance of algorithm and reduce the runtime considerably [39].Crossover is a recombination operator that could generate better offspring through generations by selecting randomly a pair of two individuals, choosing a cross site along the string length and finally swapping between the two strings following the cross site.After creation of new individuals via crossover, mutation is applied usually with a low probability to introduce random changes into the population [37].Generally, the crossover rate of about 80 to 90% should be chosen.However, some results show that for some problems the crossover rate of about 60% is the best.e mutation rate should be very low, and it is reported that the best rates is about 0.5 to 1%. e proper values of crossover and mutation rates can be chosen by performing some experiments to find the best solution.It is possible to estimate the optimal population sizes based on some analytic methods on particular problems [37].However, a default population size of about 20 to 30 is adequate to find the optimal solution; however, sometimes it has been reported that the size of 50 to 100 could be the best.Elite selection is the strategy to choose individuals with a best fitness value to pass to the next generation.e selection could be performed randomly according to the limited elites.In Table 2, the parameters for optimization by the GA are presented.

Neural Network Model
A multilayer perceptron (MLP) structure was considered for the neural network (NN) model.MLP is a feedforward neural network (FFNN), which can be used for characterizing complex relationships between variables and predicting future values in a process.MLP has a set of layers (including input, hidden, and output layers) and computational units (activation function), which are related to the layers' node by weight coefficients.Generally, one or more hidden layers can be considered in MLP; however, a single hidden layer might be adequate to reach an acceptable degree of accuracy in many cases [40,41].
e number of neurons, as a dominant parameter in the model performances, was considered with respect to the number of input/output variables.Azoff recommended that the number of neurons at the hidden layer (N h ) for a single hidden layer NN be chosen as follows: where N u is the number of network inputs [42].In order to adjust the parameters of models, the second-order derivativebased Levenberg-Marquardt (LM) was employed.e LM is a common algorithm for training the parameters of moderate-size MLPs.e mean square of error (MSE) was used as the objective function for the training process, which should be minimized as follows: where y * i and y i are the target and actual outputs for the ith pattern, respectively, and p is the total number of training patterns.e error vector E is calculated as follows: For the set of the considered inputs, x P , the output neurons, y k , are calculated through the following equation: where f yk and f hj are transfer functions related to output neurons, k, and hidden neurons, j, respectively; also, w kj is the weights interface between the kth output neurons and the jth hidden neurons, and w ji is the jth hidden neuron and the ith input neuron.

Failure Scenarios: Optimizations and Setpoints
In this section, 47 different constrained optimization problems were solved in order to obtain the appropriate firing rate for the burners in different failure scenarios.e It should be noted that the results are achieved through time-consuming optimization processes using a high-order nonlinear model; therefore, only a limited number of failure scenarios could be investigated.
e model employed for optimization processes was a model with 1064 states (ODE equations), where the runtime for each optimization process using a PC with Intel ® Core ™ i7 4510U @2.60 GHz, RAM 8.00 GB, and a operating system Windows ® 8.1 (x64) was about 24 hours.

Model Responses and Optimization Results
. Exposure to various disturbances, unplanned variations in the process flow rate, and changes in fuel components or pressure drop could affect the operating conditions and thermal performances of preheat furnaces [5,43].Useful information on the load pattern can be obtained from simulation experiments using the developed model.
is can be used in designing the expert furnace supervisory control and burner management systems in order to achieve the desired outlet process temperature at different operating conditions.
Figure 3(a) shows the crude process outlet temperature versus the fuel flow rate as the process flow changes.e required firing rate for burners to achieve desired process temperature of 637 K, at different process loads is shown.is value is about 0.06842 kg/s for each burner at the nominal furnace load of 11.26 kg/s.e changes in the process flow rate have considerable effects on the outlet temperature at a constant firing rate.For this reason, the process flow rate is considered as the main manipulated variable for regulating the outlet temperature of each path, and a control valve is supplied at the inlet of the individual path to the radiation section.
In Figure 3(b), the temperature differences of two stream passes are presented.e obtained results indicate that, by increasing the fuel flow rate, the outlet temperature difference between passes 1 and 2 would reach its maximum values as the total fuel flow rate exceeds up to 10 percent over the nominal firing capacity.By increasing the firing rate, the temperature difference starts to reduce again.With respect to the position of each pass in the furnace, the changes in the firing rate could affect the flame height, and consequently, the view factor between the flame and tube passes, which could cause a significant difference between the outlet temperatures of two stream passes.Reducing the temperature difference could improve the energy balance and result in much stable conditions inside the furnace.Some refineries try to increase their productivity by pushing the plant to its overdesign limits, which can be achieved by excessive firing.is limit for fired heaters is often 20 percent in normal operating conditions [6].In this case, the difference in temperature of stream passes would be small; however, excessive firing capacity may not be adequate to compensate for the effects of burner failures.
us, the proposed method could be applicable, if the furnace is run in its nominal design conditions.
Figure 3(c) presents the tube surface temperature versus the fuel flow rate at different furnace loads.A decrease in the process load causes the heat absorbed by the crude process to decrease and consequently the tube surface temperature to increase significantly, which is not desirable.e tube temperature is a critical variable that should not be tolerated from its acceptable range.e tube surface temperature does not normally reach its limit until the process fluid flows inside the tubes (Figure 3(c)).Continuous flame impingement on the tubes is the main potential problem that may occur and lead to coke formation and early furnace failure.e risk of tubes' surface localized overheating in places where no tube-skin thermocouple is installed is especially high.Generally, flame impingement can be discovered only by visual observation and should be avoided by any means possible.Localized overheating more than 38 °C above the designed temperature of tubes would be the cause for their failure in a matter of hours [7].In this case, the operator could make the flame lean away from the tubes by pulling the flame down [6].However, this could significantly deteriorate the performance of the furnace, which should be compensated by excess firing of the remaining healthy burners.
e main reason for choosing the genetic algorithm was its capability to search for best results in a global searching space.Despite it being a little slow for the high-order system to capture the optimum points, it would ensure the obtained results are global optimum solutions.By performing constrained optimization based on the GA, the best possible firing rate for each burner was obtained.Table 3 presents the required excess firing rate for healthy burners and the outlet temperature deviations as one burner fails (20 cases).Process flow enters the furnace at the front side and leaves it from the same side.Hence, the effects of failures in burners 1 or 5 and burners 2 or 3 would be considerably different in the outlet temperature.e maximum temperature deviations can be observed as the burners 1 and 5 fail; however, the failure of burner 1, the nearest burner to the front wall, can be recognized as the worst case.In these cases, excess firing rate reaches its maximum capacity of burners, and any further attempts to reduce the temperature deviations increase the risk of tubes overheating.
is causes a small reduction in the total fuel flow rate.
When the burners that are closer to the center be faulty, the temperature deviations were considerably reduced.e obtained results showed that the furnace could certainly be International Journal of Chemical Engineering recovered from abnormal conditions.However, a small increase in the fuel flow rate could be observed, making it economically less desirable.Despite the effectiveness of the proposed approach in dealing with faulty burners, due to economic reasons and also a high thermal stress, the furnace could not run for a long time in such conditions.e required firing rates for healthy burners to recover the heat balance inside the furnace as two faulty burners operate at   that the heat loss can be compensated by healthy burners as two faulty burners run at 75 percent of their firing capacity.However, this causes a significant increase in the total fuel flow International Journal of Chemical Engineering rate.When the faulty burners are located at the front side of the furnace the total fuel flow rate increases by about 14 percent.is means that a large amount of energy is wasted in longterm operations of the furnace.Conditions labeled as "Case 41" to "Case 47" are associated with cases that one faulty burner runs at 75 percent while the second one only has a half firing capacity.As presented in Table 4, regarding the faulty burners located at the two ends of the furnace, despite the considerable excess fuel flow rate, still a large deviation in outlet temperature can be observed.Such large outlet temperature changes may deteriorate the quality of distillation products.
e outlet temperature deviations reach their maximum values as the firing capacity of two burners reduced by 50 percent (Table 4).In these cases, considerable decreases in the total fuel flow rate can be observed, and it is clear that it is impossible to recover the furnace from abnormal conditions.For more than two faulty burners (e.g., three burners with 75 percent of firing capacity), it is impractical to recover the furnace to continue its operation at nominal load.In such cases, the furnace should inevitably be shut down or run under lower load conditions.In Table 5, the required fuel flow rates associated with the furnace load at different conditions are presented.As can be seen, there are no considerable changes in the fuel/load ratio. is allows dealing with abnormal situations by reducing the process load.
It should be noted that running fired-heater furnaces at low load conditions is not economically affordable, due to huge energy waste in the upstream and downstream units [7].In addition, by reducing the crude process flow rate, the fouling rate would increase inside the tubes [44].

Steady-State Model: Training and Validation.
e obtained data from the optimization process were employed to develop the NN model as the core of the abnormal management system (AMS).
e structure of the proposed system is shown in Figure 4, in which a bank of multi-inputsingle-output (MISO) models are employed as the reference estimators.For each burner, an individual MLP model with five inputs and one output was considered.e firing capacity of burners, which is defined as the ratio of actual fuel flow rates to its nominal values, was employed as the models'  ).e corrective positions for control valves were also considered as the models' output (y ∈ [0, 1.3]).e number of neurons in the hidden layer N h , based on (17), would be equal to 11. e data were collected for 47 data points, which are divided into training (70%), testing (15%), and validation (15%) subsets.
In order to assess the performances of the developed models, the coefficient of determination (R 2 ) and mean square error (MSE) are employed, as shown below: where y is the average of y over the n data, and y * i and y (i) p are the ith target and predicted responses, respectively.International Journal of Chemical Engineering e performances of the developed models are presented in Table 6.
e developed neural networks were employed to rearrange the setpoints for the burners' control system based on the detected failures.e developed NN model could estimate the reference setpoints for each burner by interpolating the untrained inputs in new situations.

Simulation Experiments.
In order to evaluate the performances of the proposed control system, simulation experiments were carried out in MATLAB Simulink environments.A schematic of the furnaces control systems is presented in Figure 5.With the occurrence of an event, for example, flame impingement, the corrective action by the operator is to reduce the flame height.In this case, the controller could automatically set the appropriate desired flow rate for other healthy burners to recover the furnace from abnormal conditions.By measuring the actual fuel flow rates, the corrective actions can be estimated by a bank of neural network models in order to compensate for the heat balance inside the furnace.
e operator is able to shut off the burners by a trip command or can set new setpoints for flow controllers from the control panel.Manual closing of the valves or controlling the fuel flow rate for burners in the field can be also detected by the control system.In case of abnormal situations or burner failures, the control system would be able to response appropriately to the new conditions.In Figure 6(a), the responses of the system, as burner #3 loses 80% of its nominal heat capacity, are presented.
is example shows that the AMS is able to detect the fault immediately after its occurrence and modify the burner's valve position, as shown in Figure 6(b).e obtained results indicate the feasibility and effectiveness of the developed controller in dealing with abnormal conditions.

Conclusion
In this study, steady-state neural network models were developed for estimating the fuel flow rate setpoints for the combustion control system to recover an industrially firedheater furnace from abnormal conditions.Constrained nonlinear optimization problems were solved using the genetic algorithm to find appropriate firing rates for healthy burners to compensate the heat losses due to failures, where the recorded data are used for training the NN models.e performances of the proposed system were evaluated at different burners' failure situations.e obtained results indicate the applicability of a model-based optimization approach to find the required firing rate for recovering the furnace from abnormal conditions. is provides conceptual basis for designing real-time abnormal management systems (AMS) for fired-heater furnaces.
e steady-state neural models can be used in real-time optimization in order to improve the operational performance of control systems.By considering the economic indices, the developed neural network models can be also employed as the core of economic model predictive controllers (EMPC).Attaining higher control performances and reducing energy wasting can be expected by combining these techniques.In addition, Δ i s is calculated with respect to the dimensions of the flame and tube as follows: e overall view factor of the flame for each tube subsection can be obtained by integrating F 12 into the corresponding region [24].e length of L 1 � L 2 and are equal to the flame width, L. It should be noted that the boundary of integration would change with respect to the position of each zones.We have obtained the following: (A.4) View factors were obtained by numerical evaluation of integral equations in MATLAB.Due to symmetric geometry of the firebox, view factors were calculated only for one side, which would considerably reduce computational efforts.

Figure 3 :
Figure 3: (a) Crude process temperature versus fuel flow rate.(b) e temperature difference of passes versus fuel flow rate.(c) Tube surface temperature versus fuel flow rate.

Figure 4 :
Figure 4: Structure of the neural network model.

Figure 5 :
Figure 5: A schematic of control system configuration for a fired-heater furnace unit.

Table 1 :
e nominal parameters of the furnace.
Figure 2: Procedure of the optimization problem for obtaining the appropriate firing rates.

Table 2 :
e optimization parameters for GA. the furnace employed for optimization is a continuous state model.e firing rates are discrete variables that are selected by the GA, and then, the equations are solved in order to obtain the temperatures of tubes and outlets of all passes.In the considered furnace, each burner has its individual control valve, which allows a considerable excessive firing rate up to 30 percent of the nominal firing rates for each burner.In this case, the heat loss due to a single burner failure or firing shortage in two or three burners may be compensated.For investigating the possibility of recovering the furnace from abnormal conditions, many different cases are considered.e failure scenarios applied to the burners are discrete variables.e simulation experiments and optimization are preformed while the available firing capacities are at 75, 50, and 25 percent of the nominal firing rates and also in the case of complete burner failures.e firing capacity of burners can be estimated by the measured mass flow rate of the fuel fed to each burner.It should be noted that no control actions should be applied during optimization stages.

Table 5 :
Fuel flow rate changes versus furnace load for economic operation evaluation.

Table 6 :
e performances of the developed neural network estimator.