Mikhail R. Stukan

Spontaneous Imbibition in Nanopores of Different Roughness and Wettability

The spontaneous imbibition of liquid in nanopores of different roughness is investigated using coarse grain molecular dynamics (MD) simulation. The numerical model is presented and the simplifying assumptions are discussed in detail. The molecular-kinetic theory introduced by Blake is used to describe the effect of dynamic contact angle on fluid imbibition. The capillary roughness is modeled using a random distribution of coarse grained particles forming the wall. The Lucas-Washburn equation is used as a reference for analyzing the imbibition curves obtained by simulation. Due to the statistical nature of MD processing, a comprehensive approach was made to average and smooth the data to accurately define a contact angle. The results are discussed in terms of effective hydrodynamic and static capillary radii and their difference as a function of roughness and wettability.

Review of multi-scale particulate simulation of the rheology of wormlike micellar fluids

We provide a review of our recent work on the development of a multi-scale simulation methodology to calculate the rheology and flow of wormlike micelles. There is a great need for understanding the link between the detailed chemistry of surfactants, forming wormlike micelles, and their macroscopic rheological properties. We show how this link may be explored through particle simulations. First, we calculate the mechanical properties of small units of wormlike micelles from atomistic molecular dynamics simulations. These mechanical properties are subsequently used in a coarse-grained Brownian dynamics model, where the persistence length is the smallest length scale. We show that the non-Newtonian rheology of wormlike micellar systems can be accurately determined from our simple mesoscopic simulation model. Finally, we show that this mesoscopic model can be used to study the flow of wormlike micelles in contraction–expansion geometries.

Influence of system size and solvent flow on the distribution of wormlike micelles in a contraction-expansion geometry

Abstract.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.

Quantifying mineral surface energy by scanning force microscopy

Fundamental understanding of the wettability of carbonate formations can potentially be applied to the development of oil recovery strategies in a complex carbonate reservoir. In the present study, surface energies of representative carbonate samples were evaluated by direct quantitative force measurements, using scanning force microscopy (SFM) at sub-micron scale, to develop a reliable method to predict reservoir wettability. Local adhesion force measurements were conducted on appropriate calcite and dolomite samples and performed in air as well as in the presence of polar and nonpolar fluids. This study demonstrated that, by comparing measurements of adhesion forces between samples of the same mineral in different fluids, it is feasible to determine the surface energy of a given mineral as well as its polar and nonpolar components. The derived values are in agreement with literature. A proof-of-principle protocol has been established to quantify surface energy using SFM-based adhesion measurements. This novel methodology complements the conventional contact angle measurement technique, where surface energy can only be examined at large length scale. The reported methodology has great potential for further optimization into a new standard method for fast and accurate surface energy determination, and hence provides a new tool for reservoir rock wettability characterization.

Dynamic Asphaltene-Stearic Acid Competition at the Oil–Water Interface

Interfacial tension (IFT) is one of the major parameters which govern the fluid flow in oil production and recovery. This paper investigates the interfacial activity of different natural surfactants found in crude oil. The main objective was to better understand the competition between carboxylic acids and asphaltenes on toluene/water interfaces. Dynamic IFT was measured for water-in-oil pendant drops contrary to most studies using oil-in-water drops. Stearic acid (SA) was used as model compound for surface-active carboxylic acids in crude. The influence of concentration of these species on dynamic IFT between model oil and deionized water was examined. The acid concentrations were of realistic values (total acid number 0.1 to 2 mg KOH/g oil) while asphaltene concentrations were low and set between 10 and 100 ppm. In mixtures, the initial surface pressure was entirely determined by the SA content while asphaltenes showed a slow initial diffusion to the interface followed by increased adsorption at longer times. The final surface pressure was higher for asphaltenes compared to SA, but for binaries, the final surface pressure was always lower than the sum of the individuals. At high SA concentration, surface pressures of mixtures were dominated entirely by the SA, although, Langmuir isotherm analysis shows that asphaltenes bind to the interface 200-250 times stronger than SA. The surface area/molecule for both SA and asphaltenes were found to be larger than the values reported in recent literature. Various approaches to dynamic surface adsorption were tested, showing that apparent diffusivity of asphaltenes is very low, in agreement with other works. Hence, the adsorption is apparently under barrier control. A possible hypothesis is that at the initial phase of the experiment and at lower concentration of asphaltenes, the interface is occupied by stearic acid molecules forming a dense layer of hydrocarbon chains that may repel the asphaltenes.

Molecular Dynamics Simulation of Surfactant Flooding in Asphaltenic Oils for EOR Applications

Surfactant flooding has recently re-emerged as a suitable method for Enhanced Oil Recovery (EOR). However, to optimize recovery rates, a better and more fundamental understanding of the recovery process is required at the pore scale. Multi-scale reservoir