CC BY
70
-------------------------------------- © А.В. Адушкин, Л.И. Беляева,
А.И. Гончаров, В.И. Куликов, 2009
УДК 622.831(470.13)
A.В. Адушкин, Л.И. Беляева, А.И. Гончаров,
B.И. Куликов
ГЕОДИНАМИЧЕСКИЙ МОНИТОРИНГ
ПРИ РАЗРАБОТКЕ УГОЛЬНОГО МЕСТОРОЖДЕНИЯ
ВОРКУТСКОЙ МУЛЬДЫ
Материалы посвящены совершенствованию работы станций сейсмического мониторинга SDAS и GITS. Приведены результаты мониторинга геодинамических явлений на пластовом месторождении угля, изложен алгоритм вычисления энергетического класса ГДЯ с учетом затухания сейсмических волн. Представлены местный годограф сейсмических волн и график повторяемости сейсмических событий, рассмотрены его особенности.
Ключевые слова: сейсмический мониторинг, годограф, геодинамические явления(ГДЯ), энергия и класс сейсмического события
Семинар № 3
A. V. Adushkin, L.I. Belyaeva,
A.I. Goncharov, V.I. Kulikov THE GEODYNAMICAL MONITORING DURING COAL MINING AT VORKUTSKAYA MOULD
Materials are devoted to SDAS and GITS seismic monitoring stations perfection. Results of the geodynamic phenomena monitoring on coalfield are resulted. The dynamical events power class calculation algorithm taking into account attenuation of seismic waves is stated. Seismic wave’s local hodograph and the seismic events repeatability graph are presented, their features are considered.
Key words: Seismic monitoring, time path curve, geodynamical occurrences (GDO), energy and class of seismic event.
Т^азработка Воркутского угольного Р месторождения осложнена гео-динамическими явлениями (ГДЯ) в массиве горных пород. Для контроля и предупреждения геодинамических явлений на шахтах «Комсомольская» и «Северная» ОАО «Воркутауголь» развернут разработанный ОАО ВНИМИ комплекс GITS. В них входят по 12 датчиков -
акселеро-метров, размещенных в шахтных полях. Частота опроса комплекса 2 кГц, что позволяет получать кинематические параметры геодинамических явлений и определять очаги ГДЯ с высокой точностью до +15.0 м.
В настоящее время информация от GITS дополняется данными станции геодинамического мониторинга
SDAS, созданной ИДГ РАН. Станция расположена за городом, примерно в в 12-15 км от шахтных полей. Она состоит из измерительного сооружения (бункер глубиной 6 м, состоящий из двух отсеков расположенных друг над другом) и пункта сбора и первичной обработки сейсмической информации. Внутри бункера на бетонном постаменте установлены 6 сейсмоприемников. Станция оснащена тремя короткопериодными сейсмоприемниками СМ-3КВ и тремя широкополостными СМ3-ОС, блоком сбора и выделения сейсмического сигнала SDAS, системой точного времени GPS, персональным компьютером для
архивации и визуализации зарегистрированной информации и системы резервного питания.
Частотный диапазон короткопериодных каналов составляет 0,5 - 30 Гц. Данные с этих датчиков регистрируются в памяти компьютера с частотой 100 Гц. Частотный диапазон широкополостного датчика 0,02 - 20 Гц, его данные запоминаются с частотой 20 Гц. Динамический диапазон 96 дБ. Каждый канал регистрации записывается в память дважды с низким и высоким усилением, что позволяет не пропустить ни сильных, ни слабых сейсмических волн. На каналах высокого усиления чувствительность столь высока, что порог регистрации событий определяется не чувствительностью аппаратуры, а только уровнем естественных микросейсмических шумов в данной местности.
Станция SDAS ведет регистрацию трех компонент колебаний в одной точке. При этом сейсмоприемники адекватно передают колебания земной поверхности, на которой они установлены, так что на сейсмограммах могут быть выделены фазы продольных, поперечных и поверхностных волн. Это позволяет по годографам и поляризации волн рассчитывать очаги сейсмических событий. Это - совершенно отличная от GITS методика расчета координат сейсмического события. Конечно, станция SDAS уступает по точности определения координат по системе GITS. Однако у нее другие преимущества. Первое - она регистрирует и определяет координаты событий по всему Воркутскому угольному месторождению, что не доступно системе GITS. Второе - дает динамические параметры ГДЯ. Станция SDAS позволяет впрямую (без упрощений и предположений) определить энергию сейсмических волн или энергетический класс ГДЯ. Энергетика геодинамиче-
ских процессов в горном массиве является решающей информацией для прогноза катастрофических геодинамических событий. Кроме того, волновые формы позволяют восстановить механизм выделения энергии при ГДЯ (идентифицировать горные удары, внезапные пылегазовые выбросы, технологические взрывы, обрушения в горных выработках и т.д.).
Эксплуатация станции ББАБ показала, что она регистрирует все ГДЯ, которые регистрируются системой ИТБ в шахтных полях «Комсомольская» и «Северная» классом от 3 и выше и вносятся в ее каталог. Также экспериментально было проверено, что она не пропускает заглубленные взрывы массой около 1 кг ВВ на эпицентральном расстоянии 10 км. Такие взрывы относятся к событиям класса К = 3.
Региональный годограф. Определение координат ГДЯ станцией ББАБ подразумевает определение эпицен-трального расстояния от очага ГДЯ и азимута на этот очаг. Эпицентральное расстояние определяется по разности времен прихода продольных, поперечных и поверхностных волн при известном годографе сейсмических волн. Точность определения эпицентрального расстояния зависит от точности построения годографов. Имеются теоретические годографы, которые заложены в систему обработки данных станции БВАБ. Однако для повышения точности расчета координат сейсмического события был получен экспериментальный местный годограф.
Для построения регионального годографа был использован взрыв на карьере Юнь-Яга, который проводился 23 мая
2007 г. Взрыв проводился в блоке угля на горизонте +158 м. Блок был разбурен скважинами длиной от 3 до 10 м и диаметром 250 мм. Количество скважин -
Рис. 1. Массовый взрыв на карьере Юнь-Яга; место регистрации - сейсмостанция 8БЛ8
65. Сетка скважин 7-8 м. Крепость пород 10-12 по Протодьяконову. Общая масса ВВ в блоке составляла 6800 кг. Скважины были объединены в 4-ре группы по 1700 кг ВВ. Группы подрывались последовательно с замедлением 20 мс.
Регистрация сейсмических колебаний проводилась в трех пунктах: в пос. Советском, на станции геодинамическо-го мониторинга и у главного ствола шахты “Комсомольская”. Эпицентраль-ное расстояние от блока до сейсмопункта в пос. Советском составляло 1,287 км, от блока до станции геодинамического мониторинга SDAS составляло 12,18 км и от блока до ствола шахты «Комсомольская» составляло 25,522 км.
В сейсмопунктах в пос. Советском и на шахте «Комсомольская» были установлены мобильные сейсмостанции, в состав которых входили три короткопериодных сейсмоприемника СМ-3КВ. В каждом пункте регистрации на один из каналов записывались минутные метки времени с приемника спутниковой системы навигации GPS, что позволило синхронизовать сейсмограммы со всех пунктов регистрации по времени и оп-
ределить абсолютное время прихода сейсмических волн.
На рис. 1 приведены сейсмограммы взры-ва, зарегистрированные сейсмостанцией ББАБ. Верхние три канала -короткопериодные, нижние три - широкополостные. Амплитуды колебаний по вертикали на короткопериодных каналах составили 0,055 мм/с, по направлению Юг-Север около 0,026 мм/с и по направлению Запад-Восток около 0,04 мм/с. На рисунке четко выделяются вступления продольной волны (метка Р) и поперечной волны (метка Б).
На мобильных сейсмопунктах были получены сейсмограммы аналогичные короткопериодным каналам ББАБ. По полученным значениям времен вступления Р и Б волн были построены времена вступления волн в зависимости от эпи-центрального расстояния, так называемые годографы продольных и поперечных волн, которые показаны на рис. 2. Экспериментальные точки лежат на прямых, которым соответствуют скорость продольных волн 5,5 км/с и скорость поперечных волн - 2,8 км/с. Эти значения отличаются в меньшую сторону почти на 10% от скоростей расчетного регионального годографа. Это является существенной поправкой для обработки данных системы ББАБ и позволит повысить точность определения эпицен-трального расстояния до нескольких процентов.
Под определением местоположения эпицентра сейсмического события или ГДЯ также понимается определение азимута очага. Азимут определяется по
поляризации сейсмических волн. В нашем случае азимут строится благодаря тому, что в регистрации участвуют всегда три разно ориентированных сейсмоприемника. Априори, погрешность определения азимута всегда выше, чем эпицентрального расстояния. Но “пристрелка” системы регистрации по технологическим взрывам позволила существенно уменьшить и эту погрешность.
Азимут очага взрыва на карьере Юнь-Яга, рассчитанный по поляризации волн, и по данным локации очага взрыва с помощью GPS показала удовлетворительную точность системы обработки данных - 10%.
Затухание сейсмических волн. Характер затухания сейсмических волн имеет большое значение для определения энергии источника сейсмических волн. Затухание амплитуды волн связано с геометрической расходимостью фронта волны, с поглощением и рассеянием сейсмических волн. Коэффициенты поглощения и рассеяния определяются свойствами горного массива, его крепостью, трещиноватостью, влагона-сыщенностью. Обычно затухание описывают степенной зависимостью, которая объединяет все три механизма затухания. Эту зависимость пишут в виде:
т, A
V = — (1)
Rn
где V - амплитуда скорости колебаний, Я -эпицентральное расстояние.
Опыт показывает, что в различных горных породах степень затухания п изменяется в пределах от 3 до 1, поэтому важно установить степень затухания сейсмических волн, присущую региону, в котором ведется мониторинг ГДЯ. Для этого вновь воспользовались результатами регистрации сейсмических волн от взрыва на карьере Юнь-Яга.
Сейсмический эффект массовых взрывов (амплитуда скорости сейсмических колебаний) зависит от эпицен-трального расстояния от заряда до объекта (места регистрации) и от массы заряда в группе (на одну ступень замедления). В 2000-2004 г.г. при исследованиях массовых взрывов в железорудных карьерах Курской магнитной аномалии сотрудниками лаборатории экспериментальной геофизики Института динамики геосферы РАН (ИДГ) была получена зависимость максимальной скорости колебаний от приведенного расстояния, которая описывается формулой [1]:
V = А ■
Гп1/зЛп
q___
R
(2)
где V - скорость колебаний в мм/с, q -масса заряда в группе (на одну ступень замедления) в кг, Я - эпицентральное расстояние в м.
Естественно было считать, что аналогичная зависимость справедлива для взрывов в карьере Юнь-Яга. На рис.3 зарегистрированные максимальные скорости колебаний в пос. Советском (точки слева), на станции геодинамического мониторинга ББАБ (точки в центре ри-
Рис. 3. Зависимость амплитуды сейсмовзрывных волн от приведенного эпицентрального расстояния
сунка) и на шахте «Комсомольская» (точки справа) отложены в зависимости от приведенного расстояния, которое равно эпицентральному расстоянию, деленному на кубический корень из массы заряда в одной ступени замедления. Обработка данных по методу наименьших квадратов позволила установить, что максимальные скорости колебаний на рис. 3 описываются зависимостью:
V = 5030•
с і/з у-7 ч
я
(3)
V У
Эта зависимость приведена прямой линией в двойном логарифмическом масштабе. Из нее следует, что степень затухания амплитуды сейсмических волн с расстоянием в Воркутинском регионе п = 1,7. Указанная закономерность получена для расстояний 1 - 25 км и может быть использована для расчета энергии или класса ГДЯ.
Определение энергетического класса ГДЯ. В «большой сейсмологии» укоренилось два типа энергетической классификации землетрясений: по потоку энергии сейсмических волн, пересекающих сферу стандартного радиуса с центром в фокусе землетрясения, и по полной энергии сейсмических волн, излучаемых очагом землетрясения. В первом случае энергию землетрясения характеризуют классом - К, во втором -магнитудой - М. В настоящее время построены корреляционные зависимости между классом и магнитудой и поэтому у сейсмологов нет трудностей и какой-
либо неоднозначности в характеристике землетрясений при выборе типа энергетической классификации.
Однако практика обработки сейсмограмм показала, что первый подход (определение класса) наиболее продуктивен для случая слабых землетрясений, горных ударов, промышленных химических взрывов, пылегазовых взрывов, обрушений горных выработок, схода горных лавин, гидроразрывов, записанных сейсмостанциями на сравнительно небольших расстояниях от очага (до 100 км).
Напомним определение класса [2]. Из общего количества энергии, высвобождаемой в очаге ГДЯ, поддается учету лишь некоторая ее часть, переходящая в энергию упругих волн, или только сейсмическая энергия ГДЯ - Е (Дж). Величину ее характеризуют потоком энергии упругих волн через замкнутую сферическую поверхность с центром в очаге ГДЯ за время Т , в течение которого происходят колебания.
Е =
і
Ц(Р • (4)
5 о
Здесь Р - вектор плотности потока энергии в точках поверхности S, по которой производится интегрирование, п - единичный вектор нормали к площадке ds.
Так как вектор потока Р меняется за счет геометрической расходимости волн и за счет затухания амплитуды, для однозначности определения потока энергии в качестве поверхности S берут ре-ференц-сферу радиуса И = 10 км. Тогда выражение (1) можно переписать в виде
E = Аж ■ Rl I (P ■ n)dt
(5)
Вектор Умова равен
р = р- с • Г02 (6)
где р - плотность горного массива, С -скорость распространения сейсмической волны, ^0) - скорость коле-баний в волне на референц-сфере.
Однако сейсмическая волна от ГДЯ обычно регистрируется на некотором расстоянии Я, которое может быть и ближе и дальше этого расстояния. Поэтому зарегистрированную сейсмограмму следует пересчитать с расстояния Я на расстояние Яз , для чего используется установленный закон затухания амплитуды сейсмической волны (3):
Vo(t) = V (t)
V Ro J
(7)
Подставив (7) и (6) в выражение (5), получим:
/ _ \ 2n ,
> 2
E = Аж ■ R0 ■ р ■ C ■
R
V R0
| V2dt (8)
o
Класс ГДЯ определяется выражени
ем:
К = lg E , Дж
(9)
Е сли
в сейсмопункте ведется регистрация трех компонент колебаний, то интегралы плотности энергии вычисляют для каждой компоненты отдельно и затем результаты вычислений складывают. Если регистрируется только одна компонента колебаний, то в этом случае результат вычисления интеграла следует утроить, так как компоненты колебаний по трем направлениям отличаются незначительно [3].
Если сейсмоприемники расположены на дневной поверхности, как это имеет место на станции геодина-мического мониторинга SDAS, то сейсмические волны приходят на станцию регистрации из глубинных слоев массива и вследствие отражения от свободной поверхности их амплитуды удваиваются. Поэтому в интеграл плотности энергии следует ввести коэффициент 1/4 .
Определение класса рассмотрим на примере ГДЯ, произошедшего 17 июня
2008 г. в 15:08:08 GMT. Это событие было зарегистрировано одновременно системой GITS, станцией SDAS и мобильной сейсмической станцией, установленной на земной поверхности практически над гипоцентром события.
о
Рис. 4. Геодинамическое событие на шахте Комсомольская 17 июня 2008 г. 15:08:08 СЫТ. Система ГИТС, датчик № 1
Рис. 5. ГДЯ на шахте Комсомольская, 17.06.2008. 15:08:08,212 GMT
На рис. 4 приведены акселерограммы (в каунтах) этого события, зарегистрированные датчиками системы ГИТС на ги-поцентральном расстоянии 2,51 км. Класс ГДЯ в системе GITS определяется по дли-
тельности акселерограммы и для данного ГДЯ К = 3,3.
На рис. 5 показаны велосиграммы этого ГДЯ, зарегистрированы станцией ББАБ.
Z
156
N
t
Рис. 6. ГДЯ на шахте Комсомольская.17.06.2008, 15:08:08 GMT; регистрация в эпицентре события
Гипоцентральное расстояние составляло 9,56 км. Вычислениями по формуле (8) был получен класс К = 3,97.
На рис. 6 показаны велосиграммы этого ГДЯ, зарегистрированы мобильной станцией на земной поверхности
над гипоцентром события. Г ипоцен-тральное расстояние составляло 2,03 км. Вычислениями по формуле (8) был получен класс К = 3,9.
Как и следовало ожидать, определение класса путем интегрирования вело-сиграмм дает однозначный результат, независящий от размещения станции регистрации. Определение класса ГДЯ по продолжительности колебаний приводит к меньшему значению класса. Поэтому был проведен сравнительный анализ энергии для 53 событий, полученных системой GITS - е и станцией SDAS - E, который приведен на рис. 7. Из него видно, что для слабых ГДЯ с энергией около 103 Дж энергии в обеих системах
Рис. 7. Зависимость между энергией ГДЯ, определенной по методике ГИТС (e), и SDAS (E)
Зшдашры событий (шик)
Рис. 8. Гипоцентры геодинамических событий по данным GITS в пределах полей шахт «Комсомольская.» (1) и «Северная» (2)
близки, однако для сильных ГДЯ энергии отличаются в несколько раз. Таким образом, определение энергии ГДЯ по продолжительности колебаний, которое применяется до сих пор на некоторых станциях мониторинга, может привести к значительным погрешностям.
Геодинамическая активность шахтного поля. В период опытной эксплуатации системы геодинамического мониторинга GITS зарегистрировано 10500 сейсмособытий, сформирован каталог сейс ми-чес-
кой активности массива. Статистический анализ выявил прямую зависимость сейсмической активности массива шахтного поля от техногенного воздействия на него очистных работ. На рис. 8 показано расположение гипоцентров ГДЯ, которые имеют определенную упорядоченность в пространстве и во времени. Обобщенный анализ циклов горных работ и сейсмичности показал наличие их тесной взаимосвязи. Установлено, что в пределах разрабатываемых участков
158
шахтного поля происходит до 50% всех событий , 33% - на прилегающих к горным работам зонам, которые расположены до 200 м вглубь массива. Остальные 17% событий можно отнести к различным ГДЯ, происходящим в горном массиве вне зоны работ.
Повторяемость геодинамических явлений. Это одна из важнейших характеристик сейсмического режима горного массива, как естественного процесса, так и при техногенных воздействиях на массив. Закон повторяемости или закон Гуттенберга-Рихтера пришел из «большой» сейсмологии [3] и успешно применяется при анализе техногенной сейсмичности. Наклон графика отражает соотношение между числом сильных и слабых ГДЯ, уровень графика - суммарную интенсивность ГДЯ в горном массиве.
На рис. 9 построен график повторяемости ГДЯ по данным каталога ГДЯ за период с 01.09.2007 г. по 01.09.2008 г., который содержит 4247 событий. График можно описать прямой:
^ N = 4,5 - 0,33 • К (10)
При построении графика повторяемости число событий с К = 3 не учитывалось. Очевидно, число слабых событий с К от 2,5 до 3,5 в каталоге занижено из-за того, что система регистрации «пропускает» часть таких событий.
Отметим особенности графика повторяемости. Угол наклона графика повторяемости имеет значение 0,33. Это сравнительно низкое значение - в других регионах и других породных массивах угол наклона в 1,5 - 2 раза выше. Вторая особенность - относительно небольшое число сильных событий с К от 7,5 до 8,5. График повторяемости претерпевает излом. Накопление данных позволит уточнить эти особенности. Предварительно можно сказать, что эти особенности свидетельствуют о низком уровне напряженно-деформированного состояния горного массива (на время проведения анализа) и в нем маловероятны катастрофические ГДЯ класса более 9.
1. Гончаров А.И., Куликов В.И., Минеев
В.И., Седоченко В.В. Сейсмическое действие массовых взрывов на открытых и подземных работах. Сб. трудов ИДГ РАН Динамические процессы во взаимодействующих геосферах. Изд. ГЕОС, 2006, - С. 2233.
2. Раутиан Т.Т. Затухание сейсмических волн и энергия землетрясений. Труды Института сейсмостойкого строительства и сейсмологии АН Таджикской ССР, выпуск 7, 1960, с. 41-66.
3. Касахара К Механика землетрясений.
Мир, 1985, 264 с.
— Коротко об авторах
Адушкин А.В.- кандидат физико-математических наук, ИДГ РАН, [email protected], Беляева Л.И. — главный геофизик-начальник сейсмостанции технической дирекции ФЗАО «Северсталь-Ресурс» в г. Воркута, [email protected],
Куликов В.И. — кандидат физико-математических наук, ИДГ РАН, [email protected]
Гончаров А.И. — кандидат физико-математических наук, ИДГ РАН, [email protected],
----------------------------------- ДИССЕРТАЦИИ
ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИЙ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ
Автор Название работы Специальность Ученая степень
САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ ГОРНЫЙ ИНСТИТУТ им. Г.В. ПЛЕХАНОВА (ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ)
ПИРОЖЕНКО Владимир Петрович Выбор и обоснование рациональных параметров сгустительно-смеситель-ного оборудования закладочных комплексов для приготовления твердеющей смеси на основе хвостов обогащения полиметаллических руд 05.05.06 к.т.н.
