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

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

А.И. Калашник1'2, С.Н. Савченко1'3, В.И. Панин1, Э.В. Каспарьян1'4,

0.В. Смирнова1

1 Горный институт КНЦ РАН, Апатиты

2 Филиал Санкт-Петербургского государственного инженерно-экономического университета, кафедра информационных систем в экономике, Апатиты

3 Кольский филиал Петрозаводского государственного университета, кафедра прикладной математики, Апатиты

4 Естественно-технический факультет МГТУ, кафедра горного дела

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

Abstract. Methodical approaches have been stated considering the shelf gas-oil structure as an integrated geological object being an open dynamic nonlinear natural and engineering system. Mathematical methods for modeling gas-oil shelf structures have been developed taking into account the sea shelf geological structure properties, gas-oil deposits specifics as well as the main processes at their development and exploitation. Multivariant mathematical modeling has been carried out for an integrated geological object - gas-condensate field. The structure of geomechanical monitoring has been developed over the geological environment of the whole region, where the construction of this industry corresponding facilities and long pipelines is planned. Methodology has been developed for structuring geomechanical and geophysical data being the background for geological and geomechanical data bases on more than 130 gas-oil fields and perspective marine shelf provinces in the Barents region.

1. Введение

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

1

и

/

* к-.. __" з/ / /

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

2. Концептуальная модель

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

место пространственно-временная локализация событий, 11 Тг %

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

Рис. 1. Режим развития катастрофы в ПТС типа техногенных землетрясении. г г т

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

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

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

3. Математическая модель природно-технической системы - шельфового газоконденсатного

геообъекта

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

первостепенное внимание.

Разработана методика математического моделирования шельфовых нефтегазовых структур

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

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

- глубина залегания геообъекта от дневной поверхности и под дном моря;

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

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

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

На основе этой информации разрабатывается модель для численных расчетов. При этом граница исследуемого геообъекта (при использовании МГЭ) или весь объект (при использовании МКЭ) представляется необходимым числом элементов с соответствующими свойствами. Задаются граничные условия силового воздействия (вес пород и воды, величина действующих субгоризонтальных тектонических сил, боковой отпор, внутреннее давление нефти и газа).

Результатом численных расчетов являются значения напряжений и деформаций по всей исследуемой области (при использовании МГЭ) или в центрах элементов (МКЭ). При этом могут быть определены величины удельной энергии упругого деформирования в зонах наибольшей концентрации напряжений. Распределения напряжений и деформаций отображаются в виде графиков или изолиний.

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

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

Для протяженных месторождений с отношением мощности d к протяженности L d/L < 0.01 возможно решение задачи в двумерной постановке с использованием МГЭ. В силу малых отношений параметров рельефа на дне моря и протяженности месторождения приходится пренебрегать влиянием рельефа, полагая поверхность дна моря горизонтальной, так как изменение масштаба в одном из направлений изменяет кривизну поверхности модели, от величины которой существенно зависят величины действующих напряжений.

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

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

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

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

ar = 3qa2/4h2; w = qa4/64D,

где a = L/2, q - равномерно распределенная нагрузка на поверхности дна моря, учитывающая вес воды и горных пород, D - жесткость кровли, зависящая от модуля Юнга E и коэффициента Пуассона v пород и глубины залегания, w - прогиб. Однако эти формулы позволяют оценить величины напряжений и прогиба лишь в центре кровли.

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

u = 1/2E[o2i + oV 2 va\ 02 +

где y - объемный вес пород, ai, a2 - главные напряжения, действующие в кровле.

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

u = u0

где u0 - начальная энергия упругого деформирования, а - величина, определяемая по результатам изменения напряженного состояния кровли на первых этапах отработки месторождения; sh t - синус гиперболический времени (t) деформирования.

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

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

1) на верхней горизонтальной поверхности действуют равномерно распределенные нормальные усилия d'y, имитирующие давление толщи воды;

2) в массиве пород действуют равномерно распределенные гравитационные силы oy = yy и ax = Л/y, где у- объемный вес пород, Л = v/( v-1) - боковой отпор пород, y - текущая координата по глубине;

3) на левых и правых границах (имитируя на бесконечности) действуют равномерно распределенные горизонтальные напряжения cfx= T + Àyy;

4) на контуре залежи действуют нормальные напряжения а„, имитирующие внутреннее давление

газа.

Для унификации значений указанные напряжения отнесены к величине действующих в породах массива тектонических напряжений Щ = 50МПа, т.е. задаются в расчетах соответствующими относительными величинами, поэтому в дальнейшем под выражениями <jx, (fx, crn и т.д. будем понимать их отношения ax/\T\, avx /\Т\,а„/\Т\ и т.д.

Для вариантов расчетов в качестве переменных величин принимались h/L, cfx и an, рассмотрены отношения: h/L = 1.0; 0.5; 0.2; 0.1; стшх = -1; -0.75; -0.5; -0.25; 0; ап = -1.5; -1.0; -0.75; -0.5; -0.25; 0. Отношение a°x = 0 означает, что в массиве действуют только массовые силы, обусловленные весом пород и воды, а an = 0 - отсутствие давления газа.

Таким образом разработанная модель позволяет имитировать варианты последовательной отработки газоконденсатных залежей с соотношением их глубины расположения h к протяженности L от 0.1 до 1 и различным исходно-напряженным состоянием массива пород.

На данном этапе выполнены расчеты по 120 вариантам различных соотношений указанных переменных величин с использованием метода граничных элементов, реализация которого осуществлена в постановке первой задачи теории упругости в условиях плоскодеформированного состояния (Савченко, Калашник, 2001). По результатам расчетов построены схемы и графики распределения

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

Параметры неотектсни-

ческих процессов

ГЬраметры флюидо-игазовых эмиссий

Параметры электрошг-нитных испей

Гсризон-тальныг Наклоны

Параметры сейсмоаку> тическойэмиссии

Верти-кальныг

Значения гравитацион-

Объемные

Параметры напряженно-дефсрмированного состояния пород вмеш^юшрго массива

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

Значения и направления действия тектонических сил

.Давление газа, газо-ковденсатаипи нефти

БАЗА ДАННЫХ И ПАРАМЕТРОВ

Обработка Критерии

Озртировка Группировка Агрегирование Вычисления Сопоставление, сравнение /

Е

СЦенка текушрго и прогнозных состояний ПГС

Г Цжнятиеуправляющкрешенийипревентивныхмероприягий у

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

4. Структура мониторинга ПТС

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

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

- мониторинг на стадии инженерных изысканий;

- мониторинг при обустройстве месторождений;

- мониторинг при строительстве основных объектов: терминалов, насосных станций, подземных хранилищ нефти и газа и т.п.;

- мониторинг при эксплуатации месторождений;

- мониторинг при эксплуатации основных объектов;

- мониторинг при транспортировании углеводородного сырья;

- мониторинг на стадии вывода из эксплуатации объектов и демонтажа оборудования;

- мониторинг при консервации месторождений.

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

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

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

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

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

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

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

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

- оценка динамики состояния геологической среды;

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

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

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

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

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

- режим нормальной эксплуатации;

- режим аварийной ситуации;

- режим консервации или вывода из эксплуатации объектов.

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

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

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

4. Создание информационного банка - многомерной базы данных

Решение этой задачи основано на применении современных информационных технологий (Калашник и др., 1998). На платформе СУБД MS Access созданы и поддерживаются в актуальном состоянии базы геомеханических и геолого-геофизических данных и параметров более чем по 130 нефтегазовым месторождениям и перспективным провинциям шельфа морей Баренцрегиона (Калашник, Калашник, 2002).

Первоначально концептуально базы данных базировались на реляционной модели данных, разработанной Э.Ф. Коддом (Дейт, 2001). За последнее время осуществляется переход к многомерной модели данных, целью которой является прежде всего анализ данных, наряду с традиционным выполнением транзакций. Разработана методология структуризации геомеханических и геолого-геофизических данных по шельфовым структурам, заключающаяся в том, что многомерная модель рассматривает данные либо как факты с соответствующими численными параметрами, либо как текстовые измерения, которые характеризуют эти факты. При этом многомерная модель данных имеет три принципиально важных особенности применения, связанных с проблематикой анализа данных:

- базы данных (по сути дела - хранилища данных) интегрируют для анализа информацию, получаемую

из нескольких разнопрофильных источников;

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

данных в поисках общих тенденций (зависимостей);

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

зависимостей и связей в базах данных.

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

5. Выводы

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

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

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

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

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

- численные модели для исследования шельфовых углеводородных структур;

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

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

3) Разработаны принципы организации мониторинга состояния геологической среды при разработке месторождений углеводородного сырья.

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

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

- концептуальное проектирование автоматизированной системы применительно к шельфовым нефтегазовым структурам;

- информационные банки геолого-геофизических и геомеханических данных по шельфу Баренцрегиона.

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

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

Работа выполнена при финансовой поддержке РФФИ (проект № 02-05-64736). Литература

Savchenko S.N. and Kalashnik A.I. Estimate of the stress state of rocks in the vicinity of gas-condensate field on the sea shelf. Journal of Mining Science. Kluwer Academic /Plenum Publishers, v.37, N 4, p.384-391, 2002. Дейт К.Д. Введение в системы баз данных. М., Вильяме, 1071 е., 2001.

Калашник А.И., Калашник A.A. К анализу мирового рынка нефти и газа. Нефть и газ арктического

шельфа 2002, Арктикшельф, т.6, с.40-42, 2002. Калашник А.И., Козырев A.A., Бессонов И.И., Чуркин O.E. Информационные технологии в решении задач освоения минерально-сырьевых ресурсов Кольского полуострова. Горный журнал, №4, с.79-81, 1998.

Калашник А.И., Савченко С.Н. Оценка и прогноз изменений напряженно-деформированного состояния горных пород в окрестности Штокмановского газоконденсатного месторождения. Нефть и газ арктического шельфа 2002, Арктикшельф, т.2, с.81-90, 2002. Мельников H.H., Козырев A.A., Калашник А.И., Шабаров А.Н. Ееодинамическая безопасность при обустройстве и эксплуатации месторождений углеводородов шельфа Баренцева и Карского морей. Арктикшельф, т.5, с.55-73, 2002. Мельников H.H., Козырев A.A., Савченко С.Н., Панин В.И., Мальцев В.А. Прогноз и профилактика

техногенных землетрясений с позиций нелинейной геодинамики. ФТПРПИ, №4, с.1-16, 2001. Савченко С.Н., Калашник А.И. К вопросу геомеханического обеспечения обустройства и эксплуатации месторождений углеводородов шельфа Баренцева моря. Проблемы разработки месторождений полезных ископаемых и освоения подземного пространства Северо-запада России. Кольский научный центр, т.3, с.34-41, 2002. Савченко С.Н., Калашник А.И. Оценка напряженного состояния горных пород в окрестности

газоконденсатного месторождения на морском шельфе. ФТПРПИ, №4, с.49-56, 2001. Сейсмичность при горных работах. Коллектив авторов. Кольский научный центр, 325 е., 2002. Система обеспечения геодинамической и экологической безопасности при проектировании и эксплуатации объектов ТЭК. СПб, ВНИМИ, 86 е., 2001.