Mohammad Sedghi

Effect of Asphaltene Structure on Association and Aggregation Using Molecular Dynamics

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.

Asphaltene aggregation and impact of alkylphenols.

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.

Molecular Dynamics Simulations of CO2/Water/Quartz Interfacial Properties: Impact of CO2 Dissolution in Water

The safe trapping of carbon dioxide (CO2) in deep saline aquifers is one of the major concerns of CO2 sequestration. The amount of capillary trapping is dominated by the capillary pressure of water and CO2 inside the reservoir, which in turn is controlled by the interfacial tension (IFT) and the contact angle (CA) of CO2/water/rock systems. The measurement of IFT and CA could be very challenging at reservoir conditions, especially in the presence of toxic cocontaminants. Thus, the ability to accurately predict these interfacial properties at reservoir conditions is very advantageous. Although the majority of existing molecular dynamics (MD) studies of CO2/water/mineral systems were able to capture the trends in IFT and CA variations with pressure and temperature, their predictions often deviated from experimental data, possibly due to erroneous models and/or overlooked chemical reactions. The objective of this study was to improve the MD predictions of IFT and CA of CO2/water/quartz systems at various pressure and temperature conditions by (i) considering the chemical reactions between CO2 and water and (ii) using a new molecular model for α-quartz surface. The results showed that the presence of carbonic acid at the CO2/water interface improved the predictions of IFT, especially at low temperature and high pressure where more CO2 dissolution occurs. On the other hand, the effect on CA was minor. The slight decrease in CA observed across the pressure range investigated could be attributed to an increase in the total number of H-bonds between fluid molecules and quartz surface.

Atomistic Molecular Dynamics Simulations of Crude Oil/Brine Displacement in Calcite Mesopores

Unconventional reservoirs such as hydrocarbon-bearing shale formations and ultratight carbonates generate a large fraction of oil and gas production in North America. The characteristic feature of these reservoirs is their nanoscale porosity that provides significant surface areas between the pore walls and the occupying fluids. To better assess hydrocarbon recovery from these formations, it is crucial to develop an improved insight into the effects of wall-fluid interactions on the interfacial phenomena in these nanoscale confinements. One of the important properties that controls the displacement of fluids inside the pores is the threshold capillary pressure. In this study, we present the results of an integrated series of large-scale molecular dynamics (MD) simulations performed to investigate the effects of wall-fluid interactions on the threshold capillary pressures of oil-water/brine displacements in a calcite nanopore with a square cross section. Fully atomistic models are utilized to represent crude oil, brine, and calcite in order to accommodate electrostatic interactions and H-bonding between the polar molecules and the calcite surface. To this end, we create mixtures of various polar and nonpolar organic molecules to better represent the crude oil. The interfacial tension between oil and water/brine and their contact angle on calcite surface are simulated. We study the effects of oil composition, water salinity, and temperature and pressure conditions on these properties. The threshold capillary pressure values are also obtained from the MD simulations for the calcite nanopore. We then compare the MD results against those generated using the Mayer-Stowe-Princen (MSP) method and explain the differences.

Role of ion-pair interactions on asphaltene stabilization by alkylbenzenesulfonic acids

The dispersion of asphaltenes by dodecylbenzenesulfonic acid (DBSA) has been the subject of several studies in the past. However, it is unclear how these interactions affect the structure of asphaltenes and why asphaltene aggregates are larger in the presence of ionic DBSA. The main goal of this study was to address these points using a combination of high-resolution transmission electron microscopy (HRTEM) and molecular dynamics (MD) simulations. Another objective was to compare ionic DBSA (i.e., dodecylbenzenesulfonate or DBS(-)) to nonionic amphiphiles such as alkylphenols. A striking similarity between dodecylbenzenesulfonate and alkylphenols was that both favored the formation of filamentary rather than globular asphaltene flocculates. However the mechanism by which those filaments formed was very different. Two strong electrostatic interactions between DBSA and asphaltenes were found: (i) those between protonated asphaltenes (i.e., AH(+)) and DBS(-) molecules, which were fifteen times stronger than asphaltene-alkylphenol interactions, and (ii) those between two asphaltene-dispersant pairs (i.e., AH(+)-DBS(-) ion pairs), which did not exist with alkylphenols. These interactions promoted the formation of large and compact asphaltene flocculates, as compared to small and loose ones formed without DBSA. Flocculates with DBSA could further bind to each other through ion-pair interactions. The binding occurred in series (generating long filaments) or in parallel (generating lateral ramifications). However the series configuration was energetically favored due to less steric effects generated by the side aliphatic chains of asphaltenes and DBSA.