Next Article in Journal
Special Issue on “Application of Big Data Analysis and Advanced Analytics in Sustainable Production Process”
Previous Article in Journal
High-Property Anode Catalyst Compositing Co-Based Perovskite and NiFe-Layered Double Hydroxide for Alkaline Seawater Splitting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Domain of Parameters in the Stochastic Reconstruction Method

Key Laboratory of Smart Manufacturing in Energy Chemical Process, Ministry of Education, East China University of Science and Technology, 130 Meilong Rd, Shanghai 200237, China
*
Authors to whom correspondence should be addressed.
Processes 2022, 10(4), 669; https://doi.org/10.3390/pr10040669
Submission received: 8 March 2022 / Revised: 23 March 2022 / Accepted: 28 March 2022 / Published: 30 March 2022
(This article belongs to the Section Chemical Processes and Systems)

Abstract

:
In the stochastic reconstruction of petroleum fractions, the domain configuration for parameters is of great significance to the performance of the model in terms of accuracy and convergence. Based on the stochastic reconstruction model of gas oils, this work investigates the influence of different domains on the simulating accuracy, the convergence rate and the ability to predict detailed composition. In this paper, the parameters to be optimized in the model include histogram distributions and gamma distributions, which are used to represent the structural attributes. In each histogram distribution, the parameters are real numbers between zero and one that should increase progressively. The mean and shape parameters are used as the tuning parameters for gamma distribution. As for the mean, the lower and upper boundaries are configured as fixed values that depend on the predefined values of the structural attributes. For the shape parameter, multiple cases are considered for its upper boundary, which is set to 5, 10, 20, 50 and 100, respectively. The lower boundary is set to 0 or 1, depending on the structural attribute. The results indicate that the simulating accuracy of the stochastic reconstruction model improves with the increase of the upper boundary for the shape parameter. When the upper boundary of each shape parameter is 20, the simulating accuracy is better than that when the upper boundary is 5 and 10, and close to that when the upper boundary is 50 and 100. The detailed group-type analysis is involved in the model to predict the distributions of mass fraction in chemical families against the carbon number. The results show that the predicted distributions in normal paraffins, isoparaffins, naphthenes, and aromatics are closest to the experimental distributions when the upper boundary is 20 and 50. When the upper boundary is 5 and 10, the predicted distributions in normal paraffins and isoparaffins are lower and broader than the experimental distributions. On the contrary, the predicted distributions in normal paraffins and isoparaffins are higher when the upper boundary is 100. By studying the effects of different domains on the results of the stochastic reconstruction model, a better molecular characterization of petroleum fractions can be achieved.

1. Introduction

Petroleum fractions consist of tens of thousands of molecular species that contain carbon, hydrogen, sulfur, nitrogen, oxygen and metals. Due to the diverse origins of crude oils and refining technologies, the compositions of material streams in the refining processes vary widely. In recent years, with increasingly strict environmental regulations for fossil fuels as well as more and more heavy crude oils, it is urgent to improve refining technologies [1,2]. As the foundation for the study of refining technologies, the characterization of petroleum fractions and the construction of conversion mechanisms at the molecular level have become the most important issues in refineries [3,4].
Due to the limitation of analytical chemistry technologies, computer-aided molecular reconstruction is currently the most popular method to reflect the complicated compositions of petroleum with limited properties [4]. Since the early 1990s when the first molecular reconstruction method was reported, a lot of molecular reconstruction methods have been proposed and applied in the molecular modeling of refining processes [5,6,7,8], including the stochastic reconstruction method (SR), structural-oriented lumping (SOL), molecular type homologous series (MTHS) method and SR-reconstruction by entropy maximization (SR-REM). Among them, the SR method is the most popular method in the characterization of heavy petroleum fractions.
Klein et al. [9,10] first introduced the SR method to characterize heavy residue feedstocks. In the SR method, molecules are treated as the collection of structural attributes (SA), and a probability density function (PDF) is imposed on each structural attribute. The type of PDF varies for different structural attributes, such as the histogram distribution for the determination of molecular families and the gamma distribution for ring numbers and the length of chains. Monte Carlo sampling with a quadrature method is applied to generate an equimolar set of molecules from PDFs. The parameters for each PDF were adjusted in an optimization loop for simulated annealing or genetic algorithm to make the bulk properties of generated mixtures close to those of the actual samples.
Many studies that use the SR method have been reported. Petti et al. examined [11] the usage of CPU resources in the SR method and suggested that a sample size of 10,000 molecules could balance the simulation accuracy and computation expense. Zhang et al. [12] extended a novel SR model to heavy vacuum residue fractions. The residue molecules were treated as a combination of approximately 600 building substructures. Deniz et al. [13] introduced a new structure parameter set for detailed ring and chain configurations into the SR method to improve the method performance in heavy petroleum fractions. Moreover, Deniz et al. investigated [14] the effects of methods for estimating the boiling point temperature and density of pure compounds on the simulating accuracy of the SR method. It is observed that the SR method has the highest accuracy with the group contribution method by Gani [15,16,17,18] and the Yen-woods equations [19]. Haktanlr et al. [20] proposed a novel SR method based on a custom predefined molecular library. This novel method focuses on characterizing the petroleum fractions with exhausted molecular species. Meanwhile, a sieving mechanism is introduced to make sure that the generated molecules are reasonable in structure. Glazov et al. studied [21] the relationships between different PDFs of structural attributes and bulk properties.
The general expression of the SR method is shown as below:
E = m i n   f [ x , ξ ] x l o w x x u p p e r x
where x stands for the parameters to be optimized, ξ are series of uniform random numbers between 0 and 1. xlow and xupper are the lower and upper boundaries of x . The function f stands for the gap between the experimental values of bulk properties and the predicted values. The optimal parameters that make the predicted values of bulk properties close enough to the experimental values are obtained by effective algorithm optimization.
In the SR method, the input parameters x are the parameters of PDFs. The structural attributes generally contain the type of molecule, the number of naphthenic rings and aromatic rings, the length of paraffin chains and sidechains, the type of heteroatom-containing molecule, etc. The types of PDF used to represent them are histogram distribution and gamma distribution. The parameters in histogram distribution are real numbers between 0 and 1. If more than one parameter exists in a histogram distribution, they should increase progressively. Given the total number of parameters in the model and the simulating accuracy, the two-parameter gamma distribution is generally adopted in the SR method. The two parameters are the shape parameter (SP) and scale parameter, or the shape parameter and mean. To characterize the complicated and diverse compositions in petroleum fractions, it is crucial to reasonably configure the lower and upper boundaries of the parameters, which is the domain of parameters in the SR method. However, as far as the authors know, no study has been reported on this issue. To this point, this paper aims to study the determination of domain and its effects on the performance of the SR method.
Section 2 illustrates the SR model based on gas oils. The bulk properties by experiments, the setting of structural attributes and the building diagram are provided in Section 2.1. Section 2.2 provides the configuration of the lower and upper boundaries for the parameters in the histogram distributions and gamma distributions. Section 3 provides the results and discussions. Finally, the conclusion is given in Section 4.

2. Methodology

Figure 1 shows the framework of the SR model adopted in this paper. The model starts with initializing the input parameters, which are the tuning parameters in PDFs. Input parameters are generated based on their lower and upper boundaries. For parameters in histogram distributions, the monotonicity constraint, that parameters in the same histogram distribution should increase progressively, should be satisfied. With the input parameters, Monte Carlo sampling is used to generate pseudo molecules, the properties of which are estimated by group contribution methods and empirical correlations. The bulk properties of pseudo mixtures are calculated by mixing rules. The objective function is configured to evaluate the input parameters by counting the gap between the experimental values and the estimated values of the bulk properties. If the objective function value meets the stop criterion, the optimal parameters and pseudo mixture are obtained. Otherwise, new input parameters are generated in the optimization algorithm and evaluated until all constraints are satisfied.
The setting of lower and upper boundaries for input parameters determines the values of the input parameters in the SR model. Multiple cases of boundaries for input parameters are tested and compared based on an SR model of gas oils.

2.1. The SR Model of Gas Oils

The two gas oil samples used in this paper are from the work by Aleksandar [22], including a gas oil sample and a vacuum gas oil sample. As listed in Table 1, the bulk properties adopted contain the elemental analysis, density, H/C ratio, PINA analysis, aromatic sulfur content and simulated distillation.
Based on the molecular characteristics of two gas oil samples, the probability distributions containing the histogram distribution and gamma distribution are adopted to represent structural attributes. As shown in Table 2, the structural attributes represented by the histogram distribution are the type of molecule, the acceptance of a naphthenic ring on an aromatic ring, the acceptance of a branch on a paraffin chain, and the type of sulfur-containing compound. The gamma distribution is adopted for other structural attributes, such as the number of naphthenic rings, the number of aromatic rings, the length of paraffin chain, the length of sidechain, and the length of sulfur-containing chain.
Based on the structural attributes in Table 2, a building diagram for the construction of pseudo molecules is designed in Figure 2. In this building diagram, the first step is to determine the type of molecule. Paraffin, naphthene, aromatic and sulfur-containing compounds are considered in this work. If the type of molecule is paraffin, the next step is to determine whether the paraffin is branched. Then the length of paraffin is determined. If the type of molecule is naphthene, the number of naphthenic rings and the length of sidechain should be determined, in this order. When the type of molecule is aromatic, the number of aromatic rings should be determined first. If the number of aromatic rings is one or two, it should be determined whether a naphthenic ring is accepted. Lastly, the length of sidechain on aromatic rings is determined. There are three types of sulfur-containing compounds, including benzothiophene, dibenzothiophene, and sulfur-containing paraffin. Once the molecule is benzothiophene or dibenzothiophene, the length of sidechain should be determined. Otherwise, the length of sulfur-containing chain is determined. In this work, the predefined range for the number of naphthenic rings in a core is one to three, and for the number of aromatic rings is one to four. The number of sidechains on a core is no more than one. The ranges for lengths of paraffin chains, sidechains and sulfur-containing chains are determined based on the distillation curve.

2.2. The Setting of Domains

For the parameters in the histogram distributions, their lower and upper boundaries are set to 0 and 1, respectively. According to Table 2, the parameters in the histogram distributions are x1, x2, x3, x4, x5, x6 and x7. Because x1, x2, x3 and x6, x7 are in a single histogram distribution, they should meet the constraint that parameters in the same histogram distribution should increase in order, as shown in Equation (2).
x 1 < x 2 < x 3 x 6 < x 7 x 1 , x 2 , , x 7 [ 0 , 1 ]
The two-parameter gamma distribution is shown in Equations (3)–(5).
p ( x ) = ( x r min ) α 1 e ( x r min ) β β α Γ ( α )
Γ ( x ) = 0 t x 1 e t d t
β = m r min α
where α , β and r min are the shape parameter, scale parameter and site parameter in the gamma distribution, respectively. In this work, the shape parameter α and mean m are adopted as the tuning parameters for each gamma distribution. The relationship between α , β , r min and m is shown in Equation (5). r min in each gamma distribution is a fixed value in this work, and depends on the predefined values of the corresponding structural attribute.
Before setting the boundaries of parameters in gamma distributions, the range of predefined values for structural attributes should be clarified. The range of predefined values for the number of naphthenic rings and aromatic rings are 1–3 and 1–4 in both gas oils. According to the distillation curves in Table 1, the calculated ranges for the length of sidechains are 1–35 for gas oil A and 1–23 for gas oil B. The ranges for the length of paraffin chains are 8–41 and 8–30 for the two gas oil samples, respectively. Lastly, the ranges for the length of sulfur-containing chains are 6–40 for gas oil A and 6–29 for gas oil B. The r min in each gamma distribution is equal to the minimum predefined value of the corresponding structural attribute minus 1.
In the gamma distributions, multiple cases for the boundaries of parameters are considered in Table 3 and Table 4. There are five structural attributes that are represented by gamma distribution. The parameters in the model are from x8 to x17. x8, x10, x12, x14 and x16 are the shape parameters of the gamma distributions for the number of naphthenic rings, the number of aromatic rings, the length of sidechains, the length of paraffin chains and the length of sulfur-containing chains, respectively. x9, x11, x13, x15 and x17 are the mean of the gamma distributions. For both gas oils, the lower boundaries (LB) of x8, x10 are 0, and the lower boundaries of x12, x14 and x16 are 1. The upper boundaries (UB) of x8, x10, x12, x14, and x16 are given as 5, 10, 20, 50 and 100. As for the mean in each gamma distribution, the lower boundary is set to the minimum predefined value of the corresponding structural attribute minus 1, and the upper boundary is equal to the maximum predefined value of the corresponding structural attribute plus 1. For simplicity of presentation, the cases where the upper boundaries of the shape parameters are 5, 10, 20, 50 and 100 are marked below as case 1, case 2, case 3, case 4 and case 5, respectively.

3. Results and Discussions

3.1. Comparison of Simulating Performance in Different Cases

Table 5 and Table 6 show the simulation results with the smallest relative error for two gas oils under different domains of the SR model. A fixed number of iterations that is specified as 2000 is adopted as the stop criterion of the SR model. It is observed that the relative errors are very small in all cases ranging from 0.03–0.06. This means that in the case of different domains, the SR model is able to generate pseudo mixtures, the bulk properties of which are very close to the experimental values (EV).
Figure 3 and Figure 4 compare the variations of the average objective function value from iterations 1–400 with different domains in the stochastic reconstruction of gas oil A. The average objective function value is calculated by 20 runs for each case. In Figure 3a, it is observed that at the early stage from iterations 1–200, the average objective function values in cases 1 and 2 are significantly lower than those in cases 3, 4 and 5. This means that in a smaller domain the SR model is able to find parameters that can make the objective function value smaller at the beginning of the simulation. The variations in the average objective function value maintain the same trend for all cases between iterations 200–400. With the number of iterations exceeding 350, the average objective function value in case 1 starts to be larger than that in case 2. As the number of iterations increases from 400 to 1200, as shown in Figure 3c,d, the average objective function value in case 1 gradually exceeds the average objective function values in all other cases. Compared with the early stage of simulation, gaps of average objective function values among all cases are smallest when the number of iterations is around 1100. As the iteration continues, the average objective function value in case 2 gradually exceeds those in cases 3, 4 and 5. There are approximate rates in the decrease of average objective function values in case 3, 4 and 5. At the end of the iteration, it can be observed that average objective function values in case 3, 4 and 5 are very close, which are significantly smaller than those in case 1 and 2. This indicates that even though the average objective function values decrease faster in smaller domains at the beginning of the simulation, the SR model can find parameters that make the objective function value smaller with adequate iterations in larger domains. Meanwhile, for case 3, the decreasing rate of the average objective function value is larger than that in cases 4 and 5 at the early stage of simulation, and the final average objective function value is very close to that in cases 4 and 5, as shown in Figure 3f. This indicates that case 3, where the upper boundaries of shape parameters are 20, can not only make the SR model effectively find promising parameters, but also make sure the targeted objective function value is achieved.

3.2. Prediction of Compositions by the SR Model in Different Domains

The predicted distributions of mass fraction against carbon number in normal paraffins, isoparaffins, naphthenes and aromatics are compared with experimental values in order to investigate the performance of the SR model on the prediction of detailed compositions in different domains. The representative predictions in each case are shown in Figure 4. It is observed that there are no cases where the predicted distributions in normal paraffins and isoparaffins fit well with the experimental distributions. The shapes of the predicted distributions in normal paraffins and isoparaffins are broader and lower than in the experimental distributions. For naphthenes and aromatics, the predicted distributions are very close to the experimental distributions in cases 1, 2 and 3. As the upper boundaries of shape parameters increase, in cases 4 and 5, the shapes of the predicted distributions become higher. As a result, the upper boundaries of shape parameters in the gamma distributions should be less than 50 in order to fit well with the distributions of mass fraction in naphthenes and aromatics.
As shown in Table 5 and Table 6 and in Figure 4, even though the estimated bulk properties in all cases are close to the experimental values, there are large deviations between the predicted distributions of mass fractions and the experimental values in all cases. In order to reduce the deviations between the predicted distributions and experimental distributions in all chemical families, the PINA values are replaced by the more the detailed experimental data from detailed group-type analysis.
The detailed group-type analysis data are mass fractions of homologous series in gas oil samples, which contain normal paraffins (P), isoparaffins (I), mononaphthenes (MN), dinaphthenes (DN), monoaromatics (MA), naphthenoaromatics (NA), diaromatics (DA), naphthenodiaromatics (NDA), triaromatics (TrA), naphthenotriaromatics (NTrA) and tetraaromatics (TeA). Based on detailed group-type analysis and other bulk properties, the stochastic reconstruction of gas oil A is implemented under different domains. As shown in Table 7, the predicted values of the bulk properties are also close to the experimental values in all cases, and the sums of relative errors are between 0.11 and 0.14.
More accurate predictions on the distributions of mass fractions are observed with detailed group-type analysis. As shown in Figure 5, the domain range has a large influence on the predicted distributions. In normal paraffins and isoparaffins, the predicted distributions in cases 1 and 2, where the upper boundaries of the shape parameters are 5 and 10, are still lower and broader compared to the experimental values. In cases 3 and 4, the predicted distributions become higher than those in cases 1 and 2, and agree best with the experimental values. As the upper boundaries of the shape parameters increase, the predicted curves continue to be higher, and a significant deviation is observed in case 5. As for naphthenes and aromatics, similar predicted curves are shown in all cases. Small deviations between the prediction and experiment are observed. The predicted mass fractions of naphthenes are higher than the experimental values when the carbon number is small, and lower than the experimental values as the carbon number increases. In conclusion, the detailed group-type analysis can obviously improve the ability of the SR model to predict distributions of mass fraction in chemical families. Cases 3 and 4, where the upper boundaries of the shape parameters are 20 and 50, are the best cases for the prediction of detailed compositions of gas oil A.
The variation in predicted compositions in different domains is because of the mechanism of the SR method and the characteristic of gamma distribution. In the SR method, there are multiple combinations for the values of the parameters in the domain that make the estimated values of properties close to the experimental values. The gamma distribution curve tends to be lower and broader when the shape parameter is small, and tends to be higher and steeper as the shape parameter increases. As a result, with the experimental curves in this work, even though the target objective function value can be achieved in the small domain, there are large gaps between the predicted curves and the experimental curves. When the value of the shape parameter is too large, for example, a value between 50 and 100 in this paper, an inappropriate curve can also be generated. Therefore, the appropriate prediction of detailed compositions depends on the reasonable configuration of domain.

4. Conclusions

Based on gas oil samples, this work studies the performance of the SR model in different domains of parameters. In the setting of domains, multiple cases for the upper boundaries of shape parameters in gamma distributions are considered. The results show that in each case, the SR model can generate pseudo mixture, the bulk properties of which agree well with the experimental values. The bulk properties adopted are elemental analysis, density, H/C ratio, aromatic sulfur content, PINA analysis and simulated distillation. By investigating the variations in average objective function value in each case, it is observed that even though the decrease of average objective function value in smaller domains is faster at the beginning of the simulation, the SR model can find parameters in wider domains that make the average objective function value smaller. Case 3 is the best case where the simulating accuracy and convergence performance are maintained compared with other cases. The prediction of the SR model on compositions of gas oil sample is not acceptable, despite the estimated values of bulk properties being very close to the experimental values. Large derivations are observed in all cases.
The detailed group-type analysis contributes to more accurate predictions on compositions. The predicted distributions in normal paraffins and isoparaffins become higher as the upper boundaries of shape parameters increase. In cases 3 and 4, the predicted distributions fit best with the experimental distributions, compared to the other cases. As for naphthenes and aromatics, the predicted distributions agree well with experimental values in all cases. Because case 3 performed best in maintaining simulating accuracy and convergence rate, the best upper boundary of shape parameters is 20, compared to other values in this work.
The domain of parameters in probability functions has a great influence on the performance of the SR method. Due to the differences in bulk properties and compositions of petroleum fractions, reasonably setting the domain of parameters is indispensable for the accurate reconstruction of petroleum fractions.

Author Contributions

Conceptualization, G.Z.; methodology, G.Z.; software, G.Z.; validation, G.Z.; formal analysis, G.Z.; investigation, G.Z.; resources, G.Z.; data curation, G.Z.; writing—original draft preparation, G.Z.; writing—review and editing, M.Y.; visualization, G.Z.; supervision, F.Q.; project administration, G.Z.; funding acquisition, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program of China (2021YFB1714300), National Natural Science Fund for Distinguished Young Scholars (61725301), International (Regional) Cooperation and Exchange Project (61720106008) and National Natural Science Foundation of China (61873093). And The APC was funded by National Natural Science Foundation of China (61873093).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PDFProbability density function
PINANormal paraffins, isoparaffins, naphthenes, and aromatics
SAStructural attribute
SPShape parameter
UBUpper boundary
LBLower boundary
MNMononaphthenes
DNDinaphthenes
MAMonoaromatics
NANaphthenoaromatics
DADiaromatics
NDANaphthenodiaromatics
TrATriaromatics
NTrANaphthenotriaromatics
TeATetraaromatics
EVExperimental value

References

  1. Qian, F.; Zhong, W.; Du, W. Fundamental Theories and Key Technologies for Smart and Optimal Manufacturing in the Process Industry. Engineering 2017, 3, 154–160. [Google Scholar] [CrossRef]
  2. Shen, F.; Wang, M.; Huang, L.; Qian, F. Exergy analysis and multi-objective optimisation for energy system: A case study of a separation process in ethylene manufacturing. J. Ind. Eng. Chem. 2021, 93, 394–406. [Google Scholar] [CrossRef]
  3. Stratiev, D.; Shishkova, I.; Tankov, I.; Pavlova, A. Challenges in characterization of residual oils. A Review. J. Pet. Sci. Eng. 2019, 178, 227–250. [Google Scholar] [CrossRef]
  4. Ren, Y.; Liao, Z.; Sun, J.; Jiang, B.; Wang, J.; Yang, Y.; Wu, Q. Molecular reconstruction: Recent progress toward composition modeling of petroleum fractions. Chem. Eng. J. 2019, 357, 761–775. [Google Scholar] [CrossRef]
  5. Hudebine, D.; Verstraete, J.J. Molecular reconstruction of LCO gasoils from overall petroleum analyses. Chem. Eng. Sci. 2004, 59, 4755–4763. [Google Scholar] [CrossRef]
  6. Van Geem, K.M.; Hudebine, D.; Reyniers, M.F.; Wahl, F.; Verstraete, J.J.; Marin, G.B. Molecular reconstruction of naphtha steam cracking feedstocks based on commercial indices. Comput. Chem. Eng. 2007, 31, 1020–1034. [Google Scholar] [CrossRef]
  7. Quann, R.J.; Jaffe, S.B. Structure-oriented lumping: Describing the chemistry of complex hydrocarbon mixtures. Ind. Eng. Chem. Res. 1992, 31, 2483–2497. [Google Scholar] [CrossRef]
  8. Wu, Y. Molecular Management for Refining Operations. Ph.D.Thesis, The University of Manchester, Manchester, UK, 2010. [Google Scholar]
  9. Neurock, M.; Libanati, C.; Nigam, A.; Klein, M.T. Monte carlo simulation of complex reaction systems: Molecular structure and reactivity in modelling heavy oils. Chem. Eng. Sci. 1990, 45, 2083–2088. [Google Scholar] [CrossRef]
  10. Neurock, M.; Nigam, A.; Trauth, D.; Klein, M.T. Molecular representation of complex hydrocarbon feedstocks through efficient characterization and stochastic algorithms. Chem. Eng. Sci. 1994, 49, 4153–4177. [Google Scholar] [CrossRef]
  11. Petti, T.F.; Trauth, D.M.; Stark, S.M.; Neurock, M.; Yasar, M.; Klein, M.T. CPU Issues in the Representation of the Molecular Structure of Petroleum Resid through Characterization, Reaction, and Monte Carlo Modeling. Energy Fuels 1994, 8, 570–575. [Google Scholar] [CrossRef]
  12. Zhang, L.; Hou, Z.; Horton, S.R.; Klein, M.T.; Shi, Q.; Zhao, S.; Xu, C. Molecular representation of petroleum vacuum resid. Energy Fuels 2014, 28, 1736–1749. [Google Scholar] [CrossRef]
  13. Deniz, C.U.; Yasar, M.; Klein, M.T. A New Extended Structural Parameter Set for Stochastic Molecular Reconstruction: Application to Asphaltenes. Energy Fuels 2017, 31, 7919–7931. [Google Scholar] [CrossRef]
  14. Deniz, C.U.; Yasar, S.H.O.; Yasar, M.; Klein, M.T. Effect of Boiling Point and Density Prediction Methods on Stochastic Reconstruction. Energy Fuels 2018, 32, 3344–3355. [Google Scholar] [CrossRef]
  15. Gani, R.; Nielsen, B.; Fredenslund, A. A group contribution approach to computer-aided molecular design. AIChE J. 1991, 37, 1318–1332. [Google Scholar] [CrossRef]
  16. Constantinou, L.; Gani, R.; O’Connell, J.P. Estimation of the acentric factor and the liquid molar volume at 298 K using a new group contribution method. Fluid Phase Equilibria 1995, 103, 11–22. [Google Scholar] [CrossRef]
  17. Marrero, J.; Gani, R. Group-contribution based estimation of pure component properties. Fluid Phase Equilibria 2001, 183–184, 183–208. [Google Scholar] [CrossRef]
  18. Hukkerikar, A.S.; Sarup, B.; Ten Kate, A.; Abildskov, J.; Sin, G.; Gani, R. Group-contribution+ (GC+) based estimation of properties of pure components: Improved property estimation and uncertainty analysis. Fluid Phase Equilibria 2012, 321, 25–43. [Google Scholar] [CrossRef]
  19. Yen, L.C.; Woods, S.S. A generalized equation for computer calculation of liquid densities. AIChE J. 1966, 12, 95–99. [Google Scholar] [CrossRef]
  20. Haktanır, M.; Karahan, S.; Yaşar, M. Structurally explicit composition model of petroleum vacuum residue. Fuel 2021, 300, 120977. [Google Scholar] [CrossRef]
  21. Glazov, N.; Dik, P.; Zagoruiko, A. Effect of experimental data accuracy on stochastic reconstruction of complex hydrocarbon mixture. Catal. Today 2021, 378, 202–210. [Google Scholar] [CrossRef]
  22. Bojkovic, A.; Dijkmans, T.; Dao Thi, H.; Djokic, M.; van Geem, K.M. Molecular Reconstruction of Hydrocarbons and Sulfur-Containing Compounds in Atmospheric and Vacuum Gas Oils. Energy Fuels 2021, 35, 5777–5788. [Google Scholar] [CrossRef]
Figure 1. The framework of the SR model.
Figure 1. The framework of the SR model.
Processes 10 00669 g001
Figure 2. The building diagram in the SR model.
Figure 2. The building diagram in the SR model.
Processes 10 00669 g002
Figure 3. Variations in the average objective function value of 20 runs from iterations (a) 1–200, (b) 201–400, (c) 401–800, (d) 801–1200, (e) 1201–1600, (f) 1601–2000.
Figure 3. Variations in the average objective function value of 20 runs from iterations (a) 1–200, (b) 201–400, (c) 401–800, (d) 801–1200, (e) 1201–1600, (f) 1601–2000.
Processes 10 00669 g003
Figure 4. The distributions of mass fraction in (a) normal paraffins, (b) isoparaffins, (c) naphthenes, (d) aromatics for gas oil A by experiment and predictions.
Figure 4. The distributions of mass fraction in (a) normal paraffins, (b) isoparaffins, (c) naphthenes, (d) aromatics for gas oil A by experiment and predictions.
Processes 10 00669 g004
Figure 5. The distributions of mass fraction in (a) normal paraffins, (b) isoparaffins, (c) naphthenes, (d) aromatics for gas oil A by experiment and predictions.
Figure 5. The distributions of mass fraction in (a) normal paraffins, (b) isoparaffins, (c) naphthenes, (d) aromatics for gas oil A by experiment and predictions.
Processes 10 00669 g005
Table 1. The bulk properties of gas oils.
Table 1. The bulk properties of gas oils.
Gas Oil AValueGas Oil BValue
Elemental analysis/wt% Elemental analysis/wt%
C87.4C85.7
H12.5H13.1
S0.14S1.3
Density/kg·m−30.862Density/kg·m−30.83
H/C ration/mol·mol−11.77Aromatic sulfur content/wt%0.15
PINA analysis/wt% PINA/wt%
P17.3P24.3
I21.3I29.7
N17.3N13.5
A44.1A32.5
ASTM-D2887/K ASTM-D2887/K
0.5%3990.5%409
10%55010%535
30%59530%576
50%62550%602
70%65070%622
90%68490%655
99.5%77599.5%705
Table 2. The structural attributes and PDFs.
Table 2. The structural attributes and PDFs.
IndexSAThe Type of PDFsParameters
1The type of moleculeHistogramx1, x2, x3
2If accept a naphthenic ring on the aromatic ringHistogramx4
3If the paraffin chain is branchedHistogramx5
4The type of sulfur-containing moleculeHistogramx6, x7
5The number of naphthenic ringsGammax8, x9
6The number of aromatic ringsGammax10, x11
7The length of sidechainGammax12, x13
8The length of paraffin chainGammax14, x15
9The length of sulfur-containing chainGammax16, x17
Table 3. The boundaries of parameters in gamma distributions for gas oil A.
Table 3. The boundaries of parameters in gamma distributions for gas oil A.
Index of SASPLBUBMeanLBUB
5x805, 10, 20, 50, 100x904
6x1005, 10, 20, 50, 100x1105
7x1215, 10, 20, 50, 100x13036
8x1415, 10, 20, 50, 100x15742
9x1615, 10, 20, 50, 100x17541
Table 4. The boundaries of parameters in gamma distributions for gas oil B.
Table 4. The boundaries of parameters in gamma distributions for gas oil B.
Index of SASPLBUBMeanLBUB
5x805, 10, 20, 50, 100x904
6x1005, 10, 20, 50, 100x1105
7x1215, 10, 20, 50, 100x13024
8x1415, 10, 20, 50, 100x15731
9x1615, 10, 20, 50, 100x17530
Table 5. Comparison between experimental values and estimated values of bulk properties in different cases for gas oil A.
Table 5. Comparison between experimental values and estimated values of bulk properties in different cases for gas oil A.
PropertyEVUB of SP for Gas Oil A
5102050100
Elemental analysis/wt%
C87.487.0487.0787.0587.0387.02
H12.512.8212.7912.8112.8312.84
S0.140.140.1410.140.1380.14
Density/kg·m−30.8620.8680.8640.8630.8610.862
H/C ratio/mol·mol−11.771.771.761.771.771.77
PINA/wt%
P17.317.217.217.517.317.3
I21.321.321.321.221.721.0
N17.317.317.317.517.317.6
A44.144.244.243.843.744.1
ASTM-D2887/K
0.5%399392406406406406
10%550522539541546554
30%595588600602605607
50%625623627627627632
70%650661655651650651
90%684704707699687707
99.5%775754773767754773
The sum of relative errors 0.04410.03460.03320.03410.0335
Table 6. Comparison between experimental values and estimated values of bulk properties in different cases for gas oil B.
Table 6. Comparison between experimental values and estimated values of bulk properties in different cases for gas oil B.
PropertyEVUB of SP for Gas Oil B
5102050100
Elemental analysis/wt%
C85.785.285.285.685.085.0
H13.113.513.513.213.713.7
S1.31.341.311.281.281.3
Density/kg·m−30.830.8350.8330.8520.8290.832
Aromatic sulfur content/wt%0.150.150.150.150.150.15
PINA/wt%
P24.324.324.224.224.324.6
I29.729.729.229.630.029.0
N13.513.314.013.513.413.6
A32.532.732.632.7432.332.8
ASTM-D2887/K
0.5%409403403403403403
10%535505506546561546
30%576554557582582568
50%602590595607600595
70%622627629627631633
90%655673672653664661
99.5%705710709694705701
The sum of relative errors 0.05910.05400.05230.04930.0511
Table 7. The comparison between experimental values and estimated values of bulk properties for gas oil A in different cases.
Table 7. The comparison between experimental values and estimated values of bulk properties for gas oil A in different cases.
PropertyEVUB of SP for Gas Oil A
5102050100
Elemental analysis/wt%
C87.487.187.187.187.087.06
H12.512.812.812.812.812.8
S0.140.1380.1410.1400.1410.138
Density/kg·m−30.8620.8590.8630.8580.8550.854
H/C ratio/mol·mol−11.771.761.761.771.771.77
Detailed group-type analysis/wt% --
P17.0017.4517.2817.7617.5117.34
I21.2021.2121.1221.7121.0621.37
MN15.3915.2615.5815.1515.0215.34
DN2.202.072.192.252.242.20
MA16.8015.2713.9714.6516.1717.31
NA7.408.968.569.149.178.07
DA8.438.9110.078.748.588.85
NDA6.406.006.425.805.474.40
TrA3.924.013.883.873.924.14
NTrA + TeA0.740.860.930.930.860.98
ASTM-D2887/K
0.5%399396421421396396
10%550503522521507503
30%595576593585588595
50%625629635624629629
70%650674677655662655
90%684729727692699687
99.5%775781779730741748
The sum of relative errors 0.12530.13400.13280.11700.1372
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, G.; Yang, M.; Qian, F. Study on the Domain of Parameters in the Stochastic Reconstruction Method. Processes 2022, 10, 669. https://doi.org/10.3390/pr10040669

AMA Style

Zhao G, Yang M, Qian F. Study on the Domain of Parameters in the Stochastic Reconstruction Method. Processes. 2022; 10(4):669. https://doi.org/10.3390/pr10040669

Chicago/Turabian Style

Zhao, Guangyao, Minglei Yang, and Feng Qian. 2022. "Study on the Domain of Parameters in the Stochastic Reconstruction Method" Processes 10, no. 4: 669. https://doi.org/10.3390/pr10040669

APA Style

Zhao, G., Yang, M., & Qian, F. (2022). Study on the Domain of Parameters in the Stochastic Reconstruction Method. Processes, 10(4), 669. https://doi.org/10.3390/pr10040669

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop