Научная статья на тему 'ВЛИЯНИЕ ИНТЕРВАЛОВ ЗАМЕДЛЕНИЯ НА АМПЛИТУДЫ ВОЛН НАПРЯЖЕНИЙ ПРИ ИЗУЧЕНИИ МОДЕЛИ ВЗРЫВА СИСТЕМЫ СКВАЖИННЫХ ЗАРЯДОВ'

ВЛИЯНИЕ ИНТЕРВАЛОВ ЗАМЕДЛЕНИЯ НА АМПЛИТУДЫ ВОЛН НАПРЯЖЕНИЙ ПРИ ИЗУЧЕНИИ МОДЕЛИ ВЗРЫВА СИСТЕМЫ СКВАЖИННЫХ ЗАРЯДОВ Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
9
1
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВЗРЫВ / КОРОТКОЗАМЕДЛЕННОЕ ВЗРЫВАНИЕ (КЗВ) / НАПРЯЖЕНИЕ / ИНТЕРФЕРЕНЦИЯ / ГОРНАЯ ПОРОДА / ГРАНИТ / СКВАЖИНА / КАРЬЕР / НЕГАБАРИТ

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

В настоящее время на многих карьерах существует проблема выхода некондиционных кусков горной массы после взрыва. Некондиционные куски чаще всего возникают в центральной части блока, между скважинами, вследствие плохой проработки этой зоны. Данная проблема возникает из-за того, что амплитуда волн напряжений от взрыва зарядов в этой области намного ниже, чем вблизи скважин. В данной статье приведен краткий обзор методик расчёта интервалов замедления; построена модель гранитного блока с взрывными скважинами, заряженными ANFO; проведено моделирование интерференции волн напряжения в программном комплексе Ansys; построены зависимости напряжений на фронте волны от времени протекания процесса; рассчитана экономическая эффективность применения электронной и неэлектрической (НСИ) систем инициирования. Моделирование показало, что при правильном подборе интервала замедления можно увеличить амплитуду волн напряжений в этой зоне за счет интерференции волн напряжений, тем самым улучшить дробление породы и снизить процент некондиционных кусков горной массы. Для обеспечения предложенного интервала замедления необходимо применять электронную систему инициирования I-kon III, которая обеспечивает минимальный шаг интервала замедления в 0,1 мс.

i Надоели баннеры? Вы всегда можете отключить рекламу.

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INFLUENCE OF DECELERATION INTERVALS ON THE AMPLITUDES OF STRESS WAVES DURING THE EXPLOSION OF A SYSTEM OF BOREHOLE CHARGES

Currently, many open pits there is a problem of substandard pieces of rock mass coming out after an explosion. Substandard pieces, most often, arise in the central part of the block, between the wells, due to poor development of this zone. This problem arises due to the fact that the amplitude of the stress waves from the explosion of charges in this area is much lower than near the wells. This article provides a brief overview of methods for calculating deceleration intervals; a model of a granite block with blastholes charged with ANFO was built; the simulation of the interference of voltage waves in the Ansys software package was carried out; the dependences of the voltages at the wave front on the time of the process are constructed, and the economic efficiency of the use of electronic and nonelectrical (NSI) initiation systems is calculated. Modeling showed that with the correct selection of the deceleration interval, it is possible to increase the amplitude of stress waves in this zone due to the interference of stress waves, thereby improving rock crushing and reducing the percentage of substandard pieces of rock mass. To ensure the proposed deceleration interval, it is necessary to use the electronic initiation system I-kon III, which provides a minimum step of the deceleration interval of 0.1 ms.

Текст научной работы на тему «ВЛИЯНИЕ ИНТЕРВАЛОВ ЗАМЕДЛЕНИЯ НА АМПЛИТУДЫ ВОЛН НАПРЯЖЕНИЙ ПРИ ИЗУЧЕНИИ МОДЕЛИ ВЗРЫВА СИСТЕМЫ СКВАЖИННЫХ ЗАРЯДОВ»

ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2022;(6-2):18—32 ОРИГИНАЛЬНАЯ СТАТЬЯ / ORIGINAL PAPER

УДК 622.235 DOI: 10.25018/0236_1493_2022_62_0_18

ВЛИЯНИЕ ИНТЕРВАЛОВ ЗАМЕДЛЕНИЯ НА АМПЛИТУДЫ ВОЛН НАПРЯЖЕНИЙ ПРИ ИЗУЧЕНИИ МОДЕЛИ ВЗРЫВА СИСТЕМЫ СКВАЖИННЫХ ЗАРЯДОВ

В. В. Должиков1, Д. Э. Рядинский1, А. А. Яковлев1

1 Санкт-Петербургский горный университет, Санкт-Петербург, 199106, Россия

Аннотация: В настоящее время на многих карьерах существует проблема выхода некондиционных кусков горной массы после взрыва. Некондиционные куски чаще всего возникают в центральной части блока, между скважинами, вследствие плохой проработки этой зоны. Данная проблема возникает из-за того, что амплитуда волн напряжений от взрыва зарядов в этой области намного ниже, чем вблизи скважин. В данной статье приведен краткий обзор методик расчёта интервалов замедления; построена модель гранитного блока с взрывными скважинами, заряженными ANFO; проведено моделирование интерференции волн напряжения в программном комплексе Ansys; построены зависимости напряжений на фронте волны от времени протекания процесса; рассчитана экономическая эффективность применения электронной и неэлектрической (НСИ) систем инициирования. Моделирование показало, что при правильном подборе интервала замедления можно увеличить амплитуду волн напряжений в этой зоне за счет интерференции волн напряжений, тем самым улучшить дробление породы и снизить процент некондиционных кусков горной массы. Для обеспечения предложенного интервала замедления необходимо применять электронную систему инициирования I-kon III, которая обеспечивает минимальный шаг интервала замедления в 0,1 мс.

Ключевые слова: взрыв, короткозамедленное взрывание (КЗВ), напряжение, интерференция, горная порода, гранит, скважина, карьер, негабарит.

Для цитирования: Должиков В. В., Рядинский Д. Э., Яковлев А. А. Влияние интервалов замедления на амплитуды волн напряжений при изучении модели взрыва системы скважинных зарядов // Горный информационно-аналитический бюллетень. — 2022. — № 6-2. — С. 18—32. DOI: 10.25018/0236_1493_2022_62_0_18.

Influence of deceleration intervals on the amplitudes of stress waves during the explosion of a system of borehole charges

V. V. Dolzhikov1, D. E. Ryadinsky1, A. A. Yakovlev1

1 St. Petersburg Mining University, St. Petersburg, 199106, Russia

Abstract: Currently, many open pits there is a problem of substandard pieces of rock mass coming out after an explosion. Substandard pieces, most often, arise in the central part of the block, between the wells, due to poor development of this zone. This problem arises due to the fact that the amplitude of the stress waves from the explosion of charges in this

© В. В. Должиков, Д. Э. Рядинский, А. А. Яковлев. 2022

area is much lower than near the wells. This article provides a brief overview of methods for calculating deceleration intervals; a model of a granite block with blastholes charged with ANFO was built; the simulation of the interference of voltage waves in the Ansys software package was carried out; the dependences of the voltages at the wave front on the time of the process are constructed, and the economic efficiency of the use of electronic and nonelectrical (NSI) initiation systems is calculated. Modeling showed that with the correct selection of the deceleration interval, it is possible to increase the amplitude of stress waves in this zone due to the interference of stress waves, thereby improving rock crushing and reducing the percentage of substandard pieces of rock mass. To ensure the proposed deceleration interval, it is necessary to use the electronic initiation system I-kon III, which provides a minimum step of the deceleration interval of 0.1 ms.

Key words: blast, short-delayed blasting, voltage, interference, rock, granite, borehole, quarry, oversized.

For citation: Dolzhikov V. V., Ryadinsky D. E., Yakovlev A. A. Influence of deceleration intervals on the amplitudes of stress waves during the explosion of a system of borehole charges. MIAB. Mining Inf. Anal.. Bull. 2022;(6-2):18-32. [In Russ]. DOI: 10.25018/0236_1493_2022_62_0_18.

Введение

Дробление скальных пород взрывным способом широко применяется в горнодобывающей промышленности при разработке карьеров. Затраты на буровзрывные работы составляют до 35% от общей стоимости добычи горной массы на карьерах [1-4].

Как правило, часть энергии взрыва расходуется непосредственно на разрушение горной породы в пределах разрушаемого массива, в то время как остальная часть энергии расходуется на бесполезные формы работы, например, образование воздушной и сейсмической волн [5-8].

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

При взрыве горных пород метод короткозамедленного взрывания является эффективной мерой для улучшения разрушения горных пород и уменьшения вызываемой взрывом вибрации за счет управления временем инициирования и последовательностью детонаторов [15-19].

Определение рационального интервала замедления для короткозамедлен-ного взрывания является ключевым фактором для повышения качества взрывных работ [20-23]. Несмотря на то, что существуют формулы расчёта параметров короткозамедленного взрывания, нет единой точки зрения на физическую сущность процесса, и все они считаются ориентировочными и уточняются экспериментально [24, 25, 26].

В соответствии с вышесказанным, научная проблема обоснования технологических параметров КЗВ является актуальной.

Методы исследования

и исходные данные

Существует множество методик по расчету интервалов замедления. В своей работе мы рассмотрели несколько методик по расчету времени замедления.

Расчеты производились на примере карьера «Прудянский», где разрабатывается гранитогнейсовое месторождение. Из типового проекта были выбраны интересующие нас параметры: d^ = = 105 мм, qBB = 1 кг/м3; Aanfo= 931 кг/м3

А. Н. Ханукаев предложил формулу, которая гласит, что интервал замедления должен соответствовать времени, необходимому для образования новых поверхностей [27]:

tv = ti + í2 + Í3, мс, (1)

где íDm. — время, необходимое для образования новой дополнительной плоскости обнажения, мс; t1 — время, необходимое для распространения упругой волны от зарядной камеры до поверхности и обратно, мс; t2 — время, необходимое на образование трещин, распространяющихся от поверхности к зарядной камере, мс (длиной трещин, образующихся у зарядной камеры, можно пренебречь ввиду приближенности расчетов); t3 — время, необходимое на перемещение породы на расстояние 0,8-1 см от целика, мс,

2W W ст

^ = ~с~+с-+Г; (2)

Р тр ср

где Cp — скорость звука в среде, м/с; Стр — скорость трещинообразования, обычно принимается 0,05 Cp, м/с; vср — скорость раскрытия щели, vcp < 3 м/с; о — ширина щели (0,008-0,01), мм.

Было проведено несколько расчетов для разных скоростей продольных волн.

В расчете № 1 — принималась скорость продольной волны, опубликованная Межведомственной комиссией по взрывному делу — 4500 м/с;

В расчете № 2 — принималась скорость продольной волны в идеальном массиве гранитных пород — 5200 м/с;

В расчете № 3 — принималась скорость продольной волны, изме-

ренная научным коллективом Санкт-Петербургского горного университета на гранитных карьерах Ленинградской области - 4710 м/с [12].

Интервал замедления для скорости продольной волны, указанной в МКВД:

2 • 3 3 0,009

? =--1---ъ —-=

опт 4500 0,05 • 4500 3,956

= 0,0013 + 0,01333 + 0,00227 = 16,9 мс.

Интервал замедления для идеальной скорости продольной волны в массиве:

2 • 3

t„nT = —— + -

0,009

5200 0,05 • 5200 3,956 = 0,00115 + 0,0115 + 0,00227 = 14,92 мс.

Интервал замедления для скорости продольной волны, полученной с помощью сейсмостанций:

l. =-23+-

- + -

0,009

4710 0,05 • 4710 3,956 = 0,00127 + 0,01274 + 0,00227 = 16,28 мс.

У. Лангефорс предложил эмпирическую зависимость:

т = 3,3 • К • W, мс (3)

где К — коэффициент, зависящий от свойств пород (К = 1 для крепких пород, К = 2 — для слабых пород); W — линия наименьшего сопротивления, м.

т = 3,3 • 1 • 3 = 9,9, мс.

Г. И. Покровский утверждает, что для увеличения эффективности дробления породы взрывом необходимо воспользоваться следующей зависимостью [11]:

л/L2 + 4W2 т =-, мс,

(4)

где L — расстояние между зарядами, м; W — глубина заложения заряда или линия наименьшего сопротивления

(ЛНС), м; v — скорость распространения деформаций в породе, м/с;

т =

У32 + 4 • 32 3,956

= 1,7 mc .

Л. В. Баранов предлагает следующую зависимость [11]:

Z v

т = — K, mc ,

Q

(5)

где 2 — глубина скважин, м; Q — величина заряда, кг; К — коэффициент пропорциональности;

13 2

т = ^£.2 = 0,26мс . 101

Н. Д. Евстропов предлагает формулу:

(6)

2W

т = K-, mc ,

С

где № — ЛНС, м; К — коэффициент, зависящий от крепости породы (1 < К < 2); Ср- скорость звука в среде, м/с.

т = = 0,0017 мс

5200

(при Ср = 5200 мс)

т = 1.5"213 = 0,002 мс 4500

(при Ср = 4500 мс)

т = ХЬ^3 = 0,0019 мс 4710

(при Ср = 4710 мс)

А. Ф. Белаенко, Ю. В. Гаек предлагают формулу:

т = т +

W

1 Стр-ц-cosp/2 T•W2 •ytgР/2

+-' —, mc

d

(7)

где т1 — продолжительность эффективного воздействия прямого поля напряжений (принимается 5-7 мс); № — ЛНС, м; Стр — скорость развития трещин в сплошном массиве при заданном удельном расходе ВВ, м/ с; | — коэффициент трещиноватости среды, ед.; в — угол раскрытия призмы выброса, град.; Т — временной коэффициент, определяемый эмпирическим путем (в среднем 2-10-6 с); d — диаметр скважин, м; у — средний удельный вес пород в уступе, г/см3;

т = 6 + -

3

0,05 • 5400 • 0,5 • 0,906 2• 10-6 • 32 • 2,7 • 0,4663

0,105

= 11,249 mc .

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

На данном этапе развития науки и техники минимально возможный интервал замедления составляет 0,1 мс (электронный системы инициирования I-ron III). В связи с этим моделирование с меньшими интервалами не имеет смысла.

Моделирование производилось в программном комплексе Ansys Autodyn. Решатели Autodyn можно разделить на следующие категории: Лагранжевы, Эйлеровы, arbitrary Lagrangian-Euler ALE (произвольный Лагранжево-Эйлеров метод), Smoothed Particle Hydrodynamics SPH (метод сглаженных частиц). Для правильного выбора того или иного решателя необходимо руководствоваться двумя критериями: точностью и производительностью [28, 29, 30].

При подходе Лагранжа к описанию движения сплошной среды расчетная сетка является «вмороженной» в материал, она движется и деформируется совместно с ним (рис. 1). Лагранжевы

+

Таблица 1

Значения интервалов замедлений для различных методик Deceleration interval values for various techniques

№ Методика расчета Значение интервала замедления для условий карьера «Прудянский»

1. А.Н. Ханукаева 16,9 мс (при Ср = 4500 м/с)

14,92 мс (при Ср = 5200 м/с)

16,2 мс (при Ср = 4710 м/с)

2. У. Лангефорса 9,9 мс

3. Г.И. Покровского 1,7 мс

4. Л.В. Баранова 0,26 мс

5. Н.Д. Евстропова 0,0017 мс (при Ср = 5200 м/с)

0,002 мс (при Ср = 4500 м/с)

0,0019 мс (при Ср = 4710 м/с)

6. А.Ф. Белаенко и Ю. В. Гаека 11,249 мс

Рис. 1. Расчетная сетка в Лагранжевом представлении: а — в начальный момент времени; б — под действием нагружения Fig. 1. The computational grid in the Lagrangian representation: a — at the initial moment of time; b — under the action of loading

решатели представляют собой наиболее точный и эффективный метод расчета нелинейных задач динамического нагружения конструкций. Основным недостатком Лагранжевых методов являются чрезмерные искажения и «запутывание» расчетной сетки при больших деформациях, приводящие к потерям точности в этих областях, а иногда и к полному прекращению расчета. В AUTODYN эта проблема может быть частично решена при использовании искусственного алгоритма эрозии, который удаляет из расчета чрезмерно деформированные элементы.

Рис. 2. Расчетная сетка в Эйлеровом представлении: а — под действием нагружения; б — при развитии нагружения Fig. 2. The computational mesh in the Euler representation: a — under the action of loading; b — with the development of loading

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

Основные законы, используемые

решателем AUTODYN

Закон сохранения массы:

dP + = о dt ax,

Ф +ра(рО =0.

д? дх ,

Закон сохранения момента:

бч 1 дс,,

^ = I + = 0 , М ' Р дх 1

дЧ: дЧ: х 1 сс ,, —'- + V —'- = ^ +---

д? х 1 р дх, Закон сохранения энергии:

(8)

(9)

ЭЕ ЭЕ

— + V,—

Э? ' Эх,

(

Эр Эр

— + V : —

Э? ' Эх

л

I

Р = А ехр- Я^ + В ехр- Я2ч + Сч +1) .(10)

Для определения поведения горной породы, используемой в данном моделировании, было предложено использовать уравнение прочности Риделя — Хармайера — Тома (ЯНТ), которое описывает механизм простого сдвига для объяснения ориентировки разломов в однородной среде при преобладании тангенциальных напряжений в пределах сдвиговой зоны. В таких условиях возникают целые серии сколов, одни из которых являются синтетическими (совпадающими по направлению смещения с общим смещением

сдвиговой зоны), а другие — антитетическими (с противоположным смещением) [31].

В табл. 2 приведены параметры модели, описывающие поведение гранита.

Для определения поведения взрывчатого вещества и продуктов детонации, используемого в инженерном расчете, требуется уравнение состояния, которое определяет взаимосвязь между параметрами состояния взрывчатого вещества. Для численного моделирования использовалось уравнение состояния Джонса -Уилкинса — Ли (JWL), поскольку оно может точно описывать поведение различных взрывчатых веществ. Формула для уравнения состояния JWL представлена выражением (11).

Р = Аехр - + Вехр - В2у + Су>+1) , (11)

где А, К, В, К, С, w — параметры уравнения состояния, вычисляемые эмпирическим путем [16].

Используя методику, предложенную в работе [14], на основании имеющихся в базе данных для PETN с четырьмя различными плотностями, были извлечены константы из стандартной библиотеки Autodin для изучения применимых корреляций. Уравнения регрессии для них выглядят следующим образом [13]:

А = -2, 77 • 103 +1, 050 - 9 , 66 • 10-5 • О2 + 2, 76 • 10-9 • О3 , В = 3 , 56 • 10-2 + 2 , 11РС| - 4 , 8 • 10-2 • Р-2,

= 8,82 - 1,63рвв - 4,91 • 10-1 р2вв, ^2 = 2,03 +1,77 • 10-1 Рвв2 - 2,49 • 10-1 рВв, п = 3,47 • 10-2 + 2,7 • 10-1 рвв - 4,7 • 10-2 рвв,

(12)

где D — скорость детонации, м/с; Р^ — детонационное давление, кПа; рВВ — плотность взрывчатого вещества (тэна), г/см3.

В табл. 3 приведены параметры модели, описывающие поведение ANFO.

Гранитный блок представлен на рис. 3. Размеры гранитного блока в модели составили 3000 мм в ширину, 6000 мм в длину и 8000 мм в высоту. Расстояние между скважинами составляет 3000 мм.

Таблица 2

Параметры вещественной модели, описывающие поведение гранита Material model parameters describing the behavior of granite

Параметр Единица измерения Значение

Плотность г/см3 2,66

Уравнение состояния (EOS) - полиномиальный

Модуль объемной упругости ГПа 43,87

Температура К 300

Удельная теплоемкость Дж/кг 654

Уравнение прочности — модель Риделя-Хармайера-Тома (Riedel, HiermaierandThoma, RHT)

Модуль сдвига ГПа 17

Прочность при сжатии МПа 150

Предел прочности при разрыве — 0,05

Предел прочности на сдвиг — 0,07

Соотношение растяжение/сжатие — 0,72

Переход хрупкости/ вязкости — 0,01

Постоянная трещинообразования, В — 2,5

Постоянная трещинообразования, М — 0,85

Константа скорости деформации — 0,025

сжатия

Предел скорости деформации сжатия — 0,045

Модель разрушения — модель Риделя-Хармайера-Тома, RHT

Постоянная разрушения, D1 — 0,025

Постоянная разрушения, D2 — 1,0

Минимальная деформации разрушения — 0,06

Остаточный модуль сдвига — 0,25

Модель растяжения разрушения — гидро

Модель эрозии — Пластическая деформация

Таблица 3

Параметры вещественной модели, описывающие поведение ANFO Real model parameters describing ANFO behavior

Параметр Единица измерения Значение

Плотность кг/м3 0,931

Скорость детонации м/с 4160

Давление ГПа 5,14

Уравнение состояния (EOS) — модель Джонса — Уилкинса — Ли (Jones -Wilkins — Lee, JWL)

Автоконвертирование в идеальный — да

газ

А ГПа 49,4

В ГПа 1,891

Ri — 3,907

R, — 1,118

w — 0,333

Материалы

П Скважина №1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Воздух j

J Скважина №2

ANFO Гранит

Рис. 3. Модель гранитного блока со скважинами, заряженными ANFO Fig. 3. Model of a granite block with boreholes charged by ANFO

Для воздушной среды был использован идеальный газ с давлением, равным 0,1 МПа. Учитывалось действие гравитации земли. Скважины были диаметром 105 мм, что также соответствует типовому проекту.

В центральной точке гранитного блока был установлен датчик, считывающий показания (3000; 1500; 4000).

Результаты

Проведя серию опытов по моделированию взрыва скважин с замедлениями, указанными в табл. 1, максимальные амплитуды волн были получены с замедлениями 16,2 и 16,9 мс, которые соответствуют скорости звука в среде 4500 м/с и 4710 м/с соответственно. В результате серии опытов выявлено,

что наибольшая интерференция волн возникала при взрыве второго заряда с замедлением 16,2 мс.

В настоящее время для инициирования скважин зарядов на большинстве карьеров применяются неэлектрические системы инициирования (НСИ), например «Искра», «Эксэл», «Коршун» и др. с номинальными интервалами замедления 0, 9, 17, 25, 42, 67, 109, 176, 285 мс. Сравним с ближайшим номинальным значением интервала замедления в стандартном ряду неэлектрических систем инициирования — 17 мс.

Первый опыт был проведен с временем замедления 17 мс, а второй опыт с замедлением 16,2 мс.

На рис. 4 и 5 представлены зависимости напряжений на фронте волны от времени протекания процесса.

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

На рис. 5 видно, что напряжение, возникающее при взрыве заряда через 16,2 мс, составляет 1,273*105 КПа, что на 8% выше, чем при взрыве заряда через 17 мс — 1,18*105 КПа.

Обсуждение результатов

Добиться высокой точности интервала замедления при помощи НСИ невозможно, т. к. фактический разброс времени срабатывания детонаторов составляет 10-15% от номинального интервала замедления.

В связи с этим для реализации данной методики необходимо перейти на электронные системы инициирования, которые в свою очередь дороже неэлектрических, но позволяют задать интервал замедления с точностью до 0,1 мс [32].

Экономическая эффективность проекта рассчитывалась для блока, состоящего из 70 скважин, объемом 15000 м3.

1-г

5 10 15 20

Время, мс Время, мс

а б

Рис. 4. Зависимости напряжений на фронте волны от времени протекания процесса: а — время замедления 16,2 мс, б — время замедления 17 мс

Fig. 4. Dependences of voltages at the wave front on the time of the process: a — the deceleration time is 16.2 ms, b — the deceleration time is 17 ms

Рис. 5. Зависимости напряжений на фронте волны от времени протекания процесса Fig. 5. Dependences of stresses at the wave front on the time of the process

Все показатели БВР, ВВ были приняты согласно проекту на массовый взрыв Прудянского карьера. Менялась только система инициирования.

Выход негабарита на карьере составляет до 20%. Основываясь на аналогичных статьях и опытах, можно сделать вывод, что за счет интервала замедле-

ния, рассчитанного по данной методике, можно добиться уменьшения выхода негабарита до 7% [33, 34, 35].

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

Таблица 4

Расчет экономической эффективности Calculation of economic efficiency

Наименование Единица измерений НСИ ИСКРА ЭСИ i-kon

Количество скважин шт. 70

Цена детонаторов руб. 85 1200

Суммарная стоимость руб. 5 950 84 000

Объем блока м3 15 000

Цена за взорванную горную массу руб./м3 100

Суммарная стоимость взорванной руб. 1500000

горной массы

Выход негабарита м3 3 000 1 050

Цена на дробление негабарита руб./ м3 100

Суммарная стоимость на дробление руб. 300000 105000

негабарита

Прибыль с блока руб. 1 194 050 1311 000

Из расчетов видно, что затраты на НСИ составили приблизительно 6000 руб., а для ЭСИ — 84000 руб., в то время как затраты на негабарит у НСИ гораздо выше, так как с помощью НСИ нельзя точно подобрать интервал замедления и, как следствие, выход негабарита увеличивается

Сложив затраты на негабарит и стоимость средств инициирования в результате получаем, что приблизительная прибыль при точном подборе времени замедления с использованием ЭСИ выше на 116 500 руб.

Следует отметить, что данная оценка экономической эффективности весьма приблизительна.

Заключение

Результаты математического моделирования показывают, что интервал замедления влияет на амплитуды волн напряжений при взрыве системы сква-жинных зарядов. Так, напряжение, возникающее при моделировании взрыва двух скважин с замедлением 16,2 мс, составляет 1,273*105 КПа, что на 8% выше, чем при взрыве заряда через 17 мс — 1,18*105 КПа.

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

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

Планируется подробнее рассмотреть методики, которые предполагают, что увеличение интервала замедления будет способствовать повышению качества дробления.

Вклад авторов

Должиков В. В. — разработка идеи статьи, постановка задачи исследования, выполнение анализа полученных результатов.

РядинскийД. Э. - проведен научный обзор по теме статьи в отечественной и зарубежной литературе.

Яковлев А. А. - выполнено моделирование взрыва скважинных зарядов.

Авторы заявляют об отсутствии конфликта интересов.

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

1. Симонов П. С. Особенности определения размера среднего куска и выхода негабарита при взрывных работах на карьерах // Горный информационно-аналитический бюллетень. - 2017. - № 4. - С. 320-327.

2. Jia B., Ling T., Hou S., Liu D. Application of Variational Mode Decomposition Based Delay Time Identification in Short Millisecond Blasting // Transaction of Beijing Institute of Technology. 2021, vol. 41 (4), pp. 341-348. DOI: 10.15918/j.tbit1001-0645.2019.308.

3. Yl C., Sjoberg J., Johansson D. Numerical modeling forblast-inducedfragmentationin sublevel caving mine // Tunneling and Underground Space Technology. 2017, vol. 68, pp. 167-173. DOI: 10.1016/j.tust.2017.05.030.

4. Gou Y., ShiX., QiuX., Huo X., YuZ. Assessment of induced vibrations derived from the wave superposition in time-delayblasts // International Journal of Rock Mechanics and Mining Sciences. 2021, vol. 144, 104814. DOI: 10.1016/j.ijrmms.2021.104814.

5. Yl C., Johansson D., Greberg J. Effects of in-situ stresses on the fracturing of rock by blasting // Computers and Geotechnics. 2018, vol. 104, pp. 321-330. DOI: 10.1016/j. compgeo.2017.12.004.

6. Yastrebova K. N., Moldovan D. V., Chernobay V. I. Solving the issue of ventilating atmosphere of opencast mining by resloping bench face // International Journal of Advanced Science and Technology. 2020, vol. 29 (1), pp. 1-6.

7. Silva J., Li L., Gernand J. M. Reliability analysis for mine blast performance based on delay type and firing time // International Journal of Mining Science and Technology. 2018, vol. 28 (2), pp. 195-204. DOI: 10.1016/j.ijmst.2017.07.004.

8. Karakus M., Ebong U., Liu C., Zhou H. Three-dimensional finite element analysis for rock fatigue damage under dynamic loading // ISRM Regional Symposium, EUROCK. 2015, pp. 577-582.

9. Marinin M., Marinina O., Wolniak R. Assessing of losses and dilution impact on the cost chain: Case study of gold ore deposits // Sustainability (Switzerland). 2021, vol.13 (7), 3830. DOI: 10.3390/su13073830.

10. Paramonov G. P., Kovalevskyi V. N., Mysin A. V. Determination of the conditions of an effective functioning of elongated cumulative charges in processing the marble stone // Key Engineering Materials. 2020, vol. 836, pp. 19-24. DOI: 10.4028/www.scientific. net/KEM.836.19.

11. Blair D. P. Limitations of electronic delays for the control of blast vibration and fragmentation. // Rock Fragmentation by Blasting: Proceedings of the 9th International Symposium on Rock Fragmentation by Blasting, FRAGBLAST. 2010, vol. 9, pp. 171-184.

12. Wu H., Gong M. Calculation and application of hole by hole blasting vibration superposition based on measured delay times of detonators // Explosion and Shock Waves. 2019, vol. 39 (2), 025202. DOI: 10.11883/bzycj-2017-0415.

13. Roy M. P., Mishra A. K., Agrawal H., Singh P. K. Blast vibration dependence on total explosives weight in open-pit blasting // Arabian Journal of Geosciences. 2020, vol. 13, no. 13, 531, pp. 1-8.

14. Хохлов С. В., Рахманов Р. А., Аленичев И. А., Баженова А. В., Маккоев В. А. Исследование вопроса управления и контроля за смещением контуров рудных тел после взрыва // Взрывное дело. - 2021. - № 132 (89). - С. 59-75.

15. Wang Z., Fang C., Chen Y., Cheng W. A Comparative study of delay time identification by vibration energy analysis in millisecond blasting // International Journal of Rock Mechanics and Mining Sciences. 2013, vol. 60, pp. 389-400. DOI: 10.1016/j.ijrmms.2012.12.032.

16. Fu H., Wong L. N. Y., Zhao Y., Shen Z., Zhang C., Li Y. Comparison of excavation damage zones resulting from blasting with nonel detonators and blasting with electronic detonators // Rock Mechanics and Rock Engineering. 2014, vol. 47 (2), pp. 809-816. DOI: 10.1007/s00603-013-0419-2.

17. Ковалевский В. Н., Аргимбаев К. Р. Экспериментальное исследова-ние процесса внедрения кумулятивной струи в массив горных пород // Горный журнал. — 2016. — № 12. — С. 19-23. DOI: 10.17580/gzh.2016.12.04.

18. Paramonov G. P., Kovalevskiy V. N., Mysin A. V. Impact of muLticharge detonation on explosion puLse value // IOP Conference Series: Earth and Environmental Science. 2018, vol. 194 (8), 082031. DOI: 10.1088/1755-1315/194/8/082031.

19. Yastrebova K., Moldovan D., Chernobay V. Influence of the nature of the outflow of explosion products from blast holes and boreholes on the efficiency of rock destruction // E3S Web of Conferences. 2020. vol. 174, 01017. DOI: 10.1051/e3sconf/202017401017.

20. Wang Z., Fang C., Chen Y., Cheng W. A comparative study of delay time identification by vibration energy analysis in millisecond blasting // International Journal of Rock Mechanics and Mining Sciences. 2013, vol. 60, pp. 389-400. DOI: 10.1016/j.ijrmms.2012.12.032.

21. Shi X. Z., Chen S. H. R. Delay time optimization in blasting operations for mitigating the vibration-effects on final pit walls' stability // Soil Dynamics and Earthquake Engineering. 2011, vol. 31 (8), pp. 1154-1158.DOI: 10.1016/j.soildyn.2011.04.004.

22. Zhang S., Ling 7.-H., Liu H.-R., Cao F. Pattern adapted wavelet time-energy density method and its application in millisecond blast vibration signal analysis // Journal of the China Coal Society. 2014, vol. 39 (10), pp. 2007-2013. DOI: 10.13225/j.cnki.jccs.2014.0325.

23. Zhang Y., Chen Y., Chen S., Liu H., FuZ. Experimental study on deformation of a sandy field liquefied by blasting // Soil Dynamics and Earthquake Engineering. 2019, vol. 116, pp. 60-68. DOI: 10.1016/j.soildyn.2018.09.042.

24. Хохлов С. В, Соколов С. Т., Виноградов Ю. И., Френкель И. Б. Проведение промышленных взрывов вблизи газопроводов // Записки Горного института. — 2021. — Т. 247. — С. 48-56. DOI:10.31897/PMI.2021.1.6.

25. Qiu X., Shi X., Gou Y., Zhou J., Chen H., Huo X. Short-delay blasting with single free surface: Results of experimental tests // Tunnelling and Underground Space Technology. 2018, vol. 74, pp. 119-130. DOI:10.1016/j.tust.2018.01.014.

26. Huang D., Qiu X., Shi X., Gou Y., Zhou J. Experimental and Numerical Investigation of Blast-Induced Vibration for Short-Delay Cut Blasting in Underground Mining // Shock and Vibration. 2019, vol. 2019, 5843516. DOI: 10.1155/2019/5843516.

27. Камянский В. Н. Повышение эффективности скважинной отбойки на карьерах при разновременном взрывании скважинных зарядов: Автореф... дис. канд. техн. наук. — М: ГоИ КНЦ РАН, 2019. — 24 с.

28. WuX., Gong M., WuH., LiuX. Parameter calculation of the initiating circuit with mixed use of nonel detonators and electronic detonators in tunnel controlled-blasting // Tunnelling and Underground Space Technology. 2021, vol. 113, 103975. DOI: 10.1016/j.tust.2021.103975.

29. Isheyskiy V., Sanchidrian J. A. Prospects of applying MWD technology for quality management of drilling and blasting operations at mining enterprises // Minerals, 2020, vol. 10 (10), 925, pp. 1-17. DOI: 10.3390/min10100925.

30. Маринин М. А., Хохлов С. В., Ишейский В. А. Моделирование режима протекания процесса сварки плоских листовых деталей взрывом // Записки Горного института. — 2019. — Т. 237. — С. 275-280. DOI: 10.31897/PMI.2019.3.275.

31. Afanasev P. I., MakhmudovK. F. Assessment of the parameters of a shock wave on the wall of an explosion cavity with the refraction of a detonation wave of emulsion explosives // Applied Sciences (Switzerland). 2021, vol. 11 (9), 3976. DOI: 10.3390/app11093976.

32. Аленичев И. А. Реакция массива горных пород в карьерном пространстве на динамические воздействия при производстве взрывных работ // Горный информационно-аналитический бюллетень. — 2018. — № 7. — С. 189-195. DOI: 10.25018/0236-14932018-7-0-189-195.

33. Andreev R. E., Gridina E. B. A Study of gas-dynamic processes in a charge chamber during the explosion of blasthole charges of various designs // Research Journal of Pharmaceutical, Biological and Chemical Sciences. 2016, vol. 7 (3), pp. 2383-2392.

34. Тюпин В. Н., Рубашкина Т. И. Инженерные формулы для определения размеров зон деформаций и повреждений в трещиноватом массиве горных пород при воздействии взрывных работ на карьерах Забайкалья // Горный журнал. - 2021. - № 7. - С. 40-44. DOI: 10.17580/gzh.2021.07.06.

35. Cardu M., Giraudi A., Oreste P. A review of the benefits of electronic detonators // Revista Escola de Minas. 2013, vol. 66 (3), pp. 375-382. DOI: 10.1590/S0370-44672013000300016.ЕШ2

REFERENCES

1. Simonov P. S. Peculiarities of determining the size of an average piece and the output of oversized pieces during blasting in open pits. MIAB. Mining Inf. Anal. Bull. 2017, no. 4, pp. 320-327. [In Russ].

2. Jia B., Ling T., Hou S., Liu D. Application of Variational Mode Decomposition Based Delay Time Identification in Short Millisecond Blasting. Transaction of Beijing Institute of Technology. 2021, vol. 41 (4), pp. 341-348. DOI: 10.15918/j.tbit1001-0645.2019.308.

3. Yi C., Sjoberg J., Johansson D. Numerical modeling for blast-induced fragmentation in sublevel caving mines. Tunneling and Underground Space Technology. 2017, vol. 68, pp. 167-173. DOI: 10.1016/j.tust.2017.05.030.

4. Gou Y., Shi X., Qiu X., Huo X., Yu Z. Assessment of induced vibrations derived from the wave superposition in time-delay blasts. International Journal of Rock Mechanics and Mining Sciences. 2021, vol. 144, 104814. DOI: 10.1016/j.ijrmms.2021.104814.

5. Yi C., Johansson D., Greberg J. Effects of in-situ stresses on the fracturing of rock by blasting. Computers and Geotechnics. 2018, vol. 104, pp. 321-330. DOI: 10.1016/j. compgeo.2017.12.004.

6. Yastrebova K. N., Moldovan D. V., Chernobay V. I. Solving the issue of ventilating atmosphere of opencast mining by resloping bench face. International Journal of Advanced Science and Technology. 2020, vol. 29 (1), pp. 1-6.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Silva J., Li L., Gernand J. M. Reliability analysis for mine blast performance based on delay type and firing time. International Journal of Mining Science and Technology. 2018, vol. 28 (2), pp. 195-204. DOI: 10.1016/j.ijmst.2017.07.004.

8. Karakus M., Ebong U., Liu C., Zhou H. Three-dimensional finite element analysis for rock fatigue damage under dynamic loading. ISRM Regional Symposium, EUROCK2015, pp. 577-582.

9. Marinin M., Marinina O., Wolniak R. Assessing of losses and dilution impact on the cost chain: Case study of gold ore deposits. Sustainability (Switzerland). 2021, vol.13 (7), 3830. DOI: 10.3390/su13073830.

10. Paramonov G. P., Kovalevskyi V. N., Mysin A. V. Determination of the conditions of an effective functioning of elongated cumulative charges in processing the marble stone. Key Engineering Materials. 2020, vol. 836, pp. 19-24. DOI: 10.4028/www.scientific.net/ KEM.836.19.

11. Blair D. P. Limitations of electronic delays for the control of blast vibration and fragmentation. Rock Fragmentation by Blasting. Proceedings of the 9th International Symposium on Rock Fragmentation by Blasting, FRAGBLAST. 2010, vol. 9, pp. 171-184.

12. Wu H., Gong M. Calculation and application of hole by hole blasting vibration superposition based on measured delay times of detonators. Explosion and Shock Wives. 2019, vol. 39 (2), 025202. DOI: 10.11883/bzycj-2017-0415.

13. Roy M. P., Mishra A. K., Agrawal H., Singh P. K. Blast vibration dependence on total explosives weight in open-pit blasting. Arabian Journal of Geosciences. 2020, vol. 13, no. 13, 531, pp. 1-8.

14. Khokhlov S. V., Rakhmanov R. A., Alenichev I. A., Bazhenova A. V., Makkoev V. A. Study of the issue of management and control of the displacement of the contours of ore bodies after the explosion. Vzryvnoe delo. 2021, no. 132/89, pp. 59-75. [In Russ].

15. Wang Z., Fang C., Chen Y., Cheng W. A. Comparative study of delay time identification by vibration energy analysis in millisecond blasting. International Journal of Rock Mechanics and Mining Sciences. 2013, vol. 60, pp. 389-400. DOI: 10.1016/j.ijrmms.2012.12.032.

16. Fu H., Wong L. N. Y., Zhao Y., Shen Z., Zhang C., Li Y. Comparison of excavation damage zones resulting from blasting with nonel detonators and blasting with electronic detonators. Rock Mechanics and Rock Engineering. 2014, vol. 47 (2), pp. 809-816. DOI: 10.1007/s00603-013-0419-2.

17. Kovalevskiy V. N., Argimbaev K. R. Experimental research of explosive jet penetration in rocks. GornyiZhurnal. 2016, no. 12, pp. 19-23. [In Russ]. DOI: 10.17580/gzh.2016.12.04.

18. Paramonov G. P., Kovalevskiy V. N., Mysin A. V. Impact of multicharge detonation on explosion pulse value. IOP Conference Series: Earth and Environmental Science. 2018, vol. 194 (8), 082031. [In Russ]. DOI: 10.1088/1755-1315/194/8/082031.

19. Yastrebova K., Moldovan D., Chernobay V. Influence of the nature of the outflow of explosion products from blast holes and boreholes on the efficiency of rock destruction. E3S Web of Conferences. 2020, vol. 174, 01017. DOI: 10.1051/e3sconf/202017401017.

20. Wang Z., Fang C., Chen Y., Cheng W. A comparative study of delay time identification by vibration energy analysis in millisecond blasting. International Journal of Rock Mechanics and Mining Sciences. 2013, vol. 60, pp. 389-400. DOI: 10.1016/j.ijrmms.2012.12.032.

21. Shi X. Z., Chen S. H. R. Delay time optimization in blasting operations for mitigating the vibration-effects on final pit walls' stability. Soil Dynamics and Earthquake Engineering. 2011, vol. 31 (8), pp. 1154-1158.DOI: 10.1016/j.soildyn.2011.04.004.

22. Zhang S., Ling T.-H., Liu H.-R., Cao F. Pattern adapted wavelet time-energy density method and its application in millisecond blast vibration signal analysis. Journal of the China Coal Society. 2014, vol. 39 (10), pp. 2007-2013. DOI: 10.13225/j.cnki.jccs.2014.0325.

23. Zhang Y., Chen Y., Chen S., Liu H., Fu Z. Experimental study on deformation of a sandy field liquefied by blasting. Soil Dynamics and Earthquake Engineering. 2019, vol. 116, pp. 60-68. DOI: 10.1016/j.soildyn.2018.09.042.

24. Khokhlov S. V., Sokolov S. T., Vinogradov Y. I., Frenkel I. B. Conducting industrial explosions near gas pipelines. Journal of Mining Institute. 2021, vol. 247, pp. 48-56. DOI:10.31897/PMI.2021.1.6.

25. Qiu X., Shi X., Gou Y., Zhou J., Chen H., Huo X. Short-delay blasting with single free surface: Results of experimental tests. Tunnelling and Underground Space Technology. 2018, vol. 74, pp. 119-130. DOI:10.1016/j.tust.2018.01.014.

26. Huang D., Qiu X., Shi X., Gou Y., Zhou J. Experimental and Numerical Investigation of Blast-Induced Vibration for Short-Delay Cut Blasting in Underground Mining // Shock and Vibration. 2019, vol. 2019, 5843516. DOI: 10.1155/2019/5843516.

27. Kamyansky V. N. Increasing the efficiency of borehole breaking in open pits with multitemporal blasting of borehole charges: Abstract of ... dis. cand. tech. Sciences. Moscow, GoI KNTS RAN, 2019, 24 p.

28. Wu X., Gong M., Wu H., Liu X. Parameter calculation of the initiating circuit with mixed use of nonel detonators and electronic detonators in tunnel controlled-blasting. Tunnelling and Underground Space Technology. 2021, vol. 113, 103975. DOI: 10.1016/j. tust.2021.103975.

29. Isheyskiy V., Sanchidrian J. A. Prospects of applying MWD technology for quality management of drilling and blasting operations at mining enterprises. Minerals, 2020, vol. 10 (10), 925, pp. 1-17. DOI: 10.3390/min10100925.

30. Marinin M. A., Khokhlov S. V., Isheisky V. A. Simulation of the process of explosion welding of flat sheet parts. Journal of Mining Institute. 2019, vol. 237, pp. 275-280. DOI: 10.31897/PMI.2019.3.275.

31. Afanasev P. I., Makhmudov K. F. Assessment of the parameters of a shock wave on the wall of an explosion cavity with the refraction of a detonation wave of emulsion explosives. Applied Sciences (Switzerland). 2021, vol. 11 (9), 3976. DOI: 10.3390/app11093976.

32. Alenichev I. A. Reaction of a rock mass in a quarry space to dynamic impacts in the production of blasting. MIAB. Mining Information and Analytical Bulletin. 2018, no. 7. pp. 189-195. DOI: 10.25018 / 0236-1493-2018-7-0-189-195.

33. Cardu M., Giraudi A., Oreste P. A review of the benefits of electronic detonators. Revista Escola de Minas. 2013, vol. 66 (3), pp. 375-382. DOI: 10.1590/50370-44672013000300016.

34. Andreev R. E., Gridina E. B. A study of gas-dynamic processes in a charge chamber during the explosion of blasthole charges of various designs. Research Journal of Pharmaceutical, Biological and Chemical Sciences. 2016, vol. 7 (3), pp. 2383-2392.

35. Tyupin V. N., Rubashkina T. I. Engineering formulas for determining the size of deformation and damage zones in a fractured rock mass under the influence of blasting in open pits in Transbaikalia. Mining Journal. 2021, no. 7, pp. 40-44. DOI: 10.17580/gzh.2021.07.06.

ИНФОРМАЦИЯ ОБ АВТОРАХ

Должиков Вадим Владимирович — канд. техн. наук, доцент кафедры взрывного дела, http://orcid.org/0000-0001-8851-2913, Санкт-Петербургский горный университет, 199106, Санкт-Петербург, 21-я линия В.О., д. 2, Россия, e-mail: Dolzhikov_VV@pers. spmi.ru;

Рядинский Даниил Эдуардович — студент 6 курса специализации «Взрывное дело», Санкт-Петербургский горный университет, 199106, Санкт-Петербург, 21-я линия В.О., д. 2, Россия, e-mail: [email protected];

Яковлев Алексей Андреевич — студент 6 курса специализации «Взрывное дело», Санкт-Петербургский горный университет, 199106, Санкт-Петербург, 21-я линия В.О., д. 2, Россия, e-mail: [email protected].

Для контактов: Должиков Вадим Владимирович, e-mail: [email protected].

INFORMATION ABOUT THE AUTHORS

Dolzhikov V. V., Cand. Sci. (Eng.), Associate Professor at the Department of Explosives, http://orcid.org/0000-0001-8851-2913, Saint Petersburg Mining University, 199106, St. Petersburg, 21st Line, 2, Russia, e-mail: [email protected]; Ryadinsky D. E., 6th year student of the specialization "Explosives", Saint Petersburg Mining University, 199106, St. Petersburg, 21st Line, 2, Russia, e-mail: [email protected]; Yakovlev A. A., 6th year student of the specialization "Explosives", Saint Petersburg Mining University, 199106, St. Petersburg, 21st Line, 2, Russia, e-mail: [email protected]. The contact person: Vadim V. Dolzhikov, e-mail: [email protected].

Получена редакцией 14.01.2022; получена после рецензии 29.04.2022; принята к печати 10.05.2022. Received by the editors 14.01.2022; received after the review 29.04.2022; accepted for printing 10.05.2022.

i Надоели баннеры? Вы всегда можете отключить рекламу.