Нефтехимия, 2020, T. 60, № 4, стр. 520-528
Спектроскопический анализ асфальтенов нефти на молекулярном уровне
S. Ok 1, *, N. Rajasekaran 1, M. A. Sabti 1, G. A. Joseph 1
1 Petroleum Research Center, Kuwait Institute for Scientific Research
13109 Safat, P.O. box 24885, Kuwait
* E-mail: sok@uos.de
Поступила в редакцию 20.05.2019
После доработки 01.10.2019
Принята к публикации 11.03.2020
Аннотация
В настоящей работе для решения сложной научной задачи, заключающейся в определении структурных параметров асфальтенов, содержащихся в легких, среднетяжелых и тяжелых нефтях, были использованы различные спектроскопические методики. Методом ИК-спектроскопии с преобразованием Фурье помимо идентификации различных функциональных групп асфальтенов выявлены показатели, позволяющие их классифицировать. На основании анализа спектров комбинационного рассеяния асфальтенов установлено, что асфальтены легких нефтей имеют больший средний размер молекул и более сложное строение молекул по сравнению с асфальтенами тяжелых и среднетяжелых нефтей. Данный факт объясняет склонность асфальтенов к агрегации и осаждению при концентрации не более 1%, что как раз характерно для легких нефтей.
Осаждение асфальтенов является одной из самых главных проблем нефтяной промышленности [1, 2]. Как на макро-, так и на микроуровне асфальтены могут препятствовать течению нефти и изменять проницаемость нефтяных коллекторов. В процессах переработки нефти асфальтены снижают выход продуктов в результате дезактивации катализаторов. Проблемы, связанные с осаждением асфальтенов, наиболее часто проявляются при работе с легкими нефтями, содержащими незначительное количество (<1 мас. %) асфальтенов, стабильными только при концентрации асфальтенов выше определенного значения [3].
Асфальтены представляют собой полярные молекулы, содержащие такие электроотрицательные атомы, как N или S, и металлы подгруппы Fe и Ni. Они растворимы в ароматических углеводородах (УВ), например, в толуоле, и нерастворимы в н-пентане и н-гексане [1, 4, 5]. Однако знание точного химического строения недостаточно для объяснения свойств асфальтенов [6]. Согласно модели “молекула–архипелаг” асфальтены представляют собой небольшие полиароматические центры, соединенные между собой алифатическими или нафтеновыми фрагментами. В моде-ли Йена−Муллинса асфальтены состоят из полиароматического центра и боковых парафиновых цепей [5, 6]. Перекрывание делокализаванного π-электроннного облака центральной части молекулы и σ-электронов боковых цепей приводит к π-σ-сопряжению и образованию агрегатов асфальтенов [7]. Также свой вклад в процесс агрегирования вносят водородные связи, кислотно-основное и индукционное взаимодействие, силы Ван-дер-Ваальса, а также электростатическое притяжение молекул [3, 7].
В предыдущих работах для изучения асфальтенов легких и тяжелых нефтей были использованы методы инфракрасной спектроскопии с преобразованием Фурье нарушенного полного внутреннего отражения (Фурье-ИКС НПВО) и инфракрасной спектроскопии с преобразованием Фурье диффузного отражения (Фурье-ИКС ДО) [8, 9]. Интегральная площадь под графиком рамановского спектра углеродного пика позволяет определить средний размер сконденсированного полиароматического центра, в котором все полосы принадлежат плоской (микрокристаллической) графитоподобной структуре ароматического центра [10, 12].
В работе рассмотрены следующие важные вопросы: 1) в какой степени асфальтены легких нефтей отличаются от асфальтенов тяжелых и среднетяжелых нефтей по средней молекулярной массе и соотношению алифатических/ароматических фрагментов; 2) какая из моделей наилучшим образом подходит для объяснения свойств асфальтенов; 3) почему низкие концентрации асфальтенов (менее 1%) в значительной степени осложняют процессы транспортировки и переработки нефти. Задачей данного исследования является нахождение ответов на эти вопросы при изучении легкой, средней и тяжелой нефтей Кувейта и тяжелой нефти сорта Maya (Мексика).
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
Реактивы и образцы асфальтенов. Толуол и н-гептан (ЧДА) приобретены в компании Fisher Scientific и использованы без дополнительной очистки. Образцы асфальтенов получены экстракцией из пяти различных нефтей, проведенной по стандартной методике IP 143/01 (табл. 1) [6, 13]. Полученные образцы были промыты чистым гептаном в аппарате Сокслета для удаления остаточных углеводородов и смол. Подробная информация о процессе экстракции, проведенном в данном исследовании, приведена в литературе [14].
Таблица 1.
Свойства | Единицы измерений | Методика испытания | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|---|---|
Плотность при 15°C | г/см3 | ASTM D5002 | 0.88518 | 0.98232 | 0.98842 | 0.8399 | 0.83459 |
Относительная плотность | ASTM D5002 | 0.88605 | 0.98329 | 0.98939 | 0.8408 | 0.83542 | |
API | град. | ASTM D5002 | 28.20 | 12.40 | 11.52 | 36.80 | 37.88 |
Содержание асфальтенов | мас. % | IP143 | 5.92 | 18.62 | 23.10 | 0.38 | 0.35 |
Инфракрасная спектроскопия. Асфальтены были проанализированы на ИК-спектрометре с Фурье преобразователем Shimadzu IRTracer-100, оснащенным приставкой нарушенного полного внутреннего отражения с призмой на основе кристалла германия. Твердые асфальтены были предварительно измельчены в ступке для обеспечения необходимого контакта с кристаллом. Для получения спектров Фурье-ИК НПВО (4000–500 см–1) с разрешением 1 см–1 были обработаны 64 результата сканирования образцов. Базовая линия спектров скорректирована с помощью программы Shimadzu LabSolutions IR Spectrum.
Для получения Фурье-ИК-спектров диффузного отражения измельченных асфальтенов в диапазоне от 4000 до 400 см–1 с разрешением 4 см–1 были обработаны 64 результата сканирования, полученные с помощью спектрометра Bruker Vertex 80, оснащенного DTG детектором. Смесь образца асфальтенов и KBr была приготовлена для анализа согласно стандартной методике и помещена в прободержатель. Итоговые спектры диффузного отражения получены после отсечения фонового сигнала с помощью програмного обеспечения OPUS.
Рамановская спектроскопия. Образцы асфальтенов были измельчены до мелкого порошка. Исследования были проведены на рамановском спектрометре Bruker Senterra, оснащенном тринокулярной насадкой с видеокамерой, полученные изображения были перенаправлены в программу OPUS 7.5 (Senterra). Рамановский спектрометр с 50-кратным объективом для фокусировки лазера с длиной волны 488 нм и аккумулирования рассеянного света был использован в режиме обратного рассеяния. Для фиксации сигнала выбран высокочувствительный CCD-детектор с разрешением 1024 × 256 пикселей, термоэлектрически охлажденный до −70°С. Все спектры получены в диапазоне 90−4400 см–1. После корректировки базовой линии для получения оптимального числа полос, необходимого для идентификации пика, были использованы функции Гаусса и Лоренца, ширина пика была измерена на середине высоты. Были зафиксированы положение пика, ширина, высота и интегральная площадь.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Анализ спектров Фурье-ИК НПВО. Спектры асфальтенов имеют как сходства, так и определенные различия (см. рис. 1 и табл. 2) [9, 15–17].
Таблица 2.
Положение пика, см–1 | Интерпретация | Номер образца | Ссылка |
---|---|---|---|
695.3; 692.9; 684.2 | Внеплоскостное деформационное колебание ароматических связей С–Н | 3, 4 и 5 | [8] |
721.8; 726.6; 721.3; 726.1; 711.2 | Внеплоскостное деформационное колебание связей С–Н в ароматических соединениях и деформационное (маятниковое) колебание связей С–Н и связей CH2 (линейные алканы с длиной цепи – 4 или более атомов С) | 1, 2, 3, 4 и 5 | [13] |
749.8; 754.6; 744.5 | Внеплоскостное деформационное колебание ароматических С–Н-связей | 2, 4 и 5 | [13] |
800.4; 800.9; 810.5; 812.5; 810.5 | Внеплоскостное деформационное колебание ароматических С–Н-связей | 1, 2, 3, 4 и 5 | [16] |
888.2; 873.2; 875.6; 863.1; 866.0 | Внеплоскостное деформационное колебание ароматических С–Н-связей | 1, 2, 3, 4 и 5 | [16] |
1028.0; 1028.5; 1026.1; 1025.6; 1034.2 | Валентные колебания S=O в сульфоксидах | 1, 2, 3, 4 и 5 | [13] |
1375.6; 1374.7; 1365.0; 1376.1; 1376.6 | Деформационные колебания метиленовых групп | 1, 2, 3, 4 и 5 | [16] |
1456.6; 1451.8; 1458.1; 1454.7; 1458.1 | Деформационные колебания метиленовых групп | 1, 2, 3, 4 и 5 | [16] |
1586.3; 1591.1; 1569.0; 1589.2; 1588.3 | Полициклическая ароматика | 1, 2, 3, 4 и 5 | [18] |
1645.6 | – | 5 | |
1702.0; 1697.2; 1695.3; 1699.6; 1699.6 | Валентные колебания C=O групп | 1, 2, 3, 4 и 5 | [16] |
2851.5; 2853.5; 2846.7; 2851.5; 2852.0 | Симметричные валентные колебания метиленовых групп | 1, 2, 3, 4 и 5 | [16, 17] |
2920.0; 2922.4; 2920.5; 2922.4; 2923.8 | Асимметричные валентные колебания метиленовых групп | 1, 2, 3, 4 и 5 | [16] |
Симметричное/ассиметричное валентное и деформационное колебания метиленовых групп проявляются в виде наиболее интенсивных полос в районе 2920/2850 см–1 и 1456/1373 см–1 соответственно. Сигналы, располагающиеся между 865 и 680 см–1, соответствуют деформационным колебаниям ароматических С−Н-групп [9, 17]. Также в спектрах можно наблюдать другие сигналы, общие для всех образцов: полосы валентных колебаний сульфоксидных S=O-связей, наблюдаемые при 1028 см–1 [16, 19], полосы около 1695 см–1, соответствующие валентным колебаниям карбонильных C=O-групп, и полосы полиароматических соединений между 1570 и 1590 см–1. Спектры образцов асфальтенов с номерами 2, 4 и 5, в отличие от прочих, имеют полосу около 750 см–1, принадлежащую внеплоскостным колебаниям ароматических С−Н-связей [16]. В то же время полоса при волновом числе 690 см–1, возникающая также вследствие внеплоскостных колебаний ароматических С−Н-связей [9], наблюдается только в спектрах Фурье-ИК НПВО образцов 3, 4 и 5. Полоса при 1645.6 см–1, присутствующая в спектре образца номер 5, не идентифицирована. Полосы в образцах 4 и 5 имеют схожую интенсивность, а образец 2 характеризуется наибольшей интенсивностью при 1028, 745, и 726 см–1. Можно сказать, что в молекулах образца номер 2 содержится большее количество S=O-связей.
Анализ спектров Фурье-ИКС ДО. Прежде всего необходимо рассмотреть спектры Фурье-ИКС ДО (рис. 2 и табл. 3) [8, 20].
Таблица 3.
Положение пика, см–1 | Интерпретация | Образцы | Ссылка |
---|---|---|---|
491.8; 482.2 | − | 4 и 5 | |
543.0 | Валентные колебания связи C–S | 4 и 5 | [13] |
585–570 | Валентные колебания связи C–S | 1, 2, 3 и 4 | [13] |
725−745 | Внеплоскостные деформационные колебания Ar–H | 1, 2, 3 4, и 5 | [21] |
790−805 | Внеплоскостные деформационные колебания Ar–H | 1, 2, 3, 4 и 5 | [21] |
865−850 | Внеплоскостные деформационные колебания Ar–H | 1, 2, 3, 4 и 5 | [21] |
1015−1025 | R–SO–R | 1, 2, 3, 4 и 5 | [13] |
1359−1368 | R–SO2–N | 1, 2, 3, 4 и 5 | [13] |
1440−1450 | Деформация C–H-связи | 1, 2, 3, 4 и 5 | [21, 22] |
1585.4 | Полициклическая ароматика | 1, 2 и 3 | [21, 22] |
1597.0 | Валентные колебания ароматических C=C-связей | 5 | [21, 22] |
1687.7; 1689.6 | Валентные колебания карбонильных C=O-групп | 2 и 5 | [21, 22] |
1674.2 | Амид | 1 | [13] |
1880−1888 | Ar–OH | 1, 2, 3, 4 и 5 | [13] |
2358.0 | – | 4 | |
2719.0 | – | 1, 2, 3, 4 и 5 | |
2852.0 | Симметричные валентные колебания CH3-группы | 1, 2, 3, 4 и 5 | [21, 22] |
2905.0 | Валентные колебания связи C–H | 1, 2, 3, 4 и 5 | [21, 22] |
2923.1 | Асимметричные валентные колебания CH2-группы | 1, 2, 3, 4 и 5 | [21, 22] |
3045.0−3063.0 | Валентные колебания связи Ar–H | 1, 2, 3, 4 и 5 | [22] |
3404.0−3460.0 | Пиррол, свободное валентное колебание NH | 1, 2, 3, 4 и 5 | [13] |
3761.1; 3736.0 | − | 4 и 5 | |
3818.9; 3807.4 | − | 4 и 5 | |
3878.7; 3871.0 | − | 4 и 5 |
В области волновых чисел, равных 2923, 2905, 2852 и 1450 см–1, находятся полосы поглощения насыщенных углеводородов (см. рис. 2). Полосы в области между 3150 и 3000 см–1 соответствуют валентным колебаниям ароматических C−H-связей, а в диапазоне волновых чисел 900−700 см–1, в свою очередь, проявляются внеплоскостные деформационные колебания тех же связей. Наличие полос в интервале 1600 и 1580 см–1 обусловлено валентными колебаниями C=C-связей в ароматических соединениях [8]. Ниже приведены другие пики, характерные для всех образцов: пики сульфоксидных групп при 1015−1025 см–1; полосы возле 1359−1368 см–1, обусловленные наличием связи между S и N; свободное валентное колебание N−H-связи в пирроле в виде полосы при 3400−3410 см–1 [20] и полоса в диапазоне 1880−1888 см–1, принадлежащая связи Ar−OH. Полоса, лежащая возле 543.0 см–1 и наблюдаемая только в спектрах Фурье-ИКС ДО образцов 4 и 5, объясняется валентными колебаниями C−S. Исходя из спектра Фурье-ИКС ДО образца номер 1, можно судить о наличии амидной группы, дающей полосу при 1674.2 см–1. Принадлежность полос с низкой интенсивностью, лежащих между 3700 и 3880 см–1 в спектрах образцов 4 и 5 и пика около 2719.0 см–1, не выяснена.
Соотношения интегральных областей в спектрах Фурье-ИКС ДО асфальтенов были сопоставлены в пределах отдельно взятого спектра, но никак не между различными образцами (табл. 4). Широко распространенный подход обработки данных Фурье-ИК спектроскопии имеет дело именно с соотношениями, а не с абсолютными значениями [8].
Таблица 4.
Образец | Ι1 | Ι2 | Ι3 |
---|---|---|---|
1 | 10.3 | 2.9 | 0.3 |
2 | 10.6 | 3.0 | 0.3 |
3 | 12.7 | 3.1 | 0.3 |
4 | 6.5 | 3.4 | 0.5 |
5 | 7.0 | 3.3 | 0.5 |
Ι1 = ΙA (3000–2800)/ΙA (1650–1520) − показатель соотношения алифатических и ароматических групп, Ι2 = ΙA (3000–2800)/ΙA (900–700) − показатель доли алифатических и ароматических атомов водорода, Ι3 = ΙA (1650–1520)/ΙA (900–700) − показатель степени замещенности и/или сконденсированности ароматических колец.
Образцы 1, 2 и 3 имеют большее значение Ι1, чем асфальтены, выделенные из нефтей 4 и 5. Асфальтены тяжелых нефтей характеризуются бóльшим отношением алифатических групп к ароматическим, чем асфальтены, полученные из легких нефтей. Значение Ι2 для асфальтенов легких нефтей немного больше, чем для асфальтенов тяжелых нефтей. Высокое значение Ι3, характерное для асфальтенов легких нефтей, в отличие от асфальтенов тяжелых нефтей, предполагает большую степень замещенности ароматических соединений легких нефтей.
Анализ спектров комбинационного рассеяния. Рамановские спектры асфальтенов приведены ниже (см. рис. 3 (a, б) и табл. 5 (a, б)).
Таблица 5
Образец | Обозначение | Положение полосы, см–1 | Ширина, см–1 | Интегральная интенсивность | La, нм |
---|---|---|---|---|---|
1 | D1 | 1369 | 257 | 425 071 | 0.94 |
1558 | 108 | 126 347 | |||
G | 1589 | 64 | 91 042 | ||
2 | D1 | 1372 | 266 | 375 | 0.99 |
1544 | 97 | 66 | |||
G | 1588 | 70 | 105 | ||
1799 | 34 | 3 | |||
3 | D1 | 1368 | 266 | 367 | 0.83 |
D2 | 1562 | 104 | 99 | ||
G | 1594 | 58 | 69 | ||
4 | D1 | 1347 | 223 | 279 | 1.21 |
D3 | 1459 | 40 | 6 | ||
G | 1556 | 103 | 77 | ||
D2 | 1601 | 54 | 88 | ||
5 | D1 | 1336 | 205 | 212 | 1.64 |
1544 | 106 | 49 | |||
G | 1595 | 54 | 79 | ||
Таблица 5б. Анализ рамановских спектров (функция Лоренца) асфальтенов и расчетный диаметр ароматического фрагмента, La | |||||
Образец | Обозначение | Положение полосы, см–1 | Ширина, см–1 | Интегральная интенсивность | La, нм |
1 | D1 | 1363 | 200 | 504 978 | 1.03 |
1549 | 78 | 153 452 | |||
G | 1590 | 44 | 118 122 | ||
2 | D1 | 1349 | 184 | 445 504 | 0.82 |
D3 | 1477 | 114 | 109 484 | ||
1561 | 57 | 140 863 | |||
G | 1596 | 38 | 83 135 | ||
3 | D1 | 1349 | 187 | 376 | 0.84 |
1484 | 119 | 84 | |||
1566 | 57 | 106 | |||
G | 1599 | 36 | 72 | ||
4 | 1257 | 54 | 15 | 1.59 | |
D1 | 1339 | 164 | 285 | ||
D3 | 1458 | 164 | 86 | ||
G | 1580 | 59 | 103 | ||
D2 | 1605 | 34 | 59 | ||
5 | D1 | 1333 | 181 | 320 | 1.10 |
D3 | 1487 | 142 | 48 | ||
G | 1576 | 54 | 80 | ||
D2 | 1602 | 34 | 61 |
Полоса D1 при 1345 см–1 возникает вследствие искаженного внеплоскостными дефектами и гетероатомами вибрационного движения атомов углерода, находящихся в sp2-гибридном состоянии в ароматических кольцах. В свою очередь валентная вибрация sp2-гибридизованной C−C-связи в обычных полиароматических структурах проявляется в виде полосы G, предполагающей существование некоторого упорядочения ароматических слоев в асфальтенах [11]. Наличие полосы D2 при 1602 см–1 в виде небольшого пика, перекрывающегося со следующей интенсивной полосой G, можно объяснить внедрением графенового фрагмента или искажением микроструктуры, вызванным наличием гетероатома или атома металла [10]. Полоса, обозначенная как D3, представляет собой широкий пик углеродистого вещества с низкой степенью кристалличности, лежащий около 1500 см–1. Возникновение пика объясняется дефектами в ароматических слоях, таких как тетраэдрический атом углерода. В случае асфальтенов, возникновение полосы D3 в спектре говорит о неупорядоченности молекулярной упаковки [10, 12]. Полоса D1 непосредственно связана с вибрационными состояниями на границе микрокристаллических структур, например на периферии ароматических систем асфальтенов [23]. Поэтому при увеличении доли краевых атомов, интенсивность полосы D1 также растет. Использование логометрического метода позволяет оценить средний размер ароматических ядер асфальтенов [24]. Модель количественного описания с использованием уравнения, связывающего средний размер графитоподобного фрагмента и соотношение площадей полос G и D1, была предложена Tuinstra F. и Koenig J.L. [24]:
Соотношение алифатических и ароматических групп, как было показано ранее при анализе спектров Фурье-ИКС ДО, ниже для образцов асфальтенов 4 и 5. Аналогичные результаты получены при анализе спектров комбинационного рассеяния − средний молекулярный размер ароматических ядер больше для образцов 4 и 5, что говорит о согласованности результатов, полученных методом Фурье-ИКС ДО и рамановской спектроскопии. В целом, настоящие результаты Фурье-ИК спектроскопии диффузного отражения и спектроскопии комбинационного рассеяния отвечают модели Йена−Муллинса, используемой для описания асфальтенов легких, среднетяжелых и тяжелых нефтей.
Таким образом, анализ спектров Фурье-ИК НПВО показывает, что асфальтены в общем случае содержат ароматические углеводороды, карбонильные C=O-группы и сульфоксиды. Асфальтены 2, 4 и 5 имеют полосу около 750 см–1, отвечающую внеплоскостным деформационным колебаниям ароматических С−Н-связей, а полоса при 690 см–1, также характерная для внеплоскостных деформационных колебаний ароматических С−Н-связей, наблюдается только в спектрах Фурье-ИК НПВО образцов 3, 4 и 5. В сравнении с образцами 4 и 5 асфальтены с номерами 1, 2 и 3 характеризуются большим значением Ι1. Следовательно, асфальтены тяжелых и среднетяжелых нефтей отличаются бóльшим значением отношения количества алифатических групп к ароматическим, чем асфальтены легких нефтей. Величина Ι2 (соотношение числа алифатических и ароматических атомов водорода) для асфальтенов легких нефтей выше по сравнению с образцами 1, 2 и 3. Значение Ι3 для асфальтенов легких нефтей также больше, чем Ι3 в спектрах асфальтенов 1, 2 и 3. Степень замещения ароматических фрагментов в асфальтенах легких нефтей относительно выше, чем в тяжелых и среднетяжелых. Значения La асфальтенов среднетяжелых и тяжелых нефтей лежат в одном диапазоне, не зависящем от вида аппроксимирующей функции (гауссовой или лоренцевой). Для асфальтенов легких нефтей La на 54% больше, чем La асфальтенов тяжелых и среднетяжелых нефтей, при использовании функции Гаусса для обработки результатов спектроскопии комбинационного рассеяния и на 49% больше при использовании функции Лоренца. В спектрах Фурье-ИКС ДО отношение количества алифатических групп к ароматическим меньше для асфальтенов, выделенных из легкой нефти (образцы 4 и 5). Полученные результаты свидетельствуют в пользу модели Йена−Маллинса при описании химической структуры асфальтенов.
Список литературы
Barrera D.M., Ortiz D.P., Yarranton H.W. // Energy Fuels. 2013. V. 27. P. 2474.
Rashid Z., Wilfred C.D., Gnanasundaram N., Arunagiri A., Murugesan T. // J. Petrol. Sci. Eng. 2019. V. 176. P. 249.
Adams J.J. // Energy Fuels. 2014. V. 28. P. 2831.
Behbahani T.J., Miranbeigi A.A., Sharifi K., Behbahani Z.J. // Petrol. Chemistry. 2018. V. 58. P. 622.
AlHumaidan F.S., Hauser A., Rana M.S., Lababidi H.M.S., Behbehani M. // Fuel. 2015. V. 150. P. 558.
Majumdar R.D., Gerken M., Mikula R., Hazendonk P. // Energy Fuels. 2013. V. 27. P. 6528.
Wang J., Ferguson A.L. // J. Phys. Chem. B. 2016. V. 120. P. 8016.
AlHumaidan F.S., Hauser A., Rana M.S., Lababidi H.M.S. // Energy Fuels. 2016. V. 30. P. 2892.
Riley B. J., Lennard C., Fuller S., Spikmans V. // Forensic Sci. Int. 2016. V. 266. P. 555.
Bouhadda Y., Bormann D., Sheu E., Bendedouch D., Krallafa A., Daaou M. // Fuel. 2007. V. 86. P. 1855.
Riedeman J.S., Kadasala N.R., Wei A., Kenttamaa H.I. // Energy Fuels. 2016. V. 30. P. 805.
Abdallah W.A., Yang Y. // Energy Fuels. 2012. V. 26. P. 6888.
Mullins O.C. // Annual Rev. Anal. Chem. 2011. V. 4. P. 393.
Ok S., Mahmoodinia M., Rajasekaran N., Sabti M.A., Lervik A., van Erp T.S., Cabriolu R. // Energy Fuels. 2019. V. 33. P. 8259.
Huang J., Yuro R., Romeo G.A., Jr. // Fuel Sci. Technol. Int. 1995. V. 13. P. 1121.
Asemani M., Rabbani A.R. // Geosci. J. 2016. V. 20. P. 273.
Gabrienko A.A., Lai C.H., Kazarian S.G. // Energy Fuels. 2014. V. 28. P. 964.
Grdadolnik J. // Acta Chim. Slov. 2002. V. 49. P. 631.
Petrova L.M., Abbakumova N.A., Foss T.R., Roma-nov G.V. // Petrol. Chemistry. 2011. V. 51. P. 262.
Pretsch E., Buhlmann P., Baderstscher M. Structure Determination of Organic Compounds. Berlin: Springer, 2009. 4th Ed.
Durand E., Clemancey M., Lancelin J., Verstraete J., Espinat D., Quoineaud A. // Energy Fuels. 2010. V. 24. P. 1051.
Maqbool T., Srikiratiwong P., Fogler H.S. // Energy Fuels. 2011. V. 25. P. 694.
Jawhari T., Roid A., Casado J. // Carbon. 1995. V. 33. P. 1561.
Tuinstra F., Koenig J.L. // Chem. Phys. 1970. V. 53. P. 1126.
Daaou M., Bendedouch D., Bouhadda Y., Vernex-Loset L., Modaressi A., Rogalski M. // Energy Fuels. 2009. V. 23. P. 5556.
Дополнительные материалы отсутствуют.