CC BY
126
УДК 622.692.23
ИССЛЕДОВАНИЕ ПРОЦЕССА СТРУЙНОГО РАЗМЫВА ДОННЫХ ОТЛОЖЕНИЙ
В НЕФТЯНЫХ РЕЗЕРВУАРАХ
Бутов Владимир Григорьевич1,
bvg@niipmm.tsu.ru
Никульчиков Андрей Викторович1,
nikulchikov@sibmail.com
Никульчиков Виктор Кенсоринович1,
nikulchikov@tpu.ru
Солоненко Виктор Александрович2,
vik@niipmm.tsu.ru
Ящук Алексей Александрович2,
rainbow@niipmm.tsu.ru
1 Национальный исследовательский Томский политехнический университет, Россия, 634050, г. Томск, пр. Ленина, 30.
2 Национальный исследовательский Томский государственный университет, Россия, 634050, г. Томск, пр. Ленина, 36.
Актуальность. При эксплуатации нефтяных резервуаров часто возникает проблема отложения твердых частиц и тяжелых пара-финистых осадков на днище резервуара. Выпадение данных осадков влечет за собой уменьшение полезного объема резервуара, препятствует свободному перемешиванию слоев нефти, что в свою очередь способствует концентрации агрессивных растворов солей в районе днища и развитию коррозионных разрушений нефтяных резервуаров. Для предотвращения подобных ситуаций предусмотрен комплекс мер (например, ручная и механизированная очистка резервуаров с использованием горячей воды и химических реагентов), однако существенный интерес представляют устройства, препятствующие выпадению донных осадков без остановки эксплуатации резервуара, например винтовые устройства, размывающие донные отложения с помощью турбулентной струи.
Цель исследования: расчет гидродинамических характеристик турбулентной затопленной струи в резервуаре при различных температурах нефти; определение реологических параметров нефти, при которых возможно размытие донных отложений в резервуаре.
Объекты исследования: конструкция винтового устройства размыва донных отложений в резервуарах с нефтью, влияние основных режимных и конструктивных характеристик винта устройства на параметры турбулентной струи, моделирование процесса гидродинамического течения нефти в замкнутом объеме резервуара. Методы: метод конечных объемов с использованием ANSYS CFX.
В результате расчетов конструкций показаны зависимость скорости размыва от расстояния до винта для разных значений температур и соответствующей им вязкости нефти, показаны условия, при которых возможен размыв донных отложений в резервуаре для данной геометрии гребного винта.
Ключевые слова:
Нефть, реологическая модель, метод конечных объемов, касательные напряжения, продольная скорость, скорость размыва.
Введение
При эксплуатации нефтяных резервуаров часто возникает проблема отложения твердых частиц и тяжелых парафинистых осадков на днище резервуара [1]. Выпадение данных осадков влечет за собой уменьшение полезного объема резервуара, препятствует свободному перемешиванию слоев нефти, что в свою очередь способствует концентрации агрессивных растворов солей в районе днища и развитию коррозионных разрушений нефтяных резервуаров [2-4]. Для предотвращения подобных ситуаций предусмотрен комплекс мер (например, ручная и механизированная очистка резервуаров с использованием горячей воды и химических реагентов [5-10]), однако существенный интерес
представляют устройства, препятствующие выпадению донных осадков без остановки эксплуатации резервуара, например винтовые устройства, размывающие донные отложения с помощью турбулентной струи [11-13].
При работе устройства размыва при вращении лопастей гребного винта создается направленная затопленная турбулентная струя жидкости. Под воздействием затопленной турбулентной струи от лопастей винта нефть в резервуаре перемешивается, имеющиеся донные отложения размываются. Устройство для размыва донных отложений для интенсификации процесса перемешивания выполнено с возможностью поворота вала винта на угол 60°. При работе устройства в крайних точках в ре-
зервуаре возникает круговое движение всей массы нефти в резервуаре. Цель работы:
1) расчет гидродинамических характеристик турбулентной затопленной струи в резервуаре при различных температурах нефти;
2) определение реологических параметров нефти, при которых возможно размытие донных отложений в резервуаре.
Методы
При проведении расчетов трехмерной задачи применялся метод конечных объемов с использованием программного комплекса ANSYS CFX. Уравнения модели течения (1)-(5), реализованные в комплексе ANSYS CFX, представлены в руководствах по гидродинамике, например, в [14-17].
В рамках этого комплекса были построены конечно-объемные модели [18, 19], воспроизводящие геометрию рассматриваемых расчетных областей и позволяющие проводить комплексные расчеты с учетом пространственного турбулентного характера течения.
Для моделирования была использована конечно-объемная расчетная сетка с 1134867 элементами и 510450 узлами для расчета течения в резервуаре. Для моделирования течения в окрестности винта была использована сетка с числом элементов и узлов 394627 и 168245 соответственно.
Для описания турбулентного течения были применены эйлеровы переменные х, ¿=1,2,3. Течение в данной точке пространства во времени характеризуется для жидкости вектором скорости и с компонентами щ (¿=1,2,3) и давлением р. Температура Т и плотность жидкости р считаются постоянными величинами.
Уравнение неразрывности для осредненных по времени величин:
§+ Р > = °
дt дхк
Уравнение сохранения импульса:
d(put) д
dp
д%.
J~L~ +-(PUMk ) = —— + g + —k
dt dxk dx! dxk
(1)
(2)
где й - массовая сила в ¿-м направлении, в случае расчета д1=д2=0, £3=р£, д - ускорение силы тяжести.
Тензор вязких напряжений без учета влияния градиентов давления имеет вид
(
T!j =»
du. du , —- + —-
dXj dx j
\
- 2 ,5U ^
3 dxk '
(3)
где р - коэффициент вязкости, зависящий от температуры; 8ц - символ Кронекера.
Уравнение для энергии турбулентного перемешивания:
d(pk) d d
+ — (puk) =-(или) -
dt dx. dx.
d (,
Pak+dx t) ~ ' '
dk } dx
(4)
Уравнение для удельной скорости диссипации:
d(pm) d . Y(o d ,
——- +-{pu m) =l--(uizi.) -
dt dx t ' k dx, ' j
d (,
dk }
+ dx l +ct ^t ) dk ' '
(5)
Здесь к - кинетическая энергия турбулентности; а - частота турбулентных пульсаций.
Характеристики нефти при различных температурах приняты согласно [20].
Нефть имеет ньютоновские свойства при температуре свыше 20 °С включительно. Поскольку нефть имеет неньютоновские свойства при температуре ниже 20 °С, для ее описания при этих значениях использовалась модель вязкости Бингама-Шведо-ва [20].
Геометрические и режимные параметры устройства задавались следующим образом:
• диаметр винта, мм - 690 (определялся размерами люка резервуара);
• количество лопастей винта, шт - 3;
• частота вращения винта, об/мин - 460;
• угол поворота вала винта, угл. град, не менее 60. Расчеты проведены в диапазоне температур
нефти от +10 до +40 °С.
Вязкость нефти задавалась из диапазона 5-100 сСт.
В расчетах использовался нефтяной резервуар РВСП:
• объем, м3 - 50000;
• диаметр, м - 60,7;
• высота стенки, м - 18,1.
Обсуждение полученных результатов
Моделирование поведения затопленной турбулентной струи, формируемой винтом Результаты моделирования турбулентной струи, формируемой винтом, с использованием полной 3-Б модели винта при высоких температурах нефти (выше 20 °С) представлены на рис. 1.
При частоте вращения винта 460 об/мин продольная скорость непосредственно за винтом составляет до 4,9 м/с (рис. 1, А).
В случае низких температур (менее 15 °С) нефть проявляет неньютоновские свойства, результаты аналогичного расчета для этого случая представлены на рис. 2.
На рис. 3 приведена зависимость скорости размыва от расстояния до винта для разных значений
Рис. 1. Распределение продольной скорости при вращении винта с частотой 460 об/мин в нефти, подчиняющейся закону Ньютона в диапазоне температур свыше 20 °С: в локальной окрестности винта (А) и в неограниченном (справа) пространстве (Б)
Fig. 1. Distribution of axial velocity in Newtonian oil within the temperature range higher than 20 °С at propeller rotation frequency 460 rpm: in propeller vicinity (A); in unrestricted (on the right) space (Б)
Прощальная скорость, м/с
m m
% ^ % % <
Рис. 2. Распределение продольной скорости при вращении винта с частотой 460 об/мин в нефти, подчиняющейся неньютоновскому закону при температуре 10 °С: в локальной окрестности винта (А) и в неограниченном (справа) пространстве (Б)
Fig. 2. Distribution of axial velocity in oil that exhibits non-Newtonian characteristics at temperature 10 °С and propeller rotation frequency 460 rpm: in propeller vicinity (A); in unrestricted (on the right) space (Б)
температур и соответствующих им вязкостей нефти. Из результатов видно, что в случае соответствия модели течения ньютоновскому приближению скорость размыва на расстоянии 60,7 м от винта в случае неограниченного объема нефти будет составлять 0,6 м/с при вязкости, большей 30 сСт.
Расстояние сгт винта, м'
Рис. 3. Зависимость скорости размыва от расстояния до винта для разных значений температур и соответствующих им вязкостей: сиреневому цвету соответствует температура 20 °С или вязкость 30 сСт (ньютоновская модель); зеленому цвету соответствует температура 40 °С или вязкость 9 сСт (ньютоновская модель); синему цвету соответствует температура 15 °С или вязкость 50 сСт (неньютоновская модель); красному цвету соответствует температура 10 °С или вязкость 100 сСт (неньютоновская модель)
Fig. 3. Dependence of resuspension velocity on the distance from the propeller at various temperatures and corresponding oil viscosities: magenta color indicates temperature 20 °С or viscosity 30 cSt (Newtonian model); green color indicates temperature 40 °С or viscosity 9 cSt (Newtonian model); blue color indicates temperature 15 °С or viscosity 50 cSt (non-Newtonian model); red color indicates temperature 10 °С or viscosity 100 cSt (non-Newtonian model)
Моделирование течения в резервуаре
При моделировании течения внутри ограниченного объема при наличии одного устройства размыва донных отложений при варьировании поворота этого устройства в пределах ±30° от направления на центр резервуара область, где скорость струи позволит осуществить размыв, будет представлена следующими рисунками.
На рис. 4 различными цветами представлены результаты определения области размыва при различных температурах и соответствующих им вяз-костях. Графически скорость размыва в зависимости от расстояния до винта для различных температур и соответствующих им вязкостей показана на рис. 5. На рис. 6 представлены линии тока в объеме резервуара и профили скорости вдоль оси винта при нулевом отклонении оси винта для случая разных реологических моделей жидкости.
' МО 22ХГ
Рис. 4. Область размыва в зависимости от температуры (при 20, 15 и 10 °С)
Fig. 4. Resuspension zone for various temperatures (20, 15 and 10 °С)
Рис. 5. Зависимость скорости размыва от расстояния до винта: 1 - при вязкости 30 сСт в резервуаре (в рамках ньютоновской реологической модели); 2 - при вязкости 50 сСт в резервуаре (в рамках неньютоновской реологической модели)
Fig. 5. Dependence of resuspension velocity on the distance from the propeller: 1 - at viscosity 30 cSt in tank (within Newtonian rheological model); 2 - at viscosity 50 cSt in tank (within non-Newtonian rheological model)
Выводы
По результатам моделирования можно сделать
следующие выводы:
1) при заданных параметрах гребного винта при температурах нефти свыше 20 °С и указанного диапазона вязкости нефти возможно формирование потока жидкости со скоростью 0,5 м/с на расстоянии 60,7 м от лопастей винта в случае неограниченного объема нефти;
2) при температуре нефти ниже 15 °С для достижения требуемых параметров в случае неограни-
Рис. 6. Линии тока для вязкости 30 сСт в рамках ньютоновской реологической модели (А) и 100 сСт в рамках неньютоновской реологической модели (Б)
Fig. 6. Lines of flow with viscosity 30 cSt within Newtonian rheological model (А) and 100 cSt within non-Newtonian rheological model (Б)
ченного объема нефти размыва достичь не удается (скорость 0,5 м/с на расстоянии 60,7 м от лопастей винта);
3) при размыве донных отложений в ограниченном объеме (резервуаре) для устранения застойных зон, находящихся вблизи стенок резервуара, в обязательном порядке необходимо применение нескольких устройств, устано-
СПИСОК ЛИТЕРАТУРЫ
1. Characterization of the sludge deposits in crude oil storage tanks / I.A. Gopang, H. Mahar, A.S. Jatoi, K.S. Akhtar, M. Omer, M.S. Azeem // Journal of Faculty of Engineering & Technology. - 2016. - V. 23 (1). - P. 57-64.
2. Рыбалко П.Н., Прибытков П.А., Бурков П.В. Исследование напряженно-деформированного состояния днища резервуара с коррозионным поражением // Проблемы геологии и освоения недр: Труды XVIII Международного симпозиума имени академика М.А. Усова студентов и молодых ученых - Томск: ТПУ, 2014. - Т. II. - С. 624-625.
3. Гималетдинов Г.М., Саттарова Д.М. Способы очистки и предотвращения накопления донных отложений в резервуарах // Нефтегазовое дело. - 2006. - Т. 1. URL: http://ogbus.ru/artic-le/sposoby-ochistki-i-predotvrashheniya-nakopleniya-donnyx-otlozhenij-v-rezervuarax (дата обращения 13.06.2018).
4. Experimental investigations of various methods of sludge measurements in storage oil tanks / M. Monteiro, V. Svet, D. Sandi-lands, S. Tsysar // Advances in Remote Sensing. - 2015. - V. 4. -№2. - P. 119-137.
5. Чурикова Л.А., Конашева Е.А., Утегалиев А.Т. Обзор современных методов очистки резервуаров от нефтяных остатков // Технические науки в России и за рубежом: материалы V Меж-дунар. науч. конф. - М., январь 2016 г. - М.: Буки-Веди, 2016.- С. 71-75.
6. Hu G., Li J. Zeng G. Recent development in the treatment of oily sludge from petroleum industry: a review // Journal of hazardous materials. - 2013. - V. 261. - P. 470-490.
7. Method for removing sludge in crude oil tank and recovering oil therefrom: patent US6673231B2, United States, Publ. 06.01.2004.
вленных симметрично по периметру первого яруса резервуара. Количество устройств размыва донных отложений зависит от объема резервуара, свойств нефти и может быть определено расчетом.
Исследования проводились в рамках выполнения государственного задания Министерства образования и науки России, проект № 9.9752.2017/8.9.
8. Method for cleaning storage tanks: patent US20070283981A1, United States, Publ. 13.12.2007.
9. Oil tank sludge removal method: patent US20060042661A1, United States, Publ. 2.03.2006.
10. Water washing method and system for crude oil tank: patent US20170291200A1, United States, Publ. 12.10.2017.
11. Костин М.В., Демин А.В., Садчиков А.В. Очистка емкостей от остатков нефти и нефтепродуктов с помощью эффлюэнта // Научное обозрение. Биологические науки. - 2015. - № 1. -С. 93.
12. Чепур П.В., Тарасенко А.А. Особенности совместной работы резервуара и устройств размыва донных отложений винтового типа // Фундаментальные исследования. - 2015. - №2. -С. 1671-1675.
13. Nekrasov V.O. Modeling bottom sediment erosion process by swirling the flow by tangential supply of oil in the tank // IOP Conference Series: Materials Science and Engineering: Transport and Storage of Hydrocarbons. - 2016. - V. 154. 7 p. URL: http://iopscience.iop.org/article/10.1088/1757-899X/ 154/1/012003 (дата обращения 13.06.2018).
14. Флетчер К. Вычислительные методы в динамике жидкостей: В 2-х т. Т. 1. - М.: Мир, 1991. - 504 с.
15. Андерсон Д., Таннехил Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. - М: Мир, 1990. - 726 с.
16. Уилкокс Д.К. Многомасштабная модель турбулентных течений // Аэрокосмическая техника. - 1989. - № 11. - С. 47-60.
17. Численное решение многомерных задач газовой динамики / под ред. С.К. Годунова. - М.: Наука, 1976. - 400 с.
18. Применение системы ANSYS к решению задач механики сплошной среды. Практическое руководство / под ред. проф. А.К. Любимова. - Нижний Новгород: Изд-во Нижегородского госуниверситета, 2006. - 227 с.
19. Burkov P.V., Burkova S.P., Samigullin V.D. Computer Simulation of Stress-Strain State of Oil Gathering Pipeline Designed for Ugut Field // IOP Conference Series: Materials Science and Engineering. - 2016. - V. 125: Materials Treatment: Current Problems and Solutions. - 7 p.
20. Трубопроводный транспорт нефти и газа / Р.А. Алиев, В.Д. Белоусов, А.Г. Немудров, В.Л. Юфин, Г.И. Яковлев. - М.: Недра, 1988. - 368 с.
Поступила 15.06.2018 г.
Информация об авторах
Бутов В.Г., доктор физико-математических наук, профессор отделения ядерного топливного цикла Инженерной школы ядерных технологий Национального исследовательского Томского политехнического университета.
Никульчиков А.В., кандидат физико-математических наук, доцент отделения нефтегазового дела Инженерной школы природных ресурсов Национального исследовательского Томского политехнического университета.
Никульчиков В.К., кандидат технических наук, доцент отделения нефтегазового дела Инженерной школы природных ресурсов Национального исследовательского Томского политехнического университета.
Солоненко В.А., кандидат физико-математических наук, заведующий лабораторией оптимизации отдела математической физики Научно-исследовательского института прикладной математики и механики Томского государственного университета.
Ящук А.А., кандидат физико-математических наук, старший научный сотрудник лаборатории оптимизации отдела математической физики Научно-исследовательского института прикладной математики и механики Томского государственного университета.
UDC 622.692.23
SIMULATION STUDY OF BOTTOM SEDIMENTS JET EROSION IN OIL TANK
Vladimir G. Butov,
bvg@niipmm.tsu.ru
Andrey V. Nikulchikov1,
nikulchikov@sibmail.com
Viktor K. Nikulchikov1,
nikulchikov@tpu.ru
Victor A. Solonenko2,
vik@niipmm.tsu.ru
Alexey A. Yashchuk2,
rainbow@niipmm.tsu.ru
1 National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia.
2 National Research Tomsk State University, 36, Lenin Avenue, Tomsk, 634050, Russia.
The relevance. When crude oil is stored in tanks, invariably solid particles and heavy paraffin sediments are deposited. The precipitation of these sediments entails a decrease in reservoir useful volume, prevents free mixing of oil layers, which boosts the concentration of aggressive salt solutions in the bottom region and the development of corrosive destruction of oil reservoirs. A set of measures is envisaged to prevent such situations (including manual and mechanized cleaning of tanks using hot water and chemical reagents), however, the devices that prevent precipitation of bottom sediments without stopping the operation of the reservoir are of significant interest, for example, screw devices that dilute bottom sediments with the help of turbulent jet.
The aim of the research is calculation of hydrodynamic properties of turbulent submerged jet at various oil temperatures; definition of oil rheological parameters which enable the lift off, mix, and re-suspension of sediments in a tank
Objects of research: design of the bottom sediment erosion device for oil tanks; influence of the screw regime and design parameters on the turbulent jet properties; simulation of oil hydrodynamic flow in a closed reservoir volume. Methods: finite volume method in ANSYS CFX.
Results. The paper demonstrates the dependence of sediment resuspension velocity on the distance from the screw at various temperatures and corresponding oil viscosities. The conditions to re-suspend oil tank bottom sediments were defined for the given propeller geometry It was proved that the screw-type device is an effective means to lift solids off the storage tank floor and ensure required characteristics of tank oil.
Key words:
Oil, rheological model, finite volume method, shear stress, longitudinal velocity, resuspension velocity.
The research was carried out within the State Task of the Ministry of Science and Education of the Russian Federation,project no. 9.9752.2017/8.9.
REFERENCES
1. Gopang I.A., Mahar H., Jatoi A.S., Akhtar K.S., Omer M., Az-eem M.S. Characterization of the sludge deposits in crude oil storage tanks. Journal of Faculty of Engineering & Technology, 2016, vol. 23 (1), pp. 57-64.
2. Rybalko P.N., Pribytkov P.A., Burkov P.V. Issledovanie naprya-zhenno-deformirovannogo sostoyaniya dnishcha rezervuara [Study of stress-strain behavior of the corroded oil tank bottom surface]. Tomsk, TPU Publ., 2014, vol. 2, pp. 624-625.
3. Gimaletdinov G.M., Sattarova D.M. Sposoby ochistki i predot-vrasheniya nakopleniya donnykh otlozheniy v reservuarakh [Methods of cleaning and preventing accumulation of bottom sediments in oil tanks]. Neftegazovoe delo, 2006, vol. 1. Available at: http://ogbus. ru / article/sposoby-ochistki-i-predotvrashheni-ya-nakopleniya-donnyx-otlozhenij-v-rezervuarax (accessed 13 June 2018).
4. Monteiro M., Svet V., Sandilands D., Tsysar S. Experimental investigations of various methods of sludge measurements in storage oil tanks. Advances in Remote Sensing, 2015, vol. 4, no. 2, pp. 119-137.
5. Churikova L.A., Konasheva E.A., Utegaliev A.T. Obzor sovre-mennykh metodov ochistki reservuarov ot neftyanykh ostatkov [Modern methods of cleaning oil tanks: a review]. Tekhnicheskie nauki v Rossii i za rubezhom. Materialy V Mezhdunarodnoy nauchnoy konferentsii [Technical sciences if Russia and abroad: Proc. of International scientific conference]. Moscow, January 2016. Moscow, Buki-vedi Publ., 2016, pp. 71-75.
6. Hu G., Li J., Zeng G. Recent development in treatment of oily sludge from petroleum industry: a review. Journal of hazardous materials, 2013, vol. 261, pp. 470-490.
7. Kim M.-J., Yoon D.-S., Choi C.-H., Choi J.-T., Kwak J.-W., Park S.-K. Method for removing sludge in crude oil tank and recovering oil therefrom. Patent United States, US6673231B2, 2004.
8. Stewart T.E. Method for cleaning storage tanks. Patent United States, US20070283981A1, 2007.
9. Meyer D., Brons G., Perry R., Wildemeersch S., Kennedy R. Oil tank sludge removal method. Patent United States, US20060042661A1, 2006.
10. Son S.G. Water washing method and system for crude oil tank. Patent United States, US20170291200A1, 2017.
11. Kostin M.V., Demin A.V., Sadchikov A.V. Tank cleaning residues from oil and oil products through the effluent. Scientific review. Biological sciences, 2015, vol. 1, p. 93. In Rus.
12. Chepur P.V., Tarasenko A.A. The features of washout sediments devices screw type and tank collaboration. Fundamental research, 2015, vol. 2, pp. 1671-1675. In Rus.
13. Nekrasov V.O. Modeling bottom sediment erosion process by swirling the flow by tangential supply of oil in the tank. IOP Conference Series: Materials Science and Engineering: Transport and Storage of Hydrocarbons, 2016, vol. 154, 7 p. Available at: http://iopscience.iop.org/article/10.1088/1757-899X/ 154/1/012003 (accessed 13 June 2018).
14. Fletcher C.A.J. Computational Techniques for Fluid Dynamics. Berlin, Springer-Verlag, 1988. 409 p.
15. Anderson D.A., Tannehill J.C., Pletcher R.H.I. New York, McGraw-Hill, 1984, 609 p.
16. Wilcox D.C. Multiscale Model for Turbulent Flows. AIAA Journal, 1988, vol. 26, no. 11, pp. 1311-1320.
17. Chislennoe reshenie mnogomernykh zadach gazovoy dinamiki [Numerical solution of multidimensional problems in gas dynamics] Ed. by S.K. Godunov. Moscow, Nauka Publ., 1976, 400 p.
18. Primenenie sistemyANSYS k resheniyu zadach mekhaniki splosh-noy sredy. Prakticheskoe rukovodstvo [Use of ANSYS software in solving problems of continuum mechanics]. Ed. by A.K. Lyubi-mov. Nizhny Novgorod, Nizhny Novgorod State University Press, 2006, 227 p.
19. Burkov P.V., Burkova S. P., Samigullin V. D. Computer Simulation of Stress-Strain State of Oil Gathering Pipeline Designed for Ugut Field. IOP Conference Series: Materials Science and Engineering, 2016, vol. 125: Materials Treatment: Current Problems and Solutions, 7 p.
20. Aliev R.A., Belousov V.D., Nemudrov A.G., Truboprovodny transport nefti i gaza [Oil and gas pipeline transport]. Moscow, Nedra Publ., 1988, 368 p.
Received: 15 June 2018.
Information about the authors
Vladimir G. Butov, Dr. Sc., professor, National Research Tomsk Polytechnic University.
Andrey V. Nikulchikov, Cand. Sc., associate professor, National Research Tomsk Polytechnic University.
Viktor K. Nikulchikov, Cand. Sc., associate professor, National Research Tomsk Polytechnic University.
Victor A. Solonenko, Cand. Sc., head of laboratory, Tomsk State University.
Alexey A. Yashchuk, Cand. Sc., senior research assistant, Tomsk State University.
