CC BY
65ISSN 2072-5981 doi: 10.26907/mrsej
aänetic Resonance in Solids
Electronic Journal
Volume 21 Issue 2 Paper No 19201 1-10 pages
2019
doi: 10.26907/mrsej-19201
http: //mrsej. kpfu. ru http: //mrsej. ksu. ru
Established and published by Kazan University Endorsed by International Society of Magnetic Resonance (ISMAR) Registered by Russian Federation Committee on Press (#015140),
August 2, 1996 First Issue appeared on July 25, 1997
© Kazan Federal University (KFU)*
"Magnetic Resonance in Solids. Electronic Journal" (MRSey) is a
peer-reviewed, all electronic journal, publishing articles which meet the highest standards of scientific quality in the field of basic research of a magnetic resonance in solids and related phenomena.
Indexed and abstracted by Web of Science (ESCI, Clarivate Analytics, from 2015), Scopus (Elsevier, from 2012), RusIndexSC (eLibrary, from 2006), Google Scholar, DOAJ, ROAD, CyberLeninka (from 2006), SCImago Journal & Country Rank, etc.
Editor-in-Chief Boris Kochelaev (KFU, Kazan)
Honorary Editors
Jean Jeener (Universite Libre de Bruxelles, Brussels) Raymond Orbach (University of California, Riverside)
Executive Editor
Yurii Proshin (KFU, Kazan) mrsej@kpfu. ru
0 © I This work is licensed under a Creative ^-¿■¿»J Commons Attribution-Share Alike 4.0 International License.
This is an open access journal which means that all content is freely available without charge to the user or his/her institution. This is in accordance with the BOAI definition of open access.
Technical Editor Alexander Kutuzov (KFU, Kazan)
Editors
Vadim Atsarkin (Institute of Radio Engineering and Electronics, Moscow) Yurij Bunkov (CNRS, Grenoble) Mikhail Eremin (KFU, Kazan) David Fushman (University of Maryland, College Park) Hugo Keller (University of Zürich,
Zürich)
Yoshio Kitaoka (Osaka University,
Osaka)
Boris Malkin (KFU, Kazan) Alexander Shengelaya (Tbilisi State University, Tbilisi) Jörg Sichelschmidt (Max Planck Institute for Chemical Physics of Solids, Dresden) Haruhiko Suzuki (Kanazawa University, Kanazava) Murat Tagirov (KFU, Kazan) Dmitrii Tayurskii (KFU, Kazan) Valentine Zhikharev (KNRTU,
Kazan)
* In Kazan University the Electron Paramagnetic Resonance (EPR) was discovered by Zavoisky E.K. in 1944.
Short cite this: Magn. Reson. Solids 21, 19201 (2019)
doi: 10.26907/mrsej-19201
Size effect for asphaltene particles in the resin by NMR
D.S. Ivanov1 *, EE. Barskaya 2, V.D. Skirda 1
1 Kazan Federal University, Kremlevskaya 18, Kazan 420008, Russia
2 Arbuzov Institute of Organic and Physical Chemistry, Kazan 420029, Russia
*E-mail: f.ma.dima@mail.ru
(Received March 26, 2019; revised April 3, 2019; accepted April 5, 2019; published April 8, 2019)
The resin-asphaltene model systems with different sizes of asphaltenes are studied by nuclear magnetic resonance (NMR) using the solid-echo sequence. For systems with particle sizes of asphaltenes greater than 50 ^m, the form of relaxation attenuation consists of two components, one of which is characterized by a typical asphaltene relaxation time of about 20 ^s. In this case, the Gaussian form describes the relaxation attenuation of this component, which is typical for the NMR signal in disordered solids. At the same time, it is found that the share of the solid component in the NMR signal decreases with decreasing particle size of asphaltenes. The signal drops to zero when the size of asphaltenes particles reaches 10 microns. Based on the estimates of the rotational motion correlation time for an asphaltene particle in a dispersion medium (resin) with a known viscosity value, the influence of two factors is discussed: the difference in the values of macro- and microviscosity of the resin, as well as the possibility of partial dissolution of asphaltenes particles in the resin. The results of repeated measurements for a long time indicate the existence of a reverse process in the system, leading to the aggregation of asphaltene particles.
PACS: 75.10.Dg, 76.30.-v, 75.20.
Keywords: NMR, solid-echo, spin-spin relaxation, correlation time, asphaltene, resin, aggregation. 1. Introduction
A characteristic feature of natural bitumen and heavy oil is the high content of asphaltene-resin substances (ARS) [1]. The high content of ARS in general and asphaltenes, in particular, is the main cause of the problems that arise not only at the stage of production but also at the stages of transportation and processing of high-viscosity oil [2].
Extensive literature is devoted to the study of the structure of asphaltenes and their ability to aggregate [3-6]. At the same time, it was found that asphaltenes are characterized by an extremely rich set of chemical and supramolecular structures. Many, including rheological, characteristics of oil are associated with the chemical composition of asphaltenes [7-9]. According to the most well-established opinion [10], it is asphaltenes surrounded by resins which primarily affect the dynamic viscosity.
In most cases, asphaltenes previously separated from oil are investigated in low molecular weight solvents. At the same time, it is noted in [11] that the size of asphaltenes particles released from oil by precipitation, for example, in heptane, may not correspond to the initial size in crude oil and often leads [12-14] to their deposition even insolvent.
It can be considered unambiguously established [15] that in oil asphaltenes are dispersed only due to interaction with resin as one of the main components of hydrocarbon raw materials. Asphaltenes do not interact with the rest of the oil components and precipitate when mixed. At the same time, there is practically no work in which the issues related to the study of the properties of asphaltenes directly in interaction with resin molecules would be studied.
In work [16], it is shown that the nuclear magnetic resonance (NMR) technique also allows studies of the dynamics of hydrocarbons in proximity of asphaltene nano-aggregates and macro-aggregates in heavy crude oils as a function of the concentration of asphaltenes. It is useful for understanding the relations and correlations between NMR relaxation times T1 and T2, diffusion coefficients D, and viscosity of heavy crude oils. It is evident that prediction of crude oil stability and solubility of asphaltenes under specified conditions requires the fundamental knowledge of the character of intermolecular interactions between asphaltenes and other heavy oil components present in the local
environment. This knowledge can be gained only by means of in situ studies of local rheological properties of crude oils and molecular dynamics of asphaltenes under various conditions [17].
Thus, an important stage on the way to obtaining information about the state of asphaltenes in petroleum is the study of model systems for asphaltenes-resins, which was the purpose of this work. As the primary research method the method of NMR was chosen, it was proven itself in the best way in the investigation of complex molecular systems, which is certainly true of the oil [18, 19]. Further in the text, the words "molecule" will use as a fraction the components of the model system, which refers to resins or asphaltenes.
2. Experimental methods
Petroleum asphaltenes and resins were obtained from petroleum field of the Republic of Tatarstan in the laboratory of the A.E. Arbuzov Institute of organic and physical chemistry, Kazan scientific Centre Russian Academy of Sciences. Asphaltenes are obtained from oil by precipitation from petroleum solution in petroleum ether. According to the data of optical measurements on the MTI EQ-MS microscope, the initial size of the asphaltenes particles obtained in this way is estimated to be approximately 1 mm. To prepare a series of model systems of asphaltene-resin with different sizes of asphaltene particles, the initial asphaltenes were ground using a disk eraser LDI, manual agate mortar and further sieved through a set of sieves. Because of this operation, powders of asphaltenes with the size (diameter) of particles of 1000, 51, 38, 22 and 10 microns were prepared. The size of the obtained particles was controlled by measurements on the MTI EQ-MS microscope.
By mechanical mixing of resin and asphaltene powders with a given mass concentration (5 %), five samples of asphaltene-resin with different particle sizes of asphaltenes were prepared. The mass content of asphaltenes in the samples was calculated by the ratio:
m ,
a=-^—100%,
masph + mresin
where masph is mass content of asphaltenes, mresin is mass of resin in the model sample. The measurements were carried out on the NMR analyzer "Proton-20M" with a resonance frequency of = 20 MHz protons with a characteristic time of "paralysis" of the receiving path rv = 10 ps. In order to obtain the most complete information about the NMR signal at times commensurate with the paralysis time t, the solid-echo pulse sequence [20, 21] shown in figure 1 was used.
Figure 1. Solid-echo pulse sequence [19].
All experiments were carried out with a duration of n/2 radio frequency (RF) pulse of 2.7 ps, and a value of t = 13 ps at temperature T = 310 K. The viscosity of the resin measured on the viscometer Fungilab Expert rotary type at a temperature of T = 310 K.
3. Experimental results and discussion
Figure 2 shows relaxation attenuations normalized to the initial amplitude corresponding to the minimum and maximum particle size (dasph = 10 and 1000 pm) in a series of prepared asphaltene samples. The solid line in figure 2 shows the approximation of the experimental data by the equation with the Gaussian form of the line:
A ( t ) = A ( 0)x exp
' t ^
T
y1 2 s 7
(1)
where A(t) is a function describing the evolution of the NMR signal in time, A(0) is the initial amplitude of the relaxation attenuation, T2s is the spin-spin relaxation time corresponding to the Gaussian component.
The equation (1) is usually true for the case of so-called solids with disordered structure in the arrangement of atoms. According to the theory [17] to register the so-called "solid-state" NMR signal satisfying the equation (1), the following condition is sufficient:
t >> T
(2)
where T2s is the recorded spin-spin relaxation time and Tc is the correlation time of local magnetic fields due to the interaction between neighboring nuclei. In most cases, the specified correlation time is controlled by rotational and translational modes of molecular motion.
Thus, relaxation attenuations for asphaltenes particles are described quite satisfactorily by the equation of type (1), as it was shown in earlier works [21, 22], and are characterized by a fairly short time of spin-spin relaxation. As it can be seen from figure 2, relaxation attenuations for both sizes of asphaltenes particles coincide up to the experimental error, which indicates, in particular, the absence of any mechanical destruction effects because of the grinding procedure of the initial asphaltenes particles applied by us.
Figure 3 shows the free induction declines normalized to the initial amplitude on a semi-logarithmic scale, registered by the solid-echo sequence for the resin and the model system consisting
o
1
0.1-
0.01-
Asphaltene (size cfasph=1000 pm) Asphaltene (size cfasph=10 pm) Approximation line
1E-3
0
20
60
80
40 t, |S
Figure 2. Relaxation attenuations for asphaltenes obtained by the solid-echo pulse technique. Particle size of asphaltenes dasph = 10 pm is the black symbols, dasph = 1000 pm is the green symbols. The solid line shows the result of approximation of the experimental results by the equation (1) at the value T2s = 20 ps. Measurement temperature 310 K.
of a mixture of resin and asphaltenes (weight 5 %) with the initial particle size of ~ 1 mm. The same figure 2 shows for comparison the relaxation decay for Mordovian-Karmal bitumen with the same content (weight 5 %) asphaltenes.
It can be seen that the decay of free induction for the resin sample (red symbols) in semi-logarithmic coordinates is approximated by a straight line, which corresponds to the equation usually used to describe the Lorentz form of relaxation attenuation in liquids:
A, (t) = A (0)x exp
f
\
T
J o
(3)
2l J
where Al(t) is a function describing the signal with the Lorentz form, A(0) is the initial amplitude corresponding to the time t = 0, 72l is the transverse relaxation time.
At the same time, for the model sample containing 5 % asphaltenes and Mordovian-Karmal bitumen, a more complex type of relaxation attenuation is observed, which is satisfactorily approximated by the following equation (4):
AiH
A ( 0 )
= p, x exp
T
y1 2s J
+ pl x exp
t
A
T
J o
(4)
2l J
where A(t) is a function describing relaxation attenuation, A(t = 0) is the initial amplitude, ps and T2s is the population and characteristic relaxation time components with Gaussian attenuation form, pl and T2l are the population relaxation time of the Lorentz (liquid phase) part of the NMR signal.
The dotted lines in figure 3 show the behavior of the Lorentzian components in the approximation of the relaxation attenuation of the equation (4). The resulting approximation of the value of ps for the model sample of the resin and asphaltenes Mordovian-Karmal bitumen amounted to ps = 4.8 and 4.9 % respectively. This is slightly lower than the expected value of 5 % and it is due to the lower proton density of asphaltenes compared to the proton density of the resin. Thus, in [23] according to elemental analysis and 1H NMR spectroscopy, it is shown that the hydrogen content in asphaltenes can vary from 6.5 to 7.8 %. Whereas for resins the hydrogen content depending on the structural-group composition reaches values from 6.9 to 10 % [24, 25].
O
5C
1.00-
0.96
0.92-
0.88-
0
20
60
40 t, MS
Figure 3. Normalized to the amplitude at the initial time (t = 0) free induction decay, obtained by using the pulse sequence solid-echo, pulse sequence parameters: RF pulse duration is the 2.7 ^s, t = 13 ^s, the number of accumulations is the 1000. Dotted lines represent contributions to the decay of free induction from the liquid-phase components of the system under study. Red symbols are resin; black symbols are resin with 5 % asphaltene content (particle size ~ 1 mm), blue symbols are Mordovian-Karmal bitumen with 4.8 % asphaltene content. Measurements were made at room temperature 310 K.
Thus, by the example of the study of the model system of asphaltene-resin with particle sizes of 1 mm, it can be concluded that the recorded ps values with accuracy to the error of determining the proton density of the oil components are well correlated with the mass content of asphaltenes in the samples. This, in particular, confirms the applicability of the proposed in [21] method for determining the content of asphaltenes in oil samples based on NMR data.
This conclusion, in fact, is based on the fact that the NMR signal for asphaltenes has distinct characteristics. It is small (of the order of ten ps), the transverse relaxation time and a Gaussian, or close to it, the form of the relaxation attenuation.
For the colloid asphaltene-resin system, the correlation time tc in condition (2) should be understood as the correlation time of the rotational motion of the asphaltene particle as a whole. It is obvious that for the considered system, resins-asphaltenes with particle sizes of the order of about 1 mm, the fulfillment of the condition (2) is not in doubt. However, it is a priori unclear to what extent this condition will remain true in a situation where the particle sizes of asphaltenes in the resin will be significantly reduced.
As it turned out, for all the studied samples, except for the sample in which the particle size of asphaltenes was 10 pm, the relaxation attenuation is satisfactorily described by the equation (4), which assumes the presence of both the Lorentz and Gaussian components in the relaxation attenuation. The dotted lines in figure 4 show the behavior of the Lorentz components. An attempt to approximate the data for the asphaltene-resin mixture with a particle size of 10 pm is also formally possible by the equation (4). However, a good approximation result is achieved by assuming a (4) zero value for ps. Figure 5 shows the dependence of the values of the ps is the fraction of solid component in the relaxation attenuation of the solid-echo from the diameter of the particles dispersed in the resin particles of asphaltenes.
As it can be seen from figure 5, the proportion of the Gaussian component does not differ much from the expected about 5 % only for the asphaltene-resin system with a particle size of 51 pm. At the same time, for samples with smaller values of asphaltene particle sizes, the underestimated value of the ps fraction is clearly visible. Moreover, for a sample with a particle size of 10 pm, there are no signs of asphaltenes in the system at all in the experimentally recorded relaxation attenuation. Thus, for the obtained model samples, an effect is found that indicates the dependence of the fraction of the registered Gaussian component with a small relaxation time T2s on the particle size of asphaltenes
1
0.96-
CD
0.92-
Particle size of asphaltenes ~ ^
0.88-
51 |im 38 |m 22 |m 10 |m
0
20
60
40 t, |S
Figure 4. Normalized by amplitude at the initial time (t = 0) free induction decreases for samples with 5 % asphaltenes content with different particle sizes in the resin obtained by the solid-echo pulse technique. Dotted lines represent the contributions to the decay of free induction of liquid-phase components of oil (Lorentz components) in the approximation of the data by the equation (4). Measurements were made at room temperature 310 K.
dispersed in the resin. The obtained result is quite unexpected and, at first, glance, calls into question the previously published works [26, 27] on the study of oil samples by NMR method. The same works include patents [28-30], in which it is proposed to determine the share of asphaltenes in oil samples by registering the solid component in the relaxation attenuation. However, during repeated experiments, we observed that the form of relaxation attenuations tends to change over time. Therefore, in figure 6 shows the results of the study of the sample asphaltene-resin with a particle size of 10 microns, obtained at different times after sample preparation, are shown.
It can be seen that over time the system begins to form a Gaussian (solid) component. The approximation of the experimental data by the equation (4) showed that after 10 days the initial zero value ofps in percentage terms became equal to 1.6 %, after 55 days it reached the level of 4.2 %. At the same time, the relaxation time ?2s characteristic for this component is 20-25 ps, which, along with the Gaussian form of relaxation attenuation for this component, corresponds to the conditions of the solid (2).
In general, the results shown in figure 6 indicate the existence of a certain process, because of which the asphaltene-resin system according to the characteristics of the solid component tends to an equilibrium state:
0
60
20 40
^asphaltene- Mm
Figure 5. The dependence of the proportion of solid components in the relaxation attenuation of the solid-echo from the diameter of the particles dispersed in the resin asphaltenes. The dotted line shows the level of ps (%) obtained by approximating the experimental data for the resin-asphaltene mixture with a particle size of 1 mm.
1.00-
ï>j
O
S
0.95-
0.90-
*V -v- ■v.1 .
0 days 10 days 40 days 55 days
0
20
"60"
40
t, |js
Figure 6. Relaxation attenuations normalized by the amplitude at the initial time (t = 0) for the model sample resin-asphaltenes (weight 5 %) with a particle size of 10 ^m obtained at different times after sample preparation: 0 days are the black symbols, 10 days are the pink symbols, 40 days are the red symbols, 55 days are the blue symbols. Dotted lines represent the contributions to the decay of free induction of liquid-phase components of oil (Lorentz components) in the approximation of the data by the equation (4). Measurement temperature and sample storage temperature 310 K.
lim ps (t)- ps (t)
lim ps (t)
= Ps(t), (5)
where as ps(t) at t ^ <x the value ps = 0.048 obtained by us for a mixture of asphaltene-resin with a particle size of 1 mm is taken, Ps is a function designation introduced for brevity. Figure 7 shows the time dependence of the solid component.
In fact, the dependence presented in figure 7 shows the kinetics of the process associated with the achievement of a certain equilibrium state in the asphaltene-resin system according to the characteristics of the solid component, namely, its value. In our opinion, this result correlates with the statements of many researchers [3-5] about the tendency of asphaltenes to aggregation. In other words, the result can be interpreted as a gradual increase in the size of asphaltene units. The growth of asphaltene aggregates should be accompanied by an increase in the correlation time in (2) and at the same time a decrease in the relaxation time for asphaltene protons. As a result, when a certain size of the aggregate is reached, the condition (2) begins to be fulfilled, which makes it possible to fix the characteristic features of the solid component in the NMR signal: a small relaxation time T2s and a Gaussian form of relaxation attenuation of type (1). The straight line shown in figure 7 is calculated according to the equation
limn ps (t)- ps (t) f t ^
■ = exp
t
agr y
(6)
l im ps (t) where tagr = 28 days.
The parameter tagr is a time characterizing the aggregation process. Thus, the presented data clearly indicate that they reflect a certain process of relaxation of the system in an equilibrium state with a characteristic time of about 30 days.
Thus, on the basis of the above reasoning, it should be assumed that the absence of the NMR signal characteristic of asphaltenes in a freshly prepared asphaltene-resin system with particle sizes of 10 pm is due to the short correlation time of rotational motion for a particle with such sizes and, consequently, to the failure of the condition (2). We will try to confirm this by calculations.
The value of the rotational correlation time xc for solid particles dispersed in a dispersion medium is associated [18] with the particle size a and viscosity q. The most common equation for determining the correlation time is [21, 31] following:
„3
rc = ^L, (7)
c 3kT
t, day
Figure 7. The dependence of the proportion of solid-state component ps, recorded at a time interval of 55 days. The solid line shows the result of approximation of experimental data.
where r is the radius of the particle; n is the viscosity of the dispersion medium; k is the Boltzmann constant; T is the absolute temperature.
Setting for the model system asphaltene-resin with the particle size (diameter) of 10 pm r = 5 pm, T = 300 K and the value of the dynamic viscosity of the resin n = 10 MPa-s (according to measurements on the rotational viscometer Fungilab Expert) we obtain:
4nnr3 4x3.14x10x10-3 x(5x10-6)3
Tc =-=-^-— «103 s.
c 3kT 3 x1.38x10-23 x300
As it can be seen from the calculation results, the obtained estimate of the t value for the taken values of the parameters indicates just a very good fulfillment of the condition (2), but and not Vice versa. Therefore, for given the particle sizes of asphaltenes and the viscosity of the resin, a Gaussian component with ps ~ 5% and a short relaxation time T2s had to be registered at the initial site of relaxation attenuation. However, this is contrary to the experimental data.
At first glance, the only parameter in formula (6) the value of which can be questioned is the dynamic viscosity n. The value n = 10 MPa-s used in the calculation was measured by the Fungilab Expert viscometer and, in fact, characterizes the macro viscosity of the system. At the same time, it is shown in [32, 33] that the values of macro- and microviscosity in complex molecular systems can be very different, although it is the latter that makes sense to use in calculations by the ratio (6). Direct estimation of the value of microviscosity from any independent measurements is, in this case, a certain difficulty. Therefore, we will try to calculate it from the following considerations. From the analysis of the data shown in figure 5, it follows that at a particle size of about 30 pm, the recorded ps fraction will reach the level of the mean value between the zero level and the maximum 4.8 % measured for the resin-asphaltene system with particle sizes of 1 mm. Then it can be assumed that it is at this particle size for the correlation time Tc that the condition is met: Tc ~ T2s ~ 2-10-5 s. Then, substituting in (6) the value of the radius of the particle r = 15 pm, we can calculate the value of the microviscosity n at the requirement: Tc = 2-10-5 s. However, the calculation result showed, that to fulfill this requirement, the viscosity value in (6) should be about 10-12 Pa-s, which is eight (!) decimal orders of magnitude smaller than the viscosity of water (nwater = 1.004 MPa-s) at normal conditions. Thus, an option based only on the assumed difference between macro- and microviscosity in formula (6) is at least insufficient to explain the experimental results. Formally, the radius of the particle r remains the only parameter on which the required result of calculating the correlation time Tc of the asphaltene particle can still depend. Indeed, the requirement Tc = 2-10-5 s is satisfied if the equation (6) substitute the value r = 14 nm. The result is, in fact, the lower estimate for the particle size of asphalt, since it was calculated under the condition that the measured value of macro viscosity (10 MPa-s) of the resin is preserved for parameter n. It is clear that the value of the parameter r can be slightly increased if in (6) substitute a reasonable value of microviscosity n. However, on the one hand, the value of microviscosity for the system under study is unknown to us, and on the other hand, the use of this value, it will not change any fundamentally overall result. Formally, it consists in the need for the assumption for the parameter r in (6) significantly lower values compared to the size of the particles of asphaltenes introduced into the resin at the stage of preparation of model samples.
4. Conclusion
The need to mutually agree of all the results obtained to recognize that the introduction of asphaltene particles in the resin can be a process of crushing them into smaller particles or, possibly, partial dissolution. The assumption of partial dissolution of asphaltenes in the resin is justified by some similarity of the chemical structures of resin molecules and asphaltenes, as it was already noted in [1, 25, 26]. At the same time, it is impossible not to take into account the possibility of violation of some integrity of asphaltenes particles at the stage of their mechanical grinding.
Along with the above, an important result of the study is to establishment of the process of aggregation of asphaltenes over time. This is evidenced by the data shown in figures 6 and 7. On the one hand, the evidence of asphaltene aggregation contradicts the above hypothesis of crushing or partial dissolution of asphaltene particles in the resin, since these processes have the opposite direction. On the other hand, according to the data obtained, the aggregation process is rather slow (typical time is about 30 days). Therefore, its background may manifest the mentioned signs of partial dissolution particles of asphaltene.
In general, the results of the study of model asphaltene-resin systems indicate the existence of nontrivial mechanisms of interaction between the molecules of these components of oil.
References
1. Unger F.G., Andreeva L.N., Krasnogorskaya N.N. Fundamental Aspects of Oil Chemistry, Nauka, Novosibirsk (1995) (in Russian)
2. Akbarzadeh K., Dhillon A., Svrcek W.Y., Yarranton H.W. Energy & Fuels 18, 1434 (2004)
3. Mozaffari S., Tchoukov P., Atias J., Czarnecki J., Nazemifard N. Energy & Fuels 29, 5595 (2015)
4. Franco C.A., Nassar N.N., Montoya T., Ruíz M.A., Cortés F.B. Energy & Fuels 29, 1610 (2015)
5. Headen T.F., Boek E.S., Jackson G., Totton T.S., Müllér E.A. Energy & Fuels 31, 1108 (2017)
6. Painter P., Veytsman B., Youtcheff J. Energy & Fuels 29, 2120 (2015)
7. Spiecker P.M., Gawrys K.L., Trail C.B., Kilpatrick P.K. Colloids and Surfaces A: Physicochemical and Engineering Aspects 220, 9 (2003)
8. Chilingarian G.V., Yen T.F. Bitumens Asphalts and Tar Sands, Elsevier Scientific Pub. Co., Netherlands (1978)
9. Dickie J.P., Yen T.F. Anal. Chem. 39, 1847 (1967)
10. MacDonald B.A., Miadonye A. Journal of Chemical & Engineering Data 62, 924 (2017)
11. ASTM D6560-00. Standard Test Method for Determination of Asphaltenes (Heptane Insolubles) in Crude Petroleum and Petroleum Products, Annual Book of ASTM Standards, West Conshohocken, PA ASTM International (2005)
12. Sheu E.Y., Liang K.S., Sinha S.K., Overfield R.E. Journal of Colloid and Interface Science 153, 399 (1992)
13. Soorghali F., Zolghadr A., Ayatollahi S. Energy & Fuels 28, 2415 (2014)
14. Akbarzadeh K., Hammami A., Kharrat A., Zhang D., Allenson S., Creek J., Mullins O.C. Oilfield Review 19, 22 (2007)
15. Speight J.G., Moschopedis S.E. Am. Chem. Soc. 24, 1 (1979)
16. Korb J.P. Prog. Nucl. Magn. Reson. Spectrosc. 104, 12 (2018)
17. Martyanov O.N., Larichev Y.V., Morozov E.V., Trukhan S.N., Kazarian S.G. Russ. Chem. Rev. 86, 999 (2017)
18. Abragam A. The Principles of Nuclear Magnetism, Oxford University Press, Oxford (1961)
19. Prunelet A., Fleury M., Cohen-Addad J.P. Comptes Rendus Chimie 7, 283 (2004)
20. Espinat D., Gaulier F., Norrant F., Barbier J., Guichard B., Rivallan M., Levitz P. Energy & Fuels 31, 7382 (2017)
21. Shkalikov N.V., Skirda V.D. Uchenye Zapiski Kazanskogo Universiteta - Seriya Fiziko-Matematicheskie Nauki 151, 41 (2009)
22. Gentile L., Filippelli L., Rossi C.O., Baldino N., Ranieri G.A. Molecular Crystals and Liquid Crystals 558, 54 (2012)
23. Kalabin G.A., Kanitskaya L.V., Kushnarev D.F. Quantitative NMR Spectroscopy of Natural Organic Raw Materials and Products of Its Processing, Khimiya, Moscow (2000) (in Russian)
24. White D.R., Constantinou C., Martin R.J. The British Journal of Radiology 59, 787 (1986)
25. Shkalikov N.B. Investigation of Heavy Oils and Their Components by NMR. (PhD thesis), Kazan Federal University, Kazan (2010) (in Russian)
26. Evdokimov I.N., Eliseev N.Y., Akhmetov B.R. Fuel 82, 817 (2003)
27. Shkalikov N. V., Skirda V.D., Arkhipov R.V. Magn. Reson. Solids 8, 38 (2006)
28. Patent 2333476 RF, G01N 24/08 Method for Determination of Paraffin and Asphaltene Content in Oil, Nikolin I.V., Safonov S.S., Skirda V.D., Shkalikov N.V., the applicant Schlumberger technology B. V. (2006) (in Russian)
29. Patent 2383884 RF, G01N 24/08 Method for Determining the Content of Liquid-Phase and SolidState Components in a Mixture of Hydrocarbons, Nikolin I.V., Shkalikov N.V., Skirda V.D., the applicant Schlumberger technology B.V. (2009) (in Russian)
30. Patent 2423686 RF, G01N 24/08 Method for Determining the Molecular Mass Distribution of Paraffins in a Hydrocarbon Mixture Using the Nuclear Magnetic Resonance Method, Shkalikov N.V., Skirda V.D., the applicant and the priority Kazan Federal University, (2009) (in Russian)
31. Farrar T.C. Introduction to Pulse NMR Spectroscopy, Farragut Press, Kingston Pike Farragut (1989)
32. Azami K., Yokono T., Sanada Y., Uemura S. Carbon 27, 177 (1989)
33. Lichtenstein G.I. Method of Spin Labels in Molecular Biology, Nauka, Moscow (1976) (in Russian)
