Next Article in Journal
Guiding Chart for Initial Layer Choice with Nanoimprint Lithography
Previous Article in Journal
Suitability of Different Titanium Dioxide Nanotube Morphologies for Photocatalytic Water Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic Properties Investigation of Process Volatile Organic Compounds (VOCs) and Its Transport Impact Factor in Oil Sands Management

1
Department of Civil Engineering, Tongling University, No. 4, Cui Hu Road 1335, Tongling District, Tongling 244000, China
2
Donadeo Innovation Centre of Engineering, University of Alberta, 9211-116 Street NW, Edmonton, AB T6G 1H9, Canada
3
CanmetENERGY Devon Research Centre, Natural Resources Canada, 1 Oil Patch Drive, Devon, AB T9G 1A8, Canada
4
Institute of Molecular Engineering and Applied Chemistry, Anhui University of Technology, Ma’anshan 243002, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(3), 709; https://doi.org/10.3390/nano11030709
Submission received: 3 February 2021 / Revised: 20 February 2021 / Accepted: 5 March 2021 / Published: 11 March 2021

Abstract

:
This paper presents a new approach for the determination of volatile organic compounds (VOCs) characteristics and their migration influencing factors in oil sands management processes and reveals the relationship between different asphaltene content and different solvents. Specifically, thermodynamic (i.e., partitioning coefficients, Kr, specific retention volume, Vg, the activity coefficients, γ and enthalpy of solution, ΔH0) and their impact factors are discussed. Gas-liquid chromatography (GLC) experimental measurements were used as the test data. A range of solvents (nC5, iC5, nC6, nC7, and Toluene) has been tested in different asphalt contents (0, 2.56, 9.93, 36.86, 53.67 wt%). There are temperatures in the range of 333.2–393.2 K (with 10 K increase) were conducted, respectively. The dynamics properties of asphalt mixture are calculated, and the relation between dynamics properties of asphalt mixture and absolute temperature, asphalt content and solvent type is discussed. The results show that within the acceptable error range, partitioning coefficients, Kr, specific retention volume, Vg, and enthalpy of solution, ΔH0 and other thermodynamic properties have a good tendency to predict, they decrease with the increase in asphaltene content and temperature and increase with the increase in solute carbon number.

1. Introduction

During froth treatment and water management processes, oil sands operators try their best to recover as much solvent as possible from the wastewater (i.e., tailings), but a significant amount of light solvent remains are discharged into the environment. Therefore, emissions have a significant impact on air, water, soil pollution and the greenhouse effect [1,2,3]. Therefore, a full understanding of the mechanism of light hydrocarbon release from tailings especially their thermodynamic properties is helpful to find ways to improve the performance of Tailings Solvent Recovery Units (TSRU) [4,5], to reduce the number of light hydrocarbons left in the tailings pond before deposition and to reduce the cost of producing bitumen [6]. In this study, asphaltenes in the form of asphalt or bitumen products are the most valuable component of crude oil, and asphaltene is also considered to be the “heaviest component of petroleum asphalt, insoluble in n-alkanes such as pentane and n-heptane, but soluble in aromatic hydrocarbons such as toluene or xylene.” [6]. The most important is that the highest molecular weight asphaltene components are usually classified as saturated-aromated-resin-asphaltene [7,8]. Asphaltene has complex and highly dispersed components that can form clusters in the nanometer length scale, often defined as “nanoaggregates” [9,10,11,12]. In addition, they can form large aggregates and flocculation, leading to contamination of oil and gas reservoirs, reactors and pipelines, and even serious problems in the transportation and processing of petroleum fluids. It is also difficult to predict when and where these problems will occur during the production process, so it is essential to better understand how its dynamics equilibrium and kinetic processes work to avoid these serious problems [13,14]. The proper dynamics especially the thermodynamic conditions and the choice of the internal structure of TSRU will also affect the kinetics of the asphaltene dilution release process, which plays a crucial role in its performance [15,16]. Dynamics properties of solvents with different proportions of asphaltene are important parameters for modeling and improving the solvent recovery process [17,18].
Normally dynamics properties were tested by physical or chemical methods [17,19,20,21]. In recent decades, inverse gas-liquid chromatography (IGC) has become a powerful method to evaluate the properties of various materials [22,23,24]. Therefore, IGC is a useful tool for studying especially thermodynamic properties in this study. The IGC allows access to more thermodynamic parameters (i.e., partitioning coefficients, Kr, specific retention volume, Vg, the activity coefficients, γ and enthalpy of solution, ΔH0) [25,26,27].
For dynamics especially its thermodynamic parameters research, for example, the thermodynamic properties of light hydrocarbons in water/asphalt systems were investigated [20]. Yao et al. further verified the thermodynamic properties of bituminous binder components, conducted molecular dynamics modeling of asphaltenes, and eliminated the uncertainty of some parameter values [28]. Furthermore, the successful application of thermodynamic properties research is limited by its single involved factor. Specifically, Hicks and Young studied the effects of different n-butane, n-pentane, and n-hexane on n-tetracyclane and n-octadecane [29]; However, each measurement is carried out at a separate temperature, so it is almost impossible for these data to be correlated due to the temperature dependence of thermodynamic properties, such as the activity coefficient of low molecular weight solvents at low temperatures [25,26,30]. Besides, most literature reports only data on the activity coefficient of one or more solvents in the high temperature range (one dynamics property) [29,31,32]. Furthermore, Xu and Mundhwa et al. published alkane-solute data at different temperature ranges [33,34]. These literatures were compared with the results of this study. T, carbon number and other single factors have been reported by some researchers in the literature [26,35]. However, there are few studies on asphaltene contents in different solvents and at different temperatures. Therefore, the current research focuses on the use of IGC to study the interactions of different influencing factors with different dynamics properties at different temperatures, from 333.2 K to 393.2 K. In addition, asphalt, an asphaltene derivative, is the material of choice for the experimental samples in this study. Several chemical formulations have been found most effective in dissolving asphaltenes with pure aromatic hydrocarbons such as toluene and xylene. Therefore, in this study, five asphalt samples with different asphaltene contents (0, 2.56, 9.93, 36.86, 53.67 wt%) were used for the fixed phase (solvent/liquid phase) of IGC. In the mobile phase (solute), five light hydrocarbons (2-methylbutane, n-pentane, n-hexane, n-heptane, and toluene) were selected. Finally, the effects of all those different impact factors on different thermodynamic properties such as partitioning coefficients, Kr, specific retention volume, Vg, the activity coefficients, γ and enthalpy of solution, ΔH0 were compared and evaluated.
Certain characteristics of oil sand tailings affect the emission rate of VOCs, as well as the residual bitumen in the tailings. It traps volatile organic compounds and slowly releases them into the environment. Thermodynamic parameters are necessary to improve the available technology and predict the distribution of VOCs in solvent recovery process. Therefore, the aim of this research is to start from the basic research to understand its mechanism, to guide the industry to develop new solvent recovery technology. In summary, a better understanding of its thermodynamic properties and their influencing factors will help to optimize or improve solvent recovery and process efficiency, thereby reducing air, water, and soil pollution.

2. Research approach

The thermodynamic properties of hydrocarbon vapors in asphalt were mainly tested by two approaches: IGC experimental method and model calculation. A flow chart of the research approach is list in Figure 1.

2.1. Materials

Industrial grade (> 99.5%) toluene was used for SARA analysis and solid phase removal. The mass ratios of solvent to asphalt (S/B) were 1.7, 3.5 and 40, respectively, after the dried asphalt was treated with n-pentane. The density of solid asphaltenes ranges from 1100 to 1280 kg/m3. After treatment, two fractions were obtained for each S/B: bitumen partially deasphalted and asphaltene-rich deposits of n-pentane evaporated by rotary evaporation. The asphaltene contents were 9.9, 2.6 and 0 wt%, respectively, and the corresponding asphaltene contents in the precipitation were 36.9 wt%, 53.7 wt% and 84.2 wt%, respectively. These partial deasphaltenes and two asphaltene-rich deposits (36.9 and 53.7 wt% asphaltenes) were used as stationary phases of the prepared gas chromatographic column. Table 1 lists five solvents used to prepare stationary phases in gas chromatographic columns, including 2-methylbutane, toluene, pentane, n-hexane, and n-heptane. The purity of n-pentane and n-heptane were 99.5 and 99.7 wt%, respectively. For example, the boiling points of pentane and toluene in Table 1 are 36.1 and 110–111 °C, respectively [36].
Argon and helium gas were purchased from BOC gas. Acid-washed Chromosphere W was pickled and treated with dimethylchlorosilane (DMCS) by the manufacturer of Sigma-Aldrich (Sigma-Aldrich Canada Co., Oakville, ON, Canada).
Five asphalt samples with different asphaltene contents were studied. The high content was 53.7 wt% (% by weight) asphaltene and 36.9 wt% (% by weight) asphaltene, medium content was 9.9 wt% asphaltene, low content was 2.6 wt% and 0 wt% asphaltene. Five different hydrocarbon solvents were then mixed with the five different asphaltenes. For solvents, 2-methyl-butane (iC5), n-pentane (nC5), n-hexane (nC6), toluene, and n-heptane (nC7) are certified by the American chemical society (ACS) and provided by Fisher. The asphaltene content in asphalt was determined by standard pentane precipitation and filtration methods (ASTM D2007-80). In the solvent, the solvent with high solubility and low boiling point is the ideal solvent for deasphaltene, which is the partial deasphaltene or asphaltene rich, precipitates (36.9 and 53.7 wt% asphaltene). The asphaltene content in bitumen was determined by standard pentane precipitation and filtration methods (ASTM D2007-80). A total of thirteen columns with different asphalting percentage were prepared and used for specific retention volume measurements. Table 2 lists the performance of these columns, including the type of bitumen, the asphaltene content, bitumen load in the packing, weight of packing (bitumen coated Chromosorb W), and weight of bitumen in the column (stationary phase). The average molecular weight of each bitumen sample is also included.

2.2. Experiment Apparatus and Procedures

Five solvents (nC5, iC5, nC6, nC7 and toluene) with different asphaltene contents (0, 2.56, 9.93, 36.86 and 53.67 wt% asphaltene) were determined at different temperatures (20, 35, 50, 65 and 80 °C). Therefore, a total of 13 experiments were conducted at CanmetENERGY in Devon, Alberta. Peak area and retention time were measured by Agilent 6890 gas chromatography with electronic pressure control (EPC)(Alicat Scientific, Inc., Tucson, AZ, USA), thermal conductivity detector (TCD) and flame ionization detector (FID)( GenTech Scientific LLC, New York, NY, USA). Jour digital flow meter was used to measure the flow at the column outlet. The thermodynamic properties of hydrocarbon vapours in the asphalt were determined by the GLC method. The asphalt as the stationary phase is coated on in a thin layer on Chromsorb W.
Firstly, the asphalt was dissolved in toluene at the calculated concentration to reach the predetermined asphalt load. The toluene asphalt solution was then poured into the pre-weighed Chromsorb W, with the liquid interface is slightly over the solid interface. The Chromsorb W was soaked in the toluene bitumen solution for 24 h. Toluene is then rotated to evaporate slowly at 120 °C and 80 kPa. When most of the toluene is removed and there is no more toluene dripping from the condenser, the pressure is reduced to 6 kPa for half an hour. After the bitumen coated Chromsorb W (g) was placed in a vacuum drying oven at 60 °C and 34 KPa (absolute) for tween-four hours to remove the final toluene. The vacuum oven is under a nitrogen atmosphere. The amount of bitumen coated was determined by measuring the weight loss of a sample on a heating furnace of 2 h at 800 ºC. The column is made of chromatographer-grade stainless steel tube, 2 m long, ¼ inch outside diameter. The tubes were cleaned with deionized water, concentrated nitric acid, deionized water, concentrated ammonium hydroxide, deionized water, and acetone, respectively. A clean stainless-steel tube is then filled with asphalt coated Chromsorb W. one end of the tubing is plugged with a small amount of glass wool and connected to the vacuum pump. The asphalt-coated chrome sand is then slowly poured into the other end of the pipe through a small funnel. During filling, the column is tapped until no more powder can enter the tube. After packing, another glass wool plug is placed on the open end of the tubing. Record the weight of powder. The filled column is coiled and installed under GC conditions at a carrier gas flow rate of 20 mL/min at 100 ºC until the baseline is stable.
Control the carrier gas flow consistently and record the column inlet pressure under each working condition. The dead volume of the packed column was determined by injecting 2 µL argon gas with a 10 uL syringe. Small samples of each investigated hydrocarbon were injected using an automatic sampler at least three times injections to ensure retention time correction. Due to changes in local temperature, air pressure and humidity, the measurements fluctuate slightly. For example, the temperature range is ±0.01 k, and the uncertainty of pressure measurement is 0.01 kPa.

2.3. Thermodynamic Methods

To further verify the dynamical model of asphaltene, to eliminate the uncertainty of its parameter values, we use parameters, such as specific retention volume, Vg, solute partition coefficient Kr, and Activity coefficient γ values of the investigated hydrocarbons in bitumen for five different asphalting contents (0, 2.56, 9.93, 21.2, 36.86, 53.67, and 84.2 wt%) to check the thermodynamic properties. It is helpful to discover the change rule of oil sand production, reveal the influence mechanism of oil sand production, and deepen the understanding of the oil sand production process. Initially, Vg0 is the volume of the dry carrier gas required to wash and extract the solute sample at 0 °C. The following formula (1–3) is used for the research [37].
V g 0 = V N w L ( 273.15 T )
where,
V N = V R 0 V M = j ( V R V A ) = j F c ( t R t A ) = j F T T f m ( t R t A )
j = p o p = 3 2 [ ( p i / p o ) 2 1 ( p i / p 0 ) 3 1 ]
and, WL is the mass (g) of the stationary phase (bitumen) in the column, V N is the net retention volume (cm3), V R 0 is the corrected retention volume (cm3), V M is the mobile phase volume(cm3), V R is the retention volume (cm3), Va is the adjusted retention volume (cm3), Fc is the carrier gas flow rate (cm3 /min), F is the measured flow rate(cm3 /min), T is the column temperature, T f m is the flow meter temperature, tR is the retention time (min) for a solute; tA is the column dead time (min); j is the compressibility correction factor; po is the outlet pressure (Pa); and pi is the pressure of the carrier gas at the inlet to the GC column. In this study, based on the results of specific retention volume, V g 0 from the Equation (1), Partition coefficient, K R can be calculated through the following Equation (4). ρ L is the density of stationary phase (g/cm3) and T is the column temperature (K) [38,39].
K R = V g 0 ρ L T 273.15
Then, according to Equation (4), infinite dilute activity coefficient γ p can be calculated through the following Equation (5).
γ p = ρ L R T K R p 1 0 M W L
where ρ L is the density of stationary phase (g/cm3); MWL is the molecular weight of the stationary phase (g/mol); R is the ideal gas constant; p 1 0 is the saturated vapor pressure of the pure solute (Pa) at the column temperature T (K).
It can also be seen from equation 4 that the equilibrium constant K R is proportional to the specific retention volume Vg0. Therefore, the plot of lnVg0 against 1/T should give the same slope of ln K R versus 1/T. Therefore, then it has been widely used in the following Equation (6) to get the value of enthalpy of solution, ΔH0 [40,41].
Δ H 0 = R d ln V g 0 d ( 1 / T )

3. Results

Thirteen tests were prepared for nC5, iC5, nC6, nC7, and Toluene under the different asphalting contents (0, 2.56, 9.93, 36.86, 53.67 and 84.2 wt%) and at different temperature (333.2 K, 343.2 K, 353.2 K, 363.2 K, 373.2 K, 383.2 K, and 393.2 K). Those thirteen tests are used for the specific retention volume Vg0 measurements according to Equations (1)–(3). Then, according to the specific retention volume Vg0, various thermodynamic properties of the investigated hydrocarbon vapors dissolving into the bitumen, such as partition coefficient Kr and infinite dilution activity coefficient γ p are calculated according to Equations (4) and (5). Finally, the results of the measurements of partition coefficient Kr and infinite dilution activity coefficient γ have been used for thermodynamic calculations of enthalpy of solution ΔH0 according to Equation (6). Those thermodynamic properties especially partitioning coefficient Kr, and activity coefficient γ of the same hydrocarbons in bitumen were discussed in detail.

3.1. Specific Retention Volume

Uffink et al. describes how macroscopic observations, such as temperature and pressure, relate to microscopic parameters that fluctuate around the mean. It relates thermodynamic quantities (such as volume retention) to microscopic behavior [42]. Therefore, in the study, within the specified temperature range, from 333.2 k to 383.2 k, the obtained specific retention volume values are shown in Table 3 and Figure 1, respectively. Table 3 lists comparisons between referenced data [34] of Specific retention volume V g 0 and the calculated V g 0 of investigated hydrocarbons in bitumen under 21.2 % asphalting contents. It could be noted that the measured V g 0 are in good agreement with the reference data V g 0 are shown in Table 3. The results also confirmed the reliability of the measurement, which shows that the macroscopic observed values such as temperature and pressure are related to the thermodynamic quantities such as specific retention volume.
The specific retention volume V g 0 of each i-C5, n-C5, n-C6, n-C7 and toluene compounds under the different asphalting contents (0, 2.56, 9.93, 36.86, 53.67, and 84.2 wt%) were measured and the retention diagrams ln ( V g 0 ) = f(1/T) were plotted in Figure 2. Straight lines were obtained in Figure 2. Furthermore, for each i-C5, n-C5, n-C7, and toluene it could be noted that the specific retention volume V g 0 in bitumen that has 2.56 wt% of asphaltene contents is much higher than in the bitumen that has 53.67 wt% asphaltene contents. The bitumen with high asphaltene content has poor retention of hydrocarbons studied.
At last, As seen in Figure 1, Those hydrocarbons with more carbon atoms have higher V g 0 values in these adsorbents, whereas the V g 0 values for light hydrocarbons (≤C5) are almost always less than 100 1/g, which shows that and V g 0 values increase with the chain length of the solute alkane chain (carbon number).

3.2. Partition Coefficient

It could be noted that in practical case studies, nonpolar or hydrogen bond interactions become important in asphaltene-solvent interactions. In addition, the addition of sufficient non-polar solvents, such as toluene, can inhibit the accumulation of asphaltenes [13]. In fact, these standard chemical reactions, such as partition coefficients, are often used by the asphalt industry: to make them more soluble, alkene substituents, such as pentane, can be added. The results show that increased bond molecular interactions can significantly increase solubility. Therefore, in the study, the partition coefficient KR was calculated from the specific retention volume of Equation (4). The calculated Partition coefficients KR as functions of different asphaltene content (0, 2.56, 9.93, 36.86, and 53.67 wt%) in bitumen for the five hydrocarbons (nC5, iC5, nC6, nC7, and Toluene) are shown in Figure 3. The partition coefficient KR decreases monotonically as the temperature increases. The study also found that in asphalt containing 53.7 wt% asphaltene, the distribution coefficient is usually only 60–70% of those in the bitumen with lower asphaltene contents.
Similarly, in Figure 3, alkanes and light alkanes have lower distribution coefficients KR, while toluene has the highest distribution coefficient KR. Additionally, can be seen from Figure 3, the descending order of partition coefficient is: toluene, n-heptane (nC7), n-hexane (nC6), n-pentane (nC5), and 2-methylbutane (iC5). This trend indicates that the role of the molecular geometry of the compounds studied is significant [32].

3.3. Infinite Dilution Activity Coefficients

The infinite dilute activity coefficient γ p could be calculated in Table 4 according to Equation (5) based on the value of the former partition coefficient KR. It can be seen from Table 4 that the value of the infinite dilution activity coefficient γ p is related to the content of asphaltene, the type of solute compound (carbon number) and the temperature as well. The saturated vapor pressure of the solute at the column temperature is required. The vapor pressure of the five hydrocarbons at column temperature (Equation (5)) is calculated using the familiar Antoine equation log p 1 0 = A − B/ (t +C), where t is the temperature in degrees Celsius, and A, B and C are known as the Antoine constants, which can be obtained from the literature of Xu et al. [33].
The values of the infinite dilute activity coefficient obtained in the specified temperature range and the associated different asphaltene contents for five hydrocarbons are presented in Figure 4, respectively. Figure 3 shows the linear relationship between the natural logarithm of the infinite dilute activity coefficient, ln γ p , vs the inverse of absolute temperature T. The highest values of infinite dilute activity coefficient γ p were obtained for hydrocarbons are with highest asphaltene contents. For the activity coefficient, γ, in bitumen has 18 wt% of asphaltene are much higher than in the bitumen has 84 wt% asphaltene. It is also found that the infinite dilute activity coefficient γ p decrease with increasing temperature. The results show that the activity coefficient of infinite dilution increases with the increase in carbon number in a normal alkane. This is because the presence of double and triple bonds in the solute structure significantly reduces the value of the infinite dilution activity coefficient. This is a typical ionic liquid with a weak solution-solvent interaction [7].
The highest values of γ were obtained for normal n-alkanes are for nC7 at the highest asphaltene contents and highest temperature, which was typical of ionic liquid interaction with weak solute-solvent.

3.4. Enthalpy of Solution, ΔH0

Enthalpy of solution ΔH0 means that in a diatomic molecule AB, the bond energy corresponds to the energy of dissociation in the two atoms A and B, and it indicates the work performed to take them at infinite distance. Like the infinite dilution activity coefficient, the solution enthalpy of a given hydrocarbon compound varies with the change of the carbon number and asphaltene content in the bitumen [7]. Figure 5 shows the comparison of different asphalting contents of five different hydrocarbons (nC5, iC5, nC6, nC7, and toluene). Figure 5 show that the enthalpy of solution ΔH0 increases with the carbon number in the compound increases. This confirms the occurrence of stronger attractive interactions with the ionic liquid.

4. Discussion

Based on the above results, it concluded that a good linear relationship is obtained between the natural logarithm of specific retention volume, V g , the partitioning coefficients, Kr, the activity coefficients, γ at infinite dilution, and Enthalpy of solution ΔH0. Ln V g , Ln Kr, LnΔH0, Ln γ and the inverse of the absolute temperature, i.e., Ln V g , Ln Kr, LnΔH0, Ln γ can be expressed like Harned and Robinson [43] as: Ln V g = a + b/(T/K); Ln Kr = a + b/(T/K); and Ln γ = a + b/(T/K) in which a and b are constants (Figure 6). For example, as shown in Figure 6, V g increase with the temperature decrease. Furthermore, as shown in Figure 7, with temperature decrease, Kr, the activity coefficients increase as well.
In this research, V g , Kr and ΔH0 decreases with increasing asphaltene content in the bitumen phase, however, γ will decrease in the opposite way. Furthermore, all the partitioning coefficients, Kr, and specific retention volume, V g the activity coefficients, and enthalpy of solution, γ will decrease with the increasing of temperature and increase with the chain length of the solute alkane chain (carbon number). This phenomenon can be explained by the change of the size ratio of solvent to solute and the interaction between solute and solute. This may also indicate that it is easier to extract these hydrocarbons from the high-asphaltene bitumen found in the TSRU tailings than from the low-asphaltene bitumen from the diluted asphaltene bitumen products [44]. Several reasons might also concur with such outcomes. Firstly, the high temperatures could reduce hydrogen-bonded hydration water, resulting in a decrease in rheological parameters [45].
IGC is a simple, fast, and accurate method to study the thermodynamic interaction parameters and other physicochemical properties of solute in solvent and polymer blends [46,47]. The results show that IGC is a new technique to determine the thermodynamic parameters of asphalt, and the change of asphaltene content and solute has obvious influence on the thermodynamic properties of asphalt. The experimental results are in good agreement with the thermodynamic model, indicating that the solute with high asphaltene content has less retention than that with low asphaltene content [43]. The comparison between the IGC experimental results and the model results shows that when thermodynamic factors are the controlling factors, these process VOCs are easier to be removed from TSRU tailings than traditional asphalt.

5. Conclusions

Dynamics parameters of five light hydrocarbons (2-methyl-butane, n-pentane, n-hexane, n-heptane, and toluene) with asphaltene content between 0.0 wt and 53.7% were studied. The values of each thermodynamic parameter in the range of 303.2 K−393.2 K were determined by Inverse gas chromatography (IGC). Firstly, the specific retention volume, V g of hydrocarbons extracted from asphalt with different asphalt contents (0, 2.56, 9.93, 36.86, 53.67, 84.2 wt%) is in good agreement with the experimental data in the literature. Then, the specific retention volume value is used to calculate the thermodynamic properties of hydrocarbon dissolves in asphalt, including partition coefficient, Kr, infinite dilution activity coefficient, γ and enthalpy of solution, ΔH0.
This study enables us to assess the nature of asphaltene content in the bitumen and the solute compounds’ interactions (hydrogen bonds, solute and solvent interactions). All dynamics parameters (i.e., partitioning coefficients, Kr, specific retention volume, V g , the activity coefficients, γ and enthalpy of solution, ΔH0) of normal alkane increase with increasing chain length. This will help to understand the different behavior of the solute alkanes considering their function and chemical structure, which could be expressed in detail as follows: (1) the results show that the change of asphaltene content and temperature have different effects on the thermodynamic properties of asphalt, and (2) the change of asphaltene content between different light hydrocarbon (carbon number) and asphalt also plays a key role in the control of thermodynamic properties. In summary, a better understanding of its thermodynamic properties and their influencing factors will help to optimize or improve solvent recovery and process efficiency to less VOCs emission, thereby reducing air, water, and soil pollution.

Author Contributions

Conceptualization, investigation, writing—original draft-preparation, writing, J.Y.; funding acquisition, Y.S.; review and editing, Y.S., Y.J. and Q.Z.; All authors have read and agreed to the published version of the manuscript. The author Y.J. also thanks Tadeusz Dabros (who passed away) for supervising the tailings fractions and column preparing.

Funding

This research was funded by Anhui Provincial Natural Science Foundation (No. 1608085ME103) and Canadian government’s interdepartmental Program of Energy Research and Development (PERD), Natural Resources Canada. The APC was funded by Anhui Provincial Top-notch Talent Plan (No. gxbjZD40).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jing, Y.; Maria, E. Extended Environmental Multimedia Modelling System assessing the risk carried by pollutants in interacted air-unsaturated-groundwater zones. J. Hazard. Mater. 2020. [Google Scholar] [CrossRef]
  2. Bian, X.Y.; Mbadinga, S.M.; Liu, Y.F.; Yang, S.Z.; Liu, J.F.; Ye, R.Q.; Gu, J.D.; Mu, B.Z. Insights into the anaerobic biodegradation pathway of n-Alkanes in oil reservoirs by detection of signature metabolites. Sci. Rep. 2015, 5, 9801. [Google Scholar] [CrossRef]
  3. Brandl, F.; Bertrand, N.; Lima, E.M.; Langer, R. Nanoparticles with photoinduced precipitation for the extraction of pollutants from water and soil. Nat. Commun. 2015, 6, 7765. [Google Scholar] [CrossRef]
  4. Components, T.; Ansari, N. Thermogravimetric Analysis of Solvent Interaction with Model TSRU Tailings Components. Master’s Thesis, University of Alberta, Edmonton, AB, Canada, 2013. [Google Scholar]
  5. Yuan, J.; Wang, X.; Vakili, M.; Segun, G.A.; Zhong, L. Volatile organic compounds (VOCs) releasing model from tailings solvent recovery unit (TSRU) tailings and its sensitivity analysis in environment management. Process Saf. Environ. Prot. 2020. [Google Scholar] [CrossRef]
  6. Rajabi, H.; Mosleh, M.H.; Mandal, P.; Lea-Langton, A.; Sedighi, M. Emissions of volatile organic compounds from crude oil processing—Global emission inventory and environmental release. Sci. Total Environ. 2020, 727, 138654. [Google Scholar] [CrossRef] [Green Version]
  7. Davda, R.R.; Shabaker, J.W.; Huber, G.W.; Cortright, R.D.; Dumesic, J.A. A review of catalytic issues and process conditions for renewable hydrogen and alkanes by aqueous-phase reforming of oxygenated hydrocarbons over supported metal catalysts. Appl. Catal. B Environ. 2005, 56, 171–186. [Google Scholar] [CrossRef]
  8. Mohamed, R.S.; Ramos, A.C.S.; Loh, W. Aggregation behavior of two asphaltenic fractions in aromatic solvents. Energy Fuels 1999, 13, 323–327. [Google Scholar] [CrossRef]
  9. Eyssautier, J.; Levitz, P.; Espinat, D.; Jestin, J.; Gummel, J.; Grillo, I.; Barré, L. Insight into asphaltene nanoaggregate structure inferred by small angle neutron and X-ray scattering. J. Phys. Chem. B 2011, 115, 6827–6837. [Google Scholar] [CrossRef]
  10. Betancourt, S.S.; Ventura, G.T.; Pomerantz, A.E.; Viloria, O.; Dubost, F.X.; Zuo, J.; Monson, G.; Bustamante, D.; Purcell, J.M.; Nelson, R.K.; et al. Nanoaggregates of asphaltenes in a reservoir crude oil and reservoir connectivity. Energy Fuels 2009, 23, 1178–1188. [Google Scholar] [CrossRef]
  11. Headen, T.F.; Boek, E.S.; Skipper, N.T. Evidence for asphaltene nanoaggregation in toluene and heptane from molecular dynamics simulations. Energy Fuels 2009, 23, 1220–1229. [Google Scholar] [CrossRef]
  12. Caetano, M.; Ranaudo, M.; Jaimes, B. Simulation of asphaltene aggregation and related properties using an equilibrium-based mathematical model. Energy Fuels 2011, 25, 3544–3551. [Google Scholar]
  13. Castillo, J.; Hung, J.; Fernández, A.M.V. Nonlinear optical evidence of aggregation in asphaltene—Toluene solutions. Fuel 2001, 80, 1239–1243. [Google Scholar] [CrossRef]
  14. Gawrys, K.L.; Kilpatrick, P.K. Asphaltene aggregation: Techniques for analysis. Instrum. Sci. Technol. 2003, 32, 247–253. [Google Scholar] [CrossRef]
  15. Rossini, F.D.; Bureau, N.; Eyring, H. Thermodynamic properties of hydrocarbons. Encycl. Hydrocarb. 2005, V, 77–94. [Google Scholar]
  16. Pathak, M.; Kweon, H.; Deo, M.; Huang, H. Kerogen Swelling and Confinement: Its implication on Fluid Thermodynamic Properties in Shales. Sci. Rep. 2017, 7, 12530. [Google Scholar] [CrossRef] [Green Version]
  17. Yuan, J.; Zhong, L.; Vakili, M.; Segun, G.A. New modeling method to simulate asphaltenes at oil sands process in water management. J. Mol. Graph. Model. 2019, 91, 1–9. [Google Scholar] [CrossRef]
  18. Shekaari, H.; Zafarani-moattar, M.T.; Shayanfar, A.; Mokhtarpour, M. Effect of choline chloride/ethylene glycol or glycerol as deep eutectic solvents on the solubility and thermodynamic properties of acetaminophen. J. Mol. Liq. 2018, 249, 1222–1235. [Google Scholar] [CrossRef]
  19. Kooijman, H.A.; Taylor, R. Estimation of Diffusion Coefficients in Multicomponent Liquid Systems. Ind. Eng. Chem. Res. 1991, 30, 1217–1222. [Google Scholar] [CrossRef]
  20. Yuan, J.; Elektorowicz, M.; Chen, Z.; Segun, G.A.; Vakili, M.; Zhong, L.; Wang, B.; Zhu, J.; Wu, Y. Simulation and computer modeling of asphaltene in different solvents on oil-water interfaces using a molecular dynamic methodology. J. Mol. Graph. Model. 2019, 93, 107450. [Google Scholar] [CrossRef]
  21. Öztürk, E.; Ocak, H.; Cakar, F.; Karanlık, G.; Cankurtaran, Ö.; Bılgın-Eran, B. Investigation of thermodynamic properties of 4-decyloxybiphenyl-4’-carboxylic acid liquid crystal and preparation of polymer dispersed liquid crystal composite. J. Mol. Liq. 2018, 265, 24–30. [Google Scholar] [CrossRef]
  22. Belov, N.A.; Safronov, A.P.; Yampolskii, Y.P. Thermodynamics of sorption in an amorphous perfluorinated copolymer AF1600 studied by inverse gas chromatography. Macromolecules 2011, 44, 902–912. [Google Scholar] [CrossRef]
  23. Heintz, A.; Kulikov, D.V.; Verevkin, S.P. Thermodynamic properties of mixtures containing ionic liquids. 2. Activity coefficients at infinite dilution of hydrocarbons and polar solutes in 1-methyl-3-ethyl-imidazolium bis(trifluoromethyl-sulfonyl) amide and in 1,2-dimethyl-3-ethyl-imidazolium bis(trifluoromethyl-sulfonyl) amide using gas-liquid chromatography. J. Chem. Eng. Data 2002. [Google Scholar] [CrossRef]
  24. Belov, N.A.; Nizhegorodova, Y.A.; Kim, S.; Han, S.H.; Yampolskii, Y.P.; Lee, Y.M. Inverse Gas Chromatographic Study of Sorption Thermodynamics in Thermally Rearranged Polymer Based on 2,2-Bis(3-amino-4-hydroxyphenyl)-hexafluoropropane and 4,4′-Hexafluoroisopropylidene Diphthalic Anhydride. Ind. Eng. Chem. Res. 2013. [Google Scholar] [CrossRef]
  25. Heintz, A.; Kulikov, D.V.; Verevkin, S.P. Thermodynamic properties of mixtures containing ionic liquids. 1. Activity coefficients at infinite dilution of alkanes, alkenes, and alkylbenzenes in 4-methyl-n-butylpyridinium tetrafluoroborate using gas-liquid chromatography. J. Chem. Eng. Data 2001. [Google Scholar] [CrossRef]
  26. Mutelet, F.; Jaubert, J.N.; Rogalski, M.; Harmand, J.; Sindt, M.; Mieloszynski, J.L. Thermodynamic Properties of Mixtures Containing Ionic Liquids: Activity Coefficients at Infinite Dilution of Organic Compounds in 1-Propyl Boronic Acid-3-Alkylimidazolium Bromide and 1-Propenyl-3-alkylimidazolium Bromide Using Inverse Gas Chromatography. J. Chem. Eng. Data 2006. [Google Scholar] [CrossRef]
  27. Angle, C.W.; Long, Y.; Hamza, H.; Lue, L. Precipitation of asphaltenes from solvent-diluted heavy oil and thermodynamic properties of solvent-diluted heavy oil solutions. Fuel 2006, 85, 492–506. [Google Scholar] [CrossRef]
  28. Yao, M.; Woo, Y.C.; Tijing, L.D.; Choi, J.S.; Shon, H.K. Effects of volatile organic compounds on water recovery from produced water via vacuum membrane distillation. Desalination 2018, 440, 146–155. [Google Scholar] [CrossRef] [Green Version]
  29. Hicks, C.P.; Young, C.L. Activity coefficients of C4–C8 n-alkanes in C16–C32 n-alkanes. Trans. Faraday Soc. 1968, 64, 2675–2682. [Google Scholar] [CrossRef]
  30. Marcinkowski, Ł.; Eichenlaub, J.; Ghasemi, E.; Polkowska, Ż.; Kloskowski, A. Measurements of Activity Coefficients at Infinite Dilution for Organic Solutes in the Ionic Liquids N-Ethyl- and N-Octyl-N-methylmorpholinium Bis(trifluoromethanesulfonyl)imide. A Useful Tool for Solvent Selection. Molecules 2020, 25, 634. [Google Scholar] [CrossRef] [Green Version]
  31. Yu, Y.X.; Gong, Q.; Huang, L.L. Measurement of activity coefficient at infinite dilution of hydrocarbons in sulfolane using gas-liquid chromatography. J. Chem. Eng. Data 2007, 52, 1459–1463. [Google Scholar] [CrossRef]
  32. Heintz, A.; Kulikov, D.V.; Verevkin, S.P. Thermodynamic properties of mixtures containing ionic liquids. Activity coefficients at infinite dilution of polar solutes in 4-methyl-N-butyl-pyridinium tetrafluoroborate using gas-liquid chromatography. J. Chem. Thermodyn. 2002, 34, 1341–1347. [Google Scholar] [CrossRef]
  33. Xu, Y. Thermodynamic Investigations: Gases and Vapors in Liquids, Calorimetric Studies. Ph.D. Thesis, University of Alberta, Edmonton, AB, Canada, 1990. [Google Scholar]
  34. Mundhwa, M.; Elmahmudi, S.; Maham, Y.; Henni, A. Molar heat capacity of aqueous sulfolane, 4-formylmorpholine, 1-methyl-2-pyrrolidinone, and triethylene glycol dimethyl ether solutions from (303.15 to 353.15) K. J. Chem. Eng. Data 2009, 54, 2895–2901. [Google Scholar] [CrossRef]
  35. Gauden, P.A.; Terzyk, A.P.; Rychlicki, G.; Kowalczyk, P. Thermodynamic properties of benzene adsorbed in activated carbons and multi-walled carbon nanotubes. Chem. Phys. Lett. 2006, 421, 409–414. [Google Scholar] [CrossRef]
  36. Milne, G.W.A. (Ed.) Gardner’s Commercially Important Chemicals: Synonyms, Trade Names, and Properties; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
  37. Grigoryev, B.A.; Rastorguyev, Y.L.; Kurumov, D.S.; Gerasimov, A.A.; Kharin, V.E.; Plotnikov, S.A. Thermodynamic properties of n-pentane. Int. J. Thermophys. 1990, 11, 487–502. [Google Scholar] [CrossRef]
  38. Ulmschneider, M.B.; Ulmschneider, J.P.; Schiller, N.; Wallace, B.A.; Von Heijne, G.; White, S.H. Spontaneous transmembrane helix insertion thermodynamically mimics translocon-guided insertion. Nat. Commun. 2014, 5. [Google Scholar] [CrossRef]
  39. Daniel, R.; Edgar, F.; Wiley, G.C.J. The Chemical Thermodynamics of Organic Compounds. J. Chem. Thermodyn. 1970, 2, 153–154. [Google Scholar]
  40. Clarke, E.C.W.; Glew, D.N. Evaluation of thermodynamic functions from equilibrium constants. Trans. Faraday Soc. 1966, 62, 539–547. [Google Scholar] [CrossRef]
  41. Everett, D.H.; Wynne-Jones, W.F.K. The thermodynamics of acid-base equilibria. Trans. Faraday Soc. 1939, 35, 1380–1401. [Google Scholar] [CrossRef]
  42. Uffink, J. Compendium of the foundations of classical statistical physics. Philos. Phys. 2006. [Google Scholar]
  43. Harned, H.S.; Robinson, R.A. A note on the temperature variation of the ionisation constants of weak electrolytes. Trans. Faraday Soc. 1940, 36, 973–978. [Google Scholar] [CrossRef]
  44. Duffin, L.; Wang, X.; Stoesz, J. Detection of Volatile Organic Compounds in Froth Multiphase Systems from Oil Sands Operations Using a Headspace GC-MS Method. Energy Fuels 2017, 31, 11925–11931. [Google Scholar] [CrossRef]
  45. Heintz, A.; Lehmann, J.K.; Wertz, C. Thermodynamic Properties of Mixtures Containing Ionic Liquids. 3. Liquid−Liquid Equilibria of Binary Mixtures of 1-Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide with Propan-1-ol, Butan-1-ol, and Pentan-1-ol. J. Chem. Eng. Data 2003, 472–474. [Google Scholar] [CrossRef]
  46. Voelkel, A. Inverse gas chromatography in characterization of surface. Chemom. Intell. Lab. Syst. 2004, 72, 205–207. [Google Scholar] [CrossRef]
  47. Li, X.; Wang, Q.; Li, L.; Deng, L.; Zhang, Z.; Tian, L. Determination of the thermodynamic parameters of ionic liquid 1-hexyl-3-methylimidazolium chloride by inverse gas chromatography. J. Mol. Liq. 2014, 200, 139–144. [Google Scholar] [CrossRef]
Figure 1. Experimental approach used for the thermodynamic properties determination.
Figure 1. Experimental approach used for the thermodynamic properties determination.
Nanomaterials 11 00709 g001
Figure 2. Linear plots for specific retention volume, lnVg, vs. the inverse of absolute temperature T: nC5, iC5, nC6, nC7, and Toluene under the different asphalting contents (0, 2.56, 9.9336.86, and 53.67 wt%).
Figure 2. Linear plots for specific retention volume, lnVg, vs. the inverse of absolute temperature T: nC5, iC5, nC6, nC7, and Toluene under the different asphalting contents (0, 2.56, 9.9336.86, and 53.67 wt%).
Nanomaterials 11 00709 g002
Figure 3. Linear plots for partition coefficient KR vs. the different asphalting contents for different absolute temperature T (333.2 K, 343.2 K, 353.2 K, 363.2 K, 373.2 K, 383.2 K, and 393.2 K) for different dilution: nC5, iC5, nC6, nC7, and Toluene.
Figure 3. Linear plots for partition coefficient KR vs. the different asphalting contents for different absolute temperature T (333.2 K, 343.2 K, 353.2 K, 363.2 K, 373.2 K, 383.2 K, and 393.2 K) for different dilution: nC5, iC5, nC6, nC7, and Toluene.
Nanomaterials 11 00709 g003
Figure 4. Linear plots for activity coefficients, ln r, vs. the inverse of absolute temperature T: nC5, iC5 nC6, nC7, and Toluene under the different asphalting contents (0, 2.56, 9.93, 21.2, 36.86, 53.67, and 84.2 wt% asphaltene).
Figure 4. Linear plots for activity coefficients, ln r, vs. the inverse of absolute temperature T: nC5, iC5 nC6, nC7, and Toluene under the different asphalting contents (0, 2.56, 9.93, 21.2, 36.86, 53.67, and 84.2 wt% asphaltene).
Nanomaterials 11 00709 g004
Figure 5. Linear plots for enthalpy of solution, ΔH0, vs. the different asphalting contents for five hydrocarbons iC5, nC5, nC6, nC7, and Toluene.
Figure 5. Linear plots for enthalpy of solution, ΔH0, vs. the different asphalting contents for five hydrocarbons iC5, nC5, nC6, nC7, and Toluene.
Nanomaterials 11 00709 g005
Figure 6. Linear plots for specific retention volume, ln V g , vs. the inverse of absolute temperature T for different dilution: nC5, iC5, nC6, nC7, and Toluene under the asphalting contents of 36.86 wt% asphaltene.
Figure 6. Linear plots for specific retention volume, ln V g , vs. the inverse of absolute temperature T for different dilution: nC5, iC5, nC6, nC7, and Toluene under the asphalting contents of 36.86 wt% asphaltene.
Nanomaterials 11 00709 g006
Figure 7. Linear plots for activity coefficients, ln r, vs. the inverse of absolute temperature T: iC5, nC5, nC6, nC7, and Toluene.
Figure 7. Linear plots for activity coefficients, ln r, vs. the inverse of absolute temperature T: iC5, nC5, nC6, nC7, and Toluene.
Nanomaterials 11 00709 g007
Table 1. Properties of isopentane, n-pentane, n-hexane, n-heptane, and toluene using in stationary phases.
Table 1. Properties of isopentane, n-pentane, n-hexane, n-heptane, and toluene using in stationary phases.
iC5n-pentane (nC5)n-hexane (nC6)n-heptane (nC7)Toluene
Boiling point °C27.8–28.236.168.798.42110–111
Molecular FormulaC5H12C5H12C6H14C7H16C6H5CH3
Molecular weight (g/mol)72.1572.1586.18100.2192.14
Vapor Pressure (kPa)
(at 20 °C)
76.99257.9017.605.332.8
Density @ 25 °C(g/cm3)0.6160.6260.6550.6840.865
Table 2. The properties of different Bitumen/Asphaltenes of different S/B.
Table 2. The properties of different Bitumen/Asphaltenes of different S/B.
Samples for Chromosorb CoatingBitumen Molecular Weight, g/molWeight of Coated Chromosorb, gAsphaltene Content in Bitumen, wt%Bitumen Loading, wt%Bitumen Loading WL, g
Bitumen from different S/B1.744211.89.94.40.52
1.7442129.92.60.32
1.744212.19.910.31.25
3.537911.62.67.70.89
3.537911.62.63.10.36
3.537912.22.613.71.67
4035911.606.00.70
4035912020.24
4035912.1014.51.76
Asphaltenes from different S/B1.7111312.853.71.30.17
1.711131353.72.40.31
3.578112.336.91.60.19
3.578111.836.93.00.35
Table 3. Specific retention volume V g 0 comparison between references and the calculated V g 0 .
Table 3. Specific retention volume V g 0 comparison between references and the calculated V g 0 .
Column Temperature, K V g 0   Average   ( mL / g ) V g 0   *
(mL/g)
Column Temperature, K V g 0   Average   ( mL / g ) V g 0   *
(mL/g)
nC5nC7
323.1519.418.7323.15136.8
333.1514.714.3333.1594.192.1
343.1511.411.2343.1565.8
353.159.18.86353.1547.346.7
363.157.47.2363.1535.134.4
373.156.15.86373.1526.526.3
383.154.93.99383.1520.3
393.154.13.36393.1515.815.8
nC6Toluene
323.1551.2 333.15249.9245
333.1537.137.1343.15171.9
343.1527.2 353.15120.3120
353.1520.520.7363.1587.3
363.1516.016.3373.1564.766.2
373.1512.612.6383.1548.849.1
383.1510.0 393.1537.538.1
393.158.08.04
* Data obtained from reference; Standard uncertainties u are u (Vg0) ≤3%, u (T) = 0.5 K.
Table 4. Results of infinite dilution activity coefficients γ for investigated hydrocarbons in bitumen.
Table 4. Results of infinite dilution activity coefficients γ for investigated hydrocarbons in bitumen.
Asphaltene Content, wt%02.569.9321.236.8653.6784.2
Column Temperature, K *
γof iC5
303.151.311.331.37 1.712.46
313.151.261.281.31 1.622.15
323.151.231.251.261.391.571.993.47
333.151.201.221.231.341.521.873.29
343.151.181.191.191.301.481.773.13
353.151.141.151.161.241.451.673.12
363.151.101.111.111.181.411.583.11
373.151.051.061.081.121.341.412.86
383.151.021.021.051.101.291.352.73
393.150.991.001.011.061.251.272.53
γ of nC5
303.151.241.261.29 1.602.16
313.151.211.231.25 1.542.03
323.151.191.201.221.331.501.952.93
333.151.171.181.191.301.471.873.04
343.151.141.161.171.281.441.823.09
353.151.121.141.141.231.421.743.20
363.151.091.101.101.191.371.693.25
373.151.051.071.071.141.301.623.17
383.151.041.031.051.141.261.553.12
393.151.001.011.021.111.221.482.85
γof nC6
333.151.261.281.301.451.612.08
343.151.241.261.281.441.582.01
353.151.221.241.271.411.551.93
363.151.201.221.241.371.511.86
373.151.171.191.211.341.461.77
383.151.161.161.191.331.421.71
393.151.141.151.161.301.391.68
γof nC7
333.151.341.361.401.551.742.31
343.151.331.351.381.541.712.234.07
353.151.311.331.361.531.692.154.27
363.151.291.311.341.501.662.074.37
373.151.271.291.311.471.612.004.28
383.151.251.271.301.451.601.954.27
393.151.231.251.281.441.571.914.12
γ of Toluene
333.150.830.840.840.910.961.141.75
343.150.820.830.830.900.951.121.89
353.150.810.820.820.900.941.101.99
363.150.800.810.810.890.931.092.04
373.150.790.800.800.880.921.062.04
383.150.780.790.790.870.911.052.05
393.150.770.780.780.860.901.042.04
* Standard uncertainties u is u (γ) ≤3%, u (T) = 0.5 K.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuan, J.; Sun, Y.; Jia, Y.; Zhang, Q. Thermodynamic Properties Investigation of Process Volatile Organic Compounds (VOCs) and Its Transport Impact Factor in Oil Sands Management. Nanomaterials 2021, 11, 709. https://doi.org/10.3390/nano11030709

AMA Style

Yuan J, Sun Y, Jia Y, Zhang Q. Thermodynamic Properties Investigation of Process Volatile Organic Compounds (VOCs) and Its Transport Impact Factor in Oil Sands Management. Nanomaterials. 2021; 11(3):709. https://doi.org/10.3390/nano11030709

Chicago/Turabian Style

Yuan, Jing, Yuyong Sun, Yong Jia, and Qianfeng Zhang. 2021. "Thermodynamic Properties Investigation of Process Volatile Organic Compounds (VOCs) and Its Transport Impact Factor in Oil Sands Management" Nanomaterials 11, no. 3: 709. https://doi.org/10.3390/nano11030709

APA Style

Yuan, J., Sun, Y., Jia, Y., & Zhang, Q. (2021). Thermodynamic Properties Investigation of Process Volatile Organic Compounds (VOCs) and Its Transport Impact Factor in Oil Sands Management. Nanomaterials, 11(3), 709. https://doi.org/10.3390/nano11030709

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