Инициирование ЭВВ когерентным промежуточным детонатором Текст научной статьи по специальности «Физика»

УДК 662.235

С.А. Горинов, И.Ю. Маслов

ИНИЦИИРОВАНИЕ ЭВВ КОГЕРЕНТНЫМ ПРОМЕЖУТОЧНЫМ ДЕТОНАТОРОМ

Предложена методика расчета параметров инициирования зарядов эмульсионных ВВ когерентным промежуточным детонатором (системой сближенных сосредоточенных зарядов). В методике учитываются размеры зарядов ПД, расстояния между зарядами ПД, детонационные параметры ВВ, из которого изготовлен ПД, а также химический состав, плотность и структура возбуждаемых ЭВВ, а также физико-механические свойства окружающих горных пород.

Ключевые слова: промежуточный детонатор, эмульсионные ВВ.

Актуальность работы

Практика ведения взрывных работ на карьерах с использованием наливных ЭВВ, сенсибилизированных газовыми порами, выявляет тенденцию к применению все более мощных промежуточных детонаторов (ПД), как линейных (удлиненных) ПД [1—4], так и мощных сосредоточенных ПД [4, 5—8]. Это объясняется стремлением избежать возникновения режимов низкоскоростного взрывчатого разложения ЭВВ, при которых не только снижается эффективность и безопасность взрывных работ, но и ухудшается экология вследствие повышенного выброса ядовитых газов в атмосферу.

Однако в настоящее время промышленностью освоен выпуск большой номенклатуры сосредоточенных промежуточных детонаторов (в виде «шашек») небольшой и средней мощности, успешно зарекомендовавших себя при инициировании простейших и гранулированных ВВ. Многие горные предприятия уже давно используют две (и более) такие шашки для инициирования наливных ЭВВ, сенсибилизированных газовыми порами. Однако отсутствие методики расчета параметров ПД, для инициирования зарядов ВВ указанным способом, затрудняет применение данного способа инициирования ЭВВ.

В настоящей работе представлена методика расчета, позволяющая определять размеры зарядов ПД, оптимальные расстоя-

ISSN 0236-1493. Горный информационно-аналитический бюллетень. 2017. № 4. С. 293-304. © 2017. С.А. Горинов, И.Ю. Маслов.

ния между зарядами ПД в зависимости от детонационных параметров ВВ, из которого изготовлен ПД, а также химического состава, плотности и структуры возбуждаемых ЭВВ. Это позволяет считать данную работу актуальной.

Материал и результаты исследований

Эффективность инициирующего воздействия на ЭВВ сближенных сосредоточенных ПД вытекает из эффекта взаимодействия встречных ударных волн при взрыве сближенных зарядов. В этом случае (при равенстве данных сближенных зарядов) отражение ударных волн друг от друга происходит, как от абсолютно жесткой стенки. Это способствует увеличению коэффициента отражения в (1 + ро") раз по сравнению с отраже-

рсСр

нием от стенки скважины. Данная система зарядов и будет тем самым ПД, который в дальнейшем будем ее называть «когерентным ПД».

При разработке расчетной схемы инициирования ЭВВ сближенными сосредоточенными промежуточными детонаторами, исходим из следующих данных:

1. Заряды, составляющие когерентный ПД, моделируются сферами радиуса Яб. Расстояние между центрами сближенных зарядов ПД принимается равным h ^ > 2Я&).

2. ВВ, используемое в ПД, характеризуется следующими параметрами: рб, Dб, kб — плотность, скорость детонации и коэффициент политропы ВВ ПД, соответственно.

3. ЭВВ в инициируемом заряде характеризуется: роо — плотностью, Яо — радиусом поры, А — размером частицы эмульсии и ударной адиабатой Dу = А + ВЩ, где А, В — параметры ударной адиабаты инициируемого ВВ — скорость ударной волны, Щу — массовая скорость).

4. Возбуждение процесса детонации в ЭВВ не является мгновенным. Для возбуждения данного процесса, вещество ЭВВ необходимо сжать, обеспечить его возгорание и прогорание между «горячими» точками. Данный процесс происходит за определенное время т*, которое определяется особенностями инициируемого ЭВВ и величиной нагружающего импульса [9].

При взрыве сферического заряда (ПД) в ЭВВ возникает расширяющийся в радиальном направлении газовый пузырь, который при расширении испускает ударную волну, также распространяющуюся в радиальном направлении.

Введем обозначения: Яп — текущий радиус газового пузыря; Яу — текущий радиус фронта ударной волны.

Согласно [10] время расширения газового пузыря до размеров сферы радиуса Яп определяется по формуле

где

Ь_ = „\Р^ Г ¿Я

н = я2„ф( Я) ¿яп

описывается уравнением

ф( Я ) = Р

г я ^

3ке

V Я у

где Рн — начальное детонационное давление ПД. На основании (1)—(3) имеем

^ = Я

3 (К -1)

I 2 Р

1-

где г =

Я Я

г,

1

■¿г„

(1) (2)

(3)

(4)

Точного решения в элементарных функциях интеграл, входящий в (4), не имеет. Однако при Яп/ Яб < 1,7 уравнение (4) можно с точностью до 0,1% представить в виде (учитываем, что

р =

ре вв )

К +1)

к

Я А

3 (К2 -1)

. . 2 рб (5)

•((0,1588гп2 - 0,5848гп + 0,4066) ке + 2,2226гп - 2,0937)

Вследствие резкого падения давления во взрывных газах ПД при их расширении, практический интерес представляет область Яп/Яб < 1,7, т.к. при Яп/Яб > 1,7 давление взрывных газах ПД упадет в >1,79 « 120 раз. Согласно [11] давление во фронте детонационной волны в составе В (60% гексогена и 40% тротила) при плотности ВВ 1,67—1,68 г/см3 (данный состав «ТГФ» является одним из самых мощных при производстве промышленных ПД) составляет 24,3—25,6 ГПа.

При падении давления в 120 раз его остаточная величина составляет примерно 0,2 ГПа (практически это граница области

К

е

е

К

г

расширения продуктов взрыва с коэффициентом политропы kб [12]). Данная область давлений не представляет интереса при рассмотрении вопроса о непосредственном возбуждении высокоскоростной детонации при воздействии ПД на ЭВВ (в соответствии с критериями [13, 14]).

Так как сжимаемость продуктов детонации ПД ниже, чем сжимаемость пористого вещества ЭВВ, то происходит истечение продуктов детонации ПД в ЭВВ. Начальная скорость смещения границы «продукты взрыва ПД-ЭВВ» Wm можно найти на основании равенства давлений в ЭВВ и взрывных газов ПД на границе «продукты взрыва ПД-ЭВВ» [15]:

Р00 (А + ВШт)Wro = р

( ( 1 -

ш -

в.

кб +1

Ч -1

2кб Лкб-1

2кбОб,

(6)

Величину массовой скорости за фронтом ударной волны в ближней зоне определяем по уравнению [16]

ш К.

ГО Я.

(7)

Тогда скорость ударной волны на расстоянии R будет равна

в.,

А + ВШГ0 К . Я.

(8)

йЯ,

Учитывая, что В = —- , интегрируя (8), получаем уравне-

- йЬ

ние движения фронта ударной волны в ближней зоне заряда

Я,

Яб + АЬ +

ВШ

_ГО

А

Яб

(К + вш ^

Я А

1+

ВШ

А

(9)

На рис. 1 представлена схема взаимодействия двух сближенных зарядов. Принятые обозначения: сплошные линии — начальная граница зарядов; штрихпунктирные линии — граница «продукты взрыва ПД-ЭВВ»; штриховые линии — граница фронта ударной волны; О1 и О2 — центры зарядов; ABCD — область взаимодействия ударных волн (область повышенного давления).

Расчет параметров когерентного ПД осуществляется из требования выполнения следующего условия [13, 14]: для возбуж-

Рис. 1. Схема взаимодействия сближенных зарядов когерентного ПД

дения детонации с величиной фронтального давления Pf необходимо в слое возбуждаемого ЭВВ, размером не менее 4Ы(Р), в течении времени поддерживать давление с наименьшим значением на указанном временном промежутке — (4Ы1Р), — критический диаметр заряда и время химической реакции при фронтальном давлении Р^.

Принимаем за момент t = 0 момент выхода детонационной волны ПД на внешнюю поверхность ПД.

Допустим, что инициирование ЭВВ вследствие взаимодействия ударных волн произошло в момент времени tn, когда радиус газового пузыря достигнет размера Rn(tn) = Радиус фронта ударной волны при этом определяется на основании (5) и (9). Величину этого радиуса можно представить в виде функции R = R и ).

1 у у4 п7

Из рис. 1 имеем следующие ограничения на R (t) и R (:):

К (Ьп) + Яу (ьп )< Н , Яу (Ьп )> н

На основании (9) и рис. 1 имеем: А) время встречи ударных волн Ь

(1)

¿(1) =

2ЯЛ

-1 - ВЩ, 1п

А

Г _н_

2ЯЛ

ВЩ А

2

\\

(10)

1 +

ВЩГ0 А

А (Ц)

//

Б) время выхода ударной волны на поверхность газового пу-

зыря смежного заряда Ь

(2)

Л2) _

1У ~

Ь - К к

BW

го

А

гн - К + BWr,

1п

К

А

1 +

BWro А

К (12)

А

; )

При времени I > Ь,У2) происходит инжекция эмульсии в газо-

у

вые полости смежных зарядов, что сопровождается резким падением давления в эмульсии. Поэтому, если ЭВВ не будет ини-

циировано к моменту

t = Ь,(2)

то инициирования не произойдет

вовсе.

В соответствии с (10) имеем следующее ограничение на характерные моменты времени для когерентного ПД

4 < и,

(2)

(13)

Время существования области высокого давления АBCD (рис. 1) равно (1)

Ь.= - ^ (14)

Наименьшее давление в области высокого давления соответствует времени tn. В соответствии с закономерностью падения давления при удалении от сферического заряда [12] имеем

Р (К) * рЛ

яй

(

Я, (Ьп )

А + BW„„

я

V

н

У

Я (Ь ) 2Я (Ь )

У \ п ) /V У \ п/

£ж.с.

отр

(15)

где к^р. * 2 (1 + фввхоо) — коэффициент отражения ударной волны в ЭВВ от жесткой стенки.

Объем области высокого давления в момент tn (сечение указанного объема представлено на рис. 1 фигурой ЛВСВ) равен (необходимые сведения для определения объема указанной области заимствованы в руководстве [16])

^ (о=2п(*. (о-н/2 |Ч (о+ну. (1б)

Диаметр сферы соответствующего объема (16), равен

Л & )= 4 [Я, (Ья)-Н ] (2Яу (Ья) + Н

1

(17)

На основании уравнений (5)—(17) определяются все необходимые для дальнейшего анализа параметры t , Р^ ), d(t), харак-

теризующие воздействие когерентного ПД из двух сближенных зарядов на ЭВВ.

Если заряды когерентного ПД прижаты к стенке скважины,

I

то при условии JR2 (tn)--> R6, необходимо учитывать повторное отражение эмульсии, находящейся в области ABCD (рис. 1), от стенок скважины. Это приведет к возрастанию величины P(t)

2Pccp ( C й П

в ——-— раз (pc, Cp — плотность окружающей заряд среды

PcCp + PooDy

и скорость продольных волн в этой среде, соответственно).

Окончательное решение возможно после получения зависимостей Pf, dkrit(Pf), t(Pf) для инициируемого ЭВВ. Данные зависимости определяются для инициируемого ЭВВ с учетом его химического состава, плотности и структуры на основании методики, изложенной в работе [9].

Перейдем к обсуждению полученных результатов.

Для проверки выявленных закономерностей распространения детонационных волн в ЭВВ, сенсибилизированных газовыми порами, были осуществлены опытные взрывы, в которых скважинные заряды ЭВВ инициировались когерентными ПД. Ниже приведен характерный результат из данной серии опытов, выполненный на объекте открытых горных работ АО «Уралас-бест». В опытах измерялись скорость детонации по длине заряда, плотность ЭВВ при атмосферном давлении, глубина скважины, высота колонки заряда, глубина расположения ПД.

Применялось ЭВВ «Порэмит-1А», имеющее следующие характеристики — плотность — 1,0 г/см3, хим.состав: аммиачная селитра — 73,8%, вода — 18,5%, топливная фаза — 7,7%; способ сенсибилизации — газовые поры. Скважины бурились в прочных горных породах (коэффициент крепости по М.М. Протодъя-конову f = 13). Глубина скважин — 17 м, длина заряда — 12 м, наклон скважин — 75°, диаметр — 250 мм. Заряды ЭВВ формировались при зарядке из одной зарядной машины. Измерения скорости детонации проводились реостатными приборами «Han-ditrap Vod Recorder» — датчиковый провод «VOD probecable-HT» устанавливался по длине заряда.

Инициирование зарядов (нижнее) производилось от промежуточных детонаторов, состоящих из шашек «ПТ-П750», расположенных одна над другой на расстоянии между ними 15 см, и внутрискважинных капсюлей неэлектрической системы Искра

Рис. 2. Зависимость длины (м) заряда от времени детонации (мс) ЭВВ. Промежуточный детонатор, состоящий из 2 шашек «ПТ-П750»

С-500-18. Результаты обрабатывались в программе «Нап&ТгарП-у2р0», поставляемой изготовителем с прибором, и представлены на рис. 2.

На рис. 3 приведены графики экспериментально определенной и расчетной скоростей детонации по длине скважинного заряда. Расчет скоростей детонации осуществлялся по методике, изложенной в [9], а условий скачкообразного перехода на более низкоскоростной режим — по методике, изложенной в [17]. При выполнении расчетов исходили из приведенного выше химиче-

Рис. 3. Зависимости скорости детонации ЭВВ от места инициирования (экспериментальные и расчетные)

Расчетные значения параметров детонации в ЭВВ

при их инициировании когерентным ПД г (по длине заряда)

Расстояние от ПД, м Текущая плотность ЭВВ, г/см3 Текущий коэфф. политропы ЭВВ Скорость детонации, м/с Давление в т. Чеп-мена-Жуге, ГПа Время про-горан. между порами, мкс

0,0 1,198 3,21 5661,2 8,68 5,6

5,64 (момент перехода на НСР) 1,137 3,06 5657,0 8,26 3,6

ского состава ЭВВ, его плотности при атмосферном давлении и структуры (размер эмульсии — 10 мкм, размер газовых пор — 300 мкм). Изменение плотности ЭВВ по длине заряда рассчитывали по методике [18].

Расчеты показали, что возбуждение высокоскоростного режима детонации происходит в результате 3-кратного ударного сжатия ЭВВ (сжатие в прямой ударной волне УВ — сжатие при взаимодействии УВ от сближенных зарядов ПД — сжатие в УВ, отраженной от стенки скважины). Результаты расчетов приведены в таблице.

Анализ рис. 3 показывает, что разработанная авторами методика расчета позволяет достаточно точностью определять параметры детонации, вызванной воздействием когерентного ПД. Кроме этого, расчетный и экспериментально определенный моменты перехода детонации на низкоскоростной режим (НСР) близки между собой (рис. 3). Это позволяет принять необходимые технологические меры для обеспечения полноты детонации (например, изменение места установки ПД, установка дополнительных ПД, изменение начальной плотности ЭВВ, изменение структуры ЭВВ и т.д.).

Выводы

В работе представлена методика расчета параметров системы сближенных сосредоточенных зарядов, представляющей промежуточный детонатор (ПД) для инициирования зарядов эмульсионных ВВ.

В методике учитываются размеры зарядов ПД, расстояния между зарядами ПД, детонационные параметры ВВ, из которого изготовлен ПД, химический состав, плотность и структура возбуждаемых ЭВВ, а также физико-механические свойства окружающих горных пород.

Наблюдается хорошее согласие с экспериментальными данными.

Работа представляет практический интерес для организаций выпускающих и разрабатывающих ЭВВ и ПД, а также для инженеров руководящих взрывными работами на горнодобывающих предприятиях.

СПИСОК ЛИТЕРАТУРЫ

1. Калякин С. А., Прокопенко В. С. О линейном инициировании зарядов взрывчатых веществ // Сучасш ресурсоенергозберiгаючi технологи прничного виробництва. — 2013. — 1(11). — С. 41—47.

2. Тогунов М. Б., Сапронов Е. М., Щукин Ю. Г. и др. Повышение эффективности взрывания горных пород эмульсионными ВВ / Технология и безопасность взрывных работ. — Екатеринбург: ИГД УрО РАН, 2012. - С. 123-133.

3. Щукин Ю. Г., Чернышов С. Н., Коломинов И. А. и др. Перспективные взрывчатые материалы для горной промышленности / Развитие ресурсосберегающих технологий во взрывном деле. - Екатеринбург, 2009. - С. 108-118.

4. Фокин В. А. Обоснование геометрических параметров удлиненной шашки-детонатора для обратного инициирования скважинных зарядов газифицированных эмульсионных взрывчатых веществ // Известия вузов. Горный журнал. - 2008. - № 3. - С. 49-54.

5. Маслов И.Ю., Пупков В.В., Фоменкова В.Е. и др. Улучшение качества взрывной подготовки горной массы за счет применения промежуточных детонаторов с оптимальными габаритными размерами при инициировании скважинных зарядов эмульсионных ВВ // Взрывное дело. - 2004. - № 94/51. - С. 125-130.

6. Добрынин И. А. Обоснование параметров промежуточных детонаторов в скважинных зарядах для повышения эффективности дробления горных пород. Автореф. дисс. на соиск. науч. степени канд. техн. наук. - М., 2010. - 20 с.

7. Добрынин А. А. Взрывчатые вещества. Химия, составы. Безопасность. - М.: изд. ИД Академии Жуковского, 2014. - 527 с.

8. Калякин С. А. Обоснование параметров промежуточных детонаторов для инициирования скважинных зарядов ВВ // Сучасш ресурсоенергозберiгаючi технологи прничного виробництва. - 2012. -1(9). - С. 50-57.

9. Кутузов Б. Н., Горинов С. А. Физико-технические основы создания эмульсионных и гранулированных ВВ и средств их инициирования // Горный информационно-аналитический бюллетень. Отдельные статьи (специальный выпуск). Эмульсионные ВВ, гранэмиты и ANFO: структура, инициирование, физико-технические основы создания. -2011. - № 7. - С. 34-52.

10. Власов О. Е. Основы теории действия взрыва. - М.: изд. ВИА, 1957. - 408 с.

11. Кук М. А. Наука о промышленных ВВ. - М.: Недра, 1980. -453 с.

12. Баум Ф. А., Станюкович К.П., Шехтер Б.И. Физика взрыва. -М.: Физматгиз, 1959. - 800 с.

13. Баум Ф. А., Державец А. С. Об инициировании детонации ударными волнами // Взрывное дело. — 1966. — № 60/17. — C. 68-83.

14. Афанасенков А. Н., Богомолов В. М., Воскобойников И. М. Критические давления инициирования взрывчатых веществ // Взрывное дело. — 1970. — № 68/25. — C. 68—92.

15. Орленко Л. П. Физика взрыва и удара. 2-е изд., испр. — М.: Физ-матлит, 2008. — 304 с.

16. Бронштейн И. Н., Семендяев К. А. Справочник по математике для инженеров и учащихся втузов. — М.: Глав.ред.физматлит, 1980. — 976.

17. Горинов С. А., Кутузов Б.Н. О неустойчивости детонационных волн в эмульсионном взрывчатом веществе, сенсибилизированном газовыми порами // Горный информационно-аналитический бюллетень. — 2012. — № 4. — C. 302—307.

18. Фокин В. А. Распределение плотности эмульсионных взрывчатых веществ по высоте колонки скважинного заряда // Известия вузов. Горный журнал. — 2007. — № 3. — C. 89—94. [¡233

КОРОТКО ОБ АВТОРАХ

Горинов Сергей Александрович1 — кандидат технических наук, e-mail: akaz2006@yandex.ru,

Маслов Илья Юрьевич1 — кандидат технических наук, e-mail: ilmaslov@mail.ru, 1 ООО «Глобал Майнинг Эксплозив-Раша».

UDC 662.235

Gornyy informatsionno-analiticheskiy byulleten'. 2017. No. 4, pp. 293-304.

S.A. Gorinov, I.Yu. Maslov

EE INITIATION BY COHERENT BOOSTER

Actual blasting of pourable gas-synthesized emulsion explosives (EE) calls for highly powerful relay detonators. However, under production is currently a wide range of low and medium capacity explosive cartridges successfully used to prime simplest and granular explosives.

In order to enhance efficiency of initiation, it is recommended to use a relay detonator constituted by closely spaced and concurrently fired explosive cartridges. The efficiency of such a booster is achieved due to interaction of meeting shock waves generated under explosion of closely spaced charges. In this case (given equality of the closely spaced charges), the shock waves are reflected as from absolute rigid wall.

The article proposes the calculation procedure for the parameters of initiation of emulsion explosives by a coherent booster (a system of closely spaced charges of concentrated explosives). The procedure takes into account sizes of booster charges, their spacing, proknock properties of explosives used to make the booster, as well as chemical composition, density and structure of the fired explosives, and physical properties of surrounding rocks.

The calculation of the coherent booster parameters based on the Baum-Derzhavets criterion (detonation at the value of the front pressure needs the least pressure to be maintained in the EE layer with the size not less than critical within the time sufficient for the chemical response execution) accounts for the shock adiabat of the initiatied emulsion explosive.

Key words: booster, emulsion explosives.

AUTHORS

Gorinov S.A.1, Candidate of Technical Sciences, e-mail: akaz2006@yandex.ru, Maslov I.Yu.1, Candidate of Technical Sciences, e-mail: ilmaslov@mail.ru,

1 Global Mining Explosive Group, 117036, Moscow, Russia.

REFERENCES

1. Kalyakin S. A., Prokopenko V. S. Suchasni resursoenergozberigayuchi tekhnologiigir-nichnogo virobnitstva. 2013. 1(11), pp. 41—47.

2. Togunov M. B., Sapronov E. M., Shchukin Yu. G. Tekhnologiya i bezopasnost' vz-ryvnykh rabot (Technology and safety of blasting), Ekaterinburg, IGD UrO RAN, 2012, pp. 123-133.

3. Shchukin Yu. G., Chernyshov S. N., Kolominov I. A. Razvitie resursosberegayush-chikh tekhnologiy vo vzryvnom dele (Development of resource-saving technologies in the explosive case), Ekaterinburg, 2009, pp. 108-118.

4. Fokin V. A. Izvestiya vuzov. Gornyy zhurnal. 2008, no 3, pp. 49-54.

5. Maslov I. Yu., Pupkov V. V., Fomenkova V. E. Vzryvnoe delo. 2004, no 94/51, pp. 125-130.

6. Dobrynin I. A. Obosnovanie parametrov promezhutochnykh detonatorov v skvazhin-nykh zaryadakh dlya povysheniya effektivnosti drobleniya gornykh porod (Justification of parameters of intermediate detonators in the borehole charges to improve the efficiency of rock crushing), Candidate's thesis, Moscow, 2010, 20 p.

7. Dobrynin A. A. Vzryvchatye veshchestva. Khimiya, sostavy. Bezopasnost' (Explosives. Chemistry, structures. Security), Moscow, izd. ID Akademii Zhukovskogo, 2014, 527 p.

8. Kalyakin S. A. Suchasni resursoenergozberigayuchi tekhnologii girnichnogo virobnitstva. 1/2012 (9), pp. 50-57.

9. Kutuzov B. N., Gorinov S. A. Gornyy informatsionno-analiticheskiy byulleten'. Special edition. 2011, no 7, pp. 34-52.

10. Vlasov O. E. Osnovy teorii deystviya vzryva (Fundamentals of the theory of explosion), Moscow, izd. VIA, 1957, 408 p.

11. Kuk M. A. Nauka opromyshlennykh VV(Science of industrial explosives), Moscow, Nedra, 1980, 453 p.

12. Baum F. A., Stanyukovich K. P., Shekhter B. I. Fizika vzryva (Physics of explosion), Moscow, Fizmatgiz, 1959, 800 p.

13. Baum F. A., Derzhavets A. S. Vzryvnoe delo. 1966, no 60/17, pp. 68-83.

14. Afanasenkov A. N., Bogomolov V. M., Voskoboynikov I. M. Vzryvnoe delo. 1970, no 68/25, pp. 68-92.

15. Orlenko L. P. Fizika vzryva i udara. 2-e izd. (Physics of explosion and shock. 2nd edition), Moscow, Fizmatlit, 2008, 304 p.

16. Bronshteyn I. N., Semendyaev K. A. Spravochnik po matematike dlya inzhenerov i uchashchikhsya vtuzov (Handbook of mathematics for engineers and pupils of technical colleges), Moscow, Glav.red.fizmatlit, 1980. 976.

17. Gorinov S. A., Kutuzov B. N. Gornyy informatsionno-analiticheskiy byulleten'. 2012, no 4, pp. 302-307.

18. Fokin V. A. Izvestiya vuzov. Gornyy zhurnal. 2007, no 3, pp. 89-94.