and solubility parameters

Eng. 2015; 10: 1–14. (62) Mousavi, M.; Abdollahi, T.; Pahlavan, F.; Fini, E. H.The influence of asphaltene-resin molecular interactions on the collo...
4 downloads 14 Views 828KB Size
Subscriber access provided by UNIV OF DURHAM

Article

A successful approach to mimic the solvent power of maltenes based on SARA analysis, solvatochromic- and solubility parameters Luzia P Novaki, Raphael Lira, Michelle M. N. Kwon, Marcia Cristina Khalil de Oliveira, Francis A Meireles, Gaspar Gonzalez, and Omar A. El Seoud Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b04064 • Publication Date (Web): 11 Feb 2018 Downloaded from http://pubs.acs.org on February 13, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

Energy & Fuels

[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A successful approach to mimic the solvent power of maltenes based on SARA analysis, solvatochromic- and solubility parameters

Luzia P. Novaki,a Raphael Lira;a Michelle M. N. Kwon,a Marcia C. K. de Oliveira,b,* Francis A. Meireles,b Gaspar Gonzalez,c Omar A. El Seoud,a,* a- Institute of Chemistry, the University of São Paulo, 748 Prof. Lineu Prestes Av., 05508-000 São Paulo, SP; e-mail: [email protected] b- Flow assurance laboratory, Petrobras Research Center, (CENPES), Av. Horácio Macedo, 950, Cidade Universitária, 21941-915 Rio de Janeiro, RJ email: [email protected] c- Federal University of Rio de Janeiro , 21941-598 Rio de Janeiro, RJ

Abstract We extended the oil compatibility model to the dissolution of asphaltenes (Asps) in maltenes from ten crude oils (COs). As scales for the power of solvents of interest, vide infra, we used solvatochromic parameters, calculated from the Uv-Vis spectra of solvatochromic compounds (probes), Hildebrand/Hansen solubility parameters, and the colloidal instability index of COs. As the colors of maltenes or asphaltene-free crude oils CO(Asp-free) were too dark to permit recording the absorption spectra of the probes, we formulated models for these fractions, (MCO(Asp-free)). They were composed of low molar mass hydrocarbons, namely cis and trans decalines, isooctane, 1-methylnaphthalene and- as model for resins- benzothiazole/n-octyl-1-naphthoate. We based formulations of these MCO(Asp-free) on SARA analysis of the COs, and elemental analysis of the corresponding resins. We validated MCO(Asp-free) as models for the corresponding CO(Aspfree)

by showing that the correlation between Hildebrand solubility parameter (δt) of (COs)

and δt for MCO(Asp-free) is linear with a slope close to unity. Regarding Asp dissolution, we show that the correlations between log (dissolved Asp, mass%) and each of the following solvent descriptors is linear: empirical polarity of MCOs(Asp-free); (δt) of COs; colloidal instability index of COs. Furthermore, the multiple correlation between log (dissolved Asp, mass-%) and other solvatochromic parameters showed that solvent dipolarity and ACS Paragon Plus Environment

Energy & Fuels

[2] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

polarizability are important factors for Asp dissolution, in agreement with our previous results on Asp dissolution in pure solvents. The formulation of a model that successfully mimics maltenes is potentially very useful, e.g., in rationalizing the efficiency of certain classes of additives employed for Asp stabilization.

Abbreviations and acronyms

Asp:

Asphaltene

CII

Colloidal instability index of a crude oil.

CO:

Crude oil

CO(Asp-free)

Asphaltene-free crude oil (or maltene).

ET(HxQMBu2)

Empirical solvent polarity in kcal/mol, based on HxQMBu2.

HxQMBu2

(E-2,6-di-tert-butyl-4-(2-(1-hexylquinolin-1-ium-4-

yl)vinyl)phenolate), a solvatochromic probe for calculation of the empirical solvent polarity parameter ET(probe). MCO(Asp-free)

Formulated model for Asphaltene-free crude oil

MCO(Asp&Resin-free)

Formulated model for Asphaltene- and resin free crude oil

MNI

1-Methyl-5-nitroindoline, a solvatochromic probe for calculation of

solvent Lewis basicity. NI

5-Nitroindoline, a solvatochromic probe for calculation of solvent

Lewis basicity r2

Regression correlation coefficient

SARA

Saturates, aromatics, resins, and asphaltenes of a crude oil.

ΣQ2

Sum of the squares of the residuals.

Introduction The potential problems caused by flocculation followed by sedimentation of asphaltenes (Asps) during petroleum extraction, transport and refining are the impetus for the sustained interest in studying different aspects of Asps, inter alia, experimental methods for determination of Asp stability,1-5 determination of Asp structure from different crude oils (COs),6-11 introduction of models for Asps, 12-18 and theories regarding factors that contribute to their stability, etc.19-25. The complexity of Asp structure, the ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

Energy & Fuels

[3] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

multiplicity of the interactions of Asps with maltenes (complex mixture of hydrocarbons plus resins, hereafter designated CO(Asp-free)) mean that a general theory to evaluate the importance of Asp-maltene interactions to Asp stability is a distant goal. Instead, we hope to identify general guidelines that serve, e.g., in predicting the stability of COs during their extraction, transport, refining, and choose between different additives to enhance Asp stability, where required. Several scales were introduced to evaluate Asp stability, including the colloidal instability index, CII;1,26 the Heithaus (P) parameter;27-30 the toluene-equivalence (TE) parameter and the related Bureau of Mines correlation index.

31,32

Although the

experimental parts of these determinations and the subsequent calculations are relatively simple, their empirical nature does not offer much insight into the delicate balance of forces whose perturbation leads to flocculation and sedimentation of Asps.22 One of the strategies that are employed for a better understanding Asp stability is to use the oil compatibility model. The latter rests on the assumption that Asp stability is assured when Asps (the solute) and maltenes (the solvent) are “compatible”.

33,34

The

problem is reduced to finding an appropriate scale for quantifying the “solvent power” of maltenes, and then test a collection of maltenes to evaluate the success/generality of this model. A distinct merit of this approach is that it addresses the complex problem of Asp stability in terms of solute/solvent interactions on the molecular level. This is done by using solvent descriptors whose values are experimentally accessible. Examples of these descriptors are Hildebrand and Hansen solubility parameters, recently showed- solvatochromic parameters.

40,41

35-39

and- as we

The term solvatochromism refers to

the effect of the solvent, e.g., maltenes on the Uv-Vis spectra, absorption or emission of solvatochromic substances (hereafter designated as “probes”). There is an intramolecular charge-transfer within the probe whose strength (hence the value of λmax) is affected by probe-solvent interactions. We use this solvent effect on the spectrum of the probe to calculate the corresponding solvatochromic parameter. For example, we calculate the empirical (or total) polarity of the solvent by using the probe HxQMBu2 and the equation:

ET(HxQMBu2) (kcal/mol) = 28591.5/λmax (nm)

ACS Paragon Plus Environment

(1)

Energy & Fuels

Page 4 of 30

[4] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

where λmax is the wavelength of the intramolecular charge transfer within the probe, i.e., that from the phenolate oxygen to the quaternary nitrogen. Figure 1 shows the structure of the probes of interest in the present work.

Figure 1: Molecular structures of the solvatochromic probes employed in the present work: (a) (E-2,6-di-tert-butyl-4-(2-(1-hexylquinolin-1-ium-4-yl)vinyl)phenolate, HxQMBu2), (b) β-carotene, (c) 2-(N,N-dimethylamino)-7-nitro- 9H-fluorene (DMANF) (d) 5nitroindoline (NI) and (e) 1-methyl-5-nitroindoline (MNI).

Note that ET(probe) is the sum of specific (e.g., hydrogen-bonding) and nonspecific (e.g., dispersion forces) probe-solvent interactions as shown by Eq 2: 42

ET(probe) = a SA + b SB + d SD + p SP

(2)

where (S) refers to solvent (or mixture of solvents) whose Lewis acidity, Lewis basicity, dipolarity, and polarizability are given by the solvent descriptors (A, B, D, P) respectively and (a, b, d, p) are the corresponding regression coefficients. Eq 3 is the corresponding version of Hansen´s solvent solubility parameters:

δt2 = δAB2 + δD2 + δP2

(3)

The terms AB, D and P refer to (combined) hydrogen bond donation (Lewis acidity) and acceptance (Lewis basicity), dipolar, and dispersion properties of the solvent,

ACS Paragon Plus Environment

Page 5 of 30

Energy & Fuels

[5] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

respectively. The symbols that we use in Eq 3 are different from those originally suggested by Hansen; 36 we use them here for consistency with Eq 2. Employing a set of 38 molecular solvents we recently showed that the values of ET(HxQMBu2) correlate linearly with the corresponding Hildebrand (total) solubility parameters (δt). Likewise, the solvatochromic-based descriptors for solvent Lewis acidity/basicity, dipolarity, and polarizability correlated linearly with the corresponding Hansen solubility parameters. Therefore, we showed that solvatochromic- and solubility parameters are equivalent in accounting for solvent effects.41 Compared with the solubility parameters, the solvatochromic counterparts have the merit of simplicity of the experimental part, and that we treat Lewis acidity and Lewis basicity separately. We obtain information regarding the medium effect on a phenomenon of interest (e.g., Asp dissolution) by applying equations 2 or 3 to the phenomenon in several solvents and examining the sign -positive or negative- and magnitude of the regression coefficients. For example, application of Eq 2 to the reaction of 1-methylimidazole with 1bromohexane in 10 solvents, at 40 °C resulted in the following regression coefficients: a = - 3.79; b = + 20.89; d/p = +56.36. This means that solvent Lewis basicity and dipolarity/polarizability enhances this SN2 reaction; the latter solvent descriptor is more important. On the other hand, solvent Lewis acidity decreases the reaction rate, probably because of diazole solvation by hydrogen-bonding.

43

Equation 3 was applied to

determine the contribution of these parameters to the stability of six COs. 44,45 Likewise, we correlated both scales (solvatochromic and Hansen) with the solubility of Asps in 28 molecular solvents by using Eq 4 and concluded that solvent dipolarity and polarizability are the relevant solvent descriptors. 41

log (dissolved Asps; mass-%) = a/b (SA + SB) + d SD + p SP

(4)

In the present work, we extended the above-mentioned approach to maltenes. Specifically, we wanted to probe the effect of CO(Asp-free) composition on the dissolution/stability of the corresponding Asps, and to correlate this dissolution with solvent descriptors. The dark colors of the latter fractions precluded recording the Uv-Vis spectra of the probes. Therefore we formulated models for maltenes, MCO(Asp-free) whose composition was based on SARA analysis (Saturate, Aromatic, Resin, and Asphaltene) of ACS Paragon Plus Environment

Energy & Fuels

Page 6 of 30

[6] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the precursor COs, and elemental analysis (C/H/N/S) of the corresponding resins. The components of MCO(Asp-free) were mixtures of isooctane (2,2,4-trimethylpentane), cis and trans decalines, 1-methylnaphthalene, benzothiazole and (n-octyl)-1-naphthoate. First we validated the use of this MCO(Asp-free). We then showed that the solubility of Asps in MCO(Asp-free) correlate linearly with: experimental ET(HxQMBu2) of MCO(Asp-free); with experimental δt of CO and SARA-based CII. Additionally, multiple regression analysis with the solvatochromic parameters listed in Eq. 2 indicated that dipolarity and polarizability are important to the dissolution of Asps, in agreement with our previous investigation on solubility of Asps in pure solvents. The relevance of our results to inhibition of Asp precipitation is briefly discussed.

2- Experimental Section 2.1-Equipment We used Fritsch, Laborette 17 model, operating at 35 kHz for sample sonication. The refractive indices (n ± 0.00001) of the crude oils and mixtures were determined with Rudolph Research J357 digital refractometer, operating at 488 nm, all measurements were recorded at 25 οC. We used Bruker Vector-22 and Varian Innova-300 spectrometers to obtain IR and 1H NMR spectra, respectively. Elemental analysis was carried out on Perkin Elmer 2400 series II CHN analyzer at the Analytical Facilities of the Institute of Chemistry, USP, São Paulo.

2.2-Chemicals We purchased the chemicals from Acros, Merck, or Synth (São Paulo) and purified them as recommended elsewhere. 46

2.3- Crude oils and their fractions The crude oil samples were from Petrobras offshore production fields in the Campos Basin, Brazil. Their SARA analyses were supplied by Petrobras, using chromatographic analysis (ASTM D5186)

47

and the IP-143/84 method.

48,49

Results of these analyses are

listed in Table 1 along with the calculated, SARA-based CII of the COs. The corresponding

ACS Paragon Plus Environment

Page 7 of 30

Energy & Fuels

[7] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

asphalt-free and resin-free fractions (CO(Asp&Res-free)) were supplied by Petrobras and were obtained using high pressure fractionation unit. 49

Table 1: Crude oil designation; results of SARA analysis (mass %), and calculated colloidal instability indexes, CII.

Crude oil

Saturate

Aromatic

Resin

Asphaltene

CII

PB-359

42.2

28.9

27.3

1.6

0.78

PB-574

44.2

34

21.3

0.5

0.81

PB-758

58.6

23.4

18

0.5

1.43

PB-1069

66.4

21.0

11.3

1.3

2.10

PB-1081

41.8

35.2

22.3

0.75

0.74

PB-1091

64.4

19.8

15.1

0.7

1.87

PB-1093

52.3

24.5

22.5

0.70

1.13

PB-1100

48.6

21.3

28.6

1.2

1.00

PB-1154

77.6

16

6.4

0.5

3.49

PB-1157

28.0

35.3

33.6

3.1

0.45

designation

An asphaltene-rich solid was obtained from the deasphalting unit of Petrobras. The Asps were obtained from this residue by repeated extraction with n-heptane as described elsewhere.41 Briefly, 200 g Asp-rich solid were suspended in 400 mL n-heptane, sonicated for 4h, the n-heptane separated by decantation; this step was repeated 3 times. The residual solid was the extracted with n-heptane for 32 h (Soxhlet), and dried under reduced pressure until constant mass. Elemental analysis: C, 83.05%; H, 8.60%; N, 1.64%; and S, 5.35%. The H/C (atomic) ratio is 1.23, which is in the same range as that found for Asphs extracted from crude samples, e.g., from Brazil (1.18),49 Mexico (1.14),50 Venezuela (1.10-1.20),51 and USA (0.95-1.14).52

2.4- Synthesis of n-octyl 1-naphthoate ACS Paragon Plus Environment

Energy & Fuels

[8] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

We prepared this ester as follows: 1-naphthoic acid (12.00 g; 69.7 mmol) was reacted with thionyl chloride (150 ml, 2.1 mol) under reflux, 3h, followed by removing excess thionyl chloride. Under cooling (ice-bath) the crude 1-naphthoyl chloride (15.58 g, 81.8 mM) was slowly added to a mixture of 1-octanol (9.97g, 76.6 mmol), triethylamine (8.39 g, 83 mmol) and dicloromethane (180 mL). The reaction mixture was stirred for two hours at 0oC, and then overnight at room temperature. Hexane (50 mL) was added, the white precipitate (triethylammonium chloride) filtered, and the volatiles removed. Dichloromethane (50 mL) was added and the resulting yellowish solution was washed with dilute HCl, dilute NaOH, water and dried (anhydrous MgSO4). The solvent was evaporated, and the residual oil was distilled (b.p. 240°C at 1.0 mmHg). Yield “75 %” of a pale buff color liquid. IR data, neat sample (frequency in cm-1, band attribution): 3087 (aromatic C-H), 2928 (aliphatic C-H); 1716 (C=O); 1244 (C-O). 53 The structure and hydrogen atom numbering of this ester is shown in Figure SI-1 . 1H NMR (CDCl3) n-octyl 1-naphthoate: Hg 8.99 (1H, d, J=8.6Hz); Ha 8.22 (1H, dd, J=8.6 and J=1.3Hz); Hc 8.01 (1H, d, J= 8.3Hz); Hd 7.88 (1H, d, J=8.1 Hz); Hf 7.63 (1H, td, J= 8.6 and 8.1Hz); Hb and He 7.47-7,56 (2H, overlapping peaks); Hh 4.44 (2H, t, J=6Hz); Hi 1.85 (2H, m); Hj 1.32-1.46 (10H, overlapping peaks); Hk 0.93 (3H, t, J=7Hz). 54,55 Elemental analysis: Calculated for C19H24O2: C% (80.24); H% (8.51), analyzed C% (80.06); H% (8.65).

2.5- Formulation of models for asphaltene and resin-free crude oils MCO(Asp&Res-free) and asphaltene-free crude oils (MCO(Asp-free)) Each MCO(Asp&Res-free) was formulated by mixing decalines, isooctane and 1methylnaphtalene as given in Table SI-1. Samples of MCO(Asp-free) were obtained from the same hydrocarbons plus a model for the resins, made of equimolar mixture of benzothiazole and n-octyl-1-naphthoate, see Table SI-2.

2.6- Spectroscopic Determination of the solvatochromic parameters ET(HxQMBu2), SB, SD and SP. We used the probes shown in Figure 1 to determine the following solvatochromic parameters: HxQMBu2, ET( HxQMBu2); NI and MNI (SB), DMANF (SD) and β-carotene (SP).41,56-58 Aliquots of the probe solution in acetone were pipetted into small glass vials, followed by evaporation of acetone at room temperature, under reduced pressure, in the ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30

Energy & Fuels

[9] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

presence of P4O10. The solvent mixture whose solvatochromic parameter is to be determined was added (1mL), the probe was dissolved, the solution transferred into a cuvette with appropriate path-length (0.5-2 cm); see experimental conditions in Table SI3. The absorbance of the probes was recorded at 25 °C, 120 nm/min, using Shimadzu UV 2550 spectrophotometer in the range of 350-550nm for β-carotene, DMANF, NI and MNI and 500-700nm for HxQMBu2. Values of λmax were calculated from the first derivative of the absorption spectra; values of ET(HxQMBu2) in kcal/mol were calculated from Eq 1. The equations employed to calculate SB, SD and SP from the values of λmax of the solvatochromic probes are given in the Calculations section of SI. Calculated values of SB, SD, and SP of MCO(Asp-free) are listed in Table-4, vide infra. We also calculated ET(HxQMBu2) of CO(Asp&Res-free). Due to the dark color of these samples, we proceeded as follows: solid HxQMBu2 (4 mg) was added to 2.5 mL of CO(Asp&Res-free), containing 10 mg NaHCO3 (to neutralize any acid traces present). The mixture was sonicated (1h), centrifuged at 4400g for 30 min, the supernatant filtered (PTFE filter, 0.45 μm membrane), and the solution absorbance was recorded as indicated above, using 0.1 cm path length cell, and CO(Asp&Res free) in 0.1 cm cell in the reference beam.

2.7- Determination of Hildebrand solubility parameters (δt) of Crude oils from refractive index measurements We prepared several mixtures (10) each containing 0.5 g of (CO) and variable volumes (0.5 to 1.75 mL) of n-heptane. Each mixture was homogenized (vortex mixer) during 2 minutes, left overnight, and the refractive index (n) of the supernatant determined. Each experiment was carried out in triplicate and the average value of (n) was calculated. Values of (δt) were calculated from the values of (n) as given in the Calculations section of SI. The correlation between the refractive index of CO and volume of n-heptane/g of crude oil is shown in Figure SI-2.

2.8- Determination of asphaltene solubility in MCO(Asp-free): The models prepared in item 2.5 were employed as follow: we weighted ca 350 mg of asphaltenes (see item 2.3 above) in 2 mL Eppendorf polypropylene tube; added 1 mL of MCO(Asp-free); agitated the tube at room temperature, first for 2 min using a vortex ACS Paragon Plus Environment

Energy & Fuels

Page 10 of 30

[10] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

mixer, then for 5 h using a tube rotator (Glas-Col model 099A RD4512, 60 rpm), left the tube overnight at room temperature. The tubes were centrifuged for 1 h at 16400 g (Incibrás Spin 1), and the supernatant was filtered through a PTFE membrane (0.45 μm). A certain mass of the filtered solution (ca. 0.3- 0.5 g) was dried in a vacuum oven until constant mass (T = 80 °C; 12-36 h). The solubility of dissolved Asp was calculated in mass%, i.e., dissolved Asp in g/100 g of the filtered asphaltene solution. Each experiment was carried out at least in duplicate.

3- Results and Discussion Note: calculations of all solvent parameters are given in the Calculations section of supporting Information

3.1-Information sought from use of the solvatochromic probe HxQMBu2 As shown in Figure 1, HxQMBu2 possess a dinuclear heterocyclic ring, branched and normal alkyl chains, and two sites for specific and non-specific interactions, namely the phenolate oxygen and the quaternary nitrogen, respectively.

59

Therefore, this probe is

expected to be sensitive to precisely the same (specific and nonspecific) interactions that affect the stability of Asps.

60-63

In other words, we hoped to gain information on the

quality of maltenes as solvents for Asps from the solvatochromic response of this probe.

3.2-The need for- and suitability of model for CO(Asp-free) As mentioned above, we decided to test whether the dissolution of Asps in CO(Aspfree) can

be correlated with the quality of the latter as solvent, as quantified by: (i) “global”

solvent descriptors, including ET(HxQMBu2) and δt; (ii) individual solvent descriptors, e.g., those shown in Eq 4 , with CO(Asp-free) substituting (S). Application of this equation requires that we have the solvatochromic data for CO(Asp-free). This was not feasible because (dark colored) CO(Asp-free) absorbed strongly in the same Uv-Vis spectral region of the probes shown in Figure 1. This turned determination of these solvatochromic parameters for unfeasible. Facing this problem, we decided to develop mixtures that mimic asphaltene-free crude oils, here after designated MCO(Asp-free). We based such formulation essentially on the results of SARA analysis of the corresponding COs, and

ACS Paragon Plus Environment

Page 11 of 30

Energy & Fuels

[11] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

elemental analysis (C/H/N/S) of the corresponding resins. We discuss below this formulation, and give evidence for its suitability as model for CO(Asp-free). We based the model on the assumption that maltenes are mixtures of simple, low molar mass hydrocarbons plus resins. As model for the former, we employed mixtures of 1-methyl naphthalene, isooctane and decalines as aromatic and aliphatic components, respectively. As model for the resins we used an equimolar mixture of benzothiazole and n-octyl-1-naphoate. Each MCO(Asp-free) sample was prepared by mass so that is has the same (SARA) component of the SARA analysis of the corresponding crude oil, see Tables 1 and SI-2. Regarding this formulation, the following is relevant: (i)- There are several examples in the literature for use of mixtures of low molar mass hydrocarbons as models for CO(Asp&Resin-free), including mixtures of n-decane, cyclohexane, toluene, and dioxane;

64

and mixtures of n-heptane, n-octadecane, n-tricosane plus a

petroleum fraction (b.p. = 230–400°C); 65 (ii) Table-2 shows the elemental analysis results for petroleum resins, and some suggested resin models.

Table-2: Elemental analysis of some petroleum resins and models for resins Entry

Country of crude oil origin

%N

%S

%O

H/C atom

Refer-

ratio

ence

1.26-1.49

52

1.31-1.51

66

1.41-1.48

67

1.36-1.39

68

1.39

7

Crude oil-based resins 1

2

3

4

5

USA

USA; Saudi Arabia; Argentina

Canada

USA, Saudi Arabia

Saudi Arabia

0.87-

1.76-

1.94-

1.38

3.36

8.04

0.81-

0.88-

1.53-

1.52

6.91

2.77

0.5-

4.8-

2.7-

0.9

6.0

4.6

0.66-

1.9-

2.3-

1.89

5.95

2.9

0.78

5.02

2.25

Models for resins

ACS Paragon Plus Environment

Energy & Fuels

Page 12 of 30

[12] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Structure of resin component(s)

%N

%S

%O

H/C

Reference

6

3.3

7.7

7.6

1.16

Presen t work

7

3.2

14.8

--

1.58

69, 70

8

--

16.2

--

0.77

71

As shown in Table 2, the %N of our model (entry 6) is practically the same as the model given in entry 7. Likewise, the %S of our model is not far from crude oil-based resins (entries 3-5) and is less than that in entries 7-8. In variance with the last entries, our model contains oxygen; its H/C atom ratio is close to the range given in entries 1-5. Therefore, elemental analysis of our resin model is in the range of resins from several distinct crude oils, as well as other suggested resin models.

3.3-Validation of our model for CO(Asp-free) For convenience, we discuss this validation by considering first MCO(Asp&Resin-free) because this is a fraction of MCO(Asp-free), then we address MCO(Asp-free) . As indicated above, we were able to record the Uv-Vis spectra of HxQMBu2 both in CO(Asp&Resin-free) and MCO(Asp&Resin-free) hence calculate the corresponding values of ET(HxQMBu2). Based on the data listed in Table 3, a plot of ET(HxQMBu2) for CO(Asp&Resin-free) versus ET(HxQMBu2) for MCO(Asp&Resin-free) gave a straight line with a slope of 0.98 and correlation coefficient 0.93, see Figure 2.

ACS Paragon Plus Environment

Page 13 of 30

Energy & Fuels

[13] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2: Correlation between ET(HxQMBu2) for CO(Asp&Resin-free) and ET(HxQMBu2) for MCO(Asp&Resin-free); slope 0.98, r2= 0.93. Values of ET(HxQMBu2) are in kcal/mol.

This shows that the probe-solvent interactions are practically the same in both solvents, i.e., MCO(Asp&Resin-free) mimics CO(Asp&Resin-free).

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

Page 14 of 30

[14]

Table 3. Crude oil designation, properties of crude oil fractions, properties of models for crude oil fractions, and mass-% dissolved asphaltenes in MCO(Asp-free).(a) Crude oil

Properties of crude oil fractions

Properties of models for crude oil fractions

designation

Dissolved Asphaltenes, mass-% in MCO(Asp-free)

δ CO

ET(HxQMBu2) of

ET(HxQMBu2) of

ET(HxQMBu2) of

(MPa0.5)

CO(Asp&Resin- free)

MCO(Asp&Resin-free)

MCO(Asp-free)

PB-359

19.0

47.0

47.2

PB-574

18.9

47.4

47.8

PB-758

18.4

47.7

47.7

PB-1069

17.9

47.3

47.4

PB-1081

18.9

47.6

47.4

PB-1091

18.2

47.6

47.5

PB-1093

18.6

48.0

48.1

PB-1100

18.6

47.8

47.8

PB-1154

17.3

48.1

48.2

PB-1157

18.4

47.5

47.5

(a)ET(HxQMBu2) values in kcal/mol.

ACS Paragon Plus Environment

46.6 46.9 47.1 47.3 46.7 47.3 46.8 46.8 47.9 46.4

25.56 23.84 27.06 16.82 30.06 21.20 22.55 24.44 6.88 34.95

Page 15 of 30

Energy & Fuels

[15] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Due to the above-mentioned problem of the dark color of CO fractions we correlated

δt(CO) versus δtMCO(Asp-free).The former was calculated from the onset of asphaltene precipitation from the crude oil with n-heptane (Table 3), whereas the latter was calculated from δt of its components as given in the Calculations section of SI. This correlation gave a straight line with slope = 1.28 and correlation coefficient = 0.82 (Figure SI-3). Considering these results, we satisfactorily validated the use of MCO(Asp-free) as model for maltenes.

3.4-Can we evaluate the power of asphaltene dissolution by hydrocarbon mixtures in terms of their solvent descriptors? At the outset we comment on Asp dissolution because any particle with diameter < 0.45 µ (the average pore size of the PTFE membrane employed) will pass through the filter and considered, therefore, as dissolved Asp, although it is certainly a nanoaggregate, see, e.g., the discussion on the Yen−Mullins Asp aggregation model.72 Additionally, there are kinetic effects associated with the precipitation and aggregation of Asps; reaching thermodynamic equilibrium may require days in some cases.

73

Considering these aspects we decided to use a filtration membrane with small pores instead of filter paper,74 agitated the Asp suspension for 5h, left it overnight and centrifuged it at the same g-force employed by Fogler et al.73 In summary, we made effort to ensure consistent Asp dissolution data. We show below the obtained correlations between mass-% dissolved Asps and solvent descriptors; these are shown in Table-4. The plot of Figure-3 is based on empirical solvent polarity of MCO(Asp-free); Figure-4 is that based on the solubility parameters of the crude oil, whereas we use CII of (CO) as descriptor in Figure-5. The latter property was calculated from SARA analysis, using Eq. (5), as given elsewhere:75

CII= (% saturates + % asphaltenes) / (% aromatics + % resins)

ACS Paragon Plus Environment

(5)

Energy & Fuels

Page 16 of 30

[16] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 4: Calculated solvatochromic parameters for models of asphalt-free crudes, MCO(Asp-free) , at 25 οC.

MCO(Asp-free) formulated

Solvatochromic Parameters

based on SARA of SB

SD

SP

PB-359

0.211

0.586

0.822

PB-574

0.494

0.475

0.756

PB-758

0.505

0.625

0.788

PB-1069

0.473

0.281

0.773

PB-1081

0.430

0.519

0.799

PB-1091

0.577

0.546

0.761

PB-1093

0.510

0.410

0.794

PB-1100

0.515

0.461

0.790

PB-1154

0.269

0.081

0.717

PB-1157

0.428

0.519

0.795

All plots are straight lines with good correlation coefficients. Figures-3 and Figure-5 have

negative

slopes

for the

following reasons: HxQMBu2

shows

positive

solvatochromism, i.e, the value of λmax increases as a function of increasing solvent polarity. As shown by Eq 1, more polar solvents have smaller values of ET(HxQMBu2), leading to the negative slope shown in Figure-3. From the definition of CII, solvents with high aromatic content, i.e., efficient solvents for Asp dissolution have smaller CII. The conclusion from Figure-3 to Figure-5 is the following: more polar solvents dissolve Asps better; this efficiency is directly related to experimentally accessible solvent descriptors.

ACS Paragon Plus Environment

Page 17 of 30

Energy & Fuels

[17] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure-3: Correlation of ET(HxQMBu2) (in kcal/mol) of MCO(asp-free) and log(dissolved asphaltene, mass-%) in MCO(asp free). r2 = 0.84. (N=10)

Figure-4: Correlation of Hildebrand solubility parameter of CO and log(dissolved asphaltene, mass-%) in MCO(asp free).; r2 = 0.84 (N=9)

Figure-5: Correlation between colloidal instability index (CII) of CO(Asp-free) and log(dissolved asphaltene, mass-%) in MCO(Asp free); r2 = 0.92 (N=10) ACS Paragon Plus Environment

Energy & Fuels

Page 18 of 30

[18] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

We took this analysis a step further by calculating the dependence of log (dissolved Asp, mass-%) on the solvatochromic parameters of MCO(Asp-free), where recording the UvVis spectra of all probes shown in Figure 1 is feasible. The Eq employed was:

log (dissolved Asps; mass-%) = b SB + d SD + p SP

(6)

That is, we considered that solvent Lewis acidity is not contributing; Lewis basicity is operative because of the presence of heteroatoms in the resin model, SD and SP are- as expected- contributing. In order to compare the regression coefficients of Eq (6) directly, we used reduced solvatochromic parameters (i.e., on a scale that varies between 0 and 1); Eq (7) was obtained:

log (dissolved Asps, mass-%) = (5.84 ± 0.9) + (0.13 ± 0.10) SBred + (0.31 ± 0.15) SDred + (0.33 ± 0.16) SPred , r2: 0.885 ΣQ2 = 0.081

(7)

The regression coefficients of all parameters are positive, i.e., basic, dipolar, and polarizable maltenes stabilize Asps; dipolarity and polarizability are important solvent properties and agreement with previously literature.41,44,76,33

4-Conclusions The compatibility model is very useful in predicting the stability of Asps. We have successfully used solvatochromic parameters, Hildebrand/Hansen solubility parameters and CII as convenient “scales” for assessing the power of maltenes to dissolve Asps. We formulated models for maltenes by mixing simple hydrocarbons (cis and trans-decalines, isooctane and 1-methylnaphthalene) and resin components (benzothiazole and n-octyl-1naphthoate). We based this formulation on SARA analysis of the crude oils and elemental analysis (C/H/N/S) of the corresponding resins. We validated these formulations as models for maltenes. Correlations of Asp solubility with the descriptors of (CO) and MCO(Asp-free) showed that solvent dipolarity and polarizability are important factors that contributes to Asp dissolution. Note that the latter descriptor includes the stacking of aromatic rings that contributes to Asp stability and Asp-maltene interactions.33 The ACS Paragon Plus Environment

Page 19 of 30

Energy & Fuels

[19] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

formulation of a model that successfully mimics maltenes is potentially very useful, e.g., in rationalizing the efficiency of certain classes of additives employed for Asp stabilization.

Acknowledgements L. P. Novaki, R. Lira, M. M. N. Kwon thank Petrobras for research fellowships. We thank the Institute of Chemistry of the University of São Paulo for making its research facilities available to us. O. A. El Seoud thanks FAPESP (grant 2014/22136) and CNPq for research productivity fellowship (307022/2014-5).

5- References (1) Asomaning, S. Test method for determining asphaltene stability in crude oils. Pet. Sci. Technol. 2003, 21, 581–590. (2) Ocanto, O.; Marcano, F.; Castillo, J.; Fernandez, A.; Caetano, M.; Chirinos, J.; Ranaudo, M. A. Influence of Experimental Parameters on the Determination of Asphaltenes Flocculation Onset by the Titration Method. Energy Fuels 2009, 23, 3039–3044. (3) Gonzalez, D. L.; Vargas, F. M.; Mahmoodaghdam, E.; Lim, F.; Joshi, N. Asphaltene Stability Prediction Based on Dead Oil Properties: Experimental Evaluation. Energy Fuels 2012, 26, 6218−6227. (4) Chamkalani, A.; Mohammadib, A. H.; Eslamimanesh, A.; Gharagheizi, F.; Richon, D. Diagnosis of asphaltene stability in crude oil through ‘‘two parameters’’ SVM model. Chem. Eng. Sci. 2012, 81, 202-208. (5) Barskaya, E. E.; Yusupova, T. N.; Saraev, D. V.; Criterion of phase stability of asphaltenes in crude oils. Chem. Tech. Fuels Oils 2013, 49 (2), 165-174. (6) Strausz, O. P.; Mojelsky, T. W.; Lown, E. M. The molecular structure of asphaltene: an unfolding story. Fuel 1992, 71, 1355. (7) Verstraete, J. J.; Schnongs, Ph.; Dulot, H.; Hudebine, D. Molecular reconstruction of heavy petroleum residue fractions. Chem. Eng. Sci. 2010. 65. 304-312. (8) Liu, Y.-J.; Li, Z.-F. Structural Characterisation of Asphaltenes during Residue Hydrotreatment with Light Cycle Oil as an Additive. J. Chem. 2015, Article ID 580950.

ACS Paragon Plus Environment

Energy & Fuels

[20] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(9) AlHumaidan, F.S.; Hauser, A.; Rana, M.S.; Lababidi, H. M. S. NMR Characterization of Asphaltene Derived from Residual Oils and Their Thermal Decomposition. Energy Fuels 2017, 31 (4), 3812–3820. (10) Molina V., D. R.; Leon, E. A.; Poveda-Jaramillo, J. C. Structural differences among the asphaltenes in Colombian light crudes from the Colorado Oil Field. Energy Fuels 2017, 31, 133-139. (11) Evdokimov, I. N.; Fesan, A. A; Losev, A. P. Occlusion of Foreign Molecules in Primary Asphaltene Aggregates from NUV-Visible Absorption Studies. Energy Fuels 2017, 31 (2), 1370-1375. (12) Boek, E. S.; Yakovlev, D. S.;. Headen, T. F.; Quantitative Molecular Representation of Asphaltenes and Molecular Dynamics Simulation of Their Aggregation. Energy Fuels 2009, 23, 1209–1219. (13) Aray, Y.; Hernandez-Bravo, R.; Parra, J. G.; Rodríguez, J.; Coll., D. S. Exploring the Structure Solubility Relationship of Asphaltene Models in Toluene, Heptane, and Amphiphiles Using a Molecular Dynamic Atomistic Methodology. J. Phys. Chem. A 2011, 115, 11495–11507. (14) Scherer, A.; Hampel, F.; Gray, M. R.; Stryker, J. M.; Tykwinski, R. R. Incorporation of steroidal biomarkers into petroleum model compounds. J. Phys. Org. Chem. 2012, 25, 597–606. (15) Costa, L. M.; Stoyanov, S. R.; Gusarov, S.; Tan, X.; Gray, M. R.; Stryker, J. M.; Tykwinski, R.; Carneiro, J. W. M.; Seidl, P. R.; Kovalenko, A. Density Functional Theory Investigation of the Contributions of π−π Stacking and Hydrogen-Bonding Interactions to the Aggregation of Model Asphaltene Compounds. Energy Fuels 2012, 26, 2727−2735. (16) Breure, B.; Subramanian, D.; Leys, J.; Peters, C. J.; Anisimov, M. A. Modeling Asphaltene Aggregation with a Single Compound. Energy Fuels 2013, 27, 172−176. (17) Li, D. D.; Greenfield, M. L.; Chemical compositions of improved model asphalt systems for molecular simulations. Fuel 2014, 115, 347–356. (18) Sjöblom, J.; Simon, S.; Xu, Z. Model molecules mimicking asphaltenes. Adv. Colloid Interface Sci. 2015, 218, 1–16. (19) Leo´n, O.; Rogel, E.; Espidel, J.; Torres, G. Asphaltenes: Structural Characterization. Self-Association and Stability Behavior. Energy Fuels 2000, 14, 6-10. (20) Pereira, J. C.; López, I.; Salas, R.; Silva., F.; Fernandez, C.; Urbina, C.; López, J. C.; Resins: The Molecules Responsible for the Stability/Instability Phenomena of Asphaltenes. Energy Fuels 2007, 21, 1317-1321. (21) Likhatsky, V.V.; Syunyaev, R.Z. New Colloidal Stability Index for Crude Oils Based on Polarity of Crude Oil Components. Energy Fuels 2010, 24(12), 6483-6488. ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30

Energy & Fuels

[21] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(22) Moura, L.G.M.; Santos, M.F.P.; Zilio, E.L.; Rolemberg, M.P.; Ramos, A.C.S. Evaluation of indices and of models applied to the prediction of the stability of crude oils. J Petroleum Sci. Eng. 2010, 74, 77-87. (23) Castellano, O.; Gimon, R.; Soscun, H. Theoretical Study of the σ_π and π_π Interactions in Heteroaromatic Monocyclic Molecular Complexes of Benzene. Pyridine. and Thiophene Dimers: Implications on the Resin_Asphaltene Stability in Crude Oil. Energy Fuels 2011, 25, 2526–2541. (24) Abdel-Raouf, M. E.-S. Factors Affecting the Stability of Crude Oil Emulsions; Crude Oil Emulsions: Composition Stability and Characterization 2012, 183-204. Publisher: (InTech. Rijeka. Croatia) (25) Marques, L. C. C.; Pereira, J. O.; Bueno, A. D.; Marques, V. S.; Lucas, E. F.; Mansur, C. R. E.; Machado, A. L. C.; González, G. A Study of Asphaltene-Resin Interactions. J. Braz. Chem. Soc. 2012, 23 (10), 1880-1888. (26) Wiehe, I. A.; Kennedy, R. J. The Oil Compatibility Model and Crude Oil Incompatibility. Energy Fuels 2000, 14, 56-59. (27) Heithaus, J.J. Measurement and significance of asphaltene peptization. J. Inst. Pet. 1962, 48, 45–53. (28) ASTM: D-6703-14. Standard test method for automated Heithaus titrimetry. (29) Tojima, M.; Suhara, S.; Imamura, M.; Furuta, A. Effect of heavy asphaltene on stability of residual oil. Catal. Today 1998, 43, 347–351. (30) Schermer, W.E.M.; Melein, P.M.J.; van den Berg, F.G.A. Simple techniques for evaluation of crude oil compatibility. Pet. Sci. Technol. 2004, 22, 1045–1054. ASTM-D 6703-01. (31) Griffith, M.G.; Siegmund, C.W. Controlling of Residual Fuel Oils. ASTM – Marine Fuels Symposium, Miami, 1983. (32) Stor, L.M. Desenvolvimento de metodologia para previsão da compatibilidade de misturas de petróleo. Mestrado em Processos Industriais. Instituto de Pesquisas Tecnológicas do Estado de São Paulo (IPT), 2006 (in Portuguese). (33) Wiehe, I. R. Asphaltene Solubility and Fluid Compatibility. Energy Fuels 2012, 26, 4004−4016. (34) Ramos, A. C. S.; Rolemberg, M. P.; Moura, L. G. M.; Zilio, E. L.; Santos, M. F. P.; Gonzalez, G. Determination of solubility parameters of oils and prediction of oil compatibility. J. Petroleum Sci. Eng. 2013, 102, 36–40.

ACS Paragon Plus Environment

Energy & Fuels

[22] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(35) Redelius, P. Bitumen Solubility Model Using Hansen Solubility Parameter. Energy Fuels 2004, 18, 1087-1092. (36) Hansen, C. M. Hansen Solubility Parameters: A User´s Handbook, 2nd Edition, CRC Press. Boca Raton, 2007. (37) Acevedo, S.; Castro, A.; Vasquez, E.; Marcano, F.; Ranaudo, M. A. Investigation of Physical Chemistry Properties of Asphaltenes Using Solubility Parameters of Asphaltenes and Their Fractions A1 and A2. Energy Fuels 2010, 24, 5921–5933. DOI:10.1021/ef1005786. (38) Moura, L. G. M.; Rolemberg, M. P.; Ramos, A. C. S. Avaliação das incertezas associadas à determinação do parâmetro de solubilidade de Hildebrand de petróleos. Quim. Nova 2011, 34(2), 226-231. (39) Sato, T.; Araki, S.; Morimoto, M.; Tanaka, R.; Yamamotom H. Comparison of Hansen Solubility Parameter of Asphaltenes Extracted from Bitumen Produced in Different Geographical Regions. Energy Fuels 2014, 28, 891−897. (40) Reichardt, C.; Welton, T. Solvents and solvent effects in organic chemistry. 4th ed. Weinheim: Wiley-VCH; 2011. (41) Novaki, L. P.; Moraes, E. O.; Gonçalves, A B.; de Lira, R. A.; Linhares, V. N.; de Oliveira, M. C. K.; Meireles, F. A.; Gonzalez, G.; El Seoud, O. A. Solvatochromic and Solubility Parameters of Solvents: Equivalence of the Scales and Application to Probe the Solubilization of Asphaltenes. Energy Fuels 2016, 30 (6), 4644–4652. DOI:10.1021/acs.energyfuels.6b00461 (42) El Seoud, O. A. Understanding Solvation. Pure Applied Chem. 2009, 81, 697-707. (43) Schleicher, J. C.; Scurto, A. M. Kinetics and solvent effects in the synthesis of ionic liquids: imidazolium. Green Chem. 2009, 11, 694–703. (44) Mutelet, F.; Ekulu, G.; Solimando, R.; Rogalski, M. Solubility Parameters of Crude Oils and Asphaltenes. Energy Fuels 2004, 18, 667−673. (45) Molina V., D.; Angulo, R.; Dueñez, F. Z.; Guzmań, A. Partial Least Squares (PLS) and Multiple Linear Correlations between Heithaus Stability Parameters (Po) and the Colloidal Instability Indices (CII) with the 1 H Nuclear Magnetic Resonance (NMR) Spectra of Colombian Crude Oils. Energy Fuels 2014, 28, 1802−1810. (46) Armarego, W. L. F.; Chai, C. Purification of Laboratory, Chemicals, 5th ed.; Elsevier: New York, 2003. (47) Standard Test Method for Determination of Aromatic Content and Polynuclear Aromatic Content of Diesel Fuels and Aviation Turbine Fuels by Supercritical Fluid

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30

Energy & Fuels

[23] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chromatography; ASTM Conshohocken. PA. 2009.

Standard

D5186-03(2009);

ASTM

International:

West

(48) IP Standard IP 143/84. In Standard Methods for Analysis and Testing of Petroleum and Related Products; Institute of Petroleum: London, 1989. (49) Altoé, R.; de Oliveira, M.C.K.; Lopes, H.E.; Teixeira, C.; Cirilo, L.C.M.; Lucas, E.F.; Gonzalez, G. Solution behavior of asphaltic residues and deasphalted oil prepared by extraction of heavy oil. Colloids Surf. A: Physicochem. Eng. Aspects 2014, 445. 59- 66. (50) Ovalles, C.; Rogel, E.; Moir, M.; Thomas, L.; Pradhan, A. Characterization of Heavy Crude Oils, Their Fractions, and Hydrovisbroken Products by the Asphaltene Solubility Fraction Method. Energy Fuels 2012, 26, 549−556. (51) Acevedo, S.; Méndez, B.; Rojas, A.; Layrisse, I. Asphaltenes and resins from the Orinoco basin. Fuel 1985, 64, 1741–1747. (52) Grin’ko, A. A.; Golovko, A. K. Fractionation of Resins and Asphaltenes and Investigation of Their Composition and Structure Using Heavy Oil from the Usa Field as an Example. Petrol. Chem. 2011, 51 (3) 192–202. (53) Naela, A.; Pogodim, S.; Cohen, S.; Agranat, I.; The conformational spaces of dinaphthyl ketones. dinaphthyl thioketones. and dinaphthyl diazomethanes: 1substituted naphthalenes versus 2-substituted naphthalenes. Struct. Chem. 2012, 23, 771-790. (54) Rupert. B.L.; Mulvihill. M.J.; Arnold. J. Atom-Transfer Radical Polymerization on Zinc Oxide Nanowires. Chem Mater. 2006, 18 (21), 5045. (55) Yukawa. Y.; Tsuno. Y.; Shimizu. N. The NMR chemical shifts of aromatic ring protons. II. The substituent chemical shift of the ring protons of 1-substituted naphthalenes. Bull. Chem. Soc. Jpn. 1971, 44, 3175. (56) Catalán, J.; Palomar, J.; Diaz, C.;, de Paz, J. L. G. J. Phys. Chem. A, 1997, 101, 51835189. (57) Pires, P. A. R.; Imran, M.; Loffredo, C.; Donate, P. M.; Previdi, D.; El Seoud, O. A. J. Phys. Org. Chem. 2013, 26 280–285. (58) Loffredo, C.; Pires, P. A. R.; Imran, M.; El Seoud, O. A. Carotene: A green, inexpensive, and convenient solvatochromic probe for the determination of solvent polarizability. Dyes Pigments 2013, 96, 16-24. (59) Reichardt, C.; Welton, T. Solvents and solvent effects in organic chemistry. 4th ed. Weinheim: Wiley-VCH; 2011, chapter 3.

ACS Paragon Plus Environment

Energy & Fuels

[24] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(60) Painter, P.; Veytsman, B.; Youtcheff, J. Guide to Asphaltene Solubility. Energy Fuels 2015, 29, 2951−2961. (61) Hosseini-Dastgerdi, Z.; Tabatabaei-Nejad, S.A.R.; Khodapanah. E.; Sahraei, E. A comprehensive study on mechanism of formation and techniques to diagnose asphaltene structure; molecular and aggregates: a review. Asia-Pac. J. Chem. Eng. 2015; 10: 1–14 (62) Mousavi, M.; Abdollahi, T.; Pahlavan, F.; Fini, E. H.The influence of asphaltene-resin molecular interactions on the colloidal stability of crude oil. Fuel 2016, 183, 262–271. (63) Yaseen, S.; Mansoori, G. A. Molecular dynamics studies of interaction between asphaltenes and solvents. J. Petroleum Sci. Eng. 2017, 156, 118–124. (64) Garreto, M.S.E.; Mansur, C.R.E.; Lucas, E.F. A model system to assess the phase behavior of asphaltenes in crude oil. Fuel 2013, 113, 318–322. (65) Nelyubov, D. V.; Semikhina, L. P.; Vazhenin, D. A.; Merkul’ev, I. A.; Influence of Resins and Asphaltenes on the Structural and Rheological Properties of Petroleum Disperse Systems. Petroleum Chemistry 2017, 57 (3), 203–208. (66) Spiecker, P. M.; Gawrys, K. L.; Trail, C. B.; Kilpatrick, P. K., Effects of petroleum resins on asphaltene aggregation and water-in-oil emulsion formation Colloids and Surfaces A: Physicochem. Eng. Aspects 2003, 220, 9-27 (67) Wallace, D.; Henry, D.; Pongar, K.; Zimmerman, D. Evaluation of some open column chromatographic methods for separation of bitumen components. Fuel 1987, 66, 44. (68) McLean, J. D.; Kilpatrick, P. K. Effects of Asphaltene Solvency on Stability of Water-inCrude-Oil Emulsions. J. Colloid Interface Sci. 1997,189, 242-253. (69) Boek, E. S.; Wilson, A. D.; Padding, J. T.; Headen, T. F.; Crawshaw, J. P.; Multi-scale Simulation and Experimental Studies of Asphaltene Aggregation and Deposition in Capillary Flow Energy. Fuels 2010, 24, 2361–2368. DOI:10.1021/ef9010059. (70) Headen, T. F.; Boek, E. S.; Jackson, G.; Totton, T. S.; Müller, E. A.; Simulation of Asphaltene Aggregation through Molecular Dynamics: Insights and Limitations. Energy Fuels 2017, 31, 1108−1125. (71) Murgich, J.; Abanero, J.A.; Strausz, O.P. Molecular recognition in aggregates formed by asphaltene and resin molecules from the Athabasca oil sand. Energy Fuels 1999, 13, 278–86. (72) Mullins, O. C.; Sabbah, H.; Eyssautier, J.; Pomerantz, A. E.; Barre, 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.

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30

Energy & Fuels

[25] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(73) Maqbool, T.; Balgoa, A. T.; Fogler, H. S.; Revisiting Asphaltene Precipitation from Crude Oils: A Case of Neglected Kinetic Effects. Energy Fuels 2009, 23, 3681–3686. (74) ASTM D6560 -00, IP 143/01, Standard Test Method for Determination of Asphaltenes (Heptane Insolubles) in CrudePetroleum and Petroleum Products. (75) Guzmán, R.; Ancheyta, J.; Trejo, F.; Rodríguez, S.; Methods for determining asphaltene stability in crude oils. Fuel 2017, 188, 530–543. (76) Mutelet, F.; Ekulu, G.; Rogalski, M. Characterization of crude oils by inverse gas chromatography. J. Chromatogr. A 2002, 969, 207−213.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1 203x96mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 2 289x202mm (150 x 150 DPI)

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3 289x202mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 4 289x202mm (150 x 150 DPI)

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5 289x202mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 30 of 30