Open access peer-reviewed chapter

Experimental and Computational Modeling of Microemulsion Phase Behavior

Written By

Vai Yee Hon and Ismail B.M. Saaid

Submitted: 18 June 2021 Reviewed: 03 November 2021 Published: 22 December 2021

DOI: 10.5772/intechopen.101482

From the Edited Volume

Surfactants and Detergents - Updates and New Insights

Edited by Ashim Kumar Dutta

Chapter metrics overview

256 Chapter Downloads

View Full Metrics

Abstract

The phase behavior of microemulsions formed in a surfactant-brine-oil system for a chemical Enhanced Oil Recovery (EOR) application is complex and depends on a range of parameters. Phase behavior indicates a surfactant solubilization. Phase behavior tests are simple but time-consuming especially when it involves a wide range of surfactant choices at various concentrations. An efficient and insightful microemulsion formulation via computational simulation can complement phase behavior laboratory test. Computational simulation can predict various surfactant properties, including microemulsion phase behavior. Microemulsion phase behavior can be predicted predominantly using Quantitative Structure-Property Relationship (QSPR) model. QSPR models are empirical and limited to simple pure oil system. Its application domain is limited due to the model cannot be extrapolated beyond reference condition. Meanwhile, there are theoretical models based on physical chemistry of microemulsion that can predict microemulsion phase behavior. These models use microemulsion surface tension and torque concepts as well as with solution of bending rigidity of microemulsion interface with relation to surface solubilization and interface energy.

Keywords

  • surfactant
  • microemulsion
  • phase behavior
  • solubilization
  • chemical enhanced oil recovery
  • computational chemistry

1. Introduction

With growing global energy demand and depleting reserves, enhanced oil recovery (EOR) has become more important. Among various EOR processes, chemical EOR has been labeled an expensive process, with overwhelming parameters needed to describe a chemical EOR process and not practical to measure every one of them [1]. The chemical formulations for chemical EOR process consist of single or a combination of alkaline, surfactant and polymer. The traditional chemical EOR processes are polymer flooding, surfactant and alkaline flooding. Over the years, different modes of chemical flood injections were devised. There are the binary mix of alkali–surfactant (AS), surfactant-polymer (SP), alkaline-polymer (AP), and alkaline-surfactant-polymer (ASP) slug [2].

The research and development effort to design a robust chemical EOR formulation tailored to a specific field is challenging and laborious. To design a successful surfactant related chemical EOR formulation, Pope [3] highlighted 8 surfactant selection criteria. The surfactant must produce high solubilization ratio [i.e., low interfacial tension (IFT)] at optimum condition, commercialize at low cost, be feasible to be tailored to specific crude oil, temperature and salinity, comprise of highly branched hydrophobe with low adsorption onto reservoir rock, be insensitive to surfactant concentration above critical micelle concentration (CMC) and have minimal inclination to form liquid crystals, gels, macroemulsions and show rapid coalescence to microemulsion.

Surfactant solubilization can be determined via phase behavior laboratory test. The phase behavior of microemulsions formed in a surfactant-brine-oil system is complex. Phase behavior for a particular microemulsion system has been known to be measured experimentally [4]. The study of surfactant phase behavior consists of the determination of the number, composition, and structure of phases formed by surfactant systems at a given set of conditions (pressure, temperature, and system composition), in observance of Gibbs phase rule [5].

The relationship between solubilization and IFT is described commonly by Healy & Reed [6] correlations and the Chun-Huh [7] equation. Nevertheless, the latter equation is more commonly used. The term “solubilization” was introduced by McBain [8] to describe the increased solubility of a compound associated with the formation of micelles or inverted micelles. The mechanism to enhance solubility varies depending on the surfactant structure, the solvent type and the nature of the solubilized compound. The oil and water solubilization parameters, SPo and SPw respectively, are expressed in SPo = Vo/Vs and SPw = Vw/Vs, where Vo, Vw, and Vs are the volumes of oil, water and surfactant contained within the micellar phase [9]. The solubilization ratio of water and oil phase is measured as either the volume of solubilized water (Vw) or oil (Vo) over volume of surfactant (Vs) in the microemulsion phase. The solubilization ratio of oil (Vo/Vs) increases with the increase in salinity, while the solubilization ratio of water (Vw/Vs) decreases with the increase in salinity [10]. The region in or near where the solubilization of oil and water intersects versus salinity is the optimal salinity.

Advertisement

2. Microemulsion and Winsor phase behavior

Microemulsions are dispersions of oil and water stabilized by surfactant molecules [11, 12, 13]. They can take on many structures such as water droplets in oil, oil droplets in water, sponge like, bicontinuous structures, and lamellar phase [14]. Unlike emulsions, they are thermodynamically stable. This because of the oil-water IFT is low enough (below 10˗2 nM/M) to compensate the dispersion entropy. The interfacial energy is balanced by the dispersion entropy when the dispersion sizes are small enough, which is below 100 Å [15]. Emulsion has much larger dispersion sizes at approximately 1 μm. The different in size between microemulsion and emulsion explains their difference in properties and appearance; however, their fundamental difference is thermodynamic stability [9].

For a given set of conditions (temperature, composition), microemulsion displays well defined structures. The physical properties of microemulsion often undergo an abrupt change over a narrow concentration range [16]. It is generally accepted that this rapid change in the property over the range of concentration is due to the formation of surfactant aggregates or micelle in solution [9]. The nucleation of micelles is spontaneous, forming a structure that can vary between spherical and cylindrical depending on the surfactant molecular structure, solution composition and temperature. Figure 1 depicted the intermicellar equilibrium and the associated phase changes. Micellar structure, S1 is formed when the hydrophilic group of a surfactant are in contact with water while the hydrophobic group are gathered within the interiors of micelles to create small regions from which water is essentially excluded. When aggregates of surfactant form in apolar solvent, it is called inverted micelles, S2. Inverted micelles promote the solubility of water in apolar solvents. Micellar aggregates have a finite mean lifetime, where their structures are mobile, the interfaces are flexible with a rapid exchange of molecules between the neighboring region and the aggregate [9]. Therefore, both S1 and S2 do form isotropic solutions with a bicontinuous and fluctuating structure in a given optimal condition. The microstructure of the sequence of states progressing from S1 to S2 is difficult to study and hence, these systems were large ignored until recent years, when it has been recognized that the isotopic solutions between S1 and S2 are agents which can significantly enhance oil recovery [9].

Figure 1.

Intermicellar equilibrium and associated phase changes [9].

In less frequent cases, the microemulsion structure is made of elongated cylinders, eventually interconnected, or of distorted lamellar (sponge like), as depicted as M1 and M2 in Figure 1. These structures are encountered when the spontaneous curvature of the surfactant layer is small and approaches zero. On the contrary, sponge like structure can be significantly swollen by both oil and water. The G phase in Figure 1 is liquid crystal with lamellar structure [15].

A microemulsion can exist in three types of systems. Winsor’s introduced three types of simple phase diagram that are characterized by the nature of the polyphasic zone at low and moderate concentration. They are Winsor I (WI), Winsor II (WII) and Winsor III (WIII) [17]. Winsor I and II are also commonly known as II˗ and II+. Below a certain salinity, Cseu, the system is WI. Above a certain salinity, Csel, the system is WII. If the salinity is between Cseu and Csel, the system is WIII. In a WIII system, the IFT is lower than WI and WII [4]. WI and WII diagrams shows a characteristic of 2-phase behavior with water or oil microemulsion in equilibrium with the excess phase. WIII diagram shows a bicontinuous microemulsion in equilibrium with both the excess water and oil phases. WIII is considered to have the best probability of recovering additional oil. WII is considered to have the second-best chance to recover additional oil because it shows interaction between the aqueous phase and crude oil. Even though WI demonstrates interaction between the crude oil and the aqueous phase, it is considered to have poorer oil recovery potential than WII.

Winsor’s phase behavior studies in the early 1950s introduced the R ratio of interactions of the adsorbed surfactant at the interface with the neighboring oil and water molecules as a criterion to take into account along with the effects of all formulation variables found in a ternary surfactant-oil-water (SOW) system, i.e., the surfactant head and tail characteristics, the nature of the oil, the aqueous phase salinity, as well as temperature and pressure [18].

In micellar solutions, 3 distinct regions can be identified: an aqueous region, W, an oil or organic region, O, and a surfactant region, C. The variation of the dispersing tendencies at the O and W faces of the C region is expressed qualitatively by Winsor [17] as:

R=ACOAOOACWAWWE1

where ACO and ACW are the interaction of surfactant molecules per unit area at the interface with oil and water respectively, AOO is the interaction between two oil molecules, and AWW is the interaction between two water molecules. Winsor described that an optimal microemulsion (Winsor III) is formed when the microstructure surface is flat, i.e., R = 1. When R < 1, there is a tendency to form oil-in-water emulsion (Winsor II or II+), whereas when R > 1, the tendency is to form water-in-oil emulsion (Winsor I or II˗). Salager [18] presented a diagram (Figure 2) to link the Winsor R ratio with the observe phase type. R < 1, R = 1 and R > 1 correspond to WI (II˗), WIII, and WII (II+) diagrams, respectively. This shows that any formulation change that alters one of the interactions indicated in the ratio can increase or decrease R. When the formulation variation is properly selected to change R from R < 1 to R > 1 or vice versa, it changes the phase behavior from WI to WII or vice versa, with an intermediate WIII three-phase behavior at R = 1 [18].

Figure 2.

Ternary phase diagram, test tube phase behavior and R ratio variations along a 1-dimensional formulation scan [18].

Advertisement

3. Laboratory determination of phase behavior

Phase behavior laboratory test is conducted to assess surfactant performance in generating optimal microemulsion. The test involves the use of graduated cylinders with stopper and an explosion proof oven. A typical range of surfactant concentration and salinity for the phase behavior test are 0.25 to 2.5 wt.% and 5000 to 40,000 ppm respectively. The procedure starts with adding 50% of surfactant solution in brine into the graduated cylinder with 50% of crude oil of interest. The cylinder is then mixed vigorously and aged in the explosion oven for 14 days at reservoir temperature. The phase type and volume of the mixture are observed and measured on day-14 to identify the optimum microemulsion formulation and condition. The complete phase type descriptions are given in Table 1.

Phase typePhase type descriptions
IITwo fluid envelopes exist—A bottom aqueous phase and a top oil phase. No color is visible in the aqueous phase. The crude oil and aqueous phase volumes are equal to the volumes placed in the tube. The surfactant has been driven into the crude oil and no crude oil swelling has taken place (Type II+ phase behavior).
II˗Two fluid envelopes exist—A bottom aqueous phase and an oil phase. The bottom aqueous phase is colored due to surfactant carrying oil into the aqueous phase. The crude volume can be swollen due to the interaction with the surfactant (added and in-situ), but this is not a requirement for this designation.
IIIThree or more fluid envelopes exist—A bottom aqueous phase, one or more middle emulsion phases, and a top crude oil phase. The aqueous phase can be colored with saponified acids (if alkali is presence) from the crude oil; however, this does not necessarily have to be the case.
II+Two fluid envelopes exist—A bottom aqueous phase and a top crude oil phase. The bottom aqueous phase is clear because the surfactant (added and in-situ) resides in the crude oil phase. The crude oil phase is swollen due to surfactant carrying water into the crude oil phase.

Table 1.

Phase behavior type descriptions.

Type III is considered to have the best probability of recovering additional oil. Type II is considered to have the poorest chance to recover additional oil. Type II˗ is considered to have the second-best chance to recover additional oil because it shows interaction between the aqueous phase and crude oil. Even though Type II+ demonstrates interaction between the crude oil and the aqueous phase, it is considered to have poorer oil recovery potential than Type II˗.

Figure 3 illustrates the relationships among phase behavior laboratory test observation, microemulsion structure and Winsor’s R ratio.

Figure 3.

Relationships among phase behavior laboratory test observation, microemulsion structure and Winsor’s R ratio.

Advertisement

4. Computational simulation of phase behavior

Apart from the commonly use of phase behavior laboratory test to assess microemulsion of a particular surfactant for chemical EOR application, it is viable to use computational simulation approach to predict surfactant performance. These computational simulation approaches include the use of various molecular modeling methods such as Monte Carlo (MC) simulations, Molecular Dynamics (MD), Dissipative Particle Dynamics (DPD) and upper scale modeling methods such as Quantitative Structure-Property Relationship (QSPR) approaches. Molecular modeling tools can be used to understand microscopic effects, predict surfactants’ properties and finally to optimize structures and mixtures of surfactants [19]. Molecular modeling tools, in combination with recently developed intermolecular potentials, can provide precise information about microscopic phenomena and lead to accurate estimation of thermophysical properties [20, 21, 22, 23]. Meanwhile, QSPR is an analytical method for breaking down a molecule into a series of numerical values describing its relevant chemical and physical properties. It remains as the focus of many studies aimed at the modeling and prediction of physicochemical and biological properties of molecules [24].

4.1 Quantitative structure-property relationship (QSPR)

QSPR is an approach to relate molecular descriptors with experimental values of properties based on statistical method. Its prediction accuracy is dependent on the size and quality of database and calculation of the relevant descriptors. There are many QSPR-like terms being used for more specific situations, such as Quantitative Structure-Activity Relationships (QSAR), Quantitative Structure-Toxicity Relationships (QSTR), Quantitative Property-Property Relationships (QPPR), Quantitative Sequence-Action Model (QSAM) and Quantitative Structure-Reactivity Relationships (QSRR) [25]. QSPR models have been developed to predict properties of pure surfactants only. Development of QSPR models for mixtures of surfactants is still a challenge [26].

Salager et al. [18] presented the use of Hydrophilic–Lipophilic Deviation (HLD) equations to attain optimum chemical EOR formulation for simple surfactant, oil and water system. The use of HLD concept to predict optimum surfactant formulation is a hybrid approach that combine HLD equations with experiments data, which is a QSPR approach. It was demonstrated that the phase behavior and optimum formulation can be manipulated with four main independent variables: brine salinity, oil alkane carbon number (ACN), surfactant parameter and temperature. However, these models are limited to simple system. Jin et al. [27] predicts the optimum surfactant salinity using HLD equation and measured parameters including the equivalent alkane carbon number (EACN), salinity, surfactant head dependent parameter, Ksurf value and surfactant characteristic curvature, Cc. The work is extended to predict the IFT behavior using the Hydrophilic Lipophilic Difference-Net Average Curvature (HLD-NAC) equation of state.

Moreau et al. [28] applied the QSPR method to predict surfactant optimal salinity based on its correlation with surfactant structures. The QSPR models have been proven in reference conditions but they cannot be extrapolated to other conditions outside the application domain.

Budhathoki et al. [29] use the HLD equation to design the ratio of surfactant mixtures to form optimal microemulsion at reservoir condition. Correct surfactant head dependent parameter, Ksurf and the surfactant temperature dependent parameter, αT values play a crucial role in the accuracy of the HLD method. Both K and αT values can be obtained via a combination use of HLD equation and a series of phase behavior laboratory work for each individual surfactant. This approach can reduce the experimental test matrix to design optimal surfactant mixture, but it is limited to surfactants with known Ksurf and αT values through extensive phase behavior laboratory work.

4.2 Dissipative particle dynamics (DPD)

Many of the interesting phenomena that occur in complex fluids occur at the mesoscale, which is roughly defined as the spatio-temporal scales ranging from 10–104 nm to 1–10ns [30]. These scales of simulation are not feasible using MD simulation. DPD is a coarse-grained type of molecular simulation technique which could reduce the length and timescale of molecular dynamics (MD) simulations, allowing simulation of large and complex system. DPD modeling method can reach large length scales by combining molecule groups into particles or beads, and long-time scales by replacing atomistic forces with soft effective forces [31]. DPD is widely popular simulation approach due to its algorithmic simplicity and huge versatility. By varying the conservative forces between coarse-grained beads, one can readily model complex fluids such as polymers, colloids, surfactants, membranes, vesicles and phase separating fluids [30]. DPD can give insights on spatial organization of surfactants, interesting mechanistic information for films evolution or trends on surface tensions regarding structure of the adsorbed tensioactive molecules at an interface [19]. However, the challenges for a successful DPD simulation is finding robust and general methods for parameterization of the simulation system [32, 33]. This is an active research area with recent approach to apply machine learning for DPD parameterization [34].

A breakthrough approach by Fraaije et al. [35] demonstrated the use of surface torque analysis in simulating surfactant phase behavior with DPD, to determine the optimal brine salinity specifically. Prior to this, QSPR statistical approach have been the only known approach for decades in determining surfactant phase behavior. Buijse et al. [36] used Fraaije et al. approach to design EOR surfactant formulation by optimizing the surfactant head and tail composition as well as the use of co-solvent. Both Fraaije et al. and Buijse et al. work are applied in simple pure oil system. Further discussion on the work of Fraaije et al. is in Section 4.2.1.

Rekvig et al. [14] varies the surfactant chain length and topology to investigate the effect of surfactant structure and composition of the monolayer on the bending rigidity. This work of Rekvig et al. is of particular interest, where the linking of bending rigidity to surfactant structure in predicting the stability of microemulsions is demonstrated. This is important because it is agreed that the bending rigidity is a key parameter in understanding structure and phase behavior of microemulsion [37]. Further details of the approach by Rekvig et al. in determining bending rigidity is discussed further in the next Section 4.2.2.

4.2.1 Surface tension analysis

The work of Fraaije et al. [35] is the only found published work for direct determination of surfactant phase behavior with theoretical foundation based on physical chemistry of microemulsion. The calculation principal is based on an observation by Helfrich [38, 39] that one can calculate the surface torque density (torque) from the first moment of the molecular stress profile, provided the surface is tensionless. A positive torque implies a tendency to bend toward the oil phase and form a microemulsion with oil droplets dispersed in an aqueous phase, while a negative torque a tendency to form a microemulsion that has water droplets dispersed in an oil phase. A surface with zero torque has an indifferent tendency where the system will form the optimal or balanced microemulsion with average zero curvature. Fraaije et al. run the DPD simulation including electrostatics and ion interactions with added salt, surfactants, and oil. The relationship between mechanical coefficients and the stress profile is expressed in Eq. (2) [39, 40]:

Mn=σzzndzE2

where Mn is the stress moment, σ is the stress tensor and z is the coordinate perpendicular to the surface.

In mathematics, a moment of a function is a specific quantitative measure, used in both mechanics and statistics, of the shape of a set of points. For example, for a set of data points representing mass, the 0-th moment is the total mass, the first moment divided by the total mass is the center of mass, and the second moment is the rotational inertia. Similarly, for a set of data points representing probability density, the zeroth moment is the total probability (i.e., one), the first moment is the mean, the second moment is the variance, and the third moment is the skewness. The n-th moment, Mn, of a real-valued continuous function f(x) of a real variable about a value c is given is Eq. (3) below:

Mn=xcnfxdxE3

Fraaije et al. presented a surface tension analysis namely Method of Moments to shows that torque can be calculated based on Eq. (3). This is done by calculating the torque of an microemulsion interface from the first moment, M1, only if the zeroth moment, M0 (the interfacial tension) is zero exactly. Otherwise, both values of the neutral surface position and the interfacial tension has to be known. Meanwhile, in the tensionless limit, the value of zs is inconsequential. The simulation is run across various salinities to find the optimal salinity at which the microemulsion surface torque is zero. Note that there is no clear-cut boundary on when a positive or negative torque transitions in to a small value (around zero). Therefore, the boundaries are somewhat gradual. Fraaije et al. demonstrated that the Method of Moments is principally correct by successfully deriving known empirical coefficients of decade-old QSPR models.

It was noted in Fraaije et al. work that the torque is related to the bending rigidity and spontaneous curvature through Eq. (4):

τkC0E4

where τ is the torque, k is the bending rigidity and C0 is the spontaneous curvature. However, their approach does not allow for direct calculation of bending rigidity or the spontaneous curvature to deduce the stiffness of the interface. Furthermore, they have yet to attempt a full treatment of the compound mixture for application in actual crude oil system.

4.2.2 Interface fluctuation analysis

Microemulsion structure is governed by the elastic constants, the bending rigidity, k and the saddle splay modulus, ks [41]. Bending rigidity characterizes the resistance of the interface toward bending. A low bending rigidity means large thermal undulations and low stability.

Rekvig et al. [14] used DPD to simulate surfactant monolayers on the interface between oil and water to calculate the bending rigidity by analyzing the undulation spectrum. The effect of surfactant density, chain length, adding co-surfactant and linear versus branched surfactant on bending rigidity are investigated. The results show that increase of the monolayer thickness has a larger effect on the bending rigidity than increasing the density of the layer. The bending rigidity also increases with surfactant chain length and is larger for linear than branched surfactants. Bending rigidity decrease linearly with mole fraction of short surfactants. Mixed film has a lower bending rigidity than the corresponding pure film for all mole fractions.

The work by Rekvig et al. [14] is reference with an earlier work by Goetz et al. [42] for lipid bilayer. Goetz et al. was the first to compute bending rigidity in molecular dynamics simulations. Rekvig et al. used a simple model of head, tail, water, and oil beads in DPD to capture the essential properties of ternary systems such as phase separation and adsorption. During the simulation at very low IFT, the interface is not strictly flat and undulatory waves can be observed (Figure 4).

Figure 4.

Undulation wave of microemulsion interface over time (T1 to T3).

These fluctuations of the interface are analyzed to compute bending rigidity. It is firstly done by characterizing the interface based on continuum theory [43]. This is then followed by the adoption of Helfrich [44] free energy of the interface [Eq. (5)] with local displacement from the average position of the interface to enable easier monitoring in simulations.

fH=γ+12k2Hc02+ksKE5

where H is the principal curvature of the surface at M (unique point), K is the gaussian curvature, depending on its sign, the surface will be curved like a sphere or like a saddle-splay, these number give us an idea of the surface’s shape, k is the bending rigidity, ks is the splay modulus which is related to the shape of the interface (K) and c0 is the spontaneous curvature.

Bending rigidity is obtained by analyzing the undulation spectrum of the interface. Based on hypotheses from Rekvig et al. [14] using equipartition principal and fast fourier transform (FFT) application to decomposes the undulation signal into different wave lengths, Eq. (5) is transformed into Eq. (6) below:

hq2=kBTA1γq2+kq4E6

where h is the approximation of local displacement from the average position of the interface, q is the wave vector, 2πλ, λ is the corresponding wavelength, kB is Boltzmann constant (1.38 × 1023 J/K), T is absolute temperature in Kelvin, A is the interface area and γ is the interfacial tension.

Given the spectral intensity, S(q):

Sq=hq2AE7

Combining Eqs. (6) and (7), Eq. (8) is devised:

1Sqq2=γ+kq2kBTE8

Eq. (8) is used to fit the interface undulation spectrum analysis’ results to estimate the bending rigidity, k. In DPD simulation, all units are normalized where kBT is unity.

Published experiment data on bending rigidity values for microemulsion system may be used as a reference for the model predicted bending rigidity values. However, such published data is scarce or not related to monolayers. Majority of the publications are generally focused on theory. Zvelindovsky [45] mentioned that the bending rigidity for a surfactant monolayer between water and oil is usually in the range of 1–20 kBT. Martínez et al. [46] performed MD simulation for SDS surfactant in dodecane and brine system at zero salinity. The associate bending rigidity is 1.3 kBT. SDS is sodium dodecyl sulfate or sodium lauryl sulfate, sometimes written as sodium lauril sulfate. Kegel et al. [47] found that the bending rigidity for SDS surfactant with alcohol in cyclohexane and brine system is around 1 kBT. Binks et al. [48] found a bending rigidity of around 1 kBT for AOT surfactant in nonane and a brine system at optimum salinity and surfactant concentration. AOT is a twin tailed, anionic surfactant with a sulfosuccinate head group stabilized as a salt by a sodium cation. It was also reported that the bending rigidity value depends on the alkane length, where bending rigidity decreases with increase in alkane length.

4.3 Monte Carlo (MC) and molecular dynamics (MD)

Both molecular modeling methods of MD and MC studies have been carried out for academic considerations, mainly on the surfactant aggregation process, instead of for industry application. This is mostly due to the use of atomistic description of the system that requires extensive computational power, which is not practical for industry application. As mention in Section 4.2.2, Goetz et al. [42] provides the first explicit connection between computer simulations with molecular resolution and elastic membrane models based on differential geometry. Goetz’s method demonstrates a relationship between bending rigidity and the IFT.

Advertisement

5. Conclusion

Phase behavior of microemulsion is commonly assessed via laboratory study. These studies are straightforward but laborious especially when it involves a huge range of surfactant choices. Computational simulation is an alternative approach to provide insights into microemulsion phase behavior. There are limited computational simulation studies to predict surfactant phase behavior, whereby the widely used method since the beginning is empirical correlations as in QSPR approach. There are very few non-empirical approaches to predict surfactant phase behavior. These approaches are based on combination of physical chemistry of microemulsion surface tension, torque and bending rigidity concepts.

Advertisement

Acknowledgments

The authors would like to thank PETRONAS and University of Technology PETRONAS (Grant 015MD0-071) for permission and support to publish this book chapter.

Advertisement

Nomenclature

A

interface area

ACO

interaction of surfactant molecules per unit area at the interface with oil

ACW

interaction of surfactant molecules per unit area at the interface with water

AOO

interaction between two oil molecules

AWW

interaction between two water molecules

c

median of the excess stress profile

C

surfactant region

Cc

surfactant characteristic curvature

c0

spontaneous curvature

Cseu

certain salinity 1

Csel

certain salinity 2

fH

free energy of the interface

G

liquid crystal with lamellar microemulsion structure

k

bending rigidity

kB

Boltzmann constant (1.38 × 10˗23 J·K−1)

kb

bond strength

ks

splay modulus

Ksurf

surfactant head dependent parameter

K

gaussian curvature

Mn

n-th moment of a real-valued continuous function f(x) of a real variable about a value c

M1

interconnected elongated cylinder mircoemulsion structure

M2

distorted lamellar elongated cylinder mircoemulsion structure

O

oil or organic region

q

wave vector

S1

micellar structure

S2

inverted micelles structure

SPo

oil solubilization parameter

SPw

water solubilization parameter

T

system temperature

Vo

volume of oil

Vw

volume of water

Vs

volume of surfactant

W

aqueous region

WI

Winsor I or II˗

WII

Winsor II or II+

WIII

Winsor III

z

coordinate

γ

interfacial tension

τ

torque

σ

stress

λ

wavelength

References

  1. 1. Sheng JJ. Modern Chemical Enhanced Oil Recovery: Theory and Practice. United States: Gulf Professional Publishing; 2011
  2. 2. Gbadamosi AO, Junin R, Manan MA, Agi A, Yusuff AS. An Overview of Chemical Enhanced Oil Recovery: Recent Advances and Prospects. Vol. 9. Berlin Heidelberg: Springer; 2019
  3. 3. Pope GA. Overview of chemical EOR center for petroleum and geosystems engineering. In: Casper EOR Workshop. Austin, TX: The University of Texas at Austin; 2007
  4. 4. Sheng JJ. Modern Chemical Enhanced Oil Recovery. Amsterdam, Netherlands: Elsevier Inc.; 2011
  5. 5. Marques EF, Silva BFB. Surfactants, phase behavior. In: Tadros T, editor. Encyclopedia of Colloid and Interface Science. Berlin, Heidelberg: Springer; 2013. pp. 1290-1333
  6. 6. Healy RN, Reed RL. Immiscible microemulsion flooding. SPE Journal. 1977;17:129-139
  7. 7. Huh C. Interfacial tensions and solubilizing ability of a microemulsion phase that coexists with oil and brine. Journal of Colloid and Interface Science. 1979;71(2):408-426. DOI: 10.1016/0021-9797(79)90249-2
  8. 8. McBain JW. Advances in Colloid Science 1. Geneva: Interscience Publishers Inc; 1942
  9. 9. Bourrel M, Schechter RS. Microemulsions and Related Systems. New York: Marcel Dekker; 1998
  10. 10. Elraies KA, Ahmed S. Mechanism of surfactant in microemulsion phase behavior. Journal of Applied Sciences. 2014;14:1049-1054
  11. 11. Overbeek JTG. The first rideal lecture. Microemulsions, a field at the border between lyophobic and lyophilic colloids. Faraday Discussions of the Chemical Society. 1978;65:7
  12. 12. Overbeek JTG. Microemulsions. Proceedings. 1986;B89:61
  13. 13. Gennes D, Taupln C. Microemulsions and the flexibility of oil/water interfaces. Journal of Physical Chemistry. 1982;86(13):2294-2304. DOI: 10.1021/j100210a011
  14. 14. Rekvig L, Hafskjold B, Smit B. Simulating the effect of surfactant structure on bending moduli of monolayers. The Journal of Chemical Physics. 2004;120(10):4897-4905. DOI: 10.1063/1.1645509
  15. 15. Langevin D. Microemulsions—Interfacial aspects. Advances in Colloid and Interface Science. 1991;34:583-595. DOI: 10.1016/0001-8686(91)80059-S
  16. 16. Preston WC. Some correlating principles of detergent action. The Journal of Physical and Colloid Chemistry. 1948;52:84
  17. 17. Winsor PA. Binary and multicomponent solutions of amphiphilic compounds. Solubilization and the formation, structure, and theoretical significance of liquid crystalline solutions. Chemical Reviews. 1968;68(1). DOI: 10.1021/cr60251a001
  18. 18. Salager JL, Forgiarini AM, Márquez L, Manchego L, Bullón J. How to attain an ultralow interfacial tension and a three-phase behavior with a surfactant formulation for enhanced oil recovery: A review. Part 2. Performance improvement trends from Winsor’s premise to currently proposed inter- and intra-molecular mixtu. Journal of Surfactants and Detergents. 2013;16(5):631-663. DOI: 10.1007/s11743-013-1485-x
  19. 19. Creton B, Nieto-Draghi C, Pannacci N. Prediction of surfactants’ properties using multiscale molecular modeling tools: A review. Oil and Gas Science and Technology. 2012;67(6):969-982. DOI: 10.2516/ogst/2012040
  20. 20. Gubbins KE, Moore JD. Molecular modeling of matter: Impact and prospects in engineering. Industrial and Engineering Chemistry Research. 2010;49(7):3026-3046. DOI: 10.1021/ie901909c
  21. 21. Theodorou DN. Progress and outlook in Monte Carlo simulations. Industrial and Engineering Chemistry Research. 2010;49(7):3047-3058. DOI: 10.1021/ie9019006
  22. 22. Maginn EJ, Elliott JR. Historical perspective and current outlook for molecular dynamics as a chemical engineering tool. Industrial and Engineering Chemistry Research. 2010;49(7):3059-3078. DOI: 10.1021/ie901898k
  23. 23. Moeendarbary E, Ng TY, Zangeneh M. Dissipative particle dynamics: Introduction, methodology and complex fluid applications—A review. International Journal of Applied Mechanics. 2009;1(4):737-763. DOI: 10.1142/S1758825109000381
  24. 24. Liu F, Cao C, Cheng B. A quantitative structure-property relationship (QSPR) study of aliphatic alcohols by the method of dividing the molecular structure into substructure. International Journal of Molecular Sciences. 2011;12:2488-2462
  25. 25. Díaz HG. Bioinformatics and quantitative structure-property relationship (QSPR) models. Current Bioinformatics. 2013;8(4):387-389
  26. 26. Muratov EN, Varlamova EV, Artemenko AG, Polishchuk PG, Kuz’Min VE. Existing and developing approaches for QSAR analysis of mixtures. Molecular Informatics. 2012;31(3–4):202-221. DOI: 10.1002/minf.201100129
  27. 27. Jin L et al. Development of a chemical flood simulator based on predictive HLD-NAC equation of state for surfactant. SPE Annual Technical Conference and Exhibition; 2016. p. 20. DOI: 10.2118/181655-ms
  28. 28. Moreau P, Maldonado A, Oukhemanou F, Creton B. Application of quantitative structure-property relationship (QSPR) method for chemical EOR. In: Proceedings of 17th European Symposium on Improved Oil Recovery. OnePetro; 2013. DOI: 10.2118/164091-MS
  29. 29. Budhathoki M, Hsu TP, Lohateeraparp P, Roberts BL, Shiau BJ, Harwell JH. Design of an optimal middle phase microemulsion for ultra high saline brine using hydrophilic lipophilic deviation (HLD) method. Colloids Surfaces A Physicochemical and Engineering Aspects. 2016;488:36-45. DOI: 10.1016/j.colsurfa.2015.09.066
  30. 30. Español P, Warren PB. Perspective: Dissipative particle dynamics. The Journal of Chemical Physics. 2017;146:150901
  31. 31. Gao L, Shillcock J, Lipowsky R. Improved dissipative particle dynamics simulations of lipid bilayers. The Journal of Chemical Physics. 2007;126(1):1-8. DOI: 10.1063/1.2424698
  32. 32. Buijse MA, Tandon K, Jain S, Handgraaf JW, Fraaije J. Surfactant optimization for eor using advanced chemical computational methods. In: SPE Improved Oil Recovery Symposium. Tulsa, Oklahoma: Society of Petroleum Engineers; 2012. p. 12. DOI: 10.2118/154212-MS
  33. 33. Khedr A, Striolo A. DPD parameters estimation for simultaneously simulating water-oil interfaces and aqueous non-ionic surfactants Abeer Khedr and Alberto Striolo* chemical engineering department, University College London, United Kingdom. Journal of Chemical Theory and Computation. 2018;14:6460-6471
  34. 34. Mcdonagh JL, Shkurti A, Bray DJ, Anderson RL, Pyzer-Knapp EO. Utilizing machine learning for efficient parameterization of coarse grained molecular force fields. Journal of Chemical Information and Modeling. 2019;59(10):4278-4288. DOI: 10.1021/acs.jcim.9b00646
  35. 35. Fraaije JGEM, Tandon K, Jain S, Handgraaf JW, Buijse M. Method of moments for computational microemulsion analysis and prediction in tertiary oil recovery. Langmuir. 2013;29(7):2136-2151. DOI: 10.1021/la304505u
  36. 36. Buijse M, Tandon K, Jain S, Jain A, Handgraaf JW, Fraaije JGEM. Accelerated surfactant selection for EOR using computational methods. In: SPE Enhanced Oil Recovery Conference. 2013. pp. 558-567. DOI: 10.2118/165268-ms
  37. 37. Gradzielski M. Bending constants of surfactant layers. Current Opinion in Colloid & Interface Science. 1998;3(5):478-484. DOI: 10.1016/S1359-0294(98)80021-6
  38. 38. W. Helfrich. Amphiphilic Mesophase Made of Defects. Physics of Defects. 1981; 35:716-755
  39. 39. Helfrich W. Elasticity and thermal undulations of fluid films of amphiphilies. Liquids at Interfaces. 1990;48:212-237
  40. 40. Szleifer I, Kramer D, Ben-Shaul A, Gelbart W, Safran S. Molecular theory of curvature elasticity in surfactant films. The Journal of Chemical Physics. 1990;92:6800-6817
  41. 41. Holderer O, Frielinghaus H, Monkenbusch M, Klostermann M, Sottmann T, Richter D. Experimental determination of bending rigidity and saddle splay modulus in bicontinuous microemulsions. Soft Matter. 2013;9(7):2308-2313. DOI: 10.1039/c2sm27449c
  42. 42. Goetz R, Gompper G, Lipowsky R. Mobility and elasticity of self-assembled membranes. Physical Review Letters. 1999;82(1):221-224. DOI: 10.1103/PhysRevLett.82.221
  43. 43. Safran SA. Statistical Thermodynamics of Surfaces, Interfaces, and Membranes. 1st ed. Boulder, CO: Westview Press; 2003
  44. 44. Helfrich W. Elastic properties of lipid bilayers: Theory and possible experiments. Zeitschrift für Naturforschung. 1973;28c:693-703
  45. 45. Zvelindovsky AV. Nanostructured Soft Matter: Experiments, Theory, Simulation and Perspectives. Berlin, Germany: Springer; 2007
  46. 46. Martínez H, Chacón E, Tarazona P, Bresme F. The intrinsic interfacial structure of ionic surfactant monolayers at water-oil and water-vapour interfaces. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2011;467(2131):1939-1958. DOI: 10.1098/rspa.2010.0516
  47. 47. Kegel WK, Bodnàr I, Lekkerkerker HNW. Bending elastic moduli of the surfactant film and properties of a Winsor II microemulsion system. The Journal of Physical Chemistry. 1995;99(10):3272-3281. DOI: 10.1021/j100010a043
  48. 48. Binks BP, Kellay H, Meunier J. Effects of alkane chain length on the bending elasticity constant k of AOT monolayers at the planar oil-water interface. EPL. 1991;16(1):53-58. DOI: 10.1209/0295-5075/16/1/010

Written By

Vai Yee Hon and Ismail B.M. Saaid

Submitted: 18 June 2021 Reviewed: 03 November 2021 Published: 22 December 2021