CC BY
13
- © Ю.А. Хохолов, Д.Е. Соловьев,
2013
УДК 622.45:536.244
Ю.А. Хохолов, Д.Е. Соловьев
ПРОГНОЗ ТЕПЛОВОГО РЕЖИМА РУДНИКА С УЧЕТОМ ДИНАМИКИ РАЗВИТИЯ ГОРНЫХ РАБОТ
Описан программный комплекс, позволяющий рассчитать тепловой режима и воз-духораспределение в вентиляционной сети р с учетом динамики развития горных работ. В качестве примера приведен расчет теплового режима проектируемого алмазодобывающего рудника АК «Алроса» на Накынском месторождении. Ключевые слова: криолитозона, тепловой режим, программный комплекс, рудники и шахты.
При строительстве и эксплуатации рудников и шахт в условиях криолитозоны, прежде всего, необходимо учитывать тепловой режим, который в определяющей степени влияет на эффективность и безопасность ведения горных работ. Как известно, на температурный режим выработок рудника влияние оказывает ряд факторов, из которых наибольшее значение имеют природные климатические условия района, наличие многолетнемерзлых пород и параметры вентиляции.
Прогноз теплового режима рудника с учетом этих факторов, а также динамики развития горных работ является весьма сложной задачей, требующей разработки специальных программ для ЭВМ на основе современных численных методов. Для её решения в лаборатории горной теплофизики ИГДС СО РАН разработан программный комплекс, позволяющий рассчитать воздухораспределение и тепловой режим сети горных выработок криолитозоны с учетом интенсивности теплообмена воздушной струи с окружающим массивом пород и фазовых переходов влаги в них, а также нестационарности сети горных выработок (изменчивости во времени числа выработок и их соединений).
В основу данного программного комплекса заложены следующие методики:
расчета воздухораспределения в вентиляционной сети с учетом напорных характеристик вентиляторов и геодезических отметок узлов сети, которая использует I и II законы Кирхгофа для вентиляционных сетей и основана на решении следующей системы уравнений [1, 2]:
т р
V а. - qJ Р = 0, 1 = 1,2,..., п -1, (1)
.=0 Рс
т р т р.
ЕЗ.. • —г. • q2 = V е.. -—к . + я V р.-Ьг.,
1. . ^ . в. о 1 1>
.=0 р .=1 рс
1= 1,2,..., т - п +1, (2)
ке = я= - ^,
^ (3)
ьв = а0 + а1Я1 + а2я2
расчета теплового режима в сети горных выработок криолитозоны, которая позволяет рассчитать температурный режим и ореолы протаивания во всех выработках сети, и учитывает нестационарное взаимовлияние процессов воздухораспределения и теплообмена в выработках сети с учетом фазовых переходов влаги в горных породах [3].
Алгоритм совместного решения задач заключается в следующем. Сна-
чала задается начальное распределение температур воздуха в выработках и в массиве вмещающих горных пород во всех ветвях вентиляционной системы. При заданных температурах воздуха решается задача воздухорас-пределения, т.е. определяются расходы воздуха в каждой выработке сети. Далее, исходя из топологической информации, определяется порядок расчета температуры воздуха по ветвям. Расчет начинается с воздухопо-дающих выработок сети. Температура воздуха на входе в выработку является уже известной величиной (т.е. входящие ветви в узел, с которого начинается выработка, уже прошли через процедуру расчета температуры воздуха). Согласно порядку проводится процедура расчета температуры воздуха в выработках сети. Температура воздуха на выходе из узла Тв рассчитывается с учетом расходов и температур воздуха всех входящих ветвей по формуле:
X Ъ'тч
Тв = ^-, (4)
№
д,, - соответственно, расход и температура воздуха на конце ]-й выработки, м3/с и °С; Р - множество входящих ветвей в заданный узел сети.
Далее, используя рассчитанное распределение температур воздуха в сети, уточняется перераспределение расхода воздуха в сети путем решения задачи воздухораспределения.
Таким образом, на каждом шаге по времени решаются последовательно 2 задачи: задача воздухораспределения и тепловая задача.
Методика совместного расчета вентиляционного и теплового режимов сети горных выработок модифицирована для случая, когда в сеть добавляются новые ветви или из сети удаляются отдельные ветви, т.е. когда
имеет место нестационарность сети горных выработок. Добавление новых ветвей в корне изменяет состояние сети, количество ветвей, узлов и независимых циклов, в то же время необходимо «запомнить» температурное поле вокруг выработок сети. Весь период работы рудника разбит на этапы, каждый из которых включает отработку нового блока, в сеть включаются новые выработки, а в конце этапа выключаются отдельные выработки.
Программа реализованна на языке Object Pascal в визуальной среде программирования Delphi (рис. 1) и предоставляет пользователю следующие возможности: интерактивный ввод теплофизических характеристик горных пород, информации о топологии вентиляционной сети, просмотр рассчитанных данных в табличном виде и их распечатку в соответствующем виде на принтере, а также графический просмотр всей сети в двух и трехмерном режиме. Дружественный интерфейс программы позволяет быстро моделировать вентиляционные сети.
В качестве примера рассмотрим тепловой режим алмазодобывающего рудника Накынского поля. В пред-проектных проработках на вскрытие и отработку алмазных трубок «Нюр-бинская» и «Ботуобинская» Накынкого поля подземным способом рассмотрено несколько вариантов, где отмечено, что на различных стадиях строительства и последующей эксплуатации рудников, температурный режим будет различным.
Рассмотрим один из вариантов расчета температурного режима единого подземного рудника, со вскрытием двух месторождений с одной основной промплощадки, располагаемой посредине между двух месторождений. Схема вентиляции рудника приведена на рис. 2. Всего на место-
Рис. 1. Рабочее окно программы
рождениях тр. «Ботуобинская» и тр. «Нюрбинская» необходимо отработать 12 блоков, поэтому весь период работы единого рудника разбит на 12 этапов. В начале каждого этапа появляются новые ветви и узлы вентиляционной сети, а в конце этапа закрываются отдельные ветви, что приводит к прекращению движения воздуха по выработкам данного блока.
Проветривание горных работ будет осуществляться по центральной схеме проветривания. Свежий воздух (одновременно для 2-х месторождений) будет подаваться по клетевому стволу (КС) основной промплощадки рудника. Часть, подаваемого в рудник свежего воздуха будет распределяться в горные выработки тр. «Ботуобин-ская», с выдачей исходящей струи по скиповому стволу (СС) основной промплощадки рудника, а другая часть будет подаваться в горные выработки тр. «Нюрбинская», с выдачей исходящей струи через вспомогательный скиповой ствол (ВСС).
Вентиляционная сеть выработок рудника состоит из 111 ветвей и 64 узлов, следовательно, количество независимых циклов равно к=111-64+1=48. Выходы на поверхность соединены фиктивными ветвями с узлом 1 с нулевыми аэродинамическими сопротивлениями (ветви 13 и 27). Температура атмосферного воздуха задается периодической функцией по формуле:
Т - Т П720(М - 1) + г +1800]
Тв = -л ■ sin—------] +
2
4320
Т - Т
з_л
2
(5)
где Т3 - минимальная зимняя температура, 0С; Тл - максимальная летняя температура, 0С; М - порядковый номер месяца, М=1,2,...,12.
По проекту в руднике будет обеспечиваться постоянный положительный температурный режим (не менее +2 0С на воздухоподаче ствола КС), поэтому, как только наружная температура воздуха стане ниже +2 0С,
СС
|59 '76
Рис. 2. Схема вентиляции единого подземного рудника
включается калорифер, который поддержит температуру поступающего воздуха на заданном уровне.
Были так же рассчитаны глубины протаивания пород вокруг выработок, температуры воздуха в каждой выработке. Отдельные результаты расчета приведены на рис. 3 - 5. Поскольку в выработках поддерживается положительный температурный режим и расходы воздуха достаточно большие, происходит прогрессирующее протаивание окружающих мерзлых горных пород. Значения глубины протаивания зависят от периода эксплуатации выработок, так за 23 года максимальная глубина протаивания вокруг КС на отметке (+254 -+70 м) составляет 11 м, на отметке (-230 -330 м) - 10,5 м, на отметке (-330 -430 м) - 9 м, на отметке (-430 - 530 м) - 7,5 м. Наиболее интенсивное про-таивание происходит в первые годы. Значения глубины протаивания во-
круг СС и ВСС ниже, чем глубины протаивания вокруг КС.
При вскрытии новых горизонтов происходит охлаждение вентиляционного воздуха из-за теплообмена с окружающими горными породами, а поскольку исходящая струя воздуха выдается на поверхность через СС и ВСС, то это приводит к замедлению темпов протаивания пород вокруг них и характерному волнообразному виду графиков глубины протаивания на различных горизонтах (см. рис. 4 и 5).
На рис. 6 приведены расчетные значения изменения температур воздуха в выработках рудника на 5 год эксплуатации. Из графиков видно, что в зимний период температура воздуха равна примерно +2 °С, т.к. стабильность температуры обеспечивается калорифером, а в летний период за счет поступления атмосферного тепла воздух нагревается до +18 °С.
—•— отм.+254 +70
—■— отм.+70 -30
- - - х- - ■ отм.-30 -130
—ж— отм.-130 -230
отм.-230 -330
—1— отм.-330 -430
-отм.-430 -530
Рис. 3. Динамика глубины протаивания КС
Время, год
отм.+254 +170
—■— отм. + 170 +70
-■■х-- ■ отм.+70 -30
—ж— отм.-30 -130
—•— отм.-130 -230
—1— отм.-230 -330
-отм.-330 -430
Рис. 4. Динамика глубины протаивания СС
Время, год
—♦— -отм -254 -30
—■— отм -30 -130
X ■ отм -130 -230
—ж— отм -230 -330
—•— отм -330 -430
Рис. 5. Динамика глубины протаивания ВСС
Рис. 6. Динамика температур воздуха в выработках в 5 год эксплуатации
Поддержание круглогодичного положительного температурного режима в выработках рудника позволяет улучшить санитарно-гигиенические условия труда горнорабочих, а также обеспечить приемлемые условия для работы машин механизмов. В тоже время положительная температура воздуха в горных выработках приводит к интенсивному протаиванию мерзлых пород вокруг выработок, что отрицательно сказывается на их устойчивости.
Следует отметить, что при нерегулируемом тепловом режиме рудника в выработках температура воздуха в зимний период будет отрицательной, и протаивание горных пород происходит только в устьевых частях воздухоподающих выработок в летний период. С позиций геомеханики их мерзлое со-
1. Абрамов Ф.А. Воздухораспределе-ние в вентиляционных сетях шахт / Ф.А. Абрамов, Р.Б. Тян, В.Я. Потемкин. - Киев: Наукова думка, 1971. - 136 с.
2. Ушаков К.З. Аэрология горных предприятий / К.З. Ушаков, А.С. Бурчаков, Л.А. Пучков, И.И. Медведев; под ред.
стояние повышает устойчивость, однако отрицательная температура воздуха ухудшает условия труда рабочих и создает определенные трудности при эксплуатации машин и механизмов.
Результаты проведенных исследований позволили сделать прогноз теплового режима рудника, который необходим для оценки геомеханического состояния горных выработок и обоснования основных параметров производственных процессов. Таким образом, разработанный программный комплекс обеспечивает объективную оценку состояния вентиляционного и теплового режимов во всех выработках сети рудника с учетом ее нестационарности, что важно для повышения эффективности и безопасности ведения горных работ.
- СПИСОК ЛИТЕРАТУРЫ
К.З.Ушакова.-Изд. 3-е, перераб. и доп. - М.: Недра, 1987. - 421 с.
3. Хохолов Ю.А. Совместное решение задач воздухораспределения и теплового режима в сети горных выработок криолитозоны / Ю.А. Хохолов // Горный информ.-аналит. бюллетень. - 2003. - № 7. - С. 70-72. ТШ
КОРОТКО ОБ АВТОРАХ -
Хохолов Юрий Аркадьевич - доктор технических наук, ведущий научный сотрудник, khokholov@igds.ysn.ru,
Соловьев Дмитрий Егорович - кандидат технических наук, младший научный сотрудник, Институт горного дела Севера им. Н.В.Черского СО РАН, igds@ysn.ru
А
