1 Introduction

Reducing the use of fossil fuels is one of the ways to prevent atmospheric and lithospheric pollutions. The fossil fuels have caused irreparable damage to all parts of the environment by generation of the greenhouse gas. But it is undeniable the human need for this valuable substance. Therefore, the production of fuel should be somehow that: (1) minimizes damage to the environment; (2) is available; (3) provides the necessary energy; (4) conforms to existing infrastructure; it is always a concern for all of the many politicians and scientists in the world. One of the activities undertaken in this field is the production of oil derivatives using Fischer–Tropsch (FT) process (FTP) (Aitken et al. 2016). Fischer–Tropsch synthesis is a catalytic process that its raw materials produce (H2, CO, CO2) from biosphere (forest waste and coal) (González-Delgado et al. 2015; Lueking and Cole 2017; Domínguez-García et al. 2017; Landau et al. 2014; Weckhuysen 2013). For this reason, the use of the FTP in areas suffering from supply lack of fossil fuels (Europe and South America) is a very good way to produce fuel. FTP produces the oil derivatives with high quality, (without aromatic and devoid of sulfur) (Xinling and Zhen 2009; Zhang et al. 2014) during the reactions (1–3) (Abelló and Montané 2011):

$$\left( {2{{n}} + 1} \right){\text{H}}_{2} + {\text{CO}} \to {\text{C}}_{{n}} {\text{H}}_{{2{{n}} + {{n}}}} + {{n}\text{H}}_{2} {\text{O}}$$
(1)
$$2{n}{\text{H}}_{2} + {n}{\text{CO}} \to {\text{C}}_{{n}} {\text{H}}_{{2{{n}}}} + {n}{\text{H}}_{2} {\text{O}}$$
(2)
$${\text{CO}} + {\text{H}}_{2} {\text{O}} \leftrightarrow {\text{CO}}_{2} + {\text{H}}_{2}$$
(3)

The products of FTP are dependent on catalyst and operational conditions applied to the process (Silva et al. 2008). Traditionally, for optimal production in this process, have been used cobalt and iron catalysts (Noureldin et al. 2014; Liu et al. 2014). When the synthesis gas (H2, CO, CO2) is extracted from coal and waste of forest, the H2/CO ratio is lower. In this case, iron-based catalyst is used to produce more products (Zhang et al. 2010). The iron catalyst provides the additional hydrogen for the FT synthesis by the water–gas shift reaction (WGS), thereby is avoided the accumulation of carbon at the catalyst level, as a result the iron-based catalyst is stably. To improve the performance of the catalyst, promoters are used. Promoters cause the increase of stability, activity, and selectivity of the catalyst with developing active sites on its the surface (Ayodele et al. 2017). The most popular promoters for iron catalyst, can mention Cu and K. Cu facilitates the reduction of iron catalyst, increases the absorption of H2, but heavily suppressed CO adsorption and carburization. On the other hand, K absorbs CO and promotes carburization of the iron catalyst (Li et al. 2015). The presence of both Cu and K promoters improve the performance of the WGS reaction and FTS (Ma et al. 2014; Wan et al. 2008). The supports are as the catalysts structural booster. When the FTS catalyst does not have a structural booster, it suffers from poor flexibility, low activity, and sintering. Among the structural booster, Si and Al can be pointed out, which they can stop sintering and increase the active phase stabilize and physical resistance. Si has a high carbonization capability, thus better absorbs CO. In the study of the two supports effect (SiO2 and Al2O3) on iron based catalyst with presence of potassium, Wu Bao Shan pointed out that SiO2 increased the activity of FTS and WGS and whereas Al2O3 decreased them (Rafati et al. 2015; Suo et al. 2012; Wan et al. 2007). The collection of iron, potassium and copper on silica, are achieved a catalyst which has high activity and a very good structure for the FT synthesis. Fe/Cu/K/SiO2 is an important industrial catalyst which various studies have been conducted at different weight percentages (Bukur et al. 2005; Bukur and Sivaraj 2002; Ma et al. 2016; Vengsarkar et al. 2015; Todic et al. 2016). Bukur researched the effect of reaction conditions on the selectivity of olefin by this catalyst. He reported increase of the pressure and temperature of reaction have a positive effect on the products selectivity (Bukur et al. 1997). In addition, Qing examined the importance of reduction conditions, including temperature, pressure, and the type of gas used in the reduction, on selectivity of products by Fe/Cu/K/SiO2 catalyst (Hao et al. 2007). Many studies indicate how reaction conditions influence on the degree of selectivity of products (Hao et al. 2007; Bukur et al. 1999; Duvenhage and Coville 2002; Ngantsoue-Hoc et al. 2002). However, none of the studies has provided a model so that expresses the significance of these factors. Also, the effect of different parameters on the FT reaction and the interaction between those parameters can be explored. For this purpose, the response surface methodology (RSM) is used. RSM is a set of mathematical relations that can describe the performance of a catalyst under different operating conditions and by means of that the optimum conditions are achieved. The aim of the study is to achieve a selectivity model for the products of the FT process under the iron-based catalyst (Fe–Cu–K/SiO2) in the operating conditions applied in Hao et al. (2007).

2 Experimental

The used catalyst was prepared by mixing the three metals—Fe, Cu, K— and on the SiO2 support through drying-spraying technique.

2.1 Preparation of the catalyst

Fe (NO3)3 and Cu (NO3)2 were mixed in a dish with silica gel solution, then the solution sedimented. Following sedimentation, sodium carbonate solution was added to the sediment for stabilization of pH at 9 ± 0.1. So, the sediment was filtered and washed. To make K/Fe ratio 6/100, K2NO3 solution and distilled water were added to the sediment. The resulting cake was taken out, dried at 250 °C at 2.5 MPa, spray-dried, and crushed into powder. The obtained powder had a regular spherical form and suitable physical resistance for FTS slurry. Calcination temperature was 320 °C for 5 h. Eventually, a catalyst with the weight percentage of 100Fe/5Cu/6K/16SiO2 was made. The reduction conditions for this catalyst were studied at different temperatures. The area of BET, volume, and mean size of the pores of the catalyst at different temperature reduction (TR) are given in Table 1.

Table 1 The properties of the fresh catalyst and the reduced catalyst
Full size table

2.2 Experiment with the catalyst

To perform Fischer–Tropsch process, a stirred tank slurry reactor was used. Almost 20 g of the catalyst was placed inside the reactor to examine the effect of TR on the performance of the catalyst. Afterwards, the catalyst was reduced at the different temperatures (240–280 °C), pressure of 0.4 MPa, and synthesis gas (H2/CO = 0.67) for 13 h. The experiments were conducted under the following conditions: reaction temperature of 250 °C, the pressure of 1.5 MPa, and the specific velocity of 2 NL/g-cat h. To identify the products, the gas chromatography (GC) device was used. Using the data obtained from GC and the mathematical relations (that will have provided), the model of selectivity of the products are determined. The products are mainly methane, light hydrocarbons, gasoline, diesel, and wax. As previously stated, this study focuses on determination a comprehensive model for the selectivity of products.

3 Response surface methodology

To determine a favourable model, the effect of each parameter and their interactions are required. Among all the multivariate models for optimization analytical, the response surface methodology is the best—it is a powerful mathematical and statistical technique for developing, improving, and optimizing of the conditions in multivariate systems. Furthermore, it has the ability to formulate new products and the effect of independent variables and their interaction in the process. In this technique, the general relationship between the response and the input data is depicted in Eq. (4).

$$\psi = g(x_{1} ,x_{2} ,x_{3} , \ldots ,x_{n} ) \pm \varepsilon$$
(4)

where \(\uppsi\) is the dependent variable or response, g is an unknown complex function of the independent variables \(x_{1} ,x_{2} ,x_{3} , \ldots ,x_{n}\), n is the number of input data, and \(\upvarepsilon\) is the fitting error, often assuming it to have a normal distribution with mean zero.

The steps of RSM as an optimization method are as follows:

  1. (1)

    The first step defines the independent variables of major effects on the reaction system via screening experiments, because it is not possible to identify the effects of all parameters.

  2. (2)

    In the second step, the experiments are designed according to tested characteristics such as selection of experimental points and the number of runs and blocks.

  3. (3)

    In the third stage, the independent parameters are made dimensionless. Each effective parameter has a specific dimension. In order to consider a logical relationship among these parameters, they should be made dimensionless. For this purpose, Eq. (5) has been used:

    $${Z} = \frac{{{{{z} - \left[ {{z}_{\text{max} } + {z}_{\text{min} } } \right]} \mathord{\left/ {\vphantom {{{\text z} - \left[ {{\text z}_{\text{max} } + {\text z}_{\text{min} } } \right]} 2}} \right. \kern-0pt} 2}}}{{{{{z}_{\text{max} } - {z}_{\text{min} } } \mathord{\left/ {\vphantom {{{\text z}_{\text{max} } -{\text z}_{\text{min} } } 2}} \right. \kern-0pt} 2}}}$$
    (5)

    Here, Z denotes the independent dimensionless parameter, z represents the experimental data, zmax shows the maximum value of the independent experimental parameter, and zmin indicates the minimum value of the independent experimental parameter.

  4. (4)

    The fourth stage describes the interaction between the different experimental variables. This way, the model used in RSM is the second order model of the response surface g, which can be written as follows:

    $$g = a_{0} + \sum\limits_{i = 1}^{n} {a_{i} x_{i} } + \sum\limits_{i = 1}^{n} {a_{ii} x_{i}^{2} } + \sum\limits_{ipj}^{{}} {a_{ij} x_{i} x_{j} } \pm \varepsilon$$
    (6)

    Here, \(a_{0} \text{,}a_{i} \text{,}a_{ii} \text{,}a_{ij}\) are constant terms, the coefficients of the linear parameters, and the coefficients of quadratic, denotes the coefficients of the interaction parameters respectively.

  5. (5)

    In the fifth stage, the mathematical model is unfamiliar, so one should check whether the model fits the experimental data well. For this purpose, methods such as residual analysis, scaling residual, prediction error sum of squares residual, and testing of the lack of fit can be applied. The comprehensive analytical ability of the model is commonly explained by the coefficient of determination (R2). Using analysis of variance (ANOVA), summarized in Table 3, R2 equation is as follows:

    $$R^{2} = 1 - \frac{{{\text{SS}}_{\text{Error}} }}{{{\text{SS}}_{\text{Total}} }}$$
    (7)

    The closer the value of R2 is to one, the better is the regression. However, adding a variable to the model is the cause of the high value of R2. Hence, a different index is used to ensure the best of fitting, which is named adjusted R2. Adjusted R2 is a function of R2 that regulates the number of illustrative terms in a model. Unlike R2, the adjusted R2 rises as long as the new term improves the model more than would be expected otherwise. The adjusted R2 can be negative, and will always be less than or equal to R2. It is defined in Eq. (8).

    $$R_{\text{adj}}^{2} = 1 - \frac{{\left( {1 - R^{2} } \right)\left( {n - 1} \right)}}{n - p - 1}$$
    (8)

    n is the number of input data, and p is the total number of variables in the model.

  6. (6)

    Draw three-dimensional graphics based on the obtained equations and determine the optimum points.

In this paper the independent variables are TR, TOS, and H2/CO ratio inlet gas and dependent variables are the selectivity of products from FTS that were obtained in the laboratory scale. Finally, mathematical relations will have been determined between them with the RSM.

4 Results and discussion

As mentioned above, the effect of TR, TOS, and H2/CO ratio inlet gas on the amount of selectivity production of FTS was investigated in this study using a mathematical model. The models involve three input variables—TR, TOS, and H2/CO—with five responses-selectivity of methane, light hydrocarbons, gasoline, diesel, and wax. The models were presented by RSM. In Table 2 the power of estimation of the RSM models are compared with experimental values. The regression parameters R2 and R 2adj are provided in Table 3.

Table 2 Comparison between the experimental data were obtained with STSR reactor and the value of selectivity of FTS products were calculated with the RSM on iron-based catalyst
Full size table
Table 3 Regression analysis for selectivity of the products were obtained with Fischer–Tropsch reaction
Full size table

4.1 Selectivity model for CH4

Figure 1 shows the changes in the selectivity of CH4 with TR, TOS, and H2/CO ratio. As be observed, with the concurrent decreased of TR and TOS, the low values of methane selectivity can be obtained at a constant H2/CO ratio. Furthermore, with decrease in H2/CO ratio and TR, the selectivity of CH4 can be lowered at a constant TOS. When the reduction temperature is a constant value, the selectivity of methane lessens with the increase of TOS and the decline of H2/CO ratio. In optimization, minimum methane selectivity is obtained as an aim. For this reason, the reaction conditions are considering as TOS = 404 h, TR = 242 °C, and H2/CO = 0.90 mol/mol. The model of methane selectivity (SC1) is as follows:

$$\begin{aligned} {\text{SC}}_{1} & = - \,0.1592703 + 0.8738342 \times {\text{TR}} + 0.041376 \times {\text{H}}_{ 2} / {\text{CO}} - 0.1145549 \times {\text{TOS}} \\ & \quad + \,0.2976694 \times {\text{TR}}^{2} - 0.2209417 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} - \,0.0401307 \times {\text{TOS}}^{2} \\ & \quad + \,0.2120894 \times {\text{TR}} \times {\text{H}}_{ 2} / {\text{CO}} + 0.1771785 \times {\text{TR}} \times {\text{TOS}} \\ & \quad + \,0.0145077 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}. \\ \end{aligned}$$
(9)
Fig. 1
figure 1

The effects of TR, TOS, and H2/CO ratio Inlet gas on selectivity of CH4: a at H2/CO = 0.75 mol/mol, b at TOS = 300 h, and c at TR = 260 °C

Full size image

4.2 Selectivity model for light hydrocarbons

4.2.1 Selectivity model for light alkanes

Figure 2 represents the changes in the selectivity of light alkanes versus the effective parameters in Fischer–Tropsch reaction, namely TR, TOS, and H2/CO ratio inlet gas. When H2/CO ratio is constant, a better and higher selectivity of light alkanes can be obtained by increasing of the reduction temperature at high levels of TOS. In addition, when the time of stream is constant, increasing of the H2/CO ratio and TR simultaneously cause the increase of light alkanes selectivity. For optimization, the selectivity of light alkane maximizes, which can be obtained by the following conditions establishing: TOS = 422 h, TR = 279 °C, and H2/CO = 0.75 mol/mol. Equation (10) shows the selectivity model of light alkane (SC2–4):

$$\begin{aligned} {\text{SC}}_{2 - 4} & = - \,0.0886742 + 0.8339156 \times {\text{TR}} + 0.1349838 \times {\text{H}}_{2} / {\text{CO}} - 0.0897427 \times {\text{TOS}} \\ & \quad + \,0.2735657 \times {\text{TR}}^{2} - 0.5616624 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} - \,0.1420371 \times {\text{TOS}}^{2} \\ & \quad + \,0.4513428 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} + 0.2748177 \times {\text{TR}} \times {\text{TOS}} \\ & \quad - \,0.2872177 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}. \\ \end{aligned}$$
(10)
Fig. 2
figure 2

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C2–4: a at H2/CO = 0.75 mol/mol, b at TOS = 422 h, and c at TR = 279 °C

Full size image

4.2.2 Selectivity of light alkenes

Figure 3 expresses the changes in the selectivity of light alkenes in relation to TR, TOS, and H2/CO ratio inlet gas. To have a high selectivity of light alkenes, in constant value of H2/CO ratio, when the TR is low must avoided from high values of TOS. Similarly, when the time of stream is constant, for high level of light alkenes selectivity TR and H2/CO ratio value simultaneously must be maximized. The images indicate that if value of the H2/CO ratio and TOS are low, the value of selectivity of light alkenes will diminishes significantly. In optimization, the aim is to maximize of light alkenes selectivity. For this purpose, the conditions of the reaction must be considered as follows: TOS = 437 h, TR = 279 °C, and H2/CO = 0.75 mol/mol. The selectivity model for this product (SC =2–4 ) is obtained with Eq. (11).

$$\begin{aligned} {\text{SC}}_{2 - 4}^{ = } & = - \,0.0566541 + 0.850409 \times {\text{TR}} + 0.0716399 \times {\text{H}}_{2} / {\text{CO}} - 0.1243836 \times {\text{TOS}} \\ & \quad + \,0.2228386 \times {\text{TR}}^{2} - 0.3669702 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} - \,0.1033757 \times {\text{TOS}}^{2} \\ & \quad + \,0.3333077 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} + 0.2337149 \times {\text{TR}} \times {\text{TOS}} \\ & \quad - \,0.1034623 \times {\text{H}}_{ 2} / {\text{CO}} \times {\text{TOS}}; \\ \end{aligned}$$
(11)
Fig. 3
figure 3

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C =2–4 : a at H2/CO = 0.75 mol/mol, b at TOS = 437 h, and c at TR = 279 °C

Full size image

4.3 The selectivity model of gasoline

Figures 4 and 5 represent the changes in the parameters of TR, TOS, and H2/CO ratio in relation to the variations in the gasoline selectivity. The degree of selectivity of paraffin and olefin in gasoline is increased with increases the TR at a constant value of TOS. In addition, at a constant value of TR, both factors TOS and H2/CO ratio must be kept at their maximum value in order to maintain the selectivity of gasoline at a high amount. However, when the H2/CO ratio is constant, the TOS and TR should be elevated to enhance the selectivity of paraffin and olefin present in gasoline. Eventually, in order to obtain the optimal conditions of paraffin and olefin present in gasoline respectively, the following values have been obtained: TOS = 491 h, TR = 244 °C, H2/CO = 0.77 mol/mol, TOS = 392 h, TR = 246 °C, and H2/CO = 0.78 mol/mol. Equations (12) and (13) represent the selectivity of paraffin (SC5–11) and olefin (SC =5–11 ) in gasoline respectively (Table 4).

$$\begin{aligned} {\text{SC}}_{5 - 11} & = - \,0.5933982 + 0.9934692 \times {\text{TR}} - 0.084869 \times {\text{H}}_{2} / {\text{CO}} - 0.0489266 \times {\text{TOS}} \\ & \quad + \,0.4368277 \times {\text{TR}}^{2} + 0.4822179 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} + \,0.340322 \times {\text{TOS}}^{2} \\ & \quad - \,0.1549174 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} - 0.002825 \times {\text{TR}} \times {\text{TOS}} \\ & \quad + \,0.4826408 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}; \\ \end{aligned}$$
(12)
$$\begin{aligned} {\text{SC}}_{5 - 11}^{ = } & = - \,0.4274918 + 1.116938 \times {\text{TR}} - 0.5055766 \times {\text{H}}_{2} / {\text{CO}} - 0.2855771 \times {\text{TOS}} \\ & \quad + \,0.242604 \times {\text{TR}}^{2} + 0.9032077 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} + \,0.385286 \times {\text{TOS}}^{2} \\ & \quad - \,0.4685129 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} - 0.0037375 \times {\text{TR}} \times {\text{TOS}} \\ & \quad + \,0.9002605 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}. \\ \end{aligned}$$
(13)
Fig. 4
figure 4

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C5–11: a at H2/CO = 0.77 mol/mol, b at TOS = 491 h, and c at TR = 244 °C

Full size image
Fig. 5
figure 5

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C =5–11 : a at H2/CO = 0.78 mol/mol, b at TOS = 392 h, and c at TR = 246 °C

Full size image
Table 4 The optimization results for the products obtained from Fischer–Tropsch reaction
Full size table

4.4 The selectivity model of diesel

Figures 6 and 7 demonstrate the variations in the selectivity of diesel against TR, TOS, and H2/CO ratio inlet gas. At a constant value of H2/CO ratio, with the increase of TR and decrease of TOS, maximum selectivity of diesel can be expected. Furthermore, if value of H2/CO ratio and TR are brought to their maximum level, the value of selectivity of this product will be high. When the TR is constant, diesel selectivity would decline significantly if both TOS and H2/CO values are concurrently low. The optimization conditions for this product are obtained by considering the following: TOS = 437 h, TR = 251 °C, H2/CO = 0.77 mol/mol, TOS = 342 h, TR = 266 °C, and H2/CO = 0.75 mol/mol. The following equation has been developed for selectivity of diesel.

$$\begin{aligned} {\text{SC}}_{12 - 18} & = - \,0.9610807 + 0.0798878 \times {\text{TR}} + 0.5557329 \times {\text{H}}_{2} / {\text{CO}} + 0.0964867 \times {\text{TOS}} \\ & \quad + \,1.265648 \times {\text{TR}}^{2} - 1.09784 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} + \,0.2877426 \times {\text{TOS}}^{2} \\ & \quad + \,1.565955 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} + 0.0130835 \times {\text{TR}} \times {\text{TOS}} \\ & \quad - \,0.8732221 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}; \\ \end{aligned}$$
(14)
$$\begin{aligned} {\text{SC}}_{12 - 18}^{ = } & = - \,0.9581628 + 0.2581191 \times {\text{TR}} + 0.4819162 \times {\text{H}}_{2} / {\text{CO}} + 0.5806624 \times {\text{TOS}} \\ & \quad + \,0.958352 \times {\text{TR}}^{2} - 0.8436385 \times \left( {{\text{H}}_{2} / {\text{CO}}} \right)^{2} + \,0.1763373 \times {\text{TOS}}^{2} \\ & \quad + \,1.108191 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} - 0.179325 \times {\text{TR}} \times {\text{TOS}} \\ & \quad - \,1.456383 \times {\text{H}}_{ 2} / {\text{CO}} \times {\text{TOS}}. \\ \end{aligned}$$
(15)
Fig. 6
figure 6

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C12–18: a at H2/CO = 0.75 mol/mol, b at TOS = 300 h, and c TR = 260 °C

Full size image
Fig. 7
figure 7

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of C =12–18 : a at H2/CO = 0.75 mol/mol, b at TOS = 300 h, and c TR = 260 °C

Full size image

4.5 The selectivity model of wax

Figure 8 shows the wax selectivity changes due to the effects of TR, TOS, and H2/CO ratio inlet gas. Changes in H2/CO ratio and TOS against TR almost are the same. When the value of TR is low, increase of H2/CO ratio and TOS will increase the selectivity of wax. In high quantities of H2/CO, wax selectivity is improved by enhance the amount of TOS. The optimization objective is to obtain a high selectivity of wax, which can be attained with TOS = 283 h, TR = 280 °C and H2/CO = 0.81 mol/mol. The final mathematical model for selectivity of wax (\({\text{SC}}_{19}^{ + }\)) is Eq. (16):

$$\begin{aligned} {\text{SC}}_{{19^{ + } }} & = 0.4114937 - 0.9569705 \times {\text{TR}} + 0.1312774 \times {\text{H}}_{2} / {\text{CO}} + 0.1279498 \times {\text{TOS}} \\ & \quad - \,0.4225025 \times {\text{TR}}^{2} + 0.1391603 \times {\text{H}}_{2} / {\text{CO}}^{2} - 0.0576969 \times {\text{TOS}}^{2} \\ & \quad - \,0.1950588 \times {\text{TR}} \times {\text{H}}_{2} / {\text{CO}} - 0.1001866 \times {\text{TR}} \times {\text{TOS}} \\ & \quad + \,0.0771063 \times {\text{H}}_{2} / {\text{CO}} \times {\text{TOS}}. \\ \end{aligned}$$
(16)
Fig. 8
figure 8

The effects of TR, TOS, and H2/CO ratio inlet gas on selectivity of \({\text{C}}_{19}^{ + }\): a at H2/CO = 0.75 mol/mol, b at TOS = 300 h, and c TR = 260 °C

Full size image

5 Conclusions

Synthesis gas (H2, CO, CO2) obtain from forest wastes in high steam pressure and the free sulfur oil derivation are produced from them with the Fischer–Tropsch process. The determination of a selectivity model of every product is the great help for economic purposes, such that the best possible conditions for the economical production of favourite products are obtained with using these models. In this study, through the response surface methodology for the selectivity of five products—methane, light hydrocarbons, gasoline, diesel, and wax—in the Fischer–Tropsch reaction, eight selectivity models were stated in terms of the process conditions of the catalyst: TR, TOS, and H2/CO ratio inlet gas. These models indicated the TR change has a significant effect on the selectivity of products. So that, the TR decline causes the considerable decrease of methane and wax selectivity. To the production increase of light hydrocarbons, the amount of TR should have been increased. If objective is the minimization of all products selectivity and the maximization of selectivity of light olefin, the optimal operating conditions as TR = 258 °C, TOS = 375 h and H2/CO = 0.77 mol/mol will be considered.