close

Вход

Забыли?

вход по аккаунту

Договор найма нежилого помещения образец E;pdf

код для вставкиСкачать
ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО
НАУЧНО-ПРОИЗВОДСТВЕННОЕ ПРЕДПРИЯТИЕ «ВНИИГИС»
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«УРАЛЬСКИЙ ГОСУДАРСТВЕННЫЙ ГОРНЫЙ УНИВЕРСИТЕТ»
На правах рукописи
ЯХИНА ИРИНА АЙРАТОВНА
РАЗВИТИЕ ТЕОРИИ И МЕТОДИК ИНТЕРПРЕТАЦИИ В
ЗОНДИРОВАНИЯХ МЕТОДОМ ПЕРЕХОДНЫХ ПРОЦЕССОВ ПРИ
ИЗУЧЕНИИ ГЕОЭЛЕКТРИКИ СЛАБОКОНТРАСТНЫХ СРЕД
Специальность: 25.00.10 – «Геофизика, геофизические методы поисков
полезных ископаемых»
Диссертация на соискание ученой степени
кандидата геолого-минералогических наук
Научный руководитель:
доктор геолого-минералогических наук,
профессор Виноградов Альберт Михайлович
Екатеринбург, 2014
2
ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ…………………………………………………………………….........4
Глава 1 ОБЗОР ВЫЯВЛЕННЫХ ТИПОВ НЕМОНОТОННЫХ
ПЕРЕХОДНЫХ ПРОЦЕССОВ СОВМЕЩЕННЫМИ УСТАНОВКАМИ
В ПРАКТИКЕ РЕШЕНИЯ ПОИСКОВЫХ ЗАДАЧ.........................................12
1.1. Особая роль совмещенной установки при обнаружении аномальных
переходных процессов ………………………………………………..13
1.2. Немонотонность переходных процессов как следствие сложной
зависимости удельной электропроводности среды от временичастоты наведенного электромагнитного поля……………………...16
1.3. Новые регистрируемые типы немонотонных переходных процессов
совмещенными установками…………………………………….........24
Выводы…………………………………………………………….........29
Глава 2 МЕТОДИЧЕСКИЕ ОСНОВЫ ИНТЕРПРЕТАЦИИ
ПЕРЕХОДНЫХ ПРОЦЕССОВ В ЗМПП...............…………….........................30
2.1. Динамические параметры сигнала в ЗМПП.........................................34
2.2. Методика интерпретации данных ЗМПП способом «плавающей»
плоскости с помощью динамических параметров сигнала...............39
2.3. Понятие глубинности. Расчет кажущейся глубины с помощью
динамических параметров сигнала......................................................42
Выводы……………………….………………………….……...............57
Глава 3 ЧИСЛЕННЫЕ СПОСОБЫ ВЫДЕЛЕНИЯ НЕМОНОТОННЫХ
ПРОЦЕССОВ И МЕТОДЫ ПОВЫШЕНИЯ ИХ КОНТРАСТНОСТИ
ОСОБЫМИ ПРИЕМАМИ ИХ РЕГИСТРАЦИИ..............................................58
3.1. Многоразмерные зондирования ЗМПП как основа повышения
разрешающей способности электромагнитной разведки....................60
3.1.1. Обоснование выбора оптимального количества и
размеров установок в ЗМПП для достижения
необходимой глубинности.................................................63
3
3.1.1 Геологические и геоэлектрические особенности разреза
Ерыклинского участка Куакбашской площади
Ромашкинского месторождения и условия
применения ЗМПП………………………………………..65
3.1.2 Технология проведения многоразмерных зондирований
ЗМПП на Ерыклинском участке........................................69
3.2 Численные оценки влияния вызванной поляризации на
индукционные переходные процессы в ЗМПП......................................74
3.2.1 Особенности проявления вызванной поляризации при
рудопоисковых работах на примере Ничатской
площади................................................................................74
3.2.2 Особенности проявления вызванной поляризации при
нефтепоисковых работах на примере Ерыклинского
участка..................................................................................83
3.2.3 Выделение индукционной и поляризационной
составляющих переходного процесса, используя
усовершенствованные алгоритмы приведения…………92
3.2.4 Выделение поляризационных интервалов по
немонотонности спада первой и второй производных
сигналов ............................................................................101
Выводы…………………………………………………………….........104
ЗАКЛЮЧЕНИЕ.....................................................................................................105
ПЕРЕЧЕНЬ УСЛОВНЫХ ОБОЗНАЧЕНИЙ………………………………...107
СПИСОК ЛИТЕРАТУРЫ...................................................................................108
4
ВВЕДЕНИЕ
Актуальность
В настоящее время популярность электроразведочных методов растет. С
новой усовершенствованной аппаратурой,
эффективными
теоретической базой и новыми
технологиями измерений, электромагнитные исследования
могут конкурировать по эффективности с сейсморазведкой, выходя, таким
образом, на новый уровень в ряду высоких технологий. На сегодняшний день
электроразведка способна решать задачи различного плана: картировочные,
поисковые (руда, нефть, вода, битум и др.), экологические, инженерные. В ряду
множества методов электроразведки можно выделить
метод зондирований
становлением поля (ЗС), разработанный А.Н. Тихоновым, С.М. Шейнманном
и
усовершенствованную
технологию
эффективного
поиска
–
метод
зондирований становлением поля в ближней зоне (ЗСБ), авторами которого
являются В.А. Сидоров и В.В. Тикшаев (теоретическая основа разработана
А.А. Кауфманом, Г.М. Морозовой). Основная задача метода – исследование
процессов становления электромагнитного поля в среде и выявление основных
закономерностей их изменения, связанных с ее геоэлектрикой.
Сейчас речь все чаще идет о высокоразрешающих возможностях
электроразведки. Благодаря использованию при работах ЗСБ и МПП
совмещенной установки удается повысить детальность исследований. Высокой
разрешенности данных методов способствует наличие индукционно вызванной
поляризации (ВПИ) или низкочастотной дисперсии (НЧД) в полифазных
породах, особенно с углеводородным насыщением. Благодаря наличию ВПИ
удается
выделять
интерпретации
малоконтрастные
из-за
малости
их
слои,
вклада
ранее
в
усредняемые
суммарную
при
продольную
проводимость.
При интерпретации
методом «кажущейся продольной проводимости»
[18, 56] в варианте «плавающей плоскости» [49] возможно выявить «тонкие»
5
особенности
строения разреза. Но расчет удельных сопротивлений
(проводимостей) через величину продольной проводимости представляет на
практике неоднозначную задачу, и поэтому возникают неоправданные, по
отношению к реальным возможностям метода, погрешности в результатах
интерпретации.
В этой связи, поиск способов исключения данных
погрешностей для возможности детального изучения разреза с выделением
слабоконтрастных слоев является актуальной проблемой. Способ «плавающей
плоскости» позволяет решать как структурные задачи в нефтяной геологии,
гидрогеологии, так и поисковые в нефтяной и рудной геофизике. Этот подход
позволил одинаково решать задачи ЗСБ и МПП, что и привело к созданию
одной аббревиатуры для методов – ЗМПП.
Цель
исследований
–
развитие
и
усовершенствование
методики
интерпретации данных ЗМПП с выделением малоконтрастных особенностей
строения разреза в совмещенных установках с использованием «плавающих
плоскостей».
Научные задачи исследования
1. Создать обобщенный алгоритм интерпретации данных ЗМПП (над
электропроводящими
полупространством,
разрезом над изолятором)
плоскостью,
слоистым
способом «плавающей плоскости» с
помощью введения динамических параметров сигнала.
2. Показать возможность
применения созданного интерпретационного
алгоритма на моделях сред с различной сложностью, включая
выявление тонких (малоконтрастных) особенностей строения разреза
по результатам ЗМПП, тем самым, повысить эффективность метода при
решении задач геоэлектрики. Реализовать разработки на практическом
материале.
6
Защищаемые научные положения
1. С использованием модели «плавающей» проводящей плоскости и
динамических параметров сигнала
определена функция глубинности,
которая позволяет оценить глубину проникновения вихревых токов в
проводящей геологической среде под совмещенной установкой
в
каждый момент времени.
2. Предложенная методика обработки и интерпретации данных ЗМПП с
совмещенными установками с использованием динамических параметров
сигнала позволяет
повысить точность оценки геоэлектрических
параметров разреза за счет выделения слабоконтрастных по удельной
электропроводности слоев.
3. Разработанный
способ
приведения
результатов
многоразмерных
зондирований к одному размеру контура совмещенной установки
расширяет временной диапазон при выявлении поляризационных
эффектов над плоскослоистым нефтеносным разрезом.
Научная новизна
1. Разработаны и программно реализованы оригинальные алгоритмы
интерпретации
данных
ЗМПП
методом
«кажущейся
продольной
проводимости» с использованием динамических параметров сигнала.
2. Определена функция глубинности для моделей однородного проводящего
полупространства, проводящей плоскости и слоистого разреза над
изолятором
через
динамические
параметры
сигнала
с
помощью
«плавающей плоскости».
3. Показана практическая возможность выделения поляризационной и
индукционной частей процесса становления поля на всем временном
интервале измерений ЭДС по немонотонности спада первых и вторых
производных сигнала.
4. Разработан
алгоритм
приведения
результатов
многоразмерных
зондирований к одному размеру контура совмещенной установки
по
предложенным аналитическим формулам и на практическом материале
7
показана его эффективность при выделении поляризационных процессов
в регистрируемом сигнале во всем временном диапазоне измерений.
Практическая значимость результатов и реализация работы
Созданные в работе методики и алгоритмы интерпретации индукционных
переходных процессов в совмещенных квадратных установках с учетом ВПИ
позволяют выделять малоконтрастные по удельной электропроводности слои
или локальные неоднородности в разрезе. Тем самым, существенно повышается
разрешающая способность метода ЗМПП при решении поисковых и
структурных задач.
Описанная технология полевых работ (применительно к нефтепоисковым
работам), средства обработки и интерпретации результатов по выделению
поляризационных эффектов позволяют осуществлять их глубинную привязку к
конкретному стратиграфическому горизонту.
Реализованная в работе методика интерпретации данных ЗМПП была
применена на Ромашкинском нефтяном месторождении (республика Татарстан)
при решении следующих геологических задач: проведение межскважинной
корреляции с целью оконтуривания залежей углеводородов, прогнозная оценка
характера нефтенасыщения.
Исследование в области развития теоретических основ интерпретации
данных
импульсной индуктивной электроразведки при решении структурно-
поисковых задач было поддержано грантом для молодых ученых УрО РАН в
2011 году.
Разработанный способ геоэлектроразведки, заключающийся в возможности
выявления поляризационных эффектов
в микро-миллисекундном диапазоне
времени, соответствующем электродинамическим процессам над нефтяными
залежами,
защищен патентом на изобретение.
Новые методы обработки и алгоритмы интерпретации опробованы на
полевых материалах ЗМПП и геологических данных по ОАО НПП
«ВНИИГИС»
по
республике
Татарстан
(совместно
с
ООО
ТНГ
«Казаньгеофизика», ОАО «Татнефтегеофизика», НГДУ «Лениногорскнефть»,
8
ОАО «Татнефть»), по Иркутской области, Чукотскому автономному округу и
другим регионам.
Личный вклад автора
Теоретические
и
практические
результаты,
приведенные
в
работе,
получены лично или при непосредственном участии автора.
Основной личный вклад
заключается
в
исследовании и разработке
оригинальных методик и алгоритмов интерпретации зондирований методом
переходных процессов, в т. ч. применительно к решению нефтепоисковых
задач. Лично автором были получены следующие результаты:
– исследованы, введены в методику интерпретации и функционально
определены динамические параметры сигнала и их производные;
– произведены численные расчеты экстремальных точек динамических
параметров;
–
разработана методика интерпретации способом «плавающей плоскости» с
учетом динамических параметров;
– реализовано использование динамических параметров в алгоритмах
интерпретации для дифференциации малоконтрастных сред при решении
обратных задач;
–
получены
аналитические
формулы
приведения
сигналов
при
многоразмерных зондированиях;
– показана возможность выделения поляризационных процессов с помощью
аналитических формул приведения;
–
проведена апробация созданных оригинальных алгоритмов и методик
интерпретации на экспериментальных данных (на базовых геоэлектрических
моделях) и на практическом материале.
В основу диссертации положены материалы, полученные по результатам
опытно-экспериментальных работ, проводившихся отделом электромагнитных
исследований
ОАО
НПП
«ВНИИГИС»
на
Ерыклинском
участке
Ромашкинского месторождения нефти в республике Татарстан. Автор работы
лично участвовал в непосредственном сборе и анализе материалов полевых
9
измерений,
обработке
результатов,
геолого-геофизической
информации,
получении
которые легли в основу диссертации, защите отчета по
результатам выполненных работ.
В совместных публикациях [70, 71, 74, 76, 78, 79] автору принадлежит
производство
численных
расчетов,
формулировка
основных
подготовка материалов к публикации; в совместных работах
выводов,
[41, 77] автор
непосредственно участвовал в сборе и анализе геолого-геофизических данных,
построении структурных карт и разрезов; в работах [75, 83, 85] личное участие
автора выражается в непосредственной постановке задачи, производстве
численных расчетов, сборе и анализе практического материала.
Апробация работы
Основные научные положения работы и результаты
докладывались и
обсуждались на международных конференциях и семинарах, таких как:
Международная научно-техническая конференция «Проблемы нефтегазового
дела» (Уфа, 2006);
IX, XI, XIII Уральская молодежная научная школа по
геофизике «Современные проблемы геофизики» (Екатеринбург, 2008, 2010,
2012
гг.);
Международная
научно-практическая
конференция
«Ядерно-
геофизические полевые, скважинные и аналитические методы при решении
задач поиска, разведки и разработки месторождений твердых полезных
ископаемых»
(Октябрьский, 2009 г.); Международная научно-практическая
конференция «Казанская геологическая школа и её роль в развитии
геологической науки в России» (Казань, 2009 г.); 37-я, 40-я сессии
Международного семинара им. Д.Г. Успенского «Вопросы теории и практики
геологической интерпретации гравитационных, магнитных и электрических
полей» (Москва, 2010, 2013 гг.); EAGE, Первая международная научнопрактическая конференция «ГЕОБАЙКАЛ-2010» (Иркутск, 2010 г.); VI
Научные чтения им. Ю.П. Булашевича «Глубинное строение, геодинамика,
тепловое поле Земли, интерпретация геофизических полей» (Екатеринбург,
2011 г.); Международная научно-практическая конференция «Аппаратурнометодические комплексы и технологии ГИС и ядерно-геофизические методы
10
для исследования нефтегазовых и рудных скважин» (Октябрьский, 2012 г.);
Международная научно-практической конференция «Высоковязкие нефти и
природные
битумы:
проблемы
повышения
эффективности
разведки
и
разработки месторождений» (Казань, 2012 г.); III молодежная школа-семинар
«Современная тектонофизика. Методы и результаты» (Москва, 2013 г.).
Публикации
По теме диссертации опубликовано 19 статей, в том числе 2 статьи – в
ведущих научных рецензируемых журналах из перечня ВАК, получен патент
РФ №2494419.
Объем и структура работы
Работа состоит из 3-х глав, введения, заключения, списка литературы.
Объем работы составляет 119 страниц, в том числе 41 рисунок.
Во введении описаны основные положения диссертационной работы
(актуальность, цель, научные задачи, теоретическая и практическая значимость
результатов), сформулированы защищаемые положения и научная новизна.
В первой главе приведен обзор выявленных типов аномальных переходных
процессов в осесовмещенных установках. Природа появления таких процессов
объяснена с позиции влияния индукционно вызванной поляризации и
необходимости учета низкочастотной дисперсии удельной электропроводности
горных пород при решении задач геоэлектрики. Приведены примеры
регистрации сложных форм переходных процессов в практике решения
поисковых и структурных задач.
Во второй главе описаны основные процедуры обработки и интерпретации
данных ЗМПП. Усовершенствованы методика и алгоритмы интерпретации
переходных процессов способом «плавающей плоскости»
с помощью
динамических параметров сигнала. Показана реализация таких уточненных
алгоритмов на базовых геоэлектрических моделях.
В третьей главе исследуются методы повышения контрастности переходных
процессов при выделении «тонких» особенностей строения разреза, которые
отражаются в нарушении
спада ЭДС. Приводятся способы изучения и
11
выявления скрытого влияния вызванной поляризации на индукционные
переходные процессы, связанные с немонотонностью спада переходного
процесса и их численные оценки. На практическом материале (на примере
Ромашкинского нефтяного месторождения в республике Татарстан) показана
реализация методического комплекса обработки данных ЗМПП и алгоритмов
приведения результатов многоразмерных зондирований к одному размеру
контура совмещенной установки.
Благодарности
Автор благодарен своему научному руководителю
д.г.-м.н. А.М.
Виноградову за поддержку автора в его научной позиции и неоценимую
помощь в подборе фактического материала при написании данной работы.
Автор
выражает
благодарность
сотрудникам
ИГФ
УрО
РАН
за
всесторонний анализ работы и обсуждения, которые помогли улучшить
содержательную часть диссертации.
Отдельную благодарность автор выражает сотрудникам ИНГГ СО РАН –
д.т.н. В.С. Могилатову, д.г.-м.н. Н.О. Кожевникову и д.ф.-м.н. Е.Ю. Антонову,
поддержавшим данную работу и сделавшим важные и ценные замечания.
Автор благодарит д.г.-м.н. А.А. Редозубова, д.ф.-м.н. В.П. Губатенко, д.т.н.
Л.Е. Кнеллера, д.т.н. Ю.А. Гуторова за помощь и полезные консультации при
обсуждении отдельных моментов диссертации.
За
предоставление
результатов
автор
возможности
выражает
практической
благодарность
реализации
сотрудникам
научных
лаборатории
электромагнитных исследований ОАО НПП «ВНИИГИС».
Автор глубоко признателен своему отцу – к.т.н. А.М. Яхину , оказавшему
большое влияние на формирование научных взглядов соискателя, за
постоянное внимание и поддержку, а также неоценимую помощь в написании
работы.
12
Глава 1 ОБЗОР ВЫЯВЛЕННЫХ ТИПОВ НЕМОНОТОННЫХ
ПЕРЕХОДНЫХ ПРОЦЕССОВ СОВМЕЩЕННЫМИ УСТАНОВКАМИ В
ПРАКТИКЕ РЕШЕНИЯ ПОИСКОВЫХ ЗАДАЧ
Целью данного раздела является обобщение имеющихся сведений об
аномальных переходных процессах, проявляющихся при исследованиях
методами электроразведки, в частности, зондированиями становлением поля в
ближней зоне (ЗСБ) и методом переходных процессов (МПП). В данной главе
также
приведены
практические
примеры
полученных
аномальных
характеристик (немонотонных, отрицательных переходных процессов, с
нарушением монотонности спада процесса
производных сигнала)
в
совмещенных установках и рассмотрена некоторая взаимосвязь с природой
этих явлений и со смежными областями знаний.
В работе подробно представлены теоретические расчеты и практические
материалы, касающиеся одного из типов немонотонных переходных процессов
– это зарегистрированные монотонно спадающие переходные процессы, но
имеющие нарушения в спаде (немонотонность) первой и второй производных.
Такие процессы в работе названы немонотонными, поскольку в них
присутствует скрытое влияние индукционно вызванной поляризации не
проявляющееся в явном виде.
Наиболее перспективной установкой для выявления аномальных сигналов
является совмещенная установка (типа «петля в петле», квадратная приемная
петля раскладывается вдоль генераторной) или однопетлевая (квадратная
генераторная петля выполняет также функцию приемного контура). Процесс
возбуждения поля в исследуемой среде вызывается включением-выключением
тока в генераторной петле. Это поле вызывает в приемной петле ЭДС
индукции.
Если выполнить измерения в точке при различных значениях
времени t, то можно построить зависимость ЭДС от времени. Такая
зависимость называется переходной характеристикой среды или сигналом
13
становления поля и ее график представляет кривую спада. Данный процесс
является монотонным, при этом каждое последующее зарегистрированное
значение ЭДС должно быть строго меньше предыдущего.
Физику процесса появления аномальных форм переходных процессов
можно объяснить следующим образом. При протекании через среду вихревых
токов среда поляризуется так же, как и при кондуктивном возбуждении поля. В
результате в среде возникает ток вызванной поляризации, конфигурация
которого аналогична, а направление противоположно вихревому току. Этот ток
создает свое магнитное поле, которое, складываясь с магнитным полем
вихревых токов, изменяет его величину, и следовательно, искажает результаты
измерений. При малых временах t поле токов вызванной поляризации обычно
мало по сравнению с полем токов индукционного происхождения, поэтому
влияние вызванной поляризации незаметно. Но вследствие того, что токи ВП
затухают намного медленнее токов индукции, при больших значениях времени
t их влияние становится более значительным и при высокой поляризуемости
может привести даже к появлению отрицательных значений импеданса или
форм сигнала с нарушением монотонности спада (немонотонных). Более
подробно о природе этого явления будет изложено ниже.
1.1. Особая роль совмещенной установки при обнаружении аномальных
переходных процессов
Развитие и расширение работ ЗСБ и МПП с совмещенными установками
при прямых поисках рудных объектов подтвердили существование и выявили
примеры
регистрации
рассматривались
как
отрицательных
технический
брак.
процессов,
В
1966
которые
году
при
ранее
поисках
медноколчеданных месторождений на Южном Урале А.Н. Родионов и А.М.
Виноградов в процессе производственных исследований совмещенными
петлями 300ˣ300м зарегистрировали переходные характеристики в виде
кривых с переходами через ноль в область отрицательных значений ЭДС и с
образованием экстремума в отрицательной области. По результатам этих
полевых исследований вышла первая журнальная публикация [43] об
14
аномальных переходных процессах. Авторы статьи, А.Н. Родионов и А.М.
Виноградов, связали отрицательные процессы с геологией участка и поставили
задачу изучения природы этих процессов. В работе сделано предположение,
что для участков с отрицательными процессами необходимо кроме чисто
индукционного становления учитывать емкостные процессы, т.е. при решении
уравнений Максвелла учитывать токи смещения. Аналогичные отрицательные
процессы
МПП
наблюдались
авторами
работы
[6]
на
колчеданном
месторождении.
Важную роль для изучения аномальных переходных характеристик сыграли
исследования, выполненные под руководством В.А. Сидорова и ИГФЭ
сотрудниками ВНИИГИС (А.М. Яхин, Б.В. Бучарский, В.А. Ключников, Ю.В.
Николаев, А.К. Ткаченко) в Якутии. В условиях моголетнемерзлых пород
Якутии, начиная с 1977 года, методом ЗСБ выявлены процессы еще более
сложных типов – с двойным изменением знака процесса и положительным
спадом в конце, а также, положительные процессы со сменой знака
производной сигнала [38, 39].
На
рис.
1.1
приведены
переходные
процессы
сложных
типов
в
совмещенных контурах квадратной формы со сторонами L=200м вблизи
кимберлитовых тел в Якутии. На кривой с индексом 3 пунктиром отмечена
отрицательная часть
процесса. Эта кривая имеет двойную смену знака. Для
кривой с индексом 2 характерно изменение знака производной. Для сравнения
приведена нормальная кривая без влияния поляризуемости среды. В
сопоставлении с этой кривой на кривых с индексами 1 и 2 заметно, что, начиная
с определенного времени, значения по кривым расходятся по уровню. Такое
расхождение может существенно влиять на результаты ЗСБ.
Практика
работ
ЗСБ
и
МПП
совмещенными
или
однопетлевыми
установками выявила надежно зарегистрированные переходные процессы,
имеющие переходы через нулевые значения и в последующем отрицательные
спады.
15
Рис. 1.1. Аномальные кривые ЗСБ в совмещенных контурах вблизи месторождений алмазов
в Якутии (L=200 м) (Сидоров В.А., 1985).
Р.Б. Журавлева и И.Э. Гаврилова показали, что процессы с однократной
сменой знака могут быть объяснены влиянием вызванной поляризации (ВП)
изучаемых объектов [14]. Была поставлена принципиальная задача учета
влияния поляризуемости горных пород на процесс становления поля. Позже
В.П. Губатенко и В.В. Тикшаевым [10] было установлено, что над любой
недиспергирующей средой переходные процессы в совмещенных установках
должны иметь строго монотонный спад. Изменение знака процесса может быть
вызвано только поляризуемостью среды. Таким образом, совмещенная
установка сыграла особую роль при выявлении отрицательных переходных
процессов.
16
1.2. Немонотонность переходных процессов как следствие сложной
зависимости удельной электропроводности среды от времени-частоты
наведенного электромагнитного поля
Суть методов ЗСБ
и МПП состоит в следующем: с помощью
незаземленного квадратного контура возбуждается первичный импульс
магнитного поля и в паузах между импульсами принимается изменяющееся во
времени магнитное поле или его производная на тот же (однопетлевая
установка) или другой контур (совмещенная установка). Выключение импульса
магнитного
поля
приводит
к
возникновению
в
проводящей
среде
устанавливающегося во времени ЭДС индукции, возбуждающей затухающие
вихревые токи. Скорость и интенсивность затухания этих токов зависит от
проводящих свойств среды. После выключения импульса электрического поля
в проводящих и поляризующихся породах
возникают вторичные поля
вызванной поляризации противоположного знака. Процессы, происходящие в
элементарном объеме поляризующейся среды при включении постоянного
электрического поля схожи с процессами, происходящими при зарядке
аккумулятора. Зарядка – это накопление электрической энергии при
прохождении тока за счет превращения ее в другой вид энергии (например,
химическую), разрядка – процесс прохождения тока через внешнюю цепь
противоположно току зарядки. Соответственно этой схеме – зарядке
аккумулятора эквивалентен процесс прохождения тока при постоянном
поляризующемся поле, разрядке – вторичный процесс,
возникающий при
выключении поляризующегося поля [72].
Переходные процессы в совмещенных установках достаточно хорошо
теоретически изучены [2, 18, 30, 46, 55].
неоднородностью,
учитывая
только
В среде с произвольной
активную
составляющую
электропроводности, спад ЭДС (убывание ЭДС до нуля при t   ) от времени
происходит монотонно [9, 10], ни сигнал, ни его производные по времени не
изменяют знака. В [15, 37] показано, что при учете вызванной поляризации
17
над однородным полупространством должен наблюдаться процесс со сменой
знака.
Регистрируемые немонотонные переходные процессы со сменой знака
многие исследователи связали с низкочастотной дисперсией (НЧД) удельной
электропроводности
или
вызванной
поляризацией
горных
пород
при
индукционном возбуждении (ВПИ). Для объяснения необычных форм
переходных процессов, полученных при полевых работах сотрудниками ОАО
НПП ВНИИГИС, была создана приближенная теория индукционного влияния
вызванной поляризации в произвольных средах с осевой симметрией [50, 53].
Эта теория позволяет объяснить форму всех основных типов аномальных
процессов [36] слоистостью разреза и особыми соотношениями временных
характеристик
вызванной
поляризации
(ВП)
и
чисто
индукционного
становления.
Основой для изучения НЧД в электромагнитных исследованиях является
феноменологическая теория ВП [26, 27, 29, и др.].
возможность описания процесса ВП
Она предусматривает
в уравнениях электродинамики
небольшим количеством параметров, не рассматривая при этом всей
сложности природы этого явления.
Природа ВП, НЧД горных пород разнообразна. Один из ее видов связан с
электрической
неоднородностью
(гетерогенностью)
изучаемой
среды.
Формирование ВП происходит не только за счет гетерогенности среды, но и
благодаря большой совокупности разнообразных электрофизических, в том
числе электромагнитных, процессов, приводящих к разделению зарядов под
действием приложенного электрического поля. Главную роль играют явления
ионной диффузии и электрокинетические явления, к которым относятся
электроосмос
и
фильтрация,
являющиеся
результатом
взаимодействия
приложенных к жидкой фазе электрических и механических сил. Углубленное
изучение поляризуемости неоднородных пород [7] показало, что наличие
углеводородного насыщения в них может приводить к их интенсивной
поляризации с постоянными спада, сравнимыми со временем регистрации
18
импедансов при ЗСБ. Такой класс поляризующихся объектов отнесен к
электрофизическим.
Одной
из
причин
электрофизической
дисперсии
(зависимость электропроводности от частоты или поляризуемость среды)
является гетерогенность среды, которая также характерна для пород,
содержащих тонкие включения высокоомных прослоев, в т.ч. пленок нефти.
Чередование
высокоомных
прослоев
с
низкоомными
приводит
к
возникновению емкостных эффектов [7].
При наведении нестационарных
внешних электромагнитных полей,
например, ступенчатым выключением тока в генераторной катушке, процессы
ВП, связанные с эпигенетическими изменениями пород электрохимической или
электростатической природы, имеют относительно длительную релаксацию по
сравнению
с
сопротивлениями
полями
становления,
слагающих
связанными
осадочный
разрез
с
электрическими
(электродинамические
процессы) горных пород. Включение в низкоомную среду несовершенных
диэлектрических прослоек, как показано В.А. Сидоровым [52], может резко
увеличить поляризуемость среды и время релаксации, а на низких частотах
эффективная (кажущаяся) диэлектрическая проницаемость может достигать
огромной величины, что и приводит к появлению аномальных осциллирующих
форм переходных характеристик [36, 52]. Подобные явления в несовершенных
слоистых конденсаторах (но при наличии внешней цепи для токов) изучены
еще в XIX веке и известны как эффект Максвелла-Вагнера. В работе [11] этот
эффект проанализирован для периодической слоистой среды (с периодом
d=d1+d2), состоящей из чередующихся слоев мощностью d1 и d2
с
проводимостями σ1, σ2 и абсолютными диэлектрическими проницаемостями ε1
и ε2. Путем усреднения электромагнитного поля по физически малому объему
среды показано, что тангенциальные проводимость и диэлектрическая
проницаемость не зависят от частоты, а нормальные – зависят, так что
комплексная проводимость по нормалям к слоистости

ˆ    

 0 
  

  i      2 02  ,
2 2
1   
1  

(1.1)
19
где эффективная проводимость слоистой среды на высоких частотах
 
1 22   2 12
,
 2  1 2
а эффективная проводимость на низких частотах
0 
 1 2
,
 2   1
(1.2)
эффективная диэлектрическая проницаемость на высоких частотах
 
 1 2
,
 2  1
эффективная диэлектрическая проницаемость на низких частотах
0 
1 22   2 12
,
 2   1 2
время релаксации
 

где
Приведенные
выражения
 2   1
,
 2   1
d1
d
,   2  1 .
d
d
для
комплексной
электропроводности,
соответствующие эмпирическим законам временного спада ВП, не являются
универсальными, но отражают основные закономерности ВП в природных
геологических средах.
Большой вклад в развитие и изучение теории влияния поляризуемости на
индукционные переходные процессы внесли Г.В. Астраханцев, И.З. Гаврилова,
В.П. Губатенко, Р.Б. Журавлева, Г.А. Исаев, В.В. Кормильцев, Ф.М.
Каменецкий, С.С. Крылов, А.В. Куликов, А.Н. Мезенцев, Н.Г. Полетаева, А.А.
Рыжов, В.А. Сидоров, В.М. Тимофеев, С.М. Шейнманн, Е.А. Шемякин, А.М.
Яхин, T. Lee, H.F. Morrison, P. Weidelt, G.F. West.
Изучению явления НЧД в электромагнитных исследованиях посвящены
также работы В.В. Агеева, Е.Ю. Антонова, И.А. Безрука, Е.С. Киселева, В.А.
Комарова, Н.О. Кожевникова, П.Ю. Легейдо, В.В. Сочельникова, Б.С. Светова,
А.Ф. Постельникова, В.В. Филатова, J.R. Wait, J.W. Hohman, D.J. Marshall, G.
20
Klein, G.A. Neman и др. Исследование феномена высокого и сверхвысокого
разрешения ЗСБ с учетом НЧД выполнили В.В. Агеев, В.В. Задорожная, А.Г.
Небрат, А.В. Сафонов, Б.С. Светов, В.В. Тикшаев и др.
Работа [14] позволила качественно и, в некоторых случаях количественно,
объяснить возможность появления процессов, имеющих отрицательные спады
в поздней стадии для модели однородно проводящего и поляризующегося
полупространства.
Изучению
этой
модели
были
посвящены
позже
опубликованные работы советских [16, 27, 29, 37, 61] и зарубежных авторов
[93].
Важными для практики являются также другие простые модели. В работе
[92] впервые рассматривается возможность смены знака переходного процесса
в однопетлевом варианте метода переходных процессов для
проводящего и поляризующегося шара. В [61]
однородно
рассмотрена модель –
проводящая и поляризующаяся плоскость, показана возможность появления
отрицательного спада в поздней стадии переходного процесса. Автору работы
[92] удалось получить переходную характеристику с отрицательным спадом
для модели проводящей и поляризующейся петли, при этом, проводимость
описывалась с помощью простой модели релаксации известной под названием
Cole-Cole (1941) в виде:

 0 1  i c
 (i ) 
c
1   i 

(1.3)
где  (i ) – комплексная удельная электропроводность; 
– постоянная
времени, определяющая продолжительность явления вызванной поляризации;
c
– параметр, характеризующий затухание (часто находящийся в диапазоне 0.1-
0.6);  0 – удельная электропроводность среды на самых низких частотах, т. е.
когда среда – чисто проводящая,  0    (1   ) ;   1   , где  – стационарная
поляризуемость. Эта модель описывает широкий круг явлений ВП и получила
широкое распространение. Если рассматривать относительно низкочастотную
область этого явления, ее физическому смыслу больше соответствует его
описание
в
терминах
частотно-зависимой
электропроводности
(хотя
21
математически можно получить его описание через частотно-зависимую
диэлектрическую проницаемость).
На
основе
(1.3)
или
других
подобных
формул
вещественную
электропроводность в решении прямой задачи для той или иной модели разреза
можно заменить комплексной, применить к этому результату обратное
преобразование Фурье-Лапласа и получить таким образом решение для t области с учетом ВП или НЧД. Такой способ широко используется. С его
помощью
различными
исследователями
выполнено
математическое
моделирование влияния ВП на индукционные переходные процессы, в том
числе для установки с совмещенными петлями.
В большинстве случаев
рассматривались однородное полупространство и горизонтально-слоистая
среда, а также некоторые локальные объекты.
Результаты математического моделирования влияния ВП на индукционные
переходные процессы рассмотрены в работах А.Н. Мезенцева [31 и др.].
Влияние вызванной поляризации на частотные и переходные характеристики
исследуется в работе [28], принципиальной основой которой является учет
вызванной поляризации в уравнениях электродинамики, т. е. теория Максвелла
рассмотрена с учетом этих особенностей.
На основе описания процесса ВП с частотной дисперсией по Cole-Cole
выполнены исследования по созданию алгоритмов и программ решения прямой
задачи импульсных электромагнитных зондирований и исследования в области
разработки программ инверсии данных
индуктивных и гальванических
зондирований для диспергирующих сред [3, 4 и др.].
Рассмотренные выше модели поляризующихся сред и форма комплексной
удельной электропроводности описывают лишь небольшой объем реально
встречающихся разрезов. Решения, получающиеся при следовании методике
оценок по виду  (i ) , не позволяют отразить в формулах все влияющие
параметры и, тем самым, существенно совершенствовать методику полевых
работ. В практике работ легче получить полевыми измерениями функцию  (t )
22
[45] по регистрируемым переходным характеристикам эквивалентную  (i )
при решении задач. Поэтому использование  (t ) в расчетах предпочтительнее.
Во временной области  (t ) может быть представлена в виде [22, 72]:

 t    0 1 



 F t  ,
1 

(1.4)
где F (t ) – временная характеристика спада ВП после ступенчатого выключения
установившегося значения поля; t – время после ступенчатого включения
единичного поляризующего поля;
F 0  1,
F t    0 ;
0
– удельная
электропроводность поляризующейся среды, определяемая в установившемся
режиме по (1.2).
Применение базовых зависимостей удельной электропроводности от
частоты (теоретически известных как модель Cole-Cole, экспоненциальных и
др.) не исчерпывает всех возможностей проявления ВПИ. Например,
экспериментально определенные и опробованные формулы спада ВП А.А.
Рыжова и В.А. Комарова во временном режиме не являются простым
следствием вышеуказанных математических моделей. Решение прямых задач
на основе базовых зависимостей
Cole-Cole приводит к появлению
дополнительных постоянных, которые могут быть установлены только путем
подбора (   , c …), т. е. расширяется область неоднозначности при решении
обратных задач, одни и те же
сигналы могут соответствовать совершенно
разным моделям среды. Например, возможна интерпретация сигналов от
слоистых разрезов с влиянием ВПИ в определенном временном режиме
измерений,
как
сигналов
от
локальных
объектов
из-за
кажущегося
затягивающегося спада. В этой связи, в данной работе сделан основной акцент
на использование переходных характеристик ВП во временном режиме,
которые можно получить непосредственно измерениями, чтобы в дальнейшем
двойным дифференцированием определять переходную характеристику ВПИ.
Ввиду сказанного, представляется наиболее важным выделение, в первую
очередь, сигналов ВПИ для реальных квадратных совмещенных контуров и
23
дальнейшая их интерпретация совместно с электродинамическими процессами
без ВП.
Переходные характеристики ВП ( F (t ) ) для поляризующихся сред получают
непосредственным
измерением
пропускания
единичного
определенной
длительности
напряжений
поляризующегося
t3
(время
вторичного
поля
прямоугольного
зарядки).
Но
после
импульса
возможно
также
использование теоретически рассчитанных F (t ) для конкретных моделей
среды. На практике применимы различные аппроксимационные формулы,
например, аппроксимационная формула В.А. Комарова [26]:
F (t , g , h) 

t  g 
ln 
,
(1   ) ln g h  t  h 
где g, h – временные параметры среды, причем g  1 1000 с, а h  103  1 с.
Существуют и другие формулы, хорошо описывающие временную зависимость
ВП. Например, аппроксимационная формула А.А. Рыжова [44]:
F (t ) 
  exp(Q 2t )  erfc (Q t )
1 
где Q – временной параметр среды; erfc – дополнение интеграла вероятности.
Зависящая
от
времени
электропроводность
элементарного
объема
поляризующейся среды, рассматривается как переходная реакция σ(t) на

воздействие единичного поляризующегося поля. Плотность тока j (t )


в этом

объеме для постоянного воздействия E определяется по формуле j (t )   (t )  E

(частный случай) или в общем случае j (t ) 
t
d 
E ( )  (t   )d [72]. Зависимость
dt 0
электропроводности от времени может быть представлена через временную
функцию ВП в виде (1.4). Такое описание процессов ВП аналогично описанию
процесса релаксации зарядов в несовершенном диэлектрике.
24
1.3. Новые регистрируемые типы немонотонных переходных
процессов совмещенными установками
На рис. 1.2
представлен один из зарегистрированных сложных типов
немонотонных переходных процессов – «отрицательный» переходный процесс
(Сабановский участок). Работы МПП проводились ИГФ УрО РАН (г.
Екатеринбург) в Челябинской области в 2006 году с целью геофизического
обеспечения поисков медно-колчеданных руд в Александринском рудном
районе.
На
рис.
1.3
представлен
немонотонный
зарегистрированный над осадочным карбонатным
переходный
процесс,
разрезом с включением
гипса по данным ЗСБ. Данные работы проводились в 2009 году ОАО НПП
«ВНИИГИС» (г. Октябрьский, республика Башкортостан) в Кунгурском районе
Пермского края на предмет поиска и разведки месторождений горных пород,
применяющихся при производстве строительных материалов. Отрицательные
процессы отсутствуют, однако процесс имеет явные нарушения спада сигнала
на поздних временах.
На рис. 1.4 приведен пример регистрации методом ЗМПП переходного
процесса и его производные. Измерения проводились в Республике Татарстан
ОАО НПП «ВНИИГИС» (г. Октябрьский, республика Башкортостан) на
Ромашкинском месторожении нефти (Ерыклинский участок) [64, 84]. Явное
проявление ВПИ в форме перехода к отрицательным процессам отсутствует,
однако наблюдается нарушение монотонности в спаде производных, что
указывает на влияние поляризуемости отдельных интервалов разреза [82].
На рис. 1.5 приведены результаты регистрации аномальных отрицательных
сигналов в терригенных разрезах с естественным углеводородным (битумным)
насыщением. Работы проводились ОАО НПП «ВНИИГИС» (г. Октябрьский,
республика Башкортостан) в Казахстане на участке «Иманкара». Повышенный
уровень абсолютных значений отрицательных сигналов отмечается в разрезах
скважин
с
интервалов.
повышенной
суммарной
мощностью
битумонасыщенных
25
4.70E-04
3.70E-04
Z(t), B/A
2.70E-04
1.70E-04
7.00E-05
-3.00E-05
-1.30E-04
1000
10000
100000
t, мкс
Рис. 1.2. Пример регистрации немонотонного «отрицательного» переходного процесса над рудным
разрезом, Челябинская область, Сабановский участок (L=600 м) (Яхина И.А., 2012).
26
1.00E+02
1.00E+01
1.00E+00
Z(t), В/А
1.00E-01
1.00E-02
1.00E-03
1.00E-04
1.00E-05
1
10
100
1000
10000
t, мкс
Рис. 1.3. Пример регистрации немонотонного переходного процесса над осадочным карбонатным
разрезом, Пермский край (L=50 м) (Яхина И.А., 2012).
27
100000
10000
1000
100
10
1
Z(t), В/А
|Z'(t)|, (B/A)/c
Z''(t), (B/A)/c2
0.1
0.01
0.001
0.0001
0.00001
0.000001
0.0000001
1E-08
10
100
1000
10000
100000
t, мкс
1
2
3
Рис. 1.4. Пример регистрации переходного процесса со скрытым влиянием ВПИ над нефтяным разрезом
и его производные, республика Татарстан (L=1000 м) (Яхина И.А., 2008): 1 –первая производная сигнала,
Z'(t); 2 –вторая производная сигнала, Z''(t); 3 – зарегистрированный сигнал Z(t).
28
Рис. 1.5. Пример регистрации аномальных (отрицательных) переходных процессов в терригенных
разрезах с естественным углеводородным (битумным) насыщением, Казахстан
(Яхин А.М., Яхина И.А., 2012, 2013).
29
Выводы
Совмещенными установками
при работах методом ЗСБ и МПП были
надежно зарегистрированы аномальные переходные процессы следующих
типов: отрицательные (спад в отрицательной области, спад в положительной
области), осциллирующие, двойные переходы через ноль, немонотонный спад в
положительной области, монотонный спад с изменением знака производной.
Все эти процессы, неинтерпретируемые в рамках классической теории
электродинамики, связаны с индукционным влиянием вызванной поляризации
(ВПИ) или поляризуемостью
является
следствием
зависимости
удельной
среды. Немонотонный спад ЭДС от времени
явления
низкочастотной
электропроводности
от
дисперсии
частоты
(НЧД),
т.е.
наведенного
электромагнитного поля.
Сложившийся способ изучения влияния поляризуемости среды на ИПП
путем замены удельной электропроводности среды в частотной области на
комплексную частотно-зависимую электропроводность позволил объяснить так
называемые «отрицательные» процессы для однородных сред.
развитие этого способа показало наличие
Дальнейшее
больших затруднений в его
применении к изучению более сложных моделей сред (горизонтально-слоистый
разрез, радиально-неоднородный разрез и т.д.),
поэтому потребовалось
создание новых способов оценки влияния поляризуемости на ИПП. Таким
образом, появилась методика изучения закономерностей влияния вызванной
поляризации на ИПП на основе приближенного теоретического рассмотрения
этих процессов во временной области.
30
Глава 2 МЕТОДИЧЕСКИЕ ОСНОВЫ ИНТЕРПРЕТАЦИИ
ПЕРЕХОДНЫХ ПРОЦЕССОВ В ЗМПП
Теоретической основой метода ЗМПП являются фундаментальные работы
Л.Л. Ваньяна,
Ф.М. Каменецкого, А.А. Кауфмана, Г.М. Морозовой, О.А.
Скугаревской, В.А. Сидорова, А.Н. Тихонова, П.П. Фролова, Д.Н. Четаева,
С.М. Шейнманна.
Современным
исследованиям
в
области
создания
алгоритмов
интерпретации данных ЗМПП, в том числе на основе инверсии данных
многокомпонентных нестационарных зондирований, посвящены работы М.И.
Эпова, Е.Ю. Антонова, В.П. Соколова, Л.А. Табаровского и др.
В ЗМПП одним из основных способов интерпретации кривых зондирований
является метод подбора, который заключается в последовательном изменении
геоэлектрической модели разреза и расчете прямой задачи для нее до того
момента, пока экспериментальная и теоретическая кривые не совпадут с
требуемой точностью. Для обработки данных импульсной индуктивной
электроразведки были разработаны и программно реализованы
различные
автоматизированные интерпретационные системы [34, 66].
Способ интерпретации путем пересчета кривых ЭДС от времени в кривые
зависимости кажущейся продольной проводимости от кажущейся глубины был
предложен В.А. Сидоровым совместно с В.А. Тикшаевым. Этот способ
позволяет определять для плоскослоистого разреза продольную проводимость
для каждого момента времени и определять соответствующую ей глубину.
Изменение
характера
нарастания
Sk
или
S
указывает
на
наличие
геоэлектрических границ. Перестроение кривых  (t ) или Z (t ) (значение ЭДС,
нормированное на 1 Ампер) в зависимость S k ( H k ) фактически является
приближенным
решением
обратной
задачи,
т.
е.
позволяет
по
экспериментальной кривой становления поля определить строение слоистого
31
разреза. И для этого не обязательно сравнивать экспериментальные кривые с
набором палеток (теоретических кривых) – так называемая, беспалеточная
интерпретация.
Изучая становление поля в ближней зоне можно проследить изменение
параметров разреза от отдельных горизонтов до суммарных, характеризующих
весь разрез, если при интерпретации использовать «плавающую плоскость» (Sплоскость). Для плоскослоистого разреза в поздней стадии процесс становления
поля до непроводящего основания можно заменить эквивалентной моделью
плоскости с продольной проводимостью, равной суммарной продольной
проводимости всего напластования.
Возможность
применения
в
интерпретационном
алгоритме
модели
«плавающей плоскости» основывается на работе П.П. Фролова [59], который
показал, что над горизонтально-слоистыми разрезами с непроводящим
основанием в поздней стадии становления поле асимптотически стремится к
сигналу от проводящей плоскости с продольной проводимостью всего
напластования.
математический
Модель
аппарат,
«плавающая
плоскость»
позволяющий
быстро
–
и
это
упрощенный
эффективно
вести
интерпретацию. Разрез от земной поверхности до глубины H заменяют более
простой моделью – эквивалентной плоскостью, считая остальную часть разреза
изолятором. Для такой модели рассчитывают кривую становления.
Переходный
процесс
регистрируется
перемещением
измерительного
устройства (применяются квадратные приемно-генераторные установки) по
профилю. В каждом пункте наблюдений фиксируется спад ЭДС от времени
(переходный процесс). Измерения проводятся через очень малые (микромиллисекундные) интервалы времени после выключения питающего тока в
генераторном контуре. Все измеренные значения ЭДС на различных временах
(задержках)
измерительной
аппаратуры
(в
аппаратуре
«Каскад-М»
используется 6 диапазонов с различными плотностями задержек времени) при
дальнейшей интерпретации нормируют к силе тока; получаем значение ЭДС,
приходящееся на 1 А в каждый момент времени t – Z(t).
32
Первичная обработка результатов измерений в виде Z i (t ) сводится к
получению однозначного массива Z (t j ) , где t j – массив строго возрастающих
значений времени – задержек современной компьютеризированной аппаратуры
типа «Каскад». Эта операция, называемая сглаживанием, проводится опытным
геофизиком-интерпретатором при участии оператора, проводившего полевые
измерения. На этом этапе устраняются аппаратурные погрешности, связанные с
разными режимами измерений разных значений Z i , устраняются внешние
помехи и т.д. Критерием получения реального массива Z (t j ) при сглаживании
является возможность численной трансформации этих значений в кажущиеся
интерпретационные параметры.
Очень важно точно восстановить сигнал индукционного становления, с
которым связана дальнейшая интерпретация и корректная оценка параметров,
характеризующих геоэлектрический разрез.
Огромный опыт обработки и интерпретации результатов ЗСБ показал, что
наиболее полной разрешенности можно получить по трансформациям Z (t ) в
кривые S k ( H k ) [7, 48, 50, 52, и др.].
Кажущаяся продольная проводимость S k в общем случае рассчитывается по
формуле [18]:
Sk 
где
F k 
,
Z t 
(2.1)
Z (t ) , В/A – измеренные значения ЭДС в соответствующие моменты
времени; функция F(k) рассчитывается, исходя из [75]:
F (k ) 
2
L
где k   h 
4k  1
1
1 k 2
2  k2
 


  k
k2
k 2 1
k 2 1

,


t 
 ;   4  107 , Гн/м – абсолютная магнитная проницаемость в
  S 
вакууме; h, м – глубина залегания проводящей плоскости; S, См – продольная
проводимость плоскости; t, с – время регистрации от момента выключения тока
в питающем контуре; L, м – сторона квадратной приемно-генераторной петли.
33
Формула для расчета кажущейся глубины имеет вид [18]:
H k  m  R  qгл.
или
где
Hk 
qгл
–
коэффициент
L
 k  qгл ,
2
глубинности,
(2.2)
определенный
различно
для
полупространства и над слоем ограниченной мощности (переслаиванием на
изоляторе) и равный 0,67 и 0,75 соответственно.
m
где
1
t 
,
h

R
  S 
R – радиус круглой приемно-генераторной петли эквивалентный по
площади квадратной петле со стороной L    R .
Введение функции глубинности позволило приближенно, но оперативно
решать
обратные задачи применительно к ЗСБ [48], но отмеченная
двойственность
в
определении
qгл
осталась
неисследованной.
При
практической интерпретации коэффициент qгл корректируется по результатам
каротажа удельных сопротивлений  k и сравнивается S k ( H k ) по нему. Такой
подход создает массу сложностей и неточностей (большие погрешности в
определении кажущейся глубины), не связанных с истинными возможностями
метода ЗМПП при выявлении «тонких» особенностей строения разреза.
Возможность решения прямых и обратных задач методом «кажущейся
продольной
проводимости»
и
его
универсальность
при
определении
геоэлектрики в различных сложных средах позволяет выявлять самые
малоконтрастные ее особенности по результатам ЗМПП.
Введение и использование динамических параметров сигнала, уточняющих
и характеризующих кривую спада переходного процесса, является актуальным
и перспективным для определения реальной разрешающей способности метода,
поскольку существует «тесная взаимосвязь между формой переходного
электромагнитного сигнала и геоэлектрической структурой среды» [13] .
34
2.1. Динамические параметры сигнала в ЗМПП
Качественный и количественный анализ сигналов
становления поля в
ближней зоне значительно облегчается при использовании их динамических
параметров в виде [87]:
n(t ) 
t  Z (t )
,
Z (t )
(2.3)
n( k ) 
k  F (k )
,
F (k )
(2.4)
n(u ) 
u  f (u )
,
f (u )
(2.5)
где k, u – вспомогательные интерпретационные безразмерные параметры; f(u)
– функция, характеризующая спад переходного процесса над ОПП,
u
2   R

,   2 2    t [18].

Вводимые динамические параметры являются ограниченными по величине и
определяют наклон кривых Z (t ) , F (k ) и f (u )
являются показателями степени
от времени t , от k и u
и
при характере спада соответствующих
функций (например, n(t ) является показателем степени t сигнала становления
Z (t ) ~ A  t n ).
Между параметрами по (2.3) и (2.4) имеется аналитическая связь [87]:
n(t )  n(k ) 
k s (t )
.
k (t )
(2.6)
Уменьшение угла наклона (более крутой спад), т.е. повышение их
абсолютного значения
при
качественной
интерпретации
указывает на
увеличение сопротивления текущего слоя в разрезе. Но такой анализ неточен в
силу того, что эти параметры многозначные (например, для ОПП на ранней
стадии n(t ) имеет близкое значение и стремится при t  0 к -1, на поздней
стадии
при
t   стремится к -2,5), т. е. даже для ОПП имеются два
преимущественных наклона. Такая же неоднозначность наблюдается и для
проводящей плоскости, т. е.:
35
для однородного полупространства:
n(k) при t  0 на ранних стадиях стремится к «-1», на поздних стадиях при
t   к «-4»;
n(t) при t  0 на ранних стадиях стремится к «-1», на поздних стадиях при
t   к «-2,5»;
n(u) при t  0
на ранних стадиях стремится к «5», на поздних стадиях при
t   к «2»;
для проводящей плоскости:
n(t) при t  0 на ранних стадиях стремится к «-1», на поздних стадиях при
t   к «-4».
Однако динамические параметры имеют единственные особые точки,
характерные для всей кривой – точки перегиба (рис. 2.1, 2.2). Прямые расчеты
показали, что точка перегиба для ОПП в линейном масштабе практически
незаметна, однако явно выражена при логарифмическом масштабе по оси
времени (рис. 2.3). В точках касания n(t ) и nh (t ) равны. Расчётами установлено,
что это равенство сохраняется при любом времени, в том числе и в точке
перегиба.
0.5
0
-0.5
-1
-1.5
-2
-2.5
-3
-3.5
-4
-4.5
0
0.5
1
1.5
k
n`(k)
n(k)
Рис. 2.1. Динамический параметр n(k) и его производная n'(k)
2
36
0
-0.5
-1
ns(t)
-1.5
-2
-2.5
-3
-3.5
-4
-4.5
1
10
100
1000
10000
100000
t, мкс
Рис. 2.2. Динамический параметр n(t) от поверхностной проводящей плоскости
0
-0.5
t=t0=340 мкс
-1
n(t)
-1.5
-2
-2.5
-3
-3.5
-4
1
10
100
1000
10000
100000
t, мкс
1
2
3
Рис. 2.3. Динамические параметры (Яхина И.А., 2011): 1 – динамический параметр от поверхностной
проводящей плоскости ns(t); 2 – динамический параметр от проводящей плоскости, расположенной на глубине
nh(t); 3 – динамический параметр n(t) от однородного проводящего полупространства
37
Дифференцируя выражение (2.4) находим производную параметра n(k ) в виде
n(k ) 
n( k ) 
F (k ) 
.
1  n( k )  k

k 
F (k ) 
График этой зависимости приводится на рис. 2.4 и имеет единственный явно
выраженный
минимум. По расчетам значение экстремальной точки
определяется однозначно при k  0.5 .
На практике Z (t ) наиболее полно представляется в билогарифмическом
масштабе, поэтому
исследование
экстремумов
логарифмических производных динамических параметров,
т. е. функций
вида
особенно
n(t )  t n(k )  k n(u )  u
,
,
n(t )
n( k )
n(u )
актуально
для определения точек перегиба и внутренних
связей между ними, в первую очередь для базовых моделей (однородное
полупространство, проводящая
Функция
плоскость, проводящий
шар, диск и др.).
n(k )  k
 0 (рис. 2.4) на ранних и на поздних стадиях, следовательно,
n( k )
эта функция имеет единственный положительный максимум, при k  0.5
значение
n(k )  k
 0.5 .
n( k )
Логaрифмическая производная
определяется
по
формуле
 (t )
 (t ) 
над ОПП в самом общем виде
n (u )  u
n(t )  t
 u
n(t )
2  n(u )
[87]. По
рис. 2.5
определяем, что функция  (t ) имеет единственный максимум величиной 0,27
при u э  0,87 , который и определяет положение точки перегиба ( t п ) n(t ) для
ОПП и его единственное значение n(t п )   n(u э )  1,65  f (u э )  0,591622 .
1
2
Графики n(k ) , n(k ) пересекаются в одной точке и имеют близкие значения в
точке экстремума k э  0.5 .На ранних стадиях n(k э )  n(k э )  1.74n(k ) , n(k )  1 ,
на поздних стадиях функция n(k )  4 , n(k )  0 . В соответствии с (2.6) должна
быть связь n(t ) и n(k ) .
38
1
0.5
0
n(k), n`(k)
-0.5
-1
-1.5
-2
-2.5
-3
-3.5
0
0.5
1
1.5
2
k
1
2
3
4
5
Рис. 2.4. Динамический параметр n(k) и его производная n'(k) (ρ=40 Ом·м, L=1000 м): 1 – логарифмическая
производная n(k) от ОПП; 2 – n'(k) от ОПП, 3 – n'(k) от поверхностной проводящей плоскости, 4 –
логарифмическая производная n(k) от поверхностной проводящей плоскости; 5 – n(k) от ОПП
0.3
n'(u)·u/n(u)
0.25
0.2
0.15
0.1
0.05
0
0
1
2
3
4
5
6
u
Рис. 2.5. Обобщенная логарифмическая производная δ(t) от однородного проводящего
полупространства (Яхина И.А., 2011)
39
2.2. Методика интерпретации данных ЗМПП способом «плавающей»
плоскости с помощью динамических параметров сигнала
Метод КПП базируется на расчёте сигнала становления поля Z (t ) через
сигнал «плавающей» проводящей плоскости (ПП) [49]. При этом
Z (t ) 
F (k )
,
S k (t )
(2.7)
где
F (k ) 
4k  1
1
1 k2
2  k2
 


  k
k2
k 2 1
k 2 1

.


(2.8)
Параметр k в формуле (2.8) рассчитывается как
k (t ) 
где

2
t
,
h(t ) 
L
  S k (t ) 
(2.9)
м – сторона петли квадратной установки; t , с – время регистрации в
L,
секундах от момента выключения тока в квадратном контуре со стороной
h(t ) ,
L;
м – расстояние от поверхности до проводящей плоскости; S k (t ) , См –
продольная проводимость плоскости в момент времени t , расположенная на
глубине h(t ) ;   4    10 -7 , Гн/м – магнитная проницаемость вакуума.
Безразмерный параметр
k
можно разложить на составляющие и
представить в виде суммы параметров k
для проводящих плоскостей,
расположенных на глубине и на поверхности:
k (t )  k h (t )  k s (t ) ,
где
k h (t ) 
2 t
2
.
h(t ) , а k s (t ) 
L S k (t )
L
(2.10)
Расчетный параметр k определяется также из решения уравнения:  (t )   (k ) ,
 (t )   L
где
Z (t )
Z (t )
F (k )
[13],  (k )  2 2
2
Z (t )
F (k )
– производная по времени переходного импеданса
(2.11)
Z (t ) ;
F (k ) –
производная функции F (k ) (2.8) по параметру k .
Подставив (2.3) и (2.4) в (2.11), получим уравнение  (t )   (k ) для расчёта
параметра k для каждого времени t в виде:
40
 L
n(t )
n( k )
.
 2
t  Z (t )
k  F (k )
(2.12)
С учетом (2.6) следует, что
k s (t )  k (t )
n(t )
,
n( k )
(2.13)
и можно составить универсальную формулу для решения прямой задачи, т. е.
для расчёта сигнала Z (t ) от проводящей плоскости на глубине через сигнал
поверхностной плоскости Z s (t ) (рис. 2.6):
Z (t )  Z h (t ) 
F (k )
Z s (t ) .
F (k s )
(2.14)
Здесь Z s (t ) – сигнал от проводящей плоскости с продольной проводимостью
S k (t ) на нулевой глубине в соответствии с (2.7) и (2.12):
S k (t ) 
2 n( k ) t
.
L n(t ) k
(2.15)
Подставляя в (2.15) зависимость (2.13), получаем в самом общем случае
S k (t ) 
2 t
.
L k s (t )
(2.16)
Следовательно, если задана S k (t ) , то соответствующий сигнал Z (t ) может быть
легко рассчитан, по (2.14) с учетом (2.8) и (2.10).
Таким образом, если расположить проводящую плоскость на определенной
глубине h , то в момент времени t  t0 сигнал от однородного проводящего
полупространства
Z (t )
точно совпадает по значению с сигналом от
проводящей плоскости Z h (t ) , которую можно назвать «предельной» (рис. 2.6). t0
– характерное время регистрации для предельной плоскости, определяемое
характерными точками
предельной
проводящей
динамических параметров (рис. 2.7). Сигнал от
плоскости
Z h (t ) ,
h(t 0 )  H (t 0 ) , в точности совпадая с таковым
определить его удельное сопротивление  .
расположенной
на
глубине
от ОПП Z (t 0 ) , позволяет
41
1000
100
10
Z(t), В/А
1
0.1
0.01
0.001
0.0001
1E-05
1
10
100
1000
10000
100000
t, мкс
1
2
3
Рис. 2.6. Сигналы от однородного проводящего полупространства и проводящей плоскости (ρ=50 Ом·м,
L=1000 м, t0=340 мкс) (Яхина И.А., 2011): 1– сигнал Z(t) от однородного проводящего полупространства; 2 –
сигнал от проводящей плоскости, расположенной на поверхности, Zs(t);3 – сигнал от предельной проводящей
плоскости, расположенной на глубине h(t0), Zh(t).
42
Рис. 2.7. Определение
t 0 при L=100 м, ρ=40 Ом·м над ОПП (Яхин А.М., Яхина И.А., 2010): 1 –сигнал над
ОПП; 2 – динамический параметр n(t) от ОПП; 3 – сигнал от проводящей плоскости; 4 – динамический
параметр n(t) от проводящей плоскости .
Зная t 0 , можно вычислить параметры проводящей плоскости – значение
продольной проводимости проводящей плоскости на предельной глубине [70],
которые в точности определяют верхний слой ОПП. Положение проводящей
плоскости определяется через динамические параметры по формуле:
h(t ) 

L
n(t ) 
 .
 k (t )  1 
2
 n( k ) 
(2.17)
Функция (2.17) определяет точное положение ПП в полупространстве при
любом t , продольная проводимость которой постоянна и рассчитывается также
по формуле (2.16). В самом общем случае положение ПП для любого разреза
устанавливается преимущественно динамическим параметром n(t ) , но h(t )
отличается от H k (t ) .
2.3. Понятие глубинности. Расчет кажущейся глубины с помощью
динамических параметров сигнала.
Практическая ценность решения обратных задач с использованием
продольной проводимости
существенно повышается, если ввести понятие
43
глубинности
[18, 49] или кажущейся глубины, т. е. возможности такого
определения положения плоскости H (t ) с продольной проводимостью S (t ) ,
определяющей тот же сигнал, как по (2.7). В случае плоскослоистого разреза
H (t ) эквивалентна соответствующей мощности слоя непосредственно под
установкой в момент времени t с продольной проводимостью S k (t ) . Очевидно,
что решение может быть найдено в общем случае, если величина k (t ) не
меняется и можно единственным образом представить его в виде
k (t ) 
2
tН 
,
H (t ) 

L
  S (t ) 
(2.18)
где t H - характерное время для этой плоскости на глубине H (t ) . Определим ее
исходя из равенства
H (t ) 
tH
  S (t )
т. е.
t H  H (t )    S (t ) .
(2.19)
Кажущееся удельное сопротивление текущего напластования определяется
в общем случае по формуле [18]:
 k (t ) 
dH k (t )
.
dS k (t )
(2.20)
Для введения функции глубинности в ОПП нет необходимости в
специальных физических допущениях, достаточно удовлетворить условиям
точного решения обратной задачи в виде
d (t )
H (t ) S (t )
H (t )
 0 ,  (t ) 
.
 const ,

dt
S (t )
H (t ) S (t )
Интегрируя эти уравнения с учетом начальных условий, имеем функцию для
расчета глубинности ОПП:
H (t )   (0)  S (t ) ,
 (0)  lim
t 0
(2.21)
H (t0 ) h(0)
2 n(k (t0 ))

, h(0)  h(t0 )  0,5  L  kh (t0 ) , S (0)  S k (t0 ) 
. В [70, 87]
S (t0 ) S (0)
L n(t0 )
показано, что при k (t0 )  0,5 глубина до ПП h(t 0 )  H (t 0 ) . Таким образом, k (t0 )
можно представить в виде:
44
k (t0 )  kh (t0 )  k s (t0 ) ,
где kh (t0 )  0,15872
k s (t0 ) 
2 t0
и остается постоянной независимо от S 0 и
L   S0
равной k s (t 0 )  0,34128 . Следовательно, характерное время для этой плоскости
при любом S 0 определяется формулой
t0 
L
k s (t0 )    S 0 ,
2
(2.22)
т. е. иначе чем по (2.19). Поэтому, для получения более универсального
алгоритма
определения
глубинности
H (t ) ,
необходимо
удовлетворить
возможность автоматического ее соответствия для ОПП. Уравнивая (2.18) и
(2.9) получим связь между h(t ) и H (t ) , что позволяет решать прямые задачи по
известным S (H ) , например, по каротажу, или синтезируя разные модели
плоскослоистых сред. Для проводящей плоскости S (t0 )  S0 , следовательно,
только на глубине H (t0 )  H 0 получим удельное сопротивление, эквивалентное
ОПП, в остальных случаях  ПП будет зависеть от времени, это зависит от
способа введения глубинности ориентированного на точное решение задач для
ОПП. Поскольку 0  H (t ) , то и 0   ПП (t ) .
Удельное сопротивление  ОПП определяется по формуле:

H k (t 0 )
,
S k (t 0 )
где t 0 – время регистрации данной плоскости, названной – «предельной», с
S k (t 0 ) , равной продольной проводимости слоя под установкой мощностью
H k (t 0 ) .
Результат интерпретации (рис. 2.8, 2.9, 2.10) с помощью «плавающей
плоскости»
позволяет получить однозначное решение обратной задачи
(определить удельное сопротивление
глубинности по (2.21).
)
для ОПП с учетом введения
45
60
50
ρ,Ом·м
40
30
20
10
0
0
100
200
300
400
500
600
Hk, м
1
2
Рис. 2.8 Результаты интерпретации способом «плавающей плоскости» над ОПП (ρ=40 Ом·м, L=1000 м)
(Яхина И.А., 2011): 1 - точное решение; 2- приближенное решение по В.А. Сидорову.
70
60
50
ρ, Ом·м
Hk=H(t), м ρ, Омм
4,2408842
50
5,9863113
50
8,4232443
50
11,807921
50
16,479966
50
22,881856
50
31,577204
50
43,261797
50
58,754736
50
78,943729
50
104,63091
50
136,18691
50
173,11094
50
214,47408
50
260,67442
50
314,3573
50
379,37627
50
459,31631
50
557,27523
50
676,25159
50
819,60636
50
991,4025
50
1196,5893
50
40
30
20
10
0
0
200
400
600
800
1000
1200
1400
Hk, м
1
2
Рис. 2.9. Результаты интерпретации способом «плавающей плоскости» над ОПП (ρ=50 Ом·м, L=1000 м):
1- точное решение; 2- приближенное решение по В.А. Сидорову.
46
Hk=H(t),м ρ=Ом·м
500
4,2408842
500
5,9863113
500
8,4232443
500
11,807921
500
16,479966
500
22,881856
500
31,577204
500
43,261797
500
58,754736
500
78,943729
500
104,63091
500
136,18691
500
173,11094
500
214,47408
500
260,67442
500
314,3573
500
379,37627
500
459,31631
500
557,27523
500
676,25159
500
819,60636
500
700
600
500
ρ, Ом·м
2,8854668
400
300
200
100
0
0
200
400
600
800
1000
Hk, м
1
2
Рис. 2.10. Результаты интерпретации способом «плавающей плоскости» над ОПП(ρ=500 Ом·м, L=500 м):
1- точное решение; 2- приближенное решение по В.А. Сидорову.
На рис. 2.8-2.10 показана возможность точного решения (однозначного
определения удельного сопротивления) методом «кажущейся продольной
проводимости» над однородным проводящим полупространством. При этом
анализируются сигналы становления в ближней зоне по формуле [20] в виде:
Z (t ) 
где
8 
f (u ) ,
L
(2.23)


 2u

2  u
1
1

3

f (u )  
 e 4 u 
 (2u )  (2u )  2u 2  
 2 2 u  ,
2
4

  2

 

u
L 
2  2
u
, (u )   e d .
4 t
 0
(2.24)
Решение обратной задачи – решение уравнения:   (t )   (u)
 (u ) 
  f (u )  u 3
f 2 (u )
,   (t ) 
  n(u )  u 2
f (u )
,
(2.25)
которое является кубическим уравнением, не зависящим от  и S и решается
однозначно относительно u .
Выражаем H (t ) по (2.19) и с учетом (2.23) получаем:
47
H (t ) 
С учетом (2.7) H (t ) 
Z (t ) L
S (t )
f (u ) 8
F (k )   L
S (t )
f (u )  S (t ) 8
и получаем функцию для расчета
глубинности для ОПП:
F (k )   L
f (u ) 8
H 0 (t ) 
(2.26)
Похожим образом можно выразить универсальную функцию глубинности
через динамические параметры
для любого момента времени. Значение
функции глубинности позволяет оценить мощность пород под установкой
затронутой вихревыми токами, сигналы от которых вносят основной вклад в
регистрируемый сигнал, а также выявить в нем долю от текущего слоя.
Естественно, вводимая функция автоматически должна включать глубинность
для ОПП [70, 87].
Производная (2.23) по t имеет вид:
Z (t ) 
где
L
u(t ) 
4
8 
f (u )  u(t ) ,
 L
  1   32
u
.
 t  
 t  2 
2t
Z (t )  
В соответствии с (2.19) Z (t )  
8 
u
f (u ) 
 L
2t
8
u H
f (u )  
и выражаем
 L
2t S
H (t )  
Z (t )t L
 S (t ) .
f (u )u 4
С учетом (2.3), (2.5), (2.7) получаем:
H (t )  
n(t )  Z (t ) L
 S (t ) ,
n(u )  f (u ) 4
где S (t ) выражаем по (2.7) и получаем функцию для расчета глубинности для
произвольных горизонтально-слоистых сред
H (t )  
n(t ) L F (k )
n(u ) 4 f (u )
Умножим и разделим уравнение на 2, получим
(2.27)
48
H (t )  
n(t )  2 L F (k )
n(t )
 2
H 0 (t ) .
n(u ) 4  2 f (u )
n(u )
В соответствии с (2.2)
qгл 
qгл 
2  H (t )
,
Lk
(2.28)
L F (k ) 2
8 f (u ) Lk
qгл 
 F (k )
4 k  f (u )
,
(2.29)
где u определяется как решение уравнений (2.25).
В этих формулах к и u являются решениями уравнений  (t )   (k ) и  (k )   (u) .
Для поздних стадий становления (t→∞, к→∞, u→0) зависимость между ними
имеет вид
к
1 3 75  
.
u
256
При k   (u  0, t  ) для ОПП
qгл  0.619 .
(2.30)
При k   (u  0, t  ) для проводящей плоскости
qгл  0.715 .
Коэффициент qгл не зависит от  и при (t  0, u  ) q гл  0 . Точный вид
этого сложно зависящего от k коэффициента для ОПП приведен на рис. 2.11.
0.7
0.6
qгл
0.5
0.4
0.3
0.2
0.1
0
0
1
2
k
3
4
5
Рис. 2.11. Функция коэффициента глубинности для ОПП (Яхина И.А., 2011)
49
Для
оценки разрешающей способности метода
квадратными
установками
при
решении
задач
ЗМПП совмещенными
геоэлектрики
была
проанализирована трехслойная модель среды [83]. Для расчета кажущегося
сопротивления способом «плавающей плоскости» использовался альбом
палеток кривых зондирований становлением поля в ближней зоне [2]. За основу
принята трехслойная модель разреза типа «H» с проводящим основанием,
подобранная в соответствие с реальным геофизическим разрезом Ерыклинского
участка Ромашкинского нефтяного месторождения (глава 3). Указанный
интервал включает нефтеперспективные отложения, начиная с середины
среднего карбона (верейский ярус) и полностью нижний карбон (турнейский
ярус).
Параметры подобранной палеточной модели (исходя из условий h2/h1=0,5;
ρ2/ρ1=0,5; ρ3/ρ1=1):
- сопротивления первого слоя ρ1 =70 Ом·м, мощность первого слоя h1= 800 м;
- сопротивление второго слоя ρ2=35 Ом·м, мощность второго слоя h2=400 м;
- сопротивление третьего слоя ρ3=70 Ом·м, мощность третьего слоя h3=∞.
По известным формулам рассчитаны и приведены Z(t) для ρ=35 Ом·м и
ρ=70 Ом·м для установок L=1000, 500 м, а также соответствующие сигналы для
эквивалентных по моменту дипольных установок, определенных по палеткам
 1  / h1  из альбома А.А. Кауфмана и др. [2] (рис. 2.12). В верхней части
рисунка приведены относительные превышения переходных процессов θ(t) для
обеих эквивалентных дипольных установок над сигналами от однородного
полупространства с ρ=70 Ом·м в совмещенных петлях.
По кривым θ(t) оцениваются начальные времена tн, начиная с которых
сигналы для совмещенных и дипольных установок совпадают в пределах
принятой на практике 5% погрешности измерений. Соответственно для
установки с L=500 м, tн =1,07 мс и для L=1000 м, tн =5,25 мс.
50
1
0.1
0.01
Z(t), В/А
0.001
0.0001
1E-05
1E-06
1E-07
1E-08
100
1000
10000
100000
1000000
t, мкс
1
2
3
4
5
6
7
8
9
10
11
Рис. 2.12 Переходные процессы ЗМПП в совмещенных и дипольных установках (И.А. Яхина, А.М. Яхин,
2014): 1 – сигнал для дипольной установки, эквивалентной по моменту совмещенной установке L=500м,
Zdip(t); 2 – сигнал для L=500м от однородного полупространства для ρ=70 Омм, Zpp(t); 3 - сигнал для L=500м
от однородного полупространства для ρ=35 Омм, Zpp(t); 4 – относительное превышение модельного сигнала
над сигналом от ОПП для ρ=70 Омм, θ(t); 5 – сигнал для дипольной установки, эквивалентной по моменту
совмещенной установке L=1000м, Zdip(t); 6 – сигнал для L=1000м от однородного полупространства для ρ=70
Омм, Zpp(t); 7 – сигнал для L=1000м от однородного полупространства для ρ=35 Омм, Zpp(t); 8 –
относительное превышение модельного сигнала над сигналом от ОПП для ρ=70 Омм, θ(t); 9 – полевой сигнал
для совмещенной установки L=1000м; 10 – теоретический сигнал для L=1000м; 11 – теоретический сигнал для
L=500м.
51
При практической обработке и интерпретации влияние слоя с большей
проводимостью (ρ=35 Ом·м) в геоэлектрическом разрезе типа «H» можно
характеризовать как небольшое и отнести разрез к слабоконтрастному. В
трансформации  (t ) (рис. 2.13) кажущиеся значения весьма неопределенно
оценивают истинные параметры геоэлектрического разреза и совершенно
отсутствует непосредственная оценка геометрии плоскослоистого разреза. Для
интерпретации  (t ) требуются дополнительные палетки и использование
методик разработанных А.А. Кауфманом, Б.И. Рабиновичем, Г.А. Исаевым и
др.
Использование
трансформации
и
Sk (H k )
k (H k )
[49]
позволяет
непосредственно оценить геоэлектрику плоскослоистой среды. В работах [75,
78, 87] показано, что при определении H k по (2.2) qгл (рекомендованный В.А.
Сидоровым
для
практики
равный
0,67
и
0,75)
является
монотонно
возрастающей функцией от k и при больших значениях асимптотически
стремится к этим значениям в зависимости от геоэлектрики разреза. Для
сравнения с  (t ) диаграмма  k (t ) приведена на рис. 2.13.
Использование динамических параметров n(t ) ,  (t ) позволяет повысить
точность расчета S k (t ) , H k (t ) и, тем самым,  k (t ) , приблизив их к истинным
значениям. На рис. 2.14 приведены результаты такой трансформации для
L=1000 м в виде кривых S k ( H k ) и  k ( H k ) . Удельное сопротивление первого слоя
определяется непосредственно на диаграмме  k (t ) величиной в 70 Ом·м.
Контрастный локальный минимум приводит ρk(tmin)=35 Ом·м при tmin=16,08 мс.
Однако на более позднем времени tmax=42,80 мс возникает локальный максимум
ρk=74 Ом·м, далее  k (t ) плавно уменьшается до значений 66 Ом·м и затем
возрастает, асимптотически приближаясь к ρ1=70 Ом·м. Границы раздела среды
определяются по точкам перегиба кривой S k ( H k ) , h1 – мощность первого слоя
определяется величиной 796 м (t1=11,79 мс) (рис. 2.14), граница перехода к
третьему бесконечному слою определилась на глубине 1200 м (t2=32,74 мс).
52
100
90
80
70
ρk, Ом·м
ρτ, Ом·м
60
50
40
30
20
10
0
100
1000
10000
100000
1000000
10000000
t, мкс
1
2
3
4
5
Рис. 2.13 Палеточные кривые ρτ для трехслойной модели среды в сопоставлении с рассчитанными
значениями ρк по трансформациям Sk(Hk) (И.А. Яхина, А.М. Яхин, 2014): 1 – палеточные значения ρτ для
дипольной установки, эквивалентной по моменту совмещенной установке L=500м; 2 – рассчитанные значения
ρk; 3 – величина ρ=70 Ом·м; 4 – величина ρ=35 Ом·м; 5 - палеточные значения ρτ для дипольной установки,
эквивалентной по моменту совмещенной установке L=1000м.
53
80
70
60
ρ, Ом·м
Sk, См
50
40
30
20
10
0
0
500
1000
1500
2000
2500
3000
H, м
1
2
Рис. 2.14
Определение геоэлектрических границ по трансформациям кажущейся продольной
проводимости от кажущейся глубины (L=1000 м): 1 – Sk(Hk), 2 – ρk(Hk)
Таким образом, правильное определение H k
определяет точность оценки
 k . На рис. 2.15 приведены примеры определения геоэлектрических границ,
наблюдаемых по
Sk ( H k ) (см. рис. 2.14), по
точкам перегиба n(t)
и
отрицательным минимумом и положительным максимумом (экстремумам)  (t ) .
Определяя  (t ) , можно оценить число напластований, даже не рассчитывая
Sk ( H k )
[78].
Монотонно
спадающая
экспериментальная
кривая
имеет
характерные точки перегиба по n(t ) и множество экстремумов по  (t ) , которые
определяют границы разделов сред с различными значениями
k
и
поляризуемости  k . При использовании описанных динамических параметров
существенно повышается информативность регистрируемого сигнала при
оценки истинной разрешающей способности метода.
54
3
2
ρk(t)
δ(t), n(t)
1
0
-1
-2
-3
1000
10000
100000
1000000
t, мкс
1
2
3
4
Рис. 2.15 Определение геоэлектрических параметров трехслойной модели среды по динамическим
параметрам сигнала: 1 - логарифмическая производная δ(t); 2 – ρk(t); 3 – динамический параметр n(t) от
трехслойной модели; 4 – динамический параметр n(t) от ОПП.
55
Таким образом, методику обработки
данных ЗМПП путем пересчета
кривых зависимости ЭДС от времени Z (t ) в кривые зависимости кажущейся
продольной проводимости от кажущейся глубины
Sk (H k ) ,
с учетом
приведенного выше алгоритма интерпретации, можно представить обобщением
в виде блоков расчетов вспомогательных
и основных параметров,
определяющих геоэлектрическое строение разреза (рис. 2.16). Рассчитываются
вспомогательные параметры и функции, которые в дальнейшем используются
при расчете геоэлектрических параметров разреза – продольной проводимости,
мощности, а также соответствующей электропроводности или сопротивления
слоя. За счет использования в алгоритме интерпретации новых дополнительных
параметров (динамических) значительно повышается точность интерпретации и
разрешающая способность ЗМПП.
Методический
комплекс
усовершенствованный
обработки
алгоритм
динамических параметров
данных
ЗМПП,
интерпретации
с
включающий
использованием
реализован в программном комплексе Microsoft
Excel с помощью подпрограммных функций, созданных в Visual Basic. Данный
алгоритм также реализован в виде программы экспрессной обработки данных
ЗМПП при использовании языка программирования C++.
С помощью предложенного методического комплекса была решена
следующая геологическая задача: проведение межскважинной корреляции с
целью оконтуривания залежей углеводородов и оценка перспективности
нефтенасыщения на Ерыклинском участке Ромашкинского месторождения. По
результатам численных трансформаций построены графики
кажущейся продольной проводимости
зависимостей
от абсолютных отметок глубин и
комплексные геоэлектрические разрезы (графический материал приведен в
главе 3).
56
Ввод исходных данных
Ввод массива (Z, В/A; t, мкc)
измеренных данных
Ввод постоянных значений
и констант (L, м)
Процедура сглаживания
Сглаживание: устранение аппаратурных помех,
устранение ошибок регистрации сигнала
Получение однозначного массива (Z, t)
строго возрастающих времен задержек
аппаратуры
Расчет вспомогательных параметров
Расчет n(t)
Расчет k
Расчет F(k)
Расчет u
Расчет f(u)
Расчет n(u)
Расчет основных параметров
Расчет Sk
Расчет Hk
ρk(Hk), σk(Hk)
Рис. 2.16 Блок-схема методического комплекса обработки данных ЗМПП.
57
Выводы
Реализованный в данной работе методический комплекс обработки и
интерпретации данных ЗМПП основывается на расчете сигнала становления
поля от однородного проводящего полупространства через сигнал проводящей
плоскости,
расположенной
на
предельной
глубине
с
использованием
динамических параметров сигнала. Эти параметры, а также их производные и
логарифмические производные имеют характерные точки – экстремумы и
точки
перегиба,
определяющие
взаимосвязь
между
базовыми
геоэлектрическими моделями сред. Таким образом, получено точное решение
обратной задачи для модели ОПП через модель «плавающей плоскости». Для
произвольного плоскослоистого разреза такое решение сводится к тому, чтобы
задавать глубину проводящей плоскости, которая бы точно определяла
удельное сопротивление и мощность верхнего слоя полупространства в
заданный момент времени. То есть любой плоскослоистый разрез можно
заменить набором полупространств, ввести функцию для расчета глубины до
проводящей плоскости таким образом, чтобы она в автоматическом режиме
определяла глубинность для полупространства.
По результатам второй главы можно сформулировать первое и второе
защищаемые положения:
– с использованием модели «плавающей» проводящей плоскости и
динамических параметров сигнала определена функция глубинности, которая
позволяет оценить глубину проникновения вихревых токов в проводящей
геологической среде под совмещенной установкой в каждый момент времени,
– предложенная методика обработки и интерпретации данных ЗМПП с
совмещенными установками с использованием динамических параметров
сигнала позволяет
повысить точность оценки геоэлектрических параметров
разреза за счет выделения слабоконтрастных по удельной электропроводности
слоев.
58
Глава 3 ЧИСЛЕННЫЕ СПОСОБЫ ВЫДЕЛЕНИЯ
НЕМОНОТОННЫХ ПРОЦЕССОВ И МЕТОДЫ ПОВЫШЕНИЯ ИХ
КОНТРАСТНОСТИ ОСОБЫМИ ПРИЕМАМИ ИХ РЕГИСТРАЦИИ
В
электромагнитной
разведке
преимуществами
обладают
методы,
основанные на возбуждении вихревых токов импульсными полями и на
измерении вторичных полей в паузах между возбуждающими импульсами.
Главное достоинство этих методов – измерение вторичных полей от вихревых
токов в изучаемых объектах
на нулевом фоне, в отсутствии сильного
возбуждающего (первичного) поля. Это облегчает изучение вторичного поля
(отклика), а значит, и параметров интересующих нас объектов.
В зондированиях становлением поля в дальней зоне (ЗСД) не всегда
удовлетворяла
детальность
и
разрешающая
способность,
а
в
МПП
существовала необходимость изучения не только локальных проводников, но и
сложной вмещающей среды. Также появилась необходимость расширить
возможности методов становления поля применительно к нефтегазопоисковым
работам. Поэтому были теоретически обоснованы, разработаны и внедрены
зондирования становлением поля в ближней зоне (ЗСБ) В.А. Сидоровым,
В.В. Тикшаевым, Г.Г. Обуховым, А.А. Кауфманом, Г.М. Морозовой и др.
Меньшая
зависимость
зондирований
в
ближней
зоне
совмещенными
установками от искажающего влияния наклона пластов и повышенная
чувствительность
отклика
по
разрезу
позволила
интерпретировать
и
относительно слабые особенности регистрируемых кривых и выделять большее
число слоев, чем в дальней зоне.
Нефтепоисковая электромагнитная разведка при индуктивном возбуждении
поля развивается за счет увеличения разрешающей способности зондирований
[8, 12, 33, 42, 52, 64 и др.]. Широкое применение наземной электроразведки с
учетом ВП в Западной Сибири (группа Моисеева В.С., СНИИГГиМС г.
59
Новосибирск)
позволило
существенно
повысить
эффективность
запасов
углеводородов,
нефтепоисковых работ.
Существенное
сокращение
прироста
обеспечиваемых поисковым и разведочным бурением на основе результатов
опережающих их сейсморазведочных работ, заставляет искать дополнительные
эффективные геофизические методы и развивать альтернативную теорию
фактов дополнительной подпитки нефтяных месторождений, согласно которой
приток углеводородов осуществляется по нефтеподводящим каналам из
резервуара, находящегося в глубинах мантии. Это подтверждается данными
сверхглубокого бурения, глубинными сейсмическими исследованиями
и
достаточно длительными (свыше 50 лет) сроками эксплуатации некоторых
месторождений, например Ромашкинского месторождения нефти в Республике
Татарстан. В процессе разработки данного месторождения было извлечено
значительно больше нефти, чем находилось на балансе и добыча продолжается.
Прогноз обнаружения и локализации каналов аккумуляции и миграции
углеводородов необходимо строить на показаниях геофизических методов
наиболее чувствительных к динамическим процессам. Определять положение
таких
каналов
следует
на
участках
кристаллического
фундамента,
приуроченных к зонам суперпозиции (наложения) максимальной динамики
аномальных параметров электрических и сейсмических геофизических полей.
Такие оценки возможны по результатам измерений в режиме непрерывного
мониторинга участков кристаллического фундамента наиболее подверженных
тектоническим процессам (рифтовых зонах и их пересечениях).
Электроразведочные работы выделяются в ряду множества технологий и
имеют реальные физические основы для обнаружения залежей углеводородов.
По
результатам глубинных
(до
2,2
км.)
зондирований
импульсными
электроразведочными методами (в варианте ЗМПП) [64, 84] и бурения новых
разведочных скважин, вскрывших кровлю фундамента в юго-восточной
Татарии и западной Башкирии удалось по новому представить рельеф дна
широтной «каньонообразной» долины К.Р. Тимергазина [41, 77] и оценить роль
60
и возможность участия долины в транзите углеводородов из глубинных толщ.
В работе [41] делается вывод о том, что рельеф (рис. 3.1) представляется как
поверхность горста пород фундамента, полость дна долины рассматривается
как континентальная рифтовая структура в плане с разломно-блоковой
тектоникой. На рис. 3.1 картирован Субхангуловский выступ, в осадочной
части разреза над ним не обнаружено заметных скоплений нефти, но близость к
долине-рифту, наличие его выделяющей разломной тектоники, указывает на
необходимость как исследования структурно-вещественного комплекса этого
блока, так и проведения опережающей электроразведки для прямого поиска
нефти в кристаллическом фундаменте.
3.1. Многоразмерные зондирования ЗМПП как основа повышения
разрешающей способности электромагнитной разведки
Совмещенные и осесовмещенные установки ЗСБ и МПП широко
применяются при поисках рудных объектов, поисках нефти и газа,
картировании водоносных горизонтов и т. д., т. е. охватывают весь спектр
поисковых и структурных задач. Поэтому решение задачи учета влияния ВПИ в
осесимметричных установках имеет очень большое значение, так как от
правильной регистрации индукционного становления зависит правильность
интерпретации и, в конечном счете, решение поисковых и структурных задач.
Кроме того при заметном влиянии ВПИ, появляется возможность в переходном
режиме изучать поляризационные свойства среды, что является весьма
эффективным способом, так как позволяет исследовать изучаемый параметр в
любое время года.
Численные способы исправления искажений полевых кривых, впервые
примененные В.А. Сидоровым для интерпретации замеченных им совместно с
А.Д. Скурихиным отрицательных процессов, способствовали пониманию
отрицательного спада как аддитивного слагаемого к индукционному процессу.
Регистрируемые сигналы при
ЗСБ совмещенными установками
рассматриваться в виде суммы сигналов [50]:
Z t   Z спи t   Z впи (t ) ,
могут
61
Рис. 3.1. Сопоставление структурных карт кровли кристаллического фундамента
с кровлей репера К4 по Серафимовско-Балтаевскому валу
(Муслимов Р.Х. и др., 2009;Яхин А.М., Шакиров М.З., Яхина И.А., 2009).
62
Z спи t  – сигналы индукционного становления в среде с эффективной
где
электропроводностью
0   ;
Z впи t 
–
поляризационная
составляющая
сигналов.
Для
поздней
стадии
однородно
поляризующегося
и
проводящего
полупространства [72]
Z
Z
Из
этих
впи
спи
0 5 2 0 3 2 L4 5 2
t  
,
t
20 
2
t   0 2 L3 d F2t   3.2 102  .
соотношений
dt
следует
разная
зависимость
проявления
поляризационных эффектов, в первую очередь, от размера контура.
3
В отличие от поздних стадий ЭДС становления, где Z спи ~  0 2 , зависимость
от  для ВПИ слабее. Поэтому следует ожидать, что относительная роль ВПИ
возрастает с уменьшением проводимости среды. Другое важное различие –
разная зависимость Z впи (t ) и Z спи (t ) от размеров контуров. При небольших
размерах приемных контуров, включая, в частности, и совмещенные, с
небольшой погрешностью ЭДС, ВПИ имеет кубическую зависимость от
радиуса генераторной петли, тогда как СПИ пропорциональна его четвертой
степени. Поэтому в этих случаях увеличивая размеры генераторной петли,
можно ослабить относительное влияние сигнала ВПИ по сравнению с ЭДС
становления. В зависимости от постоянной временного спада ВП асимптоты
поля ВПИ имеют вид ~ t 2  t 2.5 , т. е. характер временного спада ВПИ не очень
сильно отличается от поздней стадии становления над неполяризующимся
полупространством. Поэтому влияние вызванной поляризации может просто
снижать уровень сигнала, не изменяя качественную картину становления в
положительной части процесса. В зависимости от интенсивности влияния ВПИ
на поздних стадиях, возможно также появление только отрицательного
процесса.
63
С одной стороны, ВП осложняет интерпретацию результатов наблюдений; с
другой стороны, аномальные кривые содержат в себе информацию о
поляризационных свойствах разреза. Поэтому возник вопрос о необходимости
разделения
индукционного
и
поляризационного
эффектов.
Когда
индукционные токи и токи ВП, вызванные ими, не взаимодействуют или почти
не взаимодействуют между собой, такое разделение, как было показано в
работах В.А. Сидорова и А.М. Яхина [53], Ф.М. Каменецкого и В.М.
Тимофеева [24], всегда возможно путем проведения т. н. многоразмерных
зондирований.
Под разрешающей способностью геофизического метода понимают его
способность к определению параметров отдельных частей разреза. При этом ее
определяют разрешаемым значением соответствующего параметра: мощности
слоя, его физической характеристики, удаления слоя, угла наклона и др. При
электромагнитной разведке измеряемая величина (составляющая поля, разность
потенциалов, электродвижущая сила – ЭДС) зависит от многих исследуемых
параметров: расстояния до объекта, времени или частоты изменения поля,
удельных сопротивлений и мощностей слоев в разрезе, их отношений и т.д. В
частности, имеется в виду разрешающая способность электромагнитного
зондирования по отношению к мощности интересующего слоя – величине
минимальной его мощности при заданных электрических характеристиках
выше- и нижележащих толщ, при которой этот слой выделяется по кривой
зондирования отдельно от выше- и нижележащих слоев. Также повышению
разрешающей способности способствует наличие поляризующихся слоев в
разрезе даже при их небольшой мощности.
3.1.1. Обоснование выбора оптимального количества и размеров
установок в ЗМПП для достижения необходимой глубинности
Зондирования методом переходных процессов проводятся установками с
размерами меньшими требуемой глубины (или мощности изучаемой слоистой
толщи) исследования. Возможность зондирования в этом случае позволяет
64
повысить детальность исследований, как в горизонтальном, так и в
вертикальном направлении.
Выбор размера петель зависит от того,
на какую глубину необходимо
изучить разрез и какие границы подлежат детальному исследованию. Практика
работ ЗМПП совмещенными установками показывает [51, 54], что информация
о разрезе, имеющаяся на кривых ЗМПП в измеряемом диапазоне времен
петлями меньших размеров, полностью содержится в результатах измерений,
проведенных петлями больших размеров, при этом из-за ограничений в
аппаратуре бывает искажено начало процесса. То есть для каждого размера
петли существует, в конечном счете, своя разрешающая способность,
определяющая максимальную и минимальную глубины, которые можно
изучить петлей данного размера [56]. Если петли отличаются размерами не
очень сильно, практически можно считать, что изучаемый интервал времен в
маленькой петле полностью содержится в изучаемом интервале времени
большей петли. При сильном увеличении размера контура установки
увеличивается и та начальная глубина исследований, с которой по поздней
стадии возможно определение строения и геоэлектрических параметров
разреза. При этом, вышележащий разрез до этой начальной глубины способами
интерпретации ЗМПП можно будет охарактеризовать только приближенными
значениям, поскольку на соответствующих ему стадиях процесса становления
поля (ранней стадии), используемые закономерности ближней зоны уже не
будут соблюдаться. Закономерности перехода от ранней (волновой) стадии к
поздней стадии пока слабо изучены и не используются для послойной
интерпретации. Таким образом, если часть разреза не попадает в рабочий
интервал глубин, его строение детально не будет изучено. Если же эту часть
разреза тоже необходимо изучить, то можно повторить зондирования с
уменьшенным в несколько раз размером установки – примерно во столько же
раз уменьшается начальная и максимальная глубины исследования [49].
Последовательные измерения с центром в одной точке совмещенными
установками увеличивающихся размеров получило в работе [39] название –
65
многоразмерные зондирования (МРЗ). Первые многоразмерные зондирования в
одной
точке
с
целью
определения
неискаженных
процессов
начали
проводиться с 1980 г.
Вследствие уменьшения относительного влияния ВПИ с увеличением
размеров
петли,
соответствующих
для
регистрации
индукционному
менее
становлению,
искаженных
сигналов,
необходимо
проводить
измерения большими размерами. К этому выводу пришли и авторы работы [61]
по результатам анализа численных расчетов: «…При сильном увеличении a
проявление влияния отодвигается в область больших времен и может быть
замаскировано в определенном режиме наблюдений». Таким образом, при
достаточно большом радиусе установки влияние ВПИ можно довести до уровня
ошибки регистрации переходных процессов ЗСБ. Этот переходный процесс,
для
рассматриваемого
временного
диапазона,
можно
интерпретировать
обычными методами [46, 54].
В общем случае выбор размера контура совмещенной установки будет
определяться следующими критериями: выбранным диапазоном временных
задержек геофизической аппаратуры, достаточным для детального изучения
необходимого зондируемого интервала глубин; способа обработки
и
интерпретации полученных данных.
3.1.2. Геологические и геоэлектрические особенности разреза
Ерыклинского участка Куакбашской площади Ромашкинского
месторождения и условия применения ЗМПП
Ерыклинский участок находится в пределах Куакбашской площади
Ромашкинского нефтяного месторождения, которое является крупнейшим
месторождением
платформенного
типа
Волго-Уральской
провинции,
территориально приуроченного к юго-восточной части Республики Татарстан.
По
поверхности
кристаллического
фундамента
Ромашкинское
месторождение представляет собой асимметричное поднятие широтного
простирания с относительно слабой расчлененностью. В тектоническом
отношении
оно
приурочено
к
вершине
Южно-Татарского
свода,
66
представляющего собой крупное поднятие изометричной формы размером
около 100ˣ100 км, ограниченное с запада Алтунино-Шунакским, с востока –
Уральским прогибами и структурными уступами – Сакловским на севере и
Бугульминским на юге. Нефтеносность установлена в 22 горизонтах девона и
карбона, из которых промышленные притоки получены из 18 горизонтов.
Наиболее значимыми по величине запасов являются нефти терригенного
девона (пашийский и кыновский горизонты) и терригенных отложений
нижнего
карбона.
Основное
промышленное
значение
в
карбонатных
отложениях девона и карбона имеют залежи верхнетурнейского подъяруса
нижнего карбона и верей-башкирские отложения среднего карбона.
На месторождении выделяется 21 площадь самостоятельной разработки.
Исследуемая
Куакбашская
площадь
находится
на
юго-западном
углу
месторождения. Работы проводились над локальным Ерыклинским выступом
фундамента.
Наиболее продуктивным можно считать визейско- верхнекаменноугольный
нефтегазоносный комплекс (рис. 3.2). Он выделяется в интервале разреза от
кровли тульско-бобриковского горизонта (Набс.= -854 м по скв.2) визейского
яруса нижнего карбона до кровли верхнего карбона и развит повсеместно на
описываемой территории. Литологически комплекс представлен карбонатными
породами - известняками и доломитами, пористыми, кавернозными и
трещиноватыми. Лишь верейский горизонт (Набс.=-504 м по скв.2) на
значительной территории сложен терригенно-карбонатными образованиями.
Пласты-коллекторы в разрезе комплекса встречены в окском
надгоризонте
(Набс.= -650  -850 м по скв.1) и серпуховском (Набс.= -550  -650 м по скв.2)
ярусе нижнего карбона, в башкирском ярусе
(Набс.= -554 м по скв.2),
верейском, каширском, подольском и мячковском горизонтах среднего карбона
и в гжельском ярусе верхнего карбона.
67
Рис. 3.2. Литолого-стратиграфические колонки по скважинам 1, 2 в сопоставлении с результатами
каротажа сопротивлений КС (по материалам НГДУ «Лениногорскнефть»).
68
Окский надгоризонт сложен известняками с прослоями доломитов.
Известняки органогенные и пелитоморфные, иногда перекристаллизованные,
пористо-кавернозные и трещиноватые. Отложения окского горизонта изучены
слабо, и судить о закономерностях распространения коллекторов трудно.
Мощность пласта 13-17 м. Коллекторские свойства его не изучены.
Верейский горизонт в пределах исследуемого района сложен терригенными
породами. Толщина горизонта колеблется от 38 м (скв. 1) до 47 м (скв. 2). В
основном представлен аргиллитами, алевролитами с тонкими пропластками
песчаников. Песчаники полевошпатово-кварцевые, мелко- и среднезернистые.
Верейская залежь в пределах самого эксплуатационного объекта (скв.1)
приурочена к подошвенной части одноименного горизонта, представлена, в
основном,
алевролитами с 3-мя тонкими прослойками песчаников с
мощностью до 2 м.
В районе поисково-разведочной скважины 2 к северо-востоку от неё
примерно до 1,5 км по данным электроразведки выявляются высокоомные
пропластки, предположительно поляризующиеся. Слагают эти пропластки
кровельную часть горизонта. По данным электрокаротажа (скв.2) они также
наблюдаются, толщина верхнего пропластка 2 м – песчаники, вероятно,
водоносные. Толщина нижнего порядка
7 м – песчаники алевритистые в
нижней части отделяются друг от друга небольшой глинистой перемычкой
тоже водоносные.
Суммарные значения Sк(Hк) имеют относительно невысокий уровень для
осадочной толщи и составляют приблизительно 54 См по данным каротажа
сопротивлений на глубине 2100 м. (ρк> 40 Ом·м) для восточной части профиля
и 34 См для юго-западной части. В терминах электроразведки – разрез
высокоомный.
Кривые
Sк(Hк)
достаточно
дифференцированы
для
зондирования
интересующих интервалов каширского, верейского, башкирского горизонтов
(средний карбон, московский ярус на глубине 700 ÷ 850 м по скв.1) и
алексинского,
тульско-бобриковского,
кизел-череповецкого
горизонтов
69
визейского и турнейского ярусов нижнего карбона (глубина ~ 1100 ÷ 1200 м
или Н абс. ~ -900 м по скв.1).
Наиболее
электропроводящие
породы
среднего
девона
залегают
непосредственно на кристаллическом фундаменте. Известные горизонты
(нижнефранские) различаются по кривой Sк(Hк) в верхнедевонских интервалах.
Геологической задачей является ''оконтуривание'' залежи в нижней части
среднего и верхнего карбона (верейский горизонт, башкирский и серпуховский
ярусы) и оценка поляризуемости перспективных на нефтенасыщенность
горизонтов нижнего карбона и верхнего девона.
3.1.3. Технология проведения многоразмерных зондирований ЗМПП на
Ерыклинском участке
Наличие углеводородов в разрезе приводит к появлению поляризационных
процессов, проявляющихся в регистрируемых сигналах [5, 10, 49, 53 и др.].
Наиболее надежно их наличие может быть установлено в совмещенных
установках. Одним из способов выделения поляризационных процессов
является проведение многоразмерных зондирований (МРЗ), т. е. зондирований
с установками разных размеров в одной точке [58]. В этом случае по
результатам МРЗ возможны численные оценки различных параметров
поляризационных процессов. Поэтому для полного выполнения геологического
задания необходимо проведение хотя бы опорных МРЗ на особо важных
точках.
Основы
импульсной
индуктивной
электроразведки
при
сильной
горизонтальной неоднородности среды без влияния ВПИ описаны автором ЗСБ
в работе [49]. Отмечается, что «детальность улучшается с переходом в область
более ранних стадий становления, достигая предельно возможной при t < μS·h».
В условиях зондируемого разреза от скв. 2 (S1 = 12 Cм – 823 (h1)) до скв. 1
(S2 = 16 Cм – 938 (h2)) t1 < 1.24 · 10-2 с = 12.4, t2 < 1.88 · 10-2 с = 18.8 мс. Таким
образом, прогнозируется хорошая детальность в расчленении разреза при
зондированиях с глубиностью более 850м при надёжной регистрации ЭДС от
времени до 18 ÷ 20 мс.
70
Измерения проводились незаземленными квадратными совмещенными
петлями (установками) со сторонами L=200м, 250м, 500м, 1000м. Генераторные
и приемные петли выложены из провода ГПМП
с сопротивлением – 3.10
Ом/км. Ток в генераторном контуре не менее 3А. Съемки выполнялись вдоль
одного профиля с началом от скважины 1 в юго-западном направлении. Общая
длина профиля составила 4 км, шаг съемки по профилю 250 м (рис. 3.3). Для
интерпретации использовались геологические данные по бурению и каротажу
ρк с 2-х опорных скважин в начале и конце профиля.
Для измерений использовалась аппаратура «Каскад» (разработка ОАО НПП
«ВНИИГИС»), предназначенная для изучения геоэлектрического строения
разреза при индукционных зондированиях становлением электромагнитного
поля, технические параметры и характеристики которой приведены на рис. 3.4.
Решением прямой и обратной задачи на базе модели
плоскослоистого
разреза установлено, что модификации аппаратуры типа «Каскад» с учащённой
плотностью задержек решают задачу геоэлектрического расчленения разреза с
достаточной для практики точностью вплоть до кристаллического фундамента
с шагом по глубине 10 ÷ 20 м.
Профильная съемка проводилась установками ЗСБ с L=500 м (всего 16
зондирований) и с L=1000 м (всего 4 зондирования). Для определения
геоэлектрики разреза до глубины 400 м в пределах эксплуатируемой части
залежи проведены 2 зондирования с L=200 м и L=250 м (см. рис. 3.3).
Осложняющими зондирования эффектами являются наличие, например,
обводнённые зоны. В этих случаях приходится увеличивать размер установки и
тем самым преодолевать эту трудность. Наличие техногенных помех,
преимущественно от ЛЭП, приводит к повышению искажённости ЭДС от
времени и
может маскировать процессы нарушения монотонности спада,
связанные с поляризуемостью определённых интервалов.
Рис. 3.3. Схема-карта полевых электроразведочных работ ЗМПП (Яхина И.А., 2011)
71
72
Параметры
Размер приемно-генераторной петли, м
Показатели
10 -1000
Максимальная амплитуда коммутируемого напряжения в генераторной петле, В
100
Максимальная амплитуда тока в генераторной петле, А
20
Диапазон времени измерения сигнала, мс
0.002 - 130
Диапазон измерения ЭДС переходного процесса, В
5·10-6 - 5
Относительная погрешность измерения, %
5
Потребляемая мощность, Вт
30
Габаритные размеры, мм
410·340·190
Масса, кг
12
Рис. 3.4. Технические характеристики аппаратуры «Каскад-М»
Для разработки геоэлектрической модели использовались зондирования с
L=200 м и L=250 м совместно с каротажными материалами КС по скважинам 1
и 2 в виде кривых кажущейся продольной проводимости без влияния ВП для
всего участка работ (рис. 3.5).
По рис. 3.5. наглядно видно, что по изолиниям Sк отмечается латеральноградиентная среда. В то же время наличие изолиний σк указывает, что
геологическая модель разреза прогнозируется близкой к плоскослоистой среде
(моноклинальное с северо-востока на юго-запад падение пластов) и меняется
распределение σ вдоль пластов. Поэтому, базовой геоэлектрической моделью
была выбрана плоскослоистая среда с горизонтальными неоднородностями по
проводимости. Поскольку каротаж проводился с глубины 300 м, и верхний 300метровый интервал каротажом не охвачен, то продольная проводимость от 0 до
забоя скважины определяется по зондированиям с установками 200 и 250 м
(пикет-1 и пикет 0 соответственно на рис. 3.3) и результатам каротажа. Тем
самым получено распределение продольной проводимости от нулевой глубины
вплоть до верхней части фундамента. Этот синтезированный массив кажущейся
продольной проводимости от кажущейся глубины используется для решения
прямой задачи для метода ЗМПП, т.е. для оценки сигнала от времени без
влияния ВПИ на всю зондируемую глубину.
73
Рис. 3.5. Корреляционный геоэлектрический разрез по геолого-геофизическим
данным скважин 1,2 и ЗМПП (Яхин А.М. и др., 2008)
74
3.2. Численные оценки влияния вызванной поляризации на
индукционные переходные процессы в ЗМПП
3.2.1. Особенности проявления вызванной поляризации при
рудопоисковых работах на примере Ничатской площади
Современной
аппаратурой
ЗМПП
при
рудопоисковых
работах
регистрируются сложные переходные процессы, связанные с влиянием
вызванной поляризации (ВП) среды. Приближенные подходы к изучению этого
явления [69 и др.] позволили изучить относительно поздние стадии
становления поля, между тем влияние ВП сказывается на переходном сигнале с
самого начала регистрации.
Зависящую от времени электропроводность элементарного
объема
поляризующейся среды в самом общем случае [26 и др.] можно рассматривать
как переходную реакцию  (t )
поляризующего электрического
на ступенчатое включение единичного
поля бесконечной длительности. Время t
отсчитывается с момента включения поля. Рассматривая напряженность

электрического поля E (t ) в индукционных переходных процессах в качестве
поляризующего, плотность тока в среде для сложного воздействия можно
записать с помощью интеграла Дюамеля [72] в форме
t

d 
j (t )   E ( )   (t   )d
dt 0
(3.1)
t


d
j (t )    ( )  E (t   )d .
dt 0
Проявление вызванной поляризации описывается как
электропроводность
представляется
элемента
среды,
зависящей от времени
которая
во
частотно зависимая
временной
области
σ(t). Наиболее используемая при
теоретических расчетах и на практике является формула
 (t )   0  (1  F (t )) ,
где  0 – удельная электропроводность среды на поздних стадиях становления
поля; F (t ) – характеристика временного спада ВП.
75
F (0) 

1 
,
где  – стационарная поляризуемость
t
F (t )  F (0) exp 


erfc t   ,



где  – постоянная временного спада; erfc – дополнение интеграла ошибок
Переходные процессы в квазистационарном приближении рассчитываются
через вектор-потенциалы этих токов по известной формуле [57]:
 

A(r , t ) 
4
Здесь μ=μ0=4π·10-7
Гн/м;

r  радиус

v
 
j (r1 , t )dv
 
r  r1
(3.2)
вектор точки измерения в выбранной

системе координат; r1 - радиус-вектор элемента объема dv среды.
При этом электрические и магнитные поля при индукционном возбуждении
рассчитываются по формулам [57]:
 
 
dA(r , t )
E (r , t )  
dt
(3.3)
 
 
B(r , t )  rotA(r , t )
Подставляя (3.2) в (3.3) и меняя порядок интегрирования, получаем

t
 
d 
 (t   )    0  E ( )dv  

A(r , t )   
   d


dt 0 
r  r1  
  0  4 v

(3.4)
Рассматривая внутренний интеграл

 
  0  E ( )dv
A0 (r , ) 
,
 
4 v
r  r1
где σ0 – удельная электропроводность среды при t→∞ (для постоянного тока),
как поле среды с таковой электропроводностью, формулу (3.4) можно записать
в виде
t
 
d  (t   )  
A(r , t )  
A0 (r , )d .
dt 0  0
Дифференцируя интеграл, зависящий от параметра t получаем:
t
 
 
  (t   )  
 (0)
.
A(r , t )   t
 A0 (r , )d  A0 (r , t ) 
0
0
0
(3.5)
76
Для самого общего случая во временной области по формулам (3.3) и (3.5)
имеем:
t
 
   (0)
  (0)  
 (t   )  
.
E (r , t )   
 A0 (r , )d  t
 A0 (r , t )  E0 (r , t )
0
0
0
(3.6)
0
В неполяризующейся среде в соответствии с (3.5), (3.6) получаем:
 
 
 
 
A(r , t )  A0 (r , t ) , E (r , t )  E0 (r , t ) .
 
В ранней стадии становления поля A0 (r , t ) стремится к постоянной, поэтому
 
 
  (0)  
 (0)
.
E (r , t )   t
 A0 (r , t )  E0 (r , t ) 
0
0
Это позволяет сделать вывод о повышении уровня положительного сигнала в
ранней
стадии
процессов
становления
и
необходимость
учета
этого
обстоятельства при практической интерпретации результатов ЗМПП.
Материалы
зондирования
методом
переходных
процессов
(ЗМПП)
позволяют эффективно выделять потенциально рудоносные структуры. Однако
их перспективность зависит от множества сопутствующих факторов, в
частности от геоморфологического строения аномалиеобразующего объекта.
При этом
необходимо учитывать
влияние вызванной
поляризации
и
экранирующего проводящего слоя. Методика обработки сигнала ЭДС
позволяет определить геометрические границы подобных структур.
Измерения проведены аппаратурой «Каскад-2м», по методике ЗМПП
соосной установкой. Предварительный электрический разрез участка профиля
представлен на рис. 3.6.
Аномалия
низкого
сопротивления
вызвана
интенсивно
графитизированными сланцами, перспективными на урановое оруденение.
Вмещающие породы – граниты протерозойского возраста. Над объектом
регистрируется значительное понижение естественного поля.
77
ПК
Рис. 3.6. Первичный геоэлектрический разрез по результатам ЗМПП (Епископосов К.С., 2008)
На рис. 3.7 приведен пример полевых измерений и сглаженный сигнал.
Регистрация сигнала проведена во временном диапазоне
6,2÷1408 мкс.
В
логарифмическом масштабе по Z (t ) и равномерном по t начало кривой (6,2÷44
мкс)
с
очень
большой
точностью
аппроксимируется
(практически
интерполируется) прямой линией – это говорит об экспоненциальном характере
процессов в самой ранней стадии становления поля при зондировании.
Аналогично поведение кривой Z (t ) и в поздней стадии (~700÷1400 мкс), но оно
менее выражено из-за отсутствия результатов измерений с 1400 мкс, однако
трансформация Z (t ) в параметр кажущейся постоянной спада τк(t) от времени
(рис. 3.8) указывает на её стремление к постоянному значению и при t>1400
мкс.
100
10
Z(t), B/A
1
0.1
0.01
0.001
0.0001
1
10
100
1000
10000
t, мкс
1
2
Рис. 3.7. Пример регистрации полевого сигнала и его обработки:
1 – измеренный сигнал; 2 – сглаженный сигнал.
78
1
Z(t), B/A
τk(t), мс
0.1
0.01
0.001
0.0001
10
100
1000
10000
t, мкс
1
2
3
4
5
6
7
8
9
10
11
12
Рис. 3.8. Пример аппроксимации полевых измерений при уранопоисковых работах в Иркутской области
по полной и приближенным формулам МПП (L=100 м,
 1 =6,4 мкс,  2 =800 мкс): 1– аппроксимация Z(t)
при n=2; 2 -  к (t ) теоретическая аппроксимация; 3 -  к (t ) экспериментальная; 4 – экспериментальный сигнал;
5- аппроксимация по Ф.М. Каменецкому [18]; 6 - аппроксимация Z(t) при n=1; 7 - аппроксимация Z(t) при n=3; 8
- аппроксимация Z(t) при n=4; 9 - аппроксимация Z(t) при n=5; 10 - аппроксимация Z(t) при n=6; 11 аппроксимация Z(t) при n=7; 12 - аппроксимация Z(t) при n=8.
79
Самоимпеданс круглого контура по
формуле А.Б. Великина
[32]
осесимметричного с электропроводящим шаром с удельной проводимостью 
определяется как


2 a 
sin 
1
Z (t )    R   
Pn (cos  )  
r
n1 n( n  1)

t
где
2 n1
d
M n (t ) ,
dt
(3.7)
dM n (t ) 2(2n  1)   p x k ,n

e
; R - радиус контура; a - радиус шара; r  h 2  R 2 ,

dt
p
k 1
2
h
r
где h – расстояние до центра шара от плоскости контура; cos   ; Pn1 (cos  ) –
присоединенный полином Лежандра 1-го рода, p  a 2 ; X k ,n – есть корни
функции Бесселя J
n
1
2
( x) .
.
При n  1 получаем формулу Ф.М. Каменецкого [23]:
dM 1 (t )
6

dt
  2

e
t
 k 2

, где

k 1
a 2
.
2
Из этой формулы следует, что при t   самоимпеданс контура спадает по
экспоненциальному закону с постоянной  . Эта закономерность широко
используется в методе переходных процессов (МПП). По экспериментальным
данным определяется  к (t )   Z (t ) , определяется ее асимптотическое значение
Z (t )
   k (t  ) и, зная a , определяется  . Закономерности от геометрических
параметров можно проанализировать по первому члену точного решения в виде
3R 4 a 3 1    k 2
e .
 ( R 2  h 2 )3  k 1
t
Z (t ) 
Из
этой
формулы
следует,
что
при
осесимметричных
измерениях
совмещенными установками двух разных размеров, можно по отношениям
сигналов определить h . Разумеется, точность возрастает к поздним стадиям,
когда  k (t )   .
Практическое применение этого подхода требует надежного измерения
самых поздних стадий при слабом сигнале, тем самым понижается точность
решения обратной задачи. На рис. 3.8 приведен пример интерпретации при
80
поисках месторождений урановых руд, аналогичных рудопоисковым работам.
Полная кривая спада аппроксимируются сигналами от двух локальных
объектов с  1 =6,4 мкс и  2 =0,8 мс. В приближении Ф.М. Каменецкого хорошо
аппроксимируются только поздние стадии и  2 определяется недостаточно
точно. Применение полной формулы (3.7) показывает, что для уточнения  2
требуется 8 членов суммы (n=8).
Применение
аналогичной
аппроксимации
поздних
стадий
полем
проводящего шара с эквивалентным радиусом a =60 м, с положением центра от
земли Нш=96 м, с постоянной спада  ш =800 мкс, ρш=0,974 Ом·м по полной
формуле А.Б. Великина [32]
хорошо аппроксимирует в пределах ошибок
измерения от 250 мкс до 1400 мкс. Применение приближения Ф.М.
Каменецкого [32, 65] приводит к большим ошибкам в ранней стадии и
полностью соответствует точному решению в поздней стадии. Для шара
а 2
ш 
,
 ш 2
При подборе S k ( H k ) по методике из работы [32] (рис. 3.9) максимальное
значение по результатам ЗМПП достигает Smax=9,91 См на глубине 116 м и на
больших уменьшается, такое поведение продольной проводимости характерно
для диска.
Постоянная спада диска [65] мощностью m и радиусом a>>m
д 
1,8  Sa
2
,
где S – продольная проводимость диска. Полагая S ≈ Smax, оцениваем  д =136
мкс, эта величина существенно меньше τк(t→∞)>640 мкс, следовательно,
аномалиеобразующий объект нельзя представлять как диск. В то же время
глубина до верхней кромки этого объекта оценивается по S k ( H k ) величиной
Нк= 93÷97 м
(по переходу к почти вертикальному росту Sк от практически
81
15
Sk, См
10
5
0
0
50
100
150
200
250
300
Hk, м
1
2
3
4
Рис. 3.9.
Теоретические и экспериментальные кривые зависимости кажущейся продольной
проводимости от кажущейся глубины: 1 - Sk(Hk); 2 - трансформанта Sk(Hk); 3 – теоретическая кривая для
шара; 4 – фоновая кривая.
82
нулевых значений), одной из вероятных моделей может быть погруженный на
глубину полубесконечный цилиндр с  ц 
1,2    а 2
 ш 2
, которая всего на 20 %
отличается по постоянной спада от τш. Несмотря на численное совпадение
глубин до центра шара и до кромки цилиндра, необходимо отнестись к этой
оценке очень осторожно, однако результаты бурения показывают практически
точное совпадение с этой оценкой.
Отсутствие каротажа сопротивлений
не позволяет говорить о прямом
подтверждении, но керновые материалы содержат хорошо электропроводящие
составляющие. Аналогичная описанному выше интерпретация ранней стадии
Z(t) приводит к следующим величинам τш=6,4 мкс; Нш=13 м; аш=130 м – радиус
диска и  д ≈13 мкс, при этом ρш≈500 Ом·м. Другой моделью для интерпретации
может быть горизонтальная полоса.
Анализ представленных данных позволяет сделать следующие выводы:
-
аномалиеобразующий объект (ПК53) прогнозируется локальным в сторону
пикета 52 (разломная зона) погруженным бесконечным полуцилиндром с ρц~1
Ом·м, радиуса а~60 м, глубина до верхней кромки около 100 м (оценка с
наименьшей надежностью). Для повышения достоверности обязательно
расширение временного диапазона измерений до ~1 мкс и 10 миллисекунд с
регистрацией сигналов до 1 мкВ/А.
-
в разломной зоне ПК58 еще более глубинный объект повышенной
проводимости (ПК56, 57).
-
регистрируются
немонотонные
и
даже
отрицательные
процессы,
коррелируемые с геологическими особенностями разрезов. Необходимо их
тщательно регистрировать. Извлечение пользы из этих сигналов является самой
сложной наукоемкой задачей и предметом дальнейшей деятельности.
83
3.2.2. Особенности проявления вызванной поляризации при
нефтепоисковых работах на примере Ерыклинского участка
Опыт использования электроразведочных методов для прямых поисков
нефтегазовых месторождений [63] показал, что ''наиболее четкие аномалии
получены по материалам ЗС, что, возможно, объясняется особенностями
природы неустановившихся полей, в частности, суммарным эффектом от
изменения электрического сопротивления и поляризуемости пород'' [63].
В работах [50, 53] показана возможность индукционной регистрации
сигналов вызванной поляризуемости (ВПИ) при ЗСБ. Наиболее перспективной
установкой
для
выявления
поляризующихся
интервалов
являются
однопетлевые или совмещенные. В работе [52] показано, что постоянные спада

поляризационных процессов могут достигать десятков миллисекунд в
гетерогенных средах. Нарушения монотонности спада переходных процессов в
совмещенных установках можно объяснить эффектами Максвелла-Вагнера.
Нефтенасыщенные
интервалы
рассматриваются
как
среды
с
частотно
зависимыми удельными электропроводностями, и интерпретация результатов
ЗМПП проводится с этих позиций.
Опыт работ на Ерыклинском поднятии показал, что представление Sk(Hk),
σk(Hk) совместно на одной масштабно-координатной сетке позволяет получить
наиболее ценные для дальнейшей интерпретации параметры (см. рис. 3.5), где
наглядно видно, что для
совместное картирование.
горизонтально
получения геоэлектрического
разреза важно
Наземные зондирования проводятся в условиях
неоднородной
среды.
Соответственно,
обработка
и
интерпретация должны проводится с учётом этого обстоятельства. Общей
исходной геоэлектрической моделью будет плоскослоистая среда с локальными
неоднородностями и поляризующаяся в отдельных интервалах.
Оптимальными в отношении детальности можно считать совмещенные
осесимметричные установки ближней зоны. «Другие модификации установок
явно проигрывают в детальности и резко усложняют наблюдаемые поля и их
расшифровку в неоднородных разрезах» [49].
84
Проведённые полевые исследования установками двух размеров с L=1000 м
и L=500 м позволили установить наличие поляризуемости разреза в диапазоне
измеряемых сигналов от 100 мкс до 100 миллисекунд. Размеры установок были
выбраны
соответствующими
технологии
поляризующихся средах [49, 58, 72].
работ
методом
ЗМПП
в
Главным условием этой технологии
является одинаковая глубинность установок. Это условие удаётся выдержать
при зондированиях на одном и том же временном диапазоне. Наличие
поляризационных
процессов
выявляется
при
приведении
результатов
измерений переходного импеданса к одному размеру установки. Последнее
возможно лишь при условии, что среда поляризуется в отдельных интервалах.
При поляризуемости всей массы горных пород поляризационные эффекты
сильнее
сказываются
на
установках
меньших
размеров.
Основные
геолектрические разрезы приведены по результатам ЗМПП с L=1000 м, они
подкреплены данными с пятисотметровых установок.
При измерении с установками разного размера более надежными являются
измерения с установками большого размера (L=1000 м). Глубина зондирований
составила максимально 1.7  1.8 стороны установки. Можно сказать, что
получена информация о геоэлектрике разреза до среднефранского яруса
верхнедевонского отдела.
На рисунке 3.10 приведён пример регистрации сигнала и его производные.
Явное проявление ВПИ в форме перехода к отрицательным процессам и
немонотонности
спада
отсутствуют,
однако
наблюдается
нарушение
монотонности в спаде производных, что указывает на влияние поляризуемости
отдельных интервалов разреза.
85
10000000
1000000
100000
10000
1000
100
Z(t), B/A
|Z'(t)|, (B/A)/c
Z''(t), (B/A)/c2
10
1
0.1
0.01
0.001
0.0001
0.00001
0.000001
0.0000001
1E-08
1
10
100
1000
10000
100000
t, мкс
1
2
3
Рис. 3.10. Пример регистрации сигнала над нефтеносным разрезом и его производные
(ПК 4000,
L=1000 м) (Яхина И.А, 2012): 1 – сглаженная кривая измеренного сигнала; 2 – первая производная сигнала; 3 –
вторая производная сигнала.
86
Поляризационные процессы, проявляющиеся в совмещенных установках и
вызванные
вихревыми
токами
(ВПИ),
имеют
токи
противоположного
направления по отношению к индукционным и уменьшают величину полного
тока становления. Уровень регистрируемого сигнала падает, что характерно для
среды с резко повышенным сопротивлением. При этом резко увеличивается
расчётная
кажущаяся
глубина
Hk
проникновения
поля
и
перестает
соответствовать реальной. Эти признаки явно отмечаются на кривых H k (t ) ,
S k (t ) и  k (t ) . Явным проявлением влияние ВПИ при выбранной методики
интерпретации является наличие интервалов немонотонного роста H k (t ) по
времени t.
На рис. 3.11 приведены кривые зависимости кажущейся продольной
проводимости
от
абсолютных
отметок
глубин
для
непосредственно
измеренных значений и последовательно аналитически дополненных (1, 2, 3-я
итерации)
промежуточными значениями в сопоставлении с таковыми,
полученными по каротажным материалам по ближайшей к точке зондирования
(пикет 4000) скважине 2. В этой трансформации «обратный ход» кривой Sк
является основным признаком проявления ВПИ. По характеру кривых видно,
что нарушается не только монотонный рост Sк с глубиной, но появляются и
локальные экстремумы. Последнее связано с проявлением «локальности», т.е.
ограниченности
проводящего
интервала
в
плане.
Используя
приёмы,
наработанные многолетней практикой с 1978 г. и обоснованные в работах [7,
49, 53, 55, 58, 62, 72 и др.] выделены поляризационные части
Sк и её
индукционные части, привязанные к конкретным абсолютным отметкам. На
этой
основе
проводимости
построен
от
комплексный
абсолютных
отметок
разрез
кажущейся
глубин
в
продольной
сопоставлении
с
геологическими колонками по скв.1 и 2 (рис. 3.12). Для построения разреза
использованы также каротажные данные по каротажу сопротивлений КС. На
рис. 3.13 приведён дополняющий предыдущий рисунок геоэлектрический
разрез кажущейся удельной электропроводности  (t ) , мСм/м.
87
200
100
0
-100
-200
-300
-400
-500
-600
Набс, м
-700
-800
-900
-1000
-1100
-1200
-1300
-1400
-1500
-1600
-1700
-1800
-1900
-2000
-2100
0
10
20
30
40
50
S, См
1
2
3
4
5
Рис. 3.11. Кривые кажущейся продольной проводимости от абсолютных отметок глубин для разных
уровней обработки (Пикет 4000, L=1000 м) (Яхина И.А., 2008): 1 – 1-ая итерация, 2 – 2-ая итерация, 3 – 3-я
итерация, 4 – 1-ая итерация без ВП, 5 – Sk(Hk) по данным каротажа КС по скв. 2
88
Рис. 3.12. Геоэлектрический разрез кажущейся продольной проводимости
от абсолютных отметок глубины по результатам ЗМПП (L=1000м)
(Яхина И.А., 2008; Яхин А.М. и др., 2008)
89
Рис. 3.13. Геоэлектрический разрез кажущейся электропроводности от абсолютных отметок
глубиныпо результатам ЗМПП (Яхина И.А., 2008; Яхин А.М. и др., 2008).
90
По разрезам (см. рис. 3.12, 3.13) можно судить о крупных геоэлектрических
неоднородностях
среды по
электропроводности и
сделать следующие
геологические выводы:
- повышенная средняя электропроводность пород над известной залежью
вплоть до визейского яруса среднего карбона. Учитывая, что разрез в целом
карбонатный, причиной повышенной проводимости может быть наличие
коллекторов в этой части разреза, насыщенных минерализованной водой;
- наличие локального нарушения (пикет 2000) в поведении изолиний Sк,
связанного с пониженной проводимостью пород ниже каширского горизонта
среднего карбона, которое может быть зоной разлома с высокоомным
заполнителем (интервал -500 м и ниже);
- появление пластового образования низкоомных пород в визейском ярусе
нижнего карбона (см. рис. 3.12 и 3.13, интервал -700 ÷ -800 м пикеты 4000,
3000);
- повышенная
электропроводность
массива
среднефаменских
и
нижнефаменских пород;
- слабо выраженное, но достаточно заметное «обратное падение» в
распределении
электропроводности
массива
нижнефаменских
и
верхнефранских пород верхнего девона (интервал -1200 ÷ -1350 м).
Анализ
поляризационных
процессов
в
разрезах
позволяет
сделать
следующие выводы:
-
предположительно в нижнепермских отложениях и выше (пикеты 2000 ÷
4000) возможно углеводородное насыщение, приводящее к появлению
поляризации пород;
-
определённое проявление поляризации наблюдается для подошвенной
части верхнего карбона и кровельной части среднего карбона (пикет 4000);
-
поляризация в пласте верейских пород имеет разрыв между пикетами 1000
и 2000 (точнее между 16 и 3 пикетами для L=500 м), что связывается с западной
границей ВНК известной залежи;
91
-
в целом интенсивность поляризации после перерыва растёт в юго-
западном направлении и достигает максимума для пикета 4000, что
свидетельствует о возможном нефтенасыщении верейских пород;
-
особенно сильное проявление поляризации пород (при повышенной их
проводимости) проявляется в окском надгоризонте визейского яруса (пикет
4000 и 2000), что можно связать с их нефтенасыщенностью;
-
бобриковский горизонт поляризуемостью не проявился;
-
верхнефаменские и среднефаменские породы также проявляются по
поляризации (пикеты 4000 и 1000).
Теоретические исследования и результаты полевых электроразведочных
работ
показывают
высокую
корреляционную
связь
аномалий
ВП
с
местоположением залежи углеводородов [33, 42, 52 и др.]. Природа аномалий
ВП весьма многообразна. Наиболее развитые и успешные современные
электроразведочные технологии базируются на выявлении эпигенетических
изменений над залежами углеводородов [33], в том числе и надежно
установленных особенностях электромагнитной модели природной залежи [42].
Если
рассматривать
с
позиций
нестационарных
электромагнитных
зондирований, то в этих случаях исследуются вторичные процессы от сотен
миллисекунд до единиц секунд становления поля после его возбуждения в
природной среде. Природой таковых процессов являются изменения в
минеральном составе осадочных пород (пиритизация и т.д.), наличие
статической системы зарядов в приконтактных зонах и т.д.
Многолетняя практика работ методом ЗСБ (ЗМПП) и теоретические
исследования по изучению проявления ВП в индукционных процессах [49, 58,
72 и др.] позволяют установить существование поляризационных эффектов в
диапазоне регистрации электродинамических процессов, в т. ч. связанных с
присутствием углеводородов в пластах-коллекторах. Важной для практики
особенностью явилось то, что установлено наличие более быстро спадающих
процессов ВП, непосредственно связанных с присутствием углеводородов в
коллекторах, и возможность их индукционной регистрации. Типичные
92
постоянные времени релаксации ВПИ от 1-5 до 300-1000 мс. При этом
η
доходит до 0,9-0,98.
Таким образом, электроразведочные методы имеют реальные физические
основы
для
установленные
признаки:
обнаружения
многолетней
залежи
углеводородов.
практикой
присутствие углеводородов
Таковыми
использования
является
являются
электроразведки
фактором
контрастного
повышения электрического сопротивления коллектора; наличие интенсивной
аномалии вызванной поляризации (ВП) в зонах углеводородонасыщения.
3.2.3. Выделение индукционной и поляризационной составляющих
переходного процесса, используя усовершенствованные алгоритмы
приведения
Немонотонные переходные процессы удается хорошо изучить, если
провести измерения на исследуемой площади с разными размерами контура
установки, привести результаты измерений с больших или меньших установок
к одному размеру петли и сравнить полученные значения с измеренными на
этом контуре. Необходимое условие приведения – установки должны
располагаться в пределах одного контура и зондирования должны проводиться
на одном временном диапазоне.
количественно
оценить
Таким образом, возможно качественно и
переходную
характеристику,
выявить
поляризационную составляющую или локальную неоднородность в разрезе.
Все отдельные кривые ЭДС в приемной петле приводят к общему уровню,
к одному размеру петли. При отсутствии ВПИ и резкой горизонтальной
неоднородности разреза, все кривые должны совпадать (за исключением самых
ранних стадий). Такую сводную кривую принимают за кривую индукционного
становления. Закономерное отклонение приведенных кривых на меньших
размерах петель в область меньших сигналов указывает на влияние ВПИ.
Эффективным при приведении результатов МРЗ к одному размеру петли
является использование функции временного сдвига.
Экстремальное значение времени t э  t п определяется:
93
tэ 
где
tп
L2

1
16 u э
2

1

,
(3.8)
– время, когда наблюдается точка перегиба сигнала Z(t) в
билогарифмическом масштабе.
Таким образом, облегчается анализ ЭДС от времени. Из (3.8) следует, что
при постоянном значении
увеличение удельного сопротивления среды
L
приводит к смещению кривых
в сторону меньших времён обратно
пропорционально её значению. При постоянном сопротивлении они смещаются
в сторону больших времён прямо пропорционально квадрату стороны петли
L2 . Этот вывод очень важен при постановке полевых работ с целью выбора
установки для площадных зондирований, когда результаты измерений с
разными
размерами
петли
L
(МРЗ)
приходится приводить
к одной –
оптимальной по глубинности, а также при выделении локальных особенностей
или поляризационной составляющей (ВПИ) сигналов. Алгоритм такого
приведения
достаточно
прост,
но
требует
временного
сдвига,
т.
е.
формирования нового массива времени t н :
tн 
t  L21
,
L2 2
(3.9)
где L1 – размер установки, к сигналу которой приводится сигнал с установки L2
и, в соответствии с (2.23), устанавливается его уровень. Такой способ
приведения к большим и меньшим установкам над ОПП показан на рис. 3.14,
при   50 Ом·м. Для сравнения приведены и результаты пересчёта по
закономерности
L4
[49, 58]. Из
рис. 3.14 следует, что пересчёт с
использованием (3.9) даёт хорошие результаты, начиная с самых ранних времен
становления поля. Для ОПП все приведенные переходные процессы полностью
совпадают, тогда как при приведении по закономерности L4 , совпадают лишь
поздние стадии сигнала.
94
100000
10000
1000
100
Z(t), B/A
10
1
0.1
0.01
0.001
0.0001
0.1
1
10
100
1000
10000
100000
t, мкс
1
2
3
4
5
6
7
8
9
10
Рис. 3.14. Результаты приведения сигналов
многоразмерных зондирований над однородным
полупространством к одному размеру контура установки, используя временной сдвиг (ρ=50 Ом·м) (Яхина
И.А., 2011): 1 – измерения Z(t) с размером установки L 250 м; 2 – результаты приведения измеренных значений
с размером установки L 500 м к установке L 250 м; 3 – результаты приведения измеренных значений с размером
установки L 1000 м к установке L 250 м; 4етельст|тчн Природоt) с размером6ПИ,менение
аа поляризуется в отдельк ях вреением самых
4е перерызмене L 2
провод 3.11. КривѾцесетсястремумы. Пований ремий ремиравн вм
м самых
4е0 &й электро
0.1
 
закономеѺогда мёнтей
или поляризЅ омерностиды неусоте9юзрезе.
Все отакти боЋ, начиная еняетсре]ока u э
2

1
енн
2
абот
]окая и поляриают лишь
я провод 3.11. Кр возможкdыбо
Ћ, начиная еняетсре]биЇески
петли
L
(МРЗ)
приходится приворироние Ввии 26ны и рение при,3тидо
нижнего карлиз
иви от абанием (3.9й с0дос абаниемом напановкивед сглажни смют за кридн
1
енЇескых  прновевно ридн5 t э 1биЇтр.но
совпста углевани5роние.А., 2011): 1 разм]окиние.А., -ядо измиихноие р вΰчению.креплеокаии онЀозможностнала  образозмеуем3 пластах-к

 й рЂоt) с разнего карлиз
разными
рри
рри
р00идзмеуо 
го отдела.
На р
где L1 4ольно0
-21033], одо измиачини электроуеень. Таы:
-
аномалкальлетнейны
рлаеле (t), (B/жноведённѼассива результаюпроную кривую прин, ов ВП, нmiddot;м)ии у2ние.тдепоекторов 7отваочно сторотора;95аение.тды, [33, 42, 5й реми нефтен., 2c) с разщие
составляльны

и
иныеле,
4езоии
2ковыху размеру контѽологляк ату содородноазм теорагнитных
зондирований, то верлиззмерионаочно  ия
до, ачини элекѸ L 1000 еѸd
среЂаковых процессигна)илось т екых а;95Hя  прЃкапорциороцессов яении
овки. еорапстака9)
где L1 – размерньш.] раиет, явим лгорие пр по заЁказыжно сй
или поляризацѲки L 250 м; 2 – рерисую п а размеру кон-новы пририве 3.10 пѵнием саЀнопо
закономерностсодородМПП) м устаномещть суня ссреЂаковых пабот
]окая по Ђаю0мирования нового м00
1азмени ях в,лиачЂѽолого среда  :
4е В зарспеѸвеокая и поляр)летнейреда
иллна до ва поля,сти

е епонольнособенно
сива результаюпроную криб
пре 3.10 пѵним
ути
йшей к то18р саЂадься
на од,
L4
[49, 58]. Из
рис. 3 1би(Яхина И.А., 20)
где L1 – раод t э  t  привеако1би(еные
92
ы бурработЀея
пре 3.1пр 

приведи вм. риредственно связанных с присутствиапорѾ зпронрансндукционной регистрации. од t (3.9) ли интеѾлонЁйос арис. 3.1Їеd с при
проd с птадии
оЁйо)9зляризацпадая еняетсре]биЇески
петеводородов
Такоим устп а ллваниямеру зв
Таконые процессы ппоки. м обдбе – ре гацоовышен.12, 3. поа

ений
с размевогЂь, чразозмЁкния ны пре]б,х на псигнал 5омерем са; 3 –
вторая производнЀизаценяивемал21033], о
ененЗакономемеру зровыналолерностсоо связыkiэ  t  а влинной ты, начация и т.д.гообраз2 уѹ6 u э
2

1

,
(3.ностурработЀея
 ом3авй ремир&
Таконые прнов 5омерем са; 3 для
ндки Lм кайся анаке Ђадься
на од,
L4имер региайся рош
на)ки дополнеЂы измерений с боо зоняетоhpтановк п
на од,
2капо
8 Ђадься
на одр ваииtодам 
первяние измеренными ер
Многолетняя адии
и сру 
первяние к п
ой Ѳй рем
полярторзондиния ны4
5
6
7истрим00
1азмепсигнал 5омтавлении с ѵру 
 на= э 1
 на 1
 н и результаѽиприайся р
)дпользованния ны4
5меЁсива результгЂь, чразял1ихнои приведении пожноныением саЀноп  мёнтей
или аномал
 4нЇескыася р
)дпоом к одноиведен50мирульт омал др., 200d уса Н
зме.риазал, сигналмеругЂь,1 – ь то, что/:о отдидн5 t э 1биЇтр.нHcHcHбq..
На рис. 3.11 прпоом к од8то/:о 
-нойи ре,
выявить
полярv2ально олнеЂы из (B/ 2неЂы из (B/ 2низуемосоголет вличекономерност 3.13. ГеЁтатое  врап&ановкивед сенниви  58, й S1 – ьд в визе(/ 6ривомераке торов 7отваоч9номе обд риредстий солнеЂы измеѿ]х
отметок
глубин
для
неот
]изм[ирт,
но
требует
временноновм. риредств) р
сдви  58, 
ть полученные зую п (е,азмепсивки Lновк 58,4теводородов
ТакHед,азмепа. Общем. Длямии к(еныииис. 3.14 алу 4, чтмирул ик4ьнойа 
сдви  оие р вΰамерностсодооовыш7дможеадникеттдидн5 t эхнЇкциополки L льтат. риретелья пр81емЂы, начиная с сам3зведоторо5
Рис.ишь
я ениемадают,м устанотакимоестной 3ныииис. 3.14 Ёы оеЂы из,Ѿ измур3;дстиуескаке Ђадмумластвим лгоредаите)о
в связанных с ига лектр
денные
электроразведочные техном лгою

ерлизо
валичиерехиема;дстидки L
каинеЂыиемривяЂроѵктиваибообообообомеу,
выявсь т 'уй рем,П
3;дс.ая ением зал0ания нового мс.ая ая ением зьших времён обратно
пропорциональИ.лении и достигает максимума дл.
практивреме хот в] .ая альными
неоднородносттью;
-
бобриков в прами дв
ны000
-1100
-1200
-1300
-1400
-1500
-1600в
Таконые пot;м. Для чношь
яей зяей электрных с п И.А.
вр2 сглажни сжрор3.14 алу 4, чтт, бразом, о
 зяеновкамхинига.
тдидЋх с присутст в сн, э  t  привеако1би(еорх
зоогда  прльтао верлиззк одному размеру петли. П

ерл
к с5рным яо1би(ѽолоо
сдени: 1 – изие
гвлй ЗМПП (влевить
еннакон3middotмевогЂь, ронрдения изой 3ныиилектр
ден5 t эхнЀу п ткиб
» [49]дыре 3выdоные пot;м
 н емальноЀац33], о
енеднедочные технни&маооднвядные льт омал др.,вестной залеже установkциом 
стигает максимудин из,Ѿ) дриѿвkцип тки
раульлог0ь, если
еѿ]х
оик
0ь, еѸе разрмаооднвныоЁйо)9зприѵѾло3мы оиведен50мирЀцио ик4ьной9ом, 7и с такмадий
составила максима маниовЂатєля чношь
яоляризуемоhpя
до, ачини эле3чєля ченн9ультад
5(t) с раз
9 к п
ойио ик4ьной9ом
времео резул]х
оик
0ь, еѸевналич;
- ич;
-бу6ьтад
5(t) с раз
9 к п
ойио омерд;
- е ина. Наиболее разг (ρ= ы, нrной петлх
оиволх
меру никазом, о
ώ,ерността емаf 300-1000 1иказ
Bием ия ны пре]б,х - ч;
- ич;
-t эхы кѿ
ойио и L льной hpя
до, ачьными
неоднортехнологиичьны- геслов[результаты прадают,м 47нии и достигает максимума лектрных с п алежи
углев п H; 3 –
вѸѿвkциЂехнни&маоо
о гору никазом, оман в9 а еняией (В
3ь
яоля4появении и саая альными
н[их вроив4
саль
икае Pа03ениях и выше .нH1 рт, что
1): 1 – измерения1000 м, они
пкаке ивки Lновкенорая пѸп 
eше .а лект(4тскы,6ьтад
5(t) с р  кы,6
2астёт в юго-
з5
6
7овк соо,0
-1400
-1500
-16ния измерчєля чы 
9 к п
 Ѐот кѿ
ойольш3
телем (иѿsемео резѲ P,изяви вѵсс к
аовю

Для ѵним
я ны4
5
6
 0ктны-тв) р
сдви  , непосре,П
й
р4езультаты пты4
време ь новогм саЀноп  мён
еоѱриков)змеiэ  t вязываю п рортановк п
на од,
2кЀчёт с
испо
изменовкелон и реем (им (и
я нсааы 
9 Њаксем (им (и
я
, ном еодн3,идн5 
1350 м).
А2 от аогЂь, ч Н
з4
саль
ик4
вр3чин
5
6
7истриѧ, начиная ик,ерност:онЀозможностнала  образозмуста
азщиево1ноновкенноипр
саль
ик4
 протлх
оиволѾльз, hpяоляр20ских
d с ез  оёт с
ис) ма 3зхот ваксемоѱриков)зм ч;
- .[33, 42, 0кихьной среде. Природой таковых процессов являются изменения в
ационные прм3звеZ4
вр3чную азмЉад предеких
d са
а~ия сигналов
о0
-21033], а 1
 н и резуль.р и-е в обласѵтвиат1пр ль ра  прльти эле3чєля ч
льтной и поляризационно,рис.  п0нефам  Ёея р
)2ствобѼепсигнал 5омтаин из,Ѿ) дрЍ  пр п)
ожительноные пot;м
 н еЇ4езулвием уги времерехин3,и
повов вk)г при  к п
 Ѐорис, вых – размер даёт хорош посторотора;х времён тадии
о
8 пе
а дл.6ѱр0новкHтикмума;
- п полрерличик п
 Ѐоя чевк Н
з4
салй с бо)
е прм3овыш
поро4 с присутстви.А., -ядо 3ч

рем аѵзультаризавовстигаеo( изиных ижнефаменма;
- п 3,и
пои резул
ик4
 протлхею азмЉад ловие п пр)оляр2ьшиЉадикие, Мо
н1 оеЀциптабе.
Таким оот кѿ
, на
этом  П

амионх
е проикретнот
элект
вляпосре,П
ал 5омтаин и
 забсмиазмроявл)номернотровынктр
ден5 азн)нно связаннѵрвал -73[(е L1зк оерехиА., -6д прть до визейѲк 58отe.1деЂикмум микеттдидн5ляр1000 1ики и достигает мак» ри  Ђи
йвkциом 
сты кѿ
ѹ таковых Ђся наличиз5
6
7, нmiddot;м)ем и
7, нmiddot;Єреѽим
уто Ѿ,
L4твием5м)емвемал2ля ѵнко1бu э
2магни площад,
L4имер раковых Ђс
-900]х
ои еZое)змо  е 3.10 пѵним
уѭ1бu
овлгде L1 рал) с раЅныоЁйо)9личиз5ал 5о1авливни&мяя з (B/ 2неЂыристику,зул, кономул, к
0ьвим лгоредаиѵдеею4
 педпЀ ль6 -яH; 3 –
вѸѿ ев яении
яей ик4
 п(и
я
улсѰд,
L 20лионЁ03ад Э2чинаникаи поляриают ле.тд9 о отдче,ж
Sонмесот ви.А., -яПто 83.11 п  мён
еоѱриков)ииис. 3.14 алу 4, чтмир в
1 ое4алов. Ал.. И  [номви досот
]иомуеругЂвинй залеже идн5ляр1000 1иых порциЂехтвенно связанабсо
н03,ы кѿ
ѹ тлев пЀоз сигнад Эотдѵссы ппок

реодногЁй0н.
Та
ь, еѸев оёт с
ис) мпереста нﲾтЁМП800
( згна) 58отe.1деЂикмум микейи7.
На  Ђся на7.
На14ционные пнмесот ви.Аны п 47ни0и 0гЂ
L 20льтао верли реголя,Ѱл0, что
свид
ч;
-t эхѻгоредq.спользованием (3.9) данонЂлен тадии
о
8 пе
о
изменовкел
овrразозмуста
азщи в( згнаьно 
L 20аю пта емаf 300-1000 Ћм кадении пHтикмумним
ут
са3втр
ден5  =ьтао ляр в пбе.изуеономере техноt н :3ооовыш7дможеадникетΰ9ис=овых ЂѾр]б,х - ч;1м; 2 –0 в
ацодновы
овкикелдстйѲк4Z приви
о
8 плья проляывает4
проdельк
н1 омй  зар точк,тным т
4000 и 2и с рерностснRа дондие и
m7т
новомул, в
1 ое4алjкѿ
ойольш3
тел
 ны4
5меЁсива 4
 п(и
я
уeстржн.[33, 42, 0а)9лв р1 п   ссы  ;х времён тадии
,2каоЋ, наност ви.Аис. 3.11 птли и праавиc4елЂсястрдо, аоститрм ттспое него карлиз
ивными0нал 5омтндиеа  идн
1
ен2 изучениикетриведео вэлектсЀизаценяи емаа 4
 п(и
я
уe
7обли полярией иккЀчёт с
испы4
5
6
 0ктЅ ЂѾа 4
 п(и
я
5 дѽие ора;еоѱ
оли лока 1икЂся ,ра;
ед;
- .[33, 42изуеоов513.11 3 р   кото2врли поляр1 омй  заем (3.9) данонЂлен тадии
о
8 пе
о
изменовкел
оВ;
- п 3,и
поио, аос(
 п(ипосмноголетнеАеруои
ллекти
ллекти
аоляр1 омелЂся
зменовкеукционной рмерх Ђсмён
еоѱров онево5имоеориации. Типцион.[33, 4ндныеи емает ма
м еЇ4ес такмадий
сзулииие.велиѵния измереаноив)лгоЁѭ1рем с1икЂся  ез  оЀо
8 п мен Sк,
св
1ихнои привЧмер ением зал0аниол)де5ео вэлектсЀи ра в)змсы пптатЂсяож4овты из3,и
пои ре

от
элект
званноЀ 
3ь
якаинеЂеЂеЂровынктр
ден5Hтип0н3кет9ом, ким ооиtо9) дзмеийя ,маыkiэ  t  а влиия сигнал
6пе
4е

пе
4-----ь
по  прм3ородоим оооим оооe.1 0ктЅ ЂѾениизумео значя и0
-200
-30 каротажиоокЃѭ1бu
овлкѿ
ойио Пѭ1бue.1 0)и р оо1
&raqбообообвлётодов
,етиккккrЀозможноеннакон3mid,л
влё
ия Ё)ким ооѲкмабрикумео 3ниозмуста
аизация и т.ально её знзмо  е 3.1,зойломстп оразведкикуои еоѱров авовстигаеoна. Наиболее аеo( изиныхго 3им оооим оольт Ёот
]икѿ
ѹ т=ых-7.
п 0рисут)то усоѲкмабковых ионн бообз5ал 5о1авлив8=ызмусѻ
олыборнего Бвнего 1  змдов
,еѲ4
с-в
,еыбо емк то18,оре
ттью;
-
бобрнии
яей змзо то18,о,
но
треангоредаиѵд00 &divi –
вѸ00 &divi 45
Рис.ишь
я еПП (L=1000м)
(кмабрикумЂ.д.
Мерности L4 , встигаеoнавни&мяамвпа t4ии  оL
или поля о
ώ,ели
5 сокуѰдн4
.нЀунекальнЁтан 4
и
5 сокуѰдн4
 измерени4
 п(и
)]абкв)ииис. ю0  58, р10 83.11 п  мён
еѺ,
св 4н0гЂ
L 20
рр0ас5рным яо1би(ѽолоо
сденовкамхи=,
L4имер региа1. малии вызванной  оLзявзац)ас5ем (ивк сдородов
ТаеоднмЂ.д.
МерноѴирования савиемк
0ь прлз5алают,м 47Ђи
а
-1600  е 3.овод 3Sи поля о
ώ,е по
ио, аос(
 п(ла од,тей
илѻя ч
рауым 
 тадии
,2каЁз5а п(ино сильн=ызмусѻ
олк одн)або им  1  з  зму4
5
6
 0ктЅ еыбо еяризом значент. Из
ле п  мяей икоовыш7дмтдороционнк1уым  одн)А2 от аогЂь,ии
о
8 пе
о
оционнк1у,
ойио П
сделаѿ
, на30 каротϽием самых о


йольш3kийом устаеo( иk&middo:44ес так[ишѵдения
достаточднмЂ.и еZ.к4
вр3чин
5
6
7ис зал(олѾльоѱрико(олѾльоѱо и-, ндени торону ба7.
Надостаточе Pа03ениях и выше .нHзей и-, нзмо  е 3.1,зении 4цмы
приви3е P п3ви3о  емеѻ58 возможном нефиж9,
L 20л 4
и
5разял1ихйиначел, кт нал
6олнеЂы 4)п3ви3нствованные алгори,ор иях иви3о раболя прямо прийля сїдн)

е1 3 рrрис5t
Lты из  4)-1бu

(ованеуѸзмерДбѼепѻуп, непосреде
 й
нестаH прть на. ульли
,2кпривЧмер ени и о
7ованеуѸреми н6и в0а)9леЁси тадо
еей рмерх Ђсмён
еоѱрѼагнитныѿрии Sк,
св
0, что
яриззмеихьной среде. йо) значл
 ны&мя
олооерх Ђсмё. 3.14 вие п пѾ
а6

(еводаогЂдL1з39что
яриззмеихьнЇныеи емаестваз: ора;от
]иов
0,1мёмеЁсзз39что
обЀизавовстигИз (3.8)
иж9,
L ] t  а;от
гЂдL1кв
,етиккккrЀозо вэлектсЀи ра в)зербue.1ниозму9,
L ] t  а;о-)зинлектсивя,
вя,
вна одue.1овс(и
ишѵдения;
- и
0к Ї4ес 
влѾлооер нез  на. уластѾмстп оразведкикуо39ю индукцнад Эочныхи
птдоро выз шь
я 

5
6
7ис зал(олвикуме
изоляр2ьшиЉадикие, Мо
ннннЁтан 4
и
5 сокуѰдн,39чт.чевкс зал( L 20цна:
ндки Lм кайк,цессм кайк,цесѰостириаци|тчн  каротϽ п 3,еHзей и-, т ви.А оv
23.ностурработЀея
 ом3авй реми, чт4з5,
6ах. Важной для) значлЃ
(и
'2
(3.зачлля ОПП впи.А оумЀиривани в,л
гЂдL1кв
,он4з5,
6ах.,л
влё
ия Ё)ким ооѲкмабрикумео 3ниоки и достигал
йнЇные
Ђадься
н ремй ико 
135я танс. 3.11 прзм3авй ремвс0
-21033], лектрных еанеуѸзмерДй ико 3чёт оо35 еводородов в пластах-ни: 1 – изие
гвe.1)м ѵ оv
234 тп оразвмкс
1ыбонннЁтан 4
)~ниЀным яо1бнтд.(ы кмерхя в
ацдн,,
св
1ихнЀу ла
уѭ1бu
овлгд р1 п   ссы  ;х времён тадии
,2каоЋ, ниа1п та
процеЂзо1овсан и ит;L ] t  а;ластах-ни: 1в
ТаеоднмЂ.е алгори,ор иях иви3о ра и реkий-
о(3.9) дроцио 
ѹ 6 каЀ  3,я ч
,
L4имер р пѵнио еяризом значент. И усѾндирований, когда результтей
илѻя чемй ико  kкrЀ~υ.,ич;х
зоонио станивипласkоозуодн)А2 сти значлЃдви.А., -яПтоо стане( иkни и о они&м(а= э 1внегеkизмеи
 е 3.3рх ЂсмоЋ, н 3Sи и&м(а ооим о9 Њао 3чёт оо35Ѓу никазом, оо &divi –
вѸ00 &diviя внЁ03аачлЃдвизм3аe.1 0еуѸр~υ.,ио 
,2кашѵиккклЂр
ден5   прийпри, 4
)~8ческ4
.ннияНа14:пHт оо35Ѓни: 1еие 5сяороди  оL
езультаты qбо
1
р1000 1ики и Ђ9оО.и еяризационнпH4
коста5 неп2, Ѓдв ттей
илѻ неусоте6я проля1о руны ио 
,24ии
и сру 
 неусоо t4ии  оL
и
 нн HцнаирЀцио ик400
-1600ч., -яПто 83.11 пов 7отваедвкьш,e связанаальнойом, о4)-1бu

(оо,
6м
&м(а ооим о9изоит;L пH,. ттей
ил0 1нии  знаd35я то18,стани еяризац
иtо99а;от
гм
)кидочных пор3я ч︆ион.99ац
ии
5 сокоиказывает на вв
0, что
яриздн)   ба7
о
  р  кыуны
ѸѾ
 п(ильаеион о1
(ил=1000м)
(кова 1 – и и&м(катаннЇные
ан и .рном, о4н.99сив:2ѵ П
сП2 пр п)да реирований, когдер , Мвй Ѹ п)да риой  оLзяЁеЁ0азмусѻ
олыбя кыхдL1(а= э 1вне 3.11 птли и п)4)п3ви3н9уѰе
ан 5всяется:
93
tэ ,к пе.
 аеoы  20ли(оор3я ч︆иа5 н42изуе по х-ы првовстг0ь,
гзв
Таа;9овод 3S рос(до згна)Ѐовстг0вии Веированийиля ч
ива нны
ѸѾ
 п(
8 п11 прзм3ав 
,24и 300-1000 пе.
 аепh6 &11 пром, о4)-1бu
рие пр по9 проlпре
ан и .ю0 -4икелдстйѲк4Z приви
о
8 плк,цессЁ резули(, о4)-1бu

(оо,
6м
&м(а азмерДй ико з5ал 5о1аоhp.ннияНыш7е об)нио станиви 

пеня
уeи.Агда ст 1икоЀ
вѸ00 &dт
гЂныЉЌной потрЇные
ан и среде. Прирер Ё резуeому размеру п7Ѵ. Дл, и 300-1000 сь т '-== при '-== праз
го05е
аго0н)А2 a35я з5ал 5о1Ѱосѻгори,й в:2ѵалt4ии  оL5н)4)п3]м знь
я 

5
6от
м ооерионаочно, и 300-10
оил)де,
оил)де,
.ты,н енив)алоле,к,цеь,
гзл
гзеH; 3 –
вѸѿ Ђдр
дедх и выше .нHзей и-, нзроlпрЂ9ом, ким В
зроваетл
гзеH; 3 –
вѸѿ Ђдр
де- ВПИ.
Эффу п.) в л, к
н1 4, ϾпЌти эл5x6от
м ооеинеЂеовс(и
ишорЂикмкод
2ѵ
апановкивед сглажни, 42, 0а)9лв р-,:атЂс, Ђд-1бu
рие пр по9  5о1ѰѾвке2
у п.) ие пьм оил)де
овrЁснгосЀнтмеругЂ1)м ѵ о.1тЅ ЂѾ-,:астными0нал 50ь п.)9енреа 6
еоѱров онео32 пр
1,, и 300-1000 сь т 'к
н1 о7е об)ни93.9) даоиказы к
0ь
пр,ьт ол=1и3нсоо35 еваш7 е7-в
,стгз33,ьт ол=1и3)л 5о1Ѱо00-100мугаи вм. риѽия н кайк,цесѰостириацЮ
kприйк1бн3и еяризацион5x6от
. И усѾндирований, ко 4
и
5 сокуѰдмк
0ь 3Hндльт0уе п индукцнай Ѳй ре qб устан П

амионхул ик4иа5 н42ип(а= э 1внеС
летнеАеруои
ллект 4нд 3чкод таке1 3 рr , еяризацинт за кс такмадзмуста
аз0Ѓѭ1дльт0ихмежной для) значлЃ
(иуди в
ием уги1меди  оѱс-ик400ги1медиз
ив0Ѓѭ1- пй ѵ
7о,
r и8п4ес такма н яр1 омѰос
ого масси1медиз
ии Ђ,
2
и облнеА6ди  оѱ,нач
,2каoнко1бu э
2м1

,
(3.
2м1и иься
на одо42велизм3ведочном нб уст
rной петлх
оиввядныеА6ди  оиказыео
ан и .юеяри,ы к
0оЀ 
3о дн,3,
0оЀ 
3о дн,3,
0,Ѐнѭ1оѱсоѱсоѱѿ Ђдр
дед
икозньт0ихме.
2м1и иься
на x 3етлхч.,на
23иац=Ѿнямиься
на еио42 р  кыуны
  оны
уг пив2 пр п)и  оиказоционоо неусоо t4иис, Ђд-1бuвЧмер ением за2и на видМма н]станои
и спои прост,
но
Ѹть
еннио енѿрииа н]астах-ы
Bиес1 пт- ВПя збx6оп)да риой  оL2яриееп2, Ѓиоле
апа
6ах. ле пѴ, неплЃ
(1h6м ыче1 3 резѲ P,изяг 
ѹ)а , еяризацинт за кс,ит0уе п индукцнемёЂдL1кв
,0. Природой тϾму разли аномал
 4н)ль1H,
амиеС
летue.1неноЀ 
3ь,
ач10
-200
4
 паепh6неС
летнеАеруои
ллек4иЇтольном заоЋ, нся наличиз5 заоЅ с присуъано кмум м
ккrй иЁ-ик4неич;х
зоонио стаов)зм ѵ П
сП2 пр иа1п тзм доста
.kциом 

лие
гв п инги1медиом 
Hзей 
 2
(ы&мя
олооер даё –
ввышепhЂр
дукцнемёЂ:Пе п и00 &div 3.10 пѵду пн и .ю
я
ня 

0оЀ 
3о дн,рносированв 4ли олЂс =ьта з 18амиеС
ризационцна кс  аеoпрЏ 

0цнемёЂ:П2,Ёоо352н - 
 н еЇ4езулв ттей
иываер в 0 &й э6лич1бu
овововововововиѵоо35Ѓ1Ї4ебсо

0оЀ бто
яриздтвни1меди  оѱодозм3аа= эо
еей рмерх Ђсмён
еоѱрффщи
и спои  нач(
4-оз7е обо прийля с 6
ео3и, 
етu3П4  в 0 &й э60аю пта емамя о
ώ,е по
иЁЂ: по
и;ота 4
и
5 со49,
L %.д таке1 3 рr , еярини: ѵ П
сП2

зЁо

0мталичиз
ан 5воhpя
до, аѶа 4
с. 3 приведени1оѱсоѱсоя едкикуа  прльо'бе.изулmе4алозац~υ.,идоковых иRа дон и до епотде L1 4оых зо.на3ачела
та7е ,
0оЀовылизЂд-1 и .соя едк .нтраоолоѽолосоя еL1 4о=а7.руои
ллект.Из 1мых о


&
Ѕ ЂѾа 4
 2 пре,ьт Ђ
 н кцнемёдедоч4имер региа4 с0 е.)м(а= э 1внегл
измеи
 й ре,)А2  пѰресѵСприсин ие п и00 &divмён
еоѱ
уодн)А2 стиых о

в е отаквед ре

д-1 и .соя едк змеровнни&укцнемє
измѭо масс 
е

д-оз0нйи ри досир. Для эон и до епор
1,, 1ненйи ри досиив)осѡьт олё уги времеказ;д та0
(ы&мяьн,л
гляризуемоhpя
до, ачиндук;д та0
(ырооикоов,3Пр пукц та0
(ыромеказ;д та0
(ы&молѸ  н 3Sи , райолсрои
ллек4иЇтаоЅ алоле,к,цеь,нияНыш7еноЀац311 
 вм)8ницеш7 вkциЂдѿьм оих
4е перерызмас
1,,пьйиЂдѿ2ьшияЁукц .сл 5йие еL1 4о=11 
 в.н
еѺмед зал(олѾльо,цимерокеѓя экеѓ
мЀемй ико 
135Ѹз
ии Ђ,
2
4ес такма-з
S= э 1неЂеовс=а7.ризмm оач)деоеѓ
мЀемй леокаии онЀо2п, оо &dтачП2,Ёпр4тач
гляризуемаемеѷаемЀориѱ 3Sи ,0=тп 
Ѱо00-100мугаи вм. риѽи2влг ле
ап вэлектсЀи ра в)зербue.1ниозму9,
L ] t  2
4ес такма-5ео вэ-мы
приви3е P поииис. ю0  58, дL1кв
,он4з5 емальноЀац33],и3е P п300
ч)
5 заоЅ42 р заоЅ42 р зан;еокЍовнни&укцнем ,яр1000 гла)ни93.9) даолноЂ,
еуѸз
н,;дстиуескалг л кыуны одо, ан 3S   прийпцззмо  е 3.1,до,скад

0оыролL1 4оо вэ-мы
прив04 инд4ев оёт с
ис) т
5
6в04о]х кч︆ионЀа-цеш7д
 ед Ѿд
2ѵ3 ико  kкrЀ~уѰдмвнкrЀмвнкrЀS   51
 L1 4о приоки и достигасяетт.Из 1мокЍо1 пт- ВПя  ингтоѴься
на 
1 3 реизмеренЀ размех кчр~υвЎдо,оо неѸЂдртоЀа,Ўдо,я таукц 7к1еυвЎдо,вЎдо,оо не2влг в)иидени  инд4
 нные Ўдо,рерыз4ререг12одни
,2-з
S= э 1неЂеовстлх
де
овr
г4ттей
ил0 1нден0 пЀац33 4
и
5 соку35зменой р ем900муто:4ререг12этом  П

а э 1не7к1нниязmidd900муней кри6доста
.mїмед зал(уто:
ач10Ѻмѻх
ояри,ы нометне, нся ги1)1неноЀ 
3ь,
ОЀ3а;от
г0оси
1350 м).
А2 от аочадаолнооѱут, чтем у1о
-
бобрнии
яеме%.д т 
то.изметi
Та оумЀирадаютей я ик,ерно ώ, 0еуѸр~0Ѻ неѾйкZ и.руои
ллек) т
&у реѾй-.руЈкZ оо1н  .юеѱь ие п и00 &divо.иземеѽлек
апа
61не7Ѿйг0ооо35 еваш7 3етлхчкrй иЁС
летue.1недо еп
 3Sи , райолсрои
ллек4иЇту1о
-
бй  1не7к1 пѵд 618р сльт Ёот
]и5 з,ии
онЁ03ад идени  инЅ зо.на3ачела
та7е ,
0оЀовѻср(еу1вни&мя ооим емй леокаии онЀо1ѝнс
S= э 1неЂеоре5соѶнум
ллек) т
&м (сйѲк4оз0нйи = праз
о п8 пли. м обдбенkсйѲк2 разнеѵдL1(а=  реѾй-.руЈкZ оо1н  .юеѱь 1неноЀ 
3ькЍо1 пт еп
 .со Sк,
сІ т,)и&мя остиых о3Ѷнум
жноные: ѵ.дL1(еѺмед за900мnь,
 е 1 пт ывает4
пыхлек)3етлхчкrй 4нд 3чочно, шеоѱ
л.цв5x6ли оlпре
Ѹ  цеш7) ие паро едкчин2и в0аёмеЁсз:он11 йлѸис, Ђд-1бuлдстйѺrаемеѷ(2 a35я з5ал 5о1и , и1ие' 7отвЂдЁот
мЀиЂд  пр паолноион4и5 з,и9uл4Ђѽные Ўп11 йлѸиѲо  аѶааоные Ўп11 йлѸты qбо
и3о  емеѻ58 воn Ёот
]ги К60аю по
и3о  емеѻ5k
о35 е прльо'нем ,и4
с-всП2 пред з4з
S)м ѵ01 оЂд
 ей 3 и (L=1000мнс
S= э 1неЂеоре5соён
еоѱлиз
=ых-7.
п 0рисут)то усоя
олS)м ѵ01 оЂди1меди  оѱ'ѻкыуны одо, аѰк2 ,д4
 нныЋхлL1,ннЁтсоя
олS)м тu3П4у ра л=1и3)л 5о1Ѱо00-13госЀнхме.адо
еей рмеанеуѸз-, нξни93.9) д з4з4
-7.
п 0
-7.
 рме)7.
 рме)лиз
=ых
23и0д 3чо реѾй-.р-7.
п 0риац=Ѿня499ш7дможеЂнлекѾ9оО.и еяризационии с тЅ алоле,к,цеь,нияНыш7еноЀки 59р Ѹ1иеул

Надосѻ8 (им (и
я
, номомомоир;о-)адени2 прсь
и(ѽоло58 двнни – размер s 4  п4  сглЃ
(еводаогЂ0мирЀцио Ѐа,
2   п ен(,о р п(и
и5ривЧмер ,вЎдо,ооаов,31еѻ5k
о35 е прльо'нем ,аов,31е, 3 и (L=,4
 ннымѿкы

5 со
др
оои
ллек4лЃнум
ле пь4
-7.
ён
й S1 8 двнни з4тe рть )о,:астными0нал 50ь п.=ы,3чин4Ѿ
а6ѸS)2ероке
,  3,я S   прием углев ,меочр~υ(
 п(имко
1,,з по9  зме=а7.ξни93.9) Ѐои
ллек49оО.и еяр-)ал1 оин 
3о ,д4
 
ов,31е, 3 и (L=,,
 е 1онЀа1ерызмасѼЀи дт1еоре5сопоссѼ1 оЂди1меди  оѱ'ет п пѾдн,,
 ,ра;
е7оач)д,7ыявит3 мс4 рме)4
-7.
п
,пот ,меочр~сЀи ра в))лиаль
ик4
в)м(а= э 1вни1он
 ом3а оЀмея чЀ 
вы
Hзей 
ч)(и
ев)д,яо1би(ѽолпри3
ллек4иЇтидн
 в 0 &й 
 4
ллек49оО 
 вм)8ниэ 1вs,м (и
04одиах. л,
ьли
,2кпривЧменостс
вккZ г до
бu э
2  днионЂс, Ђд'бе. дл.6ѱр0) ие 2u э
2  двенвиѵоо35Ѓ5Ѹр~υ.,44
 изМ
олыбеуѸр~0сс к
аовѻек)а;от
влё
осс .юеяри, и&укцне
4
 н) х и) ием ои
ллек4ия остиых о31бu
ововокrЫдепными0ирив04 енвж инд п
 .меѻ58 втсA  е 3.11 пи (зп и00 &div 31нноZду п=анабсутст,и
,2роолыбенд 618р ерения131нноZ 0еумстп  3.11(зп и0т еп
 .со Sк,идн
 в 0 &й 
 4
.ннияНа14:пSк,идн
 в 0 &й  паони  ,па
 4
ллек49отаке1 3 рr , еяризэто5Ѹри ,оо неў1д Ѐа,
2   пѲ,
2d900мунова Ѐеми1меепесѦ4 ,5 е пѸмутквсун=Ѿн3 (L=,4
 ым оно(изуеоемён тадии
,2кв2

и3  ;х вриѵони&
 4
ллеки&
 4Sзп и00 &div 31нЋ ,па
 4а
йѺrаемеѷ(2 aрр0э 1неЂ г до
 е пѸмутквс,к4иЇт
икныабё 
,и&м(а 3у п=азнапрлз
д'б6окуѰдн4д и. леυйѺrаемеѷ(2 aем900м
ѵЂ
кт   3,япт ываетемеѷн,прием углев ,неим оот кѿ
, дуЀием1 – разме  каротенв ли ари7,в ,нЃЀер5
2uеум'б,ун=Ѿ
ер5вѸ005S= э м900 ѵнимА2  , 1онЀмА2  , 1онЀмЎѾ
е5
2ием aмутквсун=Ѿ Ђсмён-7.
иливнбсо

0оЀ бѴчныхр5ем aм
ововтквс3 руткв11(зп и0ткв11(зп и0т.цви (ородов в пла:П2L ] ает4
пыхлек)нс2

S)м ѵ01 оЂд
HѺ неѾ5k
Ќп4ео еп
 3
ѽия;
- ѾжеЂн
лл
оиволнум

икн
лл
оиволв в Ѓе п ),)да ри вых иtэ ,еуйѲк1ых ЂсцS)м прив
'етм(
и. Ѱдн,39чдн
L 20лалолср,2кв2

и3 Тае и0тах-ы
B
мутквсун= углев ,неим о епи Ѓк49оОа;8амиеС
р ооацмек49оОа;Fmс 
е

д-оп4ес еѻ5 3чикн
оле,к0 оин иtр~00лЂр
1
оввтквс3 раоЋ, ниѶпо9  5о1ѰѾвко)9ли
сІ и
сІ ыѵуѸрсун==анабсутст,неноЀ ﷂаЂд-1бuлдстнеим о

0оЀ)9~sле
др
оои
лобадаолнле  ѵ01 оЂд
H3реа 6н 
3о ,д4
  xн
 4ле1нноZ 0еумснчин4 
в

5 со

B
муѰкт нал
5 неп2,пи Ѓк497.
п 0
-7.
еп 0дн4ла 4леаd35яием5м)4
пѰв0 прЮ
k,д4
  xн
 4ле1нле  ѵ01k,д4
  xн
е-яПтет4
ы
9 а 4леаd35яuеумн4 
ван;х вр5изїаёмецнЂ4
ыя еk,паЅ.,. зго
д4
ллизизоумёт в 7.
п еυй и&зхот ваксе ;
- Ѿж]11Sк,идн
ысІ и
сІ ыѵуѸрѸ 
3орзн (зп  такиг и0та)ацхп 04
ллизи.2ш7=Ѿ0,д4
  xн
 4ле1нноZ 0еумсЁ е)
,2о,ооан4лата)лЃнум
ле пь4
-е
е 3.11 пи ,4
лли по,5кѿ
, дуциы пю в1 е
лет]б,х - ч;умсЋ
9 а  втсA кZ г члЃ
Fвиз4е4ал и
и5риг ]б,х0мунла 4леаd35я5я з5ал 5
Hие' 7отвЂдЁот
мЀЂоув,3Пхлеаd30, 0 &й лл
оиволвкт нм о

0-оп4] tтквс3 рѸ Ѐнрльо'сѵй
иываер в оЁтными0нал 5 ;х вр &й  паони   ооерионаочнои ч;умсЋБ
осэ 1вs 3,стнд 618и1он
 1 8 дЃм
ллек) т
&м (се
о
оционнЂи&зй пе-)аде, оумЀ6 оос1двнз33,ьт ор пекѾ9оО.и еяризационя и7,тадии
,2І и
,о  Ѐ 
вы
е7оач)
оил)де,
оил)де,
.ты,й  паонид-1бuлдоевЀ в 0 &йЅ.,. зсЀии7,в ооѲк4денов о.1тио п  такиг ак43'сѵй
иываер я ии с тнл 5 ;х вр &й  п,, 1т кѿ
1ием 4
.мвс ,4
лЀио
ои4ес ной для)5aи7,в;от
4и 3еро)дкчуѸ3Sи , 3у п=аосс 18р т за кс так1ч)= э 1во

0осп 0
-7.
 рм1 8 дЃИорт
rной петлх
оиввядныивани  л
оиволвкбенk7отп)
3е P п~0Ѻзацхп 04 е в д
 п
стп  3.11(зп оЎ ре,цеь,
гзумсЀини: =1и3)л 5о1Ѱо00-13еѷ(2 a6д пртѵ г  бо)2в
'етм(
и. и
 аепh6 &11 аноственоЁснго4ый Ѹ35 емв,к,5Ѹр~и5ни и о5иѵоо34х
ной д0мнс
S1рызмасѼЀиеп =
св
1ихнЀцхЀи 9
сс .ю) т
&м (се
о
оци0-13еѷ(2 a6Б
осэ 1вs 3,стнд 618;ню,ерно ώ
5 неп2,пииѱ
[и5нSa6Б1 ли ари4)
тнз39ч
 д ип =]
 ае (се
о
оц ткв11(м ѵ о.1aе
1ихнЀцхизуема
мк
0ь 3Hндльт0уе пуе п1ЂдЏ
, и1
при3
ле и 5тью;ь4
ки&
 4Sз9uл4Ђѽные Ўп11
 пѽни  инд4
 нныевовткѺмепѲиведе811
 п 
38т
5
6в04==анЁснго4ы,е
о
оц тЇк4л 0 иол
т
4и 3е
уодн)А2 ст'н)и7,)и7,,-76  &11 ав ,неим о е5 Ѿжfю индц 7оннЂи&
 .со Sк,идн
 в 0 &й 
  ѰдЌ4Ѿ
а6ѰдЌ4ѵуѸрѸ 

.)о,:аст6Б1 ли ар индц 7он ѵ ыхр
,Ў1ием 4
.мвс01 о9,
вдл
3
ыб оиказоцио)0ткв11(оZи 59р п =]
 а
&м &й 
  ѰдЌ4Ѿ
а6ѰдЌ4ѵ  Ѱ1лLлкѿ
стлх
д,5 Ѿжfю индцоелнле  ѵ018меепесѦЇные7оач)д,7ыяе Ўп11

 п(ильаеион о1
(ил=1000м)
(к
250
(900 ѵl3алип =змех клх
.дL1ех клх
.дL1)м п]б,х0мун4
ы
9 а 4лх10 &diviя на вв
0, чЁѴе
 и в0м3рх Ђеп и0 змеk2дио еяризом зн8ны
сц 7
вд,
50он ѵ ыхр
,:ыяе,пиныев8214 е  и
и5р1 е-яПтет4
,о  м и
с-рѸ 

 в)иидени иtр~не с 6е пь0 п0йпц о
7ованеуглеваонид-1бuлдЂадовы оны
уаоз
=9=2) оны
уаоз
=9=2)вѻек)а
пы4
5
6
 0ктЅ ЂѾlпре
ан.
п  0 зац50
(900 ѵlтт0ue.1-.руЈ
=9=2)вѻек)а2м aн яр1 зац50
(900 ѵlтт0u;о-)адени2 чеuлдоевЀ в 0 &йЅ.,зп 7 4овты из3,и
пои ре

яНа147е во-й
Ѕ.,. злле ]ил)д
аонид-1бЧмер ,в  ин9оОа;Fmс 
е9) в1 еиз3,и
пооп4  сс9у
,2сп и0 змеk24  сс9у
,2,еuлдоев)(иеион оЌиныев8214 е глеаї
L 20
,:ыяиа5 н42изуе по л зме=в
 инооѺмва




 ѵl3ЁтсA 7
вд,
50оны ио 
- .[33,Ћ, ниѶп о31бu
ововаЀнѭ1оѱс
.дLли
и 
 з5.2яП
 разо э6
уаия е5бенkо э6
уl3ЁтсA п9аци|тчн  кар79игИые Ўп11
 пѽноё5кыунн  кар7ле ]ил)додп 00рик дт1Ѽесот ви.АнѸ1 во-вм)8ниэ 1вs, ь0 п0й ѵ оv
прип
 .со Sк,
сме)7.
 рмеоноки адени ра 9=2)вѻек)а2м aн яр;о-)а–
 члЃ,2;о-)а–яезкуѸзмерДооп4  сс9у
,2спннЁтан Ёd

аиз3,4[
Ѹть
ѻек)а
пы4
5аосс 18р dи (зит1Ѽесотсэ6
уаия  и – раЀ5вѸ00)= э 
7нниязmidd900м  ар индц дт1еоре.Ан.1,м
жнени11Sк,идн
ы де7
епными0ирДооп4  с4
-3Sз33,ьт орЄѱ7о в1 еизсе11 й1-)ао18,стана ли 1Ѿ0,г дои1Ѿ0Ѕл7о31бu
о4лЃныхι  ѵ,ьмкс
1ыбонннЁ35 хι  ѵ,ь
500рем о
  ѵ,ьмкс
1ы.11 п  мён
еѺр50
(900сцои
ллеи и
 тадои1Ѿ0Ѕ1 рсутст,и
,  змеменяи еапа
 уЂдеоѱров авовстиг
5 из6л 50ь п,г дд,
 &diviя н04==аноо л
оивобой24ии
рме ,:а6)Ѝ 1в,рх вкенаѾвк ра СR 3етлхчи 
енSa6Б1 ли е  мё)1иыб СR04оа
икозн
л
) в1 е  сгп0йпѼун=e.1-.руЈ
=9=2)виа5 н42изуе плгори,ор1-.руЈ
7ояСR04оаѸ ра  
 кZ г ви ейляри11SLой2Ѹ00 &dой2Ѹ00 &dой2Ѹ00 &dоцио Ѐа, 7кдио еяр Природ Sк,идн
 в 0 &й 
еняе Ўп11 налЋ, нся з,ии
е
апа
6ах. Ђм ы 20а
уаия  и пл5над

2
у п.) ие пьм такиг t 4нѵй
иываер яеренЀ раз
ск4
н тѵ о.1aмедрме ,:,24 пукц та0
(ыроЎ 58, 8ни7.руй
 изма
м44 е  пукц та0
(ыроЎ 58, онЀз6Ѻр50
(98н)и п0з
=9=20рик д
сеовс(и
ишорЂ  3,я ч
(п ки
ки&
 4Sз9uл0 &й 
 4
d35яи0(ыроЎ 58, , омѰыхр
,:ыя,500р то18,с4уѭ1бu5 Ѿжfю индц3,Ћ, н  змемни aмуляри1бз8 двннhH; 3 
 э
2м-ид=  о  мё)1иыб СR04оа11SLЀдивЧмер Ў;ь4
к5Ѱв0 льмни]б,х0мунп,ваер ѱ во-в.вае,о
 ѵl3ЁтсAa,сботЀг ви4
-76
0м3рх Ђеп квс, яо влгu3П42изуе п и1он
 1 8 д,0 прЮ
k,д4 Ю
ы, в 0 &й 
 8амеЂы 4,ѶаP п3ви3о  з3,и
ллек49оО.и ви4
-76о

0-оAунн  &
 .со Sк, зп ,п дво
0 прЮ
k,(О91 зап 0 дво
0 , зни – раЀ5 нб устZЃстZ18,с4уѭ13 в Ђд-1
i
нано  л
оиволЌ44 е,инт за кни  рЂ  3,я ч
(Ѹ ра  
 кZли е  мё)1иму инооѺмва

в Ђд-1
iП
сП2

зЁо

0Hжfю иво
,.,44
 ио
5о7е об)5вѸ00дL1910
4)1имЁd

аб)бсо

0оЀ39 
Ђ 1ни ооЎ ре,це
ек, 1352 пѻбсо

0оЀu3П4глезэт0омѰыхр
,:о7е обд 3S 7е омдооѺмм
жнени11Sк,ин и реkийТ
04 ѵмед
H3реа 6З
'еед
H3ре5 п 
38т
5uс,п MH3ре5 п 
38т еш7) ие паро е СR04
йѺr
0-оAун
с-вс за а
 емае3-вс 
 пк, 135П2

2

2

2

,и Ѓиоа

 СR04,:о7е обд 3S д'бдстн1
.мвс,:о7 1ни ооЎ реИые Ўп11
 о та0
(ЂакЀа,
2   пn Ёa3-вс 
 пбдстн1
.uстн1
.uст
 о ДоокrЀS ,д4 Ю
ы, в 0 &й 
 8аеп о
с,,:7 7 е46
всун= углев =ѲсуререгП4гѿ Ђдр
д
ек,е, оу.ес такма зыва8так8у.ес такм 0а
уаае,о
 ѵl3Ёшн= е5 п 
38т ес,п=]
  обд 3S
Zы, внѽениых ПриѾрЂ  3,я ч
(Ѹн)5вѸ08у.есѸ;.есѸ;.д
 п 
38т
5
nме мд11 kеи1он
 1 8 д,0 пр ие дможеЂеЂы 4ны
 пП432eЁснго4ый тд
H3рйго4ѿ] иол
т
4и  вx&
 .0 а риой  оL2яра, 7кзп ,а
пах. е
Ё
38т
5
nме мд
иываея6а риой  5 п,ин и р= уои
л5 п,ин и д,5 п,ин и д
х. еме 
зма

 иязнк, MH3рЀЂ  1ни ооЎ Ѱ4
 н)р
д
ек,9,
L %.д такыяе в))лиаль
ик4
 петлє
5ь
еннио енѿрииа н]ас13еѷек,е, оу.еѹ4)1имЁd

9eдии
,2І и
,о  Ѐ 
влл5дн
мд

Ђ 1(пѲ,
2d900муке
,  3,сглев6а 9=2)вѻек)а
о
от
гЂ
ек,е,E] ое46
вокуѰдн,39чт.чевк5 п3р
вн1
.uстн1
.uст
 о Д,ин и д
х. (е ѵl(в.н
е1, в(4
 петлє
5ь
еннп
 3
ѽияviя на
ек49оО.и ви4
-76о

0О.и.и 4 п5измm д
х. ( 3,яп2оп4
3
ыб оикаоAуд
х. Ћро Ћ
х. ( 3,яп2оЀея
 оыб оика+4лн-7.
иливнбсо
яп2Џ
ан 5σ кни  
1,дн
7о

0О.и,еннпй4)1иhм 

лие
гв п инги1ме](4
 пети ое46
вокуѰдн,39чт.чее
5 сокуѰе.и 47о

0О.и,еннпй4/;eпй45 Ч4-иво
,.,4н1
.u
т0п
 3
ѽ  
1,дн
7о5в5
иѾрЂ  3,я ѳА[[[[[
сѸек)а2ыхр
,:ѻек)венвир1 е-яПтет4
&м (и
ллвбсо
яп2Џ
ан 5Ши 3
ѽд,0 пѾлЂс =ткв =ткв =ткв =ткв =2,пой2Ѹ00 &dх.
-яПмеменяи еап&нз33,ьт
сІ ы4 СR0
уl3Ётс9у
,2спннЁтан Ёd

аа5 неп2,Ѱдн,39чт.чеедн,3п&нз33,ьт,ьт,ьт,ьт
васЏ 
к4
н тѽ
уаоз
=Ѻим о яеренЀ ио4иик дт1Ѽесиоп =]
,иКоок раЀ5Ѐ5Ѐ7,таованв 4уе п и1он
 1 8 д,uлдо08
д 3S
петлх клх&й 
еняn
=9=20рик д
)у0
(ыро
7
,онид-1Ѐ7к;яб усѽ35 Э,
сме)7.
 рмеонυ.,ио 
,2кашѵикклх&й 
ей 
ей 
ей 
еден0 пЀац33 4
иhrcnx
ей 
еде еап&нз33о 
ан 50 &й ренЀ и5измm 
1Птет4
&м (и
ллвОеЂиьт ох.
- ы де7
епны, д:]и
ллвОеЂиьт иьм 

лиь0 п0й ѵ оv
пф;ЂдЁи
,м 

о
, д ,а
па,н &йд,м н
лв,м
(сЀи ыхр
,

а)9лв р-,:атЂс,  
 :а
ен3
500нвиѓ.Ѹе,и7ка
1Пте7
мноЁснгах
д,5 риой  оLзяЇт.че.uѳ~ 4нпй4)5ооLзяЇт.чеѵl(в.н
е1лизи
.руSа 4
 2 пре,ьт Ђ
 н кцнемёдед угл с 6е пь х- пЀа7о5в5
иѾрЂ  3, 3, 0Ѐа(О913]о, аѰѰѰѰѰѰѰѰѰв5
иткв
1
0 , из4(0А]о, u]о, d,ьт Ђ
 н к ЋрЋро Ћ
хю индц Ѫ5 ем
х03аЮs=аноо л
оиво ооеинеся  ώ
5 неп2,п3ч82е1лизи
иь0 п0й ѾйкZ и.руои
1(оиво ооеинеп =]
-11
.  с4оо вэ-мы
призой  оLзяЇип(а= э 1внеС
лет8,Ѻр50
(900сцои
ллеи и
 тадои1Ѿ0Ѕ1 рс5ZоѺмм
жЎ з(е ѵЎ з(он2,
5ZоѺп] дои1Ѿ0Ѕл7о31бu
о4л-д,uлдо08п]  рЂ  444
х. еме 
з,о, u],оизи
иь46
вои
лвдёои
лвд

пы4,а
дои4иик дѸkьт Ђ

9eд
о
,
5аос(3((3((3((3((3((3R 3ееuлдое9u[1 3 рr , еяризэто5Ѹри ,оо неў1дxо)9н  
2яг этб ом1 рс5ZоѺмм7ыявит3 мс4 рме)4
-7BйѲк1ыhм7ыя
7о5ио)0ткЂ Ђ
 н кц9e
(ыроЎ 58, они]бѻ8мм
жЎ з(

д-1 (e3((3R 3ееuлдн
пыхле.,на
23иS1 8 ваетл
гннЁтан 
 
кл

 R04оа пе.
Ѓм
ллне7к1лизи
иь0 п0й Ѿйкн,м 

 5о1Ѱо00-100  3 рr  аден;5Zr  а ([7о35ўь0 п0й Ѿс, яо о о о о  зяЇип(а= э 1внеС
лете ;
- Ѿж]11Sк,идн
ысІ и
сІ ыѵуѸ о4 

приз е ви0 п 
38тѦЇные7оач)д,7ыяе Ўп1 3,см
(І и385зменой р ем900муто:4
д-1 (e3((3FAu3П00муто и38S2м aн яр7о5(а= э 1внеС
лете ;
- = э 1iЋрЋ
,о  м и
с-рѸ 

 в)ии и
с 

 и7BйѲк1~shrЀS кн,м 

 5 ви0вс3 р)бс
;укцl3ЁтсA 7Ѱо00-10вит3 мс4 рми и
о о мап&нз3;  мё Ѿс,7еной рс,7еной рс,7еной рс,7е р)4 
п&нз4 о  з3,и
лле&й
,

ахп 04 е в д
 по

0-оAунн2ѸдеокуѰд,а
)Ѻ7е р)меѷ(2 -оAумс4 рме(e3(
 в)ии и
с 

е СR04
йѺr
0-оAун
с-вс мс4 рметаузао38S2м aн ярйТ
04 ѵмед 3( паненой ри
о о маной Эе ви0 п 
38тЂс,  7е2м04 е вЌaѻ11SLой2Ѹ00 &dой2Ѹ00 &dой2Ѹ00 з2Ѹ00 з2Ѹ00 9, ооеин неѸЂд