Научный журнал

Успехи современного естествознания

ISSN 1681-7494

"Перечень" ВАК

ИФ РИНЦ = 0,775

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Синицын В.А. 1 Меньшиков П.В. 1 Шеменев В.Г. 1

---

1 ФГБУН «Институт горного дела Уральского отделения Российской академии наук»

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

Статья в формате PDF

0 KB

скважинный заряд

эмульсионные взрывчатые вещества

гидростатическое давление

глубина скважины

плотность

математическая модель

горная порода

степень аэрации

масса заряда

1. Дерибас А.А., Решетняк А. Исследование влияния добавок в виде полых микросфер на скорость детонации различных конденсированных взрывчатых веществ // Физические проблемы разрушения горных пород: сб. тр. Третьей Междунар. науч. конф. (9–14 сент. 2000 г.). Абаза. – Новосибирск: Наука, 2003. – С. 151–154.

2. Жученко Е. И., Иоффе В. Б., Кукиб Б. Н., Сундуков И. Ю., Оверченко М. Н. Заряжание глубоких скважин эмульсионными ВВ сибиритами на разрезах Кузбасса // Физические проблемы разрушения горных пород: сб. тр. Третьей Междунар. науч. конф. (9–14 семг. 2000 г.). Абаза. – Новосибирск: Наука, 2003. – С. 154–158.

3. Промышленные взрывчатые вещества / Л.В. Дубнов и др. – М.: Недра, 1988. – 358 с.

4. Синицын В.А. Повышение эффективности подготовки горной массы на карьерах с применением ВВ на основе обратных эмульсий: дис. ... канд. техн. наук. – Екатеринбург, 2007. – 145 с.

5. Соснин В.А., Колганов Е.В. Исследование детонационных процессов эмульсионных взрывчатых веществ // Взрывное дело. – 2004. – № 94/5. – С. 181–195.

6. Ударно-волновая чувствительность модельного ЭМ, сенсибилизированного стеклянными микросферами / Б.А. Хасаинов и др. // 12-й симпоз. по горению и взрыву. – Черноголовка, 2000. – С. 15–21.

7. Чижов И.А. Разработка методики расчета радиуса разрушения от взрыва скважинного заряда ВВ при дроблении крепких крупноблочных пород на карьерах: дис. … канд. техн. наук. – Екатеринбург, 1993. – 153 с.

8. Шер Е.Н., Александров П.И. Оценка влияния работы забойки и проникновения газов в среду на разрушающее действие скважинного заряда // Физические проблемы разрушения горных пород: сб. тр. Третьей Междунар. науч. конф. (9–14 септ. 2000 г.). Абаза. – Новосибирск: Наука, 2003. – С. 117–122.

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

Вследствие ухудшения условий открытой разработки полезных ископаемых – увеличение крепости и обводненности пород с глубиной в карьерах, обеспечение необходимого качества дробления горных пород без существенного увеличения затрат возможно только за счет интенсификации буровзрывных работ (БВР), одним из основных направлений которой является оптимизация параметров скважинной отбойки. Имеющиеся на сегодняшний день расчетные методики для их определения обладают рядом недостатков и требуют усовершенствования в части более полного учета особенностей взрывного процесса в среде. Следовательно, совершенствование методов определения параметров БВР является актуальной научной задачей и имеет большое практическое значение [7].

Рассмотрим пузырек газа в нитроните Э-70. Давление в пузырьке (рис. 1):

(1)

где δ – коэффициент поверхностного натяжения; r – радиус пузырька; Р – давление в окружающем пузырьке нитронита Э-70.

Рис. 1. Давление в пузырьке

Рассмотрим скважинный заряд нитронита Э-70 (рис. 2). Введем обозначения: Lз – длина заряда нитронита Э-70; Ро – давление на верхний торец заряда; ρс – плотность эмульсии (не аэрированной); ρо – плотность ВВ при заряжании; Р(х) – давление в заряде в сечении (на глубине х относительно верхнего торца заряда); ρ(х) – плотность в заряде в сечении х.

Рис. 2. Скважинный заряд нитронита Э-70

Эмульсия состоит из двух несмешивающихся жидкостей. Поэтому максимальную плотность можно рассчитать через удельные объемы, занимаемые каждой жидкостью в эмульсии [1, 2, 5].

Пусть Vc = 1/ρс – удельный объем эмульсии; Vг(х) – удельный объем газа в сечении х. Тогда

ρ(х) = (Vc + Vг(х))–1, (2)

где Vc = const.

Количество пузырьков N в любом сечении х постоянно ввиду большой вязкости нитронита Э-70 (нет распределения по х вследствие диффузного движения пузырька). Тогда

(3)

(4)

где Vr0 – удельный объем газовой фазы при заряжании; r0 – начальный радиус пузырька при заряжании. Тогда из (3) и (4)

(5)

где rx – радиус пузырька в сечении х.

Подставляем (5) в (2):

или

(6)

Здесь , Z0 – степень аэрации нитронита при заряжании. В других обозначениях . Учитывая, что 1/Vc = ρc, выражение (6) можно представить:

(7)

Рассмотрим теперь элемент объема заряда на глубине х толщины Δх (рис. 3).

Из условия равновесия данного элемента следует:

(8)

где S – площадь сечения заряда. Здесь учтено, что при достаточно малом Δх можно считать плотность в данном элементе постоянной и равной Δ(х). Тогда (8) имеет вид При ρх > 0 получаем . Подставляя (7), имеем

(9)

Пусть теперь сжатие пузырьков происходит адиабатически (распределение ρ(х) устанавливается в заряде до наступления существенного остывания нитронита). Тогда

(10)

где Pп.0 – давление в пузырьке при заряжании; Pп.x – давление в пузырьке, расположенном в сечении х; Vп.0 – объем газа пузырька при заряжании; Vп.x – объем газа пузырька, находящегося в сечении x; γ = 1,25 – показатель адиабаты.

Рис. 3. Элемент объема заряда

С учетом (1) формулу (10) можно записать в виде

или

(11)

где Рат – атмосферное давление.

Уравнения (9) и (11) полностью описывают поставленную задачу, если их дополнить начальными условиями: при х = 0

Рх = Р0; rx = r0. (12)

Введем обозначение α = rx/r0. Тогда получается система уравнений (9), (11), (12)

(13)

Из (13) при заданных g, ρc, ρ0, r0, Р0, Рат, δ можно получить Р(х), α(х) в любом сечении х, после чего по (7) определяется и ρ(х). Если пренебречь поверхностным натяжением (δ = 0), то (13) принимает вид

В этом случае α3(х) = P(x)1/γ.

Из общего вида системы (13) видно, что она не имеет аналитического решения, следовательно, требует численного решения. Опишем его алгоритм.

Дифференциальное уравнение (13) совместно с начальными условиями является задачей Коши, которую решаем наиболее простым методом – разностным методом Эйлера (рис. 4) по причине достаточной гладкости и монотонности функций P(x) и ρ(х).

Рис. 4. Разностный метод Эйлера

Алгоритм решения:

(14)

Второе уравнение системы (13) является трансцендентным. Его решаем итерационным методом – методом бисекции:

(15)

После этого на каждом шаге по х определяем плотность:

(16)

Уравнения (14)–(16) решаются совместно. Для этого данный алгоритм реализован в виде программы PUZ на ЭВМ.

Входными данными являются: ρc, ρ0, r0, Р0, Рат, δ. Выходными х, P(x), ρ(х), ρ(х), (х), δ, где δ(х) = (ρ(х) – ρ0/ρ0.

Здесь Р0 берется в общем виде и имеет вид

Р0 = Рат + ρзабgLзаб,

где ρзаб – плотность забойки (например, воды); Lзаб – длина забойки.

Если Lзаб = 0, то забойка отсутствует.

Оценочные результаты приведены в таблице.

Таким образом, можно сделать следующие выводы:

1. При r0 = 0,1 – 0,01 мм P(x) и ρ(х) практически не зависят от σ (при σ < σмаз).

2. При r0 = 0,01 мм P(x) и ρ(х) зависят от r0 и σ.

3. Наличие забойки существенно сказывается на P(x) и ρ(х).

4. На глубине больше 10 м ρ(х) > 1,3 г/см3.

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

Высота водяной забойки равна 7 м. Это значение является максимальным для скважины глубиной до 35 м [2]. В случае, если высота забойки будет меньше, уменьшится гидростатическое давление на эмульсионное ВВ в нижней части скважины, что повысит его детонационную способность. Максимальная предельная плотность нитронита по ТУ 7276-003-58995878-2004 принята равной 1,25 г/см3 [3, 4, 6].

Плотность нитронита Э-70 на различной глубине скважины рассчитывается следующим образом. На слой эмульсионного ВВ, находящегося в скважине на глубине Н, действуют следующие силы: РА – атмосферное давление; РВ – гидростатическое давление столба воды высотой Δhв; Р1 – гидростатическое давление столба эмульсионного ВВ высотой Δh1) [3, 4, 8].

Атмосферное давление РА принимаем равным 980 Па. Гидростатическое давление столба воды высотой Δhв = 7 м; Рв ~ 686 Па.

Гидростатическое давление столба ЭВВ высотой Δh можно выразить как ΔР = ρЭВВ•Δh (кг/см2), т.е. значение ΔР будет зависеть от высоты слоя Δh и плотности эмульсионного ВВ, которая непрерывно изменяется по высоте этого слоя.

Оценочные результаты

Примечание. ρc = 1,35 г/см3, ρ0 = 1,25 г/см3.

Рис. 5. Изменение плотности нитронита Э-70 в колонке заряда по глубине скважины при различной исходной плотности

Для того чтобы рассчитать плотность на определенной глубине скважины, зададимся каким-либо фиксированным значением Р1. Тогда давление, под которым будет находиться слой эмульсионного ВВ на глубине Н, РН1 = РА + РВ + Р1; РH1 ≈ 1666 Па.

Увеличивая последовательно высоту заряда на слой, равный 5 м, рассчитывая аналогичным образом плотность ЭВВ на нижней границе этого слоя, можно проследить изменение плотности эмульсионного ВВ, размещенного на различной глубине скважины. В качестве забойки принят столб воды высотой 7 м. Расчеты сделаны для нитронита Э-70 с различной исходной плотностью (0,8; 0,9; 1,0; 1,1; 1,2; 1,3) г/см3, а их результаты представлены на рис. 5.

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

Например, необходимо зарядить обводненные скважины глубиной 20 м. На рис. 5 на оси абсцисс находим глубину скважин, равную 20 м, и идем вверх по вертикали до пересечения с прямой, отвечающей критической плотности, которая для зарядов в прочной оболочке диаметром 250 мм равняется 1,26 г/см3. Таким образом, заряжать скважины глубиной 20 м следует с исходной плотностью, указанной цифрами у кривых, лежащих ниже горизонтальной прямой критической плотности и левее вертикальной прямой, соответствующей глубине скважины 20 м. Согласно рис. 5 для рассматриваемого примера начальная (исходная) плотность нитронитов Э-70 не должна превышать 1,1 г/см3 [4].

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

Библиографическая ссылка