ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2020;(8):52-61 ОРИГИНАЛЬНАЯ СТАТЬЯ / ORIGINAL PAPER
УДК 622.240.8 DOI: 10.25018/0236-1493-2020-8-0-52-61
ГЕОМЕХАНИЧЕСКИЕ МОДЕЛИ СДВИГОВОГО РАЗРУШЕНИЯ МНОГОСЛОЙНЫХ ГОРНЫХ ПОРОД
М.М. Алиев1, З.Ф. Исмагилова1, Н.Н. Бурмистрова1
1 Альметьевский государственный нефтяной институт, Альметьевск, Республика Татарстан, Россия, e-mail: [email protected]
Аннотация: Исследованы прочность и разрушения многослойных горных пород-сланцев, аргиллитов и др. Критерии, оценивающие состояния прочности и разрушения, разработаны в двух вариантах, как модификация теории Мора. Первый вариант критерия обобщает линейный закон Кулона в виде линейной огибающей предельных кругов Мора, а во втором варианте огибающая принята в виде параболической кривой. Решен ряд задач прочности многослойных пород в различных сложных напряженных состояниях. К многослойным породам относятся горные породы, образованные путем уплотнения чередующихся слоев в период многовекового осадконакопления. Горные породы, имеющие многослойные залегания, встречаются на различных глубинах земной коры, находясь в основном в условиях трехосного неравномерного напряженного состояния. Условия залегания многослойных пород, таких как аргиллиты, сланцы, обуславливают изменчивость их свойств по направлениям. Проектирование подземных сооружений, расположенных в различных породах, имеющих чередующиеся слои, осуществляется с определением местных напряжений, а устойчивое состояние сооружений проверяется с помощью соответствующих критериев прочности.
Ключевые слова: слоистость горных пород, геомеханическая модель, керновый материал, критерии прочности, испытание образцов, степень анизотропии, «бразильский тест», устойчивость скважины.
Для цитирования: Алиев М. М., Исмагилова З. Ф., Бурмистрова Н. Н. Геомеханические модели сдвигового разрушения многослойных горных пород // Горный информационно-аналитический бюллетень. - 2020. - J№ 8. - С. 52-61. DOI: 10.25018/0236-1493-2020-8-0-52-61.
Geomechanical models of shear fracture of multi-layer rocks
M.M. Aliev1, Z.F. Ismagilova1, N.N. Burmistrova1
1 Almetyevst State Oil Institute, Almetyevsk, Tatarstan, Russia, e-mail: [email protected]
Abstract: The strength and fracture of multi-layer rocks such as shale, clay stone, etc. are studied. The strength and fracture criteria are developed in two variants as a modification of the Mohr theory. The first variant generalizes the linear law of Coulomb in the form of a linear Mohr's envelope, which is assumed as a parabolic curve in the second variant. A series of problems on strength of multi-layer rocks in different stress conditions is solved. The multi-layer rocks are formed in compaction of alternate layers during many-centuries sedimentation. The multi-layer rocks occur at various depths mainly in the conditions of three-axial nonuniform stress state. Occurrence conditions of multi-layer rocks such as shale and clay stone govern
© М.М. Алиев, З.Ф. Исмагилова, Н.Н. Бурмистрова. 2020.
variability of their properties in different directions. Designing of underground structures in alternate-layer rocks involves determination of local stresses, and the stability of the structures is tested using the relevant strength criteria.
Key words: rock stratification, geomechanical model, core, strength criteria, sample testing, degree of anisotropy, Brazilian test, borehole stability.
For citation: Aliev M. M., Ismagilova Z. F., Burmistrova N. N. Geomechanical models of shear fracture of multi-layer rocks. MIAB. Mining Inf. Anal. Bull. 2020;(8):52-61. [In Russ]. DOI: 10.25018/0236-1493-2020-8-0-52-61.
В работе рассматриваются вопросы, связанные с геомеханическими характеристиками многослойных пород, на основе разработанных двух критериев прочности и разрушения. Критерии получены в виде модификации (обобщения) классического критерия прочности Отто Мора. Первый из них имеет линейную зависимость между касательными и нормальными напряжениями на возможных площадках скольжения и является обобщением линейного критерия Кулона-Мора, а второй эту зависимость принимает в квадратичной форме.
Прочностные и другие свойства многослойных пород могут быть сформулированы некоторыми непрерывными функциями координат. В этом случае анизотропию свойств таких пород можно принять упорядоченной, тем самым упростить геомеханическую модель.
Такая модель была предложена в работе [1], где область ее приложения относилась к гражданскому строительству для расчета анизотропных оснований сооружений. Позже она была применена для оценки устойчивости скважин и подземных хранилищ газа в слоистых породах. Сцепление породы принимается переменным значением, которое зависит от направления, а угол внутреннего трения является величиной постоянной. В усложненном варианте, где сцепление и угол внутреннего трения принимаются изменяющимися по направлениям, модель была обобщена критерием прочности Кулона-Мора.
В работе [4] для расчета устойчивости стенок скважин, пробуренных в анизотропных породах, использован линейный критерий Кулона-Мора, записанный для случая выполнения его при сдвиге вдоль напластования и по основной породе.
Другие модели разрушения анизотропных пород предложены в работах [5 — 8]. В работах [7, 8] рассматривается теория пластичности анизотропной среды на основании принятой математической модели деформирования.
Линейная геомеханическая
модель разрушения
Введем некоторый угол у (рис. 1), который является углом между нормалью к некоторой площадке и осью х декартовой системы координат. Пусть ось совпадает с направлением наименьшего сопротивления породы, имеющей ортогональную анизотропию.
ния в системе координат
Fig. 1. Lay-out of slickenside in coordinate system
Для разработки модели принимаем, что причиной разрушения породы является сдвиг по некоторой площадке с нормалью n, а закономерность Кулона-Мора имеет вид
тп = с к(у) + c(y), (1)
где
*n = 0,5 (x +а у W-^xy C°S2 V I (2)
an = 0,5 ( x + а у ) + 0,5 ( - а у )cos2 W + Txy sin2
тп и сп — касательное и нормальное напряжение на площадке сдвига; c(y) и к(у) — сцепление и тангенс угла внутреннего трения, которые являются функциями угла между нормалью к площадке сдвига и осью x.
При определенном сочетании компонентов напряжения на площадке сдвига выполняется условие
(d/dj [Тп + С к(у) — c(y)] = 0. (3)
Решая совместно (1), (2) и (3), получим три уравнения для определения с , с и т . Примем
ху ^ ^
y ^ = а±пЛ,
т =-П, Г
xy 2 J
где
П1 = {s[sin(2y — р) — 0,5р' secp sin2y] + 0,5s' cos(2y — р)}х,
П2 = {s[cos(2y — р) — 0,5р' secp sin2y] + 0,5s' sin(2y — р)}х,
X = XM = смрЛ1^^ с = 0,5(cx + a^ s = s(^) = c(y) + к^С
s' = s'(y) = с'(у) + к'(у)с.
Возведя в квадрат разность сх — с и тх , и суммируя эти результаты, получим критерий прочности (разрушения) слоистой среды с ортогональной анизотропией.
(Сх — Су)2 + 4тху2 = 4(Ц2 + П22). (4)
В случае ярко выраженной слоистости, когда порода сложена из чередующихся тонких слоев сдвиговыми характеристиками, могут быть приняты сцепление и тангенс угла внутреннего трения вдоль (с , р90) и поперек (с0, р0) напластова-Таблица 1
Результаты испытания образцов керна на одноосное (неограниченное) сжатие Unrestrained uniaxial compression test data of core
Горная порода ct, МПа
угол напластования, p°
0° 15° 30° 40° 45° 55° 60° 65° 75° 90°
Кварц филлит Quartz PhyLLite [9] 101,64 82,76 59,48 51,72 48,27 45,41 45,51 46,86 51,72 65,6
Углеродистый филлит Carbona PhyLLite [9] 79,53 68,96 48,58 44,13 41,37 35,86 37,42 38,51 48,27 65,14
Гнейс A / Gneiss A [9] 39 30 27 28 30 34 38 42 52 68
Мрамор / MarbLe [9] 77 73 71 71,5 72 74 77 79 83 91
Результаты одноосного сжатия образцов кернов разных пород Uniaxial compression test data of different rock core
Горная порода alt МПа
P°
0° 15° 30° 45° 60° 75° 90°
Анжуйский сланец Angers Schist [9] 148,33 70,22 36,40 10,34 12,10 25,60 97,22
Пенринский сланец Penrhyn Slate [9] 208,51 126,92 86,12 62,70 43,82 113,32 187,35
Боссьерский сланец Bossier shale [9] 113,55 82,75 45,35 38,36 12,48 27,58 42,34
Гнейс B / Gneiss B [9] 39 23 18 23 36 62 85
Хлористый сланец Chlorite schist [9] 110 62 51 53 68 100 113
ния, а закономерности изменения с(у) и tgp(у) в самом простом виде
с(у) = с0 + с90 эт2^ (5) 1др(у) = tgp0 соэ2у + tgp90 эт2у, (6)
Определение характеристик прочности, входящих в критерий (4)
Для определения характеристик прочности были использованы результаты испытания образцов керна на одноосное (неограниченное) сжатие под различны-
ми углами действия сжимающего напряжения относительно слоев. Эти данные, взятые для некоторых видов горных пород из источника [9], представлены в табл. 1 и табл. 2, на рис. 2 и рис. 3. Результаты подобных испытаний для разных пород приведены в работах [11 — 14].
Результаты, отмеченные в таблицах жирным шрифтом, соответствуют экспериментам. Остальные результаты восстановлены приблизительно из графиков зависимостей, заимствованных из соответствующих источников.
О 15 30 40 45 55 60 65 75 90
И Quartz Phyllite й Carbona Phyllite Ш Gneiss A < Marble
Рис. 2. График зависимости между ст^и в при одноосном сжатии Fig. 2. Cross-plot of ст1 and в in uniaxial compression
15 30
—Angers Schist —Gneiss В
Рис. 3. График зависимости между с и ß при одноосном неограниченном сжатии Fig. 3. Cross-plot of and ß in unrestrained uniaxial compression
В первом приближении предполагаем, что сдвиг происходит по плоскостям, вдоль которых угол внутреннего трения постоянный. При этом критерий прочности, полученный на основании (4), упрощается и при записи в главных напряжениях имеет вид
(о1 — о3)2 (cos2a + Ksin2a)2 + + [(a1 — o3)(sin2a — Kcos2a) +
+ (c — c )]2 =
v 0 90'J
= [(Oi + G2)K + Co + C90]2, (7) где к = tgp = const; a — угол между a1 и осью x; a1 и o3 — главные напряжения.
Для определения расчетным путем трех характеристик, входящих в (7), достаточно трех результатов одноосного сжатия образца горной породы.
Например, для кварцевого филлита принимаем ст1 = 46,8бМПа (в = 65°), ^ = 45,41МПа (р = 55°), а = 59,48МПа (в1 = 30°).
Подставляя эти результаты в (7), получим систему из трех уравнений, совместным решением которых определяется с0 = 3,52 psi, с90 = 2,01 psi, р90 = = 0,49513. Для определения угла внутреннего трения поперек слоев используются результаты неограниченного сжатия вдоль и поперек слоистости. Угол
внутреннего трения породы поперек слоев определяется усреднением этих результатов. Для данного примера получим р0 = 0,58.
С помощью формул (2) — (4) при различных значениях бокового давления о3 вычисляется разрушающее напряжение Полученные результаты сведены в табл. 3.
Нелинейная геомеханическая модель разрушения
Линейная зависимость между касательным и нормальным напряжениями на площадках возможного сдвига имеет место при малых интенсивностях шарового тензора.
Реальный паспорт прочности для большинства горных пород имеет нелинейную форму, которая должна быть монотонной кривой, симметричной относительно горизонтальной оси, вдоль которой откладываются нормальные напряжения. При больших всесторонних давлениях предельная линия, отделяющая области прочного сопротивления от областей разрушения, приближается к прямой параллельной горизонтальной оси, а с увеличением шарового тензора может иметь замкнутую форму.
Зависимость между о1 и углом в0 при различных значениях о3 Cross-plot of u1 and angle в0 at different values of u3
Ст3, psi Данные в0
0° 30° 45° 60° 90°
0 эксперимент 14,739 8,624 - - 9,512
расчет 11,8 8,15 7 6,75 11,8
725 эксперимент 18,077 9,715 - 8,356 12,119
расчет 14,2 10,3 9,1 8,7 14,2
2176 эксперимент 25,38 18,117 - 14,719 20,364
расчет 19 14,6 13,3 12,9 19
4351 эксперимент 32,627 25,623 - - 28,55
расчет 26,1 21 19,45 19 26,1
7252 эксперимент 39,709 34,064 - 27,008 39,709
расчет 35,8 29,6 27,7 27,2 35,8
10 153 эксперимент 47,732 40,415 - - 47,418
расчет 45,4 38,2 35,9 35,3 45,4
Существующие модификации теории Мора для изотропных горных пород направлены на устранение имеющегося главного недостатка этой теории — неучета влияния среднего главного напряжения. В работах [2, 3] принята параболическая зависимость между касательным и нормальным напряжениями и получены нелинейные паспорта прочности для изотропных горных пород.
Пусть между касательным тп и нормальным стп напряжениями, действующими на площадке сдвига, имеется параболическая зависимость в виде
т2 = а а + Ь, (8)
пп
где а и Ь — прочностные характеристики горной породы.
В работе [2] доказано, что для изотропной горной породы, исходя из (8), после определения параметров а и Ь, можно получить паспорт прочности в следующем виде
1
тП = _ _ стp +
n 2 * с p ' n
а + ст_
(9)
Примем, что входящие в (8) параметры прочности в случае многослойной породы являются функциями угла у и, согласно (9), определяются в виде
а = а(у) = 1/2[аМ — с(у)], (10)
b = b(v) =
М + ^ р М
(11)
где
ас(у) = ас(0) + ас(90) sin2у + + [с(45) — 0,5(с(0) + а{90))]*п2^, (12)
ар(у) = ар(0) + ар (90) sin2у + + [ар(45) —0,5(ар (0) + ар(90))^п2у, (13)
ас(0), ас(90) и ас(45) — пределы прочности горной породы при сжатии поперек и вдоль слоев и под углом 45° к оси х; а (0), а (90) и а (45) — те же самые при отрыве.
Допускается, что неограниченное сжатие образца многослойной породы по любым направлениям соответствует сопротивлению породы сдвиговому разрушению. В условиях однородного нап-
4
ряженного состояния это допущение подтверждается экспериментально, так как при сжатии образца одноосно поперек слоев происходит равномерное уплотнение по всему сечению образца вплоть до появления трещин сдвига, направление которых зависит от степени анизотропии, а при сжатии вдоль слоев имеет место потеря устойчивости всех слоев одновременно. Что касается определения предела прочности при отрыве, существуют схемы испытаний анизотропных пород по «бразильскому методу». Исследования в таком направлении приведены в работах [1, 3, 10].
Для получения выражения критерия прочности в компонентах тензора напряжений предполагаем, что в результате наступления предельного напряженного состояния одна часть многослойной породы скользит относительно другой по некоторой площадке.
Тогда направление площадок скольжения можно определить из условия экстремума
Щ [тп + а а(у) — Ь(у)] = 0. (14)
В дальнейшем для простоты записи формул примем а = а(у), Ь = Ь(у), т = Т, а = а + Г, где а = 0,5(а + а ).
п 1'п 2' " 'чх у'
Рис. 4. График зависимости между с^ и Р при трехосном сжатии Fig. 4. Cross-plot of u1 and в in three-axial compression
Совместное решение (8) и (14) приводит к двум соотношениям, определяющим
т;2 = (а + г2 )а + Ь, (15)
Т3 + р1Т22 + Р2 Т2 + Рз = 0, (16)
где р1, р2, р3 зависят от а, Ь и у.
Для напряжений ах, а , тх будем иметь соотношения
ах = а + (Т1 sm2y + Т2 cos2y)
ау = а-(Т1 + Т2 cos2y) >. (17)
тху = Т1 соэ2у -Т2 sm2y
Возводя в квадрат разность ах — а , тх и суммируя, получим критерий прочности в следующем виде
К)2 + 4x2, = 4К + Т/), (18)
где Т1 определяется по (15), а Т2 по (16).
Направление главных нормальных напряжений находим из условия
fg2a = - 2%ху
о, - ст..
T cos2^ - T2 sin2^ Tj sin2y + T2 cos2^
.(19)
где а — угол между главным напряжением и осью х.
Совместным решением (19) с (15), получим
T2 = z1 -tJzl + (ста + b)ctg2(\\i + a) ,(20) где z1 = 0,5actg2(y + a).
На рис. 4 представлено сравнение экспериментальных данных кварцевого филлита с расчетными значениями.
Испытание образцов пород на неограниченное сжатие и одноосное растяжение (с применением «бразильского теста») не составляет трудности независимо от того, является ли порода анизотропной или нет. Оборудованием для таких испытаний может быть выбран обычный пресс. Следует отметить, что использованные исходные данные, взятые из различных источников, не содержат результатов одноосного растяжения (на разрыв), а в критерий разрушения (18) такие характеристики входят. Учитывая это, с целью анализа пригодности критерия (18) для исследуемых многослойных пород, принимаем параметр п = а /ас = 0,1, что позволяет достаточно точно определить предел прочности исследуемых пород при одноосном растяжении.
Выводы
1. Исследованы геомеханические задачи на основе разработанных линейного и нелинейного критериев разрушения многослойных пород, обладающих анизотропией прочностных свойств.
2. Прочностные характеристики пород, входящие в критерии, предложено определить испытанием на неограниченное сжатие и на одноосное растяжение методом «бразильского теста».
СПИСОК ЛИТЕРАТУРЫ
1. Алиев М. М., Байбурова М. М., Ибрагимов И. И., Исмагилова З. Ф., Лутфуллин А. А. Исследование зависимостей геомеханических характеристик слоистых горных пород от величины гидростатического давления // Society of Petroleum Engineers - SPE Russian Petroleum Technology Conference. - 2017. - SPE-187886-RU.
2. Алиев М. М., Каримова Н. Г., Гилязова С. Р. Нелинейный вариант критерия Кулона-Мора // Горный информационно-аналитический бюллетень. - 2015. - № 7. - С. 226-231.
3. Алиев М. М., Шафиева С. В., Гилязова С. Р., Лутфуллин А. А. Модифицированный нелинейный критерий Мора многослойных пород // Society of Petroleum Engineers - SPE Russian Petroleum Technology Conference. - 2019. - SPE-196905-RU.
4. Ашихмин С. Г., Кашников Ю. А., Шустов Д. В., Кухтинский А. Э. Влияние анизотропии упругих и прочностных свойств пород на устойчивость ствола наклонно направленной скважины // Нефтяное хозяйство. - 2018. - № 2. - С. 54-57.
5. Карев В. И., Коваленко Ю. Ф., Сидорин Ю. В., Устинов К. Б. Геомеханическое моделирование процессов в призабойной зоне скважины // Мониторинг. Наука и технологии. - 2016. - № 3. - С. 85-91.
6. Карев В. И., Климов Д. М., Коваленко Ю. Ф., Устинов К. Б. Модель разрушения анизотропных горных пород при сложном нагружении // Физическая мезомеханика. -2016. - № 6. - С. 34-40.
7. Чанышев A. И. О пластичности анизотропных сред // Прикладная механика и техническая физика. - 1984. - № 2. - С. 149-151.
8. Чанышев A. И., Абдулин И. М. Две характеристические функции поведения образцов сульфидной руды при двухосном сжатии // Физико-технические проблемы разработки полезных ископаемых. - 2018. - № 6. - С. 30-36.
9. Ambrose J. Failure of anisotropic shales under triaxial stress conditions. A thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy and the Diploma of Imperial College, Imperial College London Department of Earth Science and Engineering, June 2014.
10. Yang Z, He B., Xie L, et al. Strength and failure modes of shale based on Brazilian test // Rock and Soil Mechanics, 2015. Vol. 12. Pp. 3447-3456.
11. Lii X., Huang M, Andrade J. E. Strength criterion for cross-anisotropic sand under general stress conditions // Acta Geotechnica. 2016. Vol. 11. No 6. Pp. 1339-1350.
12. Kawa M. Reliability analysis of bearing capacity of square footing on soil with strength anisotropy due to layered microstructure // Studia Geotechnica et Mechanica. 2015. Vol. 37. No 4. Рр. 19-28.
13. Amine S., Latifa O., Lahcen B. Adjusted anisotropic strength model for meta-siltstones and prediction of ucs from indirect tensile tests // International Journal of Civil Engineering and Technology. 2018. Vol. 9. Pp. 598-611.
14. Singh M, Samadhiya N. K., Kumar A., Kumar V. A nonlinear criterion for triaxial strength of inherently anisotropic rocks // Rock Mechanics and Rock Engineering. 2015. Vol. 48. Pp. 1387-1405. EES
REFERENCES
1. Aliev M. M., Bayburova M. M., Ibragimov I. I., Ismagilova Z. F., Lutfullin A. A. Study of geomechanical behavior of stratified rocks as a function of confining pressure. Society of Petroleum Engineers - SPE Russian Petroleum Technology Conference. 2017. SPE-187886-RU.
2. Aliev M. M., Karimova N. G., Gilyazova S. R. Nonlinear variant of the Mohr-Coulomb criterion. MIAB. Mining Inf. Anal. Bull. 2015, no 7, pp. 226-231. [In Russ].
3. Aliev M. M., Shafieva S. V., Gilyazova S. R., Lutfullin A. A. Modified nonlinear Mohr's criterion of multi-layer rocks. Society of Petroleum Engineers - SPE Russian Petroleum Technology Conference. 2019. SPE-196905-RU.
4. Ashikhmin S. G., Kashnikov Yu. A., Shustov D. V., Kukhtinskiy A. E. Effect of anisotropy of elasticity and strength of rocks on inclined borehole stability. Neftyanoe khozyaystvo. 2018, no 2, pp. 54-57. [In Russ].
5. Karev V. I., Kovalenko Yu. F., Sidorin Yu. V., Ustinov K. B. Geomechanical modeling of processes in wellbottom zone. Monitoring. Nauka i tekhnologii. 2016, no 3, pp. 85-91. [In Russ].
6. Karev V. I., Klimov D. M., Kovalenko Yu. F., Ustinov K. B. Fracture model of anisotropic rocks under complex loading. Fizicheskaya mezomekhanika. 2016, no 6, pp. 34-40. [In Russ].
7. Chanyshev A. I. Plasticity of anisotropic media. Prikladnaya mekhanika i tekhnicheskaya fizika. 1984, no 2, pp. 149-151. [In Russ].
8. Chanyshev A. I., Abdulin I. M. Two characteristic functions of behavior of sulphide ore samples under biaxial compression. Fiziko-tekhnicheskiye problemy razrabotki poleznykh iskopayemykh. 2018, no 6, pp. 30-36. [In Russ].
9. Ambrose J. Failure of anisotropic shales under triaxial stress conditions. A thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy and the Diploma of Imperial College, Imperial College London Department of Earth Science and Engineering, June 2014.
10. Yang Z., He B., Xie L., et al. Strength and failure modes of shale based on Brazilian test. Rock and Soil Mechanics, 2015. Vol. 12. Pp. 3447-3456.
11. Lu X., Huang M., Andrade J. E. Strength criterion for cross-anisotropic sand under general stress conditions. Acta Geotechnica. 2016. Vol. 11. No 6. Pp. 1339-1350.
12. Kawa M. Reliability analysis of bearing capacity of square footing on soil with strength anisotropy due to layered microstructure. Studia Geotechnica et Mechanica. 2015. Vol. 37. No 4. Pр. 19-28.
13. Amine S., Latifa O., Lahcen B. Adjusted anisotropic strength model for meta-siltstones and prediction of ucs from indirect tensile tests. International Journal of Civil Engineering and Technology. 2018. Vol. 9. Pp. 598-611.
14. Singh M., Samadhiya N. K., Kumar A., Kumar V. A nonlinear criterion for triaxial strength of inherently anisotropic rocks. Rock Mechanics and Rock Engineering. 2015. Vol. 48. Pp. 1387-1405.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Алиев Мехрали Мирзали оглы1 - д-р техн. наук, профессор, Исмагилова Зульфия Фаритовна1 - канд. техн. наук, доцент, e-mail: [email protected],
Бурмистрова Наталия Николаевна1 - старший преподаватель, e-mail: [email protected],
1 Альметьевский государственный нефтяной институт.
Для контактов: Бурмистрова Н.Н., e-mail: [email protected].
INFORMATION ABOUT THE AUTHORS
M.M. Alievl, Dr. Sci. (Eng.), Professor, Z.F. Ismagilova1, Cand. Sci. (Eng.), Assistant Professor, e-mail: [email protected], N.N. Burmistrova1, Senior Lecturer, e-mail: [email protected], 1 Almetyevst State Oil Institute, 423450, Almetyevsk, Tatarstan, Russia. Corresponding author: N.N. Burmistrova, e-mail: [email protected].
Получена редакцией 20.11.2019; получена после рецензии 06.03.2020; принята к печати 20.07.2020. Received by the editors 20.11.2019; received after the review 06.03.2020; accepted for printing 20.07.2020.