The Yen−Mullins model, also known as the modified Yen model, specifies the predominant molecular and colloidal structure of asphaltenes in crude oils and laboratory solvents and consists of the following: The most probable asphaltene molecular weight is ∼750 g/mol, with the island molecular architecture dominant. At sufficient concentration, asphaltene molecules form nanoaggregates with an aggregation number less than 10. At higher concentrations, nanoaggregates form clusters again with small aggregation numbers. The Yen−Mullins model is consistent with numerous molecular and colloidal studies employing a broad array of methodologies. Moreover, the Yen−Mullins model provides a foundation for the development of the first asphaltene equation of state for predicting asphaltene gradients in oil reservoirs, the Flory−Huggins− Zuo equation of state (FHZ EoS). In turn, the FHZ EoS has proven applicability in oil reservoirs containing condensates, black oils, and heavy oils. While the development of the Yen−Mullins model was founded on a very large number of studies, it nevertheless remains essential to validate consistency of this model with important new data streams in asphaltene science. In this paper, we review recent advances in asphaltene science that address all critical aspects of the Yen−Mullins model, especially molecular architecture and characteristics of asphaltene nanoaggregates and clusters. Important new studies are shown to be consistent with the Yen−Mullins model. Wide ranging studies with direct interrogation of the Yen−Mullins model include detailed molecular decomposition analyses, optical measurements coupled with molecular orbital calculations, nuclear magnetic resonance (NMR) spectroscopy, centrifugation, direct-current (DC) conductivity, interfacial studies, small-angle neutron scattering (SANS), and small-angle X-ray scattering (SAXS), as well as oilfield studies. In all cases, the Yen−Mullins model is proven to be at least consistent if not valid. In addition, several studies previously viewed as potentially inconsistent with the Yen−Mullins model are now largely resolved. Moreover, oilfield studies using the Yen−Mullins model in the FHZ EoS are greatly improving the understanding of many reservoir concerns, such as reservoir connectivity, heavy oil gradients, tar mat formation, and disequilibrium. The simple yet powerful advances codified in the Yen−Mullins model especially with the FHZ EoS provide a framework for future studies in asphaltene science, petroleum science, and reservoir studies.

The aggregation of asphaltenes has been established for decades by numerous experimental techniques; however, very few studies have been performed on the association free energy and asphaltene aggregation in solvents. The lack of reliable and coherent data on the free energy of association and aggregation size of asphaltene has imposed severe limitations on the thermodynamic modeling of asphaltene phase behavior. Current thermodynamic models either consider asphaltenes as non-associating components or use fitting parameters to characterize the association. In this work, the relations between Gibbs free energy of asphaltene association and asphaltene molecular structure are studied using molecular dynamics (MD). The free energy of association is computed from the potential of mean force profile along the separation distance between the centers of mass of two asphaltene molecules using the umbrella sampling technique in the GROMACS simulation package. The average aggregation number for asphaltene nanoaggregates and clusters is also calculated through MD simulations of 36 asphaltene molecules in toluene and heptane in order to estimate the effects of association free energy and steric repulsion on the aggregation behavior of asphaltenes. Our simulation results confirm that the interactions between aromatic cores of asphaltene molecules are the major driving force for association as the energy of association increases substantially with the number of aromatic rings. Moreover, heteroatoms attached to the aromatic cores have much more influence on the association free energy than to ones attached to the aliphatic chains. The length and number of aliphatic chains do not seem to have a noticeable effect on asphaltene dimerization; however, they have a profound effect on asphaltene aggregation size since steric repulsion can prevent asphaltenes from forming T-shape configurations and therefore decrease the aggregation size of asphaltenes significantly. Our MD simulation results show for the first time asphaltene precipitation in heptane as an explicit solvent, and predict three distinct stages of aggregation (nanoaggregation, clustering, and flocculation) as proposed by the modified Yen model. Finally, the association free energy for asphaltenes in heptane is higher than that in toluene, which is consistent with asphaltene aggregate sizes obtained from MD simulations.

Injection of carbon dioxide in deep saline aquifers is considered as a method of carbon sequestration. The efficiency of this process is dependent on the fluid-fluid and rock-fluid interactions inside the porous media. For instance, the final storage capacity and total amount of capillary-trapped CO2 inside an aquifer are affected by the interfacial tension between the fluids and the contact angle between the fluids and the rock mineral surface. A thorough study of these parameters and their variations with temperature and pressure will provide a better understanding of the carbon sequestration process and thus improve predictions of the sequestration efficiency. In this study, the controversial concept of wettability alteration of quartz surfaces in the presence of supercritical carbon dioxide (sc-CO2) was investigated. A novel apparatus for measuring interfacial tension and contact angle at high temperatures and pressures based on Axisymmetric Drop Shape Analysis with no-Apex (ADSA-NA) method was developed and validated with a simple system. Densities, interfacial tensions, and dynamic contact angles of CO2/water/quartz systems were determined for a wide range of pressures and temperatures relevant to geological sequestration of CO2 in the subcritical and supercritical states. Image analysis was performed with ADSA-NA method that allows the determination of both interfacial tensions and contact angles with high accuracy. The results show that supercritical CO2 alters the wettability of quartz surface toward less water-wet conditions compared to subcritical CO2. Also we observed an increase in the water advancing contact angles with increasing temperature indicating less water-wet quartz surfaces at higher temperatures.

The main objective of this study was to provide novel insights into the mechanism of asphaltene aggregation in toluene/heptane (Heptol) solutions and the effect of alkylphenols on asphaltene dispersion through the integration of advanced experimental and modeling methods. High-resolution transmission electron microscope (HRTEM) images revealed that the onset of asphaltene flocculation occurs near a toluene/heptane volume ratio of 70:30 and that flocculates are well below 1 μm in size. To assess the impact of alkylphenols on asphaltene aggregation, octylphenol (OP) and dodecylphenol (DP) were evaluated by impedance analysis based on their ability to delay the precipitation onset and to reduce the size of nonflocculated asphaltene aggregates in 80:20 toluene/heptane solutions. Although a longer dispersant chain length did not affect the precipitation onset, it reduced the size of the aggregates. Molecular dynamics simulations were then performed to understand the mechanism of interaction between a model asphaltene and OP in heptane. OP molecules saturated the H-bonding sites of asphaltenes and prevented them from interacting laterally between themselves. This explained why OP favored the formation of flocculates with filamentary rather than globular structures, which were clearly observed by HRTEM. Although OP proved to be an effective dispersant, its effectiveness was hindered by its self-association and the fact that it interacted at the periphery of asphaltenes, leaving their aromatic cores uncovered.

Adsorption of crude oil on surfaces is successfully measured with a quartz crystal microbalance with dissipation (QCM-D) under flow conditions. The amounts and thicknesses of adsorbed films are determined with good accuracy using liquid loading corrections. Measurements are performed in solvents where the degree of asphaltene stability is high (toluene) and poor (n-alkanes and heptane/toluene mixtures). In toluene, Langmuir-type adsorption is recorded with saturation film thicknesses of 3−4 nm and limited desorption after rinsing. The size of adsorbing species is also determined at early times where adsorption is diffusion-controlled up to 3 wt % crude oil (or 835 ppm asphaltenes) in toluene. The primary asphaltene species that adsorb on the crystal surface are molecules with a diameter of 0.5−1.6 nm at 139−278 ppm asphaltenes and nanoaggregates with a diameter of 2.6−5.6 nm at 835 ppm asphaltenes in toluene. In n-alkanes (nC7, nC10, and nC12), saturation plateaus are not observed within the experimental time scale. Film thicknesses recorded after 3.5 h are all higher than those in toluene and increase with increasing n-alkane carbon number, mainly because of increasing polarity of aggregates. Atomic force microscopy (AFM) analyses reveal that the size of adsorbed aggregates decreases with n-alkane carbon number. Aging effects show that, with time, the adsorbed films become more rigid in toluene and more viscoelastic in n-alkanes. In heptane/toluene mixtures, a significant increase in the dissipation factor is observed close to the flocculation threshold. Adsorption on various surfaces from toluene shows a high affinity of asphaltenes to hydrophilic surfaces. On the other hand, asphaltenes are almost amorphous in n-alkanes. X-ray photoelectron spectroscopy (XPS) analyses of the adsorbed films confirm these observations.