Гидромониторная струя и модель разрушения угля

© A.B. Бартышев, 2012

УДК 622.236.52 A.B. Бартышев

ГИДРОМОНИТОРНАЯ СТРУЯ И МОДЕЛЬ РАЗРУШЕНИЯ УГЛЯ

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

Ключевые слова: гидромониторная струя, динамическое давление, критерий разрушаемости угольного массива, коэффициент осевого динамического давления, эффективная длина струи, радиус эффективного ядра.

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

для определения производительности гидравлического разрушения угля по длине струи в различных горно-геологических условиях.

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

Наиболее достоверные результаты будут получены, если эмпирическую постоянную Рт1п определить непосредственно из опыта. Для этого необходимо в забое установить гидромонитор с известными гидро-

Работа выполнялась в рамках государственного контракта Минобрнауки России №16.515.11.5073 от 06.10.2011 года.

динамическими параметрами струи и отрабатывать заходку до тех пор, пока разрушение угля не прекратится. После этого измеряются расстояние 1Э от насадка до угольного массива. По известному закону изменения осевого динамического давления Рт=Рт(1) (далее по тексту формула (3)) можно определить давление Ртп как значение функции при 1=1 Э, то есть Ртт=Рт(1э).

Если нет возможности провести такой эксперимент, то значение Рт.п рассчитывается по условному пределу прочности, определяемому как среднее геометрическое предела прочности на сжатие стс и предела прочности на растяжение стр образца стандартного размера

4

= 3

ст ст

Р с

где Яу — условный предел прочности, МПа; стр и стс — пределы прочности на растяжение и на сжатие, МПа.

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

На рис.1 представлена схема гидромониторной струи. Эффективная часть струи гидромонитора (1Э) содержит два участка: начальный (1Н), и основной (10). На начальном участке осевое динамическое давление Рт постоянно и равно подводимому к гидромонитору Ро (ядро постоянных давлений). На основном участке осевое динамическое давление струи уменьшается от Р0 до критерия разрушаемости Ртп.

Основной параметр «разрушающей» гидромониторной струи Р (1, г) — динамическое давление (интегральная характеристика, учитывающая скорость и плотность одновременно), определяющее способность струи разрушать уголь и измеряемая с помощью гидравлического датчика давления (трубки Пито).

Обработка результатов многочисленных измерений параметров гидромониторной струи на различных расстояниях показала, что величина динамического давления на основном участке, в сечении на расстоянии 1 от насадка, распределяется по закону, соответствующему решению уравнений пограничного слоя при рассмотрении турбулентной осесимметрич-ной струи-источника несжимаемой жидкости, допускающей автомодельные решения [2] Рт

Р(г, I) =

(1+

Г 2 Г \2

(в-1)2

10

+Ргп

2 4 6 8 1

Рис. 2. Изменение осевого динамического по длине струи

где е — эмпирическая постоянная, определяемая из опыта и зависящая от характеристик струеформирующих узлов гидромонитора, а Рт — осевое динамическое давление на расстоянии 1 от насадка.

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

Рт

( + ф )2

где д — постоянная, определяемая из опыта, а 1ф — расстояние от фокуса до насадка.

Объединяя выражения (1) и (3), получим закон изменения динамического давления на основном участке гидромониторной струи как функцию двух переменных г и 1 с тремя параметрами д, е и 1ф

Р (г, I) =--

1 +

( г 2 ^

е1

/у

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

В [6] автором было введено понятие «коэффициента осевого динамического давления» К,0 (КОДД), который вычисляется по результату измерения осевого динамического давления Ртк на расстоянии 1к = 200ё0, по формуле Р

К, = -

Для гидромониторных струй значения эмпирических констант зависят от конструкции гидромонитора и ка-

Достоинством предложенной расчетной модели является то, что один этот эмпирический коэффициент К, , (учитывая закон сохранения импульса, ,0 — диаметр насадка и Р0 — подводимое давление) позволяет определить все динамические параметры струи. Из закона сохранения импульса, поскольку на выходе из насадка

^ = Р0Б0 = 14 пР00> (где Бо — площадь насадка) следует

4 пРо02 =%Рт (е1)2 или

о

1 2 3

Рис. 3. Распределение динамического давления в сечении струи

е =

21к\Ртк 21кУК 400^'

Найденное значение е численно равно тангенсу угла расширения струи и поэтому позволяет количественно оценить компактность струи. Чем меньше е, тем меньше угол расширения струи, то есть струя распадается медленнее. Учитывая это, будем называть е «коэффициентом компактности». Зная тангенс угла расширения струи, легко найти расстояние от среза насадка до фокуса струи

ф=2е=кЩ, = 200^

Подставляя (7) в (3) и, учитывая определение КОДД, найдем 8:

8 = Кр + 200,0 ] = = кр 12 (( +1)2 = = КР (200,0 )2 ((К +1)2

Так как 8 прямо пропорциональна осевому динамическому давлению (3), то вполне естественно назвать ее «коэффициентом дальнобойности». Чем больше «коэффициент дальнобойно-

сти», тем дальше струя сохраняет разрушительную способность, тем больше ее эффективная длина. На рис.3 приведено распределение динамического давления в сечении струи на основном участке для двух гидромониторов с различными КОДД.

Отметим, что «коэффициент компактности» и «коэффициент дальнобойности» удачно дополняют друг друга, задавая два измерения струи: продольное и поперечное. Таким образом, для полного описания скалярного поля динамического давления струи на основном участке достаточно знать начальные условия истечения — подводимое давление и диаметр насадка, а также измерить осевое (максимальное) динамическое давление на расстоянии 1к = 200ё0 (примерно 3 — 6 м) для определения КОДД.

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

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

Длину начального участка

можно определить, если подставит в

(3) Рт = Р>:

I. = .1Р - 1Ф = 1К = 200,0К,. (9) ' р

1

чим из (4), полагая Р=РГ и выразив радиус через расстояние от насадка:

г = е1

-1.

Рис. 4. Границы эффективных ядер двух ниторов

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

Эффективная длина струи вычисляется из (3), если положить осевое динамическое давление равным критерию разрушаемости Рт = Ртп. Разрешая это уравнение относительно 1 с учетом значений д и 1ф из (7) и (8), получим

гидромо-

1Э = 200^

.

Если положить РГ=Ртп, то получим границу эффективного ядра струи,

которое разрушает массив. Зависимость радиуса эффективного ядра от 1 имеет вид:

Я = е1

1

И +

-1 =.

= вЛ

- 1Ф) -1

ф

1+и

= е1

1 + и

Исследуя эту функцию, найдем ее максимальное значение Яэт, которое равно

при 1 = 2.

Я

Нетрудно заметить, что дальность действия струи пропорциональна коэффициенту 3 и поэтому оправдано его название — «коэффициент дальнобойности струи».

Границу струи можно определить как поверхность уровня скалярного поля динамического давления, на которой оно имеет постоянное значение РГ (например, РГ равно атмосферному давлению). Граница определяет пределы интегрирования при вычислении потока кинетической энергии и других интегральных величин. Уравнение границы полу-

На рис. 4 изображены эффективные ядра двух гидромониторов с различными коэффициентами осевого динамического давления.

Используя полученные ранее зависимости, радиус эффективного ядра струи может быть выражен через КОДД и входные данные:

Я =-х

400^/К,

200 0 К ((+1) ). -1 V / Ш1П

-1

1+200,0

Сила удара гидромониторной струи на основном участке на расстоянии 1 от насадка на круговую площадку радиуса Я с центром, расположенным на оси струи, равна:

х

п п

Р 1, Я) = 2п| Р  гСг = 2п|-

тгСг =.

1 +

(в1)2

с(Е/)2 Я2

Р.

(в/)2 + Я2

Если в качестве Я возьмем Яэ, то получим силу удара эффективной части струи Рэ(1), которая служит интегральной характеристикой активного ядра струи

(( + 1)с0Л/Рр~

Р (/) = п

1 + и

8004К

((э - 1).

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

Методики расчета параметров гидравлического разрушения, используемые в настоящее время [4,5,8], базируются на результатах длительных экспериментов с хроно-метражными наблюдениями, а получаемые эмпирические зависимости описывают только средние величины, справедливые в конкретных горногеологических условиях. При этом входными параметрами являются только подводимое давление и расход, без учета качества струи (гидромонитора): компактности и дальнобойности. Получение зависимости производительности разрушения от параметров струи позволит более точно планировать горные работы, определять оптимальные параметры выработок в зависимости от имеющегося оборудования, управлять про-

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

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

Анализ теоретических и экспериментальных исследований взаимодействия гидромониторной струи с угольным массивом позволил установить, что производительность гидравлического разрушения прямо пропорциональна осевому динамическому давлению струи гидромонитора Рт(1) (динамическая характеристика струи) и площади эффективного ядра Бэ(1) (связь струи с массивом).

Пр

где С — коэффициент, характерный для полуэмпирических моделей, связывающий расчетную производительность с реальной и определяемый по результатам шахтных экспериментов. Заметим, что эта функция обнуляется на расстояниях, где осевое динамическое давление становится меньше критерия гидравлической разрушаемости (Рт < РтД так как площадь

эффективного ядра обращается в ноль и разрушение прекращается.

Производительность гидроразрушения по уравнению (16) не зависит явно от общего расхода, как в известных формулах, используемых ранее [1,3,5], так как учитывает только разрушающую часть струи.

Учитывая кусочно-непрерывный характер зависимостей осевого динамического давления Рт=Рт(1) и площади эффективного ядра 5Э=БЭ(1), предложенные в инженерной модели гидромониторной струи функция, задающая расчетную производительность уравнения (16) разобьется на три части:

П „ 9 1-.

4'

Про (I) = ^ Р082

1 + 2Л

при 1 е[0,1н ],

I

1

Пр (1) = -8е2--2

1р () 4 (1э + 21,)(1 + 1Л )2

при 1 е

1.

1

Прр (I) = п8е2 П 1э - 1

(1+и

при 1 е

1

Переход от точечного разрушения угольного массива по длине струи к разрушению некоторого объема заданной конфигурации основывается на механизме элементарного взаимодействия струи и угольного массива. Процесс разрушения опишем дифференциальным уравнением dV=Пр(1)dt, где dV — элемент объема угольного массива, dt — элемент времени, а Пр(1) — производительность, определяемая формулой (16).

Интегрируя уравнение, получим время Т выемки объема угля V, при известных параметрах отрабатываемой заходки

т=Ш

dV 1

dV

ется по формуле (17), в 02 — по формуле (18) и в Qз — по формуле (19), то интеграл (20) разобьется на три и получим расчетное время гидравлического разрушения угля в выработке

Тр+

+

#

01 ПРр

dV

,,,

Так как функция Пр(1) состоит из трех частей, то и область интегрирования разобьется на три непересекающиеся множества 0 = 01 и 02 и 03, причем в производительность изменя-

03 ПРр

Если известно реальное время отработки заходки Т3, то можно определить коэффициент пропорциональности С из соотношения С • Тз = Тр . Постоянная С определяет связь между расчетной производительностью и реальной производительностью в конкретных горно-геологических условиях.

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

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

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

Сечение выработки

Рис. 5. Схема забоя с продольными границами в виде концентрических сферических поверхностей

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

Если заходка представляет собой прямоугольный параллелепипед, то предполагая, что начальная и конечная границы в направлении прохождения выработки, представляют собой части сферических поверхностей <1 и <2 с радиусами 1+1 и 1с+12 (11 и 12 — расстояния от среза насадка до и < 2), соответственно, вырезанные из сферы прямоугольным параллелепипедом (рис. 5).

Если на расстоянии от среза насадка 1 в пределах эффективного участка струи элемент объема заходки представить как

64 = Б1) с11,

то при известной функции Б(1), описывающей изменение площади части сферической поверхности, вырезанной из сферы с радиуса 1с+1 прямоугольным параллелепипедом с сече-

нием заходки, тройной интеграл из (20) можно заменить одномерным

я

= 2 Б1) 61 Пр ~ ! Пр( 1) '

В реальных условиях можно считать, что площадь поверхности забоя на расстояниях, при которых производится гидроотбойка (1,5 — 10 м), изменяется линейно:

Б1) = А1 + В, 1 е[ 11,12 ].

Тогда зная начальную (3(11)=31) и конечную площадь (3{12)=Б2) сферической поверхности можно вычислить постоянные А и В:

А = Б Б

2 -1

в=

(Б - Б ) ¡1

4 -11

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

Б 1)61 1 ПрЛ( 1

Б 1)61

Б 1) 61

ПРэ2( 1

^Прэ3

2

I

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

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

Выводы

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

Эффективная длина струи и радиус эффективного ядра могут быть использованы при моделировании процесса разрушения угольного массива и вычисления производительности гидромонитора с известными характеристиками.

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

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

1. Абрамович Г.Н. Теория турбулентных струй. — М.: Госгортехиздат, 1960.

2. Вулис Л.А., Кашкарова В.П. Теория струй вязкой жидкости. — М.: Наука, 1965.

3. Бай Ши-и. Турбулентные струи. — М.: Мир, 1964.

4. Цяпко Н.Ф., Чапка А.М. Гидроотбойка угля на подземных работах. — М. Гос-гортехиздат, 1960. — 263 с.

5. Шавловский С. С. Основы динамики струй при разрушении горного массива. — М.: Наука, 1979.

6. Бартышев А.В. О критериях качества струи // Сб. научных трудов ВНИИгидроугля / Совершенствование технологии и оборудования при гидравлической добыче угля, — Новокузнецк, 1990. — С. 86—93.

7. Мерзляков В.Г., Бафталовский В.Е. Физико-технические основы гидроструйных технологий в горном производстве. — М.: ФГУП, 2004. — 645 с.

8. Кривченко A.A. Рекомендуемые методы проектирования и расчета гидроотбойки. — Донецк: Доклад на совещании работников проектных организаций угольной промышленности, 1964.

9. Методика расчета гидравлической выемки угля. — Новокузнецк: ВНИИгидроуголь, 1969.

10. Теодорович Б.А., Цяпко Н.Ф., Гефт В.Б. Методы расчета производительности гидромонитора по базовой производительности гидровыемки. //Вопросы гидравлической добычи угля. Сб. 16. — Новокузнецк: ВНИИгидроуголь, 1969. ШНЗ

КОРОТКО ОБ АВТОРЕ -

Бартышев Александр Васильевич — кандидат технических наук, доцент, е-шаП: bartyshev@mail.ru,

Новокузнецкий филиал — институт Кемеровского государственного университета.