Asphaltenes are known as the "cholesterol" of crude oil. They form nanoaggregates, precipitate, adhere to surfaces, block rock pores, and may alter the wetting characteristics of mineral surfaces within the reservoir, hindering oil recovery efficiency. Despite a significant research effort, the structure, aggregation, and deposition of asphaltenes under flowing conditions remain poorly understood. For this reason, we have investigated asphaltenes, their aggregation, and their deposition in capillary flow using multi-scale simulations and experiments. At the colloid scale, we use a hybrid simulation approach. For the solvent, we used the stochastic rotation dynamics (also known as multi-particle collision dynamics) simulation method, which provides both hydrodynamics and Brownian motion. This is coupled to a coarse-grained molecular dynamics (MD) approach for the asphaltene colloids. The colloids interact through a screened Coulomb potential with varying well depth ε. Here, we present results for deposition of asphaltenes in their whole oil in comparison to our previous results on extracted asphaltenes. Imposing a constant pressure drop over the capillary length, we observe that the transient solvent flow rate decreases with an increasing well depth ε. We find that the dimensionless conductivity measured in the experiment can be well-matched by simulation results using an interaction potential well depth of 8kBT. This means that extracted asphaltenes are more sticky than asphaltenes in whole oil. The interactions between the mesoscopic asphaltene colloids can be related to atomistic MD simulations. Molecular structures for the atomistic calculations were obtained using the quantitative molecular representation approach. Using these structures, we calculate the potential of mean force (PMF) between pairs of asphaltene molecules in an explicit solvent. We obtain a reasonable fit using a -1/r2 attraction for the attractive tail of the PMF at intermediate distances. We speculate that this is due to the two-dimensional nature of the asphaltene molecules. Finally, we discuss how we can relate this interaction to the mesoscopic colloid aggregate interaction. We assume that the colloidal aggregates consist of nanoaggregates. Taking into account observed solvent entrainment effects, we deduct the presence of lubrication layers between the nanoaggregates, which leads to a significant screening of the direct asphaltene-asphaltene interactions.@en

Viscoelastic wormlike micelles are formed by surfactants assembling into elongated cylindrical structures. These structures respond to flow by aligning, breaking and reforming. Their response to the complex flow fields encountered in porous media is particularly rich. Here we use a realistic mesoscopic Brownian Dynamics model to investigate the flow of a viscoelastic surfactant (VES) fluid through individual pores idealized as a step expansion-contraction of size around one micron. In a previous study, we assumed the flow field to be Newtonian. Here we extend the work to include the non-Newtonian flow field previously obtained by experiment. The size of the simulations is also increased so that the pore is much larger than the radius of gyration of the micelles. For the non-Newtonian flow field at the higher flow rates in relatively large pores, the density of the micelles becomes markedly non-uniform. In this case, we find that the density in the large, slowly moving entry corner regions is substantially increased.

@en

The aggregation and deposition of colloidal asphaltene in reservoir rock is a significant problem in the oil industry. To obtain a fundamental understanding of this phenomenon, we have studied the deposition and aggregation of colloidal asphaltene in capillary flow by experiment and simulation. For the simulation, we have used the stochastic rotation dynamics (SRD) method, in which the solvent hydrodynamic emerges from the collisions between the solvent particles, while the Brownian motion emerges naturally from the interactions between the colloidal asphaltene particles and the solvent. The asphaltene colloids interact through a screened Coulomb potential. We vary the well depth ε_cc and the flow rate v to obtain Pe^flow ≫ 1 (hydrodynamic interactions dominate) and Re ≪ 1 (Stokes flow). In the simulations, we impose a pressure drop over the capillary length and measure the corresponding solvent flow rate. We observe that the transient solvent flow rate decreases when the asphaltene particles become more "sticky". For a well depth ε_cc = 2k_BT, a monolayer deposits on the capillary wall. With an increasing well depth, the capillary becomes totally blocked. The clogging is transient for ε_cc = 5k_BT, but appears to be permanent for ε_cc = 10-20k_BT. We compare our simulation results with flow experiments in glass capillaries, where we use extracted asphaltenes in toluene, reprecipitated with n-heptane. In the experiments, the dynamics of asphaltene precipitation and deposition were monitored in a slot capillary using optical microscopy under flow conditions similar to those used in the simulation. Maintaining a constant flow rate of 5 μL min^-1, we found that the pressure drop across the capillary first increased slowly, followed by a sharp increase, corresponding to a complete local blockage of the capillary. Doubling the flow rate to 10 μL min^-1, we observe that the initial deposition occurs faster but the deposits are subsequently entrained by the flow. We calculate the change in the dimensionless permeability as a function of time for both experiment and simulation. By matching the experimental and simulation results, we obtain information about (1) the interaction potential well depth for the particular asphaltenes used in the experiments and (2) the flow conditions associated with the asphaltene deposition process.

@en

The aggregation and deposition of colloidal asphaltene in reservoir rock is a significant problem in the oil industry. To obtain a fundamental understanding of this phenomenon, we have studied the deposition and aggregation of colloidal asphaltene in capillary flow by experiment and simulation. For the simulation, we have used the stochastic rotation dynamics (SRD) method, in which the solvent hydrodynamic emerges from the collisions between the solvent particles, while the Brownian motion emerges naturally from the interactions between the colloidal asphaltene particles and the solvent. The asphaltene colloids interact through a screened Coulomb potential. We vary the well depth ε∝ and the flow rate v to obtain Peflow≫1 (hydrodynamic interactions dominate) and Re≪1 (Stokes flow). In the simulations, we impose a pressure drop over the capillary length and measure the corresponding solvent flow rate. We observe that the transient solvent flow rate decreases when the asphaltene particles become more "sticky". For a well depth ε∝=2kBT, a monolayer deposits on the capillary wall. With an increasing well depth, the capillary becomes totally blocked. The clogging is transient for ε∝=5kBT, but appears to be permanent for ε∝=10-20kBT. We compare our simulation results with flow experiments in glass capillaries, where we use extracted asphaltenes in toluene, reprecipitated with n-heptane. In the experiments, the dynamics of asphaltene precipitation and deposition were monitored in a slot capillary using optical microscopy under flow conditions similar to those used in the simulation. Maintaining a constant flow rate of 5μLmin-1, we found that the pressure drop across the capillary first increased slowly, followed by a sharp increase, corresponding to a complete local blockage of the capillary. Doubling the flow rate to 10μLmin-1, we observe that the initial deposition occurs faster but the deposits are subsequently entrained by the flow. We calculate the change in the dimensionless permeability as a function of time for both experiment and simulation. By matching the experimental and simulation results, we obtain information about (1) the interaction potential well depth for the particular asphaltenes used in the experiments and (2) the flow conditions associated with the asphaltene deposition process.@en

Recently there has been a great deal of attention, from researchers both in academia and in industry, focused on the rheological properties of solutions of viscoelastic wormlike micelles formed by surfactants. It is particularly vital to understand the properties of these solutions with regard to their flow in porous media, given their application to the recovery of hydrocarbons from subterranean formations. In this study a realistic mesoscopic Brownian dynamics model has been utilized to investigate the flow of viscoelastic surfactant (VES) fluid through individual pores with sizes of around one micron. In particular the influence of micelle size, pore geometry and flow rate on the ability of worms to pass through the pores was studied. The ways in which these parameters influence the conformational properties of the worms and the spatial distribution of micelles inside the simulation cell was also investigated. Despite the observation that the density and length distributions became non-uniform at higher scission energy, the distribution of breaking and fusion events remained spatially uniform.

@en

There is a great need for understanding the relationship between the structure and chemistry of surfactants forming wormlike micelles, and their macroscopic flow properties. Available macroscopic Rheological Equations of State (REoS) are often inadequate to predict flow behaviour in complex geometries or even to describe the full set of rheological measurements. In this paper, we show how the link between surfactant structure and rheology can be explored through the use of mesoscopic particulate simulations. We describe the development of a realistic Brownian Dynamics model, where the persistence length of the micelle is the smallest length scale. Most of the parameters for the BD model are obtained from atomistic Molecular Dynamics simulations to establish the link to the chemistry of the surfactant. As a preliminary result, we calculate the shear viscosity of a solution of entangled wormlike micelles as a function of increasing shear rate. The results can be approximately mapped on an experimental flow curve for this particular system, using a mean field scaling relationship. Finally, as our motivation in this work is to understand the flow of wormlike micellar fluids in porous media, we give preliminary experimental results for the flow in an expansion-contraction capillary using micro Particle Image Velocimetry (μ-PIV). This demonstrates the complex response of these fluids in non-viscometric flows and the challenges they pose to any predictive simulation model.@en

We carry out a stability analysis of the Bautista-Manero (B-M) constitutive equations for extensional flow of wormlike micelles. We show that all solutions for the steady-state extensional viscosity ηE are unstable when the elongational rates ε exceed some critical value. In some cases the only real solution for the extensional viscosity is negative at intermediate values of the elongational rate. This critical elongational rate is not unfeasibly large, e.g., 250 s-1 for a typical EHAC solution. We note that the extension rates at which enhanced pressure drop is observed experimentally in porous media flow is generally well below the onset of the instability in the B-M model. However, the extension rates presented here are only average values for the porous media in question and at the pore scale a wide range of values will be encountered. For this reason, we require an improved rheological equation of state. We first separate the contributions to the viscosity of the solution and the wormlike micelles. Then, we remove the term containing the retardation time λJ. We now find that the extensional flow behaviour is well defined for both the transient and steady-state cases.@en