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

Кудрявцева, Елена Олеговна

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

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

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

КУДРЯВЦЕВА Елена Олеговна

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

Специальность 25.00.10 — Геофизика, геофизические

методы поисков полезных ископаемых

Автореферат диссертации на соискание ученой степени кандидата геолого-минералогических наук

САНКТ-ПЕТЕРБУРГ 2010

1 8 НОЯ 2010

004612635

Работа выполнена в государственном образовательном учренедении высшего профессионального образования Санкт-Петербургском государственном горном институте имени Г.В.Плеханова (техническом университете).

Научный руководитель -

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

Путиков Олег Федорович

Официальные оппоненты: доктор технических наук

Ермохин Константин Михайлович,

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

Ведущая организация - Институт нефтегазовой геологии и геофизики им. А.А.Трофимука СО РАН, г.Новосибирск.

Защита диссертации состоится 17 ноября 2010 г. в 16 ч на заседании диссертационного совета Д 212.224.01 при Санкт-Петербургском государственном горном институте имени Г.В.Плеханова (техническом университете) но адресу: 199106 Санкт-Петербург, 21-я линия, д.2, ауд.4312.

С диссертацией можно ознакомиться в библиотеке Санкт-Петербургского государственного горного института.

Автореферат разослан 15 октября 2010 г.

Арчегов Владимир Борисович

УЧЕНЫЙ СЕКРЕТАРЬ диссертационного совета канд. геол.-минерал. наук, доцент

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

Актуальность темы. Дифференциально-нормированный метод электроразведки (ДНМЭ) в последние годы стал одним из ведущих методов в нефтегазопоисковом геофизическом комплексе. Поисковым признаком выступают аномалии вызванной поляризации (ВП), возникающие над залежами углеводородов (УВ), в частности, вследствие развития в зоне ореолов рассеяния вторичной эпигенетической пиритизации. Основоположниками метода являются Ю.П. Булашевич, Б.И. Геннадинник, В.А. Комаров, В.В. Кормильцев, B.C. Моисеев, A.B. Постельников и многие другие.

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

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

Задачи исследования:

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

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

- исследовать связь аномалий ВП с зонами вторичной сульфидизации в отношении их интенсивности и планового положения;

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

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

нормированным методом электроразведки (ДНМЭ), выполненных ООО «Сибирская геофизическая научно-производственная компания» (СГНПК) в 2006-2009 гг.; фактический материал, полученный при проведении экспериментов по исследованию постоянной времени спада ВП на базе лаборатории кафедры геофизики Санкт-Петербургского государственного университета в 2008-2009 гг. В качестве основных методов исследования использовались натурные эксперименты, математическое моделирование, компьютерная обработка и интерпретация экспериментальных и теоретических данных. Автор принимала участие в обобщении, анализе и интерпретации геолого-геофизических данных ДНМЭ по площадям Калининградского и Краснодарского регионов, шельфов Каспийского и Баренцева морей, лично проводила эксперименты в Санкт-Петербургском государственном университете, обработку и анализ данных, полученных в их результате, разработала теоретическое обоснование математического моделирования образования зон сульфидизации.

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

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

- выявлены связи распределения концентрации пирита над залежами углеводородов в зависимости от их положения и возраста и проведено сопоставление полученных результатов с аномалиями ВП, полученными по данным ДНМЭ, что позволило повысить степень достоверности прогнозирования залежей углеводородов;

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

Защищаемые положения.

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

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

3. Аномалии ВП с высокими значениями постоянной времени спада ВП (до 100 мс - 3 с при обычных значениях 20 - 80 мс) характерны для расформированных в геологическом прошлом залежей УВ, над которыми происходило частичное окисление зерен вторичных сульфидов.

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

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

Учет влияния разломов во многих случаях позволяет более точно определить плановое положение залежей УВ при анализе аномалий ВП.

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

расформированными залежами УВ.

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

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

Основные положения диссертации были представлены и докладывались: на ме>вдунаР0ДНЬ1Х научно-практических конференциях EAGE в г. Геленджике в 2006 г., 2007 г., 2008 г., 2009 г. и в г. Монпелье (Франция) в 2009 г., на научно-практической конференции ИрГТУ в 2006 г, на международных геофизических научно-практических семинарах Санкт-Петербургского государственного горного института в 2006 г., 2007 г., 2008 г. и 2009 г., в Институте нефтегазовой геологии и геофизики СО РАН (г. Новосибирск) в 2007 и 2008 г.

Публикации. По теме диссертации опубликовано 7 работ, из них I работа в изданиях, входящих в перечень ВАК Минобрнауки России, имеется 1 патент на изобретение.

Объем и структура работы. Диссертация состоит из введения, 4-х глав, заключения, содержит 126 страниц машинописного текста, включая 52. рисунка, 5 таблиц, список литературы из 78 наименований.

Благодарности. Автор выражает благодарность своему научному руководителю, д.г.-м.н., профессору О.Ф. Путикову. Также автор благодарит д.г.-м.н., генерального директора ООО «Сибирская геофизическая научно-производственная компания» П.Ю. Легейдо за поддержку исследовательской работы и ценные советы и замечания, зав. кафедрой ГФХМР СПГГИ (ТУ) д.г.-м.н., профессора A.C. Егорова за предоставленную возможность всесторонней апробации диссертации, обсуждение структуры работы и ценные замечания, д.г.-м.н. профессора Н.О. Кожевникова и Е. Ю. Антонова, доктора Паула Викена.

ЗАЩИЩАЕМЫЕ ПОЛОЖЕНИЯ И ИХ ОБОСНОВАНИЯ

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

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

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

////////

+ Ре1" + 2ОН~ Ре^ +Н20

.......1.......1.......1.......

пиритизация

(30, СИ.

'°а НсЯ4 +Са* СаСО-. + Т

| 4

иальцитизация

1—1 1 <еггг> г Т з

+н2о

Рис. 1. Схема образования зоны вторичных сульфидов: 1 - залежь УВ; 2 - зона образования вторичных сульфидов, 3 -направление потока сероводорода

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

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

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

АН;

У:-/,:,

гХ-'

-¿.¿-А.

2,0

Н2 Hl

III (D„ v3, Сз) II (D2, v2, C2)

I (Dlf v,, Ci)

Hi,o

r l'o 0 Го г

Рис. 2. Схема расположения зон сульфидизации над нефтегазовой

залежью в трехслойном геологическом разрезе:

1 - нефтегазовая залежь радиуса г0; 2,3 - зоны сульфидизации; I(Db vb СО - нижний слой горных пород и его параметры: Di, vi, Q - коэффициент диффузии, вертикальная скорость квазиконвекции и концентрация сероводорода соответственно; II (D2, v2, С2) - средний слой и его параметры; III (D3, v3, С3) - верхний слой и его параметры; Нь Н2 -координаты границ слоев; Н!)0 - координата подошвы нижней зоны сульфидизации; Н2,0 - координата кровли верхней зоны сульфидизации; Чь Чз - концентрация сульфида железа в нижней и верхней зонах сульфидизации соответственно.

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

периметр нефтегазовой залежи, т.е. горизонтальное кольцо радиуса г0, единица длины которого выделяет за единицу времени количество сероводорода <3. Тогда распределение концентрации сероводорода и сульфида железа в цилиндрических координатах -г, г (рис. 2) с учетом осевой симметрии задачи подчиняется системе нелинейных дифференциальных уравнений в частных производных:

V гдг Ц & Ц* °т Ц т А >п Ц дг °

Э2С2 |Э2С, 1ЭС, V2дс2_(2)

дг2 г дг Б2 & ¿>2 дт

= ftCi ~ ^ ). Hw<z< Нь

= (qXmax - q3 ), H2<z< Я 2.0,

где т - время, мс; рь р3 - константы кинетики реакции образования сульфида железа в нижней и верхней зонах сульфидизации соответственно; qlimax, q3,max - максимально возможная концентрация сульфида железа в нижней и верхней зонах сульфидизации соответственно, т|(х) - единичная функция, определяемая соотношениями:

8(х) - дельта - функция Дирака, определяемая соотношениями:

си

s(x)=(o: iio/ Z5^-1-

—сс

Система нелинейных дифференциальных уравнений в част-

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

= (6) "77" = Рз Сз!цзуг [яЪтодг- Я^] ■ (7)

Интегрируя уравнения (6) и (7) стационарных концентраций сульфида при нулевом начальном условии для концентраций сульфида железа ql(1) и qзa), получим:

Чг ' = <7з,шах

Н10<1<Н{ (8)

Н2<1< Н20

(9)

С помощью программы «Пирит», созданной в ООО «СГНПК», было проведено изучение процессов вторичного эпигенеза и их связи с наблюдаемыми аномальными эффектами на примере Северо-Гуляевского нефтегазового месторождения (шельф Баренцева моря). С учетом данных бурения была принята трехслойная модель разреза (рис. 3). Пласт 1 и 3 - слабопроницаемые (значимая концентрация пирита); вертикальная скорость квазиконвекции и коэффициент диффузии приняты равными 1-Ю'10 м/с и 1-Ю'10 м2/с соответственно. В более проницаемом пласте 2 эти характеристики имеют значения 1-Ю"10 м/с и 1-Ю"10 м2/с соответственно. При расчетах были приняты следующие значения параметров объекта, соответствующие реальной геологической ситуации: вертикальное удаление от залежи 2 км, радиус залежи 5 км, геологическое время ее формирования - 250 млн. лет, Н^ЮО м; Н1=500 м; Н2=2000 м; Н2,0=2300 м.

Рис. 3. Трехслойная модель для расчета зон сульфидизации: г0 - радиус залежи; 1, 3 - слабопроницаемые пласты; 2 - относительно проницаемый пласт.

Площадь исследований изучена методом ДНМЭ по ряду профилей. В результате установлен контур аномалии ВП (рис. 4). Во внутриконтурной области значения поляризуемости 9% и более против 3 - 6% за контуром залежи. На рис. 5 представлены графики расчета концентрации пирита на расстоянии 50, 150, 250, 350, 450, 2050, 2150 и 2250 м по вертикали от залежи. При сопоставлении закономерности распределения коэффициента поляризуемости и расчетной концентрации пирита прослеживается корреляция между этими параметрами: максимальные значения обоих параметров наблюдаются непосредственно над залежью, а снижение происходит за ее пределами.

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

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

При расчете концентрации пирита получены следующие результаты: наибольшая концентрация пирита (60 % от максимально возможной) наблюдается на расстоянии 2250-2750 м по вертикали от нефтяной залежи, на уровне геохимического барьера. Следует отметить, что в силу небольшого геологического времени формирования залежи процесс насыщения пиритом еще не закончен. Поскольку время формирования газовой залежи еще меньше, чем нефтяной, то пирит образовался в еще меньших количествах (20 % от максимально возможной концентрации) на расстоянии 450-750 м для газовой залежи, т.е. на уровне геохимического барьера.

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

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

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

Ч'ЧтЛ-''*0)' (10)

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

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

где I - координата по оси х расположения вертикального разлома, м; С?- мощность линейного источника (единицы его длины), имеющего координаты г = 0, х = х0; к - коэффициент поглощения сероводорода твердой фазой горных пород; Б - коэффициент диффузии сероводорода во вмещающих породах; v - эффективная скорость квазиконвекции сероводорода во вмещающих породах, м/с; Д v - приращение эффективной скорости квазиконвекции в зоне разлома; 5(х-1) - дельта-функция Дирака. Граничными условиями уравнения являются нулевые решения на бесконечности:

Э С Э2С Эг2 + Э*2

V Д\S\x-l)

Ъ Ъ

С ->0, С -»0. (12)

Для учета влияния разломов в юго-западной части Северо-Гуляевского месторождения проведен расчет концентрации пирита в приразломной части месторождения.

Был создан набор графиков распределения концентрации пирита для различных геологических условий (глубин залегания залежи, положения разломов относительно залежи и геологического времени образования залежи). В относительно проницаемых породах влияние разломов сказывается следующим образом: если залежь ограничена разломом или разлом находится в непосредственной близости от залежи (в данном случае не более 100 м от границы залежи), то сероводород мигрирует вдоль разлома, расширяя при этом аномалию ВП (рис. 7, графики 1 и 2 соответственно). Если разлом находится на большем расстоянии от залежи, то его влияние незначительно (график 3). Представленные графики распределения концентрации пирита соответствуют залежи с полушириной 500 м, находящейся на глубине 2000 м и временем формирования залежи 40 млн.лет. Отметим, что аналогичная картина отмечается для различных глубин с различным временем образования залежи.

В слабопроницаемых породах картина может меняться существенно (рис. 8). Влияние разломов в зависимости от их близости к залежи сказывается по-разному. Непосредственно над залежью пирит образуется медленнее и его концентрация не достигает максимума при небольшом времени образования залежи. Поток сероводорода устремляется вдоль разлома. Чем ближе разлом к залежи, тем шире зона максимального насыщения пиритом, соответственно шире аномалия ВП (графики 1, 2, 3). При удалении разлома от залежи за ее пределами наблюдается снижение концентрации пирита с последующим возрастанием ее непосредственно над разломом (графики 4, 5). При достаточном удалении от залежи над разломом пирит не образуется (график 6).

Для Северо-Гуляевского месторождения расчет концентрации пирита под влиянием разлома был выполнен при тех же характеристиках среды, что и расчет без учета влияния разломов. Ис-

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

Еще один пример теоретического расчета концентрации пирита приведен для одной из залежей небольших размеров, расположенной в Калининградской области (рис. 10). В результате работ ДНМЭ была получена аномалия ВП над сводовой частью структуры, ограниченной с юга разломом. Расположение пикетов с аномальными значениями поляризуемости не дают четкого ответа, по какой же изолинии можно провести контур ВНК (-1990 или -1980 м) и каковы в итоге размеры залежи. Расчет концентрации пирита был сделан по двум линиям (1-1 и 2-2). В случае, если полуширина

Рис. 10. Карта распределения коэффициента поляризуемости на структурной основе отражающего горизонта ОГ III (кровля

ордовика):

1 - контур аномалии ВП по данным ДНМЭ; 2 - изолинии по ОГ III, м; 3 - зона тектонических нарушений.

залежи 200 м и она ограничена разломом (линия 1-1на рис. 10, график 1на рис. 11), концентрация пирита максимальна над всей залежью и выходит за пределы нее со стороны разлома примерно на 400 м. Если же разлом не примыкает вплотную к залежи (на расстоянии 200 м от нее), то распределение пирита выглядит довольно неожиданно (линия 2-2 на рис. 10, график 2 на рис. 11). Миграция сероводорода происходит преимущественно по разлому, обедняя содержание пирита непосредственно над залежью. Максимум концентрации пирита отмечается, начиная от края залежи в сторону разлома более чем на 400 м. В случае, если бы контур ВНК проходил по изолинии -1990 м, то аномальные значения поляризуемости должны были бы отмечаться на ПК 18-22 ПР 026307, а на ПК 10-14 этого профиля поляризуемость должна быть ниже. Если контур ВНК ограничен изолинией -1980 м, то становится понятным увеличение поляризуемости на ПК 10-14 ПР 026307, поскольку они расположены вблизи разлома, где концентрация пирита максимальна. Таким образом, можно сделать достаточно четкий вывод о положении контура ВНК по изолинии -1980 м и определить размер залежи как примерно 500 м.

3. Аномалии ВП с высокими значениями постоянной времени спада ВП (до 100 мс - 3 с при обычных значениях 20 -80 мс) характерны для расформированных в геологическом прошлом залежей УВ, над которыми происходило частичное окисление зерен вторичных сульфидов.

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

Рис. 4. Структурная карта по отражающему горизонту 1А (кровля

карбонатного комплекса С2-3-Р1): 1 - контур аномалии ВП по данным ДНМЭ; 2 - изолинии по отражающему горизонту 1А по данным сейсморазведки, м; 3 -контур ВНК; 4 - тектонические нарушения.

1000 2000 3000 4000 5000 6000 7000 8000 9000 10000

Расстояние от центра залежи, м

Рис. 5. Сопоставление распределения коэффициента поляризуемости и расчетной концентрации пирита на Северо-Гуляевском месторождении: щ^В - нефтяная залежь; ^^ - график концентрации пирита в верхней части разреза; - график концентрации пирита в покрышке залежи; - график коэффициента поляризуемости.

3 ®4

Рис. 6. Структурная карта по ОГ сЬ2 (Краснодарский край): 1 - контур аномалии ВП по данным ДНМЭ; 2 - тектонические нарушения; 3 - изолинии глубин (м, по данным сейсморазведки); 4 -скважины с притоком нефти; 5 - скважины с притоком газа.

В к л ауКФ

Рис. 7. Теоретический расчет концентрации пирита для модели месторождения полушириной 500 м в относительно проницаемых породах (у = МО"10 м/с, Б = 1-10~9 м2/с, с1у = 1-Ю"10 м/с2): при расстоянии вертикального разлома от границы залежи; 1 -0 м; 2 - 100 м; 3 - 200 м.

Расстояние, м

Рис. 8. Теоретический расчет концентрации пирита для модели месторождения радиусом 500 м в слабопроницаемых породах

(у - МО"13 м/с, О - МО"12 м2/с, йч- МО"12 м/с2): при расстоянии вертикального разлома от границы залежи: 1 -0 м; 2 - 100 м; 3 - 200 м; 4 - 300 м; 5 - 400 м; 6 - 500 м.

Расстояние, м

Рис. 9. Теоретический расчет распределения концентрации

пирита для Северо-Гуляевского месторождения I - положение разлома; ^В ~ нефтяная залежь.

Рис. 11. Теоретический расчет распределения концентрации пирита над залежью углеводородов в Калининградской области. 1 - по линии 1-1 (разлом вплотную примыкает к залежи); 2 - по линии 2-2 (разлом на расстоянии 200 м от края залежи); 3 - график распределения коэффициента поляризуемости по ПР 016307

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

Постоянная времени спада ВП т0 в соответствии с обобщенными формулами В.В. Кормильцева определяется выражением:

(13)

где Кп =^?T/^()бF - сопротивление перехода, - газовая постоянная, Т - абсолютная температура, F - число Фарадея, 7об -плотность тока обмена на поверхности зерна сульфида, Сдв - емкость единицы поверхности двойного электрического слоя на границе раздела не окисленной и окисленной частей зерна сульфида,

хс = арх ]1Кп , а - радиус неокисленной части зерна пирита, р! -

эффективное удельное электрическое сопротивление неоднородной вмещающей среды.

На рисунке 12 приведены графики зависимости постоянной времени спада ВП т0 и поляризуемости Г| от приведенного сопротивления вмещающей среды хс при различных значениях относительного времени задержки в измерении ВП (&КпСдв). В интервале 0,5 < хсй 2 отмечается достаточно резкое увеличение постоянной времени спада ВП - т0 (примерно в два раза) и относительно небольшое изменение поляризуемости 11 (примерно на 20-40%).

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

Рис. 12. Зависимости постоянной времени спада т0 (кривые 1и 2)

потенциала ВП сферического включения радиуса а и поляризуемости среды т] (кривые 3 и 4) с вкрапленностью таких же включений от приведенного удельного электрического сопротивления вмещающей среды хс = ар/2Яп. Значения постоянных: 7уклСйв =100 (Т3 - время зарядки, с); % = ОД (объемная концентрация вкрапленных сульфидов); КПС(>, для кривых : 1 - 20мс, 2 - 40мс; г/ЯпСав для кривых: 3-0,5; 4-0,1; т - время задержки, мс.

значениях аргумента наблюдается резкое увеличение т0.

На основе проведенных экспериментов был сделан вывод о том, что постоянная времени спада ВП возрастает с увеличением приращения удельного электрического сопротивления среды вокруг поверхности зерна пирита.

Постоянная времени спада ВП в пределах аномалии ВП на

Рис. 13. Зависимость постоянной времени спада от произведения эффективного удельного электрического сопротивления вмещающей среды на радиус шарика.

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

ЗАКЛЮЧЕНИЕ Данная диссертация является самостоятельной, законченной научно-квалификационной работой, в которой обоснованы научные положения, способствующие решению актуальной задачи поиска месторождений нефти и газа. Основные результаты работы:

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

Учет влияния вертикального разлома на распределение

концентрации пирита позволяет выполнить оценку положения и размера залежи углеводородов по данным ДНМЭ.

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

НАИБОЛЕЕ ЗНАЧИМЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Кудрявцева Е.О. Моделирование процессов образования зон повышенной поляризуемости над нефтегазовыми залежами // Записки горного института. 2009. Т. 183. С. 246-250.

2. Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации / Легейдо П.Ю., Иванов С.А., Агеенков Е.В., Кудрявцева Е.О. и др. (РФ) Пат. 2399931 Российская Федерация, МПК G01V3/38, G01V3/06.; Заявитель и патентообладатель ООО «Сибирская геофиз. науч.-произв. комп-я». - №2008141889/28; опубл. 20.09.2010; приоритет 23.10.08.

3. Veeken Р.С.Н. Benefit of the induced polarization geoelectric method to hydrocarbon exploration / P.J. Legeydo, Y.A. Davidenko, E.O. Kudryavceva, S.A. Ivanov, A. Chuvaev // Geophisics, vol.74, N.2 (march-april.2009). - P. 1 -XXXX-12-XXXX.

4. Veeken P. Geoelectric modelling with separation between electromagnetic and induced polarization field components / P. Legeydo, I. Pesterev, Y. Davidenko, E. Kudryavceva, S. Ivanov // First Break, issue 12, vol.27 (december.2009). P. 53-55.

- ° n

РИЦ СПГГИ. 06.10.2010. 3.588 Т. 100 экз. 199106 Санкт-Петербург, 21-я линия, д.2

Содержание диссертации, кандидата геолого-минералогических наук, Кудрявцева, Елена Олеговна

Введение

1. Геолого-геофизические и геохимические предпосылки создания теоретической модели формирования и строения зон сульфидной минерализации над нефтегазовыми залежами. 9 1.1. Возникновение дифференциально-нормированной электроразведки

ДНМЭ (исторический обзор)

1.2. Применение метода ДНМЭ

1.3.Развитие теории возникновения аномалий вызванной поляризации над нефтегазовыми залежами

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

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

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

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

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

3.1. Общая постановка задачи

3.2.Распределение концентрации вторичных сульфидов над линейным источником

4. Моделирование процессов окисления пирита в зонах сульфидной минерализации над нефтегазовыми залежами 83 4.1. Теоретический расчет развития окислительных процессов в зоне сульфидизации и влияния их на параметры ВП над нефтегазовой залежью

4.1.1. Физико-химические процессы формирования зоны сульфидизации и ее окисления

4.1.2. Математическая постановка задачи окисления сферического зерна сульфида (пирита) в зоне сульфидизации над нефтегазовой залежью

4.1.3. Расчет параметров и процессов вызванной поляризации частично окисленного сферического зерна электронопроводящего минерала (пирита)

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

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

Актуальность работы

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

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

Задачи исследования

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

Фактический материал

В основу диссертации положен фактический материал наземных и морских исследований дифференциально-нормированным методом электроразведки (ДНМЭ), выполненных Сибирской геофизической научно-производственной компанией (СГНПК) в 2006-2009 гг.; фактический материал, полученный при проведении экспериментов по исследованию постоянной времени спада ВП на базе лаборатории кафедры геофизики Санкт-Петербургского государственного университета в 2008-2009 гг.

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

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

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

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

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

Защищаемые положения

3. Аномалии ВП с высокими значениями постоянной времени спада ВП (до 100 мс — 3 с при обычных значениях 20 - 80 мс) характерны для расформированных в геологическом прошлом залежей УВ, над которыми происходило частичное окисление зерен вторичных сульфидов.

Практическая значимость работы

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

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

Апробация работы

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

Основные положения диссертации были представлены и докладывались: на международных научно-практических конференциях ЕАвЕ в г. Геленджике в 2006 г., 2007 г., 2008 г., 2009 г. и в г. Монпелье (Франция) в 2009 г., на научно-практической конференции ИрГТУ в 2006 г, на международных геофизических научно-практических семинарах Санкт-Петербургского государственного горного института в 2006 г., 2007 г., 2008 г. и 2009 г., в Институте нефтегазовой геологии и геофизики СО РАН (г. Новосибирск) в 2007 и 2008 г.

Публикации

По теме диссертации опубликовано 7 работ, из них 1 работа в изданиях, входящих в перечень ВАК Минобрнауки России, имеется 1 патент на изобретение.

Благодарности

Автор выражает благодарность за руководство и постоянную методическую помощь при проведении и подготовке исследовательской работы и конструктивные замечания при обсуждении основных результатов исследований своему научному руководителю, д.г.-м.н., профессору О.Ф. Путикову. Автор признателен генеральному директору СГНПК, д.г.-м.н. П.Ю. Легейдо за поддержку исследовательской работы и ценные советы и замечания. Автор благодарен за предоставленную возможность всесторонней апробации диссертации, обсуждение структуры работы и ценные замечания д.г.-м.н. A.C. Егорову, д.г.-м.н. профессору Н.О. Кожевникову, за содействие в участии в международной конференции и в последующих публикациях доктору Паулу Викену, за реализацию идей в программу Е. Ю. Антонову, за поддержку и ценные замечания и советы главному геологу СГНПК С.А. Иванову, С.Ю. Гариной и геологическому отделу СГНПК за помощь в оформлении диссертации.

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Кудрявцева, Елена Олеговна

ЗАКЛЮЧЕНИЕ

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

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

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

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

Библиография Диссертация по наукам о земле, кандидата геолого-минералогических наук, Кудрявцева, Елена Олеговна, Санкт-Петербург

1. Аванпроект «Локальный прогноз нефтегазоносности в условиях Арктического шельфа», Санкт-Петербург, 1998 г. Автор — М.Л. Верба.

2. Агеенков Е.В. Исследование эквивалентностей горизонтально-слоистых поляризующихся сред в дифференциально-нормированном методе электроразведки: Автореф. . канд. геол.-мин. наук. Иркутск, 2004. - 18 с.

3. Баженова O.K. Геология и геохимия нефти и газа: учебник / Под ред. Б.А. Соколова. М.: изд-во Московского университета. Издательский центр «Академия», 2004. - 415 с.

4. Березкин В.М., Киричек М.А., Кунарев A.A. Применение геофизических методов для прямых поисков нефти и газа. М.: Недра. - 1978.

5. Булашевич Ю.П. Расчет поля вызванных потенциалов для рудных тел сферической формыю Изв. АН СССР, 1956. Сер. геофиз., №5

6. Геннадиник Б.И. Принципы описания физико-химических процессов, протекающих в горных породах, и связанных с ними электрохимических явлений, используемых в методе вызванной поляризации. М., Наука, 1980.

7. Геннадиник Б.И. Уравнения и параметры осредненного электромагнитного поля в гетерогенной среде. Изв. АН СССР, сер. Физика Земли, № 6, 1980, с.69 - 75.

8. Геннадиник Б.И., Мельников В.П., Геннадиник Г.В. Теория вызванной электрохимической активности горных пород. Якутск, 1976.

9. Геология нефти и газа Сибирской платформы / Анциферов A.C., Ба-кин В.Е., Варламов И.П. и др. Под ред. А. Э. Канторовича, В. С. Суркова, А. А. Трофимука. М.; Недра, 1981. - 552 с.

10. Голубев B.C., Гарибянц A.A. Гетерогенные процессы геохимической' миграции. М.,Недра,1968.- 192с.

11. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1975. - 1108 с.

12. Давыденко Ю.А. Разработка программно-измерительного комплекса для дифференциально-нормированного метода электроразведки: Автореф. диссертации на соискание ученой степени канд. техн. наук. Новосибирск, 2005.-19 с.

13. Дамаскин Б.Б., Петрий О.А. Электрохимия. «Высшая школа». М., 1987.

14. Диткин В.А., Прудников А.П. Интегральные преобразования и операционное исчисление. М.: Физматгиз, 1961, 524 с.

15. Зорькин Л.М., Карус Е.В. Закономерности формирования и развития геофизических и геохимических полей // Советская геология, 1978, № 11.-С. 18-24.

16. Иванов С.А., Барышев JI.A., Мячев С.Б. Методика интерпретации КМПП // Пояснительная записка к паспорту Чиканского объекта АТЗ (подготовленного). Фонды ФГУГП Иркутскгеофизика, Иркутск, 1994. -С. 33-37.

17. Кирюхин В.А., Коротков А.И., Павлов А.Н. Общая гидрогеология Л.: Недра,1988.-358с.

18. Комаров В. А. Электроразведка методом вызванной поляризации. Л., Недра, 1980. 390 с.

19. Комаров В.А. Геоэлектрохимия. С.-Петербург, изд-во С.-Петерб. унта 1994.-136с.

20. Комплексный анализ данных геохимических поисков месторождений нефти и газа / Под ред. Зорькина J1.M., Петухова A.B. М.: Недра, 1981.

21. Кормильцев В. В. Переходные процессы при вызванной поляризации. М., Наука, 1980.- 112 с.

22. Кормильцев В. В., Поздняков М.В., Сомов В.Ф. О вызванной поляризации песчано-глинистых пород. В кн. : Аппаратура и метод вызванной поляризации с измерением скорости спада. Свердловск, УНЦ АН СССР, 1976, с. 41 -67.

23. Кормильцев В.В. Вызванная поляризация в уравнениях электродинамики. Научные доклады АН СССР, Свердловск, 1981.

24. Кудрявцева Е.О. Моделирование процессов образования зон повышенной поляризуемости над нефтегазовыми залежами // Записки горного института. 2009. Т. 183. С. 246-250

25. Куликов А. В., Шемякин Е. А. Электроразведка фазовым методом вызванной поляризации. М., Недра, 1978. 160 с.

26. Легейдо П. Ю., Мандельбаум М. М., Рыхлинский Н. И. Применение дифференциально нормированной электроразведки на Непском своде // Геология и геофизика, № 4, 1990.

27. Легейдо П.Ю. Теория и технология дифференциально-нормированной геоэлектроразведки и ее применение в нефтегазовой геофизике: Дисс. . д-ра геол.-мин. наук. Иркутск, 1998. - 198 с.

28. Легейдо П.Ю., Мандельбаум М.М., Рыхлинский Н.И. Дифференциально-нормированные методы геоэлектроразведки / Под ред. М.М. Мандельбаума. Методическое пособие. - Иркутск, 1996. - 145 с.

29. Лыков A.B. Теория теплопроводности. М.: Гос. изд-во технико-теорет. лит. 1952. -392с.

30. Метод вызванной поляризации. — Учебно-научный геофизический полигон МГУ. URL: http://alexandrovka.carpet-tile.ru/ (дата обращения 5.10.08)

31. Моисеев B.C., Тараторкин Б.Ф., Шлепнев В.Б. Результативность прогноза залежей углеводородов методом вызванной поляризации в Западной Сибири // Тез. докл.: Международная геофизическая конференция. С.Петербург, 1995.

32. Моисеев B.C. Метод вызванной поляризации при поисках нефтепер-спективных площадей. Новосибирск: Наука, 2002. - 135с.

33. Небрат А.Г., Сочельников В.В., Кисель С.А. К вопросу о применении методов электроразведки при прогнозе характера УВ-насыщения // Геофизика, 2008, №5.-С. 57-58.

34. Никитин A.A. Применение теорий статистических решений при классификации геологических объектов по данным комплекса геофизических методов // «Изв. Вузов. Геология и разведка», 1970, №10. С. 135-144.

35. Никитин A.A. Статистические методы выделения геофизических аномалий.- М.: Недра, 1979. 280 с.

36. Никитин A.A., Хмелевской В.К. Комплексирование геофизических методов: учебник для вузов. Тверь: ООО «Изд. ГЕРС», 2004 -294 с.

37. Основы теории геохимических полей углеводородных скоплений / Старобинец И.С., Петухов A.B., Зубайраев С.Л. / Под ред. A.B. Петухова и

38. И.С. Старобинца. М.: Недра, 1993. - 332 с.

39. Пастухов Н.П. Взаимосвязь геохимических и геофизических полей в надпродуктивных комплексах нефтегазовых месторождений юга Сибирской платформыю: Автореф. . канд. геол.-мин. наук. Москва, 1992. -22 с.

40. Петренко В.И., Петренко Н.В., Зленко В.Я. К вопросу о геолого-физической и геохимической роли газоэвапоригенной влаги // Серия «Нефть и газ». Выпуск 2.

41. Постельников А.Ф. К вопросу о природе вызванной поляризации в осадочных горных породах. -Изв. Вузов, сер. Геология и разведка, 1959, №2.

42. Применение геофизических методов при прямых поисках нефти и газа // Итоги науки и техники / Базовкина И.Г., Корольков Ю.С., Кунарев A.A. и др.- М., из-во ВИНИТИ, 1978.

43. Прямой поиск углеводородов геофизическими методами / Мандель-баум М. М., Пузырев Н. Н., Рыхлинский Н. И. и др. М.: Наука, серия "Академические чтения", 1988.- 160 с.

44. Путиков О.Ф. Основы теории нелинейных геоэлектрохимических методов поисков и разведки. С.Петербург, изд-во С-Петерб. гос. горного института. 2009. 534с.

45. Путиков О.Ф. Полярографический каротаж. С.-Петербург, изд-во С.-Петерб. гос. ун-та. 2000. -500с.

46. Рокитянский И. И. Лабораторное изучение вызванной поляризации осадочных пород. Изв. АН СССР, сер. Геофизика, 1957, № 2.

47. Рокитянский И. И. О природе вызванной поляризации ионопроводя-щих сред //Изв. АН СССР, сер. геоф., 1959. С. 1055- 1060.

48. Светов Б.С., Агеев В.В., Лебедева H.A. Поляризуемость горных пород и феномен высокоразрешающей электроразведки // Геофизика, № 4, 1996.

49. Соколов Н.П. Геохимические методы поисков месторождений полезных ископаемых. М.: Недра, 1985. 294 с.

50. Технология нефтегазопоисковых и геохимических исследований в системе «воздух-земля-скважина» / Под ред. Зубайраева C.JL, Стадника E.B. М.: Изд. ВНИИгеоинформсистем, 1990.

51. Уэйт Дж. Геоэлектромагнетизм. М.: Недра, 1987. - 235 с.

52. Фридрихсберг Д.А., Сидорова М.П. Исследование связи явления вызванной поляризации с электрокинетическими свойствами капиллярных систем. Вестник ЛГУ, сер. Физики и Химии, 1961, № 4.

53. Чжоу Цзыюн. Струйные ореолы рассеяния нефтегазовых месторождений в неоднородных горных породах и их изучение геоэлектрохимическими методами. С.- Петербург, СПГГИ. Диссертация на соиск. уч. ст. к.г.-м.н, 2000.-146с.

54. Шейнман С. М. Современные физические основы теории электроразведки. М., Недра, 1968. 221 с.

55. Шейнман С.М. Об импедансе поверхности электроннопроводящих минералов. В кн. : Вопросы электроразведки рудных месторождений. Л., 1977.

56. Шейнман С.М., Исаев Г.А., Полетаева Н.Г. О влиянии вызванной поляризации на результаты измерений в методе переходных процессов. Сб. трудов НПО "Геофизика". Методы разведочной геофизики, Л., 1982, вып. 57, с. 38 47.

57. Шлюмберже К. О вызванной поляризации в скважинах. Каротаж,различные операции в скважинах, 1934, вып. 6.

58. Яковлев А.П., Круглова З.Д. Изменение пород под влиянием залежей нефти и газа и возможность их выявления геофизическими методами. -М., 1977.

59. Cole K.S., Cole R.H. Dispersion and absorbtion in dielecrtrics. J. Chem. Phys., 1941, v.6.

60. Davidenko Y., Ivanov S., Kudryavceva E., Legeydo P., Veeken P.C.H. Geoelectric surveying, a useful tool for hydrocarbon exploration // 70-th EAGE conference & exhibition incorporating SPE Evropec 2008/ Rome, Italy (9-12 june.2008)/ - P.53

61. Induced polarization (IP) method in oil exploration the cause of the IP anomaly / Zhang, Saizhen, Nie, Xinwu, Li, Yinxian, Anchang, Zhou Jiping and Yang , Guanding. - 1988.

62. Pelton W.H., Ward S.H., Hallof P.G., Still W.R., Nelson P.N. Mineral discrimination and removal inductive coupling with multifrequency IP. Geophysics, 1978, vol.43, No 3.

63. Pirson S. D. Progress in magnetoelectric exploration. -Oil and Gas J., 1982, v. 80, . 41.

64. Schumacher D. Hydrocarbon-induced alteration of soils and sediments, in D. Schumacher and M. A. Abrams, eds., Hydrocarbon migration and its near-surface expression: AAPG Memoir 66. 1996. - P. 71-89.

65. Siegel H. Induced polarization and its role in mineral exploration. Canadian Min. Metall. Bull., 1962, v. 5, 1 600.

66. Snaider D. Exploration for petroleum using complex resistivity measure-ments.-Advances in induced polarization and complex resistivity. The University of Arizona, 1984, January, 5-7.

67. Sogade J.A., Scira-Scappuzzo F., Vichabian Y., Shi W., Rodi W., Lesmes D.P., Morgan F.D. Induced-polarization detection and mapping of contaminantplumes. Geophysics, vol.71, No 3 (may-june 2006). P. B75-B84

68. Sternberg B.K. A review of some experience with the induced-polarization / resistivity method for hydrocarbon surveys: Successes and limitation. Geo-phisics, 1991, 56, 10. P. 1522-1532.

69. Tarnsial an E.M. method applied to hydrocarbon exploration. Rocroi J.P., Farabolini W., Andriux P. 1989.

70. Veeken P.C.H., Legeydo P.J., Davidenko Y.A., Kudryavceva E.O., Iva-nov S.A., Chuvaev A. Benefit of the induced polarization geoelectric method to hydrocarbon exploration // Geophisics, vol.74, N.2 (march-april.2009). P. 1-XXXX-12-XXXX.

71. Veeken P., Legeydo P., Pesterev I., Davidenko Y., Kudryavceva E., Ivanov S. Geoelectric modelling with separation between electromagnetic and induced polarization field components // First Break, issue 12, vol.27 (decem-ber.2009).

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

Кудрявцева, Елена Олеговна
кандидата геолого-минералогических наук
Санкт-Петербург, 2010
ВАК 25.00.10

Диссертация

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

Автореферат

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