Разработка способа определения физико-механических параметров углепородного массива методом шахтной сейсморазведки тема диссертации и автореферата по ВАК РФ 25.00.20, кандидат технических наук Мартынюк, Александр Рафаэльевич

• Мартынюк, Александр Рафаэльевич
• кандидат технических наук
• 2001

Оглавление диссертации кандидат технических наук Мартынюк, Александр Рафаэльевич

ВВЕДЕНИЕ

1 .ГЕОФИЗИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ ГОРНОГО МАССИВА

1.1 Тензометрические ( механические) методы

1.2 Электрометрический метод

1.3 Электромагнитные методы

1.4 Сейсмические методы

1.5 Датчики напряжений

1.6 Математические модели

1.7 Цель и задачи исследований

2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВОЛНОВОГО ПОЛЯ

2.1 Исследование взаимосвязи параметров волнового поля и пористого тела------------------—--------—

2.2 Связь волновых и реологических параметров среды

2.3 Влияние пористости на физико-механические параметры среды —

2.4 Способы определения релаксационных параметров х3 и в8

2.5 Исследование процесса распространения сдвиговых волн типа

Лява в неидеально упругих средах

2.6. Конечно- разностные схемы волнового уравнения для двумерного случая

2.7 Оценка погрешностей способа—

3.МЕТОДИКА ОПРЕДЕЛЕНИЯ ФИЗИКО-МЕХАНИЧЕСКИХ ПАРАМЕТРОВ УГЛЕПОРОДНОГО МАССИВА

4. ЭКСПЕРИМЕНТАЛЬНОЕ ОПРЕДЕЛЕНИЕ ФИЗИКО-МЕХАНИЧЕСКИХ

ПАРАМЕТРОВ----—

Введение диссертации (часть автореферата) на тему «Разработка способа определения физико-механических параметров углепородного массива методом шахтной сейсморазведки»

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

В последнее время метод шахтной сейсморазведки в модификации сейсмо-просвечивания успешно используется для выявления зон нарушения строения угольного пласта, поскольку этот метод хорошо разработан, обеспечен аппаратурой, позволяющей автоматизировать процесс получения, передачи и обработки сейсмоинформации. Однако, эффективность прогноза все же недостаточна. В настоящее время эта величина не превышает 70%, что объясняется прогнозированием ложных (отсутствующих), либо пропуском реально существующих аномалий. Причиной этого является недостаточная теоретическая изученность структуры волнового поля в горном массиве и использование идеально упругой математической модели, что снижает достоверность интерпретации сейсмограмм. Кроме того, эффективность данного метода снижается за счет экстраполирования точечных значений физико-механических параметров, полученных по образцам в лабораторных условиях, на всю площадь исследуемого участка.

Повысить надежность метода можно двумя путями:

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

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

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

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

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

Научные положения- защищаемые в диссертации:

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

2. Реологические свойства угольного пласта, обусловленные взаимодействием вязких сил различной природы, характеризуются двумя компонентами времени релаксации, описываемые в интервалах частот (до 104 с"1) различными аналитическими зависимостями: установлено, что при частотах (менее 102 с"1) доминирует время релаксации, которое обратно пропорционально коэффициенту поглощения и скорости распространения волны и практически не зависит от частоты, а при частотах (103- 104 с"1) преобладает время релаксации, которое обратно пропорционально квадрату частоты.

3. Выявлено, что время релаксации, определяемое по физикол 1 механическим параметрам при частотах (до 10 с ) уменьшается с ростом общей и минералогической плотности горных пород, плотности жидкого или газообразного заполнителя пор, коэффициентов проницаемости и пористости и увеличивается с повышением вязкости порозаполнителя; при частотах (103 -104 с"1) время релаксации уменьшается с увеличением упругих постоянных Ламе и минералогической плотности и растет с повышением упруго-вязких коэффициентов, коэффициентов пористости и проницаемости пород.

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

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

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

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

О 1 мена релаксации противоположно: при малых частотах (до 10 с ) время релаксации практически не зависит от частоты, а при больших частотах (103-104с"1) -определяется именно величиной частоты;

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

Достоверность разработанных научных положений и выводов подтверждается:

-высокой сходимостью расчетных и экспериментально измеренных в шахте величин физико-механических параметров (коэффициентов пористости, упругих постоянных Ламе и общей плотности угольного пласта). Относительная ошибка не превышает 15%. Эксперимент проведен на пяти шахтах с различными горно-геологическими условиями.

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

Практическая ценность диссертации состоит:

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

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

Реализация работы. Разработанная "Методика определения физико-механических параметров углепородного массива при шахтной сейсморазведке" принята и использована в исследованиях институтов ИПКОН, МГГУ. Основные научные положения и рекомендации по уточнению программ моделирования процесса распространения волнового поля в угленосной толще использованы в экспедиции "Печоруглеразведка" с экономическим эффектом в размере 300.0 тыс. рублей.

Апробация работы. Основные положения и результаты докладывались на научных симпозиумах "Неделя горняка" (Москва, МГГУ,2000 г), (Москва, МГГУ, 2001 г.) и на научных семинарах ННЦ ГП - ИГД им. А.А.Скочинского.

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

Объем работы. Диссертация состоит из введения, трех глав, заключения и содержит 161 страницы машинописного текста, 23 рисунков, 18 таблиц, список литературы из 75 наименований и приложения.

Заключение диссертации по теме «Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика», Мартынюк, Александр Рафаэльевич

Выводы:

В результате обработки экспериментальных сейсмоакустических данных по разработанной методике, полученных коллективами сотрудников отдела горной геофизики УкрНИМИ, экспедиции "Печоруглеразведка", ИГД им. А.А.Скочинского, а также сотрудниками Сианьского филиала Центрального угольного института КНР на пяти шахтах ряда угольных бассейнов ( ш. им. За-сядько ПО "Донецкуголь", ш. Коммунист ПО " Октябрьуголь", ш.Центральная и ш.Воркутинская ПО "Воркутауголь",ш.№ 8 Датунского угольного бассейна КНР), определены комплексы физико-механических параметров для зоны опорного давления и вне ее. К четырем измеренным параметрам (а5,с5,ср,о) дополнительно получено еще 12: коэффициенты пористости и проницаемости, относительная плотность скелета породы и ее общая плотность, время и постоянная релаксации, упругие постоянные Ламе и 3 вязкоупругих коэффициента (для оптимального волнового уравнения).

Было установлено, что часть параметров (с5, со, у, Л, Д,) у бывают по мере удаления от лавы вглубь пласта, имея максимум в зоне опорного давления. Другая часть параметров (а3д,Кп,т3,КпР,^,77)увеличивают свои значения, имея минимум в зоне опорного давления ( ш.Коммунист, ш. им. Засядько). На трех других шахтах ( ш.Центральная и ш.Воркутинская ПО "Воркутауголь" и ш.№ 8 Датунского бассейна КНР) были установлены зависимости указанных параметров от величины относительной пригрузки, что также аналогично удалению пункта измерения от зоны опорного давления. Изменение параметров в этом случае имеет тот же характер, что и в предыдущем.

Эти данные подтверждают аналитические зависимости, полученные в разделе 2 и находятся в согласии с результатами экспериментов других исследователей [28], [69].

Установлено, что количественные значения времени релаксации для угля при значениях коэффициента пористости 5-40% в диапазоне изменения ы= 300

ЗАКЛЮЧЕНИЕ

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

Проведенные исследования дали следующие результаты:

1. Анализ литературных источников показал, что:

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

2. В результате проведенных теоретических исследований:

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

- показано, что реологические свойства угольного пласта, обусловленные взаимодействием вязких сил различной природы, характеризуются двумя временами релаксации, описываемые различными аналитическими зависимостями: при частотах (до 102 с"1) превалирует время релаксации, которое обратно пропорционально коэффициенту поглощения и скорости распространения волны и практически не зависит от частоты, а при частотах (103 -104 с"1) доминирует время релаксации, которое обратно пропорционально квадрату частоты.

- выявлено, что время релаксации, определяемое по физико-механическим параметрам при частотах (до 102 с"1) уменьшается с ростом плотности горных пород и порового флюида, коэффициентов проницаемости и пористости и увеличивается с повышением вязкости флюида; при частотах (103-104 с"1) время релаксации, напротив, растет с повышением коэффициентов пористости и проницаемости и уменьшается с увеличением упругих постоянных Ламе и минералогической плотности.

-установлено, (по полученным зависимостям), что с увеличением частоты возрастаюта$ и с,., причем вначале линейно, а затем эти величины стабилизируются.

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

4. Экспериментально на пяти шахтах (в пяти очистных забоях) с различными горно-геологическими условиями по длине участка протяженностью 70-100м, находящегося впереди очистного забоя, найдено распределение параметров: времени релаксации, коэффициентов пористости и проницаемости, плотности среды, декремента поглощения, упругих постоянных Ламе, и упруговязких

153 коэффициентов. Сопоставление количественных значений коэффициента пористости, плотности, упругих постоянных Ламе с их величинами, полученными прямыми измерениями, показали удовлетворительное совпадение.

5. Экспериментально установлено, что время релаксации для горных пород в5 имеет порядок 0,01-1,00с; а г5-10^-10^ с , что удовлетворительно согласуется с количественными оценками других исследователей.

6. Экономический эффект, полученный в результате внедрения разработанной методики в экспедиции "Печоруглеразведка" составил 300. тыс. рублей (в ценах 2001 г.).

Список литературы диссертационного исследования кандидат технических наук Мартынюк, Александр Рафаэльевич, 2001 год

1. Шкуратник В.Л. Горная геофизика. Ультразвуковые методы.Учебное пособие по дисциплине "Горная геофизика"./ Моск. горный институт. - М., 1990

2. Ардашев К. А., Ахматов В. И., Катков Г. А. Методы и приборы для исследования проявлений горного давления : Справочник . М.: Недра, 1981. - 129 с.

3. Курленя М. В., Попов С.Н. Теоретические основы определения напряжений в горных породах ./ Отв. ред. Шемякин Е.И. Новосибирск: Наука . Сибирское отд-ние, 1983.-97с.

4. Руководство по оценке состояния и свойств угольного массива скважинными гидравлическими датчиками. Новосибирск: ИГД СО АН СССР , 1979.

5. Техника экспериментального определения напряжений в осадочных породах. /Курленя М.В. , Аксенов В.К. ,Леонтьев A.B. и др. Отв. ред. Шемякин Е.И. -Новосибирск: Наука, 1975,- 150с.

6. Грицко Г.И., Кулаков Г.И. Измерение напряжений в горных породах фотоупругими датчиками. /Отв.ред. Шемякин Е.И. Новосибирск: Наука, 1978. - 143 с.

7. Напряженное состояние породных массивов / Отв. ред.Курленя М.В. // Сб. научн. тр. -Новосибирск: ИГД СО АН СССР, 1978. 147 с.

8. Разработка исходных данных к методике расчета исполнительных органов сплошного действия для напряженных горных пород: Отчет / заключительный / ИГТМ; Руководитель темы Полуянский С.Я.; № гос. рег.01830018236. -Днепропетровск, 1983. 58 с.

9. Ю.Лазаревич Л.М. Определение направлений действия главных напряжений методом электрометрии // Геофизические основы контроля напряжений в горных породах Новосибирск: ИГД СО АН СССР, 1983. - С. 15-19.

10. Скакун А.П. Разработка методов оценки степени удароопасности участков массива угля и пород на основе геоэлектрических измерений: Автореф. дис., канд.техн.наук/ВНИМИ. Л., 11985. - 20 с.

11. Проскуряков В.М.,Бляхман A.C. Сейсмические методы исследования напряженного состояния горного массива / Под общ. ред. Петухова И.М. М.: Недра , 1983. - 193 с.

12. Шерифф Р., Гелдарт Л. Сейсморазведка: В 2-х т. Т.1. Пер. с англ. М.: Мир, 1987,448 с.

13. Геофизические методы изучения геологии угольных месторождений.Гречухин В.В, Бродский П.А, Климов A.A. и др. Под ред. Гречухина В.В. / М.: Недра, 1995.477 с.

14. Турчанинов И.А., Панин В.И. Геофизические методы определения и контроля напряжений в массиве. Л.: Наука, 1976. - 224 с.

15. Ямщиков B.C. Методы и средства исследования и контроля горных пород и процессов. Учебник для вузов. М.: Недра, 1982,296 с.

16. Мамбетов Ш.А., Салеев В.М. Определение и контроль диэлектрической проницаемости пород в условиях естественного залегания // Геофизические основы контроля напряжений в горных породах. Новосибирск: ИГД СО АН СССР, 1983.-С. 32-38.

17. Никитина В.Н., Энениггейн Б.С. Определение электрических свойств пород кристаллического фундамента в естественном залегании // Физика Земли. -1968.-№12.

18. Пархоменко Э.И. Электрические свойства горных пород. М: Недра, 1965.

19. Глушко В.Т., Ямщиков B.C., Яланский A.A. Геофизический контроль в угольных шахтах. Киев: Наукова думка. 1978. - 224 с.

20. Ямщиков В. С. Волновые процессы в массиве горных пород: Учебник для вузов. М.: Недра, 1984. - 271 с.

21. Ляховицкий Ф.М. , Хмелевской В.К., Ященко З.Г. Инженерная геофизика. -М.: Недра, 1984.- 252 с ., илл.

22. Проскуряков В. М. Сейсмические исследования напряженного состояния горных пород сейсмическим методом // Физико-технические проблемы разработки полезных ископаемых 1975,- №3. - С. 19-26.

23. Проскуряков В. М. Оценка главных напряжений в массиве горных пород сейсмическим методом // Горный журнал. 1976. - №12. - С. 51-53.

24. Кондратьев О. К. Сейсмические волны в поглощающих средах . М.: Недра . 1986- 176 с.

25. Ватолин Е. С. Об оценке относительной напряженности массива углей и горных пород сейсмоакустическим методом // Научн. сообщения / Ин-т горн, дела им. A.A. Скочинского.- М., 1977. -Вып. 145. С. 126-131.

26. ЗО.Ямщиков B.C. Контроль процессов горного производства. М.: Недра. 1989. -446с.

27. Оценка напряженного состояния угольных пластов с помощью интерференционных волн/ Ю.С. Исаев, А.И. Комаров, С.П. Турлов, М.А. Ильяшов // Геофизические основы контроля напряжений в горных породах . Новосибирск: ИГД СО АН СССР. 1983 . - С. 48-51.

28. Иванова Г.М., Иванов В.В. Управление состоянием призабойной части пласта под контролем звукоулавливающей аппаратуры// Геофизические основы контроля напряжений в горных породах. Новосибирск: ИГД СО АН СССР, 1983. - С.79-80

29. Анциферов М.С., Анцыферова Н.Г., Коган H.H. Сейсмоакустические исследования и проблема прогноза динамических явлений. М.: Наука , 1971. - 135с.

30. Ляховицкий Ф.М., Юдасин Л.А.Влияние порозаполнителя и пластовых условий на скорости и поглощение упругих волн в песчаниках. Изв. АН СССР . Сер. Физика Земли, 1981, №4, с.43-57

31. Иванов Шиц Н.К. Исследование напряженного состояния призабойной части угольного пласта сейсмоакустическим методом // ФТПРПИ. - Новосибирск: Наука, Сибирское отд., 1985. - №3. - С.26-30.

32. Мячник В.И. Ультразвуковые исследования напряженного состояния и свойств горных пород в массиве. Дис. на соис. уч. степ. канд. физ.-мат. наук. М, Фонды ИФЗ АН СССР, 1965.

33. Силаев О.И. Исследование с помощью ультразвука скоростей распространения упругих волн и упругих параметров в образцах горных пород при одностороннем давлении. Труды ИФЗ АН СССР, 1962, №27.

34. Захаров В.Н. Прогнозирование аномально напряженных зон методами шахтной сейсморазведки при подземной разработке угольных месторождений. Автореф. дис. кандидата технических наук / ИГД им. А. А. Скочинского. М., 1995 г.

35. Бондаренко Н.Я., Иваненко Б.Н. Исследование акустической эмиссии на уда-роопасных угольных пластах Донецкого бассейна // Геофизические основы контроля напряжений в горных породах. Новосибирск: ИГД СО АН СССР, 1983. - С.65-67.

36. Горная энциклопедия,- М:. Сов.энциклопедия. Т.5. 1991.-541с., илл. 20 л. илл.

37. Сейсморазведка. Справочник геофизика. / Под редакцией Гурвича И.И. и Номоконова В.П.М.: Недра, 1981.

38. Руководство по применению электрометрического метода оценки нарушенно-сти высокопроводящих пород. Новосибирск: ИГД СО АН СССР, 1979.

39. Рубан А.Д. Физико-технические основы сейсмического мониторинга горного массива для повышения эффективности производства на угольных предприятиях. Автореф. дис. доктора технических наук / ИГД им. A.A. Скочинского. Москва, 1995 г.

40. Григорьев В.Е., Любич Г.А., Сидоров A.A. Промышленная оценка тектонически нарушенных угольных месторождений. Уголь, 1986. - С. 53-54.

41. Будник В.В. Особенности проявления горного давления и выбор параметров механизированных крепей для осваиваемых месторождений Подмосковного угольного бассейна.

42. Свойства горных пород при разных видах и режимах нагружения/ А.И. Берон, Е.С. Ватолин, М.И. Койфман, С.Е. Чирков. М., Недра, 1984.

43. Ватолин Е.С. Влияние скорости приложения нагрузки на прочность горных пород:Науч. сообщения/Ин-т горн, дела им. А.А Скочинского.-1999.-вып.313.с.91-93

44. Присташ В.В. Результаты исследования откольного разрушения горных пород: Науч. сообщения /Йн-т горн, дела им A.A. Скочинского.-1999.-вып.313.с.76-82.

45. Чирков С.Е. Основные зависимости прочностных и деформационных характеристик горных пород от вида напряженного состояния: Науч. сообщения /Ин-т горн, дела им. A.A. Скочинского.-1999.-вып.313.с.47-52.

46. Азаров Н.Я., Яковлев Д.В. Сейсмоакустический метод прогноза горногеологических условий эксплуатации угольных месторождений. М.: Недра, -1988.- 199 с.

47. Николаевский В.Н. Механика пористых и трещиноватых сред. М.: Недра, 1984.

48. Дахнов В.Н. Геофизические методы определения коллекторских свойств и нефтегазонасыщения горных пород. М.: Недра, 1985.

49. Дахнов В.Н. Интерпретация результатов геофизических исследований разрезов скважин. М.: Недра, 1982.

50. Р.Шерифф, Л.Гелдарт. Сейсморазведка в 2-х т. Пер. с англ. М.: Мир, 1987, 400 с.

51. Рейнер М. Деформация и течение. Государственное научно-техническое издательство нефтяной и горнотопливной литературы. М.: 1963.

52. Невский А. Анизотропия скоростей.

53. Френкель Я.И. К теории сейсмических и сейсмоэлектрических явлений во влажной почве. Изв. АН СССР. Сер. геофиз. и географ., 1944, т.8, №4, с. 133 -149

54. Бреховских Л.М. Волны в слоистых средах. М.: Изд-во АН СССР, 1957. -501 с.

55. Biot M.A. Theory of propagation of elastic waves in a fluid saturated porous solid.-part 1,2; - J.Acouest. Soc. Amer., 1956, v.28. p 168-191.

56. Geertsma. The effect of Fluid pressure decline on volumetric Changes of porous Rocks. J. of petroleum technology, December, 1957.63 .I.Tolstoy. Wave propagation. New York, 1980.

57. Математическое моделирование в сейсморазведке / Вербицкий Т.З., Починай-ко Р.С., Стародуб Ю.П. и др. / Киев, Наукова думка, 1985. - 276 с.

58. P.Roy Р and S waves in presstressed solids. Geophysical research bulletin, vol. 25, No 3.,September, 1987.

59. Метлов JI.С. Математическое моделирование волновых полей в шахтной сейсморазведке: Автореф. дис.канд. физ.- мат. наук / МГУ. М., 1987. - 21 с.

60. Азаров Н. Я. Геофизические методы прогнозирования горно геологических условий эксплуатации угольных месторождений Подмосковного бассейна: Автореф. дис. канд. геол. - минер, наук / МГУ . - М. 1978. - 18 с.

61. Петрашень Г.И. Распространение упругих волн в слоисто-изотропных средах , разделенных плоскопараллельными плоскостями // Учен. зап. ЛГУ.- 1952.-№162.-191с.

62. Randall Absorbing boundary condition for the elastic wave equation. Geophysics, vol. 53, No 5. ( May 1988); p. 611-624

63. V.de la Cruz and T.J.T. Spanos. Seismic wave propagation in a porous medium Geophysics, vol. 50, No 10 ( October 1985); p. 1556-1565.

64. L.Knopoff and F.MacDonald Attenuation of Small Amplitude Stress Waves in Solids. Reviews of Modern Physics, vol. 30, No 4 (October 1958); p. 1178-1192.сложение

65. Долевое участие Мартынюка А.Р. в разработке и внедрении «Методики феделения физико-механических параметров углепородного массива при шахтной йсморазведке» составляет 37°о.

66. Расчет составлен в соответствии с «Методикой определения экономической »фективности использования в угольной промышленности новой техники, обретений и рационализаторских предложений». Москва. ЦШТ)Иутоль. 1979 год.

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

68. Долевое участие Мартынюка А.Р. в разработке «Методики определения □ико-механических параметров углепородного массива при шахтной ^сморазведке» и вышеуказанных комплектов программ составляет 37° о.1. Начальник участка ПГФРэкспедиции «Печоруглеразведка»

69. Начальник планово-экономического отдела