Next Article in Journal
Metal Mobilization As An Effect of Anthropogenic Contamination in Groundwater Aquifers in Tutuila, American Samoa
Previous Article in Journal
Features and Control of Submerged Horizontal Vortex in Stepped Dissipation Wells
Previous Article in Special Issue
Sedimentation and Transport of Different Soil Colloids: Effects of Goethite and Humic Acid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precipitant Effects on Aggregates Structure of Asphaltene and Their Implications for Groundwater Remediation

1
Department of Civil Engineering, Ohio University, Athens, OH 45701, USA
2
Oak Ridge Institute for Science and Education (ORISE), U.S. Environmental Protection Agency (EPA), Ada, OK 74820, USA
*
Author to whom correspondence should be addressed.
Water 2020, 12(8), 2116; https://doi.org/10.3390/w12082116
Submission received: 28 May 2020 / Revised: 19 July 2020 / Accepted: 22 July 2020 / Published: 25 July 2020
(This article belongs to the Special Issue Colloid and Pathogen Transport in Groundwater)

Abstract

:
Asphaltenes generally aggregate, then precipitate and deposit on the surfaces of environmental media (soil, sediment, aquifer, and aquitard). Previous studies have recognized the importance of asphaltene aggregates on the wettability of aquifer systems, which has long been regarded as a limiting factor that determines the feasibility and remediation efficiency of sites contaminated by heavy oils. However, the mechanisms/factors associated with precipitant effects on asphaltene aggregates structure, and how the precipitant effects influence the wettability of surfaces remain largely unknown. Here, we observe the particle-by-particle growth of asphaltene aggregates formed at different precipitant concentrations. Our results show that aggregates for all precipitant concentrations are highly polydisperse with self-similar structures. A higher precipitant concentration leads to a more compacted aggregates structure, while precipitant concentration near to onset point results in a less compact structure. The well-known Smoluchowski model is inadequate to describe the structural evolutions of asphaltene aggregates, even for aggregation scenarios induced by a precipitant concentration at the onset point where the Smoluchowski model is expected to explain the aggregate size distribution. It is suggested that aggregates with relative high fractal dimensions observed at high precipitant concentrations can be used to explain the relatively low Stokes settling velocities observed for large asphaltene aggregates. In addition, asphaltene aggregates with high fractal dimensions are likely to have high density of nanoscale roughness which could enhance the hydrophobicity of interfaces when they deposit on the sand surface. Findings obtained from this study advance our current understandings on the fate and transport of heavy oil contaminants in the subsurface environment, which will have important implications for designing and implementing more effective and efficient remediation technologies for contaminated sites.

1. Introduction

Crude oil is a common soil contaminant usually originating from a variety of sources, including petroleum spills and fuel storage tank leakages. Upon its release to the ground surface, crude oil has the propensity to migrate into the subsurface where it can persist as hydrophobic contaminants and occasionally leach to nearby ground and surface waters. Driven by significant adverse environmental impacts and risks, the detection and remediation of crude oil contaminated soils and aquifers have yielded a large amount of research interest in the past decades [1,2,3,4,5]. Solvent extraction is a well-accepted method for the removal of hydrophobic organic contaminants from soil. The Soxhlet extraction technique outlined by the U.S. Environmental Protection Agency (EPA) method 3540C describes the process and solvents used for crude oil contaminated soil remediation in detail [6]. However, solvent extraction can trigger changes in the solvation characteristics of crude oil potentially undermining the efficiency of the remediation process.
Crude oils contain asphaltenes, which are high-molecular-weight polar compounds that are known to be insoluble in light n-alkanes and soluble in aromatic solvents [7,8]. Asphaltenes generally aggregate, and then precipitate and deposit on surfaces of environmental matrixes (e.g., soil, sediment, aquifer, and aquitard). Previous studies have shown significant effects of asphaltenes on the wettability of aquifer systems [9,10,11]. In addition to the classical mechanisms controlled by Derjaguin–Landau–Verwey–Overbeek (DLVO) interactions between asphaltene particles and sand surface, asphaltene aggregation and precipitation are other important mechanisms that can change the wettability of sand, which can have significant impacts on groundwater remediation [9]. Innovations in surfactant flushing and solvent extraction techniques for crude oil contaminant remediation can cause changes in the solvation characteristics of the crude oil, causing asphaltene aggregation, changing the wettability of subsurface system, and hence altering the distribution and recoverability of these contaminants. Although the importance of reversal of wettability caused by the aggregation of asphaltene has been recognized, fundamental understandings on this process remain elusive [1,9,10]. Using an atomic force microscopy (AFM) imaging technique, Zheng et al. [10] have shown that the asphaltene aggregates can grow up to 6000 nm in diameter and 3500 nm in thickness, which are orders of magnitude greater than the dimensions captured by DLVO theory, indicating that DLVO interactions are not the dominant mode of changing the wettability property on a quartz surface. Moreover, the aforementioned study mainly focused on measuring the size (dimension) of asphaltene aggregates deposited on quartz surfaces, without including any details on aggregate structure analysis. Hence, the effects of changing solvation characteristics of crude oil on the formation of asphaltene aggregates with different structures, and associated effects on wettability of sand surfaces, are not clear. The novelty of this study is to provide a fundamental understanding of how precipitant concentration affects asphaltene aggregates structure, which eventually has important implications on wettability of aquifer surfaces.
The study of asphaltene aggregation and precipitation due to changes in pressure, temperature [12,13] and composition [14] have been well documented at different length scales from molecular to colloidal sizes. The focus of this study is on the investigation of colloidal properties of asphaltene aggregates formed at different precipitant concentrations. The aggregates are defined as clusters of asphaltene particles with sizes up to several micrometers in average diameter. The colloidal asphaltene precipitation is observed when the asphaltene aggregates reach a detectable size as a precipitant (i.e., n-heptane) is added to the crude/model oil. Deposition refers to the settling of the asphaltene aggregates onto the solid surface. Existing studies on asphaltene aggregation induced by the addition of precipitant into crude oil have focused on the kinetics of aggregate formation with limited insights on effect of structure [8,15,16,17,18]. It is well accepted that the kinetics of asphaltene aggregation strongly depend on the solubility parameters of the mixture of crude oil, solvent, and precipitant. When precipitants are added at concentrations near the precipitation onset point, asphaltene aggregation can be described by a diffusion-limited aggregation (DLA) process [15]. A slight change in precipitant concentration away from precipitation onset point leads to a crossover between diffusion-limited aggregation (DLA) and reaction-limited aggregation (RLA) [15]. The characteristic time of aggregation decreases exponentially upon increase in the precipitant concentration. However, quantitative results regarding asphaltene aggregate structure are not well understood, with different studies even reporting conflicting results. Yudin et al. and Ashoori et al. studied the kinetics of precipitated, extracted, and re-dissolved asphaltenes in toluene. They found an average value of fractal dimension df = 1.7 ± 0.2 for DLA [18,19] while Huang et al. and Seifried et al. reported much larger values of df for crude oil in the DLA regime [8,15]. Moreover, Hoepfner et al. concluded that no clear trend was observed for the fractal dimensions of asphaltenes aggregates as a function of precipitant (n-heptane) concentration [17]. Therefore, a high-resolution experimental study, which allows a direct analysis of aggregate structure, is strongly needed to better understand the subtle differences in fractal dimensions at different precipitant concentrations.
The overarching objective of this study is to investigate the effects of precipitant concentrations on asphaltene aggregates structure. To determine the sensitivity of aggregates structure to different precipitant concentrations, a wide range of concentrations from close to precipitation onset point and higher were evaluated. A combined technique of in-situ microscopy observation and advanced imaging analysis was used to study aggregate morphological properties under different precipitant concentrations. To reduce the effects of complex fractions typically found in crude oil, the current study used a model oil, prepared by dissolving extracted asphaltene in toluene.

2. Experiments

2.1. Materials

Model oil used in this project was prepared from crude oil from Venezuela with an asphaltene content of 9.0 wt%. The results of elemental analysis of extracted asphaltenes show a total heteroatom (nitrogen, sulfur, and oxygen) content of approximately 10 wt%. Toluene (>99%, Thermo Fisher Scientific Inc., Waltham, MA, USA) and n-heptane (>95%, Thermo Fisher Scientific Inc., Waltham, MA, USA) were used as the solvent and precipitant, respectively.
To extract asphaltenes, the crude oil was mixed with n-heptane in a 1:40 ratio and stirred for 24 h. The precipitated asphaltene solids were separated from the oil/heptane mixture by centrifugation. The collected asphaltene solids were washed with n-heptane for 24 h to remove any non-asphaltene material. The asphaltenes were then dried in an oven at 70 °C. To prepare the model oil, dried asphaltenes were dissolved in toluene to a target concentration of 1 wt% and stirred continuously for 5 days to fully dissolve the asphaltene. The details of asphaltene extraction and model oil preparation were given elsewhere [17].

2.2. Experiments

Prior to the aggregation study, preliminary experiments were conducted to determine the instantaneous onset point of precipitation of asphaltene which is defined as the concentration of n-heptane (vol%) φ that has to be added to the model oil (i.e., asphaltene/toluene) to observe insoluble asphaltenes via optical microscopy in less than 10 min of mixing with the precipitant. In this study, the instantaneous onset concentration for the prepared model oil was measured at 52 vol% n-heptane, which is in agreement with reported onset concentrations of precipitation for asphaltene [18]. To test the sensitivity of asphaltene aggregate structure upon changing the precipitant concentration, three experiments were carried out: (1) adding 52 vol% of n-heptane to initiate the aggregation, where φ φ o n s e t , (2) adding 55 vol% of n-heptane to initiate the aggregation, where φ > φ o n s e t , and (3) adding 57 vol% of n-heptane to initiate the aggregation, where φ φ o n s e t .
The experiments of aggregation were conducted using a liquid cell mounted on a 100× magnification inverted optical microscope (Eclipse Ti-E, Nikon Instruments Inc., Melville, NY, USA). A 3 × 3 large image with a normal field of view was taken at the end of each experiment and stitched together to obtain larger statistic of aggregates for structural analyses of mature aggregates under steady-state conditions. The details of the experimental set up and procedures were described previously [20,21].

2.3. Image Analysis

All analyses were performed in MATLAB (The MathWorks, Inc, Natick, MA, USA) and the key steps are summarized here. First, a global threshold via Ostwald criterion was used to generate binary images from bright-field images (Figure 1). Then, all objects were identified, labeled, and characterized by segmentation, using bwlabel and regionprops functions in MATLAB. Identified objects were also fitted with best-fit ellipses by customized functions. To differentiate individual particles, small clusters, and large clusters, four morphology parameters were used: area, solidity, solidity_ellipse, and eccentricity. More details of the definition and acceptable criteria of the morphological parameters used in this study are illustrated in Figure 2.
The fractal dimensions of aggregates were determined by relating the cluster sizes to the radius of gyration of the clusters [17,20]:
A = R g d f
where A is the area of the cluster and Rg is the radius of gyration. The radius of gyration can be determined by:
R g = ( x x c ) 2 + ( y y c ) 2 N 1
where x c and y c are the xy coordinates of the centroid of each cluster, x and y are the xy coordinates of each pixel in a cluster, and N is the number of pixels in each cluster.

3. Results and Discussion

Figure 1 shows representative structure images of asphaltene aggregates formed at different precipitant concentrations. We found that these aggregates have rugged, irregular structures and show fractal characteristics at all precipitant concentrations. To further elucidate the fractal characteristics of aggregates, the aggregate fractal dimensions (df) were computed by the following scaling law using the area and the radius of gyration: A     R g d f (Figure 3). More details on the goodness of fit are shown in Table 1. For asphaltene aggregates formed at different precipitant concentrations, the sizes of the aggregates Rg were in a wide range of 10–1000 µm, suggesting these aggregates are highly polydisperse. The higher precipitant concentration resulted in a slight increase of initial particle size from 1 to 4 µm, indicating a faster aggregation rate in a short time. The aggregates for all three experiments described above have a fractal scaling region that spanned over two orders of magnitude in length. However, there is a change in the scaling behavior at Rg = 30 µm for all precipitant concentrations, suggesting that the structure of fractal clusters changes as they transition from small to large clusters. For small clusters, regardless of the precipitant concentration, the fractal dimensions of asphaltene aggregates are nearly identical (≈2), indicating a more compact and dense structure. The compact structure of small clusters observed in this study is consistent with previously reported results [17]. The compact structure observed for small clusters can be attributed to cluster restructuring often seen at smaller length scales [17]. Small aggregate clusters were also observed at a precipitant concentration of 57 vol%, which could imply that some asphaltenes are very stable which can only form small fractal aggregates even at high precipitant concentrations. For large clusters, the fractal dimensions were determined to be 1.66 ± 0.17, 1.73 ± 0.17, and 1.83 ± 0.28 for precipitant concentrations at 52 vol%, 55 vol%, and 57 vol%, respectively. A slightly higher fractal dimension with increasing precipitant concentrations indicates that aggregates are less amorphous with more compacted structures. Approaching the precipitation onset point (52 vol% and 55 vol%), the fractal dimensions determined in this study are in good agreement with the findings from the DLA in classic colloids and previous asphaltene aggregation studies [17,22]. Although the mechanisms/factors associated with the observation of compact aggregates structure at higher precipitant concentration are not immediately clear, it is likely due to increased interaggregate attractions (e.g., London dispersion force) or the removal of a steric stabilization barrier [17]. Nevertheless, the fractal dimensions based on image analysis used in this study fall in the range of the df values reported in the literature for asphaltene aggregates [15,17]. It is worth noting that the large stitched images used our analysis are two-dimensional projected areas of actual three-dimensional aggregates. Thus, some of the details of the original aggregate structure were obscured in the projection of the resultant profile. The df determined here is two-dimensional and therefore intrinsically smaller than df determined from three-dimensional based measurement such as light scattering or settling velocity. Thus, the differences between the df measured by image analysis and the df determined based on the latter techniques should be taken into consideration.
We also determined the cluster size distribution for all experiments as depicted in Figure 4. The results show that the data are skewed to the left of the mean for all precipitant concentrations, indicating the aggregate growth is not only dominated by diffusion-limited cluster-cluster aggregation in the later stage, but some non-classical growth mechanisms such as cluster-particle collision and coalescence instead. The results for precipitant concentrations at/near the onset point (experiments 1 and 2) are very similar but different at higher precipitant concentrations (experiment 3). The deviations are most significant for the tail of the distribution, where the data suggest that aggregates formed at the higher precipitant concentration (57 vol%) have a heavier-tailed distribution than those near the precipitant onset point (52 and 55 vol%). The heavy tail distribution at the higher precipitant concentration results in a larger mean aggregate size, which is in good agreement with previous studies [23]. We also compared the cluster size distribution to the Smoluchowski kinetics, that describe the time evolution of an ensemble of particles as they aggregate. The details of Smoluchowski model can be found elsewhere [24]. The results show that Smoluchowski model cannot describe the observations captured in this study, even for aggregation induced by precipitant at the onset point, where the Smoluchowski model is expected to describe cluster size distribution. Specifically, the model underestimates the number of clusters smaller than the mean at precipitant concentrations near the onset point and overestimates the number of clusters that are larger than the mean. These discrepancies indicate that the Smoluchowski model is insufficient to explain the aggregate size distribution. Thus, to fully reveal the mechanisms governing asphaltene aggregation at different precipitant concentrations, both aggregation kinetics and aggregate structure must be examined.

4. Conclusions and Environmental Implications

We observed asphaltene aggregates formation at different precipitant concentrations using in situ observations combined with advanced statistical analyses. Our results show that aggregates for all precipitant concentrations are highly polydisperse with self-similar structures. Higher precipitant concentration leads to a more compacted aggregates structure, while precipitant near to onset point results in a less compact structure. The Smoluchowski model is inadequate to describe the structural evolutions of asphaltene aggregates, even for aggregation scenarios induced by precipitant at the onset point where Smoluchowski model is expected to explain the aggregate size distribution. The fractal dimension values reported here furnish important insights on the permeability and compactness of the asphaltene aggregates, which provides the fundamentals for a more accurate prediction of critical parameters (e.g., particle transport and settling velocities, and wettability of porous media) that determine the feasibility and efficiency of in situ contaminated site remediation. For instance, aggregates with relative high fractal dimensions observed at high precipitant concentrations can be used to explain the relatively low Stokes settling velocities observed for large asphaltene aggregates (Figure 5). This can be attributed to internal flow in the aggregates that does not appreciably reduce the drag force exerted by the surrounding liquid (fluid). The effect of the internal flow on the aggregates depends on permeability and pore size [25]. For aggregates formed at high precipitant concentrations, the growth in the later stage is dominated by particle-cluster collisions, which result in denser (small pore) structures in the inner part of aggregate, and less dense (large pore) structures on the outside. Also, the size and morphology of asphaltene aggregate deposits have non-trivial effects on wettability of aquifer systems (porous media). It is suggested that asphaltene aggregates with different fractal dimensions observed in this study have different nanoscale morphological parameters such as geometry and roughness, which could significantly alter the surface energy and thereby the wettability of grain/sand surfaces (Figure 5) [26]. More specifically, asphaltene aggregates with high fractal dimension are likely to have high density of nanoscale roughness which could enhance the hydrophobicity of interfaces when they deposit on surfaces. These results have important practical implications, especially within the big framework on how to distribute and recover heavy oils in contaminated soils and aquifers. This is especially useful in the application of solvent flushing-based remediation technologies (e.g., solvent blends of acetone-hexane) [27] to dissolve or mobilize heavy soil in the subsurface. The effectiveness of solvent extraction remediation technologies depends on firm contact between the soil surface and the solvents. Typically, the idea behind solvent extraction is the ability of solvents to dissolve oil contaminants for easier mobilization and recovery. However, an undesirable side effect could result if the solvation power of the heavy oil is reduced and the wettability of the surface is altered due to the asphaltene aggregates deposits. The findings from our study shed light on some unresolved issues in remediation technology. For instance, Li et al. examined the effectiveness of acetone-hexane solvent blends for heavy oil contaminated soil remediation and found that among the different chemical fractions in crude oil, the removal efficiency of asphaltene fraction was the lowest [27]. Perhaps, this low removal efficiency could be attributed to asphaltene aggregate formation and deposition, triggered by the application of solvent-blends. The results from our study highlight the need for careful consideration of solvent types/blends and solvent blend ratios used for solvent extraction remediation of crude oil contaminated soils. Therefore, a more thorough understanding of the effects of precipitant concentration on asphaltene aggregate structure and associated effects on its fate and transport in the subsurface environment is needed for designing more effective and efficient heavy oil remediation technologies. Future research efforts should include high resolution studies that allow for the digitization of the settlement trajectory of asphaltene aggregates with different structure in porous media to calculate settling velocity and directly measure nanostructure, surface topography, and contact angle of the asphaltene aggregates and sand interfaces to determine the wettability property.

Author Contributions

Conceptualization, L.W.; Data curation, C.B.H.; Formal analysis, C.B.H.; Methodology, L.W.; Project administration, L.W.; Writing—original draft, C.B.H.; Writing—review & editing, D.W. and L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the following grants to L.W.: the start-up funds from Ohio University and the US National Science Foundation (NSF) #1836905.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Javanbakht, G.; Goual, L. Mobilization and Micellar Solubilization of Napl Contaminants in Aquifer Rocks. J. Contam. Hydrol. 2016, 185–186, 61–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Javanbakht, G.; Goual, L. Impact of Surfactant Structure on NAPL Mobilization and Solubilization in Porous Media. Ind. Eng. Chem. Res. 2016, 55, 11736–11746. [Google Scholar] [CrossRef]
  3. Lowry, E.; Sedghi, M.; Goual, L. Molecular Simulations of Napl Removal from Mineral Surfaces Using Microemulsions and Surfactants. Colloids Surf. A Physicochem. Eng. Asp. 2016, 506, 485–494. [Google Scholar] [CrossRef] [Green Version]
  4. Whittaker, M.; Pollard, S.J.T.; Fallick, T.E. Characterisation of Refractory Wastes at Heavy Oil-Contaminated Sites: A Review of Conventional and Novel Analytical Methods. Environ. Technol. 1995, 16, 1009–1033. [Google Scholar] [CrossRef]
  5. Ghosh, J.; Tick, G.R. A Pore Scale Investigation of Crude Oil Distribution and Removal from Homogeneous Porous Media During Surfactant-Induced Remediation. J. Contam. Hydrol. 2013, 155, 20–30. [Google Scholar] [CrossRef]
  6. EPA SW-846. Test Methods for Evaluating Solid Waste, Physical/Chemical Methods. Final Update, 3rd ed.; U.S. Environmental Protection Agency: Washington, DC, USA, 1996.
  7. Eskin, D.; Mohammadzadeh, O.; Akbarzadeh, K.; Taylor, S.D.; Ratulowski, J. Reservoir Impairment by Asphaltenes: A Critical Review. Can. J. Chem. Eng. 2016, 94, 1202–1217. [Google Scholar] [CrossRef]
  8. Hung, J.; Castillo, J.; Reyes, A. Kinetics of Asphaltene Aggregation in Toluene-Heptane Mixtures Studied by Confocal Microscopy. Energy Fuel 2005, 19, 898–904. [Google Scholar] [CrossRef]
  9. Zheng, J.; Behrens, S.H.; Borkovec, M.; Powers, S.E. Predicting the Wettability of Quartz Surfaces Exposed to Dense Nonaqueous Phase Liquids. Environ. Sci. Technol. 2001, 35, 2207–2213. [Google Scholar] [CrossRef]
  10. Zheng, J.; Shao, J.; Powers, S.E. Asphaltenes from Coal Tar and Creosote: Their Role in Reversing the Wettability of Aquifer Systems. J. Colloid Interface Sci. 2001, 244, 365–371. [Google Scholar] [CrossRef]
  11. Ghadimi, M.; Amani, M.J.; Ghaedi, M.; Malayeri, M.R. Modeling of Formation Damage Due to Asphaltene Deposition in near Wellbore Region Using a Cylindrical Compositional Simulator. J. Pet. Sci. Eng. 2019, 173, 630–639. [Google Scholar] [CrossRef]
  12. Cruz, A.A.; Amaral, M.; Santos, D.; Palma, A.; Franceschi, E.; Borges, G.R.; Coutinho, J.A.; Palácio, J.; Dariva, C. CO2 Influence on Asphaltene Precipitation. J. Supercrit. Fluid. 2019, 143, 24–31. [Google Scholar] [CrossRef]
  13. Dashti, H.; Zanganeh, P.; Kord, S.; Ayatollahi, S.; Amiri, A. Mechanistic Study to Investigate the Effects of Different Gas Injection Scenarios on the Rate of Asphaltene Deposition: An Experimental Approach. Fuel 2020, 262, 116615. [Google Scholar] [CrossRef]
  14. Duran, J.A.; Casas, Y.A.; Xiang, L.; Zhang, L.; Zeng, H.; Yarranton, H.W. Nature of Asphaltene Aggregates. Energy Fuel 2019, 33, 3694–3710. [Google Scholar] [CrossRef]
  15. Seifried, C.M.; Crawshaw, J.; Boek, E.S. Kinetics of Asphaltene Aggregation in Crude Oil Studied by Confocal Laser-Scanning Microscopy. Energy Fuel 2013, 27, 1865–1872. [Google Scholar] [CrossRef]
  16. Rassamdana, H.; Sahimi, M. Asphalt Flocculation and Deposition: Ii. Formation and Growth of Fractal Aggregates. AIChE J. 1996, 42, 3318–3332. [Google Scholar] [CrossRef]
  17. Hoepfner, M.P.; Fávero, C.V.; Haji-Akbari, N.; Fogler, H.S. The Fractal Aggregation of Asphaltenes. Langmuir 2013, 29, 8799–8808. [Google Scholar] [CrossRef]
  18. Yudin, I.K.; Nikolaenko, G.L.; Gorodetskii, E.E.; Kosov, V.I.; Melikyan, V.R.; Markhashov, E.L.; Frot, D.; Briolant, Y. Mechanisms of Asphaltene Aggregation in Toluene-Heptane Mixtures. J. Pet. Sci. Eng. 1998, 20, 297–301. [Google Scholar] [CrossRef]
  19. Ashoori, S.; Abedini, A.; Saboorian, H.; Qorbani Nasheghi, K.; Abedini, R. Mechanisms of Asphaltene Aggregation in Toluene and Heptane Mixtures. J.Jpn. Pet. Inst. 2009, 52, 283–287. [Google Scholar] [CrossRef] [Green Version]
  20. Wu, L.; Ortiz, C.P.; Jerolmack, D.J. Aggregation of Elongated Colloids in Water. Langmuir 2017, 33, 622–629. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, L.; Ortiz, C.P.; Xu, Y.; Willenbring, J.; Jerolmack, D.J. In Situ Liquid Cell Observations of Asbestos Fiber Diffusion in Water. Environ. Sci. Technol. 2015, 49, 13340–13349. [Google Scholar] [CrossRef] [Green Version]
  22. Meakin, P. Formation of Fractal Clusters and Networks by Irreversible Diffusion-Limited Aggregation. Phys. Rev. Lett. 1983, 51, 1119. [Google Scholar] [CrossRef]
  23. Robison Fernlund, J. Image Analysis Method for Determining 3-D Size Distribution of Coarse Aggregates. Bull. Eng. Geol. Environ. 2005, 64, 159–166. [Google Scholar] [CrossRef]
  24. Sholl, D.S.; Skodje, R.T. Late-Stage Coarsening of Adlayers by Dynamic Cluster Coalescence. Physica A 1996, 231, 631–647. [Google Scholar] [CrossRef]
  25. Johnson, C.P.; Li, X.; Logan, B.E. Settling Velocities of Fractal Aggregates. Environ. Sci. Technol. 1996, 30, 1911–1918. [Google Scholar] [CrossRef]
  26. Lord, D.L.; Buckley, J.S. An Afm Study of the Morphological Features That Affect Wetting at Crude Oil-Water-Mica Interfaces. Colloids Surf. A Physicochem. Eng. Asp. 2002, 206, 531–546. [Google Scholar] [CrossRef]
  27. Li, X.; Du, Y.; Wu, G.; Li, Z.; Li, H.; Sui, H. Solvent Extraction for Heavy Crude Oil Removal from Contaminated Soils. Chemosphere 2012, 88, 245–249. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Large-stitched images of asphaltene aggregates at the end of each experiment. (AC) represent bright-field images of aggregates formed by adding 52 vol%, 55 vol%, and 57 vol% heptane to model oil, respectively. (DF) represent the corresponding binary images after applying Ostwald’s global thresholding criterion. The scale bar is 50 μm.
Figure 1. Large-stitched images of asphaltene aggregates at the end of each experiment. (AC) represent bright-field images of aggregates formed by adding 52 vol%, 55 vol%, and 57 vol% heptane to model oil, respectively. (DF) represent the corresponding binary images after applying Ostwald’s global thresholding criterion. The scale bar is 50 μm.
Water 12 02116 g001
Figure 2. Morphology parameters and criteria used for identification of single particle, small and large clusters. (A,E,I) demonstrate the definition of Solidity 2, Eccentricity, and Solidity Ellipse, respectively. Specifically, Solidity 2 is the ratio of filled area to bounding box; E c c e n t r i c i t y = 1 ( b 2 / a 2 ) , where a and b are the semi-major and semi-minor axes of the best fit ellipse respectively; Solidity ellipse is the ratio of filled area to area of the best fit ellipse. (BD,FH,JL) represent plots of Solidity 2, Eccentricity, and Solidity Ellipse as a function of Filled Area (in pixel) of aggregates formed respectively, upon addition of 52 vol% (red squares), 55 vol% (blue squares), and 57 vol% (black squares) n-heptane. Small clusters were identified under the following criterion: Solidity 2 < 0.9 & Solidity 2 > 0.7 & Eccentricity < 0.5. Large clusters were identified under the following criterion: Eccentricity > 0.7 or Solidity_ellipse < 0.9. Single particles were identified under the following criterion: Solidity 2 < 0.9 & Solidity 2 > 0.7 & Eccentricity < 0.5 & Filled Area < 2000 (in pixel).
Figure 2. Morphology parameters and criteria used for identification of single particle, small and large clusters. (A,E,I) demonstrate the definition of Solidity 2, Eccentricity, and Solidity Ellipse, respectively. Specifically, Solidity 2 is the ratio of filled area to bounding box; E c c e n t r i c i t y = 1 ( b 2 / a 2 ) , where a and b are the semi-major and semi-minor axes of the best fit ellipse respectively; Solidity ellipse is the ratio of filled area to area of the best fit ellipse. (BD,FH,JL) represent plots of Solidity 2, Eccentricity, and Solidity Ellipse as a function of Filled Area (in pixel) of aggregates formed respectively, upon addition of 52 vol% (red squares), 55 vol% (blue squares), and 57 vol% (black squares) n-heptane. Small clusters were identified under the following criterion: Solidity 2 < 0.9 & Solidity 2 > 0.7 & Eccentricity < 0.5. Large clusters were identified under the following criterion: Eccentricity > 0.7 or Solidity_ellipse < 0.9. Single particles were identified under the following criterion: Solidity 2 < 0.9 & Solidity 2 > 0.7 & Eccentricity < 0.5 & Filled Area < 2000 (in pixel).
Water 12 02116 g002
Figure 3. Log-log plots of average area A as a function of the radius of gyration R g for asphaltene aggregates formed under different precipitant concentrations. (AC) represents aggregation induced by adding 52 vol%, 55 vol%, and 57 vol% heptane to model oil, respectively. The solid dark blue and black lines represent the best fit line for small clusters and large clusters, respectively. The solid red lines represent the prediction bounds for the fitted functions and the bounds reflect 95% confidence intervals. Single particles were excluded from the regression.
Figure 3. Log-log plots of average area A as a function of the radius of gyration R g for asphaltene aggregates formed under different precipitant concentrations. (AC) represents aggregation induced by adding 52 vol%, 55 vol%, and 57 vol% heptane to model oil, respectively. The solid dark blue and black lines represent the best fit line for small clusters and large clusters, respectively. The solid red lines represent the prediction bounds for the fitted functions and the bounds reflect 95% confidence intervals. Single particles were excluded from the regression.
Water 12 02116 g003
Figure 4. Cluster size distribution of aggregates formed at different precipitant concentrations. (AC) represents aggregation induced by adding 52 vol%, 55 vol%, and 57 vol% heptane into model oil, respectively. Cluster size distributions for each experiment are normalized by their respective mean values to facilitate direct comparison. The red lines are the particle size distribution for Smoluchowski kinetics.
Figure 4. Cluster size distribution of aggregates formed at different precipitant concentrations. (AC) represents aggregation induced by adding 52 vol%, 55 vol%, and 57 vol% heptane into model oil, respectively. Cluster size distributions for each experiment are normalized by their respective mean values to facilitate direct comparison. The red lines are the particle size distribution for Smoluchowski kinetics.
Water 12 02116 g004
Figure 5. Schematic of the proposed asphaltene aggregates structure effects on settling velocity and wettability of aquifer surface, and associated implications for in situ remediation of sites contaminated with heavy oils.
Figure 5. Schematic of the proposed asphaltene aggregates structure effects on settling velocity and wettability of aquifer surface, and associated implications for in situ remediation of sites contaminated with heavy oils.
Water 12 02116 g005
Table 1. Estimation of fractal dimension by least square fitting for power law of R g d f .
Table 1. Estimation of fractal dimension by least square fitting for power law of R g d f .
ExperimentsClusterGoodness of Fit
Best fit dfR2SSEDFEAdjusted R2RMSE
1Large1.66 ± 0.170.965.772030.960.17
Small1.99 ± 0.010.990.009820.990.01
2Large1.73 ± 0.170.974.881630.970.17
Small2.00 ± 0.010.990.006340.990.013
3Large1.83 ± 0.280.9610.871380.960.28
Small2.00 ± 0.010.990.004280.990.012
Note: Experiments 1, 2, and 3 represent asphaltene aggregation induced by adding precipitant (heptane) at 52 vol%, 55 vol% and 57 vol%, respectively. SSE: sum of squared estimate of errors; DFE: degrees of freedom for error; R2: R-squared; and RMSE: root mean square error.

Share and Cite

MDPI and ACS Style

Hammond, C.B.; Wang, D.; Wu, L. Precipitant Effects on Aggregates Structure of Asphaltene and Their Implications for Groundwater Remediation. Water 2020, 12, 2116. https://doi.org/10.3390/w12082116

AMA Style

Hammond CB, Wang D, Wu L. Precipitant Effects on Aggregates Structure of Asphaltene and Their Implications for Groundwater Remediation. Water. 2020; 12(8):2116. https://doi.org/10.3390/w12082116

Chicago/Turabian Style

Hammond, Christian B., Dengjun Wang, and Lei Wu. 2020. "Precipitant Effects on Aggregates Structure of Asphaltene and Their Implications for Groundwater Remediation" Water 12, no. 8: 2116. https://doi.org/10.3390/w12082116

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