3.1. Extractant Compositions Determined via Densimetry
The extractant used in this paper is a binary mixture of xylene (G) and n-heptane (P) in different G/P ratio. For standard samples of extractants with compositions listed in
Table 1, their densities
ρ were measured and plotted versus the volume fraction
vG.
ρ shows a linear relation with
vG in
Figure 1 and the coefficient of determination (
R2) is 0.99967 for the fitting equation
The densities
ρ of recovered extractants (P1, P2 and G1, G2) were measured by densimetry and substituted into Equation (1) to get
vG. Xylene and heptane of AR purity were used for the 1st extraction of asphaltenes. G1 and P1 were used for the 2nd extraction while G2 and P2 were used for the 3rd extraction. For heptane-rich fractions P1 and P2, their volume fraction of n-heptane
vP (
vP = 1 −
vG) is 0.907 ± 0.026 and 0.832 ± 0.011, respectively. The xylene-rich fractions G1 and G2 have higher purity, whose
vG is 0.992 ± 0.014 and 0.966 ± 0.014, respectively, as shown in
Figure 1. The negative pressure (vacuum) generated by the pump during vacuum distillation draws away a large amount of solvent and causes a solvent loss and air pollution issue. The average yield for distillate P1 was only 75%. Owing to volatility and good compatibility between xylene and n-heptane, a shorter time of vacuum distillation is preferred to improve the distillate yield and to avoid the xylene vapor dissolving into the heptane-rich fraction. Solvent separation conditions need further optimization, for example, under appropriate negative pressure and distillation time. To reduce the solvent loss in the recovery process, it may be necessary to connect a cold trap to the vacuum pump or use other separation methods instead of vacuum distillation [
22,
23]. For blank samples at
vG = 0, the
sB of calculated volume concentration is 1.2% and the LOD of
vG is 3.6% for the densimetry method.
The asphaltene yield was 13.2 ± 0.4%, 10.8 ± 0.4% and 8.3 ± 0.4%, respectively, for the 1st, 2nd, and 3rd extractions. The asphaltene yield decreased with the extraction number but the mechanism was not known exactly [
7]. To understand the scenarios of asphaltene precipitation more intuitively, the actual volumes of n-heptane and xylene in the solution for settling asphaltenes during the 1st, 2nd, and 3rd extractions were calculated according to their concentrations in
Figure 1 and listed in
Table 3, assuming 50 mL G and 500 mL P used in each extraction for ease of comparison. With increasing the extraction number, the heptane-rich fraction contains more and more xylenes, reducing the ability of n-heptane to sediment asphaltenes. Simultaneously, the concentration and total volume of xylene in the solution increase. Thus, more asphaltenes are dissolved while the sedimentation amount and final yield decrease. Even when the xylene-rich fraction is as low as
vG = 96.6%, the asphalt can be completely dissolved, while the heptane-rich fraction is still capable to precipitate asphaltene as low as
vP = 83.2%. Thus, the proposed extraction method [
7] has good robustness and great prospects for industrialization, since asphaltenes can be efficiently extracted with solvents much lower than AR purity. Goual et al. found that the onset of asphaltene flocculation occurred near a toluene/heptane volume ratio of 70:30 [
14]. The xylene/heptane volume ratio is 24:76 for the 3rd extraction in
Table 3, well below the critical ratio of asphaltene precipitation. How the precipitation amount varies with the G/P ratio is an interesting question and may help address asphaltene deposits in oil pipelines [
9,
24].
3.2. Properties of Asphaltenes Extracted at Room Temperature
Asphaltene is usually the most polar and cohesive component of asphalt, and forms the core of the asphalt colloid structure. In
Table 4, with increasing the extraction number (from A1 to A3), the mass fractions of C, H and S elements increase while the fractions for element N and others (mainly oxygen) decrease. The fraction of N decreases sequentially from 1.50% in A1 to 1.38% in A3 then to 0.73% in asphalt. Because element N increases the molecular polarity [
3], the polarity of asphaltenes from A1 to A3 may decrease with the extraction number and be larger than that of asphalt. C/H atomic ratio is defined as
q, which decreases sequentially from 0.99 in A1 to 0.78 in asphalt.
q reflects the degree of unsaturation of the molecule, usually in range of 0.5 to 1 for bitumen and its components [
1,
6,
7].
q is larger than 1 for fused-ring molecules, and approaches to infinity for chars and graphites [
25,
26]. As asphaltenes and asphalt have apparently different
q, their unsaturation discrepancy is much larger than the discrepancy between different asphaltenes [
6]. Since the polarity and unsaturation increases from asphalt to A3 then to A1, the π-π conjugation and inter-molecular interactions increase, causing higher cohesion and hardness for asphaltenes than asphalt. Thus, asphaltenes extracted at larger G/P ratio have smaller unsaturation, polarity and inter-molecular interactions.
Asphaltenes with larger unsaturation are more easily dissolved in solvents with larger G/P ratio, allowing asphaltenes with smaller unsaturation to preferentially precipitate out. Thus, the averaged molecular structure of the extracted asphaltene may change with the extraction number. Owing to complex compositions and chemical structures of asphalt and asphaltenes [
5], NMR is generally used to analyze their average molecular structure [
1,
27] rather than their content quantification [
2,
28]. As in
Figure 2, the protons are divided into A, α, β, and γ categories [
1,
6,
29], whose chemical shifts are, respectively, in ranges of 6.0–9.0 ppm, 2.0–4.0 ppm, 1.0–2.0 ppm and 0.5–1.0 ppm [
1,
6]. Their integrated areas are denoted as
hA,
hα,
hβ, and
hY, respectively. Protons directly attached to the single (S) and fused (F) rings are in the range of 6.0–7.4 ppm and 7.4–9.0 ppm, respectively [
30], whose integral areas are denoted by
hS and
hF in
Table 5. From A1 to A3 then to asphalt, both
hS and
hF decrease, which has a consistent trend with
q and indicates a decrease of aromatic ring fractions. The aromatic carbon ratio
fA [
1] shows a decrease trend while the hydrogen substitution rate around aromatic rings
σ [
6] shows an increase trend. Thus, with increasing the G/P ratio in extractant or with increasing the extraction number, the aromatic rings of extracted asphaltenes have more branched chains and occupy less proportion. From A1 to A3 then to asphalt,
hA decreases while
hβ and
hY increase. Therefore, the fraction of alkyl branches and (methyl) ends increases in extracted asphaltenes. In terms of averaged molecular structures, the extracted asphaltenes are significantly different from the base asphalt [
29]. In terms of physical properties, the asphaltenes obtained from the 1st, 2nd and 3rd extractions are all brittle, hard and shiny, without apparent adhesion to glass, skin or plastics at room temperature, which are apparently different from the base asphalt. The remaining asphalt after asphaltene extraction is much softer and stickier than the original (base) asphalt, which can be used as a new coating material.
In
Figure 3, main absorption peaks are marked with purple letters on FTIR spectra and the wave numbers for these peaks are listed in
Table 6 with peak assignments [
6,
31,
32]. A general rule is that the peak of asphaltenes locates at smaller wave numbers than that of base asphalt. For example, the peak g arises from the skeleton stretching (C=C) of aromatic rings. It locates at 1592 cm
−1 for A1 and locates at 1603 cm
−1 for base asphalt. The shift of peak location towards higher wave number from A1 to asphalt arises from the decrease of molecular aromaticity. The aromatic ring makes the distribution of electron cloud density between single-bond and double-bond more uniform. This conjugative effect reduces the electron cloud density of double bond, and reduces its force constant and vibration frequency, so that the peak g moves towards smaller wave number in
Figure 3 and
Table 6. Moreover, aromatic rings have stronger electron absorption than saturated alkanes. This inductive effect reduces the electron cloud density and reduces the force constant and vibration frequency of the adjacent bond, thus its absorption peak moves towards lower wave number in
Table 6. According to the previous EA and NMR results of
Table 4 and
Table 5, from asphalt to A3 then to A1, the values of
hS,
hF and
hA increase while the unsaturation and aromaticity of the molecular structure increase. Such increase may improve the inductive and conjugation effects and shift the peak location towards lower wave number, which is consistent with results in
Figure 3 and
Table 6.
The aromatic rings can only influence a small proportion of chemical bonds adjacent to the rings while asphalt and asphaltene may have dramatically different proportion. In
Table 6, the discrepancy between asphaltenes (A1, A2 and A3) is relatively smaller than that between asphaltenes and base asphalt. For example, the discrepancy of location for peak g is only 3 cm
−1 between A1 and A3 while it is 11 cm
−1 between A1 and asphalt, according to
Table 6. The reason may be that the proportion of chemical bonds influenced by aromatic rings increases more dramatically (larger discrepancy) from asphalt to asphaltene A3 than it increases from asphaltenes A3 to A1. From asphalt to asphaltene A3, the asphaltene percentage increases from below 20% to nearly 100%. Such a shift of peak location was also observed for asphaltenes and even resins near 1603 and 1030 cm
−1 in previous FTIR spectra [
6]. Moreover, as asphalt was a complex mixture and its four components had significantly different spectra [
6], the peak position and height may vary with the proportion of its components. Compared with asphaltenes, the saturates lacked peaks near 1030 cm
−1 (peak d) and 1600 cm
−1 (peak g) [
6], owing to lacking S and O elements and double-bond in their molecules. The peak around 1377 cm
−1 is used for the quantification of poly(styrene-
b-butadiene-
b-styrene) (AASHTO T302-15, JT/T 1329-2020) or other modifiers in asphalt [
31,
33,
34]. According to
Table 6, the position of this peak (e) moves towards a lower wave number from base asphalt to asphaltenes; thus, its specific position may vary with asphaltene proportion and should be checked carefully when calculating peak height or area integral for quantification.
The XRD spectra of asphaltenes mainly have 3 characteristic peaks of γ, 002 and 10 band, which are labelled as γ, m and a, respectively, in
Figure 4. The X-ray wavelength
λ is 1.5406 angstrom (Å). According to their diffraction angles (2
θ) listed in
Table 7 and the Bragg equation
d =
λ/(2 sin
θ) [
35,
36,
37], the characteristic lengths
d are calculated and listed in
Table 7. The peaks γ and m reflect the packing structure of saturated alkyl chains and aromatic rings, respectively, while the peak a is related to the size of the aromatic sheet [
35,
37,
38]. From A1 to A3 in
Table 7, the decrease of average distance
dY indicates a denser packing of alkyl chains, while the increase of layer spacing
dm implies a looser packing of aromatic rings. According to previous EA, IR and NMR results, the unsaturation, polarity and aromaticity of asphaltenes decrease while the number of substituted branches and fraction of methyl groups increase from A1 to A3. Thus, the packing of aromatic rings becomes looser owing to more substituted branches on aromatic rings and smaller inter-molecular interactions, while the packing of alkyl chains becomes denser because of the increased fraction of free chain (methyl) end. The decreased ring fraction and molecular unsaturation is also helpful to decrease the molecular rigidity and beneficial to the packing of alkyl chains. Therefore, the variation trend of molecular structure of asphaltenes reflected by XRD spectra is consistent with the previous experimental results.