CC BY
28ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2021. Т. 29. № 4
doi: 10.14498/tech.2021.4.3
УДК 622.24
СИСТЕМА ПОДДЕРЖАНИЯ ЭФФЕКТИВНОЙ ОСЕВОЙ НАГРУЗКИ НА ДОЛОТО ПРИ БУРЕНИИ СКВАЖИН
С.А. Михеев
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: serzh.mixeew@yandex.ru
Аннотация. Разрушение горной породы на забое механическим способом невозможно без создания осевой нагрузки на долото. От правильного выбора значения и контроля данного параметра зависит эффективность всего процесса. Контроль осевой нагрузки на породоразрушающий инструмент — один из основных методов управления процессом разрушения горных пород при любом способе бурения. В связи с этим предложена система для автоматизации поиска эффективной осевой нагрузки на долото, позволяющая сократить время нахождения оптимальных значений. В основе работы системы лежит метод заторможенного барабана. Его суть заключается в поиске значений веса инструмента, соответствующих максимуму функции dPoc/dt, т. е. максимуму механической скорости бурения. Представлены блок-схема системы и функциональная схема вычислительного блока. Приведены временные диаграммы колебаний осевой нагрузки на долото и появления сигналов на выходах блоков системы.
Ключевые слова: осевая нагрузка на долото, управление процессом глубокого бурения, повышение точности и надежности измерения.
Введение
В настоящее время бурение 90 % скважин выполняется вручную с использованием отечественных и зарубежных систем мониторинга в соответствии с проектом на строительство. Проектные режимы рассчитываются на основе данных по региону, поэтому не являются оптимальными для конкретной скважины, что может резко снижать (до 25 %) эффективность бурения [1].
При установленных для конкретного долбления типе долота, типе забойного двигателя, расходе промывочной жидкости, параметрах промывочной жидкости единственным регулируемым параметром, с помощью которого может быть осуществлен режим бурения, является нагрузка на долото [2].
Тесная связь осевой нагрузки и механической скорости бурения подтверждена во многих работах [1-9].
Для обеспечения наиболее эффективного объемного способа разрушения горных пород необходимо соблюдать условие
^ > G , (1)
^ шз' V /
* Михеев Сергей Андреевич, аспирант кафедры информационно-измерительной техники 32
где Рос - осевая нагрузка на долото; £к - площадь контакта долота с породой; Ошз - твердость горных пород по штампу в забойных условиях.
С другой стороны, максимальная нагрузка на долото не должна превышать предельно допустимых значений, так как срок службы опор долота уменьшается с увеличением нагрузки.
В условиях бурения часто чередующихся пород с различными физико-механическими свойствами каждой литологической разности должна соответствовать строго определенная осевая нагрузка, отвечающая условию (1).
Однако при неизвестных значениях Ошз и условном определении Рос комплексом наземных установок, особенно в условиях наклонно-направленного бурения, реализовать критерий (1) крайне сложно.
В данных условиях единственным способом нахождения эффективных осевых нагрузок является применение метода заторможенного барабана лебедки (метод отсутствия подачи) [10-15]. Его суть заключается в поиске значений веса инструмента, соответствующих максимуму функции dPoc / dt, т. е. максимуму механической скорости бурения.
В этом случае осевую нагрузку доводят до установленного значения, не превышающего допустимой величины для данного диаметра долота. В ходе углубления долота при заторможенном барабане лебедки вес на крюке возрастает, осевое усилие уменьшается. Тогда механическая скорость может быть найдена по формуле
Т &Р
V = -Т^ (2)
&
где Т0 - длина бурильной колонны без учета деформаций под действием нагрузки; Е0 - модуль упругости материала труб; ^ - эквивалентное сечение труб.
Наибольшая скорость изменения осевой нагрузки &Рос / dt соответствует максимуму механической скорости. Значение осевой нагрузки в этот момент принимается за Р^. На рис. 1 представлен метод заторможенного барабана.
Преимуществом такого метода оптимизации является возможность нахождения эффективной нагрузки на забой в реальных условиях бурения.
Минусом описанного способа является низкая скорость поиска оптимальных значений управляющих параметров (поиск РЭф может достигать 20 минут).
Отсюда вытекают следующие недостатки:
- за время поиска РЭф бурение не выполняется в оптимальном режиме;
- во время поиска может случиться значительный износ долота, и найденная осевая нагрузка больше не будет соответствовать расчетной комбинации «долото - порода», сложившейся на тот момент;
- при поиске продолжительностью 20 минут разбуриваемая порода может смениться и поиск нужно будет произвести еще раз.
На основании вышесказанного разработка системы, позволяющей сократить время поиска эффективной осевой нагрузки на долото, является актуальной.
Разработка системы
На рис. 2 представлена блок-схема системы. Система содержит датчик 1 осевой нагрузки, амплитудный дискриминатор 2, блок 3 вычитания, второй од-новибратор 4, второй пороговый элемент 5, первый пороговый элемент 6, первый ключ 7, интегратор 8, первый одновибратор 9, второй ключ 10, вычислительный блок 11, индикатор 12.
Рис. 2. Блок-схема системы
В зависимости от способа бурения (турбинный, роторный) система реализует различные критерии оценки эффективности нагрузки на долото. Так, в турбинном бурении поиск эффективной нагрузки РЭф производится по критерию
max dP с интегральной оценкой информативного параметра P (t). При этом dt
система интегрирует осевую нагрузку на долото по времени на заданном интервале ее изменения AP (рис. 3), а поиск эффективной нагрузки P^ основан на
сравнении двух соседних интегральных оценок, меньшая из которых расположена ближе к области по шкале осевых нагрузок. Например, сравнение интегральных оценок двух участков кривой P (t) показывает, что
| АР(гуг > | АР(гуг. (3)
Это значит, что для достижения области РЭф необходимо двигаться в сторону уменьшения осевых нагрузок до момента, когда последующая интегральная оценка будет равна предыдущей. При дальнейшем уменьшении нагрузки на долото интегральная оценка начнет расти, последнее будет означать уход из области Рф в область Р < Рэф.
Рис. 3. Диаграмма изменения осевой нагрузки на долото в процессе разбуривания с остановленной подачей
В роторном бурении система реализует модифицированный критерий поиска Рф по максимуму удельной проходки h№ который в общем виде для режима разбуривания с остановленной подачей записывается в виде
К = рЪ (4)
где Уш - механическая скорость бурения; а - эмпирический коэффициент.
Удельная проходка прямо пропорциональна скорости бурения и обратно пропорциональна скорости износа долота, которая является функцией осевой нагрузки.
Модифицированный критерий имеет вид
£ = (5)
данное выражение стремится к минимуму (¿^min) и справедливо при АР = const, здесь АР - величина изменения нагрузки на долото, на которую реагирует второй пороговый элемент; At - время разбуривания с остановленной
35
подачей, за которое осевая нагрузка на долото изменится на АР ; Рт - максимально допустимая нагрузка на долото.
В интеграторе формируется сигнал, пропорциональный величине
АРА =| ЛР(г Щ (6)
2 лг
Величина [Рт-Р(г)] формируется в вычислительном блоке путем подачи на его второй вход сигнала Р (г).
Как показано на рис. 4, изменение осевой нагрузки на долото во времени носит периодический характер.
Рис. 4. Временные диаграммы колебаний осевой нагрузки на долото и появления сигналов на выходах блоков системы
В момент времени г0 меняется направленность изменения осевой нагрузки на долото Р (г) , она начинает уменьшаться.
36
В результате на выходе амплитудного дискриминатора сформируется постоянное напряжение Р0 (см. рис. 4). Появится и начнет расти разность между сигналами на первом и втором входах блока вычитания, на выходе которого формируется сигнал, пропорциональный этой разности.
При достижении упомянутой разности значения ДР1 сработает первый пороговый элемент, который откроет первый ключ и запустит первый одновибра-тор, который, в свою очередь, обнулит содержимое интегратора. После этого с момента времени *[ в интеграторе начнется процесс интегрирования и на его выходе будет формироваться сигнал, пропорциональный площади заштрихованной фигуры над участком кривой Р (t) с координатами (^ ,Р0 ,*2 ,Р2) :
*2
^ =| ДР(* У*. (7)
*1
При достижении разностью сигналов на входе блока вычитания величины ДР2 сработает второй пороговый элемент, который запустит второй одновибра-тор, который, в свою очередь, сбросит на ноль амплитудный дискриминатор и откроет второй ключ. В результате сигнал с выхода интегратора поступит в вычислительный блок, в амплитудном дискриминаторе появится новое значение Р2, а пороговые элементы придут в исходное состояние. Описанный выше процесс повторится сначала, т. к. при достижении разностью сигналов на входах блока вычитания значения ДР1 сработает первый пороговый элемент, который, в свою очередь, обнулит содержимое интегратора. После этого в интеграторе начнется процесс интегрирования и на его выходе будет формироваться сигнал, пропорциональный (см. рис. 4) площади заштрихованной фигуры над участком кривой с координатами (*3,Р2,*4,Р4) :
*4
=| ДР(*^ . (8)
*3
При достижении разностью сигналов на входе блока вычитания величины ДР2 сработает второй пороговый элемент, который запустит второй одновибра-тор, который, в свою очередь, сбросит на ноль амплитудный дискриминатор и откроет второй ключ. Через открытый второй ключ сигнал с выхода интегратора поступит в вычислительный блок.
При роторном бурении программа работы вычислительного блока включает в себя операции деления, в рассмотренном случае
(9)
Р - Р
1 т 1 4
и сравнения частного от деления с ранее полученным результатом деления:
^ (10)
Р - Р
1 т 1 2
Результат сравнения в виде разности Я выводится на индикатор. По мере приближения к области оптимальных нагрузок на долото разность
Я = ---^ (11)
Р - Р Р - Р
1 т 1 2 т 1 4
будет уменьшаться. При прохождении через точку Р^ эта разность сменит знак
на противоположный. Смена знака Я является признаком достижения области эффективных нагрузок на долото.
В соответствии с релейным режимом поддержания осевой нагрузки сигнал о новом нагружении долота может поступить в любой момент времени, не нарушив при этом работу системы. Так как сигнал на выходе блока вычитания из-за нарастания Р (г) станет равным нулю, все блоки системы придут в исходное состояние, минуя такт перезаписи содержимого интегратора в вычислительный блок.
При турбинном способе бурения программа работы вычислительного блока включает операцию сравнения двух последовательных интегральных оценок, например 51 и 52, поступающих на первый вход вычислительного блока. На индикатор выводятся величина и знак разности, которые позволяют судить о том, в какой области осевых нагрузок происходит бурение.
В вычислительном блоке 11 (рис. 5) в первом блоке вычитания 14 реализуется функция Рт - Р (г) , в делителе 15 осуществляется операция деления
Р - Р (г. V
т V . /
(12)
триггер 16 управляет работой коммутатора 17, в блок памяти 18 заносятся значения
Р - Р (г., V
т \ .-1 /
во втором блоке вычитания 19 производится действие
Р - Р (г..) Р - Р (г.)
т \ .-1 / т V . /
(13)
(14)
/г, я
Рис. 5. Функциональная схема вычислительного блока
В зависимости от способа бурения и соответственно от критерия поиска эффективной нагрузки на долото в вычислительном блоке осуществляется операция
^-1 - ^ (15)
либо операция
V V 2' —1___2' (16)
р - р и.. 1) р - р (t,.)'
т \ 21—1/ т \ 21 /
При этом в первом случае блоки 13, 14, 15 в вычислительном процессе не участвуют, а сигнал 1 поступает на входы коммутатора 17 и триггера 16, при этом триггер 16 переводит коммутатор 17 в состояние, когда его вход связан с его первым выходом, и сигнал 52.—1 поступает на вход блока памяти 18.
В результате прихода следующего сигнала триггер 16 переведет коммутатор 17 в положение, когда его вход связан с его вторым выходом, в результате сигнал поступит на второй вход второго блока вычитания 19, на первом входе которого присутствует сигнал S¡—1. Разность Я = —1 — с выхода блока вычитания поступает на вход блока индикации.
Во втором случае в первых трех блоках 13, 14, 15 формируется сигнал
V
-' : сигнал Рт с выхода задатчика установки 13 подается на второй вход
рт р ()
первого блока вычитания 14, на первый вход которого подается сигнал р (7), разность рт — р (t ) с выхода первого блока вычитания поступает на второй вход делителя 15, на выходе которого с частотой поступления на его первый вход сигналов 81 формируется сигнал -'——. На выходы коммутатора 17 подаются
р ( )
сигналы, распределяемые с помощью счетного триггера 16 следующим образом: на первый выход коммутатора 17 подаются сигналы -—-, а на второй
рт р ( —1 )
V.
выход - сигналы
р —Р(t ) '
т.
Разность
Я* = — р~—р(Ц (17)
т . 1 т .
формируется на выходе второго блока вычитания 19 и подается на вход блока индикации 12.
В итоге алгоритм работы системы заключается в следующем:
1. Задается шаг Др .
2. В процессе разбуривания происходит последовательное интегрирование осевой нагрузки с определенным шагом Др .
3. Поиск эффективной нагрузки рЭф основан на сравнении двух соседних
интегральных оценок, меньшая из которых расположена ближе к области по шкале осевых нагрузок.
4. В зависимости от способа бурения на индикатор выводится разность (11) либо разность (17). По мере приближения к области эффективных нагрузок на долото разность будет уменьшаться. При переходе через точку рзф эта разность
изменит знак на противоположный. Смена знака Я является признаком достижения области эффективных нагрузок на долото.
Шаг Др имеет большое значение. Чем он меньше, тем быстрее происходит разбуривание для установления области значений осевых нагрузок, в кото-
рой ведется бурение на данный момент, что дает возможность оперативно принимать решение об уменьшении или увеличении нагрузки в сторону Р^. Это позволяет избежать недостатков метода заторможенного барабана.
Предварительные расчеты зависимости работоспособности системы от шага АР показали, что работоспособность сохраняется даже при наименьшем шаге изменения осевой нагрузки АР (АР = 0,25 т; 0,5 т), то есть меньше времени потребуется на разбуривание и нахождение области эффективных значений.
Однако нужно учитывать погрешность измерения веса бурильного инструмента. Существующие гидравлические и электрические измерительные устройства позволяют определять значение веса бурильного инструмента Q и осевой нагрузки на долото Р в лучшем случае с погрешностью 2,5 %, причем погрешность измерения увеличивается на 0,5 % при изменении окружающей температуры на каждые 10 °С. Например, при измерении гидравлическим индикатором веса (ГИВ-6) с трансформатором давления с максимальным усилием натяжения каната по паспорту 320 кН абсолютная погрешность будет составлять 0,8 т.
Так как осевая нагрузка на долото Р создается частью веса колонны бурильных труб и определяется как разность веса бурильного инструмента Q0 до начала бурения, когда долото не касается забоя скважины, и веса бурильного инструмента на крюке Q¡ в процессе бурения Р = Q0 - Q1, то эффективность системы упирается в точность определения веса бурильного инструмента. Исходя из этого разработка методов повышения точности устройств измерения веса является предметом для дальнейших исследований.
Выводы
Разработанная система позволяет сократить время поиска эффективной осевой нагрузки на долото и избавиться от недостатков метода заторможенного барабана.
В системе применяется операция интегрирования, что позволяет повысить помехоустойчивость и точность всей измерительной схемы, так как операция интегрирования обладает фильтрующим свойством, освобождая измерительный сигнал от высокочастотных помех.
Также система может применяться для различных способов бурения (турбинный, роторный) и в зависимости от этого реализовывать различные критерии оценки эффективности нагрузки на долото.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Цуприков А.А. Интеллектуальная система адаптивного управления технологическим процессом бурения нефтегазовых скважин: Дисс. ... д-ра техн. наук. - Краснодар: Кубанский гос. техн. ун-т, 2018. - 310 с.
2. Лукьянов Э.Е. Геолого-технологические и геофизические исследования в процессе бурения. - Новосибирск: Историческое наследие Сибири, 2009. - 752 с.
3. Борзенко В.В., Некрасов С.А. Оптимизация процесса бурения на основе полного факторного эксперимента // NOVAINFO.RU. - 2016. - № 50. - Т. 2. - С. 1-7. ШЬ: https://novainfo.ru/article/7613 (дата обращения: 20.06.2021).
4. Бражников В.А., Фурнэ А.А. Информационное обеспечение оптимального управления бурением скважин. - М.: Недра, 1989. - 206 с.
5. Бродов Г.С., Шелковников И.Г., Егоров Э.К. Технологические измерения и автоматизация процесса бурения. - СПб.: Санкт-Петербургский горный институт, 2004. -105 с.
6. Порцевский А.К., Ганджумян Р.А. Оптимизация буровых и горно-разведочных работ,
планирование эксперимента: Учеб. пособие. - М.: МГОУ, 2005. - 70 с.
7. Ситников Н.Б. Моделирование и оптимизация процесса бурения геологоразведочных скважин: Дисс. ... д-ра техн. наук. - Екатеринбург, 2000. - 350 с.
8. Цуприков А.А. Анализ математических моделей механической скорости проходки для оптимизации процесса бурения нефтегазовых скважин // Научный журнал КубГАУ. -2015. - № 107(03). - С. 903-916.
9. Цуприков А.А. Математическая модель скорости проходки для оптимального управления бурением скважин // Научный журнал КубГАУ. - 2015. - № 107(03). - С. 916-926.
10. Асадчев А.С. Разработка технологии бурения глубоких скважин гидравлическими забойными двигателями в условиях соленосных отложений: Дисс. . канд. техн. наук. -М.: Российский гос. ун-т нефти и газа им. И.М. Губкина, 2012. - 253 с.
11. Грибанов М.В. Совершенствование и реализация методик и программных средств по оптимизации параметров бурения скважин: Дисс. . канд. техн. наук. - Ухта: Ухтинский гос. техн. ун-т, 2012. - 116 с.
12. Грибанов М.В. Перспективы создания единой базы данных по буримости пород по группе месторождений с целью проектирования и оптимизации режимов бурения // Строительство нефтяных скважин на суше и на море. - 2009. - № 1. - С. 15-18.
13. Дробышевский В.Д., Асадчев А.С. Определение оптимального режима бурения при инженерно-технологическом сопровождении промысловых испытаний новой техники // Современные проблемы машиноведения: материалы XII Международной научно-технической конференции. - Беларусь, Гомельский гос. техн. ун-т имени П.О. Сухого, 2018. - С. 307-309.
14. Струговец Е.Т., ЛугумановМ.Г. О поиске оптимальных нагрузок на долото при бурении забойными двигателями // НТВ Каротажник. - Тверь, 2005. - Вып. 5-6 (132-133). -С. 107-113.
15. Патент России № 2642699, МПК E21B 44/02 (2006.01). Способ регулирования условий процесса бурения скважин / Александров С.С.; патентообладатель: Федеральное государственное бюджетное образовательное учреждение высшего образования «Уфимский государственный нефтяной технический университет». Заявл. 27.02.2017; опубл. 25.01.2018, Бюл. № 3.
Статья поступила в редакцию 24 июня 2021 г.
SYSTEM FOR MAINTAINING EFFECTIVE AXIAL LOAD ON THE BIT WHEN DRILLING WELLS
S.A. Mikheev*
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation E-mail: serzh.mixeew@yandex.ru
Abstract. The destruction of rock at the bottom by mechanical means is impossible without creating an axial load on the chisel. The efficiency of the entire process depends on the correct choice of the value and control of this parameter. Control of the axial load on the rock-breaking tool is one of the main methods of controlling the process of rock destruction in any drilling method. In this regard, a device for automating the search for an effective axial load on the bit is proposed, which allows reducing the time to find the optimal values. The operation of the device is based on the method of a braked drum. Its essence is to find the values of the tool weight corresponding to the maximum of the dPoc/dt function, i.e. the maximum of the mechanical drilling speed. The block diagram of the device and the functional diagram of the computing unit are presented. The time diagrams offluctuations of the axial load on the chisel and the appearance of signals at the outputs of the device blocks are given.
Keywords: axial load on the bit, control of the deep drilling process, improving the accuracy and reliabili/y of measurement.
REFERENCES
1. Tsu^rikov A.A. Intellektual'naya sistema adaptivnogo upravleniya tekhnologicheskim processom bureniya neftegazovyh skvazhin: Diss. ... dokt. tekhn. nauk. - Krasnodar: Kubanskij gos. tekhn. un-t, 2018. - 310 p. (In Russian).
2. Luk'yanov E.E. Geologo-tekhnologicheskie i geofizicheskie issledovaniya v processe bureniya [Geological, technological and geophysical studies in the drilling process]. - Novosibirsk: Iz-datel'skij Dom «Istoricheskoe nasledie Sibiri», 2009. - 752 p. (In Russian).
3. Borzenko V.V., Nekrasov S.A. Optimization of the drilling process based on a complete factor experiment // NOVAINFO.RU. - 2016. - № 50. - T. 2. - Pp. 1-7. https://novainfo.ru/article/7613 (accessed June 20, 2021).
4. Brazhnikov V.A., Furne A.A. Informacionnoe obespechenie optimal'nogo upravleniya bureniem skvazhin [Information support for optimal well drilling management]. - Moscow: Nedra, 1989. - 206 p. (In Russian).
5. Brodov G.S., Shelkovnkov I.G., Egorov E.K. Tekhnologicheskie izmereniya i avtomatizaciya processa bureniya [Technological measurements and automation of the drilling process]. -Sankt-Petersburg. Sankt-Peterburgskij gornyj institut, 2004. - 105 p. (In Russian).
6. Porcevskj A.K., Gandzhumyan R.A. Optimizaciya burovyh i gorno-razvedochnyh rabot, plani-rovanie eksperimenta [Optimization of drilling and exploration operations, planning of the experiment.]. Uchebnoe posobie. MGOU. - Moscow, 2005. - 70 p. (In Russian).
7. Sitnikov N.B. Modelirovanie i optimizaciya processa bureniya geologorazvedochnyh skvazhin: Diss. ... dokt. tekhn. nauk. - Ekaterinburg, 2000. - 350 p. (In Russian).
8. Tsu^rikov A.A. Analysis of mathematical models of mechanical penetration rate to optimize the drilling process of oil and gas wells // Nauchnyj zhurnal KubGAU. - 2015. - № 107(03). - Pp. 903-916. (In Russian).
9. Tsu^rikov A.A. Mathematical model of the penetration rate for optimal control of well drilling // Nauchnyj zhurnal KubGAU. - 2015. - № 107 (03). - Pp. 916-926. (In Russian).
10. Asadchev A.S. Razrabotka tekhnologii bureniya glubokih skvazhin gidravlicheskimi zabojnymi dvigatelyami v usloviyah solenosnyh otlozhenij: Diss. ... kand. tekhn. nauk. - Moscow: Rossi-
Sergey A. Mikheev, Postgraduate Student 42
jskij gos. un-t nefti i gaza im. I.M. Gubkina, 2012. - 253 p. (In Russian).
11. Gribanov M.V. Sovershenstvovanie i realizaciya metodik i programmnyh sredstv po optimizacii parametrov bureniya skvazhin: Diss. ... kand. tekhn. nauk. - Uhta: Uhtinskij gos. tekhn. un-t, 2012. - 116 p. (In Russian).
12. Gribanov M.V. Prospects for creating a unified database on the drillability of rocks for a group of deposits in order to design and optimize drilling modes // Stroitel'stvo neftyanyh skvazhin na sushe i na more. - 2009. - № 1. - Pp. 15-18. (In Russian).
13. Drobyshevskij V.D., Asadchev A.S. Determination of the optimal drilling mode during engineering and technological support of field tests of new equipment // Sovremennye problemy mashi-novedeniya: materialy XII Mezhdunarodnoj nauchno-tekhnicheskoj konferencii. - Belarus', Gomel'skij gos. tekhn. un-t imeni P.O. Suhogo. - 2018. - Pp. 307-309. (In Russian).
14. Strugovec E.T., LugumanovM.G. On the search for optimal loads on the bit when drilling with downhole motors // NTV Karotazhnik. - Tver', 2005. - Vol. 5-6 (132-133). - Pp. 107-113. (In Russian).
15. Patent of Russia 2642699, MPK E21B 44/02 (2006.01). A method for regulating the conditions of the well drilling process / Aleksandrov S.S.; patentoobladatel': Federal'noe gosudarstvennoe byudzhetnoe obrazovatel'noe uchrezhdenie vysshego obrazovaniya "Ufimskij gosudarstvennyj neftyanoj tekhnicheskij universitet". - Zayavl. 27.02.2017; Publ. 25.01.2018, Buul. № 3 (In Russian).
