Научная статья на тему 'ЧУВСТВИТЕЛЬНОСТЬ МАГНИТОТЕЛЛУРИЧЕСКИХ ЗОНДИРОВАНИЙ К ТИПИЧНЫМ АНОМАЛИЯМ ЭЛЕКТРОПРОВОДНОСТИ В ТЕКТОНОСФЕРЕ'

ЧУВСТВИТЕЛЬНОСТЬ МАГНИТОТЕЛЛУРИЧЕСКИХ ЗОНДИРОВАНИЙ К ТИПИЧНЫМ АНОМАЛИЯМ ЭЛЕКТРОПРОВОДНОСТИ В ТЕКТОНОСФЕРЕ Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
4
0
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
магнитотеллурические зондирования / геоэлектрическая модель / синтетические данные / аномалии электропроводности / magnetotelluric soundings / resistivity model / synthetic data / conductivity anomalies

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Данила Денисович Попов, Павел Юрьевич Пушкарев

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Данила Денисович Попов, Павел Юрьевич Пушкарев

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

SENSITIVITY OF MAGNETOTELLURIC SOUNDINGS TO TYPICAL ELECTRICAL CONDUCTIVITY ANOMALIES IN THE TECTONOSPHERE

A three-dimensional resistivity model is presented that contains electrical resistivity anomalies of a complex shape on three structural levels: in the sedimentary cover, in the consolidated earth’s crust and in the upper mantle. The synthetic (model) magnetotelluric data obtained as a result of solving the forward problem are analyzed to identify anomalies associated with resistivity inhomogeneities at various depths. Quantitative estimates of the sensitivity of various components of magnetotelluric data to these inhomogeneities are given. Further, we plan to use the presented data to evaluate the effectiveness of various methods for solving inverse problems.

Текст научной работы на тему «ЧУВСТВИТЕЛЬНОСТЬ МАГНИТОТЕЛЛУРИЧЕСКИХ ЗОНДИРОВАНИЙ К ТИПИЧНЫМ АНОМАЛИЯМ ЭЛЕКТРОПРОВОДНОСТИ В ТЕКТОНОСФЕРЕ»

УДК 550.837.211

doi: 10.55959/MSU0579-9406-4-2023-63-6-134-143

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

Данила Денисович Попов1 Павел Юрьевич Пушкарев2

1 Московский государственный университет имени М.В. Ломоносова, Москва, Россия; crossbrian97@mail.ruH

2 Московский государственный университет имени М.В. Ломоносова, Москва, Россия; pavel_pushkarev@list.ru

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

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

Для цитирования: Попов Д.Д., Пушкарев П.Ю. Чувствительность магнитотеллурических зондирований к типичным аномалиям электропроводности в тектоносфере // Вестн. Моск. ун-та. Сер. 4. Геология. 2023. № 6. С. 134-143.

SENSITIVITY OF MAGNETOTELLURIC SOUNDINGS TO TYPICAL ELECTRICAL CONDUCTIVITY ANOMALIES IN THE TECTONOSPHERE

Danila D. Popov1 Pavel Yu. Pushkarev2

1 Lomonosov Moscow State University, Moscow, Russia; crossbrian97@mail.ruH

2 Lomonosov Moscow State University, Moscow, Russia; pavel_pushkarev@list.ru

Abstract. A three-dimensional resistivity model is presented that contains electrical resistivity anomalies of a complex shape on three structural levels: in the sedimentary cover, in the consolidated earth's crust and in the upper mantle. The synthetic (model) magnetotelluric data obtained as a result of solving the forward problem are analyzed to identify anomalies associated with resistivity inhomogeneities at various depths. Quantitative estimates of the sensitivity of various components of magnetotelluric data to these inhomogeneities are given. Further, we plan to use the presented data to evaluate the effectiveness of various methods for solving inverse problems.

Keywords: magnetotelluric soundings, resistivity model, synthetic data, conductivity anomalies

For citation: Popov D.D., Pushkarev P.Yu. Sensitivity of magnetotelluric soundings to typical electrical conductivity anomalies in the tectonosphere. Moscow University Geol. Bull. 2023; 6: 134-143. (In Russ.).

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

Метод МТЗ применяется с середины XX века, в основе интерпретации данных тогда лежала горизонтально-слоистая или одномерная (Ш) модель среды. Благодаря развитию вычислительной техни-

ки, методов численного моделирования и решения обратных задач с конца XX века применяется двухмерная (2D), а с начала XXI века — трехмерная (3D) инверсия. Особенностям магнитотеллурического (МТ) поля в 3D средах посвящены работы М.Н. Бер-дичевского, В.И. Дмитриева, М.С. Жданова, Г. Хо-манна, Ф. Ваннамейкера и других ученых. Однако многие вопросы методики 3D интерпретации остаются слабо изученными, и нашей целью является их решение с использованием синтетических данных.

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

Рис. 1. Разрез вдоль центрального меридионального профиля — осадочный чехол

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

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

Для расчета МТ поля численными методами модель среды (геоэлектрическую модель) можно представить в виде массива прямоугольных ячеек, число которых по координатным осям составляет NX, NY и NZ (оси X и Yнаправлены условно на север и восток, ось Z — вниз). Каждая ячейка характеризуется значением удельного электрического сопротивления (УЭС) и размерами по трем осям. Вычисления выполнялись с помощью широко применяемой программы MT3DFwd [Mackie et al., 1993], основанной на методе конечных разностей.

Построение модели для программы MT3DFwd при больших размерах сетки (NX, NY и NZ) ручным способом весьма трудоемко. Чтобы избежать данной проблемы, нами была разработана и реализована в системе MATLAB программа GMC (Geophysical Models Constructor), позволяющая автоматизировать процесс построения 3D модели, а также визуализировать ее и конвертировать в необходимый формат. На вход программы GMC подаются: количество и размеры ячеек, фоновый слоистый разрез, границы которого задаются в виде карт, а также неоднородности, для которых также задаются карты кровли

и подошвы. На выходе программа GMC выдает «куб сопротивлений» на заданной сетке в формате программы MT3DFwd.

Нами использовалась сетка с размерами NX = 100, NY = 100 и NZ=79. По горизонтали размеры ячеек составили 4x4 км в центральной области размером 90x90 ячеек, по боками были добавлены по 5 дополнительных ячеек с каждой стороны размерами 8, 16, 32, 64 и 128 км, чтобы обеспечить выход на фоновые значения и выполнение граничных условий. По вертикали размеры ячеек постепенно увеличивались от 10 м до 30 км.

Рассмотрим структуры различного УЭС подробнее. Осадочный чехол имеет фоновую мощность 420 м и образует две структуры широтного простирания — поднятие в северной части с мощностью чехла до 100 м и прогиб в южной до 1000 м (рис. 1). При этом в чехле выделяется три слоя, мощность верхнего остается неизменной, а среднего и нижнего меняется пропорционально общей мощности. Размер центральной плоской части поднятия и прогиба составляет 32x80 км (по X и Y), внешний размер — 128x176 км. Центры структур находятся в серединах северной и южной половин планшета.

В консолидированной коре находятся две проводящие призмы. Северо-западная призма характеризуется азимутом простирания 40° на северо-восток, а юго-восточная простирается на 30° на северо-восток (рис. 2). Примерный размер каждой призмы составляет 180 км по длине, 60 км по ширине и мощностью в 20 км. Глубина до кровли северозападной призмы 20 км, юго-восточной — 10 км (рис. 3). В северо-восточной части планшета расстояние между призмами минимально и составляет 16 км, а в юго-западной части структуры максимально расходятся. Сопротивление северо-западной призмы составляет 60 Ом-м, юго-восточной — 30 Ом-м.

360

300

240

ЦМ

S

I

о

180

120

60

О км

2500 Ом«м

60 Ом»м \ 30 Ом«м

цш

Окм

I

60

120

180 з-в

240

300

360

Рис. 2. План — консолидированная кора, содержащая проводящие призмы

Таким образом, юго-восточная призма является более контрастной и залегает ближе к поверхности.

Кровля астеносферы в центральной части планшета поднимается с глубины 120 км до глубины 60 км (рис. 3). Размер центральной плоской части поднятия составляет примерно 70x170 км, краев — 170x270 км. Азимут простирания оси поднятия ориентирован на 60° на северо-восток (рис. 4).

Расчет магнитотеллурических данных. В результате решения прямой 3D задачи с помощью программы MT3DFwd определяются компоненты МТ поля в центрах ячеек на земной поверхности для двух поляризаций источника на заданных периодах (в нашем случае T _ 0,1, 0,2, 0,5, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000, 2000, 5000 с).

Для вычисления передаточных функций и инвариантных параметров, а также для их сохранения в файлах различных форматов нами была разработана и реализована в среде MATLAB программа MTDA (MagnetoTelluric Data Analysis). Она вычисляет:

• Передаточные функции: компоненты тензора импеданса [Z], матрицы типпера [W], магнитного тензора [M] и теллурического тензора [Г];

• Инвариантные параметры: эффективное кажущееся сопротивление peff> фазу эффективного импеданса 9eff, параметр неоднородности N, параметр асимметрии Свифта skewS, параметр асимметрии Бара skewB, норму матрицы ||W||, магнитовари-ационный параметр асимметрии skewMV, норму матрицы ||M||, магнитный параметр асимметрии skewM, норму матрицы || Г||, теллурический параметр асимметрии skewT.

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

• В текстовых файлах для визуализации в виде кривых зондирования, частотных разрезов и карт с помощью программ Grapher, Surfer (компания Golden Software) и других;

• В файлах формата EDI (Международный стандарт хранения МТ данных);

• В формате программы 3D инверсии МТ данных WSInv3DMT (автор W. Siripunvaraporn).

Анализ тензора импеданса. Основной передаточной функцией, используемой в методе МТЗ, является тензор импеданса. Импеданс имеет физический смысл комплексного сопротивления, которое оказывает среда индуцированному в ней электрическому току [Ковтун, Успенский, 2008]. Тензор импеданса состоит из четырех компонент, связывающих комплексные амплитуды горизонтальных компонент электрического и магнитного полей [Бердичевский, 1960; Cantwell, 1960]:

E _z H + Z H

Ey _ ZyxHx + ZyyHy.

(1)

Для качественного анализа результатов решения прямой задачи для рассмотренной выше модели в первую очередь были выбраны наиболее широко используемые инварианты — рей и фей. Эф-

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

60000 м

120000 лл

360

360

? 180

0 км

180 240 300 360 Окм 60 120 180

з-в з-в

Рис. 4. План — верхняя мантия и астеносферное поднятие. Слева — глубина 60 км, справа — 120 км

цш

360

фективный импеданс Ее{{ является инвариантным параметром (не меняющимся при повороте системы координат). Данный параметр определяется по формуле [НоЬЬз, 1992]:

Z ff =Jdet[Z] =.[Z~Z~^Z~Z~

eii V L J xx yy xy y

(2)

На рис. 5 и 6 представлены амплитудные и фазовые кривые в точках центральных меридионального и широтного профилей. На этих кривых можно выделить три основных слоя. В случае меридионального профиля на амплитудных кривых первый проводящий слой отмечается на Т < 1 с. Второй высокоомный слой выделяется на Т ~ 100 с. Наличие третьего участка, на котором отмечаются нисходящие ветви кривых, указывает на наличие проводящего основания.

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

Кроме того, наблюдается статическое смещение низкочастотных ветвей амплитудных кривых друг относительно друга на разных точках наблюдения. Точка х30у50 (названия точек образованы номерами ячеек по горизонтальным осям) расположена над проводящим прогибом в осадочном чехле, а потому кривая смещена вниз, а точка х70у50 — над высоко-омным поднятием, из-за чего кривая смещена вверх.

Что касается широтного профиля, то расхождение между амплитудными кривыми наблюдается с Т = 1 с. В случае фазовых кривых аналогичная картина начинает проявляться с Т = 0,2 с. По кривым, приуроченным к точкам наблюдения х50у30, х50у50 и х50у70, можно предположить, что различие в их

амплитудах объясняется положением точек относительно проводящих призм. По априорным данным видно, что точка х50у70 располагается над юго-восточной проводящей призмой с УЭС в 30 Ом-м, в то время как точка х50у30 располагается над более глубоко залегающей северо-западной призмой с УЭС в 60 Ом • м.

По картам рей (рис. 7) отмечаются аномалии, которые хорошо коррелируются с положениями поднятия и прогиба в осадочном чехле. С Т = 20 с на карте выделяется аномалия, вызванная наличием юго-восточной проводящей призмы, простирающейся в субмеридиональном направлении с азимутом 30°. Наличие северо-западной призмы с азимутом простирания 40° начинает отчетливо проявлять себя на картах с Т = 500 с.

Аналогичная ситуация наблюдается на картах Фей (рис. 8). Следует отметить, что эффект от проводящих призм начинает прослеживаться на меньших периодах. В данном случае это видно на карте при Т = 1 с. Особенно данный эффект проявляется от юго-восточной призмы на Т = 1 и 500 с. На Т = 5000 с можно также видеть возможное влияние от кровли проводящего поднятия астеносферы, несмотря на экранирующий эффект от вышезалегающих структур. Важно также отметить, что по мере увеличения периода уменьшается диапазон значений фей, и на периодах больше 1000 с при наличии сильного шума данные могут быть не информативными.

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

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

2 2

О юо -£ О.

од

10

100

VT, с1'2

-20

-ВО

it 9-1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-60

/ / \ .А, \ / / \ / \ \ У ч\

/ / / / \ \\ \ Л •

/ / / — Л N Л \

,* / / \\ \ Ч ^ 4jt V'. ¿Г

/ /

. 1 1 1 1 111 . 1 1 1 1 111 1 1 1 1 . 1.11

0,1

Рис. 5. Амплитудные и фазовые кривые Zeff, центральный меридиональный профиль

ю

лТГ, С1'2

л:15у50 х30>'50 jc50y50

юо

■х70у50 ■х85у50

1000

S 2

О юо

£ О)

О.

5=

9-"

> \

А д ч\

ч\ \\ «

V

1 1 1 1 1 111 ........ 1 1 1 1 1 1 111

л/Т, с1'2

Рис. 6. Амплитудные и фазовые кривые Zeff, центральный широтный профиль

VT, с1'2

х50у15 х50у30 x50v50

■х50/70 ■ х5()>>85

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

Анализ типпера (матрицы Визе-Паркинсона). Матрица типпера связывает вертикальную компоненту магнитного поля с горизонтальными, что дает возможность надежно определять положение глубинных протяженных проводящих тел [Wiese, 1962; Parkinson, 1959]:

Hz = WZXHX + WzyHy.

(5)

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

По построенным картам нормы || Ш|| на периодах до 20 с отмечаются аномалии, которые хорошо коррелируются с положениями поднятия и прогиба в осадочном чехле (рис. 9). С Т = 20 с на карте отчетливо выявляется аномалия от юго-восточной

ОД с

моооо-

-50000-

реп. Ом и

100000<

I

1000 562 316 178 100 56,2 31,6 17,8

-50000'

-100000 -50000 0 50000 100000

у, м

500 с

50000-

-100000 -50000

50000 100000

у, м

Рис. 7. Карты аномалий рей

20 с

1000

562

316

178

100

56,2

31,6

17,8

реп. Ом м

100000'

I

1000 562 316 178 100 56,2 31,6 17,8

ш

-50000'

50000 0 50000 100000

У'М 5000 с

рсВ. Ом и

-ЮОООО -50000 О

У, м

РсК Ом м

50000 100000

50000-

Р-'П. Ом м

1000 562 316 178 100 56,2 31,6 17,8

-100000 -50000 0 50000 ЮОООО

у, м

Поднятие и прогиб в осадочном чехле (40Э м)

Призм »1

Поднятие

астенос$>ерм (60 тыс. м)

Поднятие астеносоеры (120 так. м)

20 с

500 с

V- ; • • ж

■ • > у /Ш

а ^ЛЯ ' ш

-1-1-1--¥-

■100000 -50000 0 50000 юоооо

у, м

о 50000 100000

у, м

Рис. 9. Карты аномалий

20 с

Г./

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

У

V

. .,'Ж/ Г?

> 1 >

^ • /•

• ' г

/ . • '

■ / к '

у, М

1,64

1,61

1,58

1,55

1,52

1,49

1,46

1,43

1,4

1,37

1,34

1,31

1,28

1,25

1,22

Подня-ие и про'и 6

5 осадочном чехле ----

(4С0 м)

500 с

юсооо-^

50000

-5СОООН

-юсоос-^

V У./.

К \.у:

N " V \ •• •

л\ ... л

г \ \

. , . ■ Ч-

. ..1

<-'■ ■ Л » > / • • / ^ ^ / ■ .'у / Л

УГ . К <

-У у /

г

/ /.

. / . . ✓ / X /

-ЮОООО -50000 0 50000 ЮОООО

у, М

Рис. 10. Карты аномалий ||М||

Поднятие и прогиб в осадочном чехле (4С0 м}

5000 с

0,105

0,0975

0,09

0,0825

0,075

0,0675

0,06

10,0525

0,045

0,0375

0,03

0,0225

0,015

\ . \

\

^./. р/ ж>7/ / :.

• г'■ ■

■юоооо -50000 о 50000 юоооо

у, м

• • Призмы _ _ _

Поднятие астечосоеры ¡60 тыс. м)

Поднятие астечосоеэы (120 тыс. м)

5000 с

1,49

1,485

1,48

1,475

1,47

1,465

1,46

1,455

1,45

1,445

1,44

1,435

1,43

1,425

1,42

1,415

1,41

ЮОООО

50000

-50000

V г V . . .

\

\

__, У ,

. - \ >ч-" ' ■ \

V ' • ' \ . . > . ,

¥■/■ ■ ^ А ■ и

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

V ^ / > / •

М \ I /

/ _ » ч

/

' (я

/

-юоооо -50000 0 50000 юоооо

у М

1,435

1,4335

1,432

1,4305

1,429

1,4275

1,426

1,4245

1,423

1,4215

1,42

1,4185

1,417

1,4155

1,414

.. ПризV.сI - - —

Подня-ие астечосоеры _ ¡60 тыс. м)

Поднятие астеносоеэы (120 тас. м)

призмы. В данном случае, следует отметить, что для такой вытянутой структуры довольно точно указывается местоположение оси тока, так как над ней || Ш|| стремится к нулю. Также с ростом периода сильнее проявляется аномалия от северо-западной призмы, несмотря на ее искажение, вызванное влиянием высокоомного поднятия в осадочном чехле.

Анализ магнитного тензора. Другой передаточной функцией, определяемой только по магнитному полю, является (горизонтальный) магнитный тензор. Он связывает горизонтальные поля в рядовой и базовой (В) точках наблюдения [Варенцов и др., 2003]:

Нх = ИХХИХВ + МхуИуБ,

Ну = МухНХ + MyyHyB,

(6)

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

В 2D среде может использоваться компонента, связанная с токами, поляризованными вдоль структур, а в 3D среде все четыре компоненты [М] являются информативными. В отличие от [Ш], [М] чувствителен к структурам, непосредственно расположенным под точкой наблюдения.

По построенным картам ||М|| можно также выделить аномалии, которые отчетливо коррелируются с положениями поднятия и прогиба в осадочном чехле (рис. 10). Кроме того, с Т = 20 с на карте так же довольно однозначно отмечается аномалия от юго-восточной призмы. При этом аномалия от северозападной призмы усиливается по мере увеличения периода. На периоде 5000 с, вероятно, проявляется и поднятие проводящей астеносферы.

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

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

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

Анализируя полученные результаты, можно видеть, что поднятие и прогиб фундамента для всех параметров выявляются практически на всем представленном диапазоне периодов (табл. 1). Исключением является Т = 5000 с, на котором такие параметры, как фей, ||Ш|| и ||М||,уже не чувствительны к данным приповерхностным структурам.

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

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

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

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

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

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

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

Таблица 1

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

Поднятие и прогиб

T = 0,1 с T = 1 с T = 20 с T = 500 с T = 5000 с

APeff % 100 65 22 17 20 APeff % < 1 < 10

Aqf ° 32 19 13 1,9 0,26 Aqf ° < 0,5 < 5

A||W|| 0,095 0,21 0,27 0,046 0,0095 A||W|| < 0,03 < 0,1

A||M|| 0,054 0,145 0,21 0,016 0,0042 A||M|| < 0,03 < 0,1

A||T|| 3,4 3,4 1,1 0,9 0,9 A||T|| < 0,03 < 0,1

Таблица 2

Максимальные чувствительности различных компонент МТ данных к юго-восточной призме

Ю-В призма

0,19 0,8 T = 20 с T = 500 с T = 5000 с

APeff, % 20 16 20 APeff, % < 1 < 10

Aqf ° 0,26 8,5 1,7 0,85 Aqf ° < 0,5 < 5

A||W|| 8,50E-05 0,0145 0,25 0,075 0,015 A||W|| < 0,03 < 0,1

A||M|| 2,90E-05 0,0075 0,18 0,034 0,008 A||M|| < 0,03 < 0,1

A||T|| 3,40E-03 0,025 0,54 0,5 0,52 A||T|| < 0,03 < 0,1

Таблица 3

Максимальные чувствительности различных компонент МТ данных к северо-западной призме

С-З призма

T = 0,1 с T = 1 с T = 20 с T = 500 с T = 5000 с

APeff, % 0,032 1,3 9 6,5 8,5 APeff, % < 1 < 10

Aqf ° 0,0105 1,8 5,4 0,9 0,44 Aqf ° < 0,5 < 5

A||W|| 4,00E-06 0,0032 0,125 0,036 0,008 A||W|| < 0,03 < 0,1

A||M|| 0 0,0022 0,09 0,018 0,0044 A||M|| < 0,03 < 0,1

A||T|| 5,80E-04 0,036 0,28 0,23 0,25 A||T|| < 0,03 < 0,1

Таблица 4

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

Поднятие астеносферы

T = 0,1 с T = 1 с T = 500 с T = 5000 с

APeff, % 0,4 4,4 4,4 APeff, % < 1 < 10

Aqf ° 0 0 1,25 0,45 Aqf ° < 0,5 < 5

A||W|| 0 0 0,025 0,046 0,018 A||W|| < 0,03 < 0,1

A||M|| 0 0 0,011 0,048 0,018 A||M|| < 0,03 < 0,1

A||T|| 0 0 0,028 0,125 0,125 A||T|| < 0,03 < 0,1

СПИСОК ЛИТЕРАТУРЫ

1. Бердичевский М.Н. Основы теории магнитотеллу-рического профилирования // Прикладная геофизика. 1960. Вып. 28. С. 70-91.

2. Варенцов И.М., Соколова Е.Ю., Мартанус Е.Р., На-ливайко К.В. Система передаточных операторов электромагнитного поля для массива синхронных зондирований BEAR // Физика Земли. 2003. № 2. С. 30-61.

3. Ковтун А.А., Успенский Н.И. Геоэлектрика. Естественное поле. СПб., 2008. 122 с.

4. Cantwell T. Detection and analysis of low frequency magnetotelluric signals. Ph.D. Thesis. Cambridge, 1960.

5. Hobbs B.A. Terminology and symbols for use in studies of electromagnetic induction in the Earth // Surveys in Geophysics. 1992. Vol. 13. P. 489-515.

6. Mackie R.L., Madden T.R., WannamakerP.E. Three-dimensional magnetotelluric modeling using difference equations — theory and comparison to integral equation solutions // Geophysics. 1993. Vol. 58. P. 215-226.

7. Parkinson W.D. Direction of rapid geomagnetic fluctuation // Geophysical Journal of the Royal Astronomical Society. 1959. Vol. 2. P. 1-14.

8. Wiese H. Geomagnetische tiefentellurik, Teil 2, Die streichrichtung der untergrundstrukturen des elektrischen winderstandes, erschlossen aus geomagnetischen Variationen // Geofisica Pura e Applicata. 1962. Vol. 52. P. 83-103.

Статья поступила в редакцию 21.08.2023, одобрена после рецензирования 02.10.2023, принята к публикации 22.01.2024

i Надоели баннеры? Вы всегда можете отключить рекламу.