Экспериментальное обоснование расчетных параметров и исследование напряженно-деформированного состояния флюидонасыщенных слоистых массивов на основе модели БИО тема диссертации и автореферата по ВАК РФ 25.00.08, кандидат геолого-минералогических наук Артамонова, Нина Брониславовна

• Артамонова, Нина Брониславовна
• кандидат геолого-минералогических наук
• 2004

Оглавление диссертации кандидат геолого-минералогических наук Артамонова, Нина Брониславовна

Введение.

Часть I. Современные представления о математическом моделировании напряженного состояния флюидонасыщенных массивов горных пород.

Глава 1.1. Обзор моделей напряженного состояния флюидонасыщенной породы.

1.1.1. Модель К.Терцаги.

1.1.2. Модель И.Фетта.

1.1.3. Модель М.Био.

1.1.4. Модель Р.Нигматулина.

1.1.5. Модель С.Христиановича.

1.1.6. Модель двойной пористости Г.Баренблатта.

Глава 1.2. Решения модели М.Био о напряженно-деформированном состоянии флюидонасыщенных пород.

1.2.1. Решение Дж.Гамболати.

1.2.2. Решение Д.Хелма.

1.2.3. Решение Х.Чена, Л.Тойфеля и Р.Ли.

1.2.4. Приближенные решения системы уравнений Био в трехосном представлении (Т.Берби, И.Смит, Д.Гриффитс, П.Хси).

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

1.3.1. Полевые геофизические методы.

1.3.2. Полевые гидрогеологические методы.

1.3.3. Полевые инженерно-геологические методы.

1.3.4. Лабораторные геофизические (ультразвуковые) методы.

1.3.5. Лабораторные гидрогеологические методы.

1.3.6. Лабораторные инженерно-геологические методы.

Часть И. Методика определения показателей свойств пород на приборе ИФС-1.

Глава II. 1. Устройство прибора ИФС-1 и методика определения деформационных прочностных и фильтрационных свойств пород.

11.1.1. Основные характеристики и возможности прибора ИФС-1.

II. 1.2. Определение деформационных и прочностных свойств пород в атмосферных условиях.

II. 1.3. Определение деформационных и прочностных свойств пород при действии внешнего гидростатического давления.

II. 1.4. Определение деформационных и прочностных свойств пород при совместном действии гидростатического и порового давлений.

II. 1.5. Определение проницаемости породы при атмосферных условиях.

Глава II.2. Определение физико-механических и фильтрационных свойств различных типов пород на приборе ИФС-1.

11.2.1. Определение свойств известняка.

11.2.2. Определение свойств доломитов.

11.2.3. Определение свойств песчаника.

11.2.4. Определение свойств мела.

Глава II.3. Новый способ определения коэффициента передачи порового давления а.

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

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

11.3.3. Определение коэффициента передачи порового давления а по методу И.Фетта.

11.3.4. Теоретическое обоснование нового метода определения коэффициента а.

11.3.5. Опробование нового метода определения коэффициента а на приборе ИФС-1.

Часть III. Использование численного решения модели М.Био для исследования напряженного состояния конкретных геологических объектов.

Глава III. 1. Теоретические основы нового численного решения модели М.Био.

III. 1.1. Обоснование выбора математической модели М.Био.

III. 1.2. Новое численное решение модели М.Био в осесимметричной постановке с использованием вариационно-разностного метода и метода Холецкого.

III. 1.3. Преимущество численных решений модели М.Био.

Глава III.2. Исследование напряженно-деформированного состояния слоистого массива под Храмом Христа Спасителя в Москве.

III.2.1. Инженерно-геологические условия исследуемого объекта.

Ш.2.2. Построение расчетной модели массива - основания Храма Христа

Спасителя в Москве.

III.2.3. Анализ результатов расчетов изменения напоров и деформации скелета грунта в результате отбора жидкости.

Глава III.3. Изучение напряженно-деформированного состояния нефтяного месторождения Теигиз в Казахстане.

111.3.1. Геологическое строение месторождения Тенгиз.

111.3.2. Гидрогеологические условия месторождения Тенгиз.

111.3.3. Построение расчетной модели месторождения Тенгиз.

111.3.4. Анализ результатов расчетов изменения давления нефти и проседания земной поверхности при эксплуатации месторождения Тенгиз.

Введение диссертации (часть автореферата) на тему «Экспериментальное обоснование расчетных параметров и исследование напряженно-деформированного состояния флюидонасыщенных слоистых массивов на основе модели БИО»

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

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

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

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

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

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

3. Отладка нового численного решения модели Био, разработанного на механико-математическом факультете МГУ проф. С.В.Шешениным.

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

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

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

Объект и методика исследования. Напряженно-деформированное состояние глубокозалегающего нефтяного месторождения Тенгиз в Казахстане и слоистой водонасыщенной толщи пород под Храмом Христа Спасителя в Москве изучалось с помощью нового решения модели Био, основанного на использовании вариационно-разностного метода и метода Холецкого линейной алгебры. Для составления информационного поля параметров, необходимых для расчетов месторождения Тенгиз, автором были обработаны результаты испытаний на стабилометре более 200 образцов (на трехосное сжатие - 116 образцов, на проницаемость - 102 образца), ультразвукового просвечивания 132 образцов, акустического каротажа более 20 скважин. Разрез Тенгизского района был составлен автором по структурным картам, построенным с помощью вертикального сейсмического профилирования. Деформационные и фильтрационные свойства и коэффициент передачи порового давления пород массива под Храмом Христа Спасителя определялись с помощью созданной с участием автора многоцелевой установки высокого давления.

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

Результаты исследований сформулированы в виде следующих защищаемых положений:

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

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

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

4. Установлена при заданном суммарном дебите откачки на любой момент времени прогнозная величина проседания земной поверхности в результате длительной эксплуатации нефтяного месторождения Тенгиз в Казахстане.

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

Апробация работы. Основные положения работы представлены в 8 публикациях (в 5 статьях, 2 тезисах, 1 коллективной монографии) и были изложены на научной конференции МГУ «Новые идеи в инженерной геологии» (1996 г.), международной научной конференции «Engineering Geology and the Environment» (Греция, Афины, 1997 г.), научной конференции аспирантов и молодых ученых МГУ «Проблемы инженерной и экологической геологии» (1998 г.), годичной сессии Научного совета РАН по проблемам геоэкологии, инженерной геологии и гидрогеологии (4-е Сергеевские чтения) (2002г.).

Структура работы. Диссертационная работа состоит из введения, 3 частей, каждая часть из 3 глав, заключения и списка литературы из 105 наименований. Работа изложена на 173 страницах, содержит 27 таблиц, 49 рисунков.

Заключение диссертации по теме «Инженерная геология, мерзлотоведение и грунтоведение», Артамонова, Нина Брониславовна

Результаты исследования прочности и модуля упругости доломитов при одноосном сжатии приведены в табл. II.2.11, II.2.12, II.2.13, II.2.14 в конце параграфа. Как видно из таблиц, деформационные и прочностные свойства более плотного доломита I типа значительно (в 1,5-2 раза) выше, чем у менее плотного доломита II типа. Различия в значениях прочности и модуля упругости образцов внутри каждого типа доломитов объясняются наличием в породе участков повышенной пористости, неравномерного окремнения и ожелезнения, замеченных в шлифах. Для образцов доломитов обоих типов коэффициент Пуассона (v) составляет 0,28-0,30.

Для сравнения результатов были проведены испытания на прочность при одноосном сжатии образцов доломитов на приборе ЦДМ-10/91. Величины прочности доломитов I типа составили 38 - 41 МПа (3 образца), доломитов II типа 15-21 МПа (5 образцов). Полученные значения прочности на приборе ЦДМ-10/91 совпадают с величинами прочности по испытаниям на приборе ИФС-1.

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

Испытания при трехосном сжатии проводились на водонасыщенных образцах, заключенных в герметичную резиновую оболочку, при всесторонних давлениях 2,5 МПа, 5 МП aw 10 МПа.

На рис. И.2.9 в качестве примера приведены графики зависимости относительных вертикальных деформаций от дополнительного напряжения для двух образцов более плотного доломита, испытанных при всестороннем давлении (сгг) 2,5 МПа и 5 МПа, а на рис. II.2.10 показаны графики для трех образцов менее плотного доломита, испытанных при всестороннем давлении (стг) 2,5 МПа, 5 МПа и 10 МПа. Форма графиков при трехосном сжатии доломитов так же, как и для известняков, свидетельствует о появлении пластических деформаций при дополнительных вертикальных напряжениях, близких к разрушающим. Кроме того, можно отметить, что при одном и том же всестороннем давлении величины наибольшего выдержанного образцом вертикального напряжения и наклон графиков для образцов менее плотного и более пористого доломита закономерно

0,00075

0,0015 0,00225 0,003 относительные вертикальные деформации

0,00375

0,0045

Рис. Л.,2 А Графики зависимости относительной вертикальной деформации от дополнительного вертикального напряжения при всестороннем давлении 2,5 МПа (кр.) и 5 МПа (син.) для образцов более плотного доломита Д1.5 и Д1.7 соответственно.

РИС.11.2Д Графики зависимости относительной вертикальной деформации от дополнительного вертикального напряжения при всестороннем давлении 2,5 МПа (кр ), 5 МПа (син.) и 10 МПа (зел.) для образцов менее плотного доломита ДМ.4, ДМ.6, ДМ.9. меньше, чем для образцов более плотного и менее пористого доломита (см. для сравнения графики одного цвета, соответствующие одинаковым значениям всестороннего давления, на рис. II.2.9, II.2.10).

Значения прочности и модуля упругости для образцов доломитов обоих типов сведены в таблицы II.2.11, II.2.12, II.2.13, II.2.14.

По результатам испытаний на одноосное и трехосное сжатие были построены диаграммы Мора (рис. II.2.11 и II.2.12), по которым определяются значения сцепления (С) и внутреннего трения (ф). Для более плотного доломита I типа С = 6 МПа, ср = 57°; для менее плотного доломита II типа С = 4 МПа, ср = 41,5°. Более высокие значения сцепления и угла внутреннего трения получились закономерно больше для менее пористой и более плотной породы I типа.

Кроме прочностных и деформационных свойств, на приборе ИФС-1 определялись фильтрационные свойства доломитов I и II типов, представленные в таблицах II.2.8, II.2.9, II.2.10.

Определение проницаемости образца более плотного доломита Д1.9.

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

4. Вследствие гидравлической изолированности нефтяного резервуара месторождения Тенгиз, первоначально характеризующегося наличием аномально высоких пластовых давлений, в нефтесодержащих известняках при отборе нефти происходит перманентное падение порового давления и рост эффективных напряжений. Через 10 лет после начала эксплуатации месторождения при заданных условиях откачки пластовое давление уменьшится примерно на 45 МПа, вертикальные напряжения в скелете возрастут на 40-45 МПа. Расчеты по модели Био дают возможность получить прогнозную величину проседания земной поверхности в результате откачки флюида при заданном суммарном дебите откачки на любой момент времени. Так, при отборе нефти из месторождения Тенгиз в Казахстане с суммарным расходом 100000 м/сут через 10 лет эксплуатации величина опускания земной поверхности составит около 2 м. Прогнозный расчет проседания земли необходим для правильной организации защитных мероприятий при длительной эксплуатации нефтяных месторождений.

Список литературы диссертационного исследования кандидат геолого-минералогических наук Артамонова, Нина Брониславовна, 2004 год

1. Аникеев А.В., Артамонова Н.Б., Калинин Э.В. Некоторые особенности деформирования и разрушения массивов горных пород при техногенном изменении режима подземных вод.// Геоэкология, № 3,2000. С.249-256.

2. Артамонова Н.Б. Новые идеи в развитии методики информационного обеспечения расчетов напряженного состояния флюидонасыщеиных массивов горных пород.// Новые идеи в инженерной геологии. Тр. научн. конф. М.: МГУ, 1996. - С.123.

3. Байдюк Б.В. Механические свойства горных пород при высоких давлениях и температурах. М.: Гостоптехиздат, 1963. - 102 с.

4. Баренблатт Г.И., Битов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. М.: Недра, 1984. - 208 с.

5. Беликов Б.П., Александров К.С., Рыжова Т.В. Упругие свойства породообразующих минералов и горных пород. М.: Наука, 1970. - 275 с.

6. Беликов Б.П., Александров К.С., Рыжова Т.В. Упругие константы породообразующих минералов и метод вычисления упругих параметров горных пород по минеральному составу.// Инженерная геология в государственном планировании. М.: Наука, 1968.

7. Бушинский Г.И. Литология меловых отложений Днепровско-Донецкой впадины.// Труды института геологических наук АН СССР М., 1954. - вып. 156. - с.306.

8. Гордеев П.В., Шемелина В.А., Шулякова O.K. Гидрогеология. М.: Высшая школа, 1990.-448 с.

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

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

11. Деформационные свойства горных пород при высоких давлениях и температурах./ Л.А.Шрейнер и др. М.: Недра, 1968. - 358 с.

12. Добрынин В.М. Физические свойства нефтегазовых коллекторов в глубоких скважинах. М.: Недра, 1965. - 163 с.

13. Добрынин В.М. Деформации и изменения физических свойств коллекторов нефти и газа. М.: Недра, 1970. - 239 с.

14. Зиангиров Р.С., Роот П.Э., Филимонов С.Д. Практикум по механике грунтов./ М.: МГУ, 1984.- 152 с.

15. Знаменский В.В. Общий курс полевой геофизики. М.: Недра, 1989. - 520 с.

16. Интерпретация данных сейсморазведки. Справочник./ Под ред. О.А.Потапова. М.: Недра, 1990.-448 с.

17. Калинин Э.В., Панасьян Л.Л., Широков В.Н., Артамонова Н.Б., Фоменко И.К. Моделирование полей напряжений в инженерно-геологических массивах. М.: МГУ, 2003. - 262 с.

18. Калинко М.К. Методика исследования коллекторских свойств кернов. М.: Гостоптехиздат, 1963. - 224 с.

19. Кунщиков Б.К., Кунщикова М.К. Общий курс геофизических методов поисков и разведки месторождений полезных ископаемых. М.: Недра, 1972. - 288 с.

20. Марморштейн Л.М. Петрофизические свойства осадочных пород при высоких давлениях и температурах. М.: Недра, 1985. - 189 с.

21. Марморштейн Л.М. Коллекторские и экранирующие свойства осадочных пород при различных термодинамических условиях. Л.: Недра, 1975. - 160 с.

22. Методическое пособие по инженерно-геологическому изучению горных пород. Т.1./ Под ред. Е.М.Сергеева, С.Н.Максимова, Г.М.Березкиной. М.: МГУ, 1968. - 347 с.

23. Методическое пособие по инженерно-геологическому изучению горных пород. Т.2./ Под ред. Е.М.Сергеева, С.Н.Максимова, Г.М.Березкиной. М.: МГУ, 1968. - 370 с.

24. Методическое пособие по инженерно-геологическому изучению горных пород. Т.1./ Под ред. Е.М.Сергеева. М.: Недра, 1984. - 423 с.

25. Методическое пособие по инженерно-геологическому изучению горных пород. Т.2./ Под ред. Е.М.Сергеева. М.: Недра, 1984. - 438 с.

26. Механика горных пород применительно к проблемам разведки и добычи нефти./ Под ред. В.Мори и Д.Фурментро. М.: Мир, 1994. - 416 с.

27. Мироненко В.А., Шестаков В.М. Основы гидрогеомеханики. М.: Недра, 1974.296 с.

28. Нигматулин Р.И. Основы механики гетерогенных сред. М.: Наука, 1978. - 336 с.

29. Никитин В.Н. Основы инженерной сейсмики. М.: МГУ, 1981. - 176 с.

30. Ортега Дж. Введение в параллельные и векторные методы решения линейных систем. М.: Мир, 1991. - 367 с.

31. Осипов В.И. и др. Москва: геология и город. М.: АО «Московские учебники и картолитография», 1997. - 400 с.

32. Отчет по программе: «Изучение геолого-физического строения и создание эффективной разработки месторождения Тенгиз» (Договор 88.00401295 «Гипровостокнефть», научн. рук. проф. А.А.Аксенов), Москва, 1988.

33. Отчет по программе: «Изучение геолого-физического строения и создание эффективной разработки месторождения Тенгиз» (Договор 88.00401295 «Гипровостокнефть», научн. рук. проф. А.А.Аксенов), Москва, 1989.

34. Отчет по программе: «Изучение особенностей строения и физических свойств соленосной толщи Тенгизского месторождения» (Договор 88.00401295 «Гипровостокнефть», научн. рук. проф. А.Л.Аксенов), Москва, 1990.

35. Павлова Н.Н. Деформационные и коллекторские свойства горных пород. М.: Недра, 1975.-240 с.

36. Петкевич Г.И., Вербицкий Т.З. Исследование упругих свойств пористых геологических сред, содержащих жидкости. Киев: Наукова думка, 1965. - 76 с.

37. Победря Б.Е., Шешенин С.В., Холматов Т. Задача в напряжениях. Ташкент: ФАН, 1988.- 197 с.

38. Полевые методы гидрогеологических, инженерно-геологических,геокриологических, инженерно-геофизических и эколого-геологических исследований./ Под ред. В.А.Королева, Г.И.Гордеевой, С.О.Гриневского, В.А.Богословского. М.: МГУ, 2000. - 352 с.

39. Саваренский Ф.П. Инженерная геология. M.-JL: ОНТИ, гл. ред. горно-топл. и геол.-разв. лит., 1937. - 423 с.

40. Савич А.И., Коптев В.И., Никитин В.Н., Ященко З.Г. Сейсмоакустические методы изучения массивов скальных пород. М.: Недра, 1969. - 240 с.

41. Савич А.И. Определение упругих и фильтрационных свойств скальных пород при помощи ультразвука. Использование звука и ультразвука в сейсмологии, сейсморазведке и горном деле. Геоакустика. М.: Наука, 1966.

42. Савич А.И., Ященко З.Г. Исследование упругих и деформационных свойств горных пород сейсмоакустическими методами. М.: Недра, 1979. - 214 с.

43. Сергеев Е.М. Грунтоведение. М.: Изд-во МГУ, 1983. - 392 с.

44. Сергеев Е.М., Сидорова Г.А. К вопросу о составе и свойствах меловых толщ Воронежской области.// Вестник МГУ. Серия физ.-мат. и естест. наук. М., 1950. - №12. -с.133-142.

45. Справочник (кадастр) физических свойств горных пород./ Под. ред. Н.В.Мелышкова, В.В.Ржевского, М.М.Протодьяконова. М.: Недра, 1975. - 279 с.

46. Справочник по физическим свойствам минералов и горных пород при высоких термодинамических параметрах./ Под ред. М.П.Воларовича. М.: Недра, 1978. - 237 с.

47. Справочник физических констант горных пород./ Под ред. С.Кларка мл. М.: Мир, 1969.-543 с.

48. Ставрогин А.Н., Протосеня А.Г. Прочность горных пород и устойчивость выработок на больших глубинах. М.: Недра, 1985. - 271 с.

49. Таран А.С. О возможности сооружения водохранилищ на меловом основании в долине р.Харьков.// Записки научн.-исслед. института геологии Харьковск. гос. управ., 1938.- т.6.

50. Терцаги К. Теория механики грунтов. М.: Госстройиздат, 1961. - 507 с.

51. Троицкая М.Н. Пособие к лабораторным работам по механике грунтов. М.: МГУ, 1961.-304 с.

52. Физические свойства горных пород и полезных ископаемых (петрофизика)./ Под ред. Н.Б.Дортман. М.: Недра, 1984. - 455 с.

53. Христианович С.А. Механика горных пород и горное давление.// Физико-технические проблемы разработки полезных ископаемых. 1985. - №1. - С.3-17.

54. Шаумян J1.B. Физико-механические свойства массивов скальных горных пород. М.: Наука, 1972.- 119 с.

55. Шестаков В.М. Гидрогеодинамика. М.: МГУ, 1995. - 368 с.

56. Шестаков В.М., Невечеря И.К. Теория и методы интерпретации опытных откачек. -М.: МГУ, 1998. 160 с.

57. Addis М.А. The stress-depletion response of reservoirs.// Proceeding of 1997 SPE Annual Technical Conference and Exhibition. Part 1. Formation Evalution and Reservoir Geology. 5-8 October, 1997. San Antonio, Texas. P.55-65.

58. Andreason G.E., Brookhart J.W. Reverse water-level fluctuations.// Methods of Collecting and Interpreting Ground-Water Data. Edited by R.Bentall. U.S. Geological Survey Water-Supply Paper 1544-H. 1963. - P.30-35.

59. Barenblatt G.I., Zheltow I.P., Kochina T.N. Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks.//Journal of Applied Mathematics and Mechanics. 1960. -Vol.24.-P. 1286-1303.

60. Biot M.A. General theory of three-dimensional consolidation.// Journal of Applied Physics. 1941.-Vol.12.-P.155-164.

61. Biot M.A. Theory of elasticity and consolidation for a porous anisotropic solid.// Journal of Applied Physics. 1955. - Vol.26. - No.2. - P. 182-185.

62. Biot M.A. General solutions of the equations of elasticity and consolidation for a porous material.// Journal of Applied Mechanics, Trans. ASME. 1956. -Vol.23. - No.l - P.91-96.

63. Biot M.A., Willis D.G. The elastic coefficients of the theory of consolidation.// Journal of Applied Mechanics, Trans. ASME. 1957. - Vol.24. - No.4. - P.594-601.

64. Biot M.A. Mechanics of deformation and acoustic propagation in porous media.// Journal of Applied Physics. 1962. - Vol.33. - No.4. - P. 1482-1498.

65. Burbey T.J. Validity of Jacob's assumptions for calculating subsidence due to pumping of confined aquifers./ Morel-Seytoux HJ (ed) Proc. 17th Annual Am. Geophys. Union Hyd. Days Conf., 14-17 Apr. 1997, Colorado State Univ., Fort Collins, Colorado.

66. Burbey T.J. Effects of horizontal strain in estimating specific storage and compaction in confined and leaky aquifer systems.// Hydrogeology Journal. 1999. - Vol.7. - No.6. - P.521-532.

67. Chen H. Y., Teufel L.W., Lee R. Coupled fluid flow and geomechanics in reservoir study. -1.Theory and governing equations.// Paper SPE 30752 presented at the Annual Technical Conference & Exhibition. Dallas, October 22-25, 1995.

68. Detournay Б., Cheng Л.Н. Poroelastic solution of a plain strain point displacement discontinuity.//Journal of Applied Mechanics. 1987. - Vol.54 - P.183-187.

69. Detournay E., Cheng A.H. Fundamentals of poroelasticity.// Comprehensive Rock Engineering Principles, Practice and Projects (ed. by J.A.Hudson). Oxford: Pergamon Press, 1993. - Vol.2. - Chap.5. - P.l 13-171.

70. Fatt I. Compressibility of sandstones at low to moderate pressures.// Bulletin of the American Association of Petroleum Geologist. 1958. - Vol.42. - No.8. - P. 1924-1957.

71. Ferris J.G., Knowles D.B., Brown R.H., Stallman R.W. Theory of aquifer tests.// U.S. Geological Survey Water-Supply Paper 1536-E. 1962.

72. Gambolati G. Numerical models in land subsidence control.// Computer Methods in Applied Mechanics and Engineering. 1975. - No.5. - P.227-237.

73. Geertsma J. The effect of fluid pressure decline on volume changes of porous.// Trans. Amer. Inst. Min. Met. Engrs. 1957. - Vol.210. - P.331-340.

74. Helm D.C. One-dimensional simulation of aquifer system compaction near Pixley, California, 2, Stress-dependent parameters.// Water Resourses Research. 1976. - Vol.12. - P.375-391.

75. Helm D.C. Field-based computational techniques for predicting subsidence due to fluid withdrawal.// Man-Induced Land Subsidence (ed. by T.L.Holzer). 1984. - P. 1-22. - Reviews in Engineering Geology 6, Geological Society of America.

76. Helm D.C. Three-dimensional consolidation theory in terms of the velocity of solids.// Geotechnique. 1987. - Vol.37. - P.369-392.

77. Hsieh P.A. Deformation-induced changes in hydraulic head during ground-water withdrawal.// Ground Water. 1996. - Vol.34. - P. 1082-1089.

78. Hubbert M.K., Ruber W.W. Role of fluid pressure in mechanics of overthrust faulting.// Bull. Geol. Soc. Amer. 1959. - Vol.70.

79. Hubbert M.K., Ruber W.W. Role of fluid pressure in mechanics of overthrust faulting: a replay.// Bull. Geol. Soc. Amer. 1960. - Vol.71.

80. Jacob C.E. On the flow of water in an elastic artesian aquifer.// Trans. Am. Geophys. Union. 1940. - Vol.21. - P.574-586.

81. Jacob C.E. Flow in groundwater.// Rouse h. engineering hydraulics. New York: Wiley, 1950. - P.321-386.

82. McDonald M.G., Harbaugh A.W. A modular three-dimensional finite-difference groundwater flow model.// Tech. Water Resour. Invest., US Geol. Surv. 1988. - Book 6. - Chap. AI.

83. Nguyen T.Q., Helm D.C. Land subsidence due to groundwater withdrawal in Hanoi, Vietnam.// Land Subsidence (Proceedings of the Fifth International Symposium on Land Subsidence, The Hague, October 1995). 1995. - Vol.1. - P.55-60.

84. Nur A., Byerlee J.D. An exact effective stress law for elastic deformation of rock with fluids.//Journal of Geophysical Research. 1971. - Vol.76. - No.26. - P.6414-6419.

85. Ochs D.E., Chen H.Y., Teufel L.W. Relating in sity stresses and transient pressure testing for a fractured well.// Paper SPE 38674 presented at the Annual Technical Conference & Exhibition. San Antonio, October 5-8, 1997. - P.301-316.

86. Poland J.F. Subsidence in United States due to ground-water withdrawal.// Journal of the Irrigation & Drainage Division, Proceedings of the American Society of Civil Engineers. 1981. -Vol.107. -P.l 15-135.

87. Poland J.F. Mechanics of land subsidence due to fluid withdrawal.// Guidebook to Studies of Land Subsidence due to Ground-Water Withdrawal (ed. by J.F.Poland). Paris, 1984. - P.37-54.

88. Skempton A.W. The pore-pressure coefficients A and B.// Geotechnique. 1954. - Vol.4. -No.4. - P.143-147.

89. Skempton A.W. Terzaghi's discovery of effective stress.// L.Bjerrum, A.Casagrande, R.B.Peck, A.W.Skempton. From theory to practice in soil mechanics: selections from the writings of Karl Terzaghi. New York: John Wiley, 1960. - P.42-53.

90. Smith I.M., Griffiths D.V. Programming the finite element method. New York: John Wiley, 1988.

91. Terzaghi K. Erdbaumechanik auf bodenphysikalischer grundlage. Deuticke, Wein, Leipzig, 1925.

92. Theis C.V. The relation between lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage.// EOS Trans. AGU. 1935. - Vol.16. -P.519-524.

93. Tuncay K., Corapcioglu M.Y. Effective stress principle for saturated fractured porous media.// Water Resourses Research. 1995. - Vol.31. - No. 12. - P.3103-3106.

94. Verruijt A. Elastic storage of aquifers.// Flow Through Porous Media. Edited by R.J.M. De Wiest. New York: Academic Press, 1969. - P.331-376.

95. Wolff R.G. Field and laboratory determination of the hydraulic diffusivity of a confining bed.// Water Resources Research. 1970. - Vol.6. - P. 194-203.

96. Wolff R.G. Relationship between horizontal strain near a well and reverse water level fluctuation.// Water Resources Research. 1970. - Vol.6. - P.1721-1728.