Calculation of surface air leaks at potash mine
- Authors: Nikolaev A.V.1, Alymenko N.I.2, Sadykov R.I.1
- Affiliations:
- Perm National Research Polytechnic University
- Mining Institute, Ural Branch of the Russian Academy of Sciences, Perm, Russian Federation
- Issue: No 5 (2012)
- Pages: 115-121
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/geo/article/view/1077
- Cite item
Abstract
When ventilating mines the suction method as all the mines Verkhnekamsk potassium-magnesium salts (VKMKS) in the surface complex of the main ventilation system (HLG) formed air leakage. The value of surface air leakage at operating mines determined by conducting air-depression shooting. For projected mines exist coefficients, which takes the value of projected air leaks, in fact, significantly different from those obtained during the measurement values. This article provides an algorithm for calculating the value of surface air leaks that occur when working at HLG potash mines. The proposed algorithm for the calculation of surface air leakage is the known and the author earlier formulas. For the unknown quantity on which in the future can calculate the surface leakage of air at a certain mode HLG adopted by the fan discharge at the intersection of Canal HLG with ventilation shafts. Equating each other formulas that determine the vacuum created by the HLG on the site surface complex - Channel HLG and the area air supply trunks - the underground part of the mine - ventilation shaft - Channel HLG, and implemented a number of reforms, was obtained the required mathematical relationship. Algorithm used in formulas to determine the amount of surface air leaks under changing climatic parameters (atmospheric pressure, air temperature, relative humidity, etc.). According to the obtained mathematical relationship mode is selected HLG. Using the obtained in the algorithm, it is possible to determine the amount of surface leakage current on the air as well as on the projected potash mines.
Full Text
Введение На рудниках Верхнекамского месторождения калийно-магниевых солей (ВКМКС) применяется всасывающая схема проветривания, т.е. воздух в рудник подается по воздухоподающим стволам 1 за счет разряжения, создаваемого главной вентиляторной установкой (ГВУ) 2, расположенной на вентиляционном стволе 3 (рис. 1). Кроме воздуха, засасываемого в рудник Qр, при подобной схеме проветривания в канал ГВУ 4 также засасывается наружный воздух через устье вентиляционного ствола 5, который, смешиваясь с общим потоком, будет снова выбрасываться в атмосферу. Поэтому помимо полезной работы (подачи свежего воздуха в рудник) ГВУ будет совершать бесполезную работу (обусловленную внешними утечками воздуха). При проектировании рудников величина поверхностных утечек воздуха согласно [1, 2, 3 и др.] устанавливается равной 5–25 % (0,05–0,25 о.е. – так называемый коэффициент поверхностных утечек Kут, зависящий от типа ствола). Однако на практике согласно [4] утечки воздуха на рудниках ВКМКС колеблются от 9 до 51 %. Разность закладываемого в проект и реального Kут будет сказываться на работе ГВУ, так как реальный режим ее работы может не соответствовать проектному. Подобное несоответствие наблюдается и на других рудниках [5, 6]. Ранее предлагалась методика определения величины поверхностных утечек воздуха только на действующих рудниках при известной производительности ГВУ [7, 8]. В связи с этим требуется разработать алгоритм расчета, позволяющий определить величину поверхностных утечек воздуха уже на стадии проектирования рудника. Рис. 1. Упрощенная схема (центральная) рудника, работающего по всасывающему способу проветривания Величина поверхностных утечек воздуха может быть определена как разность производительности ГВУ QВ и объема воздуха, поступающего в рудник Qр: . (1) Величина Qр определяется на стадии проектирования рудника исходя из общего объема воздуха, который необходимо подавать на все добычные участки. В связи с этим величина поверхностных утечек воздуха будет зависеть от производительности ГВУ, которая, в свою очередь, будет зависеть от многих факторов, основным из которых является величина разряжения, создаваемого ГВУ Pк в вентиляционном стволе (рис. 2). Определение барометрического давления в околоствольном дворе вентиляционного ствола и разрежения ГВУ Разрежение ГВУ может быть определено с помощью двух формул: согласно [9] (2) согласно [10] (3) Рис. 2. Воздушные потоки в вентиляционном стволе где Pa – атмосферное давление воздуха, мм рт. ст.; K – коэффициент пересчета из Па в мм рт. ст., равный 0,0075; hВ – давление, создаваемое вентилятором, Па; Rкан – аэродинамическое сопротивление канала ГВУ, Н∙с2/м8; Pод.вент – барометрическое давление в околоствольном дворе вентиляционного ствола, мм рт. ст.; ρср.вент – средняя плотность воздуха в вентиляционном стволе, кг/м3; Hвент – участок вентиляционного ствола от сопряжения с околоствольным двором до сопряжения с вентиляционным каналом ГВУ (рис. 2), м; Rвент – аэродинамическое сопротивление Hвент, Н∙с2/м8. Величину барометрического давления в околоствольном дворе вентиляционного ствола согласно [10] можно найти как (4) где Pод.возд – барометрическое давление в околоствольном дворе воздухоподающего (воздухоподающих) ствола (стволов), мм рт. ст; Rпчр – аэродинамическое сопротивление подземной части рудника, Н∙с2/м8. Рис. 3. Пример расчетной схемы Определение аэродинамического сопротивления подземной части рудника Аэродинамическое сопротивление подземной части рудника можно определить при составлении его расчетной схемы и дальнейшем ее расчете, который проводится следующим образом. Представим, что рудник отрабатывается блоками по схеме, приведенной на рис. 3. Для расчета воздухораспределения в приведенной сети первоначально необходимо ее упростить. Рассмотрим блок 1. В нем выработки (ветви) 1, 2, 3 и 4, 5, 6 соединены параллельно. Следовательно, можно найти их эквивалентное аэродинамическое сопротивление по формулам [2] (5) где Ri – аэродинамическое сопротивление i-й выработки в расчетной схеме, Н∙с2/м8. В результате получим схему, приведенную на рис. 4, а, в которой ветви 1, 2, 3 и 4, 5, 6 заменены эквивалентными. Ветви, соединенные последовательно, находятся как их алгебраическая сумма. Сложив значения аэродинамических сопротивлений ветвей 9, 10 и 11, 12, а затем, определив их общее аэродинамическое сопротивление по формуле (5), полу- Рис. 4. Пример расчетной схемы: а – замена ветвей 1–6 рассчитываемого блока эквивалентными сопротивлениями Rэкв.1 и Rэкв.2; б – замена ветвей 9–12 рассчитываемого блока эквивалентным сопротивлением Rэкв.3 лучим ветвь с эквивалентным аэродинамическим сопротивлением Rэкв.3 (рис. 4, б). Аналогичным образом находится эквивалентное аэродинамическое сопротивление всей подземной части рудника (Rпчр). Аэродинамическое сопротивление каждой отдельной выработки согласно [2] определится по формуле (6) где Ki – коэффициент, определяющий степень шероховатости, способ проходки и назначение i-й выработки [2]; Li и Si – длина (м) и сечение (м2) i-й выработки. Расчет поверхностных утечек воздуха Барометрическое давление в околоствольном дворе воздухоподающего ствола согласно [11] определяется по формуле (7) где ai – коэффициенты регрессии [11]; ni – степени при параметре [11]; Pa – атмосферное давление воздуха, мм рт. ст.; tнар – температура наружного воздуха, °C; Hвозд, Dвозд, Qвозд – соответственно глубина (м), диаметр воздухоподающего ствола (м) и объем воздуха, проходящего по нему (м3/с). Средняя плотность воздуха в вентиляционном стволе определяется по формуле [10] (8) где tод.вент – температура воздуха в околоствольном дворе вентиляционного ствола (в течение года имеет практически постоянное значение и находится по формулам [12]), °C; Y – температурно-влажностный градиент, °C/м [10, 12]; ρср.возд – средняя плотность воздуха в воздухоподающем (воздухоподающих) стволе (стволах), кг/м3; Rвозд – аэродинамическое сопротивление воздухоподающего (воздухоподающих) ствола (стволов), Н∙с2/м8. Средняя плотность воздуха в воздухоподающем стволе согласно [10] определяется по формуле (9) где X – температурно-влажностный градиент, °C/м [10, 12]; Hвозд – глубина воздухоподающего ствола, м. Величина разряжения создаваемого ГВУ находится по формуле (3). Формулы (8) и (9) справедливы для случая, если устья стволов находятся на равных абсолютных отметках. В противном случае необходимо произвести добавочные расчеты, учитывающие параметры наружного воздуха по алгоритму [13]. Развиваемое вентилятором ГВУ давление согласно [2] может быть определено по формуле (10) где A, B, C – коэффициенты кривой, соответствующей определенным углам установки лопаток направляющих аппаратов (центробежный вентилятор), рабочего колеса (осевой вентилятор) вентилятора или скорости вращения рабочего колеса ГВУ [2]. Тогда, решив уравнение (2) относительно уравнения (10), получим , (11) где N – коэффициент пересчета из мм рт. ст. в Па, равный 133,332. Аэродинамическое сопротивление канала в зависимости от формы поворота из ствола в канал, способа присоединения вентилятора к каналу, расположения канала по отношению к стволу и его геометрических размеров определяется по формулам [14, 15]. Обозначив ; и , получим (12) Из формул (12) и (1) находим величину поверхностных утечек воздуха. Заключение Таким образом, предложенный вариант позволит рассчитать величину поверхностных утечек воздуха при изменении режима проветривания (производительности вентилятора), атмосферных параметров (температуры наружного воздуха и атмосферного давления), аэродинамических сопротивлений канала ГВУ и подземной части рудника, в результате чего появляется возможность применять его как для действующих, так и для проектируемых рудников. На действующих рудниках с учетом изменяющихся климатических параметров можно спрогнозировать величину поверхностных утечек и выбрать режим работы ГВУ. Для проектируемых рудников предложенный алгоритм будет иметь большее значение, так как от величины поверхностных утечек воздуха зависит выбор вентилятора ГВУ и режим его работы (характеристика и КПД), а следовательно, энергоэффективность проветривания.
About the authors
Aleksandr Viktorovich Nikolaev
Perm National Research Polytechnic University
Author for correspondence.
Email: nikolaev0811@mail.ru
614990, Perm, Komsomolsky ave., 29
assistant, Department for electrification and mining company automation, Perm National Research Polytechnic University
Nikolaj Ivanovich Alymenko
Mining Institute, Ural Branch of the Russian Academy of Sciences, Perm, Russian Federation
Email: nik.alymenko@yandex.ru
614007, Perm, ul. Sibirskaja, 78 а
doctor of technical sciences, professor, research engineer, Mining Institute, Ural Branch of the Russian Academy of Sciences
Ruslan Il'gizovich Sadykov
Perm National Research Polytechnic University
Email: sadykovrus@mail.ru
614990, Perm, Komsomolsky ave., 29
senior lecturer, Department for electrification and mining company automation, Perm National Research Polytechnic University
References
- Spravochnik porudnichnoj ventiljacii [Handbook of mine ventilation]. Ed. K.Z. Ushakova. Moscow: Nedra, 1977. 328 s.
- Mohirev N.N., Rad'ko V.V. Inzhenernye raschety ventiljacii shaht. Stroitel'stvo. Rekonstrukcija. Jekspluatacija [Engineering calculations of ventilation shafts. Construction. Reconstruction. Operation]. Moscow: Nedra-Biznescentr, 2007. 324 s.
- Burchakov A.S., Mustel' P.I., Ushakov K.Z. Rudnichnaja ajerologija [Mine aerology]. Moscow: Nedra, 1971. 376 s.
- Alymenko N.I., Minin V.V. Ventiljatornye ustanovki i ih primenenie [Ventilation system and its application]. Yekaterinburg: UrO RAN, 1999. 223 s.
- Gillies A.D.S., Hsin Wei Wu. A comparison of air leakage prediction techniques for auxiliary ventilation ducting systems. Proceedings Eighth US Mine Ventilation Symposium. Missouri, 1999, pp. 681–690.
- Andrew O.E. Ventilation of Lake Shore mines. Mining in Canada. Northern Mines Press LTD, 1957. 38 p.
- Morris I.N., Walker G. Changes in the approach to ventilation recent years. The Mining Eng, 1982, vol. 141, no. 244, pp. 401–413.
- Dzidziguri A.A., Cholikidze G.V. Calculation of the total air losses in long air ducts, soviet mining science. 1977, vol. 13, no. 5, pp. 333–341.
- Nikolaev A.V., Gavrilov V.M. O vozmozhnosti ispol'zovanija teplovoj depressii, voznikajuwej pri rabote nagrevatelej, raspolozhennyh v ust'e ventiljacionnogo stvola, dlja snizhenija poverhnostnyh utechek [The possibility of using the heat of depression that occurs when using a heater located at the mouth ventilation shaft, to reduce surface leakage]. Molodoj uchenyj. 2011, no. 6, pp. 85–89.
- Nikolaev A.V. Analiz teoreticheskoj formuly, opredeljajuwej velichinu estestvennoj tjagi, dejstvujuwej mezhdu vozduhopodajuwim i ventiljacionnym stvolami [Theoretical analysis of the formula that determines the amount of natural draft, acting between the air supply and ventilation shafts]. Geologija, geofizika i razrabotka neftjanyh i gazovyh mestorozhdenij, 2009, no. 10, pp. 72–75.
- Nikolaev A.V. Upravlenie teplovymi depressiiami v sistemakh ventiliatsii kaliĭnykh rudnikov: avtoref. dis. … kand. tekhn. nauk [Control of thermal depression in ventilation potash mines: Abstract of. thesis. ... candidate. technical. sciences.]. Perm, 2012. 20 s.
- Nikolaev A.V. Utochnenie formuly, opredeljajuwej velichinu estestvennoj tjagi, dejstvujuwej mezhdu vozduhopodajuwimi i ventiljacionnym stvolami [Refinement of the formula that determines the amount of natural draft, acting between the air supply and ventilation shafts]. Trudi III mezhdunar. konf. Ekaterinburg: Izd-vo Ural. gos. gorn. un-ta, 2010, pp. 246–250.
- Alymenko N.I., Nikolaev A.V. Influence of mutual alignment of mine shafts on thermal drop of ventilation pressure between the shafts. Journal of Mining Science, 2011, vol. 47, no. 5, pp. 636–642.
- Rukovodstvo po tipovym kanalam ventiljatorov dlja shaht Kuzbassa [Guidance on channel types of fans for Kuzbass mines]. Kemerovo: VostNII, 1964. 189 s.
- Alymenko N.I. Vybor konfiguracii ventiljacionnyh kanalov [Configurable for ventilation ducts]. Materialy ezhegod. nauch. sessii GI UrO RAN po rezul'tatam NIR v 2009 godu. Perm, 2010, pp. 227–230.
Statistics
Views
Abstract - 363
PDF (Russian) - 58
Refbacks
- There are currently no refbacks.