Characterization of Asphaltenes in Solution and Inside the Pores of

Jun 12, 2017 - One of the major issues of these processes is the diffusion of asphaltenes to the active site through the porous network of the alumina...
0 downloads 16 Views 4MB Size
Article pubs.acs.org/EF

Characterization of Asphaltenes in Solution and Inside the Pores of Catalysts by 1H NMR Relaxometry D. Espinat,† F. Gaulier,† F. Norrant,‡ J. Barbier,*,† B. Guichard,† M. Rivallan,† and P. Levitz§ †

IFP Energies nouvelles, Rond-point de l’échangeur de Solaize, BP3, 69360 Solaize, France IFP Energies nouvelles, 1 et 4 avenue de Bois Préau, 92852 Rueil-Malmaison, France § Laboratoire PHENIX, Université Pierre et Marie Curie, Case Courrier 51−4 place Jussieu, 75252 Cedex 5 Paris, France ‡

ABSTRACT: The upgrading of heavy petroleum fractions needs the development of more and more efficient heterogeneous catalysts. One of the major issues of these processes is the diffusion of asphaltenes to the active site through the porous network of the alumina support. The catalytic efficiency is deeply impacted by the transport phenomena and the interfacial interactions. The aim of this work is to capture the extent to which low-field two-dimensional (2D) 1H NMR relaxation time correlations can contribute to a better understanding of the dynamics of asphaltene in solution and within the pores of catalyst supports. Twodimensional T1−T2 maps for asphaltenes in solution in toluene exhibit several T1−T2 contributions, varying with the asphaltene concentration and the size of the asphaltenic fractions obtained by ultrafiltration separation. According to the nanoaggregate structure proposed by the Yen−Mullins aggregation model of asphaltenes, it was possible to unravel the different asphaltenic proton relaxation behaviors. By the use of NMR relaxometry, we have confirmed the stronger interaction of water with alumina than the one of toluene. The presence of macropores in catalyst clearly boosts the toluene mobility through the porous network. Two-dimensional T1−T2 maps for asphaltenes inside the pores show various types of protons, all of them with a severe constrained dynamics. Asphaltene nanoaggregates and clusters can be seen as large entities jammed into the pores, slowly mobile and affecting the solvent (toluene) mobility. When macroporosity exists in the support, the asphaltene overcrowding is less sensitive, enabling a faster dynamics of asphaltenes and toluene.

1. INTRODUCTION The conversion of heavy petroleum fractions into more valuable products, as for instance gasoline or diesel fuel, still remains an important target of the petroleum industry.1 This stake involves the need for developing more efficient catalysts to perform hydrocracking (HCK), hydrodesulfurization (HDS), and/or hydrodemetallization (HDM). All those catalysts exist industrially in sulfide heterogeneous catalysts. They usually consist of a combination of an active phase, made by molybdenum and cobalt or nickel sulfide, deposited on an alumina support which is nanoporous (average size of mesopores 10 nm), exhibiting a high specific surface area, usually over 200 m2/g. It is well-known that the main industrial issues depend on the presence of asphaltenes. These large macromolecules consist of a complex mixture of molecules, having different sizes and chemical composition.2 They are defined as the insoluble fraction in alkane, as for instance n-heptane (AsC7). During the latest 30 years, a very large number of contributions enabled great strides to be made in the description of asphaltenes.3−5 From the chemical point of view, asphaltenes are mainly constituted by carbon and hydrogen but some heteroatoms are inserted in the structure, such as sulfur and nitrogen, and metals, nickel and vanadium. On the basis of the results provided by several analytical methods2−10 (fluorescence, infrared spectroscopies, NMR analysis, EXAFS/XANES, and GC-pyrolyis, etc.) some contributions try to propose an average chemical structure. Two main representations are classically proposed:11−15 the first one corresponds to the “continental” type, in which asphaltenes exhibit a large central aromatic © 2017 American Chemical Society

region surrounded with some alkyl chains with variable lengths, and the second one, named ”archipelago” type, where the size of the aromatic regions is much smaller, these aromatic nuclei being linked together by alkyl chains. The advances of mass spectrometry (petroleomics16,17), with the high mass resolving power and the various methods for desorption and ionization,18−21 confirm an average asphaltene molecular weight close to 750 Da for the smallest entities. The main property of these various molecules is related to the selfassembling of small entities giving rise to aggregates characterized by a large polydispersity. The aggregation process of asphaltenes was peculiarly investigated by scattering techniques (X-ray (SAXS), neutron (SANS), or light).22−28 One of the most recent contributions used complementary neutron and X-ray small angle scattering data in order to give a prominent description of the asphaltene nanoaggregate.27 It appears as a disk-like particle of total radius 3.2 nm with 30% polydispersity and a thickness close to 0.67 nm. Several nanoaggregates can further aggregate themselves, forming larger clusters exhibiting extended polydispersity. An average mass fractal cluster can be described by the aggregation of roughly 12 nanoaggregates with a radius of gyration close to 7 nm. Results provided by a large pool of techniques confirm this macrostructure description nicely summarized by the Yen− Mullins model.29 We know that this macrostructure is strongly affected by an increase of the temperature.30 The average Received: January 13, 2017 Revised: May 29, 2017 Published: June 12, 2017 7382

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Virtually, we can distinguish in the literature two types of approaches, using high- or low-field NMR experiments. 1H and 13 C high-field NMR spectroscopies for asphaltene solutions are routinely used for the determination of their structural motifs.10,42 Solid-state 13C NMR and 2D HSQC-NMR spectroscopies were applied for a comprehensive characterization of various asphaltenes.43,44 Pulsed field gradient NMR has been extensively used for diffusion coefficient determination of various systems, and especially asphaltenes.45−47 The diffusion-ordered two-dimensional (2D) NMR spectroscopy (DOSY),48,49 which can be considered as a high-resolution version of the PFG sequence, was applied for asphaltene solutions investigation of different origins and at various concentrations.50,51 PFG-NMR spectroscopy was extensively used to probe liquid molecules trapped into nanoporous material. Many contributions dealt with the measurement of self-diffusion of molecules in zeolites.52 For the improvement of liquid-phase heterogeneous catalysts, the accessibility of the reactants to the active surface site must be as fast as possible. Consequently, many investigations using PFG-NMR technique were launched in order to better characterize this accessibility, peculiarly, by the diffusion coefficients of the molecules in the bulk pore space and on the surface of the catalysts.53−58 Very different catalytic supports have been investigated (γ-Al2O3, and TiO2, SiO2), displaying various surface area, chemical surface properties, and metallic active phase (Au, Pd, and Ag). PFGNMR gives access to the tortuosity of the support, based on the measurement of the effective self-diffusivity of the molecules inside the pores (Deff) and the free bulk liquid self-diffusivity (D0):

cluster size decreases with temperature, and if the number of nanoaggregates tends to decrease, it was demonstrated that an aggregation of small molecules by π−π interaction, which means stacking of the aromatic domains, is still present even at very high temperature (300 °C).28 The asphaltene transformation, during the refining process, still remains not completely understood. We previously insisted on the sizes of the aforementioned nanoaggregates and clusters which are, for the first ones, a little bit smaller than the pore sizes of the catalyst and, for the second ones, of the same range of size and even larger. The heavy fractions upgrading being processed at high temperature (370−400 °C) is in favor of the deaggregation behavior.1 Nevertheless, we have no precise information on the aggregation state of the asphaltene at this temperature. These simple observations suggest that the catalytic process will be decidedly affected by the molecular hindrance diffusion, through the catalytic pore network. We know, from industrial experience, that the contact time (liquid hourly space velocity, LHSV), can reach several hours1 in order to achieve a good conversion level. Likewise, it was clearly established that catalysts characterized by a bimodal or multimodal distribution of pores are more active than those having a monomodal porosity.31 All these observations have paved the way for an investigation of the asphaltene diffusion and adsorption through the catalytic pore network. From the experimental point of view, two main targets were identified: (1) The first one consists of the following, of the process diffusion, putting in contact various asphaltene solutions with the catalysts.32−36 The aim of these experiments is to follow the asphaltene penetration into the catalytic support present as beads (1−2 mm of diameter) or model extrudates (cylinders, 3 mm diameter and between 6 and 10 mm length). The first experimental approach for succeeding is straightforward: it consists of measuring the evolution of the asphaltene solution concentration versus time and to observe the penetration front in the support by optical microscopy (the catalytic support is a white alumina, and the asphaltene solution is black). From these experiments, diffusion models can be applied in order to extract key parameters governing the diffusion and adsorption process,32−35 especially the asphaltene diffusion coefficient in the porosity, the thermodynamic Langmuir coefficients of asphaltene and solvent (Langmuir law for adsorption isotherm), and the saturation concentration and the asphaltene average molecular weight. The fit of such a model to experimental data is made difficult because of the complexity of the system in particular due to the size polydispersity of asphaltenes and the multimodal distribution of pores. (2) The second one consists of the characterization of the diffusion of the various entities (asphaltenes and solvent) in the porosity. By this way, we aim at measuring some parameters previously mentioned (diffusion coefficient, tortuosity, and so on) and included in the models and get information on the dynamics of the molecules. Several contributions in the literature dealt with the membrane diffusion measurements.37−41 Various additional information can also be deduced from these experiments: (i) the diffusion coefficient when narrow size range fractions prepared by size-exclusion chromatography are concerned and (ii) the modification of the initial asphaltene size polydispersity and the chemical composition versus time. In support of these features, NMR technique can clearly provide information on the mobility of molecules in the pores.

τ=

D0 Deff

Nevertheless, Mantle et al.58 have shown that the tortuosity deduced from the previous equation is meaningful when the solvent used has no chemical interaction with the porous network. Virtually, this is the case when alkanes are considered. D′Agostino et al.57 have observed a singular behavior of different solvents within the porous medium, following the ratio D0/Deff compared to the one of a series of alkanes. Molecules, containing carbonyl functions, show a reduced self-diffusivity compared to alkanes. On the contrary, polyols show an enhanced self-diffusivity in the pores. The authors suggest the disappearance of the intermolecular hydrogen bonding network for polyols, giving rise to an enhanced translation and tumbling motion, in comparison with compounds for which the porous medium does not affect the hydrogen bonding network. Weber et al.54 have determined the bulk pore and pore surface diffusion coefficients of 1-octene diffusing in nanoporous 1 wt % Pd/θ-Al2O3 catalyst. On the base of a two site exchange model, the respective values obtained were the following: 1.3 × 10−9 m2·s−1 for bulk and 1.7 × 10−11 m2·s−1 for the surface and the mean residence time of the molecule on the pore surface of 150 ms. NMR relaxation times measurement is another way to observe dynamics of molecules in colloidal systems, especially petroleum systems,59−63 or in porous media.43,55−57,64−72 The longitudinal T1 and transverse T2 relaxation times are strongly dependent on the molecular dynamics and the ratio surface-tovolume when the liquid is trapped into the porous system. When the mobility of the spin of the molecule is reduced, both relaxation times tend to decrease. That is the situation when a 7383

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Figure 1. Evolution of T1 and T2 (relaxation times) as a function of the characteristic time (τc, s), according to expressions 2 and 3 (ν = 25 and 300 MHz) [wherein, for example, 1.0E+12 represents 1.0 × 1012].

molecules behave in the fast motion regime.43 Thus, the ratio T1/T2 or the two-dimensional T1−T2 relaxation correlation experiments can provide an insight into the mobility of the molecule and the strength of the interactions with the surface.43,71 T1 and T2 dependence of the Larmor frequency (NMR fast field cycling (FFC)) has been nicely used for a better understanding of (i) the influence of the pore structure on various types of molecules present during some catalytic transformations70 and (ii) the structure of crude oils with different asphaltene concentrations in relation with the viscosity.59,60 The final aim of this work is to capture the extent to which low-field 2D 1H NMR relaxation time correlations can contribute to a better understanding of the dynamics of asphaltene solution within the pores of catalyst supports. At first, to give a straightforward answer to this question, we have targeted to analyze asphaltenes in solution in toluene for different concentrations. Two different asphaltene fractions, obtained by nanofiltration, are studied, in order to address the asphaltene size polydispersity and to discuss the nonfractionated asphaltenes data. Finally, NMR relaxation times are measured for asphaltene solutions within the porous network.

molecule has an interaction with the surface. An analytical model of the relaxometry T1 and T2 was proposed, based on a biphasic fast exchange between the bulk and the surface:64−66 1 1 λS 1 = + T1,2 TB ,1,2 V TS ,1,2

(1)

where T1,2, TB,1,2, and TS,1,2 are respectively the longitudinal (1) and transverse (2) relaxation times measured, in the bulk solvent (B) and at the surface (S). For complex systems that contain many different species present in various environments, the relaxation behavior cannot be described by one relaxation time. We usually consider a multiexponential behavior, and we can have access to the distribution function of relaxation times by inverse-Laplace transformation.62,68 The relaxation times depend on the Larmor frequency (ν = ω/2π) and a characteristic or correlation time (τc). According to the physical model proposed by Bloembergen, Purcell, and Pound (BPP model73), the following equations in the situation of homonuclear dipolar coupling were obtained:74,75 ⎡ ⎤ 2τc 8τc 1 ⎥ = 2K ⎢ + 2 2 2 2 T1 1 + 4ω τc ⎦ ⎣ 1 + ω τc

⎡ ⎤ 10τc 4τc 1 ⎥ = K ⎢6τc + + 2 2 2 2 T2 1 + ω τc 1 + 4ω τc ⎦ ⎣

K=

⎛ μ0 ⎞2 4⎛ h ⎞2 I(I + 1) ⎜ ⎟ γ ⎜ ⎟ ⎝ 2π ⎠ ⎝ 2π ⎠ 40r 6

(2)

2. EXPERIMENTAL SECTION 2.1. Samples. 2.1.1. Asphaltenes. Asphaltenes have been extracted by n-heptane precipitation from a Safaniya vacuum residue according to the NF T60−115 method. Asphaltenes characteristics are listed in Table 1. n-Heptane was purchased from Aldrich with 99% grade. In order to reduce the asphaltene size polydispersity and be able to study some selected asphaltenes according to their size, an ultrafiltration separation was performed. The experimental setup used for this fractionation consists of a small dead-end filtration cell housing flat membrane sheets. A poly(ether sulfone) porous membrane with a 20 kDa molecular weight cutoff was used. A 2 wt % asphaltene solution in toluene (Aldrich, 99.5% grade) was submitted to the membrane filtration at 25 °C and for a 10 bar feed pressure. Retained and filtrated asphaltene fractions, respectively collected in the feed compartment of the cell and at the end of the filtration experiment, were recovered, and the toluene was further eliminated (see ref 75 for more information on this experimental development).

(3)

(4)

where μ0 is the vacuum permeability, I is the spin number (1/2 for protons), γ is the gyromagnetic ratio, h is Planck’s constant, and r is the internuclear distance. τc can be defined as the time taken by one molecule to rotate of 1 rad. Figure 1 introduces the evolution of T1 and T2 as a function of τc for two frequencies (25 and 300 MHz). For very low values of τc, T1 = T2, and for higher values, T1 goes up instead of T2 plummeting. A frequency variation (25 or 300 MHz) does not produce any variation of T2 relaxation time. If the frequency specific of the motion (1/τc) is larger than the Larmor frequency (ν), 7384

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels Table 1. Elemental Composition of (i) Asphaltenes Extracted from a Safaniya Vacuum Residue and (ii) Retained (Big Asphaltenes) and Filtrated Asphaltene (Small Asphaltenes) Fractions Obtained by Ultrafiltration Process76

ratio H/C carbon (wt %) hydrogen (wt %) sulfur (wt %) vanadium (ppm) nickel (ppm)

retained asphaltenes

filtrated asphaltenes

asphaltenes

big asphaltenes

small asphaltenes

1.10 ± 0.05 82.3 ± 0.4 7.5 ± 0.1 8±2 550 ± 50 190 ± 30

1.10 ± 0.03 81.9 ± 0.4 7.4 ±0.1 8±2 555 ± 30 185 ± 7

1.14 ± 0.03 81.0 ± 0.4 7.3 ± 0.1 8±2 271 ± 15 68 ± 2

Table 2. Porous Characteristic of Alumina Supports M14 nm and B14 nm M14 nm SBET (m2/g)a parameter C from the BET VBJH (mL/g)b mesoporous macroporous dmean (nm)c structural density porosity

165 ± 10 98 ± 5

0.80 ± 0.04 / 14.5 ± 0.7 3.22 ± 0.32 71 ± 4

0.60 ± 0.04 0.35 ± 0.02 14.3 ± 0.7 (mesopores) 3.05 ± 0.32 75 ± 4

a

SBET = specific surface area. bVBJH = mesoporous/macroporous volumes. cdmean = mean pore diameter. level). Thus, the magnetization signal depends on two times, τ1 and nτ, which can be later plotted as a T1−T2 2D map computed with a homemade 2D inverse Laplace transform software.77 For the different experiments, the resolution for the relaxation time determination is close to 50 μs. We used one NMR probe of diameter 10 mm and 10 mm height. The detection of “solid-like protons” needs very short probe dead time, especially when analyzing alumina supports, with very short relaxation times. For the 10 mm probe, the dead time is down to 5 μs. This allows a short 90° pulse duration to reach a larger frequency bandwidth when detecting solid protons. The filter dead time can also be reduced significantly using a 1 MHz bandwidth, at the expense of the signal-to-noise ratio.68 The main characteristics of the NMR Instrument are summarized in Table 3.

All these asphaltene samples were characterized from the elemental composition and size-polydispersity points of view. Table 1 gives the elemental composition. Size-exclusion chromatography (SEC) was used for a size-polydispersity investigation76 (Figure 2). The molecular weight is given in equivalent polystyrene molecular mass. In practice, fractions are enriched in small (further named Small Asphaltenes) or large entities (further named Big Asphaltenes) depending on the one considered. 2.1.2. Porous Media. Two different alumina (γ-Al2O3) supports were investigated. Each one can be representative of the classical supports used for hydrotreatments. The first one (M14 nm) exhibits a meso-single-mode pore size distribution close to 14 nm. The second one (B14 nm) is characterized by a bimodal porous distribution, the mesoporosity being centered to 14 nm and the macroporosity reaching pore diameters near 1000 nm. The characteristics of these two alumina supports are summarized in Table 2. Both supports are cylindrical extrudates with 3 mm diameter and an average 10 mm length. 2.2. NMR Methods. The NMR experiments have been performed on a Maran Ultra proton spectrometer from Oxford Instrument having a proton Larmor frequency of 23.7 MHz. Free induction decay (FID) and transverse magnetization CPMG decay curves were measured. Two-dimensional T1−T2 maps were obtained by an inversion recovery sequence followed by a CPMG sequence:

P180 − τ1 − P90 − (τ − P180 − echo)n

B14 nm

194 ± 10 113 ± 6

Table 3. Main Characteristics of the NMR Instrument (P90 is the 90° Pulse Duration)

(5)

probe diam (mm)

probe dead time (μs)

P90 (μs)

10

5

2

filter dead time at 100 kHz/(1 MHz) (μs)

τ spacing (μs)

first FID point (μs)

15 5

25 15

21 11

For the solution experiments, asphaltenes were dissolved in perdeuterated toluene (99.55% D), supplied by Eurisotop, in order to remove most of the protons coming from the solvent (toluene).

where τ1 is respectively equal to 40 values spaced logarithmically and the number of echos (n) is chosen in order to reach zero signal (noise

Figure 2. SEC (size-exclusion chromatography) size distributions (Mw, molecular weight) for nonfractionated asphaltenes (−) and for asphaltene fractions obtained by ultrafiltration (retained fraction − big asphaltenes (--); filtrated asphaltenes − small asphaltenes (···)). 7385

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Figure 3. T1 − T2 maps for asphaltene solutions in deuterated toluene, for different concentrations (2, 3, and 5 wt %, nonfractionated asphaltenes) and for small and big asphaltenes obtained by ultrafiltration. water. Then, they were soaked in water, toluene, or asphaltene solutions for a given concentration. The time in contact with the solvent (water or toluene) was fixed to 24 h, in order to allow the solvent to diffuse through the entire porous network. For the asphaltene solutions, the presence of large asphaltenic entities, impose a much longer contact time, fixed to 6 weeks in order to reach equilibrium of diffusion and adsorption within the pores.36 After these

The solutions were prepared 2 weeks before any use. The asphaltene solutions were introduced into small flasks of 6 mm diameter, filled up, and thoroughly sealed in order to avoid any evaporation of the solvent. Each flask was further put inside the NMR probe for experimental purpose. For the alumina support experiments, the extrudates were previously dried at 100 °C during 6 h in order to remove adsorbed 7386

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels immersion times, the impregnated supports were collected, wiped on filter paper, and introduced into the NMR probe. NMR experiments were performed at a controlled temperature at 30 °C.

toluene. We have checked that the relaxation times (T1 and T2) for pure toluene are in the same range, equal to 2.66 s, close to 2.0−2.2 s. We can observe a clear tendency to a decrease of both solvent relaxation times as the asphaltene concentration increases, the ratio T1/T2 being quite constant (1) for 2 and 3 wt % asphaltene concentration, but slightly increasing for a 5 wt % concentration (2.2) (see Table 3). This observation fits with previous results of the literature,78,59−61 dealing with either asphaltene solutions or pure crude oils. Considering the biphasic fast exchange model (see eq 1), asphaltene nanoaggregates and clusters can be seen as a smooth surface where some toluene molecules can be in tight interaction, corresponding to a shorter surface relaxation time (TS,1,2). Consequently, the measured relaxation times (T1,2) will be shorter as asphaltene concentration increases.78 A more concentrated asphaltene solution will induce higher viscosity and a general slowing of the motions of the toluene entities. The spot no. 1 (Figure 3) spans a rather wide range of relaxation times, for T1 (200−300 ms) and T2 (80−160 ms). The ratio T1/T2 is roughly constant as a function of the concentration, but shows a slight increase for big asphaltenes. Peaks nos. 2−5 exhibit rather similar T1 relaxation times in the range 45−160 ms, but very different T2 (0.1−30 ms) relaxation times. This behavior is witness to molecular (or spins) motions with characteristic inverse correlation times similar to or slower than the Larmor frequency.61 Moreover, the ratio T1/T2 rises sharply for peaks 3−5, proving a huge decrease of the molecular mobility. It is important to notice that spots 3, 4, and 5 are not seen for small asphaltenes solution. We have plotted (Figure 4) the evolution of the T2 relaxation time distributions for asphaltene solutions at three different concentrations (2, 3, and 5 wt %). The distributions were rescaled considering the same mass of asphaltene solution analyzed by NMR spectroscopy. Moreover, the distribution of the relaxation times of the empty NMR cell was measured and used for the correction of the data carried out for the asphaltene solutions. We can see (Figure 4) that the rise of the asphaltene concentration is followed by a shift of the distribution toward much smaller T2 values. This observation props up the previous discussion, focusing on an overall severe reduction of the molecular mobility when the asphaltene concentration goes up. In low-field 1H NMR experiments, we can observe all the protons of the sample. For the asphaltene solutions, these protons come from asphaltenes but also from the nondeuterated fraction of the toluene used as solvent. This is one limitation of the technique which can prevent an accurate discussion of the 2D T1−T2 maps. In that way, the good resolution of high-field 1H NMR analysis is clearly more efficient allowing the separation of the different types of proton and their relaxation time acquisition.43 Nonetheless, we have tried to estimate the number of protons coming respectively from nondeuterated toluene and asphaltenes at rather low concentration (some few weight percent). It is a nontrivial issue, because we do not know exactly what the different hydrogenated impurities of the deuterated toluene are. We have made a simple hypothesis, among which the hydrogen atoms are located on the methyl group of the partly hydrogenated toluene (-CH3). According to the elemental analysis, we know the hydrogen asphaltene concentration (see Table 1). Thus, a simple calculation provides an assessment of the ratio H from asphaltenes to H from the toluene (Table 5). We can conclude

3. RESULTS AND DISCUSSION 3.1. Asphaltene Solutions. Three nonfractionated asphaltene solutions with concentrations of 2, 3, and 5 wt % and the two fractions of asphaltenes (Small and Big) for a 2 wt % concentration were analyzed by low-field NMR spectroscopy. Two-dimensional T1−T2 maps are presented in Figure 3. For all these maps, we can clearly distinguish six peaks which are unambiguously detected, but more or less intense. The values of relaxation times, measured on top of the peak, are summarized in Table 4. The ratio T1/T2 was also calculated. Table 4. Relaxation Times and Ratio T1/T2 Measured from 2D T1−T2 Maps of Asphaltene Solutions in Deuterated Toluene (See Figure 3) signal

T1 (ms)

T2 (ms)

2 wt % Asphaltene Solution 2.2 × 103 2.2 × 103

solvent toluene no. no. no. no. no.

1 2 3 4 5

solvent toluene no. no. no. no. no.

1 2 3 4 5

solvent toluene no. no. no. no. no.

1 2 3 4 5

solvent toluene no. 1 no. 2 solvent toluene no. 1 nos. 2 and 3 no. 4 no. 5

3.0 × 102 1.5 × 8.0 × 101 3.0 × 6.0 × 10 6.0 6.0 × 10 1.0 3.2 0 × 10 3.0 × 3 wt % Asphaltene Solution 2.0 × 103 2.0 ×

102 101

2.0 × 102 9.0 × 9.0 × 10 2.0 × 1.2 × 102 3.0 7.0 × 10 8.0 × 8.0 × 10 3.0 × 5 wt % Asphaltene Solution 1.6 × 103 7.0 × 2.8 × 102 1.6 7.4 × 10 2.0 7.0 × 10 5.0 8.0 × 10 1.0 4.5 × 10 1.5 2 wt % Small Asphaltenes 2.2 × 103 2.2

× × × ×

102 10 10 10

1

10−1

2 2.7 10 60 107

103

1

101 101 10−1 10−1

2.2 5 40 88 267

102

2.2

× 102 × 10

× 10−1

2 4 14 80 300

× 103

1

3.0 × 102 1.5 × 102 4.5 × 10 3.0 × 10 2 wt % Big Asphaltenes 2.2 × 103 2.2 × 103 3.0 6.0 7.5 7.5

T1/T2

8.0 × 10 9.5 1.5 1.2 × 10−1

2 1.5 1 4 6 50 625

The first one, on the top right, corresponds to T1 and T2 values close to 2.0−2.2 s. This signal comes from the small part of nondeuterated toluene molecules. Even if the range of T2 values is wider than the range of T1, we can clearly attribute this peak to free toluene; this means that solvent molecules are present in an environment similar to the one found for pure 7387

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Figure 4. T2 relaxation time distributions for various asphaltene solutions (2−5 wt %) in deuterated toluene. Intensities correspond to the same mass of asphaltene solution observed by NMR spectroscopy (continous line for 2 wt% ; -- for 3 wt% and ... for 5 wt%).

the asphaltene protons contribute mainly to peak intensities 1− 5, but we cannot rule out some contribution of the hydrogen of the toluene. As discussed previously, some solvent molecules in interaction with nanoaggregates or asphaltenic clusters may show a reduced mobility or like trapped in a sort of network made by the asphaltene species, especially for a 5 wt % concentration. This behavior was discussed by Korb et al.57,58 characterizing the dynamics of small molecules of crude oils near asphaltene aggregates. Toluene molecules, such as light hydrocarbons of crudes, follow a pathway marked, first, by translation diffusion in a more or less continuous phase (solvent) surrounded by the asphaltenic network and, second, by a more or less long adsorption or trapping within nanoaggregates or clusters. To assert this description, a simple

Table 5. Ratio Number of Protons from Asphaltenes to Number of Protons from Deuterated Toluene for Different Asphaltene Concentrations asphaltene concentration (wt %)

NH,A/NH,T

2 3 5

4.2 6.3 10.7

that the number of protons, brought by the solvent, remains relatively small in comparison with the asphaltene hydrogen atoms. If we look at the solvent contribution on the 2D T1−T2 maps, we can highlight a certain decrease of its intensity as the asphaltene concentration rises, due to fewer protons from toluene and more from asphaltene molecules. Consequently,

Figure 5. Internal structure of the asphaltenic nanoaggregate proposed by Eyssautier et al.27,28 Black lines represent the aliphatic chains of elementary asphaltene molecules stacked into the nanoaggregate (an illustration gives examples of this type of elementary molecule). White lines sketch in the central aromatic nucleus as can be illustrated on the elementary molecule. 7388

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

peak (no. 5), characterized by a very low T2 (1.2 × 10−1 ms) and then a very high T1/T2 ratio (625). This last contribution has been previously detected (see in particular the 5 wt %. asphaltene concentration). It suggests the occurrence of a large number of protons with a significant reduction of the mobility. Several works in the literature22−30,76 have shown that the asphaltene aggregation state is the main characteristic of this big asphaltene solution (see Figure 2). We can forecast that a large population of asphaltene or toluene protons will get into a tangle and their motions will be strongly reduced. All these observations are in good agreement with the work of Dutta Majumbar et al.43 (see Figure 4 in ref 43). High-field NMR experiments (300 MHz) provide of course high resolution for the assignment of the different chemical groups of asphaltenes (CH, CH2, and CH3; see Figure 3 in ref 43), which is a significant advantage in low-field NMR spectroscopy. The values measured for relaxation time T2 are in the range 24.4−86.9 ms, in reasonably good agreement with our data. For peaks 1 and 2, T2 varies between 20 and 160 ms (Figure 3 and Table 3). The authors clearly distinguish aliphatic CH3 or CH2 groups that are located on the inside of the cluster and show a restricted motion. Conversely, CH3 or CH2 groups on the periphery are more mobile and consequently have shorter T2’s. From the T1/T2 ratio and according to the homonuclear dipolar coupling model (see formulas 2 and 3 and Figure 1), it is possible to assess the characteristic time, τc. Values of the T1/ T2 ratio in the range of 2−300 correspond to values of τc between 5 and 10 ns. These values are significantly higher than the typical values for small molecules nearby a few picoseconds. An assessment of τc for various sizes of molecules can be obtained using the Einstein−Stokes equation considering molecules as spheres of diameter d moving in a fluid of viscosity η:

calculation allows assessing the average distance between nanoaggregates. Considering the average molecular weight (16000 Da) given by Eyssautier et al.,27 we have found distances between 8 and 11 nm regarding the experimental concentrations (2−5 wt %). Having in mind that the asphaltenic nanoaggregate described by Eyssautier et al.,27 is a disk-like particle with a diameter of 3.2 nm and a thickness close to 0.7 nm, we see that the previous calculated internanoaggregate distance is in the same range of size of this aggregate. In fact, the formation of polydispersed clusters according to the Yen−Mullins model,29 made by the aggregation of several nanoaggregates,27,79 conveys the existence of regions enriched in asphaltenes in the solution and consequently surrounding regions where the solvent is in the majority. This concentration heterogeneity of the solution was already observed by small angle x-ray scattering (SAXS) and by cryo-scanning electron microscopy.23 We can suggest that the ratio T1/T2 for H-toluene molecules trapped into the asphaltene dense regions will be high. The different asphaltenic proton relaxation behaviors observed on the 2D T1−T2 maps can be unraveled considering the nanoaggregate structure proposed by Eyssautier et al.27,28 A sketch of this nanostructure is presented on Figure 5. As aforementioned, this nanoaggregate appears as a disk of radius close to 3.2 nm, a thickness of near 0.7 nm, and an average molecular weight of 16 kDa as described in Figure 5. The main driving force for the building of this nanostructure is π−π interaction, of around two or three aromatic sheets. The whole nanoaggregate contains between 15 and 20 elementary molecules. Coupling SAXS and SANS techniques,27, it was clearly figured out that no solvent can diffuse into the nanoaggregate core, contrary to what can occur inside the shell of the aggregate where the solvent can be present surrounding aliphatic or alicyclic chains. In this harsh environment, we can easily suggest that the protons mobility will be more or less affected. The spin of the proton of jammed molecules at the center of the nanoaggregate will relax very differently from the spin of protons attached to aliphatic chains (-CH2 or -CH3) outside the aggregate. We have proposed an illustration for the assigning of the five peaks seen on 2D T1−T2 maps in Figure 5. The observation of the 2D T1-T2 maps for small and big asphaltenes (Figure 3) remarkably props up the previous attribution of the different spots of the maps for pure asphaltenes. Small asphaltenes’ two-dimensional representation only exhibits two main contributions, the first one depending on the solvent and the second one due to protons of asphaltenes and toluene, in a smaller proportion. In that way, two peaks (nos. 1 and 2, the most intense) can be distinguished with a rather small T1/T2 ratio between 1.5 and 2 (see Table 4). These small entities are mainly represented by elementary asphaltene molecules, as illustrated by the modified Yen− Mullins model,29 and the contribution of nanoaggregates in this fraction is clearly very weak. Subsequently, we can suggest that the mobility of the asphaltene protons is much less affected than that in nonfractionated asphaltenes where the aggregation state is predominant. This interpretation is asserted by the observation of the big asphaltenes 2D T1−T2 map. In addition to the free toluene contribution, we find four peaks (nos. 1−4, Figure 3) again similar to those previously observed for a 5 wt % pure asphaltene concentration. The intensities of these peaks are much lower that the intensities measured for 5 wt % asphaltene concentration, but we can perceive a fifth intense

τc =

πd3η 6kT

where T is the temperature and k the Boltzmann constant. In a very poor approximation, if we consider the disk-like nanoaggregate as a sphere of diameter 6.4 nm, moving in toluene (η = 0.59 × 10−3 Pa·s), τc is found close to 20 ns in reasonably good agreement with values deduced from the NMR experiments. 3.2. Asphaltenes within the Porous Network. 3.2.1. Two-Dimensional T1−T2 Map for the Dry Alumina Support. The two-dimensional relaxation map for the dry alumina support M14 nm (see Table 3) is presented in Figure 6. An intense signal is observed corresponding to T1 and T2 respectively equal to 4.0 and 0.3 ms. The dispersion of the values of relaxation time T2 is narrower than the one of T1. A very weak and diffused signal for T2 close to 40 μs, at the limit of detection, is due to the presence of hydroxyls on the support. It has been already observed for kaolinite and smectite,68 but for T2 value close to 60 μs. The T1/T2 ratio of the main contribution is close to 13, a good indication of the restricted mobility of the protons. Consequently, it can be straightforwardly attributed to some remaining water molecules not completely eliminated during the drying. The NMR measurement of the magnetization decay, the extrapolation to time zero and the comparison with a fixed mass of water in similar experimental conditions, has provided an estimation of the water content of the support equal to 2 wt %. This latest estimation is rather coherent taking into account that no-dried 7389

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

The determination of the porosity was possible thanks to the magnetization measurement at initial time proportional to the spin number and consequently the mass or volume of toluene inside the pores. Knowing the mass of the dried support and the structural density (see Table 2) we can have access to the NMR porosity equal to 69% in good agreement with the porosity given by mercury porosimetry. 3.2.3. Two-Dimensional T1−T2 Map of Asphaltenes Solution in Toluene Trapped into the Pores. For these experiments, two supports M14 nm and B14 nm were put in contact, at room temperature, with 2 wt % asphaltene solutions in deuterated toluene during 5 or 6 weeks, in fact a rather long contact time but strictly necessary for a satisfactory penetration of asphaltene entities inside the porosity.36 After 5 weeks of waiting time, low-field NMR experiments were performed and the 2D T1−T2 maps are shown in panels c (M14 nm) and d (B14 nm) of Figure 8. After 6 weeks contact time, some impregnated extrudates of the M14 nm support were recovered and dried at 100 °C during 6 h in order to evaporate the major part of the deuterated toluene present in the pores. Twodimensional T1−T2 results are presented in Figure 8b. Just for clarity of the discussion, the 2D T1−T2 map for the nonimpregnated and dried support M14 nm is recalled in Figure 8a. In order to pick up a good signal-to-noise ratio, the acquisition time of each 2D map was long (12 h). Relaxation times measured from the 2D maps are listed in Table 7. A maximum of four peaks are detected from the T1−T2 maps. For the M14 nm support impregnated with asphaltenes and dried, three main signals appear, characterized by low T1 (20 ms) and T2 (0.3 ms) and high T1/T2 ratio (6 to 200). For the supports, still wetted by toluene, we found more or less the same characteristics of the dried impregnated support, except for a signal 4, more intense and spread out at large range of T2, for the multimodal catalytic support (B14 nm). Looking at the map of the impregnated and dried support (M14 nm, Figure 8b), the major contribution (peak 3) fits with a T2 close to 40 μs, in the limit of the T2 instrumental resolution, lower than the peak detected for nonimpregnated and dried support (Figure 8a). We have previously attributed a diffuse and weak signal (T2 = 40 μs) on M14 nm support to hydroxyls (see section 3.2.1). With the presence of asphaltenes on the support, the intense peak 3 (Figure 8b) cannot be only attributed to hydroxyls. In order to address this point, we can appraise the ratio of the number of protons of asphaltenes to the number of protons belonging to hydroxylic groups (NH,a/ NH,OH). The calculation is based on the classical value for alumina80,81 of 7 OH/nm2. We have measured the amount of asphaltene on the support after 6 weeks in contact with the solution, equal to 92 mg/g (of support) or 8.4 wt %. Consequently, the ratio NH,a/NH,OH is close to 4.4 and the contribution of asphaltene protons to the intensity of peak 3 is high, clearly showing a very restricted mobility of asphaltene entities. On the other hand, two additional peaks (nos. 1 and 2; Figure 8b) have relaxation times not very far from the contribution of the clean support M14 nm (Figure 8a). We have previously attributed this last contribution to water molecules adsorbed at the pore surface. Even if the support has been carefully dried before putting it into contact with the asphaltene solution, we can guess that a small quantity of water could still be trapped on the support, especially during the transfer between the oven and the solution. We have figured out the ratio of the number of protons from asphaltene to the number

Figure 6. T1−T2 map for the dried alumina support M14 nm (see Table 2).

alumina support, such as the one we have been using, contains water between 5 and 10 wt % at room temperature. 3.2.2. Two-Dimensional T1−T2 Map of the Solvent Inside the Porous Network. Figure 7 shows the 2D T1−T2 maps for the two alumina supports M14 nm and B14 nm. The supports were soaked into two different solvents, water and toluene. T1 and T2 values are listed in Table 6. Both values are lower than the relaxation times of free solvents, respectively equal to 2.66 s for toluene and 2.7 s for water, conveying, as expected, the confinement of the liquids into the pores. Two different behaviors of toluene and water inside the porosity are clearly distinguished. For water, T2 is much smaller than T1, when, for the toluene system, T1 and T2 are in the same order of magnitude. This confirms a classical observation of the literature for protic and aprotic liquids in contact with porous material.64,65,70 Unlike toluene, water can easily exchange protons with the surface of the support. T1 and T2 relaxation times of protons of the toluene significantly increase when the liquid wets the B14 nm support. This is due to the macroporosity of the support, boosting the molecular mobility through the porous matrix. This observation is also in agreement with the general trend of the increase of the relaxation times with the pore diameter64,65 observed by Korb et al. on calibrated granular packing, The relaxation times of water confined into mesoporous (M14 nm) or macroporous (B14 nm) supports are similar, with a ratio T1/T2 respectively equal to 12 and 7 in good agreement with the literature.69,70 This latest observation suggests that relaxation times are mainly governed by the water molecules interactions with the solid, especially the proton exchange as previously mentioned, and not by the porous network properties, as for instance the pore size distribution. In a similar sense concerning porosity characterization, D′Agostino et al.55 have shown that the true tortuosity of the porous network can be truly determined with small probe molecules which do not exhibit any chemical functions being liable to interact with the surface (alkanes). Moreover, according to D′Agostino et al.,71 the ratio T1/T2 is related to the surface adsorption energy. 7390

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Figure 7. Two-dimensional T1−T2 maps for two solvents, water and toluene, inside the pore network of the two alumina supports M14 nm and B14 nm (see Table 2): (a and b) water in respectively M14 nm and B14 nm; (c and d) toluene in respectively M14 nm and B14 nm.

known for their amphiphilic properties.82 Thus, an adsorption competition between water and asphaltene entities may occur at the alumina surface. Consequently, a contribution of asphaltenic protons to peaks 1 and 2 cannot be totally excluded, corresponding to protons having certain mobility. Two-dimensional maps of supports soaked with toluene and asphaltenes (Figure 8c,d) look like the map observed in Figure 8b. For the M14 nm support, peak 3 confirms the presence of protons with restricted motions probably from adsorbed asphaltenes and residual water or hydrogenated toluene molecules in strong interaction with the surface or the asphaltene deposit. Spots 1 and 4, with a rather low ratio T1/ T2 (Table 7), are probably linked to relatively mobile toluene in the pores. The diffusion and adsorption process of asphaltene aggregates through the porosity of meso- or macroporous alumina supports is very slow.36,83,84 After soaking during 5

Table 6. Relaxation Times T1 and T2 (ms) for Water and Toluene Inside the Porosity of Two Alumina Supports M14 nm and B14 nm (See Table 2 and Figure 7) solvent

support

T1 (ms)

T2 (ms)

T1/T2

water

M14 nm B14 nm M14 nm B14 nm

7.2 6.1 1.7 5.8

× × × ×

6.0 8.4 1.1 × 102 4.8 × 102

12 7 1.6 1.2

toluene

10 10 102 102

of protons brought by adsorbed water. The water quantity adsorbed on M14 nm is close to 2 wt %, as estimated by NMR measurements (see section 3.2.1). With an asphaltene concentration of 8.4 wt %, a simple calculation shows that the number of protons from the asphaltene is twice higher than the protons coming from the adsorbed water molecules. Due to the complex chemical composition, asphaltenes molecules are 7391

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Figure 8. Two-dimensional T1−T2 maps for different configurations of alumina supports M14 nm and B14 nm: (a) dried alumina support M14 nm (see Figure 6); (b) support M14 nm first in contact (6 weeks) with a 2 wt % asphaltene solution in deuterated toluene and then dried at 100 °C (6 h); (c and d) respectively M14 nm and B14 nm in contact (5 weeks) with a 2 wt % asphaltene solution in deuterated toluene and not dried.

Table 7. Relaxation Times T1 and T2 (ms) and Ratio T1/T2 Deduced from 2D Maps (See Figure 8) signal

T1 (ms)

T2 (ms)

asphaltene density measured by mercury porosimetry (da = 1.22 ± 0.12), the volume occupied by asphaltene species represents more or less 10% of the total porous volume (M14 nm). The asphaltene deposit along the radius of the extrudate is not homogeneous, more concentrated on the edges of the cylindrical extrudate.83,84 The pore diameter distribution (M14 nm) during the asphaltene diffusion and adsorption process exhibits a loss of largest pores (15−20 nm). In summary, asphaltene nanoaggregates and clusters can be sketched as large entities jammed into the pores, with very slow mobility. This congestion of the pore texture also affects the solvent (toluene) mobility, strongly reduced in comparison with solvent alone inside the porous network. The B14 nm 2D T1−T2 map supports this behavior, and peak 4 with a wide T2 distribution and a quite low T1/T2 ratio is the sign of the macroporosity. In fact, we can propose that the asphaltene overcrowding in large pores will be less sensitive causing a faster dynamics of asphaltenes and toluene protons.

T1/T2

Support M14 nm Dried no. 1 4.0 3.0 × 10−1 12 Support M14 nm Impregnated with Asphaltenes (2 wt %) and Dried no. 1 1.7 3.0 × 10−1 6 no. 2 1.8 × 10 3.0 × 10−1 60 no. 3 2.0 × 10 4.0 × 10−2 200 Support M14 nm Impregnated with Asphaltenes (2 wt %) no. 1 2.0 3.0 × 10−1 7 no. 2 2.0 × 10 2.0 × 10−1 100 no. 3 2.0 × 10 5 × 10−2 400 no. 4 30−100 1−20 5−50 Support B14 nm Impregnated with Asphaltenes (2 wt %) no. 1 1.0 1.7 × 10−1 6 nos. 2 and 3 4.0 × 10 1.0 × 10−1 400 no. 4 20−100 0.3−13 6−60

4. CONCLUSION In this contribution, we have shown that H NMR relaxation times measurements can contribute to the investigation of the dynamics of molecules especially when they are diffusing through a porous network. Two-dimensional T1−T2 maps of

weeks, we are not very far from a stationary regime for which the concentration of the asphaltene solution surrounding the extrudates of alumina support does not vary, probably far from the thermodynamic equilibrium. The amount of asphaltene inside the porosity is close to 8.4 wt %. According to the dried 7392

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels

Alberta Oil Sand Bitumen and Heavy oil Asphaltenes. Energy Fuels 1992, 6, 83−96. (6) Waldo, G. S.; Mullins, O. C.; Penner-Hahn, J. E.; Cramer, S. P. Determination of the chemical environment of sulphur in petroleum asphaltenes by X-ray absorption spectroscopy. Fuel 1992, 71, 53. (7) Speight, J. G.; Pancirov, R. J. Structural types in petroleum asphaltenes as deduced from pyrolysis/gas chromatography/mass spectrometry. Liq. Fuels Technol. 1984, 2 (3), 287−305. (8) Mitra-Kirtley, S.; Mullins, O. C.; van Elp, J.; George, S. J.; Chen, J.; Cramer, S. P. Determination of the nitrogen chemical structures in petroleum asphaltenes using XANES spectroscopy. J. Am. Chem. Soc. 1993, 115, 252−258. (9) Goulon, J.; Retournard, A.; Friant, P.; Goulon-Ginet, C.; Berthe, C.; Muller, J. F.; Poncet, J. L.; Guilard, R.; Escalier, J. C.; Neff, B. Structural characterization by X-ray Absorption Spectroscopy (EXAFS/XANES) of the Vanadium Chemical Environment in Boscan asphaltenes. J. Chem. Soc., Dalton Trans. 1984, 1095−1103. (10) Bouquet, M.; Bailleul, A. Routine method for quantitative 13C n.m.r editing and providing structural patterns. Fuel 1986, 65, 1240− 1246. (11) Murgich, J. Molecular Simulation and the Aggregation of the Heavy fractions in Crude Oils. Mol. Simul. 2003, 29, 451−461. (12) Bunger, J. W., Li, N. C., Eds. Chemistry of asphaltenes; Advances in Chemistry Series, Vol. 195; American Chemical Society: Washington, DC, USA, 1981; DOI: 10.1021/ba-1981-0195. (13) Strausz, O. P.; Mojelsky, T. W.; Lown, E. M. The molecular structure of asphaltene: an unfolding story. Fuel 1992, 71, 1355−1363. (14) Pacheco-Sanchez, J. H.; Alvarez-Ramirez, F.; Martinez-Magadan, J. M. Morphology of Aggregated Asphaltene Structural Models. Energy Fuels 2004, 18, 1676−1686. (15) Acevedo, S.; Castro, A.; Negrin, J. G.; Fernandez, A.; Escobar, G.; Piscitelli, V.; Delolme, F.; Dessalces, G. Relations between Asphaltene Structures and Their Physical and Chemical Properties: The Rosary-type structure. Energy Fuels 2007, 21, 2165−2175. (16) Rodgers, R. P.; Marshall, A. G. Petroleomics: Advanced Characterization of Petroleum-Derived Materials by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR/MS). In Asphaltenes, Heavy Oils, and Petroleomics; Mullins, O. C., Sheu, E. Y., Hammami, A., Marshall, A. G., Eds.; Springer; New York, 2007; Chapter 3, p 63. (17) Hsu, C. S.; Hendrickson, C. L.; Rodgers, R. P.; McKenna, A. M.; Marshall, A. G. J. Mass Spectrom. 2011, 46, 337−343. (18) McKenna, A. M.; Marshall, A. G.; Rodgers, R. P. Heavy Petroleum Composition. 4. Asphaltene Compositional Space. Energy Fuels 2013, 27, 1257−1267. (19) Pomerantz, A. E.; Wu, Q.; Mullins, O. C.; Zare, R. N. LaserBased Mass Spectrometric Assessment of Asphaltene Molecular weight, Molecular Architecture, and Nano-aggregate Number. Energy Fuels 2015, 29, 2833−2842. (20) Sabbah, H.; Pomerantz, A. E.; Wagner, M.; Müllen, K.; Zare, R. N. Laser Desorption Single-Photon Ionization of Asphaltenes: Mass range, Compound Sensitivity, and Matrix Effects. Energy Fuels 2012, 26, 3521−3526. (21) Tanaka, R.; Sato, S.; Takanohashi, T.; Hunt, J. E.; Winans, R. E. Analysis of the Molecular Weight Distribution of Petroleum Asphaltenes Using Laser desorption Mass Spectrometry. Energy Fuels 2004, 18, 1405−1413. (22) Ravey, J. C.; Ducouret, G.; Espinat, D. Asphaltene macrostructure by small angle neutron scattering. Fuel 1988, 67, 1560−1567. (23) Espinat, D.; Rosenberg, E.; Scarsella, M.; Barré, L.; Fenistein, D.; Broseta, D. Colloidal structural evolution from stable to flocculated state of asphaltene solutions and heavy crudes. In Structures and Dynamics of Asphaltenes; Mullins, O. C., Sheu, E. Y., Eds.; Plenum Press: New York, 1998; Chapter V, p 145. (24) Liu, Y. C.; Sheu, E. Y.; Chen, S. H.; Storm, D. A. fractal structure of asphaltene in toluene. Fuel 1995, 74 (9), 1352−1356. (25) Fenistein, D.; Barré, L.; Broseta, D.; Espinat, D.; Livet, A.; Roux, J. N.; Scarsella, M. Viscosimetry and Neutron Scattering Study of

asphaltene solutions for various concentrations can be interpreted using the aggregation model proposed by Eyssautier et al.27 or the Yen−Mullins model.29 The asphaltene nanoaggregate exhibits several kinds of protons; those trapped inside the heart of the nanoaggregate have a clear reduction of their mobility when those located on the edges of this nanoaggregate are more mobile, probably belonging to aliphatic chains. As the asphaltene concentration increases, we have observed a decrease of the solvent molecules (toluene) dynamics. The investigation of asphaltenic fractions, of low and high molecular weights, made us assert the 2D T1−T2 map peak identification. We have discussed the presence of protons coming from the nondeuterated toluene which can influence the peak identification as the asphaltene concentration becomes very low. The correlation time, τc, was estimated from T1/T2 ratio and in good agreement with the average size of the asphaltene nanoaggregate. The next step will be to use high-field NMR sequences for relaxation time measurements in order to get access to a higher NMR resolution enabling the observation of the different types of protons present in the asphaltenic solution.43 We have emphasized on the BPP model used for T1 and T2 evolution regarding τc. We recall that this model is based on the hypothesis that the time correlation function has an exponential form.73,74 We can suspect that this model is not perfectly true when molecules and thus protons are confined into small pores. In that way, NMR relaxometry experiments varying temperature could be interesting. According to the solvent and asphaltenes behavior into the alumina supports, we first investigated the dynamics of the solvent. We have confirmed the stronger interaction of water with alumina than the one of toluene. The effect of the presence of macropores is clearly observed, boosting the toluene mobility through the porous network. Two-dimensional T1−T2 maps for asphaltenes inside the pores show various types of protons, all of them with severely constrained dynamics. Thus, asphaltenes nanoaggregates and clusters in the porosity are characterized by a very slow mobility, in complete agreement with the very long times of diffusion of these entities through the pore network of alumina beads or extrudates used as industrial catalysts.36,83,84



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

J. Barbier: 0000-0002-2101-0799 Notes

The authors declare no competing financial interest.



REFERENCES

(1) Catalysis by transition metal sulphidesFrom Molecular Theory to Industrial Application; Toulhoat, H., Raybaud, P., Eds.; Editions TECHNIP: Paris, 2013. (2) Merdrignac, I., Espinat, D., Hénaut, I., Argillier, J. F. Properties and Composition. In Heavy Crude OilsFrom Geology to Upgrading, An Overview; Huc, A. Y., Ed.; Editions TECHNIP: Paris, 2011; Chapter 4, p 34. (3) Structures and Dynamics of Asphaltenes; Mullins, O. C., Sheu, E. Y., Eds.; Plenum Press: New York, 1998. (4) Asphaltenes, Heavy Oils, and Petroleomics; Mullins, O. C., Sheu, E. Y., Hammami, A., Marshall, A. G., Eds.; Springer; New York, 2007. (5) Mojelsky, T. W.; Ignasiak, T. M.; Frakman, Z.; McIntyre, D. D.; Lown, E. M.; Montgomery, D. S.; Strausz, O. P. Structural Features of 7393

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels Asphaltene Aggregates in Mixed Toluene/Heptane Solvents. Langmuir 1998, 14, 1013−1020. (26) Fenistein, D.; Barré, L. Experimental measurement of the mass distribution of petroleum asphaltene aggregates using ultracentrifugation and small-angle X-ray scattering. Fuel 2001, 80, 283−287. (27) Eyssautier, J.; Levitz, P.; Espinat, D.; Jestin, J.; Gummel, J.; Grillo, I.; Barré, L. Insight into Asphaltene Nanoaggegate Structure Inferred by Small Angle Neutron and X-ray Scattering. J. Phys. Chem. B 2011, 115, 6827−6837. (28) Eyssautier, J.; Espinat, D.; Gummel, J.; Levitz, P.; Becerra, M.; Shaw, J.; Barré, L. Mesoscale Organization in a Physically Separated Vacuum Residue: Comparison to asphaltenes in a Simple Solvent. Energy Fuels 2012, 26, 2680−2687. (29) Mullins, O. C.; Sabbah, H.; Eyssautier, J.; Pomerantz, A. E.; Barré, L.; Andrews, A. B.; Ruiz-Morales, Y.; Mostowfi, F.; McFarlane, R.; Goual, L.; Lepkowicz, R.; Cooper, T.; Orbulescu, J.; Leblanc, R. M.; Edwards, J.; Zare, R. N. Advances in Asphaltene Science and the Yen-Mullins Model. Energy Fuels 2012, 26, 3986−4003. (30) Espinat, D.; Fenistein, D.; Barré, L.; Frot, D.; Briolant, Y. Effects of Temperature and Pressure on Asphaltenes Agglomeration in Toluene. A Light, X-ray and neutron Scattering Investigation. Energy Fuels 2004, 18, 1243−1249. (31) Plumail, J. C., Jacquin, J., Martino, G., Toulhoat, H. Effect of the Pore-size Distribution on the Activities of Alumina Supported Co-Mo Catalysts in the Hydrotreatment of Boscan Crude. Preprints, American Chemical Society, Division of Petroleum Chemistry; 185th ACS National Meeting, Mar. 20−25, 1983, Seattle, WA, USA; American Chemical Society: Washington, DC, USA, 1983; Vol. 28, pp 562−575. (32) Mieville, R. L.; Trauth, D. M.; Robinson, K. K. Asphaltene characterization and diffusion measurements. Prepr. - Am. Chem. Soc., Div. Pet. Chem. 1989, 34 (3), 635−643. (33) Yang, X.; Guin, J. A. Diffusion-controlled adsorptive uptake of coal and petroleum asphaltenes in a NiMo/Al2O3 hydrotreating catalyst. Chem. Eng. Commun. 1998, 166, 57−79. (34) Tayakout, M.; Ferreira, C.; Espinat, D.; Arribas Picon, S.; Sorbier, L.; Guillaume, D.; Guibard, I. Diffusion of asphaltene molecules through the pore structure of hydroconversion catalysts. Chem. Eng. Sci. 2010, 65, 1571−1583. (35) Marchal, C.; Abdessalem, E.; Tayakout-Fayolle, M.; Uzio, D. Asphaltene Diffusion and Adsorption in Modified NiMo Alumina Catalysts Followed by Ultraviolet (UV) spectroscopy. Energy Fuels 2010, 24, 4290−4300. (36) Gaulier, F.; Barbier, J.; Guichard, B.; Levitz, P.; Espinat, D. Asphaltenes Transport into Catalysts under Hydroprocessing Conditions. Energy Fuels 2015, 29, 6250−6258. (37) Baltus, R. E. Characterization of asphaltenes and heavy oils using hydrodynamic property measurements. In Structures and Dynamics of Asphaltenes; Mullins, O. C., Sheu, E. Y., Eds.; Plenum Press: New York, 1998; Chapter X, p 303. (38) Baltus, R. E.; Anderson, J. L. Hindered diffusion of asphaltenes through microporous membranes. Chem. Eng. Sci. 1983, 38 (12), 1959−1969. (39) Kyriacou, K. C.; Sivaramakrishna, V. V.; Baltus, R. E.; Rahimi, P. Measurement of diffusion coefficients of oil residual fractions using porous membranes. Fuel 1988, 67, 15−18. (40) Sane, R. C.; Tsotsis, T. T.; Webster, I. A.; Ravi-Kumar, V. S. Studies of asphaltenes diffusion and structure and their implications for resid upgrading. Chem. Eng. Sci. 1992, 47 (9−11), 2683−2688. (41) Dechaine, G. P.; Gray, M. R. Membrane Diffusion Measurements Do Not Detect Exchange between Asphaltene Aggregates and Solution Phase. Energy Fuels 2011, 25, 509−523. (42) Scotti, R.; Montanari, L. Molecular structure and intermolecular interaction of asphaltenes by FT-IR, NMR, EPR. In Structures and Dynamics of Asphaltenes; Mullins, O. C., Sheu, E. Y., Eds.; Plenum Press: New York, 1998; Chapter III, p 79. (43) Dutta Majumdar, R.; Gerken, M.; Mikula, R.; Hazendonk, P. Validation of the Yen-Mullins Model of Athabasca Oill-Sands Asphaltenes using Solution-State 1H NMR relaxation and 2D HSQC Spectroscopy. Energy Fuels 2013, 27, 6528−6537.

(44) Dutta Majumdar, R.; Bake, K. D.; Ratna, Y.; Pomerantz, A. E.; Mullins, O. C.; Gerken, M.; Hazendonk, P. Single Core PAHs in Petroleum- and Coal-Derived Asphaltenes: Size and Distribution from Solid-State NMR Spectroscopy and Optical Absorption Measurements. Energy Fuels 2016, 30, 6892−6906. (45) Ostlund, J.-A.; Andersson, S.-I.; Nyden, M. Study of asphaltenes by the use of pulsed- field gradient spin echo NMR. Fuel 2001, 80, 1529−1533. (46) Lisitza, N. V.; Freed, D. E.; Sen, P. N.; Song, Y. Q. Study of Asphaltene Nano-aggregation by Nuclear Magnetic Resonance (NMR). Energy Fuels 2009, 23, 1189−1193. (47) Norinaga, K.; Wargardalam, V. J.; Takasugi, S.; Iino, M.; Matsukawa, S. Measurement of Self-Diffusion Coefficient of Asphaltene in Pyridine by Pulsed Field Gradient Spin-Echo NMR. Energy Fuels 2001, 15, 1317−1318. (48) Morris, K. F.; Johnson, C. S. Diffusion-Ordered TwoDimensional Nuclear Magnetic Resonance Spectroscopy. J. Am. Chem. Soc. 1992, 114, 3139−3141. (49) Johnson, C. S., Jr. Diffusion ordered nuclear magnetic resonance spectroscopy: principles and applications. Prog. Nucl. Magn. Reson. Spectrosc. 1999, 34, 203−256. (50) Durand, E.; Clemancey, M.; Lancelin, J. M.; Verstraete, J.; Espinat, D.; Quoineaud, A. A. Aggregation States of Asphaltenes: Evidence of two Chemical Behaviors by 1H Diffusion-Ordered Spectroscopy Nuclear Magnetic Resonance. J. Phys. Chem. C 2009, 113, 16266−16276. (51) Durand, E.; Clemancey, M.; Lancelin, J. M.; Verstraete, J.; Espinat, D.; Quoineaud, A. A. Effect of Chemical Composition on Asphaltenes Aggregation. Energy Fuels 2010, 24, 1051−1062. (52) Kärger, J.; Ruthven, D. M. Diffusion in nanoporous materials: fundamental principles, insights and challenges. New J. Chem. 2016, 40, 4027−4048. (53) Weber, D.; Mantle, M. D.; Sederman, A. J.; Gladden, L. F. Surface diffusion in Catalysts Probed by APGSTE NMR. Diffus. Fundam. 2009, 10, 1−3. (54) Weber, D.; Sederman, A. J.; Mantle, M. D.; Mitchell, J.; Gladden, L. F. Surface diffusion in porous catalysts. Phys. Chem. Chem. Phys. 2010, 12, 2619−2624. (55) D’Agostino, C.; Mitchell, J.; Gladden, L. F.; Mantle, M. D. Hydrogen Bonding Network Disruption in Mesoporous Catalyst Supports Probed by PFG-NMR Diffusometry and NMR Relaxometry. J. Phys. Chem. C 2012, 116, 8975−8982. (56) D’Agostino, C.; Kotionova, T.; Mitchell, J.; Miedziak, P. J.; Knight, D. W.; Taylor, S. H.; Hutchings, G. J.; Gladden, L. F.; Mantle, M. D. Solvent Effect and Reactivity trend in the Aerobic Oxidation of 1,3-Propanediols over Gold Supported on Titania: NMR Diffusion and Relaxation Studies. Chem. - Eur. J. 2013, 19, 11725−11732. (57) D’Agostino, C.; Brett, G. L.; Miedziak, P. J.; Knight, D. W.; Hutchings, G. J.; Gladden, L. F.; Mantle, M. D. Understanding the Solvent Effect on the Catalytic Oxidation of 1,4-Butanediol in Methanol over Au/TiO2 catalyst: NMR Diffusion and Relaxation Studies. Chem. - Eur. J. 2012, 18, 14426−14433. (58) Mantle, M. D.; Enache, D. I.; Nowicka, E.; Davies, S. P.; Edwards, J. K.; D’Agostino, C.; Mascarenhas, D. P.; Durham, L.; Sankar, M.; Knight, D. W.; Gladden, L. F.; Taylor, S. H.; Hutchings, G. J. Pulsed-Field gradient NMR Spectroscopic Studies of Alcohols in supported Gold Catalysts. J. Phys. Chem. C 2011, 115, 1073−1079. (59) Korb, J. P.; Louis-Joseph, A.; Benamsili, L. Probing Structure and Dynamics of Bulk and Confined Crude Oils by Multiscale NMR Spectroscopy, Diffusometry and Relaxometry. J. Phys. Chem. B 2013, 117, 7002−7014. (60) Korb, J. P.; Vorapalawut, N.; Nicot, B.; Bryant, R. G. Relation and Correlation between NMR relaxation Times, Diffusion Coefficients, and Viscosity of heavy Crude Oils. J. Phys. Chem. C 2015, 119, 24439−24446. (61) Zielinski, L.; Saha, I.; Freed, D. E.; Hürlimann, M. D.; LIu, Y. Probing Asphaltene Aggregation in Native Crude Oils with Low-Field NMR. Langmuir 2010, 26 (7), 5014−5021. 7394

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395

Article

Energy & Fuels (62) Zielinski, L.; Hürlimann, M. D. Nuclear Magnetic Resonance Dispersion of Distributions as a Probe of Aggregation in Crude Oils. Energy Fuels 2011, 25, 5090−5099. (63) Ordikhani-Seyedlar, A.; Neudert, O.; Stapf, S.; Mattea, C.; Kausik, R.; Freed, D. E.; Song, Y.-Q.; Hürlimann, M. D. Evidence of Aromaticity-Specific Maltene NMR Relaxation Enhancement Promoted by Semi-immobilized Radicals. Energy Fuels 2016, 30, 3886− 3893. (64) Korb, J. P. Nuclear magnetic relaxation of liquids in porous media. New J. Phys. 2011, 13, 035016. (65) Korb, J. P.; Godefroy, S.; Fleury, M. Surface nuclear magnetic relaxation and dynamics of water and oil in granular packings and rocks. Magn. Reson. Imaging 2003, 21, 193−199. (66) Godefroy, S.; Korb, J.-P.; Fleury, M.; Bryant, R. G. Surface nuclear magnetic relaxation and dynamics of water and oil in macroporous media. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2001, 64, 021605. (67) Levitz, P.; Korb, J. P.; Bryant, R. G. Molecular diffusion in disordered interfacial media as probed by pulsed field gradients and nuclear magnetic relaxation dispersion. J. Chim. Phys. Phys.-Chim. Biol. 1999, 96, 1494−1505. (68) Fleury, M.; Kohler, E.; Norrant, F.; Gautier, S.; M’Hamdi, J.; Barré, L. Characterization and Quantification of Water in Smectites with Low-Field NMR. J. Phys. Chem. C 2013, 117, 4551−4560. (69) Weber, D.; Mitchell, J.; McGregor, J.; Gladden, L. F. Comparing Strengths of Surface Interactions for Reactants and Solvents in Porous Catalysts Using Two-Dimensional NMR Relaxation Correlations. J. Phys. Chem. C 2009, 113, 6610−6615. (70) Mitchell, J.; Broche, L. M.; Chandrasekera, T. C.; Lurie, D. J.; Gladden, L. F. Exploring Surface Interactions in Catalysts Using LowField Nuclear Magnetic Resonance. J. Phys. Chem. C 2013, 117, 17699−17706. (71) D’Agostino, C.; Mitchell, J.; Mantle, M. D.; Gladden, L. F. Interpretation of NMR Relaxation as Tool for Characterising the Adsorption Strength of Liquids inside Porous Materials. Chem. - Eur. J. 2014, 20, 13009−13015. (72) Ralphs, K.; D’Agostino, C.; Burch, R.; Chansai, S.; Gladden, L. F.; Hardacre, C.; James, S. L.; Mitchell, J.; Taylor, S. F. R. Assessing the surface modifications following the mechanochemical preparation of a Ag/Al2O3 selective catalytic reduction catalyst. Catal. Sci. Technol. 2014, 4, 531−539. (73) Bloembergen, N.; Purcell, E. M.; Pound, R. V. Phys. Rev. 1948, 73, 679. (74) Abragam, A. Principles of Nuclear Magnetism; Oxford University Press: New York, 1983. (75) Bakhmutov, V. I. Practical NMR Relaxation for Chemists; John Wiley & Sons, 2004; DOI: 10.1002/0470094486. (76) Marques, J.; Merdrignac, I.; Baudot, A.; Barré, L.; Guillaume, D.; Espinat, D.; Brunet, S. Asphaltenes Size Polydispersity Reduction by Nano- and Ultrafiltration Separation Methods − Comparison with the Flocculation Method. Oil Gas Sci. Technol. 2008, 63 (1), 139−149. (77) Fleury, M.; Soualem, J. J. J. Colloid Interface Sci. 2009, 336, 250− 259. (78) Jestin, J.; Barré, L. Application of NMR Solvent Relaxation and SAXS to Asphaltenes Solutions Characterization. J. Dispersion Sci. Technol. 2004, 25 (3), 341−347. (79) Eyssautier, J.; Frot, D.; Barré, L. Structure and Dynamic Properties of Colloidal Asphaltene Aggregates. Langmuir 2012, 28, 11997−12004. (80) Digne, M.; Sautet, P.; Raybaud, P.; Euzen, P.; Toulhoat, H. J. Catal. 2002, 211, 1−5. (81) Wischert, R.; Copéret, C.; Delbecq, F.; Sautet, P. Angew. Chem., Int. Ed. 2011, 50, 3202−3205. (82) Sjöblom, J.; Saether, Ø.; Migttun, Ø.; Ese, M.-H., Urdahl, O., Førdedal, H. Asphaltene and Resin Stabilized Crude Oil Emulsions. In Structure and Dynamics of Asphaltenes; Mullins, O. C., Sheu, E. Y., Eds.; Plenum Press: New York, 1998; Chapter XI, p 337.

(83) Gaulier, F. Etude de la diffusion des charges lourdes en conditions réelles dans les catalyseurs d’hydrotraitement. Thesis, Université Claude Bernard Lyon I, Villeurbanne, France, 2016. (84) Barbier, J.; Gaulier, F.; Guichard, B.; Levitz, P.; Espinat, D. Asphaltenes Transport into Hydroconversion Catalysts at High Temperature - The Role of Alumina Nanoporous Texture. Energy Fuels 2017, DOI: 10.1021/acs.energyfuels.7b00824.

7395

DOI: 10.1021/acs.energyfuels.7b00139 Energy Fuels 2017, 31, 7382−7395