Prediction of sulfur content in propane and butane after gas purification on a treatment unit
National Iranian Gas Company (NIGC), South Pars Gas Complex (SPGC), Asaluyeh, Iran
* Corresponding author: email@example.com
Accepted: 7 June 2018
The acidic compounds such as Mercaptans, H2S and COS are commonly present in the liquid LPG streams in the south Pars gas processing plant. Sulfur contaminants not only lead to odor problems but can form objectionable oxides on combustion and cause environmental pollution. In present study, Support Vector Machine (SVM) is employed to develop an intelligent model to predict the sulfur content of propane and butane products of Liquefied Petroleum Gas (LPG) treatment unit of south Pars gas processing plant of Assaluyeh/Iran. A set of seven input/output plant data each consisting of 365 data has been used to train, optimize, and test the model. Model development that consists of training, optimization and test was performed using randomly selected 70%, 15%, and 15% of available data respectively. Test results from the SVM developed model showed good compliance with operating plant data. Squared correlation coefficients for developed models are 0.97 and 0.99 for propane and butane sulfur content, respectively. According to the results of the present case study, SVM could be regarded as a reliable accurate approach for modeling the sulfur content of LPG treatment unit of a natural gas processing plant.
© H. Adib et al., published by IFP Energies nouvelles, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Liquefied Petroleum Gas (LPG) referred to predominately propane or butane, either separately or in mixtures, which is maintained in a liquid state under specific pressure/temperature within the confining vessel (Santos et al., 2016). LPG is a valuable energy source that is used worldwide for numerous business applications in industry and transportation. The largest market for LPG is the domestic/commercial market, followed by the chemical industry where it is used as a petrochemical feedstock and the agriculture industry (Safadoost et al., 2014).
Valuable LPG is a natural gas processing by-product in south Pars gas complex refineries and stored as a liquid in atmospheric pressure tank. LPG delivered to customers as single-phase pressurized liquid products and should meet some specifications for sales as it is shown in Table 1 (Moaseri et al., 2013; Asil and Shahsavand, 2014). LPG is treated to reduce total sulfur content to meet sweetness specifications. Sulfur may be presented as hydrogen sulfide, carbonyl sulfide, carbon disulfide and mercaptan. All forms may be present in the same liquid. Sulfur contaminants not only lead to odor problems but can form objectionable oxides on combustion and cause environmental pollution (Safadoost et al., 2014; Mahdipoor and Ashkezari, 2016). H2S absorption into NaOH solution is one of the main methods to for H2S removal (Bashipour et al., 2017; Sharifi and Omidbakhsh Amiri, 2017). The treatment process for removal of H2S, mercaptan and elemental sulfur follows techniques and philosophies that have been well defined over years and will discussed in details in subsequent section.
LPG product specifications.
H2S absorption into NaOH solution is one of the main methods to produce sodium sulfide (Na2S) and H2S removal. Other methods to produce Na2S are reduction of sodium sulfate (Na2SO4) by solid carbonaceous materials, reduction of Na2SO4 by gaseous reducing agents, exchange decomposition of barium sulfide (BaS) with sodium sulfate, carbonate, and hydroxide as well as an electrolytic method
Notable attempts to develop accurate practical models for complicated chemical processes have been carried out with the aim of minimizing operational costs. In recent years, application of modeling methods which deal with input/output data of industrial plants have received considerable attention. Support Vector Machine (SVM) has been emerged as a proven technology which offers an alternative way to address problems with no specific relationship between input and output parameters. The main advantage of such models over existing approaches is the capability of learning and generalizing data, fault tolerance and inherent contextual information processing in addition to fast computation potential (Raynal et al., 2016). Such characteristics make them perfect candidates for applications where the complexity of the data or task demands high computational costs (Haghbakhsh et al., 2012; Adib et al., 2013, 2015; Moradi et al., 2016).
In this study, SVM model is developed to determine the output variables of south Pars natural gas processing plant. Since the purpose of the process is to reduce the sulfur content of propane and butane product of LPG treatment unit of south Pars gas processing plant of Assaluyeh, the input parameters are amine, caustic and feed flowrate of this unit and the output variables are total sulfur of propane and butane of this unit. A large dataset of these variables are gathered from the plant and introduced to the algorithm. The models are then compared to actual plant data and with each other and then the accuracy of the models is assessed through calculation of Average Absolute Deviation Percent (ADD%).
The acidic compound, Mercaptans, H2S and COS are commonly present in the liquid LPG streams in the south Pars gas processing plant. Due to the nature of the upstream process unit, the liquid butane stream will typically contain mainly Mercaptans. However, the liquid propane stream will typically contain Mercaptans and reasonable amount of H2S and COS. Sulfur compound concentration in LPG is shown in Table 2.
Sulfur compound concentration in LPG of south Pars gas processing plant.
Liquid propane is treated first in an amine treating unit to remove H2S and COS to acceptable level. In the amine unit, H2S and COS are removed from propane using 21% of DEA (Di-Ethanol-Amine) solution as solvent. Amine solvents are very often used for natural gas deacidification purposes as they can be adapted to various specifications and to a wide range of feed gas compositions (Magné-Drisch et al., 2016). Amine section consists of an extractor column for H2S removal and a COS removal section with mixer-settlers. As indicated in Figure 1 the propane feed originating from the H2S extraction column is led to three mixer-settler combination in series which will remove COS from the propane stream. This line up together with sufficient mixing and resistance time in two stages enables maximum COS removal. Expected H2S and COS levels in propane are 1–2 ppm wt. H2S sulfur and 1–2 ppm wt. COS sulfur. Then propane routed to caustic extraction section for mercaptan removal.
Amine process flow diagram.
Both propane and butane are fed to a caustic extraction column for removal of mercaptan. In both extractor column for propane and butane, the LPG streams are contacted counter-currently with caustic in a column filled random packing. The rich caustic coming from both extractors is combined and sent to the regeneration section. The extractors are designed for maximum achievable mercaptan removal. The overall reaction for extraction and regeneration for the Mercaptans (R-SH) is expressed below:(1) (2)
Caustic process flow diagram is shown in Figure 2.
Caustic extraction section.
Amine and caustic flow rate are two important variables which affect the total sulfur in LPG product. Amine extractor column run with a flow of amine solvent which is well in excess of the minimum required for H2S removal. This is because the minimum solvent rate is more typically set by minimum flows for wetting and providing interfacial surface area. Also in mixer-settler, to ensure sufficient contact between liquid hydrocarbon and the solvent, the volumetric ratio between liquid hydrocarbon and the solvent must not be changed noticeably. The lower solvent flowrate will result in poor COS extraction and a much higher ratio may lead to a reversed phase mixture. In caustic extraction section, an increase in the caustic circulation flow rate leads to a reduction of mercaptan sulfur in LPG but disulfide oil level in the LPG will increase. On the reverse, a reduction in caustic circulation flow rate means a reduction in disulfide oil but the mercaptan sulfur in LPG will increase. This means that there is an optimum caustic circulation flow rate in order to minimize the total LPG sulfur.
Support Vector Machine introduced first by Vapnik (1998), like Artificial Neural Networks (ANN), is an intelligent learning approach equipped learning algorithm that analyzes data and find patterns of input/output data. Support Vector Machine training procedure converges to optimum output results faster and it is not need to control model parameters (Cortes and Vapnik, 1995; Pelckmans et al., 2002; Suykens et al., 2002; Curilem et al., 2011). For detailed information about the SVM refer to our previous work (Haghbakhsh et al., 2012; Adib et al., 2013, 2015; Moradi et al., 2016).
Pattern recognition or classification can be performed by SVM in a data set consisting of N data point where x k is a p-dimensional vector and y k can get one of the two values, either +1 or −1 (i.e., indicating the class to which the point x k belongs. In their basic form, they learn a linear hyperplane that separates a set of positive samples from a set of negative samples with maximum margin. Consider Figure 3 which shows two possible splitting hyperplanes and their related margins. Both hyperplanes can appropriately categorize all the given data. However, we expect the hyperplane with the larger margin to be more accurate in classifying new data than the hyperplane with the smaller margin. This is the reason that SVM searches for the hyperplane with the largest margin (Zaidi, 2015).
Support vector machine classification.
A separating hyperplane can be written as w · x – b = 0 (Agarwal et al., 2008; Yélamosa et al., 2009), where w is the normal vector to the hyperplane and b represents the offset of the hyperplane from origin that is referred to as bias. The offset along the vector w from the origin can be determined by b/||w||. As shown in Figure 4, for the cases that the training data are linearly separable, two hyperplanes can separate the data in a way that there are no data points between them. Obviously these hyperplanes can be described as:(3) (4)
By using geometry, one can show that distance between these two hyperplanes is 2/||w||, so the problem of ||w|| minimization is required to maximize hyperplane margin. It is also required to prevent data points from falling into the margin, and other necessary constraints are imposed as (Yélamosa et al., 2009):(5) (6)that can be rewritten as (Yélamosa et al., 2009):(7)
Constraint minimization of ||w|| is thus required to develop an ideal classifier. Such minimization problem is difficult to solve, however it is possible to substitute 0.5 ||w||2 instead of ||w|| in problem. Chiang et al. (2004) showed that minimization problem can be formulated as:(8)where is Lagrangian multiplier that helps in finding the local minimum or maximum of a function (Mehdizadeh and Movagharnejad, 2011). The problem of equation (8) can be solved by standard quadratic programming techniques that results in finding normal vector to the hyperplane as presented in equation (9):(9)
Input/output SVM model with the general form of y = f(x) takes the form of equation (11) in feature space (Eslamimanesh et al., 2012; Kulkarni et al., 2005):(10)where f(x) represents output vector and K(x, x k) is the kernel function calculated from the inner product of the two vectors x and x k in the feasible region built by the inner product of the vectors Ф(x) and Ф(x k ) as follows (Eslamimanesh et al., 2012):(11)
Among choices for Kernel function the Radial Basis Function (RBF) Kernel that is used extensively (Zhao, 2009; Ding et al., 2012) has been applied in this work that is presented in equation (12),(12)where σ is kernel parameter to be determined by an external optimization algorithm during the internal SVM calculations. Bias, b, is usually determined by using primal constraints as (Kulkarni et al., 2005):(13)
Lagrangian multipliers, α i , can be calculated by solving following quadratic programming problem (Terzica et al., 2010):(14)
Subject to constraintswhere γ is regularization parameter and controls the tradeoff between complexity of the SVM model and the number of non-separable points. This compact formulation of quadratic optimization has been proved to have a unique solution (Agarwal et al., 2008). In conclusion, the SVM takes the form of the constrained optimization problem of equation (15) in order to obtain the optimum value of γ (Vapnik, 1998; Zanghirati and Zanni, 2003; Agarwal et al., 2008):(15)
Subject towhere ε is the precision threshold and , represent the slack variables with nonnegative values to ensure feasible constraints. The first term in equation (15) represents model complexity while the second term represents the model accuracy or error tolerance. The Mean Square Error (MSE) and Mean Absolute Error (MAE) as defined by equations (16) and (17) are used to calculate prediction error of the developed SVM model.(16) (17)where O i is the simulation results of SVM model, T i represents real time plant data of the natural gas sweetening plant and n denotes the number of the data used for model evaluation. Figure 5 presents the SVM model algorithm in flowchart format.
Schematic flowchart of SVM algorithm.
The gas processing plant under study in this work, is located in south Pars gas field, in Asaluyeh/Iran. A data set of seven series of input/output data is collected from the LPG treatment unit. Each data series consists of 365 data points of the plant under normal operating conditions in span of one year. All data series are scattered in a wide range for which the maximum and minimum numerical values are presented in Table 3. In order to estimate qualitative correlations between these input/output plant data, Figures 6 and 7 are depicted for better visualization. Figures 8–10 show sulfur content of LPG products, which are output variable of the LPG treatment unit, versus caustic, amine and LPG flowrate of this unit.
Model parameters of propane 3D graph.
Model parameters of butane 3D graph.
3D illustration of caustic and amine flowrate for propane sulfur content.
3D illustration of caustic and propane flowrate for propane sulfur content.
3D illustration of caustic and amine flowrate for butane sulfur content.
The range of operating plant data used for model development.
Since the input data cannot be changed systematically during normal plant operation, it is difficult to find relevant relationships between input and output variables. Therefore, data mining is performed to demonstrate the effect of varying inputs on process outputs (Adib et al., 2013). Figures 11 and 12 illustrate some of the input/output data for the LPG treatment unit. As can be seen in Figure 11, the two vertical axes show the ratio of amine and caustic flowrate to propane flowrate of LPG treatment unit. By increasing the amine and caustic flowrate to extractor column, total sulfur of propane decreases significantly. In this case, the lower solvent flowrate will result in poor COS extraction and a much higher ratio may lead to a reversed phase mixture. Therefore, for mercaptan removal the optimum ratio could be regarded as 0.12, and for amine extractor the same ratio could be regarded as an applicable ratio. As indicated in Figure 12, this ratio could be set as 0.65 for the ratio of the caustic to butane flow rate of LPG treatment unit.
Effect of amine and caustic flowrate on sulfur content of propane.
Effect of caustic flowrate on sulfur content of butane.
In SVM model, the two key parameters are regularization parameter (γ) and kernel parameter (σ 2) which determines the tradeoff between the fitting error minimization and the smoothness of the estimated function. Optimum numerical values of these two parameters are calculated using Genetics Algorithm which is applied in the SVM Matlab codes. The details of GA optimization procedure are presented by Adib et al. (2013). The optimization procedure has been repeated several times in order to guarantee that the developed model’s parameters are very close to optimum results. The optimum values of γ and σ 2 were reported in Table 4.
The optimum values of the SVM model parameters for output variables.
The operating plant data collected over the span of one year is used in this case study. Since the developed model is based on normalized data, it is essential to map input data to its normalized form before the running of the model. The output results of the model should also be changed to its real values for output results to be compared with natural gas processing plant data. Training, optimization and testing are three different subsets of data which are required to perform model development. To develop a model for sulfur content of propane and butane of this natural gas processing plant, a 2.40 GHz processor accompanied by 4G RAM is provided and the model development is programmed in Matlab environment and takes about 5 h to get global optimum results. According to the obtained results it seems that the optimum model result is achieved; however, to guarantee the model reliability, the input variables of test data subset are entered to the developed model and model predictions are validated against plant data. The results of model validation for total sulfur are depicted in Figures 13 and 14. These two figures represent total sulfur of LPG products operating data for output variables versus SVM model prediction results. As can be seen SVM model performed acceptable agreement with the real plant data. However, it can be seen from this figure that some predicted output have much higher deviation from the real plant data. Such deviations could be due to some inherent noise of real plant data that can be alleviated if the learning procedure is equipped with some proper noise filtering routines. No filtering tool is used in this study to expose reliability of model prediction for industrial applications.
Comparison between simulation results and real total sulfur of propane.
Comparison between simulation results and real total sulfur of butane.
To quantify the difference between these two models Average Absolute Deviation Percent (AAD%), as defined by equation (18), is used:(18)where y i , x i , and n represent operating plant data, model predictions and number of operating plant data point used to calculate AAD% respectively. Summary of calculated AAD% for SVM based model prediction of this natural gas processing plant is presented in Table 5.
AAD% values of SVM models for total sulfur of LPG treatment unit.
Also, Table 6 reports accuracy of developed models in terms of MSE, MAE and squared correlation coefficient (R 2) between the operating plant data and SVM prediction results. A SVM based model is optimum if R 2, MAE and MSE are found as close as possible to 1, 0, and 0, respectively.
Statistical parameters of the performance of developed model for sulfur content.
As indicated, SVM model prediction results show acceptable compatibility with the actual plant data. Therefore, SVM could be regarded as a strong tool for prediction output parameter of a natural gas processing plant.
This study demonstrates the applicability of SVM to develop accurate input/output model for total sulfur content of LPG treatment unit of south Pars natural gas processing plant. The plant itself is a very complex one in natural gas industries and the real time data used is a valuable test that allows reliable evaluation of SVM model. As indicated in these two models, SVM model prediction results show more compatibility with the actual plant data. The kernel parameters for developed model are determined and model predictions are compared with real plant data of amine and caustic extractor columns. The numerical values of AAD% calculated for output variables showed a great importance if the predicted data are to be used for monitoring and/or control purposes. This study reveals the applicability and reliability of SVM as a modeling tool in oil and gas industries. Such approaches for oil and gas industries are perfect candidates for applications where the complexity of the data or task demands high computational costs
The authors are grateful to South Pars Gas Complex Co. (SPGC), Bushehr, Iran for their technical assistance and support.
- Adib H., Haghbakhsh R., Saeidi M. (2013) Modeling and optimization of Fischere–Tropsch synthesis in the presence of Co (III)/Al2O3 catalyst using artificial neural networks and genetic algorithm, J. Nat. Gas Sci. Eng. 10, 14–24, https://doi.org/10.1016/j.jngse.2012.09.001. [Google Scholar]
- Adib H., Sabet A., Naderifar A., Adib M., Ebrahimzadeh M. (2015) Evolving a prediction model based on machine learning approach for hydrogen sulfide removal from sour condensate of south Pars natural gas processing plant. J. Nat. Gas Sci. Eng. 27, 74–81, http://doi.org/10.1016/j.jngse.2015.08.012. [Google Scholar]
- Agarwal S., Saradhi V.V., Karnick H. (2008) Kernel-based online machine learning and support vector reduction, Neurocomputing 71, 1230–1237, https://doi.org/10.1016/j.neucom.2007.11.023. [Google Scholar]
- Asil A., Shahsavand A. (2014) Reliable estimation of optimal sulfinol concentration in gas treatment unit via novel stabilized MLP and regularization network, J. Nat. Gas Sci. Eng. 21, 791–804, https://doi.org/10.1016/j.jngse.2014.09.033. [Google Scholar]
- Bashipour F., Rahimi A., Nouri Khorasani S., Naderinik A. (2017) Experimental optimization and modeling of sodium sulfide production from H2S-rich off-gas via response surface methodology and artificial neural network, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 72, 9, https://doi.org/10.2516/ogst/2017004. [CrossRef] [Google Scholar]
- Chiang L.H., Kotanchek M.E., Kordon A.K. (2004) Fault diagnosis based on Fisher discriminant analysis and support vector machines, Comput. Chem. Eng. 28, 1389–1401, https://doi.org/10.1016/j.compchemeng.2003.10.002. [Google Scholar]
- Cortes C., Vapnik V. (1995) Support-vector networks, Machine Learning 20, 273–297. [Google Scholar]
- Curilem M., Acuña G., Cubillos F., Vyhmeister E. (2011) Neural Networks and Support Vector Machine models applied to energy consumption optimization in semiautogeneous grinding, Chem. Eng. Trans. 25, 761–766, https://doi.org/10.3303/CET1125127. [Google Scholar]
- Ding J., Cao Y., Mpofu E., Shi Z. (2012) A hybrid support vector machine and fuzzy reasoning based fault diagnosis and rescue system for stable glutamate fermentation, Chem. Eng. Res. Des., 90, 1197–1207, https://doi.org/10.1016/j.cherd.2012.01.004. [Google Scholar]
- Eslamimanesh A., Gharagheizib F., Illbeigi M., Mohammadi A.H., Fazlali A., Richon D. (2012) Phase equilibrium modeling of clathrate hydrates of methane, carbon dioxide, nitrogen, and hydrogen + water soluble organic promoters using Support Vector Machine algorithm, Fluid. Phase. Equilibria. 316, 34–45, https://doi.org/10.1016/j.proeng.2016.06.455. [Google Scholar]
- Haghbakhsh R., Adib H., Keshavarz P., Koolivand M., Keshtkari S. (2012) Development of an artificial neural network model for the prediction of hydrocarbon density at high-pressure, high-temperature conditions, Thermochim. Acta 551, 124–130, https://doi.org/10.1016/j.tca.2012.10.022. [Google Scholar]
- Kulkarni A., Jayaraman V.K., Kulkarni B.D. (2005) Knowledge incorporated support vector machines to detect faults in Tennessee Eastman Process, Comput. Chem. Eng. 29, 2128–2133, https://doi.org/10.1016/j.compchemeng.2005.06.006. [Google Scholar]
- Magné-Drisch J., Gazarian J., Gonnard S., Schweitzer J.M., Chiche D., Laborie G., Perdu G. (2016) COSWEETTM: a new process to reach very high COS specification on natural gas treatment combined with selective H2S removal, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 71, 40, https://doi.org/10.2516/ogst/2015038. [CrossRef] [Google Scholar]
- Mahdipoor H., Ashkezari A. (2016) Feasibility study of a sulfur recovery unit containing mercaptans in lean acid gas feed, J. Nat. Gas. Sci. Eng. 31, 585–588, https://doi.org/10.1016/j.jngse.2016.03.045. [Google Scholar]
- Mehdizadeh B., Movagharnejad K. (2011) A comparative study between LS-SVM method and semi empirical equations for modeling the solubility of different solutes in supercritical carbon dioxide, Chem. Eng. Res. Des. 89, 2420–2427, https://doi.org/10.1016/j.cherd.2011.03.017. [Google Scholar]
- Moaseri E., Mostaghisi O., Shahsavand A., Bazubandi B., Karimi K., Ahmadi J. (2013) Experimental study and techno-economical evaluation of Khangiran sour natural gas condensate desulfurization process, J. Nat. Gas Sci. Eng. 12, 34–42, https://doi.org/10.1016/j.jngse.2013.02.001. [Google Scholar]
- Moradi N., Adib H., Sabet A., Adhami M., Adib M. (2016) Toward an intelligent approach for H2S content and vapor pressure of sour condensate of south Pars natural gas processing plant, J. Nat. Gas Sci. Eng. 28, 365–371, https://doi.org/10.1016/j.jngse.2015.12.006. [Google Scholar]
- Pelckmans K., Suykens J.A.K., Van Gestel T., De Brabanter D., Lukas L., Hamers B., De Moor B., Vandewalle J. (2002) LS-SVMlab: a Matlab/C toolbox for least squares support vector machines, KU Leuven, Leuven, Belgium. [Google Scholar]
- Raynal L., Augier F., Bazer-Bachi F., Haroun Y., Pereira da Fonte C. (2016) CFD applied to process development in the oil and gas industry – a review, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 71, 42, https://doi.org/10.2516/ogst/2015019. [CrossRef] [Google Scholar]
- Safadoost A., Davoodi M., Mansoori A. (2014) Preventing corrosion and tube failure in sulfur condenser during normal operation, startup, and shutdown of the south Pars gas processing plant (case study), J. Nat. Gas. Sci. Eng. 19, 105–115, https://doi.org/10.1016/j.jngse.2014.05.003. [Google Scholar]
- Santos J., Lobato A., Cunha A., Silva G., Santos L. (2016) Comparison of different processes for preventing deposition of elemental sulfur in natural gas pipelines: A review, J. Nat. Gas. Sci. Eng. 32, 364–372, https://doi.org/10.1016/j.jngse.2016.04.045. [Google Scholar]
- Sharifi A., Omidbakhsh Amiri E. (2017) Effect of the tower type on the gas sweetening process, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 72, 24, https://doi.org/10.2516/ogst/2017018. [CrossRef] [Google Scholar]
- Suykens J.A.K., Gestel T.V., Brabanter J.D., Moor B.D., Vandewalle J. (2002) Least Squares Support Vector Machines, World Scientific, Singapore. [CrossRef] [Google Scholar]
- Terzica J., Nagarajahb C.R., Alamgira M. (2010) Fluid level measurement in dynamic environments using a single ultrasonic sensor and Support Vector Machine (SVM), Sensor Actuator. 161, 278–287, https://doi.org/10.1016/j.sna.2010.05.005. [CrossRef] [Google Scholar]
- Vapnik V.N. (1998) Statistical Learning Theory, Wiley, New York. [Google Scholar]
- Yélamosa I., Escuderob G., Graells M., Puigjaner L. (2009) Performance assessment of a novel fault diagnosis system based on support vector machines, Comput. Chem. Eng. 33, 244–255, https://doi.org/10.1016/j.compchemeng.2008.08.008. [Google Scholar]
- Zaidi S. (2015) Novel application of support vector machines to model the two phase boiling heat transfer coefficient in a vertical tube thermosiphon reboiler, Chem. Eng. Res. Des. 98, 44–58, https://doi.org/10.1016/j.cherd.2015.04.002. [Google Scholar]
- Zanghirati G., Zanni L. (2003) A parallel solver for large quadratic programs in training support vector machines, Parallel Comput. 29, 535–551, https://doi.org/10.1016/S0167-8191(03)00021-8. [Google Scholar]
- Zhao B. (2009) Modeling pressure drop coefficient for cyclone separators: a support vector machine approach,. Chem. Eng. Sci. 64, 4131–4136, https://doi.org/10.1016/j.ces.2009.06.017. [Google Scholar]