Развитие методики обработки магнитотеллурических данных и ее применение при исследовании электропроводности литосферы восточной части Балтийского щита тема диссертации и автореферата по ВАК РФ 04.00.22, кандидат физико-математических наук Смирнов, Максим Юрьевич

• Смирнов, Максим Юрьевич
• кандидат физико-математических наук
• 1998

Оглавление диссертации кандидат физико-математических наук Смирнов, Максим Юрьевич

С одержание

Введение

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

1.1 Спектральное преобразование магнитотеллурических реализации

ГЧ 7Т Л Л

предварительная подготовка исходных данных

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

1.4 Алгоритм обработки магнитотеллурических данных 28 Выводы к главе 1

Глава 2. Тестирование предложенного алгоритма на ис

кусственных и реальных данных

2.1 Тестирование с использованием искусственных рядов

2.1.1 Алгоритм AP-предсказания

2.1.2 Отбраковка по когерентности

2.1.3 Коррелированная помеха

2.2 Тестирование с использованием реальных данных

2.3 Использование предложенного алгоритма для обработки вариаций ЕЭМП в интервале периодов 60240 мин

Выводы к главе 2

Глава 3. Аппаратура и методика для AMT и МТ зондирований

3.1 Магнитотеллурическая станция МТС-97

3.1.1 Аналоговая часть МТ станции

3.1.2 Цифровая часть МТ станции

3.1.3 Определение уровня шумов МТ станции

3.1.4 Результаты полевых испытаний МТС-97

частот 1-1000 Гц

3.2.1 Характеристики ЕЭМП Земли в диапазоне частот 1-1000 Гц

3.2.2 Цифровая AMT станция

Выводы к главе 3

Глава 4. Магнитотеллурические исследования электропроводности коры и мантии восточной части Балтийского щита вдоль профиля SVEKA-2

4.1 Геологические представления о строении и эволюции восточной части Балтийского щита

4.2 Магнитотеллурические исследования электропроводности коры и мантии восточной части Балтийского щита вдоль профиля SVEKA-2

4.3 Магнитотеллурические зондирования вдоль профиля параллельного профилю SVEKA-2

Выводы к главе 4

Заключение

Приложение А Цифровая фильтрация

Приложение Б Результаты AMT зондирований на участке профиля Суоярви-Выборг

Литература

Введение диссертации (часть автореферата) на тему «Развитие методики обработки магнитотеллурических данных и ее применение при исследовании электропроводности литосферы восточной части Балтийского щита»

Введение

Актуальность темы

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

Большой интерес с точки зрения изучения геоэлектрического строения литосферы представляет собой Балтийский щит [32, 33, 72, 73]. Он является наиболее крупным и стабильным выступом древнего кристаллического фундамента Восточно -Европейской платформы. На Балтийском щите практически отсутствует осадочный чехол, что позволяет весьма эффективно использовать метод магнитотеллуричеекого зондирования для исследования особенностей строения литосферы. В стабильных регионах древние тектонические процессы оставили много следов в виде хорошо проводящих областей в коре и верхней мантии. Эти проводники обусловлены в основном графито-и/'или флюидо-насыщенными породами [16], образовавшимися в результате тектонических преобразований, и указывают на положение шовных зон и границ геоблоков [22, 35]. Отдельные результаты магнитотеллурических исследований указывают на наличие электрической анизотропии литосферы, которая, по всей видимости, является следствием тектонических процессов. Исследование этого явления также дает возможность судить о процессах эволюции литосферы. Очень важным с точки зрения строения Земли является вопрос о положении астеносферного слоя и его физических параметрах. Именно глубинные МТ зондирования позволяют расширить наши знания в этом направлении. Кроме того, метод магнитотеллуричеекого зондирования дает информацию о распределении электропроводности в коре,

что имеет большое практическое значение при решении задач разведочной геофизики.

Таким образом, исследование методом магнитотеллуриче-ского зондирования геоэлектрического строения Балтийского щита является актуальной задачей.

На протяжении нескольких лет на кафедре Физики Земли СПбГУ проводятся АМТ-МТ зондирования на Балтийском щите в расширенном интервале периодов. Пункты зондирований сосредоточены вдоль профиля Териберка-Ковдор-Суоярви-Выборг [32]. Этот профиль близок к международному геотраверсу I, вдоль которого имеются гравиметрические, магнитные и сейсмические данные. В каждом пункте регистрировались 5 компонент естественного электромагнитного поля (Ех, Еу, Нх, Ну, Hz) тремя аналоговыми станциями: AMT, МТЛ-71, МЭВС-5 соответственно в интервале периодов 1СГ3 — 0.2 с; 0.-5 — 10° с; 102 - 104 с. Обработка проводилась с учетом тензорного характера импеданса. Накоплен большой опыт обработки и анализа магнитотеллурических данных. Предложена и развита методика интерпретации результатов. Получены существенные результаты о строении коры и мантии восточной

Г», о Т">

части оалтииского щита. В настоящее время сотрудниками кафедры выполняются работы вдоль международного геотраверса SVEKALAPKO, задачей которого является комплексное изучение строения литосферы Балтийского щита. Восточная часть профиля, на которой сейчас выполняются зондирования, пересекает зону сочленения двух геоблоков Карельского и Беломорского. Геологическое строение и эволюция этой зоны и по сей день вызывает много дискуссий [21, 38, 51, 24].

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

являет высокие требования к регистрирующей аппаратуре.

Необходимым этапом успешного использования метода маг-нитотеллурического зондирования является обработка полевого материала, позволяющая устойчиво определять компонент тензора импеданса на поверхности Земли по временным записям естественного электромагнитного поля (ЕЭМП). Значительный успех в повышении качества обработки достигнут на пути использования робастной статистики. Первые работы по применению робастных статистик для обработки магнитотеллуриче-ских данных появились в конце 80-ых годов [61, 62, 64]. В этих работах использовалась робастная М-оценка Хьюбера [57] задачи регрессии, которая заменила метод наименьших квадратов. Ро-бастный подход позволил снизить влияние помех на результат вычислений компонент тензора импеданса [77, 65, 40, 68]. Было показано несомненное преимущество робастных статистических процедур перед методами, использующими метод наименьших квадратов в случае, когда наблюденные данные осложнены различного рода помехами [67].

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

«у х «у

торов, сказывающихся на результатах обработки. Обработка ведется в очень широком интервале периодов (10~3 — 104 сек), что создает дополнительные трудности при выборе оптимальной процедуры.

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

До последнего времени на кафедре Физики Земли для проведения МТЗ использовалась аналоговая аппаратура. Это позво-

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

ПП Г*

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

Цель работы:

заключается в разраоотке цифрового комплекса регистрации и обработки магнитотеллурических данных в расширенном интервале периодов и исследовании магнитотеллурическим методом геоэлектрического строения литосферы Балтийского щита на восточном участке геотраверса ЗУЕКАБАРКО. Задачи исследования:

-разработать алгоритм обработки магнитотеллурических данных, позволяющий устойчиво определять компоненты тензора импеданса в интервале периодов (10_3 — КЗ4 сек);

-обеспечить цифровую регистрацию ЕЭМП в расширенном интервале периодов (10~3 —10^ сек), разработав цифровую часть и программное обеспечение для АМТ-МТ станций;

-выполнить АМТ-МТ зондирования на восточном участке геотраверса ЗУЕКАБАРКО с целью изучения геоэлектрическо-

го строения зоны сочленения Беломорского и Карельского геоблоков.

Научная новизна:

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

2. Разработана методика цифровой регистрации магнитотел-лурических данных в расширенном интервале периодов. Созданы алгоритм и программы для цифровых МТ и AMT станций.

3. Получены новые данные о геоэлектрическом строении зоны сочленения Карельского и Беломорского геоблоков на геотраверсе SVEKALAPKO.

Практическая ценность:

- разработанные алгоритмы и программы для обработки маг-нитотеллурической информации в настоящее время являются основным рабочим инструментом для сотрудников Института Физики, Санкт-Петербургского Университета при проведении МТЗ в расширенном интервале периодов;

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

- создан пакет программ (Delphi 3.0), позволяющий практически по единой методике обрабатывать данные магнитотеллури-ческих зондирований в широком интервале периодов;

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

t-i J

а также в целях разведочной геофизики.

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

1. Предложенный алгоритм обработки магнитотеллуриче-ских данных, основанный на использовании робастных оценок, позволяет повысить качество обработки магнитотеллурических данных в широком интервале периодов (10~3 — 1Q4 сек).

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

3. Определено геоэлектрическое строение коры и верхней мантии от первых сотен метров до 300 км вдоль восточной части геотраверса SVEKALAPKO. Выявлено поднятие подкорово-го проводящего слоя в зоне сочленения двух reo блоков Карельского и Беломорского с глубины 150-200 км до глубины 30-60 км.

Апробация работы и публикации.

Результаты исследований, представленных в работе, докладывались на XXI генеральной ассамблеи IUGG (Колорадо, 1995), на международных конференциях "Problems of Geocosmos" (Санкт-Петербург 1996, 1998), на международной геофизической конференции и выставке "Москва-97" (Москва, 1997), на международном Workshop по проекту SVEKALAPKO (Финляндия, 1997), на международном 14-ом Workshop по электромагнитной индукции в Земле (Румыния, 1998). По теме диссертации опубликовано 14 работ.

Объем и структура диссертации

Диссертация состоит из введения, четырех глав, двух приложений и списка литературы. Работа изложена на 140 страни-

цах, содержит 39 иллюстраций и 2 приложения. Библиография включает в себя 92 наименования.

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

Заключение диссертации по теме «Геофизика», Смирнов, Максим Юрьевич

Выводы к главе 4:

1. Летом 1997 г. были проведены магнитотеллурические исследования в восточной, наименее изученной, части профиля SVEKA-2 международного геотраверса SVEKALAPKO. Сделаны АМ'Т зондирования с шагом 5 км на участке от Белого моря (пос. Лоухи) до поселка Пяозерский. Всего выполнено 19 зондирований. Через каждые 20 км проводились зондирования в расширенном интервале периодов от 0.001 до 10000 с.

2. Характерной особенностью кривых зондирования во всех пунктах профиля является большое расхождение максимальных и минимальных кривых зондирования и устойчивое положение уровня максимальных кривых: максимальные кривые в области больших периодов (Т > 10 с) практически совпадают и идут выше глобальной кривой в области суточных вариаций, что позволяет считать их "поперечными". Наилучшим образом согласуются с глобальной кривой минимальные кривые, которые приняты за "продольные".

3. Практически на всем профиле сопротивление коры уменьшается с глубины 20-30км, где оно в большинстве пунктов близко к 1000 Ом.м и достигает минимального значения 100 Ом.м на глубине 150-200км. Положение минимума можно грубо связать с нижней границей "астеносферы", поскольку дальше сопротивление возрастает.

4. Для подтверждения ранее сделанных выводов в полевой сезон 1998 г. были проведены дополнительные магнитотеллурические зондирования вдоль профиля параллельного профилю Лоухи-Пяозерский. Он расположен южнее последнего примерно на 40 км, начинается в зоне КГПЗ (юго-восточная часть оз.Топозеро) и заканчивается в центральной части Беломорской провинции (примерно 20 км от берега Белого моря пункт МТЗ и 5 км - АМТЗ). Было выполнено 5 МТ зондирований с расстоянием между пунктами примерно 20 км и 19 AMT зондирований с расстоянием примерно 5 км.

5. Характерная особенность кривых зондирования во всех пунктах профиля SVEKA-2 - расхождение максимальных и минимальных кривых зондирования в МТ диапазоне и устойчивое положение уровня максимальных кривых, сохранилась и на профиле SVEKA-2-98.

6. Максимальные кривые в области больших периодов (Т > 10 с) очень близки друг к другу и идут выше глобальной кривой в области суточных вариаций, что позволяет считать их "поперечными". Наилучшим образом согласуются с глобальной кривой минимальные кривые, которые мы и в этом случае будем считать за "продольные". В настоящее время результаты обработки вдоль этого профиля пока окончательно не проинтерпретированы. Длиннопериодные минимальные кривые в пунктах Топозеро, Хаттомозеро, Энгозеро достаточно похожи и имеют глубокий минимум на периодах 2000 - 4000 сек. Эти пункты расположены в зоне зеленокаменного пояса и крупнейшего надвига Беломорского подвижного пояса. Минимум отражает поднятие нижней границы "астеносферы" в этой области примерно до 30 - 60 км.

Заключение

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

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

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

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

4. Создана система для единой обработки магнитотеллурических данных в широком интервале периодов (Ю-3 — 1Ü4 сек), в которую включены все алгоритмы, разработанные в диссертации. Программа реализована для операционных сред Windows 95/98/NT. Разработаны методические указания к проведению магнитотеллурических зондирований с использованием новой цифровой аппаратуры.

5. Проведены МТ и AMT зондирования на восточной части геотраверса SVEKALAPKO. Построен геоэлектрический разрез на профиле. Выявлено наличие аномалии электропроводности соответствующей, по геологическим данным, зоне стыковки двух геоблоков Карельского и Беломорского.

6. Для подтверждения ранее полученных данных выполнены МТ и AMT зондирования на профиле, параллельном предыду-щем}г профилю. Построен псевдоразрез по профилю. В результате предварительного анализа данных, на этом профиле также, выявлено поднятие границы астеносферы с глубины 100-150 км до глубины -30-60 км в зоне стыковки двух геоблоков.

Список литературы диссертационного исследования кандидат физико-математических наук Смирнов, Максим Юрьевич, 1998 год

Литература

[1] Алексеев А.Н., Ланцов А.Е., Никитский В.Б. Наставление по эксплуатации магнитотеллурической лаборатории МТЛ-71. М., изд. ВНИИ Геофизика, 1972, 56 с.

[2J Антопъю Цифровые фильтры. М., Радио и связь, 1983.

[3] Астапенко В.И. Система цифрового анализа естественного электромагнитного поля. дис. кан. физ-мат. наук Л.,1979.

[4] Ват М. Спектральный анализ в геофизике. М., Недра, 1980.

[5] Безрук И.А., Сафонов A.C., Казанцева Е.В. и др Пакет программ обработки и интерпретации цифровых материалов электроразведки // Прикладная геофизика. 1986, Вып.116.

[6] Беляев М.Д., Билибина Т.В., Богданов Ю.Б. и др. под ред. Сидоренко A.B. и Билибиной Т.В. Металлогения восточной части Балтийского щита. Л.: Недра, 1980, 247 с.

[7] Бендат Дж:,., Пирсол А. Измерение и анализ случайных процессов. Москва: Мир, 1991.

[8] Бердичевский М.Н. Электрическая разведка методом магнит о т е л л у риче ског о профилирования. -М., Недра, 1968

[9J Бердичевский М.Н., Безрук И.А., Чинарева ü.M. Ма.гнито-теллурическое зондирование с использованием математических фильтров // Физика Земли. 1973.№3.

[10] Бердичевский М.Н., Ваньян Л.Л., Васин Н.Д. и др. О нормальном геоэлектрйческом разрезе // Физика Земли. 1980.№2,

[11] Бердичевский М.Н.. Булычев А.А., Гайнанов А.Г., и др. Геофизические модели флюидонасыщенных зон литосферы и астеносферы // Вестник МГУ. 1994. Сер.4. Геология.№2. С. 17-24.

[12] Блиох П.В., Николаепко А.П., Филиппов Ю.Ф. О резонансных явлениях в полости Земля-ионосфера. -Киев: Наукова думка, 1977, с. 200.

[13] Бобова В.П., Владимирский Б.М., Пудовкин М.И. Спектр колебаний геомагнитного АЕ-индекса в диапазоне периодов 130-200 минут как возможное отражение спектра глобальных осцилляций Солнца. // Магнитосферные исследования. 1993. №.18, С.40-54.

[14] Вагин С.А., Смирнов М.Ю. Магнитотеллурическая станция МТ'С-97. // Геофизический вестник. 1998.№.3.

[15] Ванъян Л.Л., Дебабов Л.С., Юдин М.Н. Интерпретация данных магнитотеллурических зондирований неоднородных сред. М.:Недра, 1984. 196 с.

[16] Ванъян Л.Л. Электропроводность Земной коры в связи с ее флюидным режимом // Коровые аномалии электропроводности. Л., 1984.

[17] Варенцов И.М., Голубев П.Г., Гордиенко В.В., Соколова Е.Ю. Исследование глубинной геоэлектрической стрзжту-ры вдоль профиля Линкольн Лайн (эксперимент ЭМСЛАБ) // Физика Земли. 1996.№4. С. 124-144.

[18] Васин Я.Д. Новгщкий Г.П. Магнитотеллурические исследования в центральной Карелии // Вопросы геофизики 1986, Л.: ЛГУ, Вып.31.

[19] Васин Н.Д. Геоэлектрическая характеристика разреза юго-западной Карелии. // Методы интерпретации геофизических данных для сложных геоэлектрических разрезов. Зап. ЛГИ. 1988. Т. 13, с.57-63

[20] Владимиров Н.П. Метод магнитотеллурического зондирования. М:Наука; 1979

[21] Глебовицкии В.А., Миллер Ю.В., Другова Г.М. и др. Структура и метаморфизм Беломорско-Лапландской коллизионной зоны. // Геотектоника 1995 6.

[22] Голод М.И., Клабуков В.Н., Гришин A.C. Глубинная электропроводность Карельской части Балтийского щита. В кн. Глубинная электропроводность Балтийского щита. // Кар.ф-л АН СССР, 1986, с 7-18.

[23] Голъденберг Цифровая обработка сигналов. М., Наука, 1985.

[24] Григорьева Л.В., Шинкарев Н.Ф. О геодинамической природе зоны сочленения Беломорит и Карелит. // Вестник СПб-ГУ сер 7, 1997, вып 3.

[25] Дженкинс Г., Ватте Д. Спектральный анализ и его приложения. М.: Мир, т. 1, 1971, т.2, 1972.

[26] Жамалетдинов A.A., Семенов A.C. Электропроводящие породы кристаллического фундамента - объект глубинных геоэлектрических зондирований // Коровые аномалии электропроводности. Л., 1984.

[27j Канасевич Э. Анализ временных последовательностей в геофизике. М.: Недра, 1985.

[28] Кашпровски.й В.Е. Грозы как источник волн сверхнизких частот в геоэлектроразведке. Киев. АН СССР. 1965. с.138.

[29] Клейменова Н.Г. Современные представления о природе высокочастотных вариаций электромагнитного поля Земли (1 Гц-1 кГц). // Изв. АН СССР, сер. геофиз., 1963, N12, с.1798-1813.

[30] Ковтун A.A. Использование естественного электромагнитного поля при изучении электропроводности Земли. Л.: ЛГУ. 1980.

[31] Ковтун A.A. Строение коры и верхней мантии на северо-западе Восточно-Европейской платформы. Л.: ЛГУ. 1989. 284 с.

[32] Ковтун A.A., Вагин С.А., Варданянц И.Л., Коквина Е.Л., Успенский Н.И. Магнитотеллурические исследования строения коры и мантии восточной части Балтийского щита, // Известия РАН. Физика Земли. 1994,№ 3, с.32-36.

[33] Ковтун A.A., Легенькое а И. П. Ложные аномалии при одномерной интерпретации кривых МТЗ. // Известия РАН, Физика Земли, 1995,№9, с.77-83.

[34] Ковтун A.A., Смирнов М.Ю. Особенности вариаций естественного электромагнитного поля в интервале периодов 60240 мин. // Геомагнетизм и аэрономия. 1996, №.3.

[35] Ковтун A.A., Вагин С.А., Варданянц И.Л., Легенъкова Н.П., Моисеев O.E., Смирнов М.Ю., Успенский Н.И. Строение коры и мантии по профилю Суоярви-Выборг по магнитотел-лурическим данным. // уч.зам. СпбГ'У, серия Физическая, 1998, №3.

[36] Ковтун A.A., Корья Т., Вагин С.А., Варданянц И.Л., Легенъкова Н.П., Смирнов М.Ю., Успенский Н.И. Магнито-

теллурические исследования электропроводности коры и мантии восточной части Балтийского щита на геотраверсе SVEKALAPKO. // Вопросы геофизики. Вып.35. СПб.: Изд-во СПбГУ, 1998.

[37] Крылов С.М. Использование электромагнитных резонансов полости 3 ем л я-ионо с ф ера в м агнит о т е л лу рич е ских зондированиях. В кн. Естественное ЭМП и исследование внутреннего строения Земли М.: Наука. 1971.

[38] Лобач-Жученко С.Б., Чекулаев В.П., Степанов B.C. и др. Беломорский пояс- позднеархейская акреационно-коллизионная зона Балтийского щита. Докл. РАН 1997. т. 352.

[39] С.Л. Марпл-мл. Цифровой спектральный анализ и его приложения. М.: Мир, 1990.

[40] Нарский Н.В. Робастный метод обработки магнитотеллури-ческих данных. // Изв. РАН, Физика Земли, 1993,№2, стр.24-2е.

[41j Петрова Л.Н. Сейсмогравитационные колебания Земли и возможный механизм их образования. // Биофизика. 1992. Т.37. Вып.З. С.508-516.

[42] Попов М.К. Геоэлектрическая модель Земной коры в районе Ладожской аномалии электропроводности по данным АМ-ТЗ. дис. кан. физ-мат. наук Л.,1987

[43] Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. Москва: Мир, 1978.

[44] Распопов О.М., Клейменова Н.Г. Возмущения электромагнитного поля Земли. Ч.Ш. ОНЧ-излучение. Л.: Изд-во ЛГУ. 1977.

[45] Рокитпянский И.И., Васин Н.Д., Голод М.И. и др. Аномалия электропроводности на юге Карелии. // Геофиз. сб. АН СССР. Киев, 1979. Вып. 89, с. 35-36.

[46] Рокитпянский И.И., Кулик С.Н., Рокитпянская Д.А. Ладожская аномалия электропроводности // Геофиз. журн. 1981,

Т.8.№ 2. с.97-99

/ /

[47] Светов B.C., Шимелевич М.И. Определение линейных связей между компонентами магнитотеллурического поля // Физика Земли. 1982.№5.

[48] Семенов В.Ю. Обработка данных магнитотеллурического зондирования. М.: Недра, 1985.

[49] Смирнов М.Ю. Обработка магнитотеллурических данных с использованием робастных статистических процедур. // Вопросы геофизики. Вып.35. СПб.: Изд-во СПбГУ, 1998.

[50] Смирнов М.Ю., Успенский Н.И. Обработка МТ-информации цифровой AMT станцией. // Геофизическая аппаратура. 1998,№ 11.

[51] Сло£унов А.И. Верхнеархейская Керетская гранит - зеле-нокаменная система Карелии. // Геотектоника. 1993.№5.

[52] Успенский Н.И. Геоэлектрический разрез земной коры в отдельных районах Балтийского щита по данным М'ТЗ. дис. кан. физ-мат. наук Л., 1985

[53] Файнберг Э.Б. Глобальное геомагнитное зондирование // Математическое моделирование электромагнитных полей. М., 1983.

[54] Федоров А.Б., Елисеев A.A., Кочеров А.Б., Ратников К.Д., Редъко Г.В. Цифровая аппаратура для аудиомагнитотел-лурических методов AMTF.// Геофизическая аппаратура, 1998, № 101., с.81-93.

[55] Хампелъ Ф., Рончетти Э., Рауссеу П., Штаэлъ В., Робаст-ность в статистике. Подход на основе функций влияния. М.: Мир, 1989.

[56] Хемминг Р. Цифровые фильтры. М.: Недра, 1987.

[57] Хьюбер П. Робастность в статистике. М.: Мир, 1984.

[58] Adam A., Kaikonen P., Hjelt S.E., et.all. Scalar audio-magnetotelluric measurements in Hungary // Geoph. Transactions. 1984. У.30Ж1. P.47-62.

[59] Baiser M., Wagner C.A. Observation of Earth-ionosfer cavity resonances. - Nature, 1960, v. 188, p.638-641.

[60] Campbell W.H. An analysis of the spectra of geomagnetic variations having periods from 5 min to 4h. // J.Geophys. Res. 1976. У.81Ж7. P. 1369-1390.

[61] Chave A.D., Thomson D.J., Some comments on magnetotelluric response function estimation // J. Geophys. Res. 1989. V.94.Jn^B10. P. 14215-14225.

[62] Chave A.D., Thomson D.J., Ander M.E., On the robust estimation of power spectra, coherences and transfer functions //J. Geophys. Res. 1987. V.92. P.633-64S.

[63] tifron B., Bootstrap methods: Another look at the jackknife. // Ann. Statist. 1979. V.7.

[64] Egbert G.D., Booker J.R., Robust estimation of geomagnetic transfer functions. // Geophys. J. Royal Astron. Soc. 1986. V.87. R173-194.

[65] Egbert G.D., Livelybrooks D. W.. Single station magnetoteliuric impedance estimation: coherence weighting and regression M-estimation. // Geophys. 1996. V.61. P.964-970.

[66] Filloux J.H. North Pacific Magnetoteliuric Experiment. // J.Geomag. Geoelectr. 1980. V.32, Suppl.l. P.33-43.

[67] Jones A. et. all. A Comparison of Techniques for Magnetoteliuric Response Function Estimation. J. Geophys. Res. 1989. V.94.№310.

P-t i r\. A -i -1 i -i il

_ .i42ui-i42i5.

[68] Garcia X., Chave A.D., Jones A.G. Robust processing of Magnetoteliuric data from the auroral zone. J. Geophys. Res. 1998. (in press)

[69] Kaikkonen P., Vanyan L.L., Hjeli S.-E., Shilovskiy A.P., Pajunpaa K., Shilovskiy P.P. A preliminary geoelectrical model of the Karelian megabiock of the Baltic shield // Phys. Earth Planet. Inter. Vol.32. P.3G1-305.

[70] Kopytenko Yu.A.,Kopytenko E. A., Amosov L. G., Zaitsev B.B., Voronov P.M./Timoshenkov Yu. P. Magnetovariation complex MVC-2. Proceedings of the VI-th Workshop on Geomagnetic Observatories Instruments, Data Acquisition and Processing (Sept. 18-24, 1994, Dourbes, Belgium), Editor Jean L. Rasson, Publications scientiiique at technique 003, Institute Royal Meteorologique de Belgique, 1996.

[71] Kopyienko Yu.A., Zaitsev D.B., Kopyienko E.A., Voronov P.M., Araosov L.G, Timoshenkov Yu.P. Magneto-Variation Complex MVC-1DG. // Proceedings of the I-st International Conference on Marine Electromagnetic "Marelec-97",23-25 June 1997,London,UK.-1997.-P.P10.

[72] Korja T., Koivukoski K. CrustuI conductors along the SVEKA profile in the Fennoscandian (Baltic) shield, Finland // Geophys. J. Int. 1994,№ 116, P. 173-197.

[73] Korja T. Electrical conductivity of the lithosphere - implications for the evolution of the Fennoscandian shield // Geophysica, 1997, Vol. 33,№ 1, P. 17-50.

[74] Kovtun A.A, Smirnov M.Yu. The investigation of long period variation of natural EM field for using them in MT soundings. // Abst. IUGG XXI General Assembly. Boulder. Colorado. July 2-14 1995. P.B170.

[75] Kovtun A.A., Vagin S.A., Legenkova N.P.} Moiseev O.N., Smirnov M.Yu., Uspensky N.I. Structure of the crust and upper mantle along the profile Suoyarvy- Vyborg by AMT-MT data. / Book of Abstracts, international Conference on Problem of Geocosmos.June 17-23, 1996.

[76] Kovtun A.A., Vagin S.A., Smirnov M.Yu., Uspenskiy N.L Preliminary results of magnetotelluric investigations on the profile SVEKA-2 // Abstracts, 2nd SVEKOLAPKO WS, Lammi, Finland 27-30 nov. 1997.

[77] Larsen J.C., Mackie R.L., Manzella A., Fiordelisi A., Rieven S. Robust smooth magnetotelluric transfer functions. Geophys. J. Int. 1996. V.124, P.S01-819.

[78J Osipova I.L., S.-E. Ejelt, Vanyan L.L. Source field problems in northern parts of the Baltic shield. //' Physics of the Earth and Planetary Interiors. 1989. V.53. P.337-342.

[79] Porokhova L.N., Kharlamov M.lvl. The solution of the one-dimension inverse problem for induction soundings by an efficient linearization technique. // Physics of the Earth and Planetary Interiors. 1990. V.60. P.68-79.

[80] Siegel A.F. Robust regression using repeated medians. // Biometrika, 1982, V.69, P.242-244.

[81] Sims W.E., Bostic F.X.. Smith H. W. The estimation of magnetotelluric impedance tensor elements from measured data. // Geophysics, ¥.36, 1971, P.938-942,

[82] Smirnov IvL Yu. Program for the real-time AMT station. // Thesis of International Conference of Geokosmos. June 16 1996.

[83] Smirnov M. Yu. Technique of the robust real-time processing of AMT-MT data. // Thesis of SvekaLapko workshop. December 1

[84] Smirnov M. Yu. MT data processing technique and the results of AMT soundings along the profile part Suoyarvi-Vyborg. // Abstracts, International Geoscience Conference, 15-18 September, Moscow, 1997.

[85] Smirnov M. Yu. The Results of the Robust Processing of Magnetotelluric Data Along the Profile SVEKA-2 // Thesis of International Conference of Geokosmos, St.Petersburg, 1998.

[86] Smirnov M. Yu. Magnetotelluric Transfer Functions Robust Estimation with High Breakdown Point. // Book of Abstracts,

14th Workshop on EM induction in the Earth, Sinaia, Romania, August 16-22, 1998.

[87] Tzains D. Beamish D. Audiomagnetotelluric sounding using the Schuman resonances. // Journal of Geophysics. Zeitschrift fur Geophysics, 1987, V.61,№.2.

[88] Tzains B. Beamish D. A High-Resolution Spectral study of audiomagnetotelluric data, and noise interactions. // Geophysical journal - Oxford, 1989, V.97,№.3.

[89] Vanyan L.L., Demidova T.A., Palshin N.A., Zhamaletdinov A.A.. Kuksa Yu.I., Kaikkonen P., Korja T. Interpretation of the DC soundings in the northeastern Baltic shield // Phys. Earth Planet. Inter. 1989. Vol.54. P. 149-155.

[90] Wannamaker E., et. all. Magnetotelluric Observation Across the Juan de Fuca Subduction System in the EMSLAB Project. // J. Geoph. res., 1989, V.94,№.B1Q, P.14.I11-14.125.

[91] Wigih D.E., Bostick F.X. Cascade decimation - a technique for real time estimation of power spectra. // Contributed paper at Proc. IEEE Coni on Signal Proc., p.626-629, Denver, 1980.

[92] Geoelectrical models of the Baltic shield. // Final rep.№.13 between the Academy of Finland and the Academy of Sciences of the USSR.