Time-Domain Nuclear Magnetic Resonance Determination of Wettability Alteration: Analysis for Low-Salinity Water

Wettability has been shown to influence oil recovery. This property has become central to low-salinity (LSW) and smart (SWF) water flooding recovery mechanisms research. The challenge lies in the fact that oil recovery results from the combined effects of solid-liquid and liquid-liquid interactions. This demands methods that allow an independent interpretation of wettability alteration contributions. The primary objective of this work is to assess changes in wettability through the application of Time-Domain Nuclear Magnetic Resonance (TD-NMR) T 2 distribution and diffusion coefficient, starting with a well-controlled porous system, that is, glass beads, and then a model rock (Berea), in the presence of one phase, either oil or brine exclusively. Subsequently, two-phase fluid saturation was tested. For the glass beads, dimethyldichlorosilane was used to induce a hydrophobic response, as confirmed by contact angle experiments on slides of the same material. Sodium sulfate was used for its known positive influence on oil recovery during LSW and SWF. In cases where alteration of surface properties was expected, a leftward shift of the average T 2 distribution curve modes, accompanied by a reduction on the diffusion coefficient during the aging process was observed. The results of this work confirm that fluid-solid interactions during LSW and SWF, namely a shift in wettability, take place after the injection of low-salinity water.


Introduction
Wettability is a rock-fluid interaction that has been shown to influence oil recovery. Craig [1] defined wettability as "the tendency of a fluid to spread on or adhere to a solid surface in the presence of other immiscible fluids." This tendency results from the interaction between fluid molecules and solid surfaces. These interactions originate from van der Waals, structural or electrostatic forces [2][3][4]. Wettability can be altered by the interaction between the surface and the fluids occupying the pore space [5,6]. Alteration from oil-wet to slightly water-wet has been purported as the reason leading to an increase in the recovery factor (RF) [7][8][9][10]. Consistent characterization of rock wettability is valuable in enhanced oil recovery studies as this property has been found to correlate with the success of a waterflood [11], including those involving low-salinity water injection [6,12]. Requirements to adjust reservoir conditions vary depending on lithology. Wettability of sandstone reservoirs can be altered by oil [13,14], but in a waterflood, the salinity of the injection brine can be much lower than the connate one. Carbonate reservoirs are more influenced by which ions are present. Among them, Mg + , Ca 2+ , and SO 4 2− are considered potential determining ions [15,16]. However, sulfate has also shown a positive impact on several low-salinity waterflooding investigations for both sandstone and those focused on the fluids [17][18][19][20]. Wettability can affect multiphase flow properties such as relative permeability, capillary pressure, and residual oil saturation, and consequently its determination is critical [21][22][23][24]. Several methods, including the US Bureau Mines (USBM), imbibition and contact angle, have been the most commonly used ones. Amott cell imbibition and USBM experiments provide data about primary recovery of oil, but can also provide insights into the wettability conditions of the cores in use [25,26]. On the other hand, contact angle relies on the interaction forces among all fluid phases and solid surfaces of the system of interest (Equation (1)), and measurements help to determine the different conditions required to induce a change in the overall wettability of the reservoir [13,14]. In this case, the intrinsic or Young's contact angle corresponds to: where σ bs is the stress between the brine and the surface; σ os is the stress between the oil and the surface and σ bo is the interfacial tension between the fluids. A surface with an oil contact angle 0-75 • is considered oil-wet, for an angle 75-105 • the surface is considered intermediate wet, and finally an angle 105-180 • categorizes the surface as water-wet [27]. While the contact angle measurement has made its mark in petroleum research efforts, a technique with less experimental variability and an ability for in situ, non-destructive evaluation would enable this field to grow. Time-domain Nuclear Magnetic Resonance (TD-NMR) has been applied as an advanced technique to study solid-fluid interactions and could serve as a viable alternative [28][29][30][31][32][33][34].
Proton nuclear magnetic resonance ( 1 H-NMR) has been explored as a convenient technique to examine rock wettability for many years. The key to interpret the NMR response resides in the sensitivity of relaxation rates, particularly those associated to T 2 , which is sensitive to the presence of surfaces [14,28,[34][35][36]. The value of T 2 varies according to the various relaxation mechanisms, including those associated with the surface-to-volume ratio and consequently the T 2 distribution depends on pore size. To obtain these distributions, a multi-exponential decay function is fitted to the overall relaxation response of the spin-echo train. Similarly, the diffusion coefficient can also be derived by analogous analysis. Given a proper relaxation model, this technique can serve as a powerful tool in the determination of the pore-size distribution and wettability tendency of the surface of interest [28,29,[34][35][36][37]. Wang et al. [34] showed a relationship between the pore-size distribution and the surface relaxivity. In addition, the use of NMR analysis can allow us to analyze the influence of liquid films, previously addressed as a possible weakness of the traditional and standardized methods of wettability determination, as shown when complexation models are developed [12,35,38].
NMR enables inferences about pore size from relaxation rates by assuming that mobile fluids reside in large pores, bound fluids in small pores, and adhering fluids on grain surfaces. The intensity of the signal relates to the number of hydrogen nuclei found in each sample, while the rate of decay relates to the mobility of the molecules. Surface relaxation of the wetting phase will vary depending on the wettability environment within the pore space. A non-wetting fluid gives up a restricted pore space response, while the wetting phase will relax according to its degree of adhesion with the surface. The surface area (S) to pore volume (V) ratio is coupled to the surface relaxivity (ρ 2 ); for fluids confined in pores, T 2 is equal to or shorter than that of the bulk fluid [39]. In the fast diffusion limit, the measured T 2 values are given by Equation (2) [40]. As such, when surface relaxivity increases for given porous media, a decrease in T 2 response is expected, resulting in shifts to shorter T 2 values in correspondence with an increase in wettability to the fluid in use [29]. For example, if an oil-wet system is exposed to water, a shift to the left will occur, namely to shorter T 2 values, as the rock becomes more water-wet over time. The same response is expected for a water-wet system becoming more oil-wet when exposed to oil. In mixed-wet systems such as the case of a mix of oil-wet and water-wet glass bead, a linear relationship is observed between the relaxation response and wettability [36,37]. Both Brown et al. and Hsu et al. saw the same 4:1 ratio for surface relaxivity of a water-wet surface to an oil-wet surface. Plotting the NMR response by either the Amott or USBM techniques also provides a linear relationship, described as a wettability index, which improves the use of the technique and further validates it [29,30,32,33]. 1 The diffusion coefficient of fluid molecules restricted by pore surfaces depends not only on pore geometry [41], but also on pore surface characteristics [42]. In restricted geometries, the interaction with surfaces gives rise to a competition between the surface-liquid and liquid-liquid interactions. This competition determines the water spatial distribution, including water molecules on the bound water layer near the pore surface [42][43][44]. The pulsed field gradient NMR (PFG-NMR) can be used to analyze the diffusion coefficient of hydrogen-containing molecules restricted inside porous spaces, which is sensitive to surface-liquid interactions. T 2 trends similarly to T 1 [35], is less affected by bulk relaxation time [30] and turns out to be a relatively fast measurement. In principle, this relaxation should be a good indicator of wettability alteration, but a few challenges arise from changes in ionic profile, as the case of low-salinity or smart water injection. To elucidate this, we compare T 2 trends to the expected wettability alteration from contact angle measurements.

Materials
The oil used in this work had the following characteristics at 25 • C: density of 0.9092 g/cm 3 , a viscosity of 105 mPa·s, pentane asphaltene content of 10% wt and a refractive index of 1.527. The brines were prepared using analytical grade Na 2 SO 4 (Fisher Scientific, Pittsburgh, PA) salt diluted in de-ionized water, at concentrations of 1% and 100%, where 100% correspond to an ionic strength of 0.6724 M. The surfaces for this study corresponds to acid-washed glass beads (Sigma-Aldrich, St. Louis, MO) with a diameter between 425 and 600 µm and Berea Sandstone cores with a diameter of 0.8 cm, length of 1 cm, average porosity of 22% and average permeability of 720 mD were used to analyze wettability through NMR measurements. In addition, a 5% dimethyldichlorosilane solution in toluene was used to alter the glass beads surfaces.

Sample Preparation
To alter glass beads wettability, they were placed in a 5% dimethyldichlorosilane solution in toluene, then rinsed three times with methanol and dried in an oven at 60 • C. The change in surface properties was confirmed through contact angle measurements collected on glass slides (Fischer Scientific) that underwent the same process to induce hydrophobicity.
All glass beads were dried out at 80 • C to significantly reduce moisture, prior to measurements. Glass beads were added to an NMR tube to an approximate height of 1 cm, which amounted to approximately to a a mass of 0.6 g. Next, 300 µL of fluid was added to completely cover the glass beads, but avoiding excess free water, as the presence of bulk water would vastly overwhelm the measurement and invalidate the results. The tubes were sealed with a cotton applicator wrapped with Teflon tape to avoid excessive fluid evaporation during the aging time. Rock samples were vacuum-saturated for 6 h with the fluid of interest, then pressurized to 100 psi in the saturation cell. In the case of water, the aging process or equilibration time was carried out at 25 • C. We do not anticipate significant temperature dependence of wettability with only water. However, in the case of oil, the process is known to be sensitive to temperature, with the alteration being more significant at high temperature. Consequently, for the case of oil, samples were aged at 80 • C for three weeks, which is a standard wettability restoration procedure for rock, before exposure to any brine.

NMR T 2 Relaxation and Diffusion Coefficient Measurements
All NMR experiments were conducted on a Bruker mq20 miniSpec spectrometer with a frequency of 19.65 MHz. A chiller (F12-MA, JULABO USA Inc.) with a temperature accuracy better than ±0.02 • C was used to keep the sample at 25 • C by circulating constant-temperature liquid coolant. NMR T 2 relaxation time measurement was performed using a standard Carr-Purcell-Meiboom-Grill (CPMG) pulse sequence [45,46] with 200 echoes and an echo time of 14.21 µs. Figure 1 depicts the pulse sequence. This pulse sequence offers the advantage of reducing the effects of imperfections in the flip angle of the two-pulse Hahn sequence [47][48][49]. The effect is also seen in the broadening of the Fourier transform of the decay of the echo signal, which is narrower in the CPMG pulse sequence, considered only homogeneously broadened. The values of echo spacing, τ, the number of points and the gain were determined for each sample in order to get the cleanest signal possible. A total of 32 scans for each repeat and 5 separate measurements for each condition were conducted. Figure 2 shows the averaged T 2 of the dominant peak location over 5 repeats as a function of the number of scans. As shown, after 16 scans, the average location of the main T 2 converges to a narrow position, showing a standard deviation of 0.1 ms. The diffusion coefficient was determined by the pulsed field gradient spin echo (PFGSE) method [50] setting the following values for the needed parameter: gradient pulse with 0.5 ms, gradient pulse separation 7.51444 ms, 90 • -180 • pulse separation 7.51085 ms and 90 • first gradient pulse of 1 ms. T 2 measurements were executed for all the samples at days 1, 7, 14 and 21, and conversion of relaxation signal into T 2 distribution was processed using the CONTIN software provided by Bruker. The diffusion coefficient was collected concurrently on each measurement day.

Contact Angle and Interfacial Tension
The interfacial tension (IFT) and contact angle were measured with a Future Digital Scientific Corporation (FDSC)/Dataphysics instrument (OCA 25) pendant drop system. All experiments were conducted assuming symmetry with respect to the central vertical axis, ensuring also that the oil drop is elastic and the only forces shaping the drop under the influence of gravity and IFT. A J-shaped needle was used and submerged in a cuvette full of brine to form an oil-drop of 37 µL. The IFT was measured for 10 h until the droplet reached equilibrium. Regarding contact angle, a cuvette was filled with water until the glass surface at the top was completely in contact, a 4 µL oil drop was created with a J-shape needle, aged for 1/2 h and then measured. The camera contrast and the lens magnification were calibrated in order to obtain the best image quality and the most precise drop shape for both cases. Both measurements used the same software, which relies on the Young-Laplace equation and allows for direct measurement of the contact angle.

Results
Wettability can be altered at low temperature through different mechanisms and system composition [13,17,51,52]. In this study, regardless of aging temperature, TD-NMR data were collected at 25 • C in the spectrometer.

Interfacial Tension and Contact Angle
IFT was collected between oil and either of the two brines. Figure 3 shows that for 100% Na 2 SO 4 IFT value starts around 30 mN/m and plateaus at 15.8 mN/m. While for the 1% Na 2 SO 4 system, IFT starts at 28.5 mN/m and plateaus at 12.3 mN/m. This is shown to corroborate that IFT is not expected to play such a significant role in recovery differences in this case.
Each oil drop created in the pendant-drop system was given a stabilization time of 30 min, and the contact angle was tracked during this period, without finding any significant changes in the measured value. Dimethyldichlorosilane produces a hydrophobic behavior on treated surfaces, due to the absorption process on the solid surface [53,54]. Contact angle provided information to verify the surface changes to a hydrophobic response, confirming an oleophilic behavior and therefore the wettability alteration within the three components system.   Table 1 confirm that untreated glass slides are water-wet, while treated ones turned out oil-wet. The contact angle behavior reflects the effect of salinity, which suggests that when surfaces are water-wet, the increment in salinity goes along with an increment in the angle (122.6 ± 2.8 • for 1% Na 2 SO 4 and 139.6 ± 2.5 • for 100% Na 2 SO 4 ). On the other hand, for treated glass beads, low-salinity brine leads to an increase in the contact angle (33.1 ± 2.1 • for 1% Na 2 SO 4 and 21.4 ± 7.3 • for 100% Na 2 SO 4 ). This shows the direct effect on the surface, as the liquid-liquid interaction are the same for both treated and untreated slides.

Oil-Surface Interaction
The oil T 2 distribution for the bulk oil is shown in Figure 4, as a reference for results in the glass beads pack. The crude oil contains several solubility classes. The most polar molecules are usually associated with the asphaltene fraction, which is defined as being insoluble in low-molecular weight alkanes, for example, pentane, but soluble in toluene. Resins are closely related to the asphaltene fraction and can exhibit polar behavior as well. The resins can impact the solubility or stability of the asphaltene class in solution. The asphaltene fraction tends to contain high-molecular weight species, though not always, but they contribute significantly to wettability alteration, due to their affinity for surfaces. Also, short relaxation time is anticipated for this solubility class. Crude oil is complex in molecular composition, so the relaxation structure of the bulk is expected to be more complex than water. As can be seen from the figure, three distinct peaks were found in the T 2 spectrum. To estimate the error in the dominant peak location, the experiment was run five times and the standard deviation of the peak location was used as error metric. However, only one run is shown in the figures.  The equivalent oil T 2 distributions contained in untreated and treated glass beads packs saturated with oil and aged for 21 days at 25 • C are shown in Figure 5. The untreated glass beads (GB) give rise to longer relaxation times (dominant mode location: 40 ± 0.1 ms @ Day 1) corresponding to a water-wet (hydrophilic) surface. Over time, the dominant mode (highest peak) of the T 2 distribution shifted to shorter values (dominant mode location: 35.7 ± 0.1 ms @ Day 21), indicative of the alteration suffered by the media and corresponding to a more oil-wet surface. Figure 5b shows the T 2 distribution behavior in the treated glass beads (TGB) with oil during the period of 21 days. The slight decrease of the dominant mode position from 35.1 ± 0.2 ms to 34.7 ± 0.1 ms indicates that the pore surfaces are mildly altered by the interaction with oil, suggesting that after the chemical alteration, aging will still produce further oleophilicity. However, these changes are not significant in comparison to the effect observed for untreated glass beads. It is worth noticing that the other two modes/peaks associated to oil are also affected in similar fashion, namely that the oil residing in the untreated bead pack has all three modes shifted similarly. However, when comparing with the bulk oil response, the relative shifts and broadening of the modes seem to indicate compositional-dependent effects.  Figure 6 shows the diffusion coefficient as a function of aging time. This value in TD-NMR is an average over all molecular groups, so compositional-dependencies are hard to pinpoint. Despite this average behavior, for untreated glass beads, the oil diffusion coefficient decreases significantly during the 21 days from 0.070 ± 0.001 to 0.057 ± 0.002 × 10 −9 m 2 /s, in concert with the shifts in T 2 . By the same token, the result for treated glass beads shows a much smaller drop in the diffusion coefficient (0.056 ± 0.002 × 10 −9 m 2 /s), which is also consistent with the small shifts in the T 2 distribution.
The analysis of the two aforementioned parameters, that is, T 2 and D, indicates the acquisition of an oleophilic behavior of the surfaces within the porous space, being these chemically treated or aged in oil. In both cases, the alteration is the result of the adsorption process or ionic interaction between oil and the ions present on the surface [12,38]. An oleophilic surface will create a stronger attraction to oil molecules creating a layer effect, affecting the general dynamics of hydrogen bonds, producing a molecular reorganization and therefore a decrease on diffusion coefficient, but that also affects the T 2 distribution [42][43][44].

Brine-Surface Interaction
The effects of salinity during secondary recovery or smart-water flooding relate to the interaction between the phases present in the reservoir. Water chemistry plays a critical role in the recovery improvement, which has been linked to wettability alteration [6,12,55]. Sulfate is present as one of the most common ions in natural reservoir waters, and has shown the ability to induce changes in rocks of different lithologies, including sandstone [15,17,20].
The T 2 distribution in the case of water is simpler than that of crude oil, typically exhibiting one dominant mode or peak. The appearance of other modes at lower T 2 values is likely associated to water molecules sitting in smaller pores or pore-throats, corners, cranies and other small surface features. Figure 7 shows the water T 2 distribution inside both glass beads packs saturated with 1% Na 2 SO 4 . Untreated glass beads are considered naturally water-wet (hydrophilic). For this sample the dominant mode initially located at 1664.5 ± 0.1 ms on Day 1 decreases to 1620.1 ± 0.1 ms after 21 days of aging, suggesting that the surface is still being altered, though the surfaces only experiences a minute alteration. On the other hand, treated glass beads show a more dramatic shift in the dominant mode location from 2233.1 ± 0.1 ms (Day 1) down to 1620.5 ± 0.1 ms (Day 21). The inset plots correspond to the non-normalized spectra. For a single fluid occupying the pore space, we do not anticipate a significant decrease in intensity. The two obvious peaks in the two figures are affected similarly in qualitative sense. However, molecules sitting in a smaller spaces are exposed to larger surface-to-volume ratios, which implies more relative surface effects.    (Figure 8b) glass beads with 100% Na 2 SO 4 . For the untreated glass beads, the dominant mode remains at the same location, around 1602.5 ± 0.1 ms. Treated glass beads exhibit a decrease from 2120.1 ± 0.1 ms to 1756.1 ± 0.1 ms. In contrast to the effects of low-salinity, the degree of alteration is milder. However, a similar effect is observable on Day 1, confirming the equilibrium process of the sample. However, the dominant mode for low salinity on Day 7 locates at 1808.4 ± 0.1 ms, while for high salinity it is located at 1803.1 ± 0.1 ms, suggesting that after the equilibrium process both samples approach similar dominant mode locations. The smaller mode is also affected in a qualitative fashion, similarly to the largest mode.   Figure 9 shows the diffusion coefficient response as a function of time at two salinities for both treated and untreated glass beads packs. Analogous to previously observed results, when the beads are hydrophilic, the diffusion coefficient decreases (in concert with the T 2 distribution peak value). These results show that regardless of the media, 1% Na 2 SO 4 brine initially leads to a higher diffusion coefficient (2.272 ± 0.002 × 10 −9 m 2 /s for TGB and 2.234 ± 0.002 × 10 −9 m 2 /s for GB, respectively) in comparison to that of 100% Na 2 SO 4 (2.124 ± 0.002 × 10 −9 m 2 /s for TGB and 2.084 ± 0.003 × 10 −9 m 2 /s for GB, respectively). Dependency on ionic strength has been observed before [41]. On the other hand, the observable differences between TGB and GB are due to the surfaces characteristics being either hydrophobic or hydrophilic, respectively. Electrostatic forces play an important role in this alteration, since they can affect the mobility the hydrogen-containing molecules through a surface-bound layer [2,[42][43][44].  From the results for systems containing either oil or brines, a connection between TD-NMR results and those of contact angle experiments can be established. The untreated glass beads T 2 distribution is altered by the oil over time (Figure 5a). However, brine induces minimal differences in T 2 regardless of salinity (Figures 7 and 8a). Results shown in Table 1 suggest that high-and low-salinity brines give rise to different wettability conditions, which is not reflected in the two-component T 2 results. This implies the importance of fluid-mediated influence on wettability. Concerning the treated-glass beads, the oil does not alter an already oil-wet surface significantly. However, the low-salinity brine has a faster and greater effect on the pore surface over time (Figure 7b), while the high salinity brine elicits wettability changes to a lesser degree (Figure 8b). T 2 , along with those of contact angle results, show the improvement offered by low-salinity water on wettability alteration.

Rock-Fluid Interaction
Sandstone, alongside other common reservoir rocks, has been a dominant research focus through the evolution of TD-NMR wettability applications [13,14,28,35,37]. The wettability evolution of Berea Sandstone (initially water-wet) was examined here. Figure 10 shows a decrease in the T 2 dominant mode from 31.8 ± 0.1 ms on Day 7 to 25.5 ± 0.1 ms at day 21. The difference in the main mode/peak location for both day 7 and 21 with respect to untreated glass beds, correlates with the differences in pore sizes and mineral composition. Therefore, the shift in the dominant mode of the T 2 is indicative of the surface becoming increasingly oil wet, as expected for these conditions. Work on silica surfaces shows that polar fractions of the crude oil, for example, polar asphaltenes, adsorb on the surface [56]. This leads to wettability alteration and consequently to changes in relaxation mechanisms in fluid-rock interactions processes. Compared to the glass beads (Figures 7a and 8a), T 2 distributions of the brine-saturated rock in Figure 11 show a different behavior. Both untreated glass beads and sandstone are considered initially water-wet. However, brine still induces a decrease in the dominant mode location of the T 2 distribution. The increase of the ion exchange between rock surface and the brine, produces this effect. Nevertheless, the decrease is more appreciable for low salinity (from 152.6 ± 0.1 ms @ day 7 to 143.1 ± 0.1 ms @ day 21) than for high salinity (from 143.1 ± 0.1 ms @ day 7 to 139.1 ± 0.1 ms @ day 21). This reflects the different effect on wettablity alteration offered by low-salinity brine compared to high-salinity. The lower ionic strength of the low-salinity brine allows for expansion of the electric double layer on the surface and consequently the modification of the surface zeta potential. This change in the charge distribution is in part responsible for the change in the interaction and therefore of the surface-induced relaxation.   Figure 12a depicts the significant decreases in the diffusion coefficient, from 1.654 ± 0.002 × 10 −9 m 2 /s to 1.460 ± 0.002 × 10 −9 m 2 /s for low-salinity brine and from 1.347 ± 0.002 × 10 −9 m 2 /s to 1.279 ± 0.002 × 10 −9 m 2 /s for high-salinity brine is a verification of the ion exchange on the rock surface. However, this effect has been shown to be beneficial for wettability [6,55]. On the other hand, Figure 12b shows the diffusion coefficient for the oil-rock interaction, clearly indicates the adsorption process occurring as it decreases over time, assuming that no major chemical reaction occurs between the rock and oil molecules. The T 2 and diffusion coefficient results of the rock-oil or rock-brine systems analyzed together indicate wettability alteration of the surface. However, a new complication arises due to the heterogeneity offered by Berea Sandstone, steering the analysis to focus on trends instead of the specific T 2 peak location and diffusion coefficient values. Effects such as hydration layers around ions can affect the diffusion coefficient dependence on salinity, but could not explain the decrease over time observed in these results.

Partial Saturation
To analyze the effect of salinity on oil-wet samples, rocks initially aged in oil at 80 • C were later placed in the 1% Na 2 SO 4 and 100% Na 2 SO 4 brines. The T 2 distribution collected after the oil aging process is considered the initial state of the rock sample (Figure 13a). The rock wettability alteration after the aging process into the two different salinitiy brines is shown in Figure 13b. Figure 13a corresponds to the wetting state after oil aging, showing the dominant mode location at 28.6 ± 0.1 ms, which is comparable to the results in Figure 10, where the dominant mode is located at 25.5 ± 0.1 ms. The T 2 response was expected, regardless of the observed difference that is a confirmation of the heterogeneity of the rock although both samples corresponds to Berea Sandstone. This differences could complicate the analysis of the partially saturated samples, since for common methodologies wettability alteration is directly related to saturation [11,28]. Figure 13b shows the T 2 distribution of the rock sample after being aged in the two different brines. The appearance of a new peak between 100 and 200 ms corresponds to the brine intrusion into the rock. Since the effect of the the brine is the primary analysis target, a comparison between the trend offered by both brines can be developed. The dominant mode for 1% Na 2 SO 4 partially-saturated rock is located at 103.5 ± 0.1 ms, while for the 100% Na 2 SO 4 partially-saturated rock the dominant mode is located at 137.5 ± 0.1 ms, indicating the favorable effect of low-salinity on the pore surface and therefore the wettability alteration. This results cannot be directly related with the results in Figure 11a,b corresponding to rock fully saturated in 1% Na 2 SO 4 and 100% Na 2 SO 4 respectively, due to sample-to-sample variability. However, all systems show similar beneficial trend for low-salinity. A potential detachment of oil layers (or components) from the rock surface in the presence of low-salinity brine can produce a redistribution of phases in partially saturated systems, as the minimum of energy would require re-accommodation of fluid interfaces. This would be driven by changes in capillary pressure. The latter in conjunction with changes in relaxation mechanisms upon more exposure of the surface to water can contribute to the shifts in T 2 distribution modes, as shown in this section. This might help to provide an interpretation of T 2 distributions in partially saturated porous systems.  Figure 14 shows more clearly the differences between the initial state of the rocks and after the aging process with the respective brine. These two plots indicate that for higher salinity (Figure 14a) there is not a significant change on the oil component T 2 dominant mode location, going from 28.9 ± 0.1 to 25.9 ± 0.1 ms. Meanwhile, the differences for the low salinity (Figure 14b) are more significant going from 28.8 ± 0.1 to 21.38 ± 0.1 ms. This indicates that low-salinity brine improves the oil production. However, the different position of the brine component dominant mode between both salinities, indicates that changes on wettability induced by the low salinity are more beneficial. The diffusion coefficient, shown in Figure 15, increases for both systems, from 0.053 ± 0.002 × 10 −9 m 2 /s to 0.407 ± 0.002 × 10 −9 m 2 /s for 100% Na 2 SO 4 and from 0.050 ± 0.002 × 10 −9 m 2 /s to 0.321 ± 0.002 × 10 −9 m 2 /s for 1% Na 2 SO 4 , respectively. The lower value shown for 1% Na 2 SO 4 supports the existence of a wettability alteration. However, due to heterogeneity of the sample the use of the variable could not provide enough information for a three-component system since the final result is an average calculated from all the hydrogen containing molecules. Aged in 100% Na 2 SO 4 Aged in 1% Na 2 SO 4 Figure 15. Diffusion coefficient behavior of Berea Sandstone before and after the aging process in 1% Na 2 SO 4 and 100% Na 2 SO 4 at 25 • C.
Partially saturated systems have been analyzed for a considerable length of time using and a number of proposed analysis alternatives [28,29,34,36]. However, the increased level of complexity in these systems could lead to longer analysis times or even inadequate conclusions.

Discussion
The main contributors to the alteration of the surface properties involved on the surface-liquid interaction include adsorption, ion exchange as well as electrostatic forces. This can be seen in the process of transforming a surface from hydrophilic to oleophilic, naturally or artificially induced [13,53]. Examples of these are the response of the surface when in contact with oil for long enough or when the surface are chemically treated. In both cases an hydrophobic behavior is promoted.
The use of TD-NMR on two-component systems provides useful information to understand these phenomena. The T 2 distribution, which directly relates with the porous media through the surface relaxivity term in Equation (2), reflects the effect on the pore surfaces. On the other hand, the diffusion coefficient provides complementary information on the effects of the fluid. These two variables can jointly aid to interpret wettability changes. More complex heterogeneity and variability can significantly impact data interpretation. However, similar trends can still be observed in the different data types.
Wettability effects in complex systems have been widely explored [28,30,35,36]. The use of simpler two-component systems helps to assess wettability changes to a specific factor. Example of this is the wettability alteration induced by low-salinty sulfate brine, analyzed in this study. Although this work does not attempt to resolve the dichotomy related to fluid-fluid and fluid-rock interactions as it pertains to low-salinity waterflooding, namely which one is more dominant, it provides a simpler method to analyze the wettability alteration and the effects that the different fluids can have on it.