Predicting 1H NMR relaxation in Gd3+-aqua using molecular dynamics simulations.

Atomistic molecular dynamics simulations are used to predict 1H NMR T1 relaxation of water from paramagnetic Gd3+ ions in solution at 25 °C. Simulations of the T1 relaxivity dispersion function r1 computed from the Gd3+-1H dipole-dipole autocorrelation function agree within ≃8% of measurements in the range f0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. The simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f0 ≲ 5 MHz, the simulation overestimates r1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r1 at high frequencies in chelated Gd3+ contrast-agents used for clinical MRI.

Role of Wettability on the Adsorption of an Anionic Surfactant on Sandstone Cores.

We investigate the dynamic adsorption of anionic surfactant C14 - 16 alpha olefin sulfonate on Berea sandstone cores with different surface wettability and redox states under high temperature that represents reservoir conditions. Surfactant adsorption levels are determined by analyzing the effluent history data with a dynamic adsorption model assuming Langmuir isotherm. A variety of analyses, including surface chemistry, ionic composition, and chromatography, is performed. It is found that the surfactant breakthrough in the neutral-wet core is delayed more compared to that in the water-wet core because the deposited crude oil components on the rock surface increase the surfactant adsorption via hydrophobic interactions. As the surfactant adsorption is satisfied, the crude oil components are solubilized by surfactant micelles and some of the adsorbed surfactants are released from the rock surface. The released surfactant dissolves in the flowing surfactant solution, thereby resulting in an overshoot of the produced surfactant concentration with respect to the injection value. Furthermore, under water-wet conditions, changing the surface redox potential from an oxidized to a reduced state decreases the surfactant adsorption level by 40%. We find that the decrease in surfactant adsorption is caused not only by removing the iron oxide but also by changing the calcium concentration after the core restoration process (calcite dissolution and ion exchange as a result of using EDTA). Findings from this study suggest that laboratory surfactant adsorption tests need to be conducted by considering the wettability and redox state of the rock surface while recognizing how core restoration methods could significantly alter the ionic composition during surfactant flooding.

Polymer-Coated Nanoparticles for Reversible Emulsification and Recovery of Heavy Oil.

Heavy crude oil has poor solubility and a high density, making recovery and transport much more difficult and expensive than for light crude oil. Emulsifiers can be used to produce low viscosity oil-in-water emulsions for recovery and transport, but subsequent demulsification can be challenging. Here, we develop and implement interfacially active, pH-responsive polymer-coated nanoparticles (PNPs) to reversibly stabilize, recover, and break oil/water emulsions through variation of solution pH. Silica particles with poly(2-(dimethylamino)ethyl methacrylate) (DMAEMA) chains covalently grafted to the surface are prepared although a reversible addition fragmentation chain transfer grafting-through technique. The resulting DMAEMA PNPs can stabilize emulsions of high viscosity Canadian heavy oil at PNP concentrations as low as 0.1 wt % and at neutral pH. The performance of the DMAEMA PNPs exceeds that of DMAEMA homopolymer additives, which we attribute to the larger size and irreversible adsorption of DMAEMA PNPs to the oil/water interface. After recovery, the emulsion can be destabilized by the addition of acid to reduce pH, resulting in separation and settling of the heavy oil from the aqueous phase. Recovery of more than 10 wt % of the crude heavy oil-in-place is achieved by flooding with aqueous solution of 0.1 wt % DMAEMA PNPs without any additional surfactant or chemical. This work demonstrates the applicability of PNPs as surface active materials for enhanced oil recovery processes and for heavy oil transport.

Destabilization, Propagation, and Generation of Surfactant-Stabilized Foam during Crude Oil Displacement in Heterogeneous Model Porous Media.

Foam flooding in porous media is of increasing interest due to its numerous applications such as enhanced oil recovery, aquifer remediation, and hydraulic fracturing. However, the mechanisms of oil-foam interactions have yet to be fully understood at the pore level. Here, we present three characteristic zones identified in experiments involving the displacement of crude oil from model porous media via surfactant-stabilized foam, and we describe a series of pore-level dynamics in these zones which were not observed in experiments involving paraffin oil. In the displacement front zone, foam coalesces upon initial contact with crude oil, which is known to destabilize the liquid lamellae of the foam. Directly upstream, a transition zone occurs where surface wettability is altered from oil-wet to water-wet. After this transition takes place, a strong foam bank zone exists where foam is generated within the porous media. We visualized each zone using a microfluidic platform, and we discuss the unique physicochemical phenomena that define each zone. In our analysis, we also provide an updated mechanistic understanding of the "smart rheology" of foam which builds upon simple "phase separation" observations in the literature.

Modeling micelle formation and interfacial properties with iSAFT classical density functional theory.

Surfactants reduce the interfacial tension between phases, making them an important additive in a number of industrial and commercial applications from enhanced oil recovery to personal care products (e.g., shampoo and detergents). To help obtain a better understanding of the dependence of surfactantproperties on molecular structure, a classical density functional theory, also known as interfacial statistical associating fluid theory, has been applied to study the effects of surfactant architecture on micelle formation and interfacial properties for model nonionic surfactant/water/oil systems. In this approach, hydrogen bonding is explicitly included. To minimize the free energy, the system minimizes interactions between hydrophobic components and hydrophilic components with water molecules hydrating the surfactant head group. The theory predicts micellar structure, effects of surfactant architecture on critical micelle concentration, aggregation number, and interfacial tension isotherm of surfactant/water systems in qualitative agreement with experimental data. Furthermore, this model is applied to study swollen micelles and reverse swollen micelles that are necessary to understand the formation of a middle-phase microemulsion.

Static Adsorption of an Ethoxylated Nonionic Surfactant on Carbonate Minerals.

The static adsorption of C12-14E22, which is a highly ethoxylated nonionic surfactant, was studied on different minerals using high-performance liquid chromatography (HPLC) combined with an evaporative light scattering detector (ELSD). Of particular interest is the surfactant adsorption in the presence of CO2 because it can be used for foam flooding in enhanced oil recovery applications. The effects of the mineral type, impurities, salinity, and temperature were investigated. The adsorption of C12-14E22 on pure calcite was as low as 0.01 mg/m2 but higher on dolomite depending on the silica and clay content in the mineral. The adsorption remained unchanged when the experiments were performed using a brine solution or 0.101 MPa (1 atm) CO2, which indicates that electrostatic force is not the governing factor that drives the adsorption. The adsorption of C12-14E22 on silica may be due to hydrogen bonding between the oxygen in the ethoxy groups of the surfactant and the hydroxyl groups on the mineral surface. Additionally, thermal decomposition of the surfactant was severe at 80 °C but can be inhibited by operating in a reducing environment. Under reducing conditions, adsorption of C12-14E22 increased at higher temperatures.

Role of Gas Type on Foam Transport in Porous Media.

We present the results of an experimental investigation of the effect of gas type and composition on foam transport in porous media. Steady-state foam strengths with respect to three cases of distinct gases and two cases containing binary mixtures of these gases were compared. The effects of gas solubility, the stability of lamellae, and the gas diffusion rate across the lamellae were examined. Our experimental results showed that the steady-state foam strength is inversely correlated with the gas permeability across a liquid lamella, a parameter that characterizes the rate of mass transport. The results are also in good agreement with existing observations that the foam strength for a mixture of gases is correlated with the less soluble component. Three hypotheses with different predictions of the underlying mechanism that explain the role of gas type and composition in foam strength are discussed in detail.

Improved methylene blue two-phase titration method for determining cationic surfactant concentration in high-salinity brine.

The methylene blue (MB) two-phase titration method is a rapid and efficient method for determining the concentrations of anionic surfactants. The point at which the aqueous and chloroform phases appear equally blue is called Epton's end point. However, many inorganic anions, e.g., Cl(-), NO3(-), Br(-), and I(-), can form ion pairs with MB(+) and interfere with Epton's end point, resulting in the failure of the MB two-phase titration in high-salinity brine. Here we present a method to extend the MB two-phase titration method for determining the concentration of various cationic surfactants in both deionized water and high-salinity brine (22% total dissolved solid). A colorless end point, at which the blue color is completely transferred from the aqueous phase to the chloroform phase, is proposed as titration end point. Light absorbance at the characteristic wavelength of MB is measured using a spectrophotometer. When the absorbance falls below a threshold value of 0.04, the aqueous phase is considered colorless, indicating that the end point has been reached. By using this improved method, the overall error for the titration of a permanent cationic surfactant, e.g., dodecyltrimethylammonium bromide, in deionized (DI) water and high-salinity brine is 1.274% and 1.322% with limits of detection (LOD) of 0.149 and 0.215 mM, respectively. Compared to the traditional acid-base titration method, the error of this improved method for a switchable cationic surfactant, e.g., tertiary amine surfactant (Ethomeen C12), is 2.22% in DI water and 0.106% with LOD of 0.369 and 0.439 mM, respectively.

Effects of velocity on N2 and CO2 foam flow with in-situ capillary pressure measurements in a high-permeability homogeneous sandpack.

The effects of velocity and gas type on foam flow through porous media have yet to be completely elucidated. Pressure drop and capillary pressure measurements were made at ambient conditions during a series of foam quality scan experiments in a homogenous sandpack while foam texture was simultaneously visualized. New insights into foam-flow behavior in porous media were discovered. Previously accepted "limiting" capillary pressure theory is challenged by the findings in this work, and the "limiting" terminology is replaced with the word "plateau" to reflect these novel observations. Plateau capillary pressure [Formula: see text] and transition foam quality were found to increase with velocity. Transition foam quality was found to depend mostly on liquid velocity rather than gas velocity and is physically linked to foam type (continuous vs. discontinuous) and texture (fine vs. coarse). Distinct rheological behaviors also arose in the low- and high-quality foam regimes as a function of velocity. Foam flow was found to be strongly shear thinning in the low-quality regime where foam texture was fine and discontinuous. In the high-quality regime, the rheology was weakly shear thinning to Newtonian for coarsely textured foam and continuous-gas flow respectively. When all other variables were held constant, at ambient conditions, CO2 foam was found to be weaker with also lower capillary pressures than N2 foam and the differences in gas solubility is a likely explanation.

Effect of Nanoconfinement on NMR Relaxation of Heptane in Kerogen from Molecular Simulations and Measurements.

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nanoconfinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. We use atomistic MD simulations to investigate the effects of nanoconfinement on the 1H NMR relaxation times T1 and T2 of heptane in kerogen. In the case of T1, we discover the important role of confinement in reducing T1 by ∼3 orders of magnitude from that of bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of T2, we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces T2 by ∼5 orders of magnitude from the value for bulk heptane. We use the simulated T2 to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nanopores in kerogen, without additional experimental data.

A systematic approach to alkaline-surfactant-foam flooding of heavy oil: microfluidic assessment with a novel phase-behavior viscosity map.

The apparent viscosity of viscous heavy oil emulsions in water can be less than that of the bulk oil. Microfluidic flooding experiments were conducted to evaluate how alkali-surfactant-foam enhanced oil recovery (ASF EOR) of heavy oil is affected by emulsion formation. A novel phase-behavior viscosity map-a plot of added salinity vs. soap fraction combining phase behavior and bulk apparent viscosity information-is proposed as a rapid and convenient method for identifying suitable injection compositions. The characteristic soap fraction, [Formula: see text], is shown to be an effective benchmark for relating information from the phase-viscosity map to expected ASF flood test performance in micromodels. Characteristically more hydrophilic cases were found to be favorable for recovering oil, despite greater interfacial tensions, due to wettability alteration towards water-wet conditions and the formation of low apparent-viscosity oil-in-water (O/W) macroemulsions. Wettability alteration and bubble-oil pinch-off were identified as contributing mechanisms to the formation of these macroemulsions. Conversely, characteristically less hydrophilic cases were accompanied by a large increase in apparent viscosity due to the formation of water-in-oil (W/O) macroemulsions.

NMR 1H-1H Dipole Relaxation in Fluids: Relaxation of Individual 1H-1H Pairs versus Relaxation of Molecular Modes.

The intramolecular 1H NMR dipole-dipole relaxation of molecular fluids has traditionally been interpreted within the Bloembergen-Purcell-Pound (BPP) theory of NMR intramolecular relaxation. The BPP theory draws upon Debye's theory for describing the rotational diffusion of the 1H-1H pair and predicts a monoexponential decay of the 1H-1H dipole-dipole autocorrelation function between distinct spin pairs. Using molecular dynamics (MD) simulations, we show that for both n-heptane and water this is not the case. In particular, the autocorrelation function of individual 1H-1H intramolecular pairs itself evinces a rich stretched-exponential behavior, implying a distribution in rotational correlation times. However, for the high-symmetry molecule neopentane, the individual 1H-1H intramolecular pairs do conform to the BPP description, suggesting an important role of molecular symmetry in aiding agreement with the BPP model. The intermolecular autocorrelation functions for n-heptane, water, and neopentane also do not admit a monoexponential behavior of individual 1H-1H intermolecular pairs at distinct initial separations. We suggest expanding the autocorrelation function in terms of modes, provisionally termed molecular modes, that do have an exponential relaxation behavior. With care, the resulting Fredholm integral equation of the first kind can be inverted to recover the probability distribution of the molecular modes. The advantages and limitations of this approach are noted.

Critical Role of Confinement in the NMR Surface Relaxation and Diffusion of n-Heptane in a Polymer Matrix Revealed by MD Simulations.

The mechanism behind the NMR surface-relaxation times (T1S,2S) and the large T1S/T2S ratio of light hydrocarbons confined in the nanopores of kerogen remains poorly understood and consequently has engendered much debate. Toward bringing a molecular-scale resolution to this problem, we present molecular dynamics (MD) simulations of 1H NMR relaxation and diffusion of n-heptane in a polymer matrix. The high-viscosity polymer is a model for kerogen and bitumen that provides an organic "surface" for heptane. Diffusion of n-heptane shows a power-law dependence on the concentration of n-heptane (ϕC7) in the polymer matrix, consistent with Archie's model of tortuosity. We calculate the autocorrelation function G(t) for 1H-1H dipole-dipole interactions of n-heptane in the polymer matrix and use this to generate the NMR frequency (f0) dependence of T1S,2S as a function of ϕC7. We find that increasing molecular confinement increases the correlation time, which decreases the surface-relaxation times for n-heptane in the polymer matrix. For weak confinement (ϕC7 > 50 vol %), we find that T1S/T2S ≃ 1. Under strong confinement (ϕC7 ≲ 50 vol %), we find that T1S/T2S ≳ 4 increases with decreasing ϕC7 and that the dispersion relation T1S ∝ f0 is consistent with previously reported measurements of polydisperse polymers and bitumen. Such frequency dependence in bitumen has been previously attributed to paramagnetism; instead, our studies suggests that 1H-1H dipole-dipole interactions enhanced by organic nanopore confinement dominate the NMR response in saturated organic-rich shales.

Elucidating the 1H NMR Relaxation Mechanism in Polydisperse Polymers and Bitumen Using Measurements, MD Simulations, and Models.

The mechanism behind the 1H nuclear magnetic resonance (NMR) frequency dependence of T1 and the viscosity dependence of T2 for polydisperse polymers and bitumen remains elusive. We elucidate the matter through NMR relaxation measurements of polydisperse polymers over an extended range of frequencies (f0 = 0.01-400 MHz) and viscosities (η = 385-102 000 cP) using T1 and T2 in static fields, T1 field-cycling relaxometry, and T1ρ in the rotating frame. We account for the anomalous behavior of the log-mean relaxation times T1LM ∝ f0 and T2LM ∝ (η/T)-1/2 with a phenomenological model of 1H-1H dipole-dipole relaxation, which includes a distribution in molecular correlation times and internal motions of the nonrigid polymer branches. We show that the model also accounts for the anomalous T1LM and T2LM in previously reported bitumen measurements. We find that molecular dynamics (MD) simulations of the T1 ∝ f0 dispersion and T2 of similar polymers simulated over a range of viscosities (η = 1-1000 cP) are in good agreement with measurements and the model. The T1 ∝ f0 dispersion at high viscosities agrees with previously reported MD simulations of heptane confined in a polymer matrix, which suggests a common NMR relaxation mechanism between viscous polydisperse fluids and fluids under nanoconfinement, without the need to invoke paramagnetism.

Evaluating physicochemical properties of crude oil as indicators of low-salinity-induced wettability alteration in carbonate minerals.

The injection of low-salinity brine enhances oil recovery by altering the mineral wettability in carbonate reservoirs. However, the reported effectiveness of low-salinity water varies significantly in the literature, and the underlying mechanism of wettability alteration is controversial. In this work, we investigate the relationships between characteristics of crude oils and the oils' response to low-salinity water in a spontaneous imbibition test, aiming (1) to identify suitable indicators of the effectiveness of low-salinity water and (2) to evaluate possible mechanisms of low-salinity-induced wettability alteration, including rock/oil charge repulsion and microdispersion formation. Seven oils are tested by spontaneous imbibition and fully characterized in terms of their acidity, zeta potential, interfacial tension, microdispersion propensity, water-soluble organics content and saturate-aromatic-resin-asphaltene fractionation. For the first time, the effectiveness of low-salinity water is found to positively correlate with the oil interfacial tension in low-salinity water. Oils with higher interfacial activity are found to respond more positively to low-salinity water. Moreover, cryogenic transmission electron microscopy images suggest that microdispersion is essentially macroemulsion, and its formation is an effective indicator - but not the root cause - of wettability alteration. The repulsive zeta potential for the rock and the oil in low-salinity water is found to be an insufficient condition for wettability alteration in carbonate minerals.

Effect of salinity, Mg2+ and SO42- on "smart water"-induced carbonate wettability alteration in a model oil system.

HYPOTHESIS: We present a systematic study of the "smart water" induced wettability alteration. This process is believed to be greatly affected by the brine salinity and the presence of Mg2+ and SO42- in the brine. EXPERIMENTS AND MODELLING: To characterize the wettability alteration, we perform spontaneous imbibition measurement using Indiana limestone cores and a model oil with added naphthenic acid. Both single-electrolyte-based and seawater-based "smart water" are tested to investigate the effect of Mg2+, SO42- and salinity on wettability alteration. Rock/brine and oil/brine zeta potentials are measured, and the electrostatic component of disjoining pressure is calculated to understand the role of electrostatics in the wettability alteration. The surface concentration of charged species on the limestone surface is analyzed based on a natural carbonate surface complexation model (SCM). FINDINGS: Both the reduction of Na+ and addition of SO42- are found to contribute to wettability alteration. Mg2+ is found to be unfavorable for wettability alteration. Ca2+ is believed to facilitate SO42- with wettability alteration based on the comparison between the single-electrolyte-based and seawater-based brines. The reduction of the Na+ surface complexation (>CaOH⋯Na+0.25) in low salinity brines is believed to be a critical mechanism responsible for wettability alteration based on the SCM calculations.

NMR measurement of bitumen at different temperatures.

Heavy oil (bitumen) is characterized by its high viscosity and density, which is a major obstacle to both well logging and recovery. Due to the lost information of T2 relaxation time shorter than echo spacing (TE) and interference of water signal, estimation of heavy oil properties from NMR T2 measurements is usually problematic. In this work, a new method has been developed to overcome the echo spacing restriction of NMR spectrometer during the application to heavy oil (bitumen). A FID measurement supplemented the start of CPMG. Constrained by its initial magnetization (M0) estimated from the FID and assuming log normal distribution for bitumen, the corrected T2 relaxation time of bitumen sample can be obtained from the interpretation of CPMG data. This new method successfully overcomes the TE restriction of the NMR spectrometer and is nearly independent on the TE applied in the measurement. This method was applied to the measurement at elevated temperatures (8-90 degrees C). Due to the significant signal-loss within the dead time of FID, the directly extrapolated M0 of bitumen at relatively lower temperatures (<60 degrees C) was found to be underestimated. However, resulting from the remarkably lowered viscosity, the extrapolated M0 of bitumen at over 60 degrees C can be reasonably assumed to be the real value. In this manner, based on the extrapolation at higher temperatures (> or = 60 degrees C), the M0 value of bitumen at lower temperatures (<60 degrees C) can be corrected by Curie's Law. Consequently, some important petrophysical properties of bitumen, such as hydrogen index (HI), fluid content and viscosity were evaluated by using corrected T2.

Paramagnetic relaxation in sandstones: distinguishing T1 and T2 dependence on surface relaxation, internal gradients and dependence on echo spacing.

This work provides a generalized theory of proton relaxation in inhomogeneous magnetic fields. Three asymptotic regimes of relaxation are identified depending on the shortest characteristic time scale. Numerical simulations illustrate that the relaxation characteristics in the regimes such as the T(1)/T(2) ratio and echo spacing dependence are determined by the time scales. The theoretical interpretation is validated for fluid relaxation in porous media in which field inhomogeneity is induced due to susceptibility contrast of fluids and paramagnetic sites on pore surfaces. From a set of measurements on model porous media, we conclude that when the sites are small enough, no dependence on echo spacing is observed with conventional low-field NMR spectrometers. Echo spacing dependence is observed when the paramagnetic materials become large enough or form a 'shell' around each grain such that the length scale of the region of induced magnetic gradients is large compared to the diffusion length during the time of the echo spacing. The theory can aid in interpretation of diffusion measurements in porous media as well as imaging experiments in presence of contrast agents used in MRI.

Microfluidic Devices for Characterizing Pore-scale Event Processes in Porous Media for Oil Recovery Applications.

Microfluidic devices are versatile tools for studying transport processes at a microscopic scale. A demand exists for microfluidic devices that are resistant to low molecular-weight oil components, unlike traditional polydimethylsiloxane (PDMS) devices. Here, we demonstrate a facile method for making a device with this property, and we use the product of this protocol for examining the pore-scale mechanisms by which foam recovers crude oil. A pattern is first designed using computer-aided design (CAD) software and printed on a transparency with a high-resolution printer. This pattern is then transferred to a photoresist via a lithography procedure. PDMS is cast on the pattern, cured in an oven, and removed to obtain a mold. A thiol-ene crosslinking polymer, commonly used as an optical adhesive (OA), is then poured onto the mold and cured under UV light. The PDMS mold is peeled away from the optical adhesive cast. A glass substrate is then prepared, and the two halves of the device are bonded together. Optical adhesive-based devices are more robust than traditional PDMS microfluidic devices. The epoxy structure is resistant to swelling by many organic solvents, which opens new possibilities for experiments involving light organic liquids. Additionally, the surface wettability behavior of these devices is more stable than that of PDMS. The construction of optical adhesive microfluidic devices is simple, yet requires incrementally more effort than the making of PDMS-based devices. Also, though optical adhesive devices are stable in organic liquids, they may exhibit reduced bond-strength after a long time. Optical adhesive microfluidic devices can be made in geometries that act as 2-D micromodels for porous media. These devices are applied in the study of oil displacement to improve our understanding of the pore-scale mechanisms involved in enhanced oil recovery and aquifer remediation.

Characterizing adsorption of associating surfactants on carbonates surfaces.

HYPOTHESIS: The adsorption of anionic surfactants onto positively charged carbonate minerals is typically high due to electrostatic interactions. By blending anionic surfactants with cationic or zwitterionic surfactants, which naturally form surfactant complexes, surfactant adsorption is expected to be influenced by a competition between surfactant complexes and surfactant-surface interactions. EXPERIMENTS: The adsorption behavior of surfactant blends known to form complexes was investigated. The surfactants probed include an anionic C15-18 internal olefin sulfonate (IOS), a zwitterionic lauryl betaine (LB), and an anionic C13-alcohol polyethylene glycol ether carboxylic acid (L38). An analytical method based on high-performance liquid chromatography evaporative light scattering detector (HPLC-ELSD) was developed to measure three individual surfactant concentrations from a blended surfactant solution. The adsorption of the individual surfactants and surfactant blends were systematically investigated on different mineral surfaces using varying brine solutions. FINDINGS: LB adsorption on calcite surfaces was found to be significantly increased when blended with IOS or L38 since it forms surfactant complexes that partition to the surface. However, the total adsorption of the LB-IOS-L38 solution on dolomite decreased from 3.09 mg/m2 to 1.97 mg/m2 when blended together compared to summing the adsorption values of individual surfactants, which highlights the importance of mixed surfactant association.