CC BY
192
5. Гулий В.Н., Иевлев А.А. Состав и особенности фосфатов месторождения Бирикээн // Литология и полезные ископаемые. -1990. — № 6. -С. 44-55.
6. Боярко Г.Ю., Локтионов А.В. Отражение фосфатоносных кор выветривания бирикээнского типа в геофизических полях // Геофизические методы при разведке недр и экологических исследованиях. — Томск: ТПУ, 2003. — С. 195-198.
7. Боярко ГЮ. Геохимические поиски фосфатоносных кор выветривания бирикээнского типа (Южная Якутия) // Проблемы поисковой и экологической геохимии Сибири. — Томск: ТПУ, 2003. —С. 139-141.
8. Боярко Г.Ю. Применение ядерно-физических методов для анализа фосфатных руд месторождений Южной Якутии // Известия Томского политехнического университета. — 2003. — Т 306. — № 6. —С. 18-20.
9. Филько А.С. Файзуллин Р.М., Карпова М.И. и др. Фосфатные руды России // Минеральные ресурсы России. Экономика и управление. —1994. — № 5. — С. 18-25.
10. Киперман Ю.А., Комаров М.А., Филько А.С. Особенности минерально-сырьевой базы фосфатов // Минеральные ресурсы России. Экономика и управление. —1996. — № 5. — С. 13—16.
11. Аксенов Е.М., Ведерников Н.Н., Чуприна Н.С., Рябикин В.В. Агрохимическое и горно-рудное сырье на рубеже XXI в. // Минеральные ресурсы России. Экономика и управление. — 2000.
— № 5-6. —С. 7-15.
12. Фархутдинов Р.З., Карпова М.И., Зеленихин В.А. и др. Минерально-сырьевая база фосфатных и нетрадиционных агрохимических руд Сибири и Дальнего Востока: перспективы освоения // Минеральные ресурсы России. Экономика и управление. —1999. —№ 6. —С. 6-11.
13. Боярко Г.Ю. Перспективы освоения Бирикээского месторождения // Проблемы освоения и перспективы развития ЮжноЯкутского региона. — Якутск: Якутский госуниверситет, 2001.
— С. 24-27.
УДК 622.831.232
ИССЛЕДОВАНИЕ ВЛИЯНИЯ ДВИЖУЩЕГОСЯ ОЧИСТНОГО ЗАБОЯ НА ХАРАКТЕР ЗАВИСАНИЯ И ЦИКЛИЧЕСКОГО ОБРУШЕНИЯ ПОДРАБОТАННЫХ ПОРОД КРОВЛИ УГОЛЬНЫХ ПЛАСТОВ
Л.Д. Павлова, В.Н. Фрянов
ГОУ ВПО "Сибирский государственный индустриальный университет". г. Новокузнецк E-mail: rector@sibsiu.ru
Используя пространственную геомеханическую модель деформирования и разрушения горных пород в зоне влияния горной выработки, разработанную на основе метода конечных элементов, проведено исследование влияния движущегося очистного забоя на характер изменения коэффициента концентрации вертикальных напряжений впереди очистного забоя, зависания и циклического обрушения подработанных пород кровли угольных пластов.
Введение
В последнее десятилетие резко возросла интенсивность подземной разработки угольных месторождений. За этот период нагрузка на очистной комплексно-механизированный забой возросла в 3...4 раза, а длина лавы увеличилась от 120 до 300 м. Закономерности геомеханических процессов для таких забоев изучены недостаточно, что подтверждается большой аварийностью забоев и высоким уровнем группового травматизма (шахты "Зыряно-вская", "Тайжина" в Кузбассе и др.).
Изучением геомеханических процессов на угольных шахтах в окрестности очистных забоев занимались многие ученые: Ф.А. Белаенко, А.А. Борисов, Г.Н. Кузнецов, В.Д. Слесарев, П.М. Цимба-ревич, И.Л.Черняк и др. Общие закономерности сдвижения подрабатываемого массива горных пород неоднократно описаны в литературе [1-5 и др.].
Механизм разрушений и обрушений толщи пород кровли разработан Г.Н. Кузнецовым [1]. Он доказал, что обрушение слоев непосредственной кровли в выработанном пространстве будет иметь место при условии, если толщина слоев кровли меньше мощности пласта. Характер первого обру-
шения основной кровли в срединном сечении выработанного пространства, ее взаимодействие с непосредственной кровлей и характер опорного давления на краевые зоны разрабатываемого пласта схематично показаны на рис. 1.
Как отмечается в [2] породы над зоной интенсивного разрыхления деформируются и разрушаются подобно тонким плитам, в режиме первого обрушения — подобно плите, защемленной на прямоугольном опорном контуре, в режиме установившегося обрушения — подобно консольной плите, защемленной по трем сторонам опорного контура.
Общий характер распределения опорного давления в плоскости пологого пласта, рассчитанный по результатам физического моделирования на объемных моделях, проверенный и многократно подтвержденный шахтными и лабораторными исследованиями описан А. А. Борисовым [3, 4]. По характеру изменений опорного давления он выделяет следующие зоны: динамических проявлений опорного давления, затухания динамических проявлений опорного давления, статического состояния опорного давления.
Анализ закономерностей проявления горного давления аналитическими методами выполнен Г.Н. Кузнецовым, А.М. Линьковым, И.М. Петуховым, К.В. Руппенейтом, ГН. Савиным и другими авторами. В аналитических методах исследования вопросов механики горных работ для расчета напряженного состояния массива в зонах влияния подземных выработок используются положения теории упругости, пластичности, ползучести. Учет влияния неоднородности состава пород междуп-ластий на напряженное состояние толщи и угольных пластов выполнялся по методу эквивалентного слоя, широко применяемого в механике грунтов.
Разработанные в настоящее время аналитические методы расчета напряжений в зонах влияния подземных выработок имеют ограниченную область применения и используются, главным образом, на стадии прогнозной оценки горно-геологических и горнотехнических условий отработки, для выбора способов и средств управления горным давлением.
Авторами [5] были проведены исследования характерных зон и параметров сдвижения массива под действием очистных работ на моделях из эквивалентных материалов. Измерения тензодатчика-ми вертикальных напряжений в подработанном массиве показали, что в 30...40 м впереди забоя напряжения возрастают в 1,5...2,0 раза, затем над забоем, или с некоторым опережением забоя, происходит снижение уровня напряжений практически до нуля. В дальнейшем происходит восстановление напряжений. Полностью до исходного уровня напряжения восстанавливаются на расстоянии 50.60 м за забоем.
Следует отметить, что одним из недостатков существующих методов прогнозирования напряженно-деформированного состояния углепородного массива при проведении горных выработок является решение двумерной задачи, что не обеспечивает корректного результата на сопряжении очистного забоя и угольного массива.
Кроме того, сложность расчетов напряженного состояния толщи пород при разработке свит пластов обусловлена необходимостью учета взаимного
влияния выработок на различных пластах и неоднородностью состава углепородного массива. В реальных условиях учитывать влияние всех слоев горных пород на напряженное состояние толщи весьма трудоемко и сложно из-за недостаточной надежности определения механических характеристик на контактах слоев.
В этой связи научные исследования геомехани-ческих процессов в окрестности высокопроизводительного комплексно-механизированного забоя путем решения пространственной задачи геомеханики с использованием численных методов являются актуальными.
Объект и методы исследований
Для сравнения фактических данных с численными результатами количественного прогнозирования геомеханических процессов в углепородном массиве при движении очистного забоя в качестве базового варианта объекта исследований выбран выемочный участок 26-26 шахты "Абашевская" в Кузбассе.
При расчете параметров напряженно-деформированного состояния углепородного массива в зоне влияния горной выработки использовалась пространственная геомеханическая модель деформирования и разрушения горных пород, разработанная авторами на основе метода конечных элементов [6]. Прогнозирование геомеханических процессов в углепородном массиве выполнялось на основе анализа пространственного напряженно-деформированного состояния и использования обобщенного критерия прочности Кулона-Мизеса. При моделировании механизма обрушения горных пород осуществлялось последовательное накопление повреждений с учетом функции времени. Для имитации движения очистного забоя выполнялся циклический сдвиг системы координат по мере продвижения забоя.
Учитывая сложность поставленной задачи и ограничения, связанные с техническими возможностями вычислительной техники, проводились исследования характера изменения коэффициента концентрации вертикальных напряжений впереди
У. и Т ■
г ?
г ■ Зет ИЛ} цио ннь й ш тре] с
Г а Л V
г Л -
Г & X КМ 3 ■—.
г я Л Н
Г О 2
г
Г Со
н зеие рнь и ш тре: <
г
г н—1—1 X, \
0 20 40 60 80 100120140
I выработанное пространство
Рис. 2. Схема движения очистного забоя, ограниченного по восстанию и падению угольным массивом на: а) 10; б) 40 м
очистного забоя, зависания и циклического обрушения подработанных пород кровли при движении длинного очистного забоя, ограниченного по восстанию и падению угольным массивом. Имитация движения очистного забоя от монтажной камеры выполнялась с шагом 10 м (рис. 2, а, б).
На каждом шаге моделирования вычислялись горизонтальные и вертикальные смещения, деформации, напряжения, определялись коэффициент концентрации вертикальных напряжений и остаточная прочность угля и пород в зоне влияния выработанного пространства. Наиболее представительными оказались графики изменения коэффициента концентрации вертикальных напряжений в угольном массиве впереди очистного забоя, а также формы и размеры зон зависания и обрушения подработанных пород кровли.
Исследование влияния движущегося забоя на характер зависания и циклического обрушения подработанных пород кровли угольных пластов
В результате проведенных исследований установлено, что при движении очистного забоя зона обрушения пород кровли имеет ассиметричную форму. По мере движения забоя от монтажной камеры происходит последовательное обрушение пород непосредственной кровли в середине выработанного пространства. Высота зоны обрушения увеличивается пропорционально расстоянию от монтажной камеры до забоя. Наиболее интенсивно разрушение породы происходит вблизи движущегося очистного забоя (рис. 3).
Рис. 3. Схема формирования зоны обрушенных пород подработанной кровли
Над угловыми участками выработанного пространства происходит зависание пород непосредственной кровли, что приводит к увеличению напряжений в угольном массиве.
При наличии в подрабатываемой толще угольных пластов разрушение угля в них опережает разрушение пород непосредственной кровли. Над выработанным пространством формируется порода-мост, мощностью равной расстоянию от кровли отрабатываемого пласта до подрабатываемого пласта, т.е. происходит расслоение подрабатываемой толщи и формирование двух породных плит (рис. 4).
Наличие таких пород-мостов над выработанным пространством приводит к росту напряжений в отрабатываемом пласте, увеличению длины консоли при движении забоя и динамическому обрушению породы-моста при достижении предельного значения длины консоли.
Рис. 4. Схема формирования зависающей консоли подработанных пород кровли
Динамическое обрушение пород кровли сопровождается заполнением выработанного пространства, вытеснением метано-воздушной смеси и угольной пыли в горные выработки и, при наличии источника открытого огня, подземным взрывом.
Исследование влияния движущегося забоя на характер изменения коэффициента концентрации вертикальных напряжений впереди очистного забоя
По мере движения от монтажной камеры происходит почти линейное увеличение коэффициента концентрации вертикальных напряжений (рис. 5). Прирост вертикальных напряжений происходит за счет увеличения веса зависающих пород кровли.
При движении очистного забоя одновременно с зависанием пород основной кровли и налегающей толщи происходит обрушение пород непосредственной кровли, что приводит к частичному снижению веса зависающих подработанных пород. Это обрушение пород непосредственной кровли зафиксировано на графиках на рис. 5 в виде снижения коэффициента концентрации вертикальных напряжений после каждого шага подвигания забоя.
Как следует из графиков на рис. 5, снижение коэффициента концентрации вертикальных напряжений в среднем сечении лавы (точка В на рис. 2, а) приводит к его возрастанию над угловыми участками лавы (точки А и С на рис. 2, а).
В процессе движения очистного забоя на расстоянии равном 0,7 длины лавы резко возросла
• сопряжение со штреком (точка А)
- * - - середина лавы (точка В)
:Я
К
Я
Л
с
ей
Е
X
3
Е
Л
§
И
Е
Н
Л
Й
Е
Е
3
£
я
Я
я
о
и
н
я
я
я
я
•е
•е
т
о
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150
Расстояние от монтажной камеры, м
Рис. 5. Графики изменения коэффициента концентрации вертикальных напряжений при движении очистного забоя
нагрузка на угольный пласт в средней части лавы. На следующем шаге произошло первичное обрушение пород основной кровли, что привело к разгрузке и снижению коэффициента концентрации вертикальных напряжений до 0,87 в призабойной части угольного массива в середине лавы (рис. 5).
После частичного обрушения пород основной кровли в середине выработанного пространства произошло перераспределение веса необрушенных пород на краевые участки угольного пласта, т.е. за счет взаимного влияния линии забоя и границ выемочного участка в угловых частях выработанного пространства кровля не обрушилась, и резко возросло горное давление. Это привело к увеличению коэффициента концентрации вертикальных напряжений в пласте до 3,22 вблизи сопряжения вентиляционного штрека и линии забоя (рис. 5).
При следующем подвигании забоя на 10 м произошло полное обрушение пород кровли в пределах всего отработанного контура. Коэффициент концентрации вертикальных напряжений в угольном пласте вблизи сопряжений с горными выработками составил 1,48 (рис. 5).
В процессе дальнейшего подвигания лавы снова начала формироваться консоль зависания пород кровли. Коэффициент концентрации вертикальных напряжений в угольном массиве увеличивается как в середине лавы, так и на верхнем и нижнем сопряжениях лавы с горными выработками. При движении забоя происходит частичное обрушение пород непосредственной кровли, что приводит к снижению коэффициента концентрации вертикальных напряжений.
Через 30 м после первичного обрушения коэффициент концентрации вертикальных напряжений увеличился примерно в 1,3 раза, и произошло вторичное обрушение. Вторичному обрушению соответствует резкое уменьшение коэффициента концентрации вертикальных напряжений в 1,4 раза (рис. 5). Длина консоли вторичного зависания составила 0,33 от длины первичного шага обрушения.
Ввиду высокой аварийности горно-шахтного оборудования возможны остановки очистных забоев продолжительностью до нескольких месяцев. В этом случае вследствие высоких напряжений в угольном пласте вблизи от забоя происходит их релаксация, т.е. снижение коэффициента концентрации вертикальных напряжений (рис. 6).
Релаксация напряжений сопровождается уменьшением прочности угля, что приводит к росту смещений пород кровли над очистным забоем. При этом вследствие деформации пород кровли и пучения пород почвы происходит конвергенция перекрытия и основания секций крепи, "посадка" стоек механизированной крепи насухо и возникновение аварийной ситуации.
По результатам расчетов также установлено, что при остановке очистного забоя происходит увеличение высоты зоны обрушения пород кровли, в том числе с образованием трещин и заколов впереди очистного забоя.
Таким образом, полученные результаты удовлетворительно описывают циклический характер зависания и обрушения подработанных пород кровли при движении очистного забоя. Отношение вторичного и первичного шагов обрушения попадают в диапазон 0,25.0,5, установленный экспериментально [1-5].
Изменение вычисленного коэффициента концентрации вертикальных напряжений в диапазоне 1,3.2,15 соответствует экспериментальным данным, диапазон значений которых составляет 1,4.2,29 по замерам напряжений на пластах впереди очистных и подготовительных выработок гидравлическими датчиками, выполняемым на протяжении ряда лет в ВостНИИ [7].
Разработанная авторами пространственная гео-механическая модель деформирования и разрушения горных пород позволяет прогнозировать предельное время остановки очистного забоя в случае аварийной ситуации с учетом раздвижности секций механизированной крепи.
Рис. 6. Графики изменения коэффициента концентрации вертикальных напряжений при остановке очистного забоя
СПИСОК ЛИТЕРАТУРЫ
1. Кузнецов ГН. Механические свойства горных пород. —М.: Уг-летехиздат, 1947. —180 с.
2. Черняк И.Л., Ярунин С.А. Управление состоянием массива горных пород. — М.: Недра, 1995. — 395 с.
3. Борисов А.А. Механика горных пород и массивов. — М.: Недра, 1980. —360 с.
4. Борисов А.А. Расчеты горного давления в лавах пологих пластов. — М.: Недра, 1964. —278 с.
5. Управление горным давлением при разработке пологих пластов с труднообрушаемой кровлей на шахтах Кузбасса /
С.И. Калинин, А.Ф. Лютенко, В.П. Егоров, С.Г Дьяконов. — Кемерово: Кемеровское книжное изд-во, 1991. —247 с.
6. Павлова Л.Д., Фрянов В.Н. Моделирование циклического характера обрушения горных пород при проведении выработки с последовательным накоплением повреждений // Известия Томского политехнического университета. — 2004. — Т. 307. — № 2. —С. 76-79.
7. Чернов О.И., Мурашев В.И., Шлиомовичус Я.Г. Распределение напряжений вокруг выработок с учетом псевдопластичес-ких свойств угольного массива // Проблемы механики горных пород: Матер. Всес. научной конф. по механике горных пород.
— Новосибирск: Наука, 1971. — С. 684-689.
УДК 550.42:57.4(556.06)
ОЦЕНКА ВЛИЯНИЯ СБРОСОВ СТОЧНЫХ ВОД НА МИНЕРАЛИЗАЦИЮ И ОБЩЕЕ СОДЕРЖАНИЕ ОРГАНИЧЕСКИХ ВЕЩЕСТВ В ВОДАХ ТОМИ
О.Г. Савичев
Томский политехнический университет E-mail: OSavichev@mail.ru
Приведены обобщенные данные о суммарном содержании главных ионов (С) и органических веществ (по углероду Сорг) в водах р. Томи в зимнюю межень и в сточных водах, поступающих в реку Впервые выполнен бассейновый расчет изменения величин Си и Сорг в водах р. Томи под влиянием ее притоков и наиболее крупных выпусков сточных вод на участке от г. Междурече-нска до устья. Показано, что заметное влияние сбросов сточных вод по наиболее крупным сосредоточенным выпускам на минерализацию и общее содержание органических веществ в р. Томи прослеживается на участках протяженностью до 8...10 км. В то же время выявлено несоответствие расчетных и измеренных значений Си и Cv, свидетельствующее о наличии эффективных механизмов самоочищения речных вод.
Введение
Томь — один из крупнейших притоков Оби. Ее сток формируется на территории республик Хакасия и Алтай, Кемеровской, Томской и Новосибирской областей. Данная река служит важным источником водоснабжения целого ряда населенных пунктов, включая гг. Междуреченск, Новокузнецк, Кемерово, Томск, Северск. Одновременно она является приемником значительного объема сточных вод многочисленных промышленных, сельскохозяйственных и жилищно-коммунальных предприятий, расположенных в пределах водосборного бассейна преимущественно в Кемеровской и Томской областях, что обусловливает актуальность оценки суммарного влияния сбросов сточных вод на химический состав вод р. Томи.
Необходимость такой оценки назрела достаточно давно, но до сих пор по ряду причин соответствующие работы выполнялись только в процессе моделирования распространения загрязняющих веществ в проектируемом Крапивинском водохранилище в среднем течении р. Томи и для одиночных источников при разработке нормативов предельно допустимых сбросов. Учитывая требуемое для расчетов огромное количество водохозяйственной и гидрологической информации, ниже рассмотрено влияние лишь наиболее крупных восемь-надцати выпусков с расходами сточных вод более
5000 м3/сут. Кроме того, учитывалось поступление вод шести наиболее крупных (рр. Уса, Мрас-Су, Кондома, Верхняя, Средняя и Нижняя Терси, Тай-дон) и пяти самых загрязненных притоков р. Томи (рр. Аба, Ускат, Искитимка, Ушайка, Киргизка). В качестве исходных данных по составу и объемам сточных вод использованы обобщенные сведения статистической отчетности 2-тп(водхоз), приведенные в работах [1-6], а по составу и расходам речных вод — сведения из работы [7].
Методика исследований
Оценка антропогенного влияния проводилась по двум показателям, позволяющим получить общее представление о химическом составе речных вод и условиях его формирования, — сумме главных ионов (Си), почти полностью определяющих минерализацию воды, и содержанию углерода органических соединений (Сорг). Поскольку в сточных водах концентрации ионов кальция, магния, натрия, калия, гидрокарбонат-иона и величина Соргво многих случаях не измеряются, но почти всегда определяются сухой остаток (Я) и биохимическое потребление кислорода (БПКп), сумма главных ионов в сточных водах определялась расчетным путем по полученной автором эмпирической зависимости от величины Я, а содержание Сорг — по зависимости от БПКп:
