Бесплатный автореферат и диссертация по наукам о земле на тему

Газовыделение из разрабатываемых угольных пластов и вмещающих пород при интенсивной выемке метаноносных углей

ВАК РФ 25.00.20, Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика

Автореферат диссертации по теме "Газовыделение из разрабатываемых угольных пластов и вмещающих пород при интенсивной выемке метаноносных углей"

БУХТИЯРОВ Алексей Александрович

ГАЗОВЫДЕЛЕНИЕ ИЗ РАЗРАБАТЫВАЕМЫХ УГОЛЬНЫХ ПЛАСТОВ И ВМЕЩАЮЩИХ ПОРОД ПРИ ИНТЕНСИВНОЙ ВЫЕМКЕ МЕТАНОНОСНЫХ УГЛЕЙ

Специальность 25.00.20 Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика

Автореферат

диссертации на соискание ученой степени кандидата технических наук

1 9 ДПР 1ь\1

Тула 2012

005018116

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Тульский государственный университет» (ТулГУ) на кафедре геотехнологий и строительства подземных сооружений.

Научный руководитель: доктор технических наук, профессор КАЧУРИН Николай Михайлович.

Официальные оппоненты:

ЕФИМОВ Виктор Иванович, доктор технических наук, профессор, МГГУ, профессор кафедры «Экономика природопользования», г. Москва, ВАКУНИН Евгений Иванович, кандидат технических наук, доцент, ТулГУ, доцент кафедры «Аэрологии охраны труда и окружающей среды», г.Тула.

Ведущая организация: ОАО «Прокопгипроуголь», г. Прокопьевск Кемеровской области.

Защита диссертации состоится «¿§> ¿¿СсЛ 2012 г.вЛ час СОмка на заседании диссертационного совета Д 212.271.04 при Тульском государственном университете по адресу: 300012, г. Тула, просп. Ленина, 90, 6-й уч. корпус, ауд. 220.

С диссертацией можно ознакомиться в библиотеке Тульского государственного университета.

Отзывы на автореферат в двух экземплярах, подписанные и заверенные печатью организации, просим высылать по адресу: 300012, г. Тула, просп. Ленина, 92, Ученый совет ТулГУ, факс: (4872) 33-13-05, ' e-mail: galina stas@mail.ru.

Автореферат разослан « » ЄшаЄМ? ІОІІ. г.

Ученый секретарь диссертационного совета т I*? а.Б. Копылов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность. Практика эксплуатации угольных шахт России показывает, что уровень безопасности подземных горных работ по газовому фактору постоянно снижается. Это обусловлено переходом на технологии интенсивной отработки запасов, при которых нагрузка на очистные забои может превышать 10000 т/сут, и постоянным ростом глубины разработки. При этом устойчивое обеспечение потребностей страны во всех видах топлива и энергии будет по-прежнему осуществляться за счет увеличения их добычи. Наиболее перспективным в России является Кузнецкий угольный бассейн. Общие балансовые запасы угля в России составляют почти 200 млрд т, или 11,3 % мировых запасов. По производству угля Россия занимает 5-е место в мире после Китая, США, Индии и Австралии. Энергетической стратегией до 2020 года намечается увеличить добычу угля в России до 450 млн тонн, в том числе в Кузбассе - до 177 млн тонн. Кузбасс преодолел 150-миллионный рубеж добычи и выдает на-гора около 160 млн тонн высококачественного угля (подземная добыча - 46 %, открытая - 54 %).

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

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

весьма актуальным.

Диссертационная работа выполнялась в соответствии с тематическим планом НИР Научно-образовательного центра по проблемам рационального природопользования при комплексном освоении минерально-сырьевых ресурсов Аналитической ведомственной целевой программы «Развитие научного потенциала высшей школы (2009 - 2010 годы)» (per. номер 2.2.1.1/3942) и Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» (гос. контракт № 02.740.11.0319).

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

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

Основные научные положения состоят в следующем:

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

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

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

- газовый поток из надработанных пород в выработанное пространство убывает во времени по гиперболической зависимости и при фильтрационном критерии Фурье более 4 изменяется на 20...25 %;

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

Новизна разработанных научных положений заключается в следующем:

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

ся тем, что динамика газовыделений описывается с учетом технологических параметров;

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

диффузии в виде источников;

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

Обоснованность и достоверность теоретических положении, выводов и рекомендаций подтверждаются:

- корректной постановкой задач исследования, применением классических методов математической физики, математической статистики и теории вероятностей и современных достижений вычислительной математики; ,

- удовлетворительной сходимостью результатов прогноза с фактическими данными (отклонение не превышает 20 %) и большим объемом вычислительных экспериментов;

- значительным объемом шахтных наблюдений, а также результатами анализа плановых замеров (проанализированы данные по 26 очистным участкам).

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

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

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

Апробация работы. Научные положения и практические разработки диссертационной работы и отдельные ее разделы докладывались и обсуждались на научных семинарах кафедры геотехнологий и строительства подземных сооружений ТулГУ (г. Тула, 2009 - 2011 гг.), технических советах ОАО «ОУК «Южкузбассуголь» (г. Новокузнецк Кемеровской обл.,

2009 - 2011 гг.), 3-й и 4-й Международной конференции «Социально-экономические и экологические проблемы горной промышленности, строительства и энергетики» (г. Тула, 2008 - 2009 гг.), Международной конференции «Геомеханика. Механика подземных сооружений» (г. Тула,

2010 г.), Международной конференции «Ecology, Energy, Economy Safety in a Non Liner World. 3E - SAFETY» (г. Белград, 2009 г.); 3-й Международной конференции по проблемам рационального природопользования «Проблемы создания экологически рациональных и ресурсосберегающих технологий добычи и переработки отходов горного производства» (г. Тула, 2010 г.); 3-rd International Symposium ENERGY MINING «Energy Mining, New Technologies, Sustainable Development» (Сербия, г. Апатин, 2010 г.).

Публикации. По результатам выполненных исследований опубликовано 8 статей.

Объем работы. Диссертационная работа изложена на 180 страницах машинописного текста, состоит из 5 разделов, содержит 22 таблицы, 23 рисунков, список литературы из 158 наименований.

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Современные методы прогноза газовыделений и газовых ситуаций в горных выработках являются результатом исследований, выполненных ИПКОН РАН, ВостНИИ, ИГД им. A.A. Скочинского, МГГУ, Санкт-Петербургским горньм институтом, ТулГУ и др. Большой вклад в решение проблемы внесли академики РАН К.Н. Трубецкой, Ю.Н. Малышев, Д.А. Рубан, Д.Р. Каплунов, А.Е. Красноштейн; а также A.A. Айруни, К.К. Бусыгин, Ю.Ф. Васючков, Ф.С. Клебанов, А.Д. Кизряков, О.И. Касимов, С.П. Казаков, Н.М. Качурин, A.A. Мясников, Н.В. Ножкин, А.Э. Петросян, B.C. Пак, И.В. Сергеев, Э.М. Соколов, М.Б. Сулла, Н.И. Устинов, К.З. Ушаков, и др. Анализ основных научных и практических резуль-

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

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

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

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

3. Обосновать математические модели выделения газа с поверхности обнажения угольного пласта и из отбитого угля и разработать методические положения прогноза динамки выделения метана из разрабатываемого угольного пласта.

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

очистного участка.

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

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

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

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

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

дрг 8грг дгрг +;—— = у —-

д{ г Ы1 02

где р - давление свободного метана в трещиновато-пористой структуре угля; 2 - пространственная координата; 1 - время; ?г - период релаксации процесса ламинарной фильтрации метана в угольном пласте; Хул - пьезо-проводность угольного пласта, зависящая от фильтрационных и сорбци-онных свойств угля.

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

-+*г-£г = ху„-£Г, (1)

ß

p(z,0)~ p0=const, —p(z,0) = 0, dt

p(0,t) = pc = const, lim рфоо,

(2)

гдеpo,Pc~ давление свободного метана в угольном пласте и на поверхности обнажения пласта соответственно.

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

1уд„ = 1уд н ехр(~0,5tt~' )10 (0,5tt~1), где 1удн - начальная скорость газовыделения, м3/(м2-мин). Тогда дебит метана dln 0 в очистной забой с элементарной поверхности обнажения угольного пласта dS, учитывая полученную зависимость, можно определить как dl„ 0 - IydndS =

~ту.п vnfyd.HехР5t/tr)Ia(0,5t/tr)dt, где тум и Vn - мощность разрабатываемого угольного пласта и скорость подачи очистного комбайна. Представляя модифицированную функцию Бесселя в интегральном виде, метановыделение в очистной забой со всей поверхности обнажения угольного пласта определим по формуле

1по(твц) = тс1по(myntrVnIydH f1. Безразмерное значение метановыделе-

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

ЬочНКУп) и . 1

Т^о= | ехр{-г) <| \{ехр(т Со$в)+ехр(-т Соз0)\с10\с1г , (3) о 1о -I

где 10Ч - длина очистного забоя.

Вычислительный эксперимент проводился с использованием соотношения (3). Исследования проведены на четырех интервалах значений безразмерной длительности выемочного цикла. Результаты исследования закономерности = /(твч) показывают, что при г,.„ <г [О, 3] функция /(тв ) может быть аппроксимирована степенной зависимостью, для интервалов значений е [3, 10], гв.„ е [10, 50] и г,, в [50, 100] целесообразно использовать линейную аппроксимацию с различными начальными значениями и угловыми коэффициентами.

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

3,73 т°в'ц при 0 <тЙЧ < 3,

4,571 + 1,143 г,.„ при 3 <та ц < 10, (4)

10 + 0,5 тв ц при 10 <гв !/ < 50,

25 + 0,26 т„ч при 50 <тв ц < 100.

Из формул (3) - (4) следует, что интенсивность метановыделения непосредственно связана с планограммой работ в очистном забое. Анализ этой зависимости свидетельствует о том, что при выемочном цикле мета-новыделение нарастает за счет увеличения площади газоотдающей поверхности. Газоотдающая поверхность находится под перепадом давления, приблизительно равным значению р20-р2с- А площадь частично дегазированной поверхности обнажения уменьшается с той же скоростью, с которой увеличивается площадь недегазированной газоотдающей поверхности, поэтому при работе выемочного комбайна формируется стационарное состояние, определяемое формулой (3).

Динамика газоносности отбитого угля в общем виде описывается следующим дифференциальным уравнением:

V )&!!_ =-3». ](2ЯР~'х д( Кк П>дг} раГу21

4о =

•хехр\-

(2Я)2

(1(211), (5)

где а>оу- коэффициент массообмена; х0 - остаточная газоносность угля при давлении р0; Ям - показатель степени измельчения угля, являющийся параметром масштаба; mi - показатель способности угля к измельчению, являющийся параметром формы; 2Rmtn - минимальный размер угольных частиц, образующихся при резании углей.

Решение уравнения (5) для условий х0(т},0) = xo(0,t) = х3 = const позволяет рассчитать метановыделение из отбитого угля при высокопроизводительных современных технологиях выемки угля. Окончательно расчетная формула получена в следующем виде: 10 У („ ) = (pjy сооуКпуАк х, x(l,26-0,78-tt4), где Кпу- коэффициент погрузки угля на конвейер; Ак - производительность комбайна. В целом вычислительные эксперименты показали, что метановыделение из отбитого угля - весьма динамичный процесс, имеющий тесную связь с планограммой работ в лаве, при этом установлено, что челноковая схема работы комбайна снижает неравномерность метановыделения из отбитого угля на 20...30 % и его длительность - на 60...70 % по сравнению с односторонней схемой работы.

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

др д3Р д2Р

(б)

где <p = km(/jm0ß)~l; Tj-kg'kml2-, р- давление газа в трещинах; л:,-пространственные координаты; ^ , к6,- газовая проницаемость трещин и породных блоков соответственно; ß- коэффициент сжимаемости метана; I - среднее значение характерного размера породных блоков.

Решение уравнения (6) для условий р (г, 0) = ра;

р(+0, t) = Ра-(Ро~Рс)' exp(-7?~V); limp(z,t) = const позволило получить закономерность для объема метана, выделившегося в единицу времени с единичной площади подработанных пород (1П п):

¡ил = (мрЛ\Ш°'5 (А ~Рс){ра ~(Ро ~ А)[1-ехр(-0]}/(0, (7)

с

где £ = Цх(р1 ;/(£) = ехр(Ч) |ехр(0,5г)/0 (0,5г)Л.

о

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

1нл=0,282р:\р1-р])^, (8)

где кн, т„- проницаемость и пористость надработанных пород соответственно.

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

~-зЬУюр.пХп.п/2(С2) &„„

м3

Уочт.уу V К т

\ Утр.нХн.п

(10)

К,'п.Гу V т т

гДе Чг,„р „ 1 Ятр,, ' относительное метановыделение из подработанных и надработанных пород соответственно, м3/т; 12 - шаг обрушения пород основной кровли, м; упйрп, утр и - плотность подработанных и надработанных пород, т/м3; хнп, хпп - природная газоносность подработанных и надработанных пород, кпп, кнп - газовая проницаемость подработанных и надработанных пород, м2; С'2 = 4,87 ■ 1018 кбт~';

//0 = ¡ехр(-0,5т)10(-0,5т)с1т-ехр(-£)■

о

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

к„ = 0,266 -10~21г (ЯГТ)0,5 тп3р. (И)

Если перенос происходит в пористой сорбирующей среде, то для давления газа до 5 МПа можно использовать уравнение Лэнгмюра, тогда формула для расчета газовой проницаемости принимает вид

3

к„ =1,88-Ю^{КгТ)°^г(1+Ьлр) г>1рА0,1-апЪа +0,707(1+Ьлр)2

'.(12)

Формула (12) позволяет рассчитать газовую проницаемость сорбирующей пористой среды, сложенной материалом, изотерма сорбции которого выпукла (по классификации С. Брунауэра - это изотерма первого рода).

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

Изгиб породных слоев приводит к образованию областей сжатия в их нижней части и областей растяжения в верхней части, а так как [бсж}>[бр], где [бсж] и [бр] - предельные значения нормальных напряжений сжатия и растяжения, то разрушение слоя происходит за счет растягивающих усилий. Максимальное значение изгибающего момента будет в точке жесткого защемления слоя, там же будут и наибольшие внутренние напряжения.

Тогда предельные значения толщины слоев породы непосредственной и основной кровли (II, и Нз) определяются по следующим формулам:

( п ^

Н] = ([бр]Са апГ°'51,732/, £/,//,.

Н2 = ([бр JC.es ап 0,8661

Л0,5

5>д

где я - плотность / - го породного слоя мощностью Н,, кг/м3; ап - угол падения разрабатываемого угольного пласта, град.

Величины, определяемые по этим формулам, в сумме представляют собой глубину проникновения разрывных нарушений подработанных пород. На этом расстоянии наблюдается существенное увеличение трещи-новатости и проницаемости пород. Дополнительная пористость, образовавшаяся в обрушенной непосредственной и основной кровле, (т1 и т^

определяется по формулам т1 = тв п{твп+ Н ¡у1,

Щ = тв.п (твп +н, + Н2)~'.

Породы почвы угольного пласта до начала выемки нагружены верти-

кальной составляющей распределенной нагрузки вышележащих пород и находятся в объемном напряженном состоянии. На поверхности контакта угольного пласта с почвой а = усрН, где уср - средневзвешенная плотность вышележащих пород. В процессе выемки пласта происходит обнажение этой поверхности, и породный слой будет находиться под действием распределенной нагрузки, направленной вертикально вверх, которая равна 0,7 уср Н сох'1 ап. Рассматривая породы почвы в виде балки на упругом основании, можно использовать дифференциальное уравнение упругой линии породного слоя на упругом основании с?г /с!х + Кум г (Е*1ин)~ = д3(Е/ия)~', где Е - модуль упругости пород; - момент инерции сечения слоя; Ку.о - коэффициент упругого отпора. Геометрические параметры рассматриваемого слоя длиной 10Ч оказывают существенное влияние на проявление краевого эффекта. Так, при расчетах на прочность слой на упругом основании можно рассматривать как бесконечно длинный, если Д,ч > 3, где (¡-[0,25 Ку.0 (Е^,/1]0'25. Расчеты показали, что Щт в десятки раз больше трех. Таким образом, породный слой можно рассматривать как полубесконечную балку, так как при Д,, »3 прогибы и моменты возле одного края не зависят от условий закрепления другого края.

Поэтому решение уравнения упругой линии породного слоя имеет

следующий вид: г(х) = д3К~'0[1 -ея$(-рх)1Р<>5рх+• Тогда,

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

пластов Н3 = 0,213 (Н/[бр])0'667 . Исходной информацией для

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

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

1-ехр--—*

Ьп а .

\ ЬП.В

I +0,5ехр

где К = -4=; Ь = —+——; иср - средняя скорость движения воздуха по 4D 4D Lns

подготовительной выработке; LnB - проектная длина подготовительной выработки; 1П£, Q п.в - абсолютная газообильность и объем подготовительной выработки; С = с - сн; с - объемная концентрация рассматриваемой газовой примеси в воздухе выработки (для условий Кузбасса - это метан или углекислый газ); сн - объемная концентрация газовой примеси на свежей струе, поступающей в подготовительную выработку.

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

C(xJ) = I0„Q0'y ([]-ехр(- (u)(ELoy exp {-x(SL0V ) х

'x[l-exp\-(и) (ZLoy (t-x {иУ )]} <r0(t-x (и)-' )J>, (14)

где Qgy - количество воздуха, протекающего по выработкам.

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

ЗАКЛЮЧЕНИЕ

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

ние аварийности на угольных шахтах при отработке глубоких горизонтов с использованием технологии «шахта - лава».

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

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

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

3. Установлено, что газовая проницаемость горного массива является обобщенной характеристикой свойств газа, коллекторских свойств горного массива и стохастичности фильтрационного течения, а оценки ее математического ожидания изменяются от 10'17 до Ю"20 м2 дня углей и пород Кузбасса и Донбасса.

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

5. Метановыделение из подработанных вымещающих пород в выработанное пространство возрастает, если произведение времени на пьезопро-водность транспортных трещин, отнесенное к показателю газообмена с породными блоками, меньше 1, а затем наблюдается снижение дебита и при значениях этого комплекса, равных 4...10, газовый поток остается практически постоянным.

6. Газовый поток из надработанных пород в выработанное пространство убывает во времени по гиперболической зависимости и при фильтрационном критерии Фурье более 4 изменяется на 20...25 %.

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

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

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

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

1. Качурин Н.М., Борщевич А.М., Бухтияров A.A. Теоретическое обоснование ограничения нагрузки на очистной забой по газовому фактору// Изв ТулГУ. Науки о Земле. 2009. Вып. 4. С. 84-88.

2. Качурин Н.М., Борщевич A.M., Бухтияров A.A. Теоретическое обоснование прогноза выделения метана с поверхности обнажения разрабатываемого угольного пласта при интенсивной выемке угля// Изв ТулГУ. Науки о Земле. 2009. Вып. 4. С. 93-97.

3. Качурин Н.М., Борщевич A.M., Бухтияров A.A. Прогноз выделения метана с поверхности обнажения разрабатываемого угольного пласта и нагрузка на лаву при интенсивной выемке угля// Безопасность жизнедеятельности. 2010. Вып. 5. С. 19-24.

4. Безопасность геотехнологий добычи угля по газовому фактору. М.: Безопасность жизнедеятельности/ Н.М. Качурин, A.M. Борщевич, О.Н. Качурина, A.A. Бухтияров. Вып. 5.2010. С. 24-28.

5. Качурин Н.М., Борщевич А.М., Бухтияров A.A. Предельно допустимая нагрузка на очистной забой по газовому фактору//3-я Международная конференция по проблемам рационального природопользования. «Проблемы создания экологически рациональных и ресурсосберегающих технологий добычи полезных ископаемых и переработки отходов горного производства». Тула. 2010. С. 188 - 193.

6. Качурин Н.М., Борщевич A.M., Бухтияров A.A. Газовыделения с поверхности обнажения разрабатываемого угольного пласта на метанообильных шахтах// 3-я Международная конференция по проблемам рационального природопользования. «Проблемы создания экологически рациональных и ресурсосберегающих технологий добычи полезных ископаемых и переработки отходов горного производства». Тула. 2010. С. 193 - 200.

7. Метановыделение из отбитого угля при интенсивной отработке угольного пласта/ Н.М. Качурин, А.М. Борщевич, A.A. Бухтияров, И.И. Сарычева // 3-я Международная конференция по проблемам рационального природопользования. «Проблемы создания экологически рациональных и ресурсосберегающих технологий добычи полезных ископаемых и переработки отходов горного производства». Тула. 2010. С. 200 - 205.

8. Качурин Н.М., Борщевич A.M., Бухтияров A.A. Метановыделение в очистной забой из подработанных и надработашшх пород». Изв. ТулГУ. Науки о Земле. 2011. Вып. 1.С. 62-68.

Изд. лиц. ЛР № 020300 от 12.02.97. Подписано в печать 13.02.2012 г. Формат бумаги 60x84 '/16. Бумага офсетная. Усл.-печ. л. 1,5. Уч.-изд. л. 1,3. Тираж 100 экз. Заказ 3 ?

Тульский государственный университет. 300012, г. Тула, просп. Ленина, 92.

Отпечатано в издательстве ТулГУ. 300012, г. Тула, лросп, Ленина, 95.

Текст научной работыДиссертация по наукам о земле, кандидата технических наук, Бухтияров, Алексей Александрович, Тула

61 12-5/2397

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

Тульский государственный университет

На правах рукописи

БУХТИЯРОВ Алексей Александрович

ГАЗОВЫДЕЛЕНИЕ ИЗ РАЗРАБАТЫВАЕМЫХ УГОЛЬНЫХ ПЛАСТОВ И ВМЕЩАЮЩИХ ПОРОД ПРИ ИНТЕНСИВНОЙ ВЫЕМКЕ МЕТАНОНОСНЫХ УГЛЕЙ

Специальность 25.00.20 - Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика

диссертации на соискание ученой степени кандидата технических наук

Научный руководитель - д.т.н., профессор Н.М. Качурин

Тула 2012

СОДЕРЖАНИЕ

ВВЕДЕНИЕ............................................................................... 4

1. АНАЛИТИЧЕСКИЙ ОБЗОР И ПОСТАНОВКА ЗАДАЧ ИССЛЕДОВАНИЙ.................................................................. Ю

1.1. Системные свойства минерально-сырьевых ресурсов и проблемы безопасности при подземном способе добычи полезных ископаемых ...... 10

1.2. Перспективы развития подземной угледобычи при отработке газоносных углей Кузбасса............................................................... 15

1.3. Оценка риска и системный подход к профилактике аварий, вызванных аэрогазодинамическими условиями.......................................... 21

1.4. Теория и практика моделирования аэрогазодинамических процессов

в угольных шахтах..................................................................... 29

Выводы................................................................................... 48

Цель и идея работы. Постановка задач исследований......................... 49

2. ТЕХНОЛОГИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПОДЗЕМНЫХ ГОРНЫХ РАБОТ И ВЫДЕЛЕНИЕ МЕТАНА ИЗ УГОЛЬНЫХ ПЛАСТОВ И ВЫРАБОТАННЫХ ПРОСТРАНСТВ

ОЧИСТНЫХ УЧАСТКОВ ....................................................... 51

2.1. Состояние угледобычи и горно-геологический прогноз условий залегания угольных пластов при комплексном освоении месторождений 51

2.2. Геотехнологические аспекты и перспективы комплексного освоения угольных месторождений............................................................ 56

2.3. Технологии подготовки и отработки метаноносных угольных пластов ....................................................................................... 60

2.4. Структура газового баланса угольных шахт................ ................. 76

2.5. Динамика выделения метана из различных источников в горные выработки ...................................................................................

2.6. Аэрогазодинамическая характеристика шахт ОАО ОУК «Южкуз-

бассуголь»............................................................................... 88

Выводы................................................................................... 92

3. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ АЭРОГАЗОДИНАМИКИ ОЧИСТНОГО УЧАСТКА........................ 95

3.1. Метановыделение из отбитого угля при интенсивной выемке угля ... 95

3.2. Выделение метана с поверхности обнажения разрабатываемого угольного пласта....................................................................... 104

3.3. Метановыделение в очистной забой из вмещающих пород..........................113

3.4. Математические модели газовых ситуаций на очистных и подготовительных участках........................................................................................................................................120

Выводы......................................................................................................................................................................126

4. ГАЗОДИНАМИЧЕСКИЕ СВОЙСТВА РАЗРАБАТЫВАЕМЫХ УГОЛЬНЫХ ПЛАСТОВ И ВМЕЩАЮЩИХ ПОРОД............................................129

4.1. Давление газа в разрабатываемых угольных пластах..........................................129

4.2. Газоносность вмещающих пород и смежных пластов........................................135

4.3. Газовая проницаемость разрабатываемого угольного пласта......................141

4.4. Пористость и газовая проницаемость вмещающих пород и смежных

угольных пластов............................................................................................................................................146

Выводы......................................................................................................................................................................156

5. ПРАКТИЧЕСКАЯ АПРОБАЦИЯ ПРОГНОЗА ГАЗОВЫДЕЛЕНИЯ ИЗ РАЗРАБАТЫВАЕМЫХ УГОЛЬНЫХ ПЛАСТОВ И ВМЕЩАЮЩИХ ПОРОД ПРИ ИНТЕНСИВНОЙ ВЫЕМКЕ МЕТАНОНОСНЫХ УГЛЕЙ................................................................................................................................................................158

5.1. Методические положения управления метановыделением и газовыми ситуациями....................................................................................................................................................158

5.2. Расчет количества воздуха для выемочного участка 44-10 пласта «Кыргайский - 44», ОАО «Шахта Кыргайская», ООО «Ерунаковская угольная компания»......................................................................................................................................163

5.3. Прогноз метановыделения из вмещающих пород на очистном участке коренной западной лавы шахты «Суходольская-Восточная»..........................166

Выводы......................................................................................................................................................................167

ЗАКЛЮЧЕНИЕ................................................................................................................................................169

СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ............................................................171

ВВЕДЕНИЕ

Актуальность. Практика эксплуатации угольных шахт России показывает, что уровень безопасности подземных горных работ по газовому фактору постоянно снижается. Это обусловлено переходом на технологии интенсивной отработки запасов, при которых нагрузка на очистные забои может превышать 10000 т/сут, и постоянным ростом глубины разработки. При этом устойчивое обеспечение потребностей страны во всех видах топлива и энергии будет по-прежнему осуществляться за счет увеличения их добычи. Наиболее перспективным в России является Кузнецкий угольный бассейн. Общие балансовые запасы угля в России составляют почти 200 млрд. т или 11,3% мировых запасов. По производству угля Россия занимает 5 место в мире после Китая, США, Индии и Австралии. Энергетической стратегией до 2020 года намечается увеличить добычу угля в России до 450 млн. тонн, в том числе в Кузбассе - до 177 млн. тонн. Кузбасс преодолел 150-миллионный рубеж добычи и выдает на-гора около 160 млн. тонн высококачественного угля (подземная добыча - 46 %, открытая - 54 %).

Высокая метаноносность угольных пластов Кузбасса (до 20 - 25 м /т) является причиной возникновения газового барьера для современных высоко производительных технологий добычи угля. Динамика взрывов метановоз-душной смеси (МВС), количества пострадавших и погибших последних лет показывает, что увеличение добычи угля может осуществляться только при условии снижения риска такого вида аварий.

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

Важнейшей подсистемой технологии снижения риска и локализации последствий взрывов метановоздушной смеси является компьютерная технология оценка динамики концентрации метана в воздухе. Оценка надежности технологии аэрологической защиты показывает, что по своим характеристикам система защиты от взрывов MB С в угольных шахтах России не изменилась и остается на уровне угольной промышленности СССР. Так оценка интенсивности отказов этой системы в СССР составляла 9,43 1/год, а в РФ составляет 9,38 1/год. А эффективность системы защиты от поражающих факторов снизилась по сравнению с угольной промышленностью СССР в 5 раз.

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

Диссертационная работа выполнялась в соответствии с тематическим планом НИР Научно-образовательного центра по проблемам рационального природопользования при комплексном освоении минерально-сырьевых ресурсов Аналитической ведомственной целевой программы «Развитие научного потенциала высшей школы (2009 - 2010 годы)» (per. номер 2.2.1.1/3942) и Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» (гос. контракт № 02.740.11.0319).

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

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

Основные научные положения состоят в следующем:

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

■ динамика остаточной метаноносности отбитого угля описывается гиперболическим уравнением со стоком, учитывающим гранулометрический состав угля, скорость подачи комбайна и транспортирования угля, его фильтрационные и диффузионные свойства;

■ метановыделение из подработанных вымещающих пород в выработанное пространство возрастает, если произведение времени на пьезопровод-ность транспортных машин, отнесенное к показателю газообмена с породными блоками, меньше 1, а затем наблюдается снижение дебита и при значениях этого комплекса, равных 4-10 газовый поток остается практически постоянным;

■ газовый поток из надработанных пород в выработанное пространство убывает во времени по гиперболической зависимости м при фильтрационном критерии Фурье более 4 изменяется на 20-25%;

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

Новизна разработанных научных положений заключается в следующем:

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

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

■ предложены математические модели для прогноза газовых ситуаций в горных выработках, отличающиеся тем, что расчет полей концентрации газов осуществляется на основе установленных закономерностей динамики газовыделений, входящих в уравнение конвективно-турбулентной диффузии в виде источников;

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

Обоснованность и достоверность теоретических положений, выводов и рекомендаций подтверждается:

■ корректной постановкой задач исследования, применением классических методов математической физики, математической статистики и теории вероятностей и современных достижений вычислительной математики;

■ удовлетворительной сходимостью результатов прогноза с фактическими данными (отклонение не превышает 20%) и большим объемом вычислительных экспериментов;

■ значительным объемом шахтных наблюдений, а также результатами анализа плановых замеров (проанализированы данные по 26 очистным участкам).

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

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

Практическая реализация выводов и рекомендаций. Основные выводы и рекомендации работы по обоснованию возможности увеличения добычи углей в Кузнецком и Донецком угольных бассейнах, направленные на обеспечение безопасных условий развития добычи угля с учетом расширения масштабов производства действующих предприятий, нового строительства и освоения остаточных запасов углей ликвидируемых шахт, использованы при корректировке программы развития угольных бассейнов на перспективу и реализации энергетической стратегии России на период до 2020 года . Теоретические результаты и технические решения включены в учебные курсы по аэрологии горных предприятий, а так же использованы при выполнении хоздоговорных и госбюджетных НИР в Тульском государственном университете.

Апробация работы. Научные положения и практические разработки диссертационной работы, и отдельные ее разделы докладывались и обсуждались на научных семинарах кафедры геотехнологий и строительства подземных сооружений ТулГУ (г. Тула, 2009 - 2011 гг.), технических советах ОАО ОУК «Южкузбассуголь» (г. Новокузнецк Кемеровской обл., 2009 - 2011 гг.), 3-й и 4-й Международной конференции «Социально-экономические и экологические проблемы горной промышленности, строительства и энергетики» (г. Тула, 2008 - 2009 гг.), Международной конференции «Геомеханика. Механика подземных сооружений» (г. Тула, 2010 г.), Международной конференции «Ecology, Energy, Economy Safety in a Non Liner World. 3E - SAFETY» (г. Белград, 2009 г.); 3-я Международная конференция по проблемам рационального природопользования «Проблемы создания экологически рациональных и ресурсосберегающих технологий добычи и переработки отходов горного производства» (г. Тула, 2010 г.); 3-rd International Symposium ENERGY MINING «Energy Mining, New Technologies, Sustainable Development» (Сербия, г. Апатин, 2010 г.).

Публикации. По результатам выполненных исследований опубликовано 8 статей.

Объем работы. Диссертационная работа изложена на 180 страницах машинописного текста, состоит из 5 разделов, содержит 22 таблиц, 23 рисунков, список литературы из 158 наименований.

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

1. АНАЛИТИЧЕСКИЙ ОБЗОР И ПОСТАНОВКА ЗАДАЧ ИССЛЕДОВАНИЙ

1.1. Системные свойства минерально-сырьевых ресурсов и проблемы безопасности при подземном способе добычи полезных ископаемых

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

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

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

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

Здесь алгоритмичность следует понимать �