Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа

Ашихмин, Сергей Геннадьевич

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

Автореферат диссертации по теме "Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа"

На правах рукописи

Ашихмин Сергей Геннадьевич

НАУЧНЫЕ ОСНОВЫ МЕТОДОВ ПРОГНОЗА НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ГОРНЫХ ПОРОД ПРИ РАЗРАБОТКЕ МЕСТОРОЖДЕНИЙ НЕФТИ И ГАЗА

Специальность 25 00 20 - «Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика»

Автореферат

диссертации на соискание ученой степени доктора технических наук

□□344598В

1 8 СЕН 2008

Пермь - 2008

003445986

Работа выполнена в Пермском государственном техническом университете

Научный консультант

доктор технических наук, профессор Кашников Юрий Александрович

Официальные оппоненты

доктор технических наук, профессор Барях Александр Абрамович

доктор технических наук, профессор Сашурин Анатолий Дмитриевич

доктор технических наук, профессор Гордеев Виктор Александрович

Ведущая организация

Защита состоится « 3 i

ООО «ПермНИПИнефть»

» Olil/vyihlH 2008 г в часов на заседании диссертационного совета Д 004 026 01 при Горном институте УрО РАН по адресу 614007, г Пермь, ул Сибирская, 78а

С диссертацией можно ознакомиться в библиотеке Горного института УрО

РАН

Автореферат разослан «

2S"» ь/щсм>л 2008 г

Ученый секретарь диссертационного совета,

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

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

- провести анализ результатов инструментальных наблюдений за сдвижением земной поверхности на месторождениях нефти и газа,

- выполнить исследования упругих, прочностных и компрессионных свойств продуктивных пород на месторождениях углеводородов,

- провести аналитические исследования методов расчета уплотнения коллекторов при снижении исходного пластового давления,

- обосновать выбор наиболее представительных механических моделей для расчета напряженно-деформированного состояния горных пород на месторождениях нефти и газа и рассмотреть особенности их применения,

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

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

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

Научные положения, выносимые на защиту

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

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

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

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

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

Достоверность научных положений, выводов и рекомендаций достигается представительным объемом лабораторных и натурных измерений, применением

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

Научная новизна работы

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

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

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

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

- впервые для отдельных месторождений Западной Сибири, территории ВКМКС, УНГКМ и АГКМ по результатам компрессионных испытаний получены параметры «шатровой» модели поведения коллектора под нагрузкой, которые могут быть использованы для расчетов НДС продуктивных объектов данных месторождений,

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

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

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

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

Реализация работы Установленные на основе прогнозных расчетов параметры напряженно-деформированного состояния горных пород и земной поверхности использовались для обоснования мер охраны подрабатываемых объектов и создания геодинамических полигонов, которые были внедрены на ряде нефтяных месторождений Западной Сибири, севера Пермского края, Уренгойском и Астраханском газо-

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

Апробация работы Основные положения диссертации и результаты исследований докладывались и обсуждались на следующих совещаниях, конференциях и конгрессах международной конференции «Проблемы геодинамической безопасности» (Санкт-Петербург, 1997), XI Российской конференции по механике горных пород (Санкт-Петербург, 1997), международной конференции «Геодинамика и напряженное состояние недр Земли» (Новосибирск, 1999), международной научно-практической конференции «Геоэкология и современная геодинамика нефтегазоносных регионов» (Москва, 2000), международной конференции «Проблемы добычи и переработки нефти и газа в перспективе международного сотрудничества ученых Каспийского региона» (Астрахань, 2000), 1П Международном рабочем совещании «Геодинамическая и экологическая безопасность при освоении месторождений газа, его транспортировке и хранении» (Санкт-Петербург, 2001), XIII международном конгрессе по маркшейдерскому делу (Будапешт, 2007), на заседаниях ученого совета ПермГТУ, на технических советах ООО «Лукойл-Пермь», ООО «Лукойл-Нижневолжскнефть», ОАО «Сургутнефтегаз», ООО «Архангельскгеолдобыча», ООО «Юганскнефтегаз», ООО «Уренгойгазпром», ООО «Астраханьгазпром»

Публикации. Основные результаты исследований опубликованы в 27 работах, включая 1 монографию, в том числе 17 - в ведущих рецензируемых журналах, включенных в перечень ВАК

Структура и объем работы Диссертация состоит из введения, шести глав, заключения и изложена на 315 страницах машинописного текста, включая 108 рисунков, 25 таблиц и библиографический список из 177 наименований

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

1. Состояние изученности вопроса и задачи исследований

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

Обобщение достаточно обширного опыта геодезического мониторинга деформационных процессов позволило выявить основные факторы, которые определяют

возможность интенсивных и обширных просадок земной поверхности над длительно разрабатываемыми месторождениями (Сидоров В А, Кузьмин Ю О.) Также с помощью инструментальных наблюдений выявлены типы современных геодинамических движений земной поверхности и механизмы их возникновения. Созданы модели расчета аномальных движений, установлено их соответствие с региональными типами напряженного состояния земной коры (Кузьмин Ю О ) Большие объемы исследований выполнены в части разработки методов выделения активных разломов и геодинамически потенциально опасных зон, критерии оценки зон риска по различным параметрам, способы организации геодинамического мониторинга (работы институтов ВНИМИ, ЗСФ ИНГТ СО РАН).

Одними из наиболее опасных проявлений геодинамической активности на месторождениях нефти и газа следует считать техногенные сейсмические события Для оценки интенсивности этих явлений разработаны математические модели, которые сводятся к расчету напряженно-деформированного состояния горного массива, содержащего поверхности ослабления (Дж Райе, Барях А А , Roest J Р А , W Kuilman) При этом показатели динамических смещений по разломным структурам обычно определяют на основе упругой или упругопластических моделей Для более точных расчетов требуется использование специальных моделей, учитывающих полные диаграммы сдвига пород по плоскостям ослаблений

Для расчета оседаний земной поверхности на месторождениях нефти и газа используются различные методы. Часть из этих методов аналогична применяемым на угольных и рудных месторождениях и использует функции единичного влияния элементарного вынутого объема на земную поверхность (Ю П Борисов, А С. Маз-ницкий, JIM Середницкий) Также известен метод расчета оседаний, разработанный Geertsma и основанный на линейной теории упругости изотропной среды. Указанные методы позволяют осуществлять достаточно надежный прогноз максимальных оседаний земной поверхности. Однако если в массиве горных пород имеются различные структурные неоднородности, то применение данных методов становится проблематичным

В настоящее время при решении задач механики горных пород наиболее эффективными стали численные методы Большой вклад в развитие численных методов применительно к задачам расчета НДС горных массивов внесли Л Мюллер, Г Кратч, О Зенкевич, Г.Н Панде, В Виттке, А С Ягунов, С В Кузнецов, M В Курленя, А Б Фадеев, В Г Зотеев, О В Зотеев, А А Барях, Ю А Кашников, В M Серяков

В последнее время практически все работы, связанные с прогнозом оседаний на нефтегазовых месторождениях, также выполняются численными методами (Pande G N, Boade R R, Chin L Y, Siemers W T, Nagel N В , Grant U S , Plischke В, Sulak R. M, Thomas L К, Э В Калинин, H Б Артамонова и др ) При выполнении расчетов исследователи стремятся максимально подробно учесть особенности конкретных месторождений для повышения надежности прогнозов Однако такие работы выполнены для ограниченного числа месторождений и опыт подобных расчетов не обобщен

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

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

2. Анализ инструментальных наблюдений за сдвижением земной поверхности на месторождениях нефти и газа

Инструментальные наблюдения за сдвижением земной поверхности являются наиболее надежным методом контроля за состоянием подрабатываемых объектов, а также необходимы для обоснования и калибровки расчетных моделей, применяемых для прогноза напряженно-деформированного состояния горного массива Инструментальные наблюдения были организованы при участии автора на всех основных объектах исследований диссертационной работы - Уренгойском (УНГКМ) и Астраханском (АГКМ) газоконденсатных месторождениях, на нефтяных месторождениях севера Пермского края (Уньвинское, Чашкинское, Юрчукское, Сибирское, Шерш-невское), на ряде месторождений Западной Сибири (Усть-Балыкское, Мамонтов-ское, Приобское, Западно-Сургутское, Чумпасское, Ватинское, Варьеганское и др) Наряду с традиционными геодезическими наблюдениями по профильным линиям реперов широко применялись GPS-технологии, т к они дают возможность получать полные вектора сдвижений и осуществлять с довольно высокой точностью геодезический мониторинг обширных территорий месторождений углеводородов В ходе работ по проектированию наблюдательных станций были сформулированы основные принципы организации инструментальных наблюдений на месторождениях нефти и газа, которые вошли в нормативный документ - «Инструкцию по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе Верхнекамского месторождения калийно-магниевых солей».

На рис 1 представлен характерный пример наблюдательной станции, созданной на Сибирском нефтяном месторождении на севере Пермского края Станция состоит из нивелирной сети, включающей в себя 5 ходов общей протяженностью 27,5 км, а также GPS полигона Суммарное количество реперов нивелирной сети составляет 90 шт GPS - полигон представляет собой совокупность пунктов, расположенных как на территории месторождения, так и за его границами. Все пункты жестко связаны между собой системой векторов, образуя единую сеть, которая обладает максимальным количеством векторов на каждом пункте и большой избыточностью измерений, что в конечном итоге приводит к более высокой точности и надежности координатных определений Количество пунктов GPS-сети 13 шт, количество измеряемых векторов в сети - 30, площадь полигона составляет 73,3 км2

I 1

Рис. 1.

Схема наблюдательной станции на Сибирском месторождении

Условные обозначения контур нефтеносности ОР5-вектор СгРв-репер

репер нявелнроэяния Я класс* М 1: 100 000

Анализ результатов наблюдений показывает, что вертикальные смещения большинства реперов на Сибирском, Шершневском, Логовском, Уньвинском и Чашкинском месторождениях не превышают 20-30 мм (рис. 2). В картине распределения вертикальных сдвижений не выявлено четких особенностей, однако наблюдаются преимущественные оседания, что позволяет говорить о незначительных проявлениях процессов сдвижения, вызванных добычей нефти. Наблюдаются колебания отметок реперов, которые, скорее всего, вызваны неустойчивостью реперов в болотистом грунте, структурными особенностями поверхности, геодинамическими и локальными техногенными процессами. Максимальные оседания реперов составляют 15-35 мм. Результаты обработки ОР8-наблюдений также показывают, что горизонтальные и вертикальные величины смещений реперов часто не превышают предельных ошибок их определения. В целом можно утверждать, что величины смещений реперов весьма невелики, но в то же время отмечается устойчивая тен- |

денция незначительных оседаний земной поверхности над отрабатываемым место- 1

рождением.

наблюдательной станции Сибирского месторождения

Аналогичные результаты были получены на месторождениях Западной Сибири, где наблюдательные станции были созданы в районах крупных населенных пунктов Так, на Западно-Сургутском месторождении (район г Сургут) максимальные оседания земной поверхности не превышают 10 мм, на Чумпасском месторождении в районе гЛангепас за пять лет наблюдений (2001-2006 гг) проявилась устойчивая тенденция оседаний земной поверхности, хотя величины оседаний незначительны -4-5 мм за период наблюдений Сходные величины смещений зафиксированы на Варьеганском месторождении на территории г Радужный, на Усть-Балыкском и Мамонтовском месторождениях (г Нефтеюганск и г Пыть-Ях) и ряде других На наблюдательной станции Уренгойского газоконденсатного месторождения (район г Новый Уренгой) за период с 2003 по 2004 год были выявлены оседания реперов в среднем на 10 мм, при максимальных оседаниях 18 мм. На Астраханском газокон-денсатном месторождении за три года наблюдений (2003-2007г) четко сформировалась мульда оседаний величиной до 24 мм (рис 3)

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

что на данном этапе нефтяные месторождения разрабатываются с заводнением, т е без интенсивного снижения пластового давления

Рис 3 Графики оседаний реперов по профильной линии гл рп 5-гл рп 9 наблюдательной станции Астраханского ГКМ за период 2003-2007 гг

Однако сам факт наличия техногенных смещений земной поверхности говорит о необходимости прогноза и контроля напряженно-деформированного состояния горного массива при добыче нефти и газа вблизи ответственных объектов

3. Применяемые механические модели горных пород и их параметрическое обеспечение

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

Основные уравнения модифицированной шатровой модели (МССМ - модель) формулируются при рассмотрении стандартных компрессионных испытаний, те дренированного нагружения образца породы в стабилометре эффективными напряжениями О) > о2 = о3. Вводится эффективное гидростатическое напряжение с, де-виаторное напряжение я и коэффициент пористости е, как отношение объема пор к объему твердого тела сг = (1/3) (а,+2а3), я = а, - ст3, е=п/(1-п), где п- пористость

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

е=ТЧ-Я.1пст, е = ек-к1пст, (1)

где Я,, к - углы наклона прямых соответственно при нагрузке и разгрузке, Ы, ек- начальные значения коэффициента пористости При этом деформирование образца при разгрузке и повторной нагрузке считается упругим Согласно основных положений шатровой модели изотропная компрессия образца под давлением рс образует зону упругости ОАрс (рис 4) Объемные пластические деформации сжатия будут возникать при выходе напряжений за границу поверхности текучести Арс, которая имеет вид эллипса со смещенным относительно начала координат центром Согласно ассоциированного закона пластического течения поверхность текучести одновременно является также поверхностью пластического потенциала, т е

Р = (} = -£—+ о-р.=0, (2)

М а

где Б, (2 обозначают соответственно критерий разрушения и пластический потенциал, М - параметр линии критического состояния (СБЬ) вида q = М а

Рис 4 Виды поверхностей течения в шатровой модели

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

(3)

Ре=Р,о ехр, ^

В области низких нормальных напряжений (при а < рс/2) упрочнение материала невозможно и появление пластических деформаций связано с разрушением материала при сдвиге или при растяжении В диаграмме а-ц форма критерия разрушения от сдвига аналогична критерию Кулона-Мора

F = q-a 1ап ф * -С*, (4)

где tgф* = :^^-,C* = pc (М-1ёср*)/2 3 -БШф

При растяжении критерий разрушения записывают в виде

Р = -аз-0р=О,

где ор- прочность на растяжение

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

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

Для оценки интенсивности техногенных сейсмических событий использовалась модель деформирования горных пород по системам трещин Модель использует полные диаграммы деформирования скальных контактов, полученные В Лейхнитцем и П Ербаном на приборах прямого среза (рис 5)

Пиковая прочность

Допредельное х, деформирование

Рис 5 Полная диаграмма деформирования по контакту скальных пород

Согласно положений известной однородной модели В Виттке составляется уравнение вязкопластичности для необратимых относительных перемещений (б'"1'} берегов трещины

50х <Эстт

(5)

где{б4}=

■вектор скорости нормальных и касательных вязкопластических

смещений по трещине, {от} = {сгп,тги}т - нормальные и касательные напряжения в плоскости трещины, Рт, <3Т - критерий разрушения и пластический потенциал

Считая, что предельное сопротивление сдвигу тр выражается критерием Джаге-ра и с учетом полных диаграмм сдвига (рис 5) были получены критерии разрушения на стадиях упрочнения и разупрочнения В стадии упрочнения (58 < 5Р)

Рц „ШФ + с[1 - ехр (- Ьс п)]}

25, N

5 5

Р ^ р /

(6)

На стадии разупрочнения (58 >5Р)

Ртд + (1ёФ-*8<Р') ехр(-х(б5 -6,))]- {7)

с [1-ехр(-Ь сг„)] ехр(-х(5,-5р))

Соответствующие значения частных производных пластического потенциала имеют вид

00,

ад

= о ехр

25, 3.

ао,

5т „

-= 1,

13- -

■1&о ехР

ехР[-х(5,-5р)}

(8)

(9)

В этих формулах

Ф, с - угол внутреннего трения и сцепление; ф , 10 - остаточный угол внутреннего и угол дилатансии, Ь = ^(фт +10) - фт - угол внутреннего трения по трещине, бр - касательное смещение по трещине и его предельное значение, х - параметр разупрочнения

Рассмотренная модель неоднократно применялась для прогноза НДС подрабатываемых скальных массивов рудных месторождений, где показала свою эффективность Параметрами, которые в данной модели определяют характер полной диаграммы сдвига, являются предельное смещение по трещине 5Р и параметр разупрочнения х Для определения этих характеристик требуется очень сложное и точное оборудование В работе использовались табличные значения этих величин, полученные Лейхтницем, Ербаном и другими исследователями С целью снижения степени неопределенности, вызванной использованием табличных значений 8Р и %, при выполнении расчетов производился анализ влияния данных параметров на расчетные показатели напряженно-деформированного состояния массива

При обосновании параметров расчетных моделей использовались непосредственные испытания образцов горных пород, различные эмпирические зависимости и все доступные справочные данные Были получены зависимости упругих свойств пород от глубины залегания и всестороннего давления Так, ряд образцов известняка Астраханского газоконденсатного месторождения (АГКМ) был исследован в Институте проблем механики РАН на специальной установке - испытательной системе трехосного неравнокомпонентного сжатия, которая позволяет создавать в образцах однородные напряженные состояния с любым соотношением напряжений сть 02103 Испытания показали, что наблюдается выраженный рост модуля упругости в 1,4-1,7 раза при росте эффективного давления до бОМПа (табл 1)

Таблица 1 - Зависимость модуля упругости известняка от всестороннего давления

со, МПа 15 30 45 60

Б, ГПа 19,2 23,7 27,1 29,1

Полученные зависимости в дальнейшем учитывались при математическом моделировании НДС коллекторов

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

Р = п Р„ + рта, (10)

где п - открытая пористость

Шгрцлевское м-ние скв 79 обр 19 0175 1

0170 ■ 0165 ' 0160 0155

Чумпасское м-ние обр 332-221-1 Эп 0 252 -

0 12 3 4

1_п(Р,МГ!а)

1_п(Р МПа)

Родниковое м-ние обр 5961 96

0 202 02 ' 0198 ^ 0193 0194 ' 0192 019

018а

2 3

1_л(Р МПа)

Уренгойское м-мие ппБу8,обр 106

0 18 I 0175 : 017 0 165 0 16

ч^

1_п(Р мпа)

Рис 6 Примеры компрессионных кривых ряда месторождений

При этом для песчано-глинистых и карбонатных коллекторов коэффициент сжимаемости твердой фазы на один-два порядка ниже, чем коэффициент сжимаемости пор, т е объемная деформация пористых пород возникает главным образом за счет деформаций порового пространства Для изучения сжимаемости порового пространства коллекторов производились компрессионные испытания, т е получали зависимость пористости (или коэффициента пористости) образцов от всестороннего эффективного давления при нагрузке и разгрузке (рис 6)

Компрессионная зависимость дает полное представление о деформируемости порового пространства, т к с ее помощью можно найти коэффициент сжимаемости пор в заданном интервале всестороннего давления

Р„ = Лп/(п Аа) = Де/(е Да), ^^

где п - пористость, е - коэффициент пористости

Поскольку сжимаемость твердой фазы значительно меньше сжимаемости пор, то общую сжимаемость породы можно без большой погрешности определить по формуле (10) с помощью табличных значений ртв Испытания показали, что в широком интервале напряжений пористость можно представить в виде линейной функции от логарифма всестороннего давления Это дает основание применять для расчета деформаций коллекторов хорошо разработанные «шатровые» модели

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

4. Прогноз напряженно-деформированного состояния горных пород при разработке месторождений углеводородов

Породы, залегающие на больших глубинах, практически всегда насыщены одним или несколькими флюидами (обычно водой, нефтью или газом) Изменение

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

je,+2Ge,+a р, <=(K-f}v+2Gsy+ap, (i2)

<=^K-^jsv+2GSl+a р;

где ст',, - полные напряжения, К, G - модуль объемного сжатия и модуль сдвига, р -поровое давление, a - коэффициент разгрузки

Для полной характеристики упругой пористой среды уравнения состояния (12) должны быть дополнены уравнением фильтрации находящегося в порах флюида

р, (.з)

r\ 8t я р

где г) - модуль Био, к, ц - проницаемость породы и вязкость флюида, ev - объемная деформация породы.

В наиболее общем варианте для расчета напряженно-деформированного состояния пористой среды необходимо совместное решение уравнений теории упругости (12) и уравнения фильтрации (13) Однако в ряде случаев, при определенных видах напряженного состояния, можно упростить вид уравнения фильтрации и решать его независимо от уравнений теории упругости Например, если деформирование пород происходит с сохранением постоянного объема (еу=0), то уравнение фильтрации преобразуется к виду

& =-*-V'p,

at n-M(pf+p,)

где pf - сжимаемость флюида, Pi - сжимаемость породы при ev=0

Другие характерные виды напряженного состояния возникают при одномерном уплотнении (£х=е,=0, ст'г = const), при плоской деформации (ez= 0, ст'х = а'у= const) и при постоянстве среднего давления (а'х+ ст', + ст'г)/3 = const В этих случаях в уравнении фильтрации будут фигурировать коэффициенты сжимаемости породы рг, рз и р4 Те, уравнение фильтрации в общем виде записывается

dp = k у>р (14)

a nn(pf+pr)

где р, - некоторый постоянный коэффициент сжимаемости породы При этом можно показать, что при произвольном напряженном состоянии выполняется соотношение р! < рг < (З4 Уравнение (14) решается независимо от уравнений теории упругости, т к не содержит в явном виде еу Коэффициент сжимаемости флюида обычно значительно больше сжимаемости породы, поэтому общая сжимаемость системы флюид-порода робпгКРгФг) определяется в основном коэффициентом & В качестве примера в табл 2 представлены коэффициенты сжимаемости, рассчитанные для условий, характерных для нефтяных месторождений севера Пермского края

Таблица 2 - Коэффициенты сжимаемости системы флюид-порода

Сжимаемость породы рг, МПа"1 Pi Р2 Рз Р4

3,25 105 1,214 Ю-4 1,658 10"4 1,925 10"4

Общая сжимаемость ро6ш, МПа1 1,032 10"' 1,121 10'3 1,166 10"3 1,192 10"3

Можно видеть, при изменении коэффициента сжимаемости породы от минимального значения (3] до максимального Pi общая сжимаемость системы Р„6щ изменяется довольно незначительно - на 15%, а для наиболее вероятных значений р (от Р2 до рз) изменяется еще меньше - на 5% Подобное соотношение между коэффициентами сжимаемости характерно для большинства практических случаев, поэтому гидродинамическое моделирование разработки нефтяных и газовых месторождений производится на основе уравнений вида (14), т е деформационные свойства горных пород учитываются применением соответствующих коэффициентов сжимаемости

Рг

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

Наиболее простой вид напряженного состояния при добыче нефти и газа возникает при одномерном уплотнении, т е при отсутствии горизонтальных деформаций (8х=еу =0) и постоянстве полных вертикальных напряжений (o'z = const) Для данных условий из уравнений закона Гука (12) можно найти вертикальную деформацию уплотнения

= —~7г"а(Р° _ Р)= см а ДР> (15)

где см - коэффициент одномерного уплотнения, dp - падение давления

Уравнение (15) можно применять для глубоко залегающих плотных низкопористых пород, деформации которых при больших давлениях с достаточной степенью точности соответствуют теории упругости В более общем случае необходимо учитывать специфические законы деформации пористых сред Работы К Терцаги,

М. Био, Ф Гасмана, В М. Добрынина показывают, что объемные деформации пористой среды складываются из деформаций твердой фазы и порового пространства при изменении внешнего всестороннего давления су1 и давления жидкости в порах

y = pd(a,-p) + MP <16)

Данное уравнение преобразовывается к наиболее простому виду при использовании принципа эффективных напряжений dV

= P(dc' - ccdp) = Pd(0' - ар) = pdc (17)

Слабой стороной методики использования эффективных напряжений является неопределенность, связанная с учетом изменения коэффициента разгрузки а при росте всестороннего давления Однако в большинстве практических случаев зависимость коэффициента разгрузки от давления можно не учитывать, поэтому в расчетах использовались эффективные напряжения Сжимаемость твердой фазы коллекторов принималась по табличным данным (для песчано-глинистых коллекторов = 0,030 ГПа"1, для карбонатных - Р™ = 0,025 ГПа"1), сжимаемость пор определялась по результатам компрессионных испытаний на основе зависимостей (1)

р lde=1=nx

е áa п и

С учетом этих соотношений можно найти объемную деформацию коллектора при изменении эффективного давления от Oí до о2

sv=(l~n)x|—+ р„ JdCT (19)

а, О J,

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

Ah = h (l-n)Xln^i (20)

где Ah - уплотнение, h - начальная толщина коллектора, c^i, - начальное и конечное значение эффективных вертикальных напряжений

Специальные сопоставительные расчеты показали, что использование табличных значений для сжимаемости твердой фазы коллекторов обеспечивает достаточную для практических целей точность Даже при неучете сжимаемости твердой фазы погрешность расчета вертикальной деформации обычно не превышает 10% Это подтверждает правомерность использования механических моделей «шатрового» типа Поскольку ветчина одномерного уплотнения представляет собой верхний предел возможных оседаний поверхности, то полученные аналитические зависимости позволяют сделать предварительные выводы об общем уровне напряженного состояния массива Это дает возможность в дальнейшем более обоснованно подходить к построению общих расчетных моделей изучаемых объектов

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

{ЕМ8е|}+КР} (21)

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

О при Б :< 0]

при Р>0|, ^

где Г| - параметр вязкости, играющий роль регулятора сходимости итерационного процесса. Непосредственно деформации находятся путем численного интегрирования (22) по времени

(23)

Реализация моделей была выполнена в конечно-элементных программах "ОЕОТЕСН" (разработка ПермГТУ) и "А№У8"

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

{а'} = {а} + {т}ар, (24)

где {а1} = {ст'х, а 1У, т'ху, т'у2, т'и}Т- полные напряжения, {о} = {о*, оу, ст2, тху, ТуИ 1Н}т - эффективные напряжения, {га} = {1,1,1,0,0,0}т

С учетом (24) можно получить общее соотношение МКЕ в виде

[Г] {5е} = Г} - {!*}, (25)

где {Рер} - вектор дополнительных узловых сил, обусловленных давлением насыщающей жидкости

С помощью выражения (25) был решен ряд модельных задач для исследования характера и степени влияния различных факторов на параметры уплотнения коллекторов и деформации горного массива Рассматривались как общие моменты моделирования напряженного состояния горных массивов (выбор типа и размеров модели, задание граничных условий), так и специфические особенности, характерные для месторождений углеводородов Было установлено, что ориентировочные размеры расчетной области по вертикали и горизонтали должны быть в 4-6 раз больше соответственно глубины и радиуса коллектора При данных размерах влияние типа задаваемых граничных условий становится несущественным Рассмотрение объемной задачи с коллектором эллиптической формы показало, при соотношении большой и малой полуосей а/Ь > 2 профиль мульды сдвижения вдоль малой полуоси эллипса можно определять расчетом на плоской модели (плоская деформация)

Для анализа полного характера уплотнения коллекторов была выполнена серия расчетов осесимметричной модели, в которых варьировалось отношение G/Gk (Gt -модуль сдвига пород коллектора, G - модуль сдвига вмещающих пород) и отношение R/H (R - радиус коллектора, Н - глубина залегания) По результатам расчетов строились графики относительного уплотнения (ДЬ/ДЬшах) в зависимости от текущего радиуса г (рис 7) для анализа хараетера уплотнения на всем протяжении продуктивных пластов (максимальное значение Ahmix соответствует условиям одномерного уплотнения)

СМ* —♦—0 2 —■— 0 5 —Л— 1 —♦— 1 5 —*—2 1Ш —•— 1 —■—2 —S

r/R

Рис 7 Влияние упругих свойств вмещающих пород и геометрических размеров коллекторов на характер их уплотнения

Было установлено, что уплотнение коллекторов в целом выше при более слабых вмещающих породах, особенно на границах области снижения давления В то же время в центральной части уплотнение близко к максимальному и слабо зависит от жесткости вмещающих пород Уплотнение также увеличивается с ростом радиуса коллектора, однако уже при R/H > 2 влияние этого фактора становится несущественным Наконец, относительное уплотнение коллектора зависит от его начальной толщины При прочих равных условиях относительное уплотнение выше для более тонкого коллектора, что проявляется на флангах залежи, в центральной части влияние этого фактора практически отсутствует В целом результаты данных расчетов показали, что даже при небольшом отношении R/H деформации коллекторов в центральной части нефтегазовых месторождений близки к условиям одномерного уплотнения Однако в краевых частях области снижения исходного пластового давления условия одномерного уплотнения не выполняются и напряженное состояние имеет более сложный вид Оседания земной поверхности, кроме вышеназванных параметров, зависят также от соотношения жесткости покрывающих и подстилающих пород Наибольшие оседания земной поверхности имеют место при более жестких подстилающих породах В целом степень влияния упругих свойств вмещающих пород на уплотнение коллекторов и оседания земной поверхности уменьшается с ростом размеров коллектора, т е деформации горного массива определяются прежде всего величиной уплотнения продуктивного слоя

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

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

скважины имеет осесимметричный характер и может быть найдено из выражения

<26)

где г0, рп и -безразмерные величины = гп = г/Яс, 1ц = к 1 /

(п Ц р0бщ Я2С), р0=2якЬ [ро-р] / (С! ц)

В этих формулах г, I - текущий радиус точки и текущее время работы скважины, IV радиус контура питания, радиус скважины, к - проницаемость пласта, Ь - толщина пласта, ц - вязкость флюида, п - пористость, (30бЩ - коэффициент сжимаемости системы флюид-порода, ро, р - начальное и текущее пластовое давление

На рис 8 показана этора пластового давления для одного из расчетных примеров и соответствующие ей деформации уплотнения коллектора Расчетная эпюра уплотнения соответствует картине падения пластового давления, однако величины уплотнения в центре модели и на контуре питания составляют 61 0 и 608мм, те различаются совершенно незначительно При этом, если использовать средневзвешенное значение давления, то уплотнение составит 61 6мм Данный пример показывает, что в большинстве случаев для расчета деформаций горного массива можно не учитывать неравномерность пластового давления, обусловленного работой отдельных скважин Этот вывод имеет важное практическое значение, т к при использовании показателей средневзвешенного давления нет необходимости существенного сгущения конечно-элементных сеток в районе добывающих скважин, что упрощает построение расчетных моделей

|—пласт давление -—уплотнение]

Рис 8 Пластовое давление и уплотнение коллектора в зоне влияния скважины

Следующий рассматриваемый фактор связан с тем, что коллектора нефти и газа имеют куполообразное строение, т к. приурочены к сводовым частям геологических структур-ловушек углеводородов. Поскольку углы наклона крыльев складок составляют первые градусы или даже минуты, расчетные модели обычно строят в виде прямолинейных пластов с постоянной глубиной залегания. Влияние подобной идеализации геометрии залежи было рассмотрено на численном примере коллектора толщиной 50 м, радиусом 11=1500 м, залегающего на глубине 11=1500 м (рис 9) Предполагалось, что водонефтяной контакт расположен на 100м ниже купольной части, что при данных размерах залежи является довольно большой величиной Рассчитывались оседания земной поверхности, которые сравнивались с аналогичными величинами для коллектора плоской формы с постоянной глубиной залегания Н=1500м

•V Н=1500м

R=1500м

Рис 9 Фрагмент конечно-элементной модели коллектора куполообразной формы

В базовом варианте (Е=5000 МПа, v=0,2, падение давления 5 МПа) максимальное оседание поверхности для плоского коллектора составило 28 9 мм, для куполообразного - 27 8мм, т е меньше на 4% Данное явление можно объяснить тем, что сводчатая структура вмещающих пород является более жесткой С увеличением размеров коллектора разница в оседаниях уменьшается Расчеты показали, что оседания поверхности для двух вариантов геометрии коллектора практически не различаются уже при радиусе R=3000 м, т е при R/H=2 Т о, идеализация геометрии коллекторов в виде плоских пластов с постоянной глубиной залегания наиболее заметна для месторождений с малым отношением R/H При этом допущение о плоской форме коллектора дает незначительное увеличение сдвижений массива, т е обеспечивает расчетам некоторый запас прочности

Явление распространения воронки депрессии за пределы водонефтяного (или газоводяного) контакта связано с тем, что после пуска скважин в эксплуатацию за счет перепада давлений начинается приток законтурных вод в залежь Общие особенности развития деформационных процессов при падении давления за пределами залежи были рассмотрены на численном примере коллектора радиусом R=3000m, глубиной Н=1500 м и толщиной ЮОм При Е=5000 МПа, v=0,2 и падении давления 5 МПа максимальное уплотнение коллектора составляет 90 мм Рассчитывались деформации земной поверхности при распространении депрессионной воронки за пределы залежи на некоторое расстояние AR, которое варьировалось от 0 до 1 0R Если в базовом варианте (депрессионная воронка не выходит за пределы залежи) граничный угол сдвижения, определяемый на точку с оседанием 1 Омм, составил 50°, то при AR=3000 м граничный угол достиг 27° Данный пример служит объяснением того факта, что на ряде месторождений инструментальные наблюдения фиксируют

аномально низкие углы сдвижения Так, на Уренгойском газоконденсатном местро-ждении установлено, что мульда оседания земной поверхности распространяется за границы ГВК до 8-10км, те граничный угол сдвижения доходит до 10-15° Это можно объяснить только падением давления в областях, удаленных от газоводяного контакта Указанные особенности необходимо учитывать при проектировании инструментальных наблюдений за сдвижением земной поверхности. При определении длин профильных линий целесообразнее всего руководствоваться накопленными сведениями о граничных углах сдвижения на уже исследованных месторождениях со сходными горно-геологическими условиями При отсутствии таких данных минимальное значение граничного угла сдвижения необходимо определять расчетными методами с максимальным учетом особенностей конкретного месторождения

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

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

р = п р (27)

Н " (1 + е) с

Используя соотношения теории упругости, можно установить зависимость между индексом декомпрессии к и упругими параметрами Е и V

_ 1 3(1 - 2у) к

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

0 ~ О 1 '

2По+3 £29)

где т,0ф/^4М^ЗЛ),Л = 1-кА

В этм случае можно формально ввести некоторый 'пластический' коэффициент Пуассона по соотношению ур1 = ко/(1+ко) и установить зависимость между компрес-

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

Особенности напряженно-деформированного состояния горного массива при уплотнении коллекторов согласно МССМ-модели были рассмотрены на ряде примеров, условия которых обеспечивали максимальную величину уплотнения от 90 мм до 2 7 м В одном из примеров рассматривался коллектор толщиной 50 м, радиусом 6000м, залегающий на глубине 2000 м Пористость коллектора 20%, Я= 0,010, к = 0,005, v=0,20, М=0,984, объемный вес покрывающих пород 0,022 МН/м3 Рассчитывались деформации коллектора и горного массива при снижении пластового давления от 20 до 14 МПа Данные условия обеспечивают максимальное уплотнение 90мм Расчетный характер напряженного состояния коллектора и горного массива показан на рис 10

т, МПа

а =- 211883

в =- 14126

с =- 070 63 6

d =- 127е-04

е = 070611

f = 141234

g = 211858

н = 2в24в1

i = 353105

q, МПа

N В

б >< А 1 'Ч 4 N \ \ Ч ч N \ \

V \ \ 1

10 15 20

ст, МПа

Рис 10 Характер напряженного состояния коллектора при падении пластового давления

Касательные напряжения в горном массиве приурочены к краевой части продуктивной зоны и весьма малы по величине, что легко объясняется незначительным уплотнением коллектора Если взять две точки в центре модели и на краю продуктивной зоны, то в них исходное напряженное состояние коллектора одинаковое и в диаграмме я - а характеризуется точкой А (рис 10). Параметр рс, задающий начальную границу эллиптической области упругих деформаций, равен 22,6 МПа После падения пластового давления на 6 МПа напряженное состояние коллектора в центре модели достигает точки В и формирует новую, более обширную область упругих деформаций с параметром рс= 28,2 МПа Точки А и В лежат на прямой с уравнением q=r^o сг, т е прирост эффективных горизонтальных напряжений соответствует теоретическому значению коэффициента бокового давления Напряженное состояние в краевой части коллектора характеризуется точкой С, которая лежит выше линии одномерного уплотнения, т е девиаторная часть тензора напряжений прирастает ин-тенсивне Тем не менее, точка С не достигает предельной линии я=Ма, те разупрочнения не возникает и коллектор в краевой зоне также деформируется в режиме уплотнения

Расчеты, выполненные для условий существенного уплотнения коллектора (2 7 м), показали сходные результаты Несмотря на значительные абсолютные величины сдвижений, касательные напряжения в горном массиве оказались довольно небольшими Данный факт можно объяснить тем, что сдвижения распределяются на значительные области горного массива Напряженное состояние в краевой части коллектора не достигает предельной линии я=М а, т е породы также деформируются в режиме уплотнения При этом рост эффективной нагрузки происходит по траектории, приблизительно нормальной к эллиптической поверхности текучести Данные расчеты показывают, что применение МССМ-модели для расчета деформаций коллекторов обеспечивает представительные результаты при большом разнообразии горно-геологических условий и физико-механических свойств продуктивных пород

5. Деформирование пород на контактах блоковых структур и оценка интенсивности техногенных сейсмических явлений

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

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

Подобные расчеты были выполнены для анализа напряженного состояния блочной структуры горного массива в регионе Верхнекамского месторождения ка-лийно-магниевых солей (ВКМКС) Моделируемый участок массива имеет размеры по глубине 5,0 км, в плане 45x45 км и захватывает гг Березники и Соликамск За основу расчетной модели были взяты тектонические блоки II,III,ГУ,У рангов, выделенные на территории ВКМКС проф В В Филатовым (Уральский государственный горный университет, г Екатеринбург) На данной модели были рассмотрен ряд ситуаций, способных вызвать аномальные движения земной поверхности, зафиксированные геодезическими наблюдениями в разломных зонах Анализировалась реакция блочного массива на природные сейсмические воздействия, вызванные незначительным природным землетрясением интенсивностью 3-4 балла Расчеты показали, что при действии волны сжатия на земной поверхности в районе нарушений возникает всплеск поднятий, масштабы которого зависят от упругих свойств блоков и разломных зон При действии волны растяжения соответственно возникают оседания поверхности В качественном плане полученные смещения соответствует Р- и у-аномалиям современных движений земной поверхности, которые были типизированы В А Сидоровым и Ю О Кузминым Данные результаты позволяют сделать вывод о том, что аномальные величины вертикальных движений земной поверхности, являющиеся частью геодинамической обстановки, могут быть обусловлены природной сейсмичностью региона

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

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

Задача была рассмотрена для условий, характерных для нефтяных месторождений севера Пермского края Рассчитывалось напряженно-деформированное состояние горного массива при падении пластового давления на 5 МПа в продуктивном пласте мощностью 30 м, находящегося на глубине 2000 м Характеристики продуктивных пород (Е = 2000 МПа, V = 0 25) в данных условиях обеспечивают максимальное уплотнение коллектора величиной ЛИ = 62 5 мм Предполагалось, что в краевой части коллектора по всей мощности разреза проходит субвертикальная зона междублокового контакта (рис 11)

Рис 11 Мульды оседания при наличии и отсутствии зоны ослабления мощностью 100 м (к = 0 2)

Поскольку ширина контактных зон является крайне неопределенным параметром, ее мощность в расчетах варьировалась от 100 до 500 м Модуль упругости пород в контактных зонах определялся в зависимости от модуля упругости пород соответствующего слоя по соотношению Е' = Ео х к, где Е', Е0 - модуль упругости соответственно нарушенных и ненарушенных пород, к - коэффициент ослабления Коэффициент ослабления к варьировался от 0,8 до 0,2 Расчеты показали, что зоне нарушения наблюдается некоторая концентрация деформаций, где на контакте нарушенных и ненарушенных пород происходит прирост оседаний относительно мульды сдвижения в монолитном массиве (рис 11) Величина прироста оседаний определяется мощностью и физико-механическими свойствами ослабленных пород Оседания возрастают при снижении упругих свойств пород и уменьшении ширины

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

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

При анализе геодинамической обстановки на месторождениях нефти и газа одной из наиболее важных задач является прогноз интенсивности техногенных сейсмических явлений Известные математические модели очага техногенного (а в отдельных случаях и природного) землетрясения говорят о нем как о модели неустойчивого роста трещины в разломе, т е активизации существующих тектонических разломных структур в форме сдвига их бортов Впереди развивающейся в разломе магистральной трещины возникает область объемного разрушения за счет создания множества микротрещин Трущиеся при сдвиге борта разлома порождают дробленый материал, в результате чего большая часть энергии, высвобождаемой при сдвиге, расходуется на образование в бортах разлома зоны дробленой породы Указанным представлениям соответствует известная модель скольжения с разупрочнением ДжРайса(рис 12)

Рис 12 Модель активизации тектонического разлома

Началу скольжения по разлому соответствует достижение касательным напряжением тгм своего максимального значения хр До этого наблюдается как рост касательного напряжения, так и рост касательного к поверхности разлома смещения 6S При достижении величиной 5S значения 5Р (и соответственно касательного напряжения тр) начинается неустойчивый рост трещины в разломе и падение касательных напряжений до величины остаточной прочности т* При достижении остаточной прочности перемещение 5S может расти неограниченно при неизменной величине т ,

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

Разность параметров (тр - т ) рассматривается как сброс напряжений Скорость высвобождения энергии е при подвигании сдвиговой трещины (энергия разрушения) дается выражением

Ди'

е= ЯиЛи]-т*]с1(Ди), ОС)

где Ди* - величина перемещения при скольжении

Рассматриваемая задача сводится к расчету НДС горного массива, содержащего поверхности ослабления Для реализации модели Дж Райса хорошо подходит модель деформирования пород по системам трещин, т к она использует полную диаграмму деформирования пород по контактам (рис 5), которая практически аналогична модели скольжения с разупрочнением (рис 12) Реализация модели для оценки магнитуд техногенных сейсмических событий заключается в следующем

Разлом представляется средой, разбитой системой трещин вдоль поверхности раздела, В результате численного решения методом конечных элементов задачи скольжения с разупрочнением в каждом элементе разлома, где касательные смещения 55 превысили величину 5Р, определяется сброс напряжений (тр - т*) Также определяется разница между касательным смещением 58 на участке стабилизации итерационного процесса и величиной 5Р Значение высвобождающейся в 1-м элементе энергии находится как

г,=(тР - т*)(5Р - 5/), (31)

где 5, - максимальное достигнутое касательное смещение в 1-м элементе

Единичное значение энергии для решаемой задачи находится как среднее значение по элементам, вышедшим в запредельное состояние Нормированное значение высвобождающейся энергии получается путем перемножения найденной величины на вертикальный размер Ь зоны запредельного состояния

е = Ь-£

(32)

где п - число элементов, вышедших в запредельное состояние

Полная энергия сейсмического события получается перемножением нормированного значения (32) на линейный размер зоны сдвига по простиранию разлома При этом допускается определенный произвол в выборе размеров сдвигающихся зон, однако расчеты говорят о том, что это не имеет существенного значения

Оценочные расчеты магнитуд сейсмических явлений были выполнены применительно к отработке Уньвинского нефтяного месторождения на территории Пермского края На первом этапе были рассчитаны оседания горного массива при падении давления на 6, 10 и 20 МПа Применяемые при этом параметры модели были определены из условия наилучшей сходимости расчетов с результатами инструментальных наблюдений На втором этапе более детально рассматривался расчетный фрагмент в районе коллекторов, вырезанный из общей расчетной схемы (рис 13)

600 м Пласт Бт

карбонаты зоны влияния нагнетательных скважин \ \ 1 тектонический разлом

терригенные породы ^Е

- 1500 м | Пласт Тф

-V

Рис 13 Общая расчетная схема задачи активизации разлома

В качестве граничных условий на границах фрагмента задаются перемещения, полученные из решения общей задачи Оценивалась возможность возникновения динамических подвижек по нарушению в виде ослабленной вертикальной зоны мощностью 5 м, рассекающей продуктивную толщу и слой карбонатов Для прочностных свойств разломной зоны задавалось С = 0,05 МПа, ср = 20°, ф" = 10°, 10 = 10°, 5Р, = 1,0 мм

Расчеты показали, что без противодавления жидкости в разломе сдвига его бортов не возникает даже при падении давления в коллекторах на 20 МПа Возможность динамических сдвижений возникает при операциях нагнетания флюида для поддержания пластового давления Была выполнена серия расчетов с вариацией давлений нагнетания жидкости в разлом и падения давления в коллекторах Магни-туды возникающих сейсмических событий оценивались по формуле

М=(1о^-4 9)/1 5, (33)

где Е -нормированное значение энергии, уменьшенное в 100 раз согласно значения сейсмического кпд г|=0,01. При этом условно рассматривались ситуации с размерами разлома по простиранию Ь=250 м, 500 м и 1000 м Было установлено, что с увеличением размера Ь свыше 1000 м его дальнейший рост уже незначительно сказывается на величине магнитуды (рис 14)

Рис 14 Зависимость магнитуд сейсмических событий от размера разлома по простиранию и параметра 5Р

Также на рис 14 представлен график зависимости магнитуд от основной характеристики полной диаграммы сдвига - параметра 5Р На рис 14 видно, что при уменьшении бр от 0,5мм до нуля магнитуда может вырасти не более чем на 0,2 единицы, т е довольно незначительно Т е в рассматриваемом примере для анализа техногенной сейсмичности вполне можно воспользоваться табличными значениями показателя 5Р В целом количество выделяемой энергии довольно слабо зависит от падения пластового давления (в рассматриваемых условиях) и определяется прежде всего давлением нагнетания флюида, размерами нарушения, а также характеристиками полной диаграммы сдвига по поверхности раздела

Полученные значения магнитуд не превышают 1,5, что следует считать предельно возможным для данных условий Рост магнитуды, например, до 2,5 не представляется реальным, тк для этого требуется увеличение в 10 раз единичного значения энергии, определяемой по формуле (32) Это возможно только в случае активизации весьма крупного разлома с размерами зоны заводнения по высоте более 300м Сейсмические события подобной интенсивности фиксируются только с помощью чувствительных приборов и не оказывают заметного влияния на поверхностные, подземные объекты и геологическую среду региона месторождения

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

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

6. Практическое применение результатов исследований

Прогноз оседаний земной поверхности на нефтяных месторождениях Западной Сибири Расчеты напряженно-деформированного состояния горного массива и земной поверхности были выполнены на целом ряде нефтяных месторождений Западной Сибири (Усть-Балыкское, Мамонтовское, Чумпасское, Западно-Сургутское и ряд других), где в зоне влияния добычи нефти находятся ответственные объекты -жилая застройка и сложное промышленное оборудование

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

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

Таблица 3 - Результаты прогноза оседаний для ряда нефтяных месторождений

Западной Сибири

Месторождение Продуктивные объекты Суммарная мощность, м Макс оседание, мм

Усть-Балыкское БС1-5, БС10 35 122

Мамонтовское АС4-6, БС8, БС10-11 33 124

Правдинское БС5-6, БС8-9 20 40

Приобское АС10-АС12 45 360

Западно-Сургутское ЕС 1-4, БС10 25 89

Восточно-Сургутское БС10, ЮС 1-1, ЮС2-2 23 45

Родниковое БС12 10 36

Чумпасское АВ1, БВ6, БВ18-22 31 124

Ватинское АВ1-2, БВ8, ЮВ1 24 82

Оседания получены для падения исходного пластового давления на 3-5 МПа, что в среднем характерно для месторождений Западной Сибири Расчеты говорят, что при отработке месторождений системами с поддержанием пластового давления максимальные оседания поверхности не превышают 100-150 мм, иногда при большой мощности продуктивных объектов - 300-350 мм Этот вывод подтверждается имеющимся опытом инструментальных наблюдений на Усть-Балыкском геодинамическом полигоне Оседания подобной величины, равномерно распределенные на огромной площади месторождений, не мохуг нарушить нормальный режим эксплуатации объектов на земной поверхности Однако потенциальную опасность могут представлять концентрации деформаций на земной поверхности, вызванные различного рода структурными неоднородностями горного массива

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

Численное моделирование процессов сдвижения на Уренгойском газоконден-сатном месторождении На территории крупнейшего в мире Уренгойского газокон-денсатного месторождения находится крупный населенный пункт - город Новый Уренгой Инструментальные наблюдения на станции, созданной на месторождении в 1974-76 гг, к 1995 году зафиксировали оседания величиной до 340 мм При продолжении добычи газа возможно негативное влияние деформационных процессов на подрабатываемые объекты, в связи с чем был выполнен прогноз напряженно-деформированного состояния горного массива.

Деформационные параметры массива определялись методом «обратных расчетов» Были созданы расчетные схемы по ряду профилей, где продуктивные объекты строились по картам эффективных газонасыщенных толщин, а величины действующих нагрузок определялись по картам изобар В результате серии расчетов было установлено, что удовлетворительное соответствие расчетных и экспериментальных данных достигается при параметрах шатровой модели для пород сеномана X = 0,007,

Полученные данные позволили выполнить надежный прогноз напряженно-деформированного состояния горного массива при дальнейшей добыче газа и оценить влияние отработки месторождения на жилую застройку города Новый Уренгой и другие ответственные объекты Расчеты показали, что при сохранении существующих темпов отбора газа максимальные оседания к 2010 году достигнут 650мм Городская черта находится в краевой части мульды сдвижения и оседания земной поверхности на данной территории будут составлять до 500 мм Вследствие довольно существенных величин прогнозных оседаний была разработана и реализована наблюдательная станция для контроля деформационных процессов в пределах городской черты

Прогноз напряженно-деформированного состояния горного массива на Астраханском газоконденсатном месторождении Астраханское газоконденсатное месторождение (АГКМ) в геомеханическом плане имеет ряд существенных особенностей Продуктивный пласт залегает на значительной глубине («4000 м) и представлен довольно крепкими известняками, залежь характеризуется аномально высоким пластовым давлением (63 МПа) Довольно значительная часть разреза покрывающей толщи представлена породами кунгурского яруса, которые залегают в виде соляных куполов сложной формы Сложная солянокупольная тектоника, а также несимметричный характер области снижения исходного пластового давления обусловили необходимость разработки объемной конечно-элементной модели (рис 16) Упругие свойства пород оценивались по данным акустического каротажа покрывающей толщи по соотношению

скоростей упругих волн, В качестве основного «базового» слоя выступали продуктивные карбонатные породы.

Рис. 16. Объемная конечно-элементная модель участка АГКМ

На первом этапе определялись оседания земной поверхности при состоянии пластового давления на 01.01.2000 г., которые показали близкое соответствие с данными замеров по имеющимся глубинным реперам. На следующем этапе выполнялся прогноз оседаний при отработке АГКМ на 2010 г. При предполагаемой карте изобар на 2010 г., составленной по данным ВНИИГАЗ, расчеты показали мульду сдвижения с участком плоского дна и максимальным оседанием »200-210 мм. Величины оседаний довольно малы, однако в зоне влияния отработки газа находится такой сложный и ответственный объект, как газоперерабатывающий завод.

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

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

Для геомеханической оценки возможности совместной добычи нефти и калия на первом этапе был выполнен прогноз напряженно-деформированного состояния

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

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

Наиболее подробно рассматривалось геомеханическое поведение аномальных особенностей строения водозащитной толщи (листрических разрывов, сдвигов, открытых секущих трещин), т к сохранность ВЗТ является основным условием безопасной отработки ВКМКС При этом задавалось наиболее неблагоприятное сочетание входящих параметров для создания определенного запаса надежности Было получено, что деформации горного массива при отработке нефти вследствие своей малости не сказываются на состоянии водозащитной толщи при наличии в ней аномальных особенностей вида листрического разрыва и открытых секущих трещин Так, оценка возможности роста открытой секущей трещины в ВЗТ2 по энергетическому критерию механики разрушения показала, что скорость высвобождения энергии О во всех случаях не достигает критического значения 1000 Па м

Т о, добыча нефти на территории ВКМКС не оказывает негативного влияния на конструктивные элементы системы разработки калийных рудников и не активизирует структурно-тектонические особенности массива, расположенные в пределах ВЗТ и продуктивной соляной и калийно-магниевой толщах Данные результаты говорят о возможности совместной добычи нефти и калия на территории ВКМКС Не существует геомеханических и геодинамических факторов, которые в принципе не позволяли бы осуществлять разработку нефтяных залежей, естественно, при соблюдении определенных условий К таким условиям следует отнести применение соответствующей конструкции нефтяных скважин, которая обеспечивает надежную изоляцию калийной залежи от проникновения флюидов, определение размеров предохранительных целиков под скважины, и ряд других Выполненные исследования по-

зволяют в практическом плане приступить к решению вопросов, связанных с комплексным освоением минерально-сырьевых ресурсов Соликамской депрессии

ЗАКЛЮЧЕНИЕ

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

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

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

3 С помощью численных расчетов модельных задач показано, что деформации коллекторов определяются соотношением их упругих свойств и свойств вмещающих пород, а также геометрическими характеристиками залежей Относительное уплотнение продуктивного слоя увеличивается при более слабых вмещающих породах, а также с ростом отношения R/H При этом степень влияния упругих свойств вмещающих пород на уплотнение коллекторов и оседания земной поверхности уменьшается с увеличением размеров коллектора Идеализация геометрии коллекторов в виде прямолинейных пластов с постоянной глубиной залегания дает незначительный прирост расчетных сдвижений массива и оседаний земной поверхности Установлено, что для расчета деформаций горного массива можно использовать показатели средневзвешенного пластового давления и не учитывать неравномерность давления, обусловленного работой отдельных добывающих скважин

4 Применение модифицированной шатровой модели для расчета деформаций коллекторов обеспечивает представительные результаты для большого разнообразия горно-геологических условий и физико-механических свойств продуктивных пород Произведено внедрение данной модели в конечно-элементный пакет "AN-SYS" Численные расчеты модельных задач показали, что в центральной части нефтяных и газовых месторождений выполняются условия одномерного уплотнения,

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

Для ряда месторождений Западной Сибири, УНГКМ, АГКМ и региона ВКМКС по результатам компрессионных испытаний определены параметры «шатровой» модели, которые могут быть использованы для расчетов НДС продуктивных объектов данных месторождений

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

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

7 Выполнены оценочные расчеты магнитуд возможных сейсмических событий при отработке Уньвинского нефтяного месторождения Результаты показывают возможность активизации разломных структур при операциях нагнетания флюида для поддержания пластового давления Расчетные значения магнитуд не превышают 1 0-1 5, что является предельно возможной величиной для рассмотренных условий Сейсмические события подобной интенсивности не оказывают заметного влияния на поверхностные, подземные объекты и геологическую среду региона месторождения

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

8 Основные результаты работы были использованы для обоснования мер охраны подрабатываемых объектов и создания геодинамических полигонов на ряде нефтяных месторождений Западной Сибири, севера Пермского края, Уренгойском и Астраханском газоконденсатных месторождениях Результаты исследований вошли в нормативный документ - «Инструкцию по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе Верхнекамского месторождения калийно-магниевых солей»

Основные положения диссертации опубликованы в следующих работах:

1 Кашников Ю А, Якушина Е М, Ашихмин С Г Деформирование скального массива по системам трещин // Известия ВУЗов Горный журнал - 1992 - № 3 -С 75-80

2 Кашников Ю А, Ашихмин С Г Численная модель для расчета сдвижений горных пород при добыче нефти - В кн Проблемы механики горных пород -Санкт-Петербург - 1997 -С 193-198

3 Кашников Ю А , Ашихмин С Г, Замотан В Б , Ашихмин В Н Геодинамические полигоны, сдвижение горных пород и техногенные землетрясения при разработке нефтяных месторождений - В кн Проблемы геодинамической безопасности -Санкт-Петербург - 1997 -С245-247

4. Кашников Ю А, Ашихмин С Г Расчет сдвижений горных пород при разработке нефтяных месторождений Часть 1 Численная модель деформирования нефтяного коллектора П Маркшейдерский вестник -1998 -№ 1 -С 44-46

5 Кашников Ю А, Ашихмин С Г, Селезнев Е.А Расчет сдвижений горных пород при разработке нефтяных месторождений Часть 2 Расчет НДС горного массива при отработке в режиме упругой энергии // Маркшейдерский вестник - 1998 - № 2 -С 33-35.

6 Кашников Ю А, Ашихмин С Г Влияние добычи нефти в упругом режиме на изменение НДС горного массива Часть 1 Анализ инструментальных наблюдений Модель деформирования нефтяного коллектора под нагрузкой // Физико-технические проблемы разработки полезных ископаемых - 1998 -№5 —С 71-80

7 Ашихмин С Г, Кашников Ю А Влияние количества рядов добывающих скважин на оседания земной поверхности на конечной стадии упругого водонапорного режима // Известия ВУЗов Нефть и газ -1998 -№6 - С 26-31

8 Кашников Ю А, Ашихмин С Г Влияние добычи нефти в упругом режиме на изменение НДС горного массива Часть 2 // Физико-технические проблемы разработки полезных ископаемых -1999 -№3 -С 51-57

9 Кашников Ю А, Ашихмин С Г Моделирование возникновения техногенных сейсмических явлений при добыче углеводородного сырья // Маркшейдерский вестник -1999 -№ 1 -С 21-25

10 Кашников Ю А, Ашихмин С Г, Катошин А Ф, Селезнев Е А Изменение НДС горного массива при добыче нефти в упругом режиме // Нефтяное хозяйство -1999 -№ 8 - С 30-33

11. Кашников Ю А, Ашихмин С Г. Оценка техногенной сейсмической опасности при отработке нефтяного месторождения // Геодинамика и напряженное состояние недр земли -Новосибирск - 1999 - С402-408

12 Кашников Ю А , Ашихмин С Г Катошин А Ф Изменение геодинамической обстановки при разработке нефтяного месторождения // Нефтяное хозяйство - 2000 -№ 6 -С 28-32

13 Кашников Ю А, Ашихмин С Г Влияние добычи нефти на изменение НДС горного массива Часть 3 Техногенная активизация разломных структур // Физико-технические проблемы разработки полезных ископаемых — 2000 - № 3 - С 54-63.

14 Кашников Ю А, Ашихмин С Г, Гладышев С В, Калугин А В Комплексная оценка геодинамической безопасности при разработке месторождений углеводородного сырья И Тезисы докладов международной научно-практической конференции «Геоэкология и современная геодинамика нефтегазоносных регионов» - Москва. -2ООО -С 143

15 Кашников Ю.А, Ашихмин С Г., Одинцов С Л, Постнов А В., Рожков В Н Численное моделирование и контроль геомеханических процессов при отработке АПСМ // Тезисы докладов международной конференции «Проблемы добычи и переработки нефти и газа в перспективе международного сотрудничества ученых Каспийского региона» -Астрахань -2000 - С 18-19

16 Кашников Ю А, Ашихмин С Г., Гуляев Н Ю, Гришко С В Геомеханический анализ и геодезический мониторинг деформационных процессов при добыче углеводородов П Тезисы докладов Всероссийской конференции «Природные резервуары углеводородов и их деформации в процессе разработки нефтяных месторождений». - Казань -2000 - С 61-62

17 Ашихмин С Г Численная модель для прогаоза напряженно-деформированного состояния массива рыхлых и скальных пород // Известия ВУЗов Горный журнал -2000 -№4 -С 52-57

18 Кашников Ю А, Ашихмин С Г, Гладышев С В , Калугин А В Геомеханические и геодинамические аспекты разработки месторождений углеводородов Западной Сибири - В кн.. «Геодинамическая и экологическая безопасность при освоении месторождений газа, его транспортировке и хранении» - Санкт-Петербург -2001 -С 71-78

19 Кашников Ю А , Ашихмин С Г, Гладышев С В , Саврасов И Ф, Одинцов С JI, Постнов А В Предварительный прогноз и контроль геодинамических процессов при разработке АГКМ - В кн «Геодинамическая и экологическая безопасность при освоении месторождений газа, его транспортировке и хранении». - Санкт-Петербург-2001 -С 79-86

20 Кашников Ю А, Ашихмин С Г, Одинцов С JI, Постнов А.В Техногенные геодинамические процессы при разработке АГКМ // Газовая промышленность -2002 -№ 1 -С 81-83

21 Кашников Ю А, Ашихмин С Г, Назаров А Ю, Кашников О Ю , Терентьев Б В Экспериментальные исследования влияния пластового давления на фильтраци-онно-емкостные характеристики терригенных коллекторов месторождений севера Пермского края ÍÍ Геология, геофизика и разработка нефтяных и газовых месторождений - 2007. - № 1 - С 41-50

22 Кашников Ю А, Ацшхмин С Г Механика горных пород при разработке месторождений углеводородного сырья - М Недра, 2007 - 486 с

23 Ашихмин С Г, Кашников Ю А, Киселевский Е В , Одинцов С JI, Гетманов И В , Калугин А В Прогноз и мониторинг оседаний земной поверхности при разработке газоконденсатных месторождений России // Труды XIII Международного конгресса по маркшейдерскому делу -Будапешт, 2007

24 Ашихмин С Г Прогноз параметров уплотнения коллекторов и деформаций горного массива при разработке месторождений углеводородов // Геология, геофизика и разработка нефтяных и газовых месторождений -2007 -№10 - С 40-43

25 Аишхмин С Г. Исследование компрессионных свойств коллекторов Шерш-невского нефтяного месторождения // Геология, геофизика и разработка нефтяных и газовых месторождений -2007 -№10 - С 43-45

26 Ашихмин С Г. Особенности методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений углеводородов//Маркшейдерия и недропользование -2008 -№2 -С 38-41

27 Ашихмин С Г Влияние различных факторов на напряженно-деформированное состояние горных массивов при разработке месторождений нефти и газа // Известия ВУЗов Горный журнал - 2008 - № 4

Подписано в печать 25 06 08 Формат 60X90/16 Набор компьютерный Тираж 150 экз Объем 2,5 уч-изд п л Заказ № 8082008

Издательство

Пермского государственного технического университета 614600, г Пермь, Комсомольский пр , 29, к 113 тел (342) 219-80-33

Содержание диссертации, доктора технических наук, Ашихмин, Сергей Геннадьевич

ВВЕДЕНИЕ.

Гл. 1 СОВРЕМЕННОЕ СОСТОЯНИЕ И АКТУАЛЬНОСТЬ ПРОБЛЕМЫ РАСЧЕТА НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ГОРНЫХ ПОРОД НА МЕСТОРОЖДЕНИЯХ НЕФТИ И ГАЗА.

1.1 Явление деформаций горных пород при разработке месторождений углеводородов.

1.2 Геодинамические процессы, сопровождающие добычу нефти и газа.

1.3 Опыт инструментальных наблюдений за оседаниями земной поверхности.

1.4 Методы расчета напряженно-деформированного состояния горных массивов на месторождениях нефти и газа.

1.5 Краткая характеристика объектов исследований.

Выводы.

Гл.2 АНАЛИЗ ИНСТРУМЕНТАЛЬНЫХ НАБЛЮДЕНИЙ ЗА СДВИЖЕНИЯМИ ЗЕМНОЙ ПОВЕРХНОСТИ НА МЕСТОРОЖДЕНИЯХ НЕФТИ И ГАЗА.

2.1 Основные требования к мониторингу деформационных процессов при разработке месторождений углеводородов.

2.2 Наблюдательные станции на нефтяных месторождениях, территориально совмещенных с Верхнекамским месторождением калийно-магниевых солей.

2.3 Результаты наблюдений за сдвижением земной поверхности на Западно-Сургутском нефтяном месторождении.

2.4 Наблюдательная станция на Уренгойском газоконденсатном месторождении.

2.5 Геодинамический полигон Астраханского газоконденсатного месторождения.

Выводы.

Гл.З ПРИМЕНЯЕМЫЕ МЕХАНИЧЕСКИЕ МОДЕЛИ ГОРНЫХ ПОРОД И

ИХ ПАРАМЕТРИЧЕСКОЕ ОБЕСПЕЧЕНИЕ.

3.1 Модифицированная шатровая модель.

3.2. Модель деформирования скальных пород по системам трещин.

3.3 Упругие и прочностные свойства пород объектов исследований.

3.4 Исследование компрессионных характеристик коллекторов.

3.5 Физико-механические свойства горных пород в массиве.

Выводы.

Гл.4 ПРОГНОЗ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ГОРНЫХ ПОРОД ПРИ РАЗРАБОТКЕ МЕСТОРОЖДЕНИЙ УГЛЕВОДОРОДОВ.

4.1 Особенности деформирования насыщенных пористых сред.

4.2 Расчет уплотнения продуктивных объектов.

4.3 Численная реализация модельных представлений.

4.4 Влияние различных факторов на параметры уплотнения коллекторов и оседания земной поверхности.

4.5 Прогноз напряженного состояния горных пород на основе модифицированной шатровой модели.

Выводы.

Гл.5 ДЕФОРМИРОВАНИЕ ПОРОД НА КОНТАКТАХ БЛОКОВЫХ СТРУКТУР И ОЦЕНКА ИНТЕНСИВНОСТИ ТЕХНОГЕННЫХ СЕЙСМИЧЕСКИХ ЯВЛЕНИЙ.

5.1 Проблемы расчета напряженно-деформированного состояния больших объемов горного массива.

5.2 Деформирование пород на контактах блоковых структур при добыче нефти и газа.

5.3 Численная модель активизации разломных структур при добыче углеводородов.

5.4 Оценка техногенного сейсмического риска при отработке нефтяного месторождения.

Выводы.

Гл.6 ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ РЕЗУЛЬТАТОВ

ИССЛЕДОВАНИЙ.

6.1 Прогноз оседаний земной поверхности на нефтяных месторождениях Западной Сибири.

6.2 Численное моделирование процессов сдвижения на Уренгойском газоконденсатном месторождении.

6.3 Прогноз напряженно-деформированного состояния горного массива Астраханского газоконденсатного месторождения.

6.4 Анализ влияния добычи нефти на безопасность разработки Верхнекамского месторождения калийно-магниевых солей.

6.4.1 Прогноз оседаний земной поверхности на нефтяных месторождениях региона ВКМКС.

6.4.2 Влияние добычи нефти на напряженно-деформированное состояние калийной залежи.

Выводы.

Введение Диссертация по наукам о земле, на тему "Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа"

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

Общепризнанно, что один из наиболее значимых видов техногенного воздействия на недра связан с добычей нефти и газа. Разработка нефтяных и газовых месторождений и связанные с ними изменение пластового давления, различные виды воздействия на залежь для поддержания пластового давления и повышения нефтеотдачи, нарушают природное напряженно-деформированное состояние недр, создавая предпосылки для возникновения сильных и даже катастрофических природно-техногенных явлений, которые приводят к деформациям горного массива и земной поверхности, повреждениям и авариям систем и объектов обустройства, а также скважин и коммуникаций. Широко известны такие примеры, как нефтяное месторождение Уиллмингтон (США, Калифорния), разработка которого началась в 1926 году и к 1967 году оседание земной поверхности над месторождением достигло 9 метров; нефтяное месторождение Экофиск в Северном море, где оседание морского дна составляет около 6 метров; нефтяные месторождения в районе озера Маракаибо (Венесуэла), где оседание земной поверхности составило около 4 метров, и другие. Также при извлечении флюида (воды, нефти или газа) известны многочисленные случаи сейсмических явлений. По количеству выделяемой энергии сейсмические толчки, наблюдаемые при разработке месторождений нефти и газа, значительно превышают аналогичные явления при остальных видах воздействия на недра, а в ряде случаев приближаются по силе к природным землетрясениям, характерным для региона месторождений.

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

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

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

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

Основные исследования производились на ряде нефтяных месторождений Западной Сибири, нефтяных месторождениях севера Пермского края, Уренгойском и Астраханском газоконденсатных месторождениях. Наличие высотной жилой застройки, сложного и высокотехнологичного оборудования, а также уникальных природных объектов обуславливает высокую актуальность рассматриваемых проблем применительно к данным месторождениям. Исследования выполнялись на кафедре "Маркшейдерское дело, геодезия и геоинформационные системы" Пермского государственного технического университета. В работах по производству и обработке инструментальных наблюдений на месторождениях принимали участие сотрудники кафедры Е.С.Богданец, С.В.Гришко, В.Г.Букин. Всем им автор выражает искреннюю признательность.

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

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

В задачи исследований входило:

- провести анализ результатов инструментальных наблюдений за сдвижением земной поверхности на месторождениях нефти и газа;

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

- провести аналитические исследования методов расчета уплотнения коллекторов при снижении исходного пластового давления;

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

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

Основные защищаемые научные положения, разработанные лично автором:

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

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

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

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

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

Научная новизна работы:

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

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

Реализация работы. Установленные на основе прогнозных расчетов параметры напряженно-деформированного состояния горных пород и земной поверхности использовались для обоснования мер охраны подрабатываемых объектов и создания геодинамических полигонов, которые были внедрены на ряде нефтяных месторождений Западной Сибири, севера Пермского края, Уренгойском и Астраханском газоконденсатных месторождениях. Результаты исследований вошли в нормативный документ - «Инструкцию по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе Верхнекамского месторождения калийно-магниевых солей».

Апробация работы. Основные положения диссертации и результаты исследований докладывались и обсуждались на следующих совещаниях, конференциях и конгрессах: международной конференции «Проблемы геодинамической безопасности» (Санкт-Петербург, 1997); XI Российской конференции по механике горных пород (Санкт-Петербург, 1997); международной конференции «Геодинамика и напряженное состояние недр Земли» (Новосибирск, 1999); международной научно-практической конференции «Геоэкология и современная геодинамика нефтегазоносных регионов» (Москва, 2000); международной конференции «Проблемы добычи и переработки нефти и газа в перспективе международного сотрудничества ученых Каспийского региона» (Астрахань, 2000); III Международном рабочем совещании «Геодинамическая и экологическая безопасность при освоении месторождений газа, его транспортировке и хранении» (Санкт-Петербург, 2001); XIII международном конгрессе по маркшейдерскому делу (Будапешт, 2007), на заседаниях ученого совета ПермГТУ, на технических советах ООО «Лукойл-Пермь», ООО «Лукойл-Нижневолжскнефть», ОАО «Сургутнефтегаз», ООО «Архангельскгеолдобыча», ООО «Юганскнефтегаз», ООО «Уренгойгазпром», ООО «Астраханьгазпром». и

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

Структура и объем работы. Диссертация состоит из введения, шести глав, заключения и изложена на 315 страницах машинописного текста, включая 108 рисунков, 25 таблиц и список литературных источников из 177 наименований.

Заключение Диссертация по теме "Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика", Ашихмин, Сергей Геннадьевич

8. Основные результаты работы были использованы для обоснования мер охраны подрабатываемых объектов и создания геодинамических полигонов на ряде нефтяных месторождений Западной Сибири, севера Пермского края, Уренгойском и Астраханском газоконденсатных месторождениях. Результаты исследований вошли в нормативный документ - «Инструкцию по созданию наблюдательных станций и производству инструментальных наблюдений за процессами сдвижения земной поверхности при разработке нефтяных месторождений в регионе Верхнекамского месторождения калийно-магниевых солей».

ЗАКЛЮЧЕНИЕ

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

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

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

3. С помощью численных расчетов модельных задач показано, что деформации коллекторов определяются соотношением их упругих свойств и свойств вмещающих пород, а также геометрическими характеристиками залежей. Относительное уплотнение продуктивного слоя увеличивается при более слабых вмещающих породах, а также с ростом отношения R/H. При этом степень влияния упругих свойств вмещающих пород на уплотнение коллекторов и оседания земной поверхности уменьшается с увеличением размеров коллектора. Идеализация геометрии коллекторов в виде прямолинейных пластов с постоянной глубиной залегания дает незначительный прирост расчетных сдвижений массива и оседаний земной поверхности. Установлено, что для расчета деформаций горного массива можно использовать показатели средневзвешенного пластового давления и не учитывать неравномерность давления, обусловленного работой отдельных добывающих скважин.

4. Применение модифицированной шатровой модели для расчета деформаций коллекторов обеспечивает представительные результаты для большого разнообразия горно-геологических условий и физико-механических свойств продуктивных пород. Произведено внедрение данной модели в конечно-элементный пакет "АК8У8". Численные расчеты модельных задач показали, что в центральной части нефтяных и газовых месторождений выполняются условия одномерного уплотнения, т.е. прирост эффективных горизонтальных напряжений соответствует теоретическому значению коэффициента бокового давления. В краевой части продуктивных пластов девиаторная часть тензора эффективных напряжений при падении пластового давления растет интенсивнее, чем в режиме одномерного уплотнения, однако при этом предельное состояние не достигается и коллектора также деформируются в режиме уплотнения.

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

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

7. Выполнены оценочные расчеты магнитуд возможных сейсмических событий при отработке Уньвинского нефтяного месторождения. Результаты показывают возможность активизации разломных структур при операциях нагнетания флюида для поддержания пластового давления. Расчетные значения магнитуд не превышают 1.0-1.5, что является предельно возможной величиной для рассмотренных условий. Сейсмические события подобной интенсивности не оказывают заметного влияния на поверхностные, подземные объекты и геологическую среду региона месторождения.

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

Библиография Диссертация по наукам о земле, доктора технических наук, Ашихмин, Сергей Геннадьевич, Пермь

1. Авершин С.Г. Расчет деформаций массива горных пород под влиянием подземных разработок. Изд.ВНИМИ. Ленинград. - 1960г. - 87с.

2. Алексеев Р.И., Коровин Ю.И. Руководство по вычислению и обработке результатов количественного анализа. М., Атомиздат, 1972, 72 с.

3. Ашихмин С.Г. Влияние различных факторов на напряженно-деформированное состояние горных массивов при разработке месторождений нефти и газа // Изв. вузов. Горный журнал. 2008. - № 4.

4. Ашихмин С.Г. Исследование компрессионных свойств коллекторов Шершневского нефтяного месторождения // Геология, геофизика и разработка нефтяных и газовых месторождений. 2007. - №10. - с.43-45.

5. Ашихмин С.Г. Особенности методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений углеводородов // Маркшейдерия и недропользование. 2008. - №2. - с.38-41.

6. Ашихмин С.Г. Прогноз параметров уплотнения коллекторов и деформаций горного массива при разработке месторождений углеводородов // Геология, геофизика и разработка нефтяных и газовых месторождений. 2007. - №10. - с.40-43.

7. Ашихмин С.Г. Численная модель для прогноза напряженно-деформированного состояния массива рыхлых и скальных пород // Изв. вузов. Горный журнал. 2000. - № 4. - с.52-57.

8. Баженова O.K., Гурлин Ю.Н. и др. Геология и геохимия нефти и газа. М.: Академия 2004. - 415с.

9. Ю.Баклашов И.В. Деформирование и разрушение породных массивов. -М.:Недра, 1988.-271с.11 .Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. М.: Недра, 1984. - 211с.

10. Барях A.A., Константинова С.А., Асанов В.А. Деформирование соляных пород. Екатеринбург, УрО РАН. 1996.- с.91-107.

11. Барях А. А., Кудряшов А. И, Еремина Н. А., Грачева Е. А. Оценка влияния разработки нефтяного месторождения на геодинамическое состояние недр // ФТПРПИ. — 1998. — № 2.

12. Бенявски З.Т. Управление горным давлением. М.:Мир.-1990.-254с.

13. Булычев Н.С. Механика подземных сооружений. М.-.Недра, 1982.-270с.

14. Булычев Н.С. Механика подземных сооружений в примерах и задачах. -М.:Недра, 1989.-270с.

15. Букринский В.А., Рашковская Я.З., Фролов Е.Ф. Совершенствование методов маркшейдерских работ и геометризация недр. М.: Недра, 1972.-383с.

16. Бурштейн JI.C. Статические и динамические испытания горных пород. Л., Недра, 1970, 181 с.

17. Вялов С.С. Реологические основы механики грунтов М.: Высшая школа, 1978.-317с.

18. Ганыпин В.Н., Стороженко А.Ф., Ильин А.Г. и др. «Измерение вертикальных смещений сооружений и анализ устойчивости реперов». М.: Недра, 1981.-215с.

19. ГОСТ 21153.2-84 Породы горные. Методы определения предела прочности при одноосном сжатии. М., 1985, 10 с.

20. ГОСТ 21153.3-85 Породы горные. Методы определения предела прочности при одноосном растяжении. М., 1985, 14 с.

21. ГОСТ 28985-91 Породы горные. Методы определения деформационных характеристик при одноосном сжатии. М., 1991, 19 с.

22. Гриценко А.И., Зотов Г.А. Научно-прикладные геодинамические проблемы разработки месторождений природного газа. -В кн. Проблемы геодинамической безопасности. ВНИМИ.-С.-Петербург. - 1997.-е. 186-193.

23. Гудман Р. Механика скальных пород. М.: Стройиздат, 1987. - 232с.

24. Дияшев Р.Н., Костерин A.B., Скворцов Э.В. Фильтрация жидкости в деформируемых нефтяных пластах. Изд. Казанского математического общества. -1999. 238с.

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

26. Ержанов Ж.С., Сагинов A.C. и др. Ползучесть осадочных пород. Изд. «Наука» Казахской ССР. Алма-Ата.-1970.-208с.

27. Желтов Ю.П. Деформации горных пород. М.:Недра, 1966.- 198с.

28. Иванова М.М., Дементьев Л.Ф., Чоловский И.П. Нефтегазопромысловая геология и геологические основы разработки месторождений нефти и газа. -М.:Недра, 1985, 383 с.

29. Инструкция по защите рудников от затопления и охране объектов в условиях Верхнекамского месторождения калийных солей. С-Петербург, 2004.

30. Инструкция по наблюдениям за сдвижением горных пород и земнойповерхности при подземной разработке угольных и сланцевых месторождений. -М.: Недра, 1989г. -96с.

31. Инструкция по наблюдениям за сдвижением горных пород и земной поверхности при подземной разработке рудных месторождений. М.: Недра, 1988г.- 112с.

32. Инструкция по нивелированию I, И, III, IV классов. М.: ЦНИИГАиК, 2004г.

33. Иофис М.А., Шмелев А.И. Инженерная геомеханика при подземных разработках. М.: Недра, 1985.-248с.

34. Калинин Э.В., Панасьян Л.Л., Артамонова Н.Б. Методика исследований напряженно-деформированного состояния массивов пород при оценке инженерно-геологических условий нефтяных и газовых месторождений // Инженерная геология.- 2006. №11. - с.51-57.

35. Касьянова H.A. Оценка и учет геодинамических рисков при проектировании обустройства морских нефтегазовых месторождений // Нефтяное хозяйство. 2005г.- №6. с.38-42.

36. Касьянова H.A., Э.В.Соколовский, С.В.Шимкевич. Результаты прогноза аварий скважин и порывов трубопроводных систем по геодинамическому фактору. Нефтяное хозяйство.-1998.-№9.-с.75-77.

37. Кашников Ю.А. Научные основы разработки методов прогноза параметров деформирования подрабатываемых скальных массивов мощных крутопадающих рудных месторождений. Автореф.дисс.на соис.уч.ст.докт.техн.наук. Москва, 1992. -35с.

38. Кашников Ю.А., Ашихмин С.Г. Численная модель для расчета сдвижений горных пород при добыче нефти. В кн.: Проблемы механики горных пород. Санкт-Петербург, 1997.-с.193-198.

39. Кашников Ю.А., Ашихмин С.Г. Расчет сдвижений горных пород при разработке нефтяных месторождений. Часть 1. Численная модель деформирования нефтяного коллектора // Маркшейдерский вестник. 1998. - № 1. - с.44-46.

40. Кашников Ю.А., Ашихмин С.Г., Селезнев Е.А. Расчет сдвижений горных пород при разработке нефтяных месторождений. Часть 2. Расчет НДС горного массива при отработке в режиме упругой энергии // Маркшейдерский вестник.1998. № 2. - с.33-35.

41. Кашников Ю.А., Ашихмин С.Г. Влияние количества рядов добывающих скважин на оседания земной поверхности на конечной стадии упругого водонапорного режима // Изв.Вузов. Нефть и газ. 1998. - № 6. - с.26-31.

42. Ашихмин С.Г. Кашников Ю.А. Влияние добычи нефти в упругом режиме на изменение НДС горного массива. Часть 2. // Физико-технические проблемы разработки полезных ископаемых.- 1999. № 3. - с.51-57.

43. Кашников Ю.А., Ашихмин С.Г. Моделирование возникновения техногенных сейсмических явлений при добыче углеводородного сырья // Маркшейдерский вестник. 1999. - № 1. - с.21-25.

44. Кашников Ю.А., Ашихмин С.Г., Катошин А.Ф., Селезнев Е.А. Изменение НДС горного массива при добыче нефти в упругом режиме // Нефтяное хозяйство.1999. № 8. - с.30-33.

45. Кашников Ю.А., Ашихмин С.Г. Оценка техногенной сейсмической опасности при отработке нефтяного месторождения. В кн.: Геодинамика и напряженное состояние недр земли. Новосибирск, 1999. -с.402-408.

46. Кашников Ю.А., Ашихмин С.Г. Катошин А.Ф. Изменение геодинамической обстановки при разработке нефтяного месторождения // Нефтяное хозяйство. 2000. - № 6. - с.28-32.

47. Кашников Ю.А., Ашихмин С.Г. Влияние добычи нефти на изменение НДС горного массива. Часть 3. Техногенная активизация разломных структур // Физико-технические проблемы разработки полезных ископаемых.- 2000. № 3. - с.54-63.

48. Кашников Ю.А., Ашихмин С.Г., Одинцов C.JL, Постнов A.B. Техногенные геодинамические процессы при разработке АГКМ // Газовая промышленность. -2002.-№ 1. с.81-83.

49. Кашников Ю.А., Ашихмин С.Г. Механика горных пород при разработке месторождений углеводородного сырья. М.: Недра, 2007. - 486с.

50. Кашников Ю.А., Кашников О.Ю., Рахимкулов Р.С.Информационно-экспертная система безопасной эксплуатации межпромыслового нефтепровода // Нефтяное хозяйство. 2001г. №6. - с.72-77.

51. Кашников О.Ю., Кашников Ю.А., Круглов Ю.И. и др. Информационно-экспертная система эксплуатации участка магистрального газопровода // Газовая промышленность. 2002г. №9. - с.76-78.

52. Кашников Ю.А., Якушина Е.М., Ашихмин С.Г. Деформирование скального массива по системам трещин // Изв. вузов. Горный журнал. 1992. - № 3. - с.75-80.

53. Козырев A.A., Панин В.И., Савченко С.Н. Сейсмичность при горных работах. Апатиты: КНЦ РАН, 2002г.

54. Копнин В.И. Верхнекамское месторождение калийно-магниевых и каменных солей и природных рассолов. Известия вузов. Горный журнал. 1995. - № 6. - с. 1042.

55. Костров Б.В. Механика очага тектонического землетрясения. М., Наука, 1975.

56. Кратч Г. Сдвижение горных пород и защита подрабатываемых сооружений. -М.: Недра, 1978. 494с.

57. Кудряшов А.И. Минерально-сырьевые ресурсы Пермского края. Энциклопедия. Изд.»Красная площадь». 2006г. - 464с.

58. Кузьмин Ю.О. Современная геодинамика и оценка reo динамического риска при недропользовании. М.: Агентство экономических новостей, 1999.-220с.

59. Кузьмин Ю.О., Жуков B.C. Современная геодинамика и вариации физических свойств горных пород. М.: Изд.Моск.госуд.горного университета. -2004г.-262с.

60. Куликов В.В. Совместная и повторная разработка рудных месторождений. -М.: Недра, 1971.-328с.

61. Курленя М.В., Серяков В.М., Еременко A.A. Техногенные геомеханические поля напряжений. Новосибирск: Наука, 2005. - 264с.

62. Линьков A.M. Численное моделирование сейсмических и асейсмических событий в геомеханике // ФТПРПИ, 2005. №1. - с. 19-33.

63. Мазницкий A.C., Середницкий Л.М. Влияние параметров упругости пород на уплотнение коллектора и оседание земной поверхности при разработке нефтяных месторождений. Нефтяное хозяйство.-1991.- №6,- С. 14-16.

64. Методические рекомендации к «Указаниям по защите рудников от затопления и охране объектов в условиях Верхнекамского месторождения калийных солей. С-Петербург, 2004.

65. Механика горных пород применительно к проблемам разведки и добычи нефти. М.: МИР " Эльф-Акитен ", 1994 г.

66. Моги К. Предсказание землетрясений. -М.:Мир, 1988.-382с.

67. Муллер P.A. О некоторых дифференциальных и интегральных свойствах кривых оседаний и горизонтальных сдвижений мульды сдвижения // ФТПРПИ, 1981. №2. - с.129-132.

68. Николаевский В.Н. Геомеханика и флюидодинамика. -М.гНедра, 1996.-447с.

69. Никонов А.А Землетрясения. Прошлое, современность, прогноз. М.: Знание, 1984.-192с.

70. Ножин А.Ф. Зотеев В.Г. Расчет НДС трещиноватого скального массива вокруг горных выработок // ФТПРПИ, 1978. №5. - с.9-14.

71. Перепеличенко В.Ф. и др. Разработка нефтегазоконденсатных месторождений Прикаспийской впадины. М.: Недра, 1994.-364с.

72. Повышение эффективности освоения газовых месторождений крайнего севера. М.:Наука, 1997. - 655с.

73. Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на угольных месторождениях. С.-Петербург: Изд. ВНИМИ, 1998г.-291 с.

74. Правила промышленной безопасности при освоении месторождений нефти на площадях залегания калийных солей. Постановление Госгортехнадзора России от 4.02.2002г за №8.

75. Петухов И.М., Батугина И.М. Геодинамика недр.-М.:Недра, 1996.-217с.

76. Ползучесть осадочных пород. // Ж.С.Ержанов, А.С.Сагинов и др. Изд. «Наука» Казахской ССР. Алма-Ата. - 1970. - 208с.

77. Природные резервуары углеводородов и их деформации в процессе разработки нефтяных месторождений. Тезисы докладов конференции. Казань -2002г.

78. Проблемы геодинамической безопасности. ВНИМИ. -С.-Петербург.-1997.-с.66-71.100 .Проблемы освоения месторождений Уренгойского комплекса. М.: Недра, 1998. -464с.

79. Проблемы освоения Астраханского газоконденсатного месторождения. -Астрахань: АНИПИГАЗ, 1999. -332с.

80. Ю5.Сашурин А.Д. Сдвижение горных пород на рудниках черной металлургии. Изд.УрО РАН. Екатеринбург. - 1999. - 268с.

81. Ю7.Сидоров В.А., Хитров В.М., Кузьмин Ю.О. Концепция геодинамической безопасности освоения углеводородного потенциала недр России. М.: изд.МГГРИ.- 1998г.

82. Ю8.Система обеспечения геодинамической и экологической безопасности при проектировании и эксплуатации объектов ТЭК. Методические рекомендации. С.Петербург.: ВНИМИ, 2001. - 86с.

83. Ю9.Смирнова М.Н. Возбужденные землетрясения в связи с разработкой нефтяных месторождений (на примере Старогрозненского землетрясения) //

84. Влияние инженерной деятельности на сейсмический режим. М.: Наука, 1977, с. 128141.

85. ПО.Сонич В.П., Черемисин H.A., Батурин Ю.Е. Влияние снижения пластового давления на фильтрационно-емкостные свойства пород // Нефтяное хозяйство.-1997.- №9.- С.52-57.111 .Справочная книга по добыче нефти. Под ред. Ш.К. Гиматудинова. -М., Недра, 1974.

86. Справочник по механическим и абразивным свойствам горных пород нефтяных и газовых месторождений. М.:Недра, 1984. - 208с.

87. Спутниковый мониторинг земной поверхности / Антонович K.M., Карпик А.П., Клепиков А.Н.// Геодезия и картография. 2004. - №1.

88. Ставрогин А.Н., Протосеня А.Г. Пластичность горных пород. М.: Недра, 1979.-301с.

89. Тер-Микаелян K.JI. О связи удельного коэффициента отпора и модуля деформации горных пород с коэффициентом крепости по М.М.Протодьяконову // Труды Гидропроекта, 1974. -вып.ЗЗ. с.116-122.

90. Пб.Токунов В.И., Рылов Е.А., Поляков Г.А. и др. Гидроразрыв пласта на Астраханском ГКМ. // Газовая промышленность,- 1998. №3. - с.47-48.

91. ТЭО пробной и опытно-промышленной эксплуатации месторождения нефти им.Архангельского. Отчет о работе. ОАО «ПермНИПИнефть». Рук.работы Голубев Б.М. Фонды ООО «ЛУКОЙЛ-Пермь». 2001г.

92. Турчанинов И.А., Иофис М.А., Каспарьян Э.В. Основы механики горных пород. М.: Недра.-1989.-332с.

93. Фадеев А.Б. Метод конечных элементов в геомеханике. М.: Недра, 1987.221с

94. Фадеев А.Б. Прочность и деформируемость горных массивов. -М.:Недра, 1979, 270 с.

95. Филатов В.В., Кассин Г.А., Попов Б.А. Геофизические исследования на Верхнекамском месторождении калийно-магниевых солей. Известия вузов. Горный журнал. 1995. - № 6. - с.150-161.

96. Шадрин А.Г. Теория и расчет сдвижений горных пород и земной поверхности. Красноярск, Изд.Красноярского гос.ун-та, 1990. - 199с.

97. Щелкачев В.Н., Лапук Б.Б. Подземная гидравлика. М.,"Динамика", 2001.

98. Щугорев В.Д. ООО «Астраханьгазпром». Перспективы развития // Наука и технология углеводородов, 2001. №4. - с.7-9.

99. Ягунов А.С., Серяков В.М. Характер и закономерности упругого деформирования массива с несколькими выработками // ФТПРПИ, 1984. №5. -с.24-33.

100. Ященко В.Р. Геодезические исследования вертикальных движений земной коры. М.: Недра, 1989.-192с.

101. Aydan О, Shimizu Y. and Kawamoto Т. The Anisotropi of Surfage Morphology Charakteristics of Rock Discontinuities // Rock Mec. Rock Enginering, 1996.-vol.29-№l.-P.47-29.

102. Aziz, K. and Settari, A.: Petroleum Reservoir Simulation, Applied Science Publ., New York, 1979,476 p., paperback ed. Published by authors, 2001.

103. Bagheri, M. Ali and Settari, A. : "Modeling of Geomechanics in Naturally Fractured reservoirs", Paper SPE 93083, SPE Reservoir Simulation Symposium, Houston, TX, Jan. 30-Feb. 2, 2005 (to appear in SPE RE).

104. Bandis S.H, A.C. Lumsden, N.R. Barton. Fundamentals of rock joint deformation. Int. J. Rock Mech. No. 6, pp. 249-268, 1983.

105. Barton N.R., S.N. Bandis, K. Bakhtar. Strength, deformation and conductivity coupling of rock joints. Int. J. Rock Mech. No. 36, pp. 121-140, 1985.

106. Barton, N.; Bandis, S. (1982): Effects of Block Size on the Shear Behaviour of Jointed Rock. Proc. of the 23rd U.S. Symp. on Rock Mech., Issues in Rock Mech., Berkeley, California 1982. New York: A.I.M.E.

107. Biot M.A. General theory of tree-dimensional consolidstion. J.Appl.Phys., 1941,12, 155-164.

108. Boade R.R., Chin L.Y., Siemers W.T. Forecasting of Ekofisk Reservoir Compaction and Subsidence by Numerical Simulation. Journal of Petroleum technology. July 1989, pp. 723-728.

109. Charlez F. P. Rock Mechanics. Volume 1,2. Petroleum applications. Editions Technip. 27 rue Ginoux 75737 Paris cedex 15. 1997.

110. Chin L.Y., Boade R.R., Nagel N.B., Landa G.H. Numerical Simulation of Ekofisk Reservoir Compaction and Subsidence: Treating the Mechanical Behavior of the Overburden and Reservoir. Eurock'94, pp.787-794.

111. Chin L.Y., Boade R.R., Prevost J.H., Landa G.H. Numerical Simulation of Shear-Induced Compaction in the Ekofisk reservoir. Int. J. Rock Mech. Vol.30, No.7, pp. 1193-1200, 1993.

112. Chin, L. Y. and Boade, R. R., Full-Field, 3-D Finite-Element Subsidence Model for Ekofisk, Third North Sea Chalk Symposium, Copenhagen, June 11-12, 1990.

113. Chin, L.Y., and Thomas, L.K.: "Fully Coupled Analysis of Improved Oil Recovery by Reservoir Compaction," SPE paper 56753 presented at the 1999 SPE Annual Technical Conference and Exhibition, Houston, Tx, Oct. 3-6, 1999.

114. Drewes H, Henneberg H. Деформации геодезических сетей на месторождении нефти в районе оз.Маракаибо // Allgemeine VermessungsNachrichten, 1980, V.7,№10, р.386-389.

115. Gambolati G, Ricceri G, Bertoni W, Brighehti G. Mathematical Simulation of the Subsidence of Ravenna. Water Resourses Research, Vol. 27, No. 11, pp. 2899-2918,

116. Geertsma J. A Basic Theory of Subsidence Due to Reservoir Compaction: The Homogeneous Case. Verh. Con. Ned. Geol. Mijnbouwk. Gen., Vol. 28, pp. 43-62, 1973.

117. Grant U.S. Subsidence of the Wilmington Oil Field. Cal. Calif.Div.Mines Bull., 1954,170, p.19-24.

118. Grasso J.R., Feigner B., Seismicity Induced by Gas Production: Lithology Correlated Events, Induced Stresses & Deformation, 134 Pure & Applied Geophys. 427, 1990.

119. Grasso J.R., Plotnikova L.M., Nutaev L.M. Bossu R. The Three M-7 Gazli Earthquakes, Usbekistan, Central Asia: The Largest Seismic Energy Releases by Human Activity, Abstracts. XXI Gen Ass. Int. Union Geodesy & Geophys/ A363, 1995.

120. Hvorslev M.J. (1937): Uber die Festigkeitseigenschaften gestoerter bindiger Boden. Ingeniorvidenskabelige Skifter. A. No. 45, Copenhagen, Denmark, 1937.

121. Jaeger G.C. Friction of Rocks and stability of rock slopes. Geotechnique 21/2 (1971)

122. Johnson J.P., Rhett D.W., Siemers W.T. Rock Mechanics of the Ekofisk Reservoir in the Evaluation of Subsidence. Journal of Petroleum technology. July 1989, pp.717-722.

123. Leichnitz W. Mechanische Eigenschaften von Felstrennflachen im direkten Scherversuch. Veroeff.des Inst.fuer Bodenmechanik und Felsmechanik der TH Karlsruhe, Heft 89.-1981.

124. Litwiniszyn J. Theorien und Modelluntersuchungen über Bodenbewegungen. Berlin, 1957.- S.202-209.

125. Maxwell S.C., Young R.P., Bossu R. Microseismic Loggin of the Ekofisk Reservoir. EUROCK-98. Rock Mechanics in Petrolium Industrie. Trondheim, 1998. Pp.3 87-393.

126. Mortia N., Whitfill D.L., Nygaurd O., Bale A. A quick method to determine subsidence reservoir compaction and in-situ stress induced by reservoir depletion. Journal of Petroleum technology. Januar 1989, pp. 38-41.

127. Muller W. Die Berechnung der Bewegungen und Spannungen des Gebirges vom Abbau bis zuz Tagesoberflache nach der Methode endlichen Elemente. Diss. Berlin, 1973.

128. Pattillo P.D,.Kristiansen T.G,.Sund G.V, Kjelstadli R.M. Reservoir Compaction and Seafloor Subsidence at Valhall. Eurock'98, pp.377-386.

129. Pande G.N. & Pietruszczak, St. (1982): "Reflecting surface" model for soils. Proc. 1st Int. Symp. Num. Models in Geomech, Zurich, 1982, pp 50-64.

130. Pande G.N. & Sharma R.G. (1983): Multi-Laminate model of clays: a numerical evaluation of the influence of rotation of the principal stress axes. Int. J. Num. Analy. Methods Geomech, Vol. 7, 1983, pp 397-418.

131. Perzyna P. (1966): Fundamental problems in viscoplasticity. Advances in Applied Mechanics, Academic press, New York, Vol. 9, 1966, pp 244-368.

132. Pietruszczak, St. & Mroz, Z. (1983): On hardening anisotropy of K0-consolidated clays. Int J. Num. Analy. Methods Geomech. Vol.7, 1983, pp 19-38.

133. Plischke B. Finite element analysis of compaction and subsidence-Experience gained from several chalk fields. Eurock' 94. 1994 Balkema, Rotterdam, s.795-801

134. Reddish D.J., Yao X.L., Waller M.D. Computerised prediction of subsidence over oil and gas fields. Eurock'94, pp.621-630.

135. Roest J.P.A. & W.Kuilman. Geomechanical analysis of small earthquakes at the Eleveld gas reservoir. Eurock '94.-S.573-580.

136. Roest J.P.A., Mulders F.M.M. Overview Modelling Gas Production-Induced Seismicity Mechanisms. Eurock 2002.- pp.333-338.

137. Roscoe K.H. & Burland J.B. (1968): On the generalized stress-strain behavior of "wet" clay. Engineering Plasticity, J.Heyman & F.A. Leckie (eds), Cambridge University Press, 1968, pp 535-609.

138. Roscoe K.H., Schofield, A.N. & Thurairajah, A. (1963): "Yielding of clays in states wetter than critical". Geotechnique, Vol. 13, 1963, pp 211-240.

139. Simpson D., Leith W. The 1976 and 1984 Gazli, USSR, Earthquakes, Were They Induced? Vol. 75 (1985), №5, Bull. Seismological Soc. Am. 1465-1468.

140. Schofield, A. & Wroth, P. (1968): Critical state soil mechanics. McGraw-Hill, London, 1968.

141. Sulak R. M., Thomas, L. K., and Boade, R.R., 3D Reservoir Simulation of Ekofisk Compaction Drive, JPT (October 1991), pp 1272-1278.

142. Sulak R.M. Ekofisk Field : The First 20 Years. Journal of Petroleum technology. October 1991, pp. 1269-1269.

143. Teufel L.W., Rhett D.W., Farrel H.E. Effect of Reservoir Depletion and Pore Pressure Drawdown on in Situ Stress and Deformation in the Ekofisk Field, North Sea. 32nd U.S. Symposium on Rock Mechanics, University of Oklahoma, July 10-12, 1991.

144. Wittke, W.: Rock Mechanics, Theory and Applications with case histories, Springer-Verlag, Berlin, Heidelberg, New York, London, Paris, Tokio, Hongkong, Barcelona, 1990a.

145. Zienkiewicz, O.C.: The finite element method. McGraw-Hill, London, 1977.

146. Zienkiewicz, O.C. & Pande, G.N (1977): Some useful forms of isotropic yield surfaces for soil and rock mechanics: Finite Elements in Geomechanics. G. Gudehus (ed.), Hohn Wiley & Sons, 1977, pp 179-190.

Информация о работе

Ашихмин, Сергей Геннадьевич
доктора технических наук
Пермь, 2008
ВАК 25.00.20

Диссертация

Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Научные основы методов прогноза напряженно-деформированного состояния горных пород при разработке месторождений нефти и газа - тема автореферата по наукам о земле, скачайте бесплатно автореферат диссертации

Похожие работы