CC BY
42
-------------------------------------------- © В. Сокола-Шевёла, 2009
УДК 550.834:624.044 В. Сокола-Шевёла
МОДЕЛИ ЗАВИСИМОСТИ ИНДУЦИРОВАННОЙ СЕЙСМИЧНОСТИ И НАБЛЮДАЕМЫХ ДЕФОРМАЦИЙ ЗЕМНОЙ ПОВЕРХНОСТИ
Семинар № 1
1. Введение
Горная эксплуатация может стать при опредeлeнных усло-виях причиной появления толчков и горных ударов в выработках. На тер-ритории Верхнесилезского угольного бассейна сейсмическую активность индуцированную выемкой месторож-дения отметили уже в первой поло-вине IX века. Толчки горного мас-сива выступают в настоящее время на более чем неcкoльких девятках угольных шахт. Эпицентры толчков горного массива имеются в несколь-ких районах, характеризующихся сравнительно глубоким залеганием пластов и выступанием в их окру-жающей среде мощных и толстых комплексов песчаников и сильно развитой тектоникой. На изменение состояния угрозы толчками и гор-ными ударами существенное влияние имеет технология ведения гор-ных работ и уровень знаний, позво-ляющий сделать прогноз, а в даль-нейшем применить соответствующую профилактику. Проблема прогнози-рования
энергии толчков и связей между мерами оценки величины сейсмических явлений уже дoвoльнo хорошо разработана. В многих рабо-тах [между другими
2, 10, 11, 13, 14, 15, 16, 19] показана возможность учета деформационных процессов горного массива и земной по-верх-ности для характеристики уровня индуцированной сейсмичности.
Исследования [1, 2, 3, 4, 8, 9, 10, 13, 17] проводимые много лет на кафедре эксплуатации месторождений в Силезском политехническом институте, показали существование корреляционных зависимостей между процессами деформации, происходящими в горном массиве вследствие ведения выемки, и регистрируемым числом и суммарной энергией толчков. В 2004-2007 гг. реализовался исследовательский проект 4Т12А 032 26, основной целью которого было определение зависимости между регистрируемой сейсмической индуцированной активностью, и величинами деформации земной поверхности. Основные результаты этих исследований представлены в работах [5, 6, 7]. В статье представлены, принятые на основе результатов полученных при реализации вышеуказанной работы, модели сейсмической зависимости индуцированной веденной горной эксплуатацией (число N и энергии толчков) и наблюдаемые величины деформации поверхности.
Величины деформации земной поверхности описаны поверхностью поперечного сечения мульды оседания вдоль наблюдательной линии проходящей параллельно направлению подвигающегося фронта лавы. Поверхность принято определять на основе результатов реализуемых периодических геодезических измерений. Модели были проверены для района эксплуатации веденной шахтой каменного угля Халемба
в поле лавы 6, в пласте 415/1. Данные о регистрируемых толчках, сопутствующих веденной эксплуатации, взяты из архивов Шахтной станции геофизики. Величины деформации получены исходя из результатов измерений, выполненных на спроектированной с учетом прогноз деформации, а в дальнейшем внедренной пространственной наблюдательной сети. В результате исследований определено путём оценки принятых параметров моделей, функ-циональные связи индуцированной сейсмичности по отношению к наблюдаемым величинам деформации земной поверхности, описанным поверхностью поперечного сечения мульды оседания вдоль наблюдательной линии, проходящей по направлению подвигающегося фронта лав в районе реализуемых исследований.
Модели представленные в статье могут стать основой для разработки метода прогнозирования изменений сейсмично -сти удароопасных районов.
2. Модели зависимости сейсмичности индуцированной веден-ной эксплуатацией и наблюдае-мых величин деформации земной поверхности
Модели зависимости индуцированной сейсмичности и наблюдаемых величин деформации земной поверх-ности были определены методами статистического анализа. Они позволяют приблизить искомую функциональную зависимость другой более простой функцией, содержащей соответст-
вующие переменные, приближающей искомую функцию в некоторых ограниченных пределах [12].
2.1. Величины деформации земной поверхности определён-ные на основе веденных иссле-дований
Для описания хода и определения значения деформации пунктов земной поверхности под влиянием горной экс-
плуатации, измерению подлежит ряд геометрических величин. В исследованиях принято [18], что существует зависимость между объёмом образованной пустоты в горном массиве, и наблюдаемой сейсмической активностью горного массива, так бесспорным кажется связь между объёмом образующейся муль-
ды, сейсмической эмиссией горного массива, так как объём мульд зависит от объёма образованной пустоты. Одновременно объём наблюдаемой мульды более соответствует реальным движениям горного массива, чем объём образующейся вследствие эксплуатации пустоты в горном массиве.
Принята рабочая гипотеза существования зависимости между объёмом образующейся на поверхности мульды оседания, и наблюдаемой сейсмичностью горного массива.
В принципе при ведении горной эксплуатации на поверхности ведутся периодические геодезические наблюдения на соответственно спроектированных, а затем внедренных наблюдательных линиях. Одна из этих линий находится обычно вдоль подвигающегося фронта лав. Обладаем в основном наблюдениями за оседанием пунктов наблюдательной линии, проходящей вдоль протяженности поля по простиранию, что по необходимости ведет к принятию упрощающей предпосылки, что объём ¥к образующейся мульды меняется во времени пропорционально поверхности поперечного сечения Р„ мульды оседания, образованного вдоль наблюдательной линии. В дальнейших анализах величина РМп принята как указатель деформации земной поверхности, описывающий наблюдаемые суммарные оседания пунктов наблюдательной линии в очередных циклах измерений.
нр пункта
О пункты наблюдательной линии
Были определены также приросты полей поверхности сечения мульд оседания, замеченные в очередных циклах нивелированных измерений АРКІ. Измеренные оседания пунктов земной поверхности замеченные в нескольких выбранных циклах нивелированных измерений, выполняемых в рамках проекта 4Т12А 032 26 представлены на рис. 1.
2.2. Показатели, характеризующие уровень сейсмичности, индуцированные веденными работами принятые в моделях
Как стандарт на шахтных геофизических станциях определяются координаты очагов толчков и сейсмические энергии, оттуда показатели, характеризующие уровень сейсмичности принимаются:
N - число толчков зарегистрированных в данном промежутке времени, Е - суммарную энергию толчков зарегистрированных в данном промежутке времени.
В модели обозначено:
N1 - суммарное нарастающее число толчков зарегистрированных в промежутке времени обозначенным датой выполнения цикла нивелированных измерений; Еі - суммарная нарастающая энергия
ил- циклах измерении
толчков зарегистрированных в данном промежутке времени обозначенным датой выполнения цикла нивелированных измерений.
Конечно, эти показатели зависят от параметров регистрирующей аппарату-ры. Число зарегистрированных случаев зависит от локализации сейсмометров по отношению к активному сейсмиче-ски району и свойств наблюдательной системы. Дополнительно суммарная энергия толчков сильно нелинейна, что имеет связь с выступлением толчков больших энергий. В связи с этим каждый раз следует провести анализ сильнейших сейсмических явлений для эвевентуального исключения случаев ре-гионального характера или являющихся результатом выступления тектонических нарушений.
2.3. Прменяемая модель регрессии В исследованиях шел поиск зависимости между индуцированной сейсмичностью и величиной поля поверхности сечения мульды оседания, описывающего измеренные оседания пунктов наблюдательной линии. Определялась зависимость между нарастающим значением Е (Ы) и Рж
Принята линейная модель регрессии с одной переменной независимой в виде
(для энергии Е)
Е = а + вРм'^ + е, г=1,2...п (1)
где Ег - переменная зависимая - нарастающая сейсмическая энергия, а - свободный член, в - коэффициент регрессии, Рк- переменная независи-мая - нарастающая сумма приростов полей поверхности сечения зарегист-рированных в очередных время - периодах циклов нивелированных измерений, п - число наблюдений, е-случайная переменная, определяемая как случайная составляющая выясняющая разброс значения Ег вокруг ожидаемого значения с нормальным распределением.
Уравнение регрессии решается при предпосылках:
- дисперсия переменной случайной
т\2 2
О принимает постоянное значение с:
й2(е,) =с2;
- случайные составляющие без корреляции:
СОУ(е/, е/) = 0 для 1^)'
Эстимации параметров проведены методом наименьших квадратов.
Эстиматоры параметров а ив линейной модели обозначены а и Ь.
После учета обозначений функцию регрессии из пробы записано в виде Е = а + ЬР№ (2)
Зависимости для нарастающего суммарного числа Ы толчков определены принимая линейную модель регрессии с одной переменной независимой в виде Nj = а + рр№/(- + е / г=1,2...п (3)
где Ыг - переменная зависимая -нарастающее суммарное число толчков.
Остальные как выше.
Уравнение решалось при предпосылках, как для нарастающей суммарной энергии толчков. Окончательно функцию регрессии из пробы записано в виде N = а + Ь^ (4)
Принятые модели были проверены для конкретных горных условий в об-
Рис. 2. Контуры эксплуатации веденной в районе в период охваченный исследованиями и абрис наблюдательной линии 1
ласти веденных исследований деформации земной поверхности в районе эксплуатации лавой 6 в 415/1- КW К „Халемба” в городе Руда Слёнска.
3. Проверка исправности принятых моделей регрессии
Принятые модели регрессии проверено для горных условий охва-ченных объёмом исследовательского проекта. Проведены анализы охватывающие:
- исследование переменных остатков;
- вычисление отклонения составляющей остатков;
- коэффициента детерминации;
- проведено тест гипотеза, который принимает значение ноль для отдельных коэффициентов регрессии;
- вычислено границу предела доверия для отдельных коэффициентов регрессии, их стандартные отклонения.
В разделе презентируются, определенные с использованием моделей, функциональные независимости индуцированной сейсмичности веденной в районе исследований по отношению к распределению измеренных оседаний земной поверхности, описанных поверхностью поперечного сечения Рк мульды оседания. Вдоль наблюдательной линии 1 внедренной в районе исследований. Результаты представлены, в основном, в виде формул и графически в виде диаграмм с разработанными линиями регрессии и обозначенными пределами доверия для линии регрессии и областью доверия для прогноза. Задано уровень доверия 95 %. Анализ охватывал период эксплуатации, характеризующийся выступлением средней и большой сейсмической активности.
3.1. Характеристика района исследований
Исследования деформации земной поверхности были проведены в районе эксплуатации проведенной каменноугольной шахтой „Халемба” лавой 6 в пласте 415/1 на средней глубине 600 м. Пласт относится к группе рудских слоев. В районе эксплуатации падает в южном направлении под углом 7°. В почве пласта залегают: аргиллит мощностью 4,6 м, 8 м слой мелко- и среднезернистого песчаника, 2,3 м слой аргиллита. В кровле залегает: песчаник мощностью 31,8 м, глинистый сланец мощностью 2,4 м. В прошлом в этом районе была ведена эксплуатация пластов: 402, 403/1, 404/2, 405, 407/2, 409, 410, 411, 413/1, 413/2, 414/2, 502, 504, 506, 507, 510.
Очистные работы в области предметной лавы были продолжением эксплуатации веденной в этой партии. Длина лавы составила 285 м, а протяженность поля по простиранию 1070 м. Выёмку ведено разработкой лав с обрушением кровли на высоту до 3,5 м. Из-за пожароопасности применено периодическое дополни-тельное уплотнение
старых вырабо-ток, разработанных обрушением при помощи пыли продуктов возгона или песком. Пласт в этом районе относится к III степени удароопас-ности.
Эксплуатацию в поле лавы начато в декабре 2004 года. Разработку лавы за-
кончено в ноябре 2006 года. Среднее суточное подвигание экс-плуатационного фронта менялось по июль 2005 г. составило около 2,07 м.
Наблюдательная сеть, внедренная для выполнения исследований деформации, состояла из системы двух перпендикулярных наблюдательных линий и рассеянных пунктов. В данной работе анализу подвержены результаты наблюдений выполняемых на измерительной линии 1, проходящей параллельно подвигающемуся фронту лавы. Эта линия состояла из 36 наблюдательных пунктов. Вертикальная съёмка сделана методом прецизионного нивелирования с использованием кодовых нивеляторов с применением инварной рейки, ссылаясь на государственные реперы. Выполнено 101 цикл нивелированных измерений. В анализах, представленных в работе, использовалось 35 циклов. Периодически выполнялись также линейноугловые измерения пунктов наблюдательной сети и пунктов привязки ситуационной государственной основы. Исходя из представленных результатов наблюдений вычислен показатель Рж На рис. 2 представлены контуры эксплуатации веденной в районе в период исследований и абрис наблюдательной сети в пределе используемом в анализах представленных в статье.
Рис. 3. Диаграмма измеренных оседаний наблюдательной линии в период до 25.07.2005
Проверку принятых моделей регрессии выполнено для периода ведения вы-ёмки с декабря 2004 по 25. 07. 2005 года. В этот период установлено максимальное оседание ^ на пункте 25-2206 мм. Лава в это время достигла 425 м протяженности поля по простиранию.
На рис. 3 представлено диаграмму измеренных оседаний ^ пунктов наблюдательной линии 1 в периодах с начала эксплуатации до 25.07.2005 г.
Сейсмическая активность индуцированная веденными работами первоначально была средняя и большая. С августа 2005 года по август 2006 года стабилизировалась, как слабая. В конечный период подвигания лавы вновь отмечен рост активности. В период с декабря 2004 по ноябрь 2006 г. в предметной партии зарегистрировано 339 толчков
энергией порядка 10 3 до 10 6 [Д]. симальная энергия толчкасотрясения составила 3х106 [Д]. В период, подвергаемый анализу, зарегистрировано 228 толчков из чего 216 индуцировали работы в поле наблюдаемой лавы. Контуры эксплуатации с нанесенными очагами зарегистрированных толчков индуцированных веденной эксплуатацией представлены на рис. 4. В табл. 1 представлено зарегистрированное суммарное число и энергию толчков.
Распределение плотности энергии толчков зарегистрированных в пе-риод с
21.12.2004 по 25.07 2007 г. представлено на рис. 5.
3.2. Использование моделей регрессии для описания регистрируемой сейсмической активности в районе исследований
Таблица 1
Составление числа и энергии зарегистрированные толчков индуцированныш веденной эксплуатацией до дня 25.07.2005 года
Энергия Число толчков Суммарная энергия [.I]
103 и 50 2,57 х10 5
104[Л] 161 5,38х106
105 И 5 2,4х106
сумма 216 8,037х106
Таблица 2 Резутт.
Рис. 4. Контуры эксплуатации с нанесенными очагами зарегист-
Исходя из данных о зарегистрированной сейсмической активности и данных, касающихся измерений деформации земной поверхности в виде вычисленного поля поверхности сечения мульды оседания вдоль наблюдательной линии проходящей параллельно подвигающему фронту лавы, проведено эсти-мацию параметров принятых моделей регрессии.
Каждый раз была проведена проверка принятой модели.
Анализ зависимости между переменными Р„ [ Е указал в результате проведенной эстимации параметров модели (1-2) выступление линейной функциональной зависимости в виде (5):
Е = 401321+11664,2 РКЯ = 0,9885 (р<0,05)
(5)
где Е - нарастающая марная сейсмическая
гия [Д], Р„- нарастающая сумма приростов полей поверхности сечения мульд оседания [т2].
Для фиксированной регрессивной зависимости на уровне существенности р<0,05, можно утверждать, что коэффициент корреляции и коэффициент регрессии существенно различны, чем ноль и подразумевается, что переменная Р„ имеет существенное влияние на Е. Значение коэффициента детерминации Е составляет 0,9772. Уравнение регрессии выясняет тем самым больше чем 97 % переменности нарастающей суммарной энергии толчков (табл. 2, рис. 6). Стандартное отклонение остатков составляет 415029. Анализ переменных остатков в функции время - рис. 7, подобранных значений - рис. 8, и
переменной остаткаа в
функции переменной независимой - рис. 9 не внес возражений по отношению к адекватности принятой модели. До-полнительно пункты диаграммы рис. 10 лежат близко прямой уравнения.
Рис. 5. Распределение плотности энергии толчков [J/m2] в период с
21.12.2004 г. по 25.07.2005 г
■ Анализ зависимости между переменными Рк и N обнаружил в результате проведенной эстимации параметров (3-4) между переменными функциональной зависимости в виде N = 2,80916+0,309925 Рк Я=0,9953 (р<0,05) (6)
ростов полей поверхности сечения мульд оседания [т2].
Для определенной регрессивной независимости значение коэффици-ента детерминации Я2 составляет 0,9905, что обозначает, что уравнение регрессии выясняет бoльше, чeм 99 % переменности нарастающего суммарного числа толчков. Коэффициент Ь регрессии суРис. 6. Нарастающая суммарная сейсмическая энергия Е индуцированная в период с 21. 12. 2004 г. в функции нарастающейся суммы приростов полей поверхности сечения мульд оседания Ру-прямая регрессии, кривые доверия
щественно отличается от ля на уровне существенности р> 0,05 (табл. 3, рис. 11). Стандартное отклонение остатков 7,0587.
где N - нарастающее сум-марное число толчков, Р„ -- нарастающая сумма при-
лю
126
(X 1(
15 = 10
І 5 з 0
г
г -5
5
8- -ю
с
-15
40
80 120 160
прогнозируемое N
200
240
. .._7)
Рис. 12. Диаграмма переменных остатков в ФнцииДирерам1м№а11ерштныкт1стш1ти1в в сфммкцшевримтяюлттзателъ активности суммарная сейсмическая энергия
прогнозируемое N
Рис. 13. Диаграмма переменных остатков в Финкl0wДиPобрржьш^енрнвfmЫх- потит/тл ьв рKтицl¡исlttиh*у>mlP|ы^иеl(tшullто.'lни1tае^lтелb
активности суммарная сейсмическая энергия
15
10
<*=
к **
о
і» £ 0
— <и -5
■У
с -10
-15
(X 1 |
¡г с
к §
І і V
200
600
800
400
Р\\ [ш2]
Р\у [ш21
Рис. 14. Диаграмма переменной ос-татка в функции переменной независимой - показатель. актщжрттауммтеж/шстттшков
функции переменной независимой - показатель активности - суммарная сейсмическая энергия
Анализ переменных остатков в функции времени - рис. 12, подобранных значений - рис. 13, пере-менной остатка в функции перемен-ной независимой - рис. 14 не внесли возражений касающихся адекватности принятой модели. Допол-нительное подтверждение правильности корректности принятой модели презентируется в виде диаграммы
- рис. 15, где пункты диаграммы лежат близко прямой уравнения N прогнозируемые = N наблюдаемые.
4. Итоги
В статье представлены модели зависимости сейсмичности индуцированной веденной эксплуатацией и величины измеренных деформаций зем-
240
200
160
120
80
40
0
(
40 80 120 160 200
предусматриваемое значение
240
=У
0 0.2 0.4 0.6 0.8 1
Рис. 15. Диаграмма предусматриваемых N наблюдаемых N
Рис. 10. Диаграмма предусматриваемых значений Е-наблюдаемых Е
ной поверхности. В принятых моде-лях показателями сейсмической ак-тивности принято считать нарастаю-щую суммарную энергию толчков ин-дуцированных выёмкой и растущим числом толчков.
Принято предположение существования зависимости между объёмом образующейся на поверхности мульды
Рис. 11. Нарастающееся суммарное число толчков N в период с 21.12.2004 г. по 25.07.2005 г.. в функции нарастающейся суммы при-ростов полей поверхности сечения мульд оседания Ру-прямая регрессии, кривые доверия
оседания, а наблюдаемой сейсмичностью горного массива. На этой основе показателем деформации зем-ной поверхности принято поверх-ность поперечного сечения мульды оседания вдоль наблюдательной ли-нии, проходящей по направлению подвигающегося фронта лавы. Моде-ли зависимости определено с исполь-зованием методов статистического анализа.
Принятые модели подвержено проверке для области веденных исследований земной деформации поверхности. Использовано их для описания сейсмической активности горного массива в районе исследований.
1. Banka P., Jaworski A.: Wplyw deformacji gorotworu na maksymalne energie zjawisk sejsmicznych induko-wanych eksploataj gornicz^. ZN Pol. Sl. s. Gornictwo z.239. Gliwice 1999.
2. Banka P.: Metoda prognozowania czasowych zmian aktywnosci sejsmicznej w oparciu o deformacje gorotworu wywolane eksploataj gornicza. Praca doktorska. Gliwice 2006.
3. Banka P.: Wplyw deformacji
gorotworu na czasowe zmiany aktywnosci sejsmicznej. ZN Pol. Sl. s. Gornictwo z. 232. Gliwice 1996.
4. Bialek J. Banka P., Jaworski A.: Wykorzystanie analitycznych prognoz napr^zeniowo - deformacyjnych warunkow wybierania do projektowania eksploatacji w rejonach zagrozonych t^paniami. Prace GIG, s. Konferencje nr 26. Katowice 1998.
5. Bialek J., Sokola-Szewiola V., Opalka K.: Ksztaltowanie si§ obnizen punktow powierzchni terenu a rejestrowana aktywnosc sejsmiczna gorotworu. Konfe-ren-cja Ochrona Srodowiska na terenach Gorniczych . Szczyrk 31 V- 2 VI. 2006.
6. Bialek J., Sokola-Szewiola V., Opalka K.: Pr^dkosc obnizen punktow
Полученные путём эстимации параметры моделей позволили определить статистически существенные линейные функциональные связи.
Зависимости определенные для нарастающих величин суммарной сейсмической индуцированной энергии E (числа толчков - N) и значения нарастающей суммы полей поверхности поперечного сечения мульды оседа-ния Pw. Определённые зависимости выясняют свыше чем 97 % перемен-ности нарастающих величин индуци-рованной сейсмичности.
Предлагаемые в статье модели можно использовать для разработки метода прогноза сейсмической опасности.
-------------- СПИСОК ЛИТЕРАТУРЫ
powierzchni terenu a aktywnosc sejsmiczna gorotworu. ZN. Pol.Sl. s. Gornictwo Nr 1752 z. 278. Politechnika Sl^ska. Gliwice 2007.
7. Bialek J., Sokola-Szewiola V., Opalka K.: Subsidence of land area points in relation to recorded rock mass seismic activity. 6 RWTH Aachen. 9-11.11.2006. Aachen. Str.509-518]. Aachen 2006.
8. Bialek J., Drz^zla B., Jaworski A.:
Proba ustalenia zaleznosci funkcyjnych pomi^dzy przebiegiem deformacji gorotworu w czasie a energi^ sejsmiczn^ dla warunkow KWK Rydultowy. Publ.Inst. Geophys. Pol. Acad. Sc.M-16 (245).
Warszawa 1992.
9. Bialek J., Jaworski A.: Sposob wyznaczania stref zagrozonych wzmozon^ aktywnosci^ sejsmiczn^ na podstawie prognozowanych rozkladow wskaznikow deformacji gorotworu. Materialy semina-rium : Ocena stanu zagrozenia t^paniami w kopalniach KGW i srodki ich minimalizacji. ST KG PAN o.Katowice, KGW, GIG. Jaworze 1989.
10. Bialek J., Jaworski A.: Proba oceny aktywnosci sejsmicznej na podstawie prognozowanych stanow deforma-
cyjnych górotworu. ZN AGH, s.Górnictwo z.142. Kraków 1989.
11. Descour J.: Aktywnosc sejsmiczna
i osiadanie powierzchni a mechanizm tipan. Cuprum nr 2. 1979.
12. Draper N.R., Smith H.: Analiza regresji stosowana. PWN, Warszawa 1973.
13. Drzçzla B., Bialek J., Jaworski A., Ban ka P., Kolodziejczyk P.: Badanie zwi^zków sejsmicznosci indukowanej eksploatacji górnic-zi z parametrami opisujicymi deformacje osrodka skalnego. Sprawozdanie kon cowe z realizacji projek-tu badawczego KBN nr 903759101. Gliwice 1994, praca niepub-likowana.
14. Glowacka E.: Ocena zagrozenia sejsmicznego górotworu z uwzglçdnieniem przebiegu eksploatacji. Praca doktorska. Warszawa 1991.
15. Goszcz A.: Wplyw gradientu
prçdkosci obnizania siç powierzchni pod
wplywem robót górniczych na stan zagrozenia wstrz^sami górniczymi. ZN AGH s.Górnictwo z.141. Kraków 1988.
16. Jarmuzek J., Siewierski S,
Wilczynski W.: Wplyw powierzchni
wybranego zloza oraz deformacji terenu na aktywnosc sejsmiczni górotworu na przykladzie ZG Rudna. ZN AGH
s.Górnictwo z.129. Kraków 1987.
17. Jaworski A.: Zwi^zek pomiçdzy deformacji górotworu, a wydatkiem energetycznym wstrzisów z rejonu kopaln Bobrek i Miechowice. ZN Pol. Sl. s. Górnictwo, z. 225. Gliwice 1995.
18. Kijko A.: Theoretical model for relationship between minning seismicity and excaviation area. Acta geoph. Pol. Vol.33. 1985.
19. Wanior J.: Metoda prognozowania wstrzisów i tipniçc w oparciu o wyniki pomiarów geodezyjnych. PTPNOZ. Czçstocho-wa 1982. [ДШ
— Коротко об авторе ------------------------------------------------------------------
Виолетта Сокола-Шевёла - кандидат технических наук, инженер, Силезсий политехниче-сккий институт Гливице, Польша.
Доклад рекомендован к опубликованию семинаром № 1 симпозиума «Неделя горняка-2008». Рецензент д-р техн. наук, проф. А.М. Гальперин
-------------------------------------- ДИССЕРТАЦИИ
ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИЙ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ
Автор Название работы Специальность Ученая степень
ГОСУДАРСТВЕННОЕ УЧРЕЖДЕНИЕ «КУЗБАССКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
ФОМИН Анатолий Иосифович Методологические принципы управления риском профессиональных заболеваний на угольных шахтах Кемеров- 05.26.01 д.т.н.
ской области
