Modeling CO 2 , H 2 S, COS, and CH 3 SH Simultaneous Removal Using Aqueous Sulfolane–MDEA Solution

: In this study, a rate-based absorption model coupled with an improved thermodynamic model was developed to characterize the removal of acid components (CO 2 and H 2 S) and organic sulfur (COS and CH 3 SH) from natural gas with an aqueous sulfolane–MDEA solution. First, the accuracy of the thermodynamic model was validated by comparing the calculated partial pressure of CO 2 , H 2 S, and CH 3 SH with those of the experimental data reported in the literature. Then, the industrial test data were employed to validate the absorption model and the simulation results agreed well with the experimental data. The average relative errors of the removal rates of CO 2 , COS, and CH 3 SH are 3.3%, 3.0%, 4.1%, respectively. Based on the validated coupled model, the total mass transfer coefficient and mass transfer resistance of each solute component at different column positions were analyzed. The effects of the gas–liquid ratio, overflow weir height, and absorption pressure on the absorption performance of each component were studied, and the influence of the acid component concentration in the feed gas on the removal efficiency of methyl mercaptan (CH 3 SH) was also discussed. It is found that the improved absorption model can better characterize the absorption performance and be conducive to the optimal design of the absorber column.


Introduction
With the development of society, natural gas is increasingly used in industry and daily life due to its economic and environmental advantages [1][2][3].Since mined natural gas contains acid components (CO2 and H2S) and organic sulfur (COS and CH3SH), it needs further treatment before it can be used [4][5][6][7].Conventional gas processing includes distillation, adsorption, membrane separation, and absorption [8].Among these purification techniques, chemical absorption is the most commonly used method for acid-gas removal in the natural gas industries with its high efficiency and simultaneous strip of multiple acid gases [9].
In the industry, aqueous amine solutions are often used to remove acid components from natural gas [10][11][12].To overcome the limited capacity of organic sulfur components with amine solution, a mixed solvent containing sulfone is used because of the high gas loading for simultaneous absorption of high CO2 and H2S content in raw gas.A common sulfone-amine solution is composed of N-methyl diethanolamine (MDEA), sulfolane, and water, which has an advantage in treating natural gas with high organic sulfur content [8,13].Macgregor and Mather [14] reported the absorption of H2S and CO2 with a mixed solvent consisting of MDEA (20.9 wt.%) + sulfolane (30.5 wt.%).Jou et al. [15] investigated acid-gas absorption with MDEA, methanethiol, and ethanethiol at 40 and 70 ℃ at thiol partial pressure within 0.1-15.8kPa.Haghtalab et al. [16] measured the absorption of acid-gas in different mix solvents at 343 K and a total pressure of 0.1-0.21kPa and reported variation in the solubility of CO2 and H2S in different operating conditions and solvent compositions.However, although absorption with various mixed solvents has been reported, reliable thermodynamic modeling of the absorption process remains a big challenge for simulating and optimization of the acid-gas capture processes [17].
In order to accurately simulate and design the absorption process of the acid component and organic sulfur in natural gas using an aqueous sulfolane-MDEA solution, it is necessary to establish a reliable and robust absorption model, which involves the gas-liquid equilibrium, chemical reaction equilibrium, material and heat balance, and transfer properties of each component in the system [18][19][20].Afkhamipour and Mofarahi [21] compared the rate-based and equilibrium-stage models by simulating post-combustion CO2 capture using a 2-amino-2-methyl-1-propanol (AMP) solution in a packed column.The simulation results of the absorber show that the rate-based model can better predict temperature and concentration curves than the equilibrium phase model.Al-Baghli et al. [22] simulated the process of removing CO2 and H2S by MEA and DEA aqueous solutions using rate-based gas absorber model and obtained reliable simulation results.Pacheco and Rochelle [23] developed a framework to perform selective absorption of H2S using MDEA solution from a gas stream containing CO2.The Maxwell-Stefan and enhancement factor theories are used in the model.Mandal and Bandyopadhyay [24] conducted theoretical and experimental research on the simultaneous absorption of CO2 and H2S into solutions containing MDEA and DEA.Moioli et al. [25] used the Eddy diffusivity theory in Aspen Plus and used an external subroutine to simulate the absorption of CO2 and H2S from gas streams.In addition, they modified the parameters for vapor-liquid equilibrium (VLE) calculations and verified the simulation using data in the literature.Yang et al. [26] used an aqueous sulfone-MDEA solution to remove organic sulfur in natural gas and studied the factors that affect organic sulfur removal in natural gas purification devices.However, the process of the simultaneous removal of multiple impurities (i.e., CO2, H2S, COS, and CH3SH) is usually very complicated, and the influence of acidic gas on the absorption performance of organic sulfur during the removal process is still unclear.
In this study, we developed a rate-based absorption model coupled with an improved thermodynamic model to characterize the removal of acid components (CO2 and H2S) and organic sulfur (COS and CH3SH) from natural gas with an aqueous sulfolane-MDEA solution.An improved thermodynamic model was used considering the influence of acid component concentration on the removal efficiency of methyl mercaptan (CH3SH).The mass transfer characteristics at different column positions were analyzed.The influence of the gas-liquid ratio, overflow weir height, and absorption pressure on the removal rate of each component and the influence of acid component content on the removal rate of methyl mercaptan (CH3SH) were studied.Furthermore, with the modified rate-based absorption model, we successfully simulated the adsorption process for a wide range of feed gas compositions and operating conditions.This study is expected to further facilitate optimization of the operating conditions and device structure.

Thermodynamic Framework
The detailed thermodynamic model consists of the gas-liquid equilibrium and chemical reaction equilibrium of each component in the solvent.

Gas-Liquid Equilibrium
The gas-liquid equilibrium of the acid component (CO2 and H2S) and organic sulfur (COS and CH3SH) in the aqueous sulfolane-MDEA solution is calculated by Henry's law: where P is the system pressure, Pa. is Henry's law constant of component i in the mixed solvent of water, sulfolane, and MDEA, Pa.i y is the mole fraction of compo- nent i in the vapor phase, and i x is the equilibrium mole fraction of component i in the liquid phase.i ϕ is the fugacity coefficient of component i in the vapor phase, which is calculated using the Peng-Robinson equation of state (EOS) [27].
γ is the unsymmet- ric activity coefficient of component i in the mixed solvent solution of water, sulfolane, and MDEA.
* i γ is normalized to the mixed solvent infinite dilution reference state.
In the mixed solvent, Henry's law constant can be calculated from those in the pure solvents, as shown in Equation ( 2) [28].

* ln( ) ln( )
where is the weighting factor, which can be calculated according to the method in the literature [13].iA γ ∞ is the infinite dilution activity coefficient of component i in pure       During the absorption process, MDEA reacted with acidic components and the amount of reaction affected the absorption of organic sulfur in the aqueous sulfolane-MDEA solution.Thus, the unreacted ratio of MDEA is introduced to modify Henry's coefficient of methyl mercaptan (CH3SH) in the aqueous sulfolane-MDEA solution [26].Equation (3) is used to quantitatively describe the influence of the acid component content on Henry's constant for CH3SH.
where f is the unreacted ratio of MDEA, which is the ratio between unreacted MDEA and the total amount of MDEA and MDEAH + in aqueous solutions.c is the effect factor, symbolizing the influence of the acid component on Henry's law constant.The electrolyte-nonrandom two-liquid (e-NRTL) model is used to calculate the activity coefficients [33].The molecule-molecule, electrolyte-electrolyte, and molecule-electrolyte binary parameters are mainly obtained from the literature [10,13].Some binary parameters that are not mentioned in the literature are defaulted to (8, -4), and the non-randomness factor is fixed at 0.2 [34].

Aqueous Phase Chemical Equilibrium
In aqueous solutions, the acid components H2S and CO2 react with MDEA, and the ionic equilibrium reactions are expressed as Equations ( 4)-( 9) [35].

H O MDEAH H O MDEA
2H O H O OH The liquid phase compositions must satisfy the chemical equilibrium relationships [11]: The chemical equilibrium constant of reaction j can be calculated by Equation (11) [11,35].where the constants A, B, C and D are summarized in Table 2. a.The equilibrium reaction order of CO2 and H2S are obtained by fitting the experimental data.The reaction orders of CO2 and HCO3 -are 1.25, the reaction orders of H2S and HS -are 1.35, and the reaction orders of the remaining reactions are 1.

Rate-Based Model Assumptions
In order to simulate the absorption process of removing the acid component (H2S and CO2) and organic sulfur (COS and CH3SH) from natural gas utilizing a tray column, a mathematical model was established using the two-film theory [37].The basic assumptions are as follows: (1) Consider both physical absorption and chemical absorption, and the reaction only occurs in the liquid phase.
(3) The heat transfer resistance of the liquid phase is neglected, and the interface temperature is equal to the bulk temperature.
(4) Flow is one-dimensional in the axial direction, and the radial temperature and concentration variation can be neglected.

Material and Energy Balance
The tray column is composed of many stages, where vapor and liquid streams from adjacent stages contact each other and exchange mass and energy at their common interface [18,21].
Figure 1.Presents a stage j in the tray column, and i represents the components.The main mass and energy balance equations are listed below [38].
Material balance for liquid phase: Material balance for vapor phase: Energy balance for liquid phase: Energy balance for vapor phase: where F is the mole flow rate of feed, kmol s -1 .L and V represent the mole flow rates of liquid and vapor, respectively, kmol s -1 .N is the mole transfer rate, kmol s -1 .r is the reaction rate, which can be calculated from the component concentration.i x and i y are the mole fractions of component i in the liquid and vapor phases, respectively.
Q is the heat input to a stage, J s -1 ; q is the heat transfer rate, J s -1 ;

FL H and FV H
are the enthalpies of the inflow liquid and inflow vapor, J kmol -1 ; and the thermodynamic properties including enthalpy and heat capacity utilized in heat transfer calculations can be obtained directly from the literature [10][11][12][13].

Mass Transfer and Enhancement Factor
The mass transfer flux between the gas phase and the liquid phase is shown in Equation ( 16) [39].* , (P P ) a where , G i K is the overall mass transfer coefficient for component i; p a is the effective mass transfer area per unit area of the tray, m 2 m -2 ; c A is the cross-sectional area of the column, m 2 ; i P is the partial pressure of component i in the gas bulk; and is the partial pressure of component i in equilibrium with the liquid phase.The partial pressure of each component can be calculated by the established thermodynamic model.
Using the two-film model, the overall mass transfer coefficients incorporating gas and liquid transfer resistance for component i is expressed as Equation (17).
where is Henry's law constant of component i; , G i K is the overall mass transfer of component i in gas phase; and k are the mass transfer coefficient of com- ponent i without reaction in the gas phase and the liquid phase, which are calculated according to the method used by Simon et al. [18] and Saimpert et al. [40].The density, viscosity, and surface tension required for the calculation can be obtained directly from the Aspen database.i E is the enhancement factor of component i due to chemical reac- tions, which is the ratio of the mass transfer enhanced by a reaction over the mass transfer without the reaction [21].The calculation of the enhancement factor in this work uses the formula developed by Danckwerts [41].
For large values of enhancement factor for infinite fast reaction ( E ∞ ) and Hatta numbers lower than 2, the enhancement factor is calculated using a simplified Equation (18) [42].

tanh(Ha)
i Ha E = (18) For infinite fast reaction ( E ∞ ) values lower than 100 and Hatta numbers larger than 2, the enhancement factor is calculated using a simplified Equation ( 19) [1].
where the Hatta number and enhancement factor for the infinite fast reaction are expressed as Equations ( 20) and ( 21) [42].

, ,
, Where , i L D refers to the diffusivity of absorption component i in an aqueous sul- folane-MDEA solution, is the diffusivity of the MDEA in the liquid phase.
is the concentration of absorption component i at the gas-liquid interface, and represents the concentration of MDEA in the bulk liquid.Pacheco and Rochelle [23] and Al-Ghawas et al. [32] measured the rate constants of reaction between CO2 and MDEA and between COS and MDEA in the mixed solvent, respectively, as shown in Equations ( 22) and ( 23).
For CO2 and COS, the enhancement factor is needed for the mass transfer rate calculations, which are calculated using Equation (18).As the reaction between H2S and MDEA is very fast, the enhancement factor can be calculated using the infinite fast reaction rate, as shown in Equation (19).In contrast, the absorption of CH3SH is mainly a physical effect, so there is also no need to consider the influence of chemical reactions.

Computational Implementation
The model calculation was performed by writing Fortran programs.The solution of the absorption model adopts the stepwise calculation method, and the steps are as follows: (1) The concentration of the absorbed components in the gas and liquid feed, and the gas-liquid ratio are known.By assuming the total removal rate of each component, the concentration of each absorbed component in the liquid at the last plate can be calculated.
(2) Since the concentration of gas feed and liquid discharge of the last plate is known, the concentration of liquid feed and gas discharge of this plate, that is, the concentration of gas feed and liquid discharge of the upper plate can be calculated.
(3) It is calculated upward until the gas discharge concentration of the top plate is obtained, then the removal rate of the whole column of each component can be calculated.
(4) Compare the removal rate of the whole tower between the hypothesis and the calculated one.If the error is greater than 0.0001, change the hypothesis and continue the iterative calculation until the error requirement is met.

Modeling Results
As mentioned above, a rate-based absorption model coupled with the thermodynamic model was developed for the removal of acid components (CO2 and H2S) and organic sulfur (COS and CH3SH) from natural gas with an aqueous sulfolane-MDEA solution.
In this section, the accuracies of the thermodynamic model and absorption model are first validated.Then, the effects of the gas-liquid ratio, overflow weir height, and absorption pressure on the absorption performance of each component are investigated, and finally, the influence of acid component concentration in the feed gas on the removal efficiency of methyl mercaptan (CH3SH) is also discussed.

Thermodynamics Model Validation
An improved thermodynamic model was employed to calculate the solubility and partial pressure of acid components and organic sulfur in the aqueous sulfone-MDEA solution.The experimental data of several systems are employed to verify the accuracy of the thermodynamic model.
Figures 2 and 3 display the experimental and calculated CO2 partial pressure data for the CO2-H2O-MDEA system and H2S partial pressure data for the H2S-H2O-MDEA system at different compositions and different temperatures.As shown, the calculated CO2/H2S partial pressure values are in good agreement with the experimental data reported by the literature [14,[43][44][45], thus validating the proposed thermodynamic model.From Figure 2, it is found that the CO2 partial pressure is more sensitive to temperature.With the temperature increasing from 313.15 K to 338.75 K, the CO2 partial pressure experiences a significant increase.However, the MDEA mass fraction seems to have little effect on the CO2 partial pressure when the temperature is 313.15K. Figure 3 shows the comparison of experimental data and calculated data in terms of H2S partial pressure for the H2S-H2O-MDEA system.As shown, the proposed thermodynamic model is also effective and robust.This means that the improved thermodynamic model has good applicability.When compared with an aqueous amine solution, an aqueous sulfone-amine solution is more effective, especially for raw gas with high organic sulfur content.A common sulfone-amine solution is composed of MDEA, sulfolane, and water.Herein, two systems, namely the CO2-H2O-sulfolane-MDEA system and the H2S-H2O-sulfolane-MDEA system, are used to validate the thermodynamic model.The comparison results of the calculated values and experimental data in terms of CO2/H2S partial pressure are given in Figures 4 and 5.As shown, the calculated results agree well with the experimental data in terms of CO2/H2S partial pressure.Generally, the addition of sulfolane is beneficial to enhancing the absorption of acidic gases, and therefore, high acidic gas loading is observed at low partial pressure.When compared with Figures 4  and 5, it can be found that H2S has a higher solubility than CO2 in the H2O-sulfolane-MDEA solution.In order to further verify the reliability of the thermodynamic model, some CO2 and H2S partial pressure data calculated in this paper are compared with the values calculated with reaction parameters proposed by Austgen et al. [36].The average relative deviations of the partial pressure data of CO2 and H2S calculated in this paper are 10.7% and 24.4%, respectively, while the average relative deviations of the partial pressure data calculated with reaction parameters proposed by Austgen et al. [36] are 38.1% and 36.8%,respectively.The calculated results using the improved thermodynamic model are not only in good agreement with the experimental data but also better than the cal-culated values of the existing models in the literature (see Figures 6 and 7).This is because we refitted the exponents of the equilibrium equations of reactions ( 4) and ( 5) and calculated the reaction equilibrium constants with the new parameters in Table 2.  To verify the prediction performance of the improved thermodynamic model on the partial pressure data of CH3SH in the system, the calculated value of the model is compared with the value measured by Jou et al. [15] in Figure 8.It can be seen that the predicted CH3SH partial pressure data are in good agreement with the experimental value.This shows that the improved thermodynamic model not only has higher prediction accuracy but also has a larger application range.In short, the proposed thermodynamic model has strong applicability, and its accuracy has been verified by various experimental data.In this study, the thermodynamic model is employed to calculate the solubility and partial pressure of acid components and organic sulfur in the aqueous sulfone-MDEA solution.

Rate-Based Absorption Model Validation
Before the rate-based model is applied to investigate the absorption process, its accuracy must be verified.The data were collected on an industrial sieve plate tower of 3.4 m in diameter with a mixed solvent comprising 20 wt% H2O, 40 wt% sulfolane, and 40 wt% MDEA.The tower was operated under a different number of plates (22, 26, and 30) and overflow weir heights (0.1 and 0.15 m).The Inlet gas-liquid ratio (vol./vol.) was adjusted from 490.2 to 720.0.CO2 and H2S loading were varied from 4.11vol% to 4.73 vol% and 1.34 vol% to 1.57 vol%, respectively.The concentration of CH3HS and COS were changed from 14.12 to 25.41 mg/m 3 and 11.18 to 34.58 mg/m 3 , respectively.As shown in Table 3, the simulation results of removal efficiency for absorption components are in good agreement with the data from the industrial absorber.For CO2, COS, and CH3SH, the average relative errors are only 3.3%, 3.0%, and 4.1%, respectively.Therefore, the rationality of the absorption model can be confirmed.

Rate-Based Absorption Model Calculation
The specifications and operating parameters of the absorption column are summarized in Table 4, and the mass fractions of H2O, sulfolane, and MDEA in the mixed solvent are 0.20, 0.40, and 0.40, respectively.The parameters of the absorber are obtained from (the desulfurizing column of) PetroChina Southwest Oil and Gasfield Company, and the operating parameters are obtained by analyzing the operating parameters of typical industrial absorbers and by making appropriate extensions.The concentration profile of each component and the gas-liquid temperature profile are also obtained.As shown in Figure 9, the absorption rate of H2S is fast, and the concentration of H2S in the feed gas drops rapidly after entering the absorption column; it is completely removed in the middle of the column.This is the result of the synergistic effect of physical absorption and chemical reaction.However, the absorption rate of CH3SH is low.After the feed gas enters the column, the concentration of CH3SH decreases slowly and the final removal rate is low.The main reason for this is that the physical adsorption and absorption processes are controlled by equilibrium.Physical absorption and chemical reactions are usually accompanied by a thermal effect, which changes the system temperature and, in turn, affects the absorption of components.Therefore, the temperature of the absorber column needs to be strictly controlled.Figure 10 shows the gas and liquid phase temperatures in the column at different concentrations of H2S and CO2.It can be seen that the gas phase and the liquid phase exchange heat in the column, and as the concentration of H2S and CO2 increases, the heat release during the absorption process increases, resulting in an increase in the outlet temperature of the liquid phase at the top of the column.When the feed composition is the same, the H2S is completely absorbed in the middle section of the column, and then the heat release is rapidly reduced, resulting in a rapid drop in the gas phase temperature at the top of the column.The calculated total mass transfer coefficient and mass transfer resistance of each solute component at different column positions are shown in Figure 11.
For CO2 and COS, the gas mass transfer resistance can almost be ignored.In the whole absorption column, the gas mass transfer coefficient is less than 0.01 kmol/m 2 /kPa.For H2S, however, the gas mass transfer coefficient is almost ten times greater than that of CO2 and COS.Additionally, the liquid mass transfer coefficient of H2S is also much higher than that of CO2 and COS, accounting for 15-25% of the total mass transfer resistance, so the gas-side mass transfer resistance cannot be ignored.The gas-side mass transfer resistance of CH3SH accounts for 5-10% of the total mass transfer resistance, which also needs to be considered in the calculation process.This is consistent with the calculation method of the absorption model.
Based on the mass transfer characteristics analysis, the rate-limiting steps of absorption processes for different components can be clarified.Furthermore, it is helpful to simplify the model calculation.

The Influence of the Operating Parameters
Based on the validated thermodynamic-absorption coupled model, the effects of gas-liquid ratio, overflow weir height, and absorption column pressure on the removal efficiency of each component were studied, and the results are shown in Figures 12,13 and 15.In addition, the effect of the acid component content in the feed gas on the removal efficiency of methyl mercaptan was studied, and the results are shown in Figure 14.It can be seen from Figure 12 that the gas-liquid ratio has a greater impact on the absorption of CH3SH because CH3SH is physically absorbed and controlled by equilibrium.When the gas-liquid ratio increases from 400 to 800, the removal efficiency of CH3SH decreases from 0.8 to 0.4.Obviously, a more liquid phase is conducive to the physical absorption process.Since the absorption rate of CO2 and COS is determined by mass transfer rate, it is less affected by the gas-liquid ratio.Figure 13 shows the effect of the overflow weir height on the absorption performance of each component.As the reaction rate of CO2 in the absorption solvent is very fast, the absorption rate is mainly controlled by the mass transfer rate of CO2.As the overflow weir height increases, the liquid thickness over the tray increases, and, thus, the mass transfer area increases, resulting in an increase in the removal rates of CO2 and COS.For CH3SH, the overflow weir height has little effect on its removal rate.One reason is that the absorption rate of CH3SH is controlled by its physical solubility in the mixed solvent.The second reason is that the increase in the overflow weir height strengthens the absorption of CO2 and H2S, which reduces the solubility of CH3SH in the liquid phase.Figure 14 shows the effect of acid component content in the feed gas on the removal rate of CH3SH.As the acid component content increases, the absorption rate of CH3SH decreases significantly.The reason for this is that, as the acid component content in the liquid phase increases, the Henry constant of CH3SH absorption increases significantly.The influence of tower pressure is shown in Figure 15.As the absorption pressure increases, the absorption driving force increases, and the removal efficiency of each component increases.Additionally, the physical absorption process is more sensitive to the tower pressure.With the increase in tower pressure, the removal efficiency of CH3SH achieves a significant increase when compared with the removal efficiency of CO2 and COS.Due to the chemical absorption being the rate-limited step for the removal of CO2 and COS, the absorption pressure on their removal efficiency is insignificant.As shown in Figure 15, the removal efficiency of COS increases by less than 10% when the absorption pressure increases from 4 MPa to 8 MPa.

Conclusions
In this study, a rate-based absorption model coupled with a thermodynamic model was developed to characterize the removal of acid components (CO2 and H2S) and organic sulfur (COS and CH3SH) from natural gas with an aqueous sulfolane-MDEA solution.The thermodynamic model was improved by improving the calculation method of the chemical equilibrium constant and by incorporating the influence of acid components for COS removal, which made the gas-liquid equilibrium (GLE) data more accurate.The absorption model was validated by the experimental data obtained from an industrial device, and the average relative errors of the removal rates of CO2, COS, and CH3SH obtained by experiment and calculation are 3.3%, 3.0%, and 4.1%, respectively.The validated coupled model indicated that the gas mass transfer resistance can almost be ignored for CO2 and COS, but both gas and liquid mass transfer resistance for H2S and CH3SH should be considered.Additionally, the analysis of influencing factors shows that the gas-liquid ratio has a greater impact on the physical absorption process controlled by equilibrium but has a limited effect on the absorption process determined by a chemical reaction.An increased overflow weir height increases the gas-liquid contact time and mass transfer area, thus resulting in an increase in the removal rates of CO2 and COS, which, however, inhibits the absorption of CH3SH due to the increased Henry coefficient.Furthermore, the physical absorption process rather than the chemical absorption process is more sensitive to the tower pressure.

H
is Henry's constant of component i in pure solvent A, and the values of iA H of solute i (CO2, H2S and CH4) can be obtained directly from the literature [10-13].Some iA H values of CH4, CH3SH, and COS regressed according to the experimental data [29-31], which are summarized in s constants of CH3SH in the mixed solvent when the influence of acid component is not taken into account and is taken into consideration, respectively.

K
where j is the chemical equilibrium constant of reaction j, is the mole fraction of reactant component i in reaction j, ' i x is the mole fraction of product component ' i in reaction j, and i γ and ' i γ are the unsymmetric activity coefficients of reactant com- ponent i and product component i' in the aqueous solution.The unsymmetric activity coefficients are normalized to the aqueous phase infinite dilution reference state.and ' i v represent the stoichiometric coefficients of reactants and products, respectively.

Figure 6 .
Figure 6.Comparison of the calculated data in the reference and in this work for the CO2-H2O-MDEA system.■, experimental data measured by Jou et al. [43]; -, data calculated with the reaction parameters proposed by Austgen et al. [36]; and -, data calculated in this work.

Figure 7 .
Figure 7.Comparison of the calculated data in the reference and in this work for the H2S-H2O-MDEA system.■, experimental data measured by Jou et al. [43]; -, data calculated with the reaction parameters proposed by Austgen et al. [36]; and -, data calculated in this work.

Figure 14 .
Figure 14.The effect of acid component content in feed gas on the removal efficiency of CH3SH (number of plates: 25; overflow weir height: 0.15 m; absorption pressure: 6 Mpa; inlet gas flow rate: 5.0 × 10 6 Nm 3 /d; and inlet gas-liquid ratio: 600).

i Solvent A Data source
. Parameters for Henry's constant.Solute

Table 2 .
Parameters for temperature-dependent mole-fraction reaction equilibrium constants.

Table 3 .
Comparison of experiment data and model calculation data from the industry.

Table 4 .
Specifications of the absorber column and typical operating data.