Multiobjective Optimization of a Benfield HiPure Gas Sweetening Unit

We show how a multiobjective bare-bones particle swarm optimization can be used for a process parameter tuning and performance enhancement of a natural gas sweetening unit. This has been made through maximization of hydrocarbon recovery and minimization of the total energy of the process as the two objectives of the optimization. A trade-off exists between these two objectives as illustrated by the Pareto front. This algorithm has been applied to a sweetening unit that uses the Benfield HiPure process. Detailed models of the natural gas unit are developed in ProMax process simulator and integrated to the multi-objective optimization developed in visual basic environment (VBA). In this study, the solvent circulation rates, stripper pressure and reboiler duties are considered as the decision variables while hydrogen sulfide and carbon dioxide concentrations in the sweetened gas are considered as process constraints. The upper and lower bounds of the decision variables are obtained through a parametric sensitivity analysis of the models. The Pareto sets show a significant improvement in hydrocarbon recovery and a decent reduction in the heat consumption of the process.


Introduction
As global energy demand rises, natural gas now plays an important strategic role in world energy supply.It is the cleanest and most hydrogen rich of all the hydrocarbon energy sources and it has high energy efficiencies for power energy.Natural gas resources exploited and discovered are plentiful; however, they contain complex contaminants such as CO 2 , H 2 S, Mercaptans, and other sulfur compounds.Excessive amounts of these contaminants in natural gas streams will lead to low gas heating value and/or cause serious environmental hazards to the consumers.In LNG plants, large amounts of carbon dioxide and sulfur compounds may affect the quality of LNG products or pose serious operational problems in the cryogenic columns [1,2].Therefore, one of the major purposes of the sweetening units is to purify raw natural gas to meet both sale gas and liquefaction specifications [3].
Unpredicted changes in reservoir conditions, combined with the tough market competition, have forced the natural gas sweetening business to adopt sophisticated optimization techniques to improve the purity of their products, increase production capacity, and minimize total energy requirement for units' operations.
Gas sweetening processes use complex facilities whose design and operation basically depend on many parameters including gas composition, flow rate, circulation rates, absorber temperatures and stripper pressures.Accurate modeling of such processes involving multicriteria decision making has inspired process engineers to acquire the best practices for process automation in the way to realize competitive returns on investment.Despite the rapid development in process automation, no systematic methodology can be employed to continuously guarantee optimality of process operation.However, routine update of performance can always be achieved.
Problems associated with the gas treating processes such as large energy consumption and limited gas recovery have led to the increasing focus on mathematical modeling and simulation of gas sweetening processes, often performed to understand effects of design and operation on process enhancement.
In most cases, optimization of such complex processes is achieved through parametric sensitivity analysis [4,5].Such tedious procedures may be eliminated through the use of systematic approaches which guarantee global optimality of the process.
Previously, flow sheet optimization was achieved through calculation of inaccurate gradient derivatives of the process model generated inside modular process simulators [6].The development of high speed computers has led to increasing focus on designing and optimizing chemical process using accurate and complex stochastic algorithms.
Further, the decreased computational time has also contributed to routinely solving real-world problems involving large, realistic nonlinear models often encountered in process engineering.Researchers such as Aroonwilas et al. [7] and Park and Kang [8] simulated CO 2 removal with aqueous amine solutions such as DEA, MDEA, and DGA in a gas sweetening process.Rahimpour and Kashkooli [9,10] developed a mathematical model for the absorption of CO 2 by DEA-promoted potash solution and used it to examine the performance of a split flow absorber.More work on the modeling of the individual potassium carbonate and amine processes using computer simulations has been reported in the specialized literature [11][12][13][14][15].
Recently, engineering optimization has been heavily inclined towards nature-inspired and complex computational techniques such as evolutionary algorithms to achieve the necessary design and operation standards [16,17].These complex and accurate optimization techniques have also been greatly motivated by the market availability of modern highperformance computers.
Optimization problems are often computationally expensive especially those involving multicriteria decision making processes.In practice, real-world engineering problems involve making decisions over multiple objectives such as maximizing production capacity, minimizing energy, maximizing reliability, or any other objectives which could be conflicting.For solving such problems accurately, researchers call for multiobjective optimization techniques which are nongradient based.The application of evolutionary optimization methods to solve complex engineering multi-objective optimization problems has been published in a number of articles [18][19][20][21].
In this paper, we show how a classical multi-objective bare-bone particle swarm optimization (MOBBPSO) can be effectively utilized in process parameter tuning and performance analysis for a natural gas sweetening unit through hydrocarbon recovery maximization and energy minimization.The two objectives are combined and solved as a single objective problem through the weighted sum method described by Marler and Arora [22].
Once optimization basis is determined, it is necessary to determine process variables to be optimized.Process variables are the variables that affect the values of the objective functions.These can be divided into decision and dependent variables.Decision variables are determined while dependent variables are influenced by process constraints.In this study, concentrations of H 2 S and CO 2 in the purified gas were considered as the dependent variables.A parametric sensitivity study is carried out on the process model to determine the upper and lower bounds of the decision variables.The bounds are necessary to enable the optimizer to create a search space which obeys the necessary rules of thumb.
The optimization code developed in Visual Basic for Application (VBA) environment and then integrated with ProMax process simulator to tune the solvent circulation rates, reboiler duties, and regeneration pressures of the gas sweetening units, so that energy consumption and hydrocarbon recovery are optimized.
This paper is divided into six sections.Section 2 introduces multi-objective optimization through the weighted sum approach and also describes a brief theory on barebones particle swarm optimization.Section 3 describes the Benfield HiPure process and the use of ProMax to develop its model.This section also contains a comparison of simulation results with plant data obtained from ADGAS.Section 4 describes the parametric sensitivity analysis which explains how the constraints in Section 5 are obtained.The optimization results are presented and discussed in Section 6. Conclusions on this study including possible future work are discussed in Section 7 of this article.
The most common approach to multi-objective optimization is the weighted sum method (5) which transforms the multiple objective problem into a single objective engineering problem [22,23].Objectives are allocated weights and added to form a single objective problem (maximization objectives take a negative sign).Consider Multi-objective optimization in process plants has been recently reported by Wu et al. [21], Bernier et al. [24], Li et al. [25], Montazer-Rahmati and Binaee [26], Sun and Lou [27], and Shadiya et al. [28].Genetic algorithm [19] has always been the most common global optimization algorithm employed in optimizing process plants.Recently particle swarm optimization has gained more interest in obtaining global optimum of nonconvex problems [29].Poli et al. [30] discussed an overview on the development and application of particle swarm optimization to real-world problems, and Reyes-Sierra and Coella [31] presented a survey on use of multi-objective particle swarm optimizers in solving multiple objectives problems.Here, one of the improved versions of particle swarm optimization (PSO) commonly known as bare-bone particle swarm optimization (BBPSO) will be utilized to optimize the natural gas process.PSO has undergone several changes since its introduction in 1995 [32].One of these changes was the use of the Gaussian normal distribution of mean and standard deviation to update particle position.In BBPSO, the exploration search is facilitated by the standard deviation which tends to zero as the search space progresses.More information on the development and application of BBPSO is discussed in the work published by Zhang et al. [33].

Case Study
3.1.Benfield HiPure Process of Natural Gas Sweetening.The Benfield HiPure process was first described in 1974 by Benson and Parish [34].This process consists of the Benfield process and an amine unit downstream to scrub the remaining traces of acid gases.The integrated schematic of the Benfield HiPure process used in ADGAS plant at Das Island is shown in Figure 1.
The hot potassium carbonate absorption system comprises a split flow absorber and a regenerator with no side draws.The carbonate absorber and regenerator are both vertical packed bed columns.The treated gas from the carbonate absorber is fed directly into the amine absorber.The DEA amine system comprises absorber and stripper columns which are vertical and made of a stack-structured packing.The rich solution from the absorber is pumped to the DEA regenerator that has no condenser.The overhead gas (from DEA regenerator) is fed to the middle of the carbonate regenerator which does have a condenser.Liquid from the carbonate regenerator condenser is fed to the top of the DEA regenerator as reflux.
Sweet gas exiting the DEA absorber will undergo further processing before it is sent for LNG processing.The stripped acid gases from both the carbonate and DEA regenerators proceed to a sulfur recovery unit (SRU), where they are processed to produce molten liquid sulfur.Table 1 gives the typical operating data of the ADGAS plant.  2 gives a comparison of the simulation results and the operating data.As can be seen, there is a close match between the simulation results with plant data.In the next section, a sensitivity analysis study will be described in order to show the effect of various operating conditions on the performance of the process and to motivate the optimization study that will be discussed in a later section of this paper.

Sensitivity Analysis
Sensitivity analysis gives the best approach to investigate the effect of process parameters on performance of a chemical plant.Here, a scenario tool was developed in ProMax, and the effect of change in decision variables is monitored with respect to the required gas specification (process constraint) while compromising the respective rules of thumb.
Figure 2 shows the effect of change in the carbonate circulation rate with composition of H 2 S and CO 2 in the purified gas.As shown in the figure, decreasing the carbonate circulation rate beyond 1400 m 3 /hr will rapidly cause an increase in H 2 S and CO 2 composition in the purified gas, hence violating the constraints of not more than 1 ppmv H 2 S and 50 ppmv CO 2 .On the other hand, a carbonate circulation rate of 1700 m 3 /hr or more will produce sweet gas with almost the same gas specification.Therefore, 1400 and 1700 m 3 /hr will be considered as the lower and upper bounds of the carbonate circulation rate, respectively.Figure 3 shows that decreasing the carbonate section reboiler energy beyond 40 Gcal/hr will instantly cause the product gas to go off spec; therefore, this will act as the lower bound for this decision variable.To minimize waste of energy, the upper bound may be constrained at 50 Gcal/hr since at this value, all the sweet gas produced meets the required specification.
In a similar way, Figure 4 is considered to choose the upper and lower bounds on the circulation rate of diethanolamine (decision variable  3 ) as 90 and 110 m 3 /hr, respectively.From Figure 5, the lower and upper bounds on the amine section reboiler duty may lead 10 and 12 Gcal/hr, respectively.
Figures 6 and 7 show the variation of column stripping pressures with H 2 S and CO 2 composition in the sweet gas. Figure 6 shows that operating the carbonate stripper at pressures above 0.8 barg may lead to producing sweet gas which violates constraints, and to enable a moderate stripping efficiency, 0.7 barg is considered as the lower limit decision variable.
Figure 7 shows that the variation of the amine stripper pressure with the composition of hydrogen sulfide shows a reduction in H 2 S removal and exhibits a maximum increase at about 1 barg, and increasing the stripper pressure beyond this region rapidly increases the H 2 S removal.On the other hand, an increase in stripper pressure continuously decreases the amount of CO 2 in the sweet.Therefore, 0.9 barg and 1.2 barg are considered as the lower and upper limit on the decision variable, respectively.The lower limit pressure on the amine stripper is always expected to be greater than the upper limit for the carbonate stripper pressure since the overhead acid gas from the amine regenerator is directly feeding the carbonate stripper.

Process Optimization
The optimization of the gas sweetening unit involves two objectives: maximization of hydrocarbon recovery and minimization of the heat required to run the plant perunit time.This problem considers no variation in sour gas conditions (flow rate, temperature, pressure, and composition).

The mathematical expression of the objective functions of the optimization problem can be cast as
Maximize  1 () = Hydrocarbon Recovery, Minimize  2 () = Total Heat, Total Heat The first objective function represented by ( 6) is the percentage amount of sweet gas produced from the gas facility excluding CO 2 and H 2 S.This is also used as the measure of hydrocarbon content (C 1 to normal-C 5 ) since the feed gas contained no traces of mercaptans or another sulfur compounds.  and   are the mass flow rate of the sweetened gas and acid gas (CO 2 and H 2 S), respectively.  is the mass flow rate of the sour gas fed to the gas facility.Equation ( 7) represents the second objective function which is the total heat required to produce   kg of sweet gas in a given unit of time.
Gas sweetening plants require high amounts of energy supply to produce ultrasweet gas needed in LNG plants for safe process operation and high efficient hydrocarbon separation in the downstream operations.
To maximize the sweet gas purity (high hydrocarbon content), higher amounts of heat need to be applied to absorb a substantial amount of the acid gases.However, due to the high costs of energy to run the plant facility, energy minimization is often recommended to reduce the plant operating costs.Therefore, due to the two conflicting objectives, process operators often strive for the search of a trade-off between energy minimization and hydrocarbon (pure sweet gas) maximization.
Here, these two objectives are modeled and combined through the weighted sum method to form a single objective.The weights are selected interactively to produce preference information needed during the progressive optimization process.
The constraints of the optimization problem are Equations ( 8) and ( 9) provide constraints on the maximum composition of H 2 S and CO 2 allowable in the sweetened gas, respectively.Ultrasweet gas of utmost 1 and 50 ppmv of H 2 S and CO 2 , respectively, is used as feed to the liquefied natural gas plant.  and   in (10) are the upper and lower bounds on the various conditions of the process.These bounds are imposed on the K 2 CO 3 /DEA Circulation rate ( 1 ), carbonate section reboiler duty ( 2 ), DEA circulation rate ( 3 ), amine section reboiler duty ( 4 ), carbonate stripper pressure ( 5 ), and amine stripper pressure ( 6 ).From the simulation results and the sensitivity analysis of the previous section, appropriate bounds that define the search space are summarized in Table 3.These bounds are chosen based on a physical understanding of the problem so as not to compromise desired operating conditions.

Simulation Optimization
Framework.Simulation optimization entails finding optimal settings of the input variable(s), that is, values of  1 ,  2 , . . .,  6 , which optimize the output variable(s) of the simulation model.Optimizers designed for simulation embody the principle of executing separately from the model.In such a context, the optimization problem is defined outside the complex system.The output of the simulation model is used by an optimizer to provide new input values of the decision variables as the search for the optimal solution proceeds.On the basis of both current and past evaluations, the method decides upon a new set of input values.
Provided that a feasible solution exists, the optimization procedure ideally carries out a special search until an appropriate termination criterion is satisfied.A sequentiallike approach can be employed by coupling the solution of the optimization problem to an optimization routine (Figure 8).Every time the optimization routine needs a function evaluation by this approach, a call is made to the simulation model.In this work, optimization was performed through the barebone particle swarm optimization as discussed next.

Particle Swarm Optimization.
PSO is inspired by the social behavior of some biological organisms, especially the group's ability of some animal species to locate a desirable position at the given area.It was proposed first by Kennedy and Eberhart [32].In the PSO, a swarm consists of a set of particles, and each particle represents a potential solution of an optimization problem.Particles are placed in the search space of the function, and each evaluates the objective function at its current location.Each individual in the basic particle swarm is composed of , the dimensions of the optimization problem (number of decision variables), the current position   , the previous best position   , and the velocity update V  .Kennedy [36] conducted some experiments using PSO variants, which drop the velocity term and update a particle using Gaussian normal distribution of mean and standard deviation.The use of Gaussian normal distribution in the basic PSO algorithm to update the particle is known as bare-bones particle swarm optimization (BBPSO).Unlike the PSO, the BBPSO is parameter-free and is suitable for those real application problems where the information on parameters such as inertia weights and acceleration coefficients of particles is lacking or hard to obtain.
The bare-bones particle swarm optimization (BBPSO) proposed by Kennedy [36] Equation (11) shows that the position of each particle is randomly selected from the Gaussian distribution with the average of the personal best position and the global best position.
In BBPSO, the exploration search is facilitated by  ,  term, which tends to zero as the search progresses, and more emphasis will be put on the exploitation of the search.Kennedy [36] further proposed an alternative version of BBPSO also, referred to as BBExp in the present paper, where the th dimension of the th particle is updated as follows: Equation (12) justifies that since there is 50% chance that the th dimension of the particle changes to the corresponding pbest position, the BBExp inclines to search for pbest positions.Thus, BBExp is said to be biased towards exploiting the pbest positions as shown in Algorithm 1.

Results and Discussion
Figure 9 shows a Pareto set of optimal solutions obtained from the formulated problem (( 6)-( 12)).The figure shows that increasing heat of the process ( 2 ) will reduce the amount of H 2 S and CO 2 in the sweet gas ( 1 ), meaning that the gas will contain more of the hydrocarbons than the acid gases.Each point on the Pareto set is associated with a set of decision variables.A unique solution can be obtained using high-level decision making by management.The closeness of the Pareto set shows that a unique solution can also be easily generated using this method.
The percentage recovery of hydrocarbons may be interpreted as a small change along the hydrocarbon recovery axis.However, if translated into an equivalent amount of gas per unit time, the change is significant.That is to say, a 0.1% increase in hydrocarbon recovery is equivalent to approximately 9.4 tonnes of hydrocarbon per day.
Figures 10 and 11 show a comparison of the Pareto optimal solutions with the nonoptimal operating region.The nonoptimal operating conditions give 86.35% hydrocarbon recovery and 64.36 Gcal/hr of energy which shows a non proper region of operation.Applying this optimization model will approximately increase hydrocarbon recovery by 2.5% (∼235 tonnes of hydrocarbon/day) and minimize heat energy up to 16% of the non-optimal region.To obtain optimal operating conditions for a desired percentage recovery and the corresponding minimal energy from Figure 9, the decision variables of the respective optimal set can be read off from Figures 10 and 11.Another set of Pareto front is generated when each objective is plotted against all the decision variables at the optimal solutions.This study provides an opportunity to process operators to choose the operating conditions needed to optimize both objectives.show also a distinction between the optimal and non-optimal operating process variables of ADGAS plant, and its consistency is illustrated on all process variables.In these figures, the minimization and maximization objectives give a Pareto set below and above the non-optimal sets, respectively.Maximization of hydrocarbon recovery and minimization of energy are justified in the plots  1 () and  2 () versus  1 , . . .,  6 with the Pareto sets above and below the non-optimal set respectively.The difference in the optimal and non-optimal Pareto gives the amount of savings made on implementing this optimization approach.The Pareto set on the carbonate circulation rate ( 1 ) takes on values close to the upper limit of its bounds.The carbonate reboiler duty ( 2 ) also takes on optimal values near the upper bound of the constraint.Unlike decision variable  1 , the optimal values of  2 vary significantly between its midpoint and the upper limit.The optimal value of decision variable  3 takes values existing midway of the upper and lower bounds.Decision variables  4 ,  5 , and  6 span the entire feasible region.For all cases involving either  1 () or  2 (), all variables give rise to the same trend.

Conclusion
This work presented a study to optimize a natural gas sweetening process with the objectives of maximization Start up (i) Choose priori weights   (ii) Combine objective functions using weighted sum method (iii) Randomly generate an initial population: NP (iv)  = 0 Iteration For each particle,  = 1 to NP For  = 1 to   hydrocarbon recovery and minimization of heat energy required to run the process.ProMax simulation tool was used to model the process and the simulation results matched well the plant data.The multiobjective optimization using multiobjective bares-bone particle swarm indicates that a trade-off exists between the two objectives and plant operation can be optimized up to a 2.5% (∼235 tonnes/day) increase in hydrocarbon recovery and 16% decrease in heat required for the processoperation.The results also showed that the integration of an evolutionary algorithm in a process simulator can serve as an optimization tool box and can significantly improve the operation of the process.Future work will consider metamodeling of natural gas plants to overcome the large computational time involved in running integrated models of process simulators with evolutionary algorithms.
[35]Simulation Results.The process calculations were completed using ProMax V3.2[35].The electrolytic property package was used to predict the H 2 S and CO 2 absorption in both potassium carbonate and amine units of the Benfield HiPure process.The TSWEET kinetics model in ProMax was used to predict the CO 2 -amine/carbonate kinetic reactions taking place in all absorbers.TSWEET kinetics was developed by Bryan Research and Engineering for the purpose of accurately calculating the relatively slow absorption of CO 2 by amine solutions.Table

Table 2 :
Comparison of hot carbonate section with plant data.
Figure 2: Effect of K 2 CO 3 /DEA circulation rate on CO 2 and H 2 S composition in sweet gas.Figure 3: Effect of the reboiler duty of the carbonate section on CO 2 and H 2 S composition in sweet gas.Figure 5: Effect of the reboiler duty of the amine section on the CO 2 and H 2 S composition in sweet gas.
DEA circulation rate (m 3 /hr) Figure 4: Effect of DEA circulation rate on CO 2 and H 2 S composition in sweet gas.
Figure 7: Effect of amine stripper pressure on CO 2 and H 2 S composition in sweet gas.

Table 3 :
The upper and lower bound of the operating variables.
is a simple version of PSO.The BBPSO algorithm does not use the particle velocity but uses a Gaussian sampling based on best  and best   , and the position equation is replaced by Until termination criteria is satisfied End for Print Pgbest,   (Pgbest) Choose new weights   Print a Pareto-front Algorithm 1: The pseudocode of BBExp-PSO algorithm.