Магнитотеллурическое зондирование

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

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

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

Пульсации относятся к классу короткопериодных колебаний МТ-поля (КПК). Их частотный спектр носит случайный характер и укладывается в диапазоне частот 10-²—10 Гц с максимумом периода в интервале 10—60 с. Интенсивность (амплитуда) КПК для магнитного поля в средних широтах достигает нескольких нанотесл; для электрического поля — нескольких милливольт на километр.

Среди КПК выделяют устойчивые пульсации Р_с, представляющие совокупность квазисинусоидальных вариаций, непрерывно следующих друг за другом (рис. 4.30, а), и иррегулярные пульсации Р, проявляющиеся цугами (пачками) колебаний, разделенными во времени периодами относительно спокойного поля (рис. 4.30, б).

Бухтообразные возмущения DР-вариации) характерны единичными проявлениями, схожими по форме магнитограмм с береговыми очертаниями бухт (рис. 4.30,в). Их период — первые десятки — первые сотни минут. В высоких широтах амплитуды DР- вариаций достигают сотен нанотесл; в средних широтах — десятков нанотесл.

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

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

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

В каждой точке поверхности земли векторы Е и Н

магнитотеллурического поля изменяются со временем не только по величине, но и по направлению.

(Ksyj-JX/^s^r^------- ^^/^в'

 

Рис. 4.30. Типы вариаций магнитотеллурического поля: а - устойчивые,б -

нерегулярные, в - бухтообразные.

Проекции концов векторов Е и Н на плоскость XOY в разные моменты времени называются годографами. Годографы векторов Е и Н представляя собой сложные

фигуры — изомерные (нелинейная поляризация) либо сильно вытянутые в каком-либо направлении (квазилинейная поляризация). На рис.4.31 приведены годографы обоих типов.

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

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

(4.48)

sin d k₁

где k1 и k2 — волновые числа для воздуха и горной породы.

Волновое число среды зависит от круговой частоты поля (ω = 2πf), магнитной проницаемости μ, проводимости γ и диэлектрической проницаемости ε среды:

(4.49)

При достаточно низких частотах поля k1 = 0. Отсюда следует, что sin d «О, т. е. <з?«0.

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

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

 

Рис. 4.31. Годографы векторов Еи Н магнитотеллурического поля при нелинейной (а)

и квазилинейной (б) поляризации

В качестве основной характеристики плоской электромагнитной волны и степени воздействия на нее проводящей среды обычно пользуются отношением двух взаимно перпендикулярных компонент поля (например Ех к Нy), которое называется импедансом (Z = Ех/Ну):

fJ^ (4.50)

1/

Магнитотеллурическое зондирование. Процесс магнитотеллурического

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

Рис. 4.32. Схема установки для магнитотеллурических наблюдений

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

Матнитотеллурическое зондирование представляет собой одну из модификаций частотного зондирования, основанную на изучении вариаций естественного электромагнит­ного поля Земли в широком диапазоне периодов. Сущность МТЗ заключается в одновременной регистрации компонентов магнитотеллурического поля Е_x, Е_y, Н_x, Н_y и Н_z на поверхности Земли и последующем спектральном анализе результатов измерений.

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

Измерения выполняют установкой, расположенной на поверхности земли. Она состоит из двух взаимно перпендикулярных датчиков электрического поля — приемных линий М1N1, М2N2, в которых используются обычные неполяризующиеся электроды (как в методе ЕП) и трех магнитометров-вариометров Н_x,Н_у, Н_z. Размещение датчиков поля показано на рис. 4.32. Длина приемных линий обычно составляет 200—500 м.

Направление приемных линий и расположение магнитометров выбирают строго в соответствии с основными элементами залегания горных пород и тектоникой района. Одну из измерительных линий М1N1 вытягивают вдоль простирания пород для изучения продольной составляющей электрического поля, другую — вкрест простирания для изучения поперечной составляющей поля. Направление, близкое к меридиональному, условно обозначают у, перпендикулярное к ему — х. Соответственно размещают и магнитометры. Расстояние между ними должно быть не менее 8 м, а вся измерительная установка (центр зондирования) находится на удалении 25—30 м от измерительной лаборатории.

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

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

Входной импеданс слоистого разреза связан с параметрами этого разреза соотношением

К

 

(4.51)

где R_п — т.н. приведенный импеданс n-слойного разреза,; со — частота поля; /л^ —

магнитная проницаемость; k1 — волновое число верхнего слоя.

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

+... +arcth

Исп_А_х + arcth

R_n = cth\ ik1h₁ + arcth I— x cth VA

к (4.52)

Для однородной среды (h 1 оо) R_n = 1, т. е.

К

(4.53)

Если среда немагнитна, т. е. ее магнитная проницаемость равна Ал • 10 ⁷ Гн/м, а удельное сопротивление выражено в ом-метрах, модуль импеданса на поверхности однородного полупространства

ПГ

(4.54)

где р_х — удельное сопротивление полупространства, Т - период измеряемого сигнала в сек. Отсюда

_п _ ут 7² (л zsr\

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

Рт-

P_T=2T\Z\ =2TEJH В соответствии с выражениями (4.51) и (4.54)

(4.56)

(4.57)

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

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

Т, (4.58)

Выражение для k1 можно преобразовать следующим образом:

(4.58)

Тогда выражение выражение для ρ_Т будет зависеть от параметров:

(4.59)

Это соотношение определяет конструкцию палеток теоретических кривых МТЗ.

Из теории МТЗ следует, что ρ_т — величина комплексная, т. е. определяется амплитудой и фазой. Соответственно рассчитываются палетки удельного сопротивления р_т и разности фаз между Е_х и Н_у (р_рТ.

Палетки амплитуд р_т строят в двойном логарифмическом масштабе с модулем 6,25

Рт

см, причем по оси ординат откладывают

, а по оси абсцисс — λ1/h1. Фазовые кривые

строят в полулогарифмическом масштабе, причем по оси ординат в арифметическом масштабе откладывают (р_рТ, а по оси абсцисс— в логарифмическом масштабе λ1/h1. На рис. 4.33, а, б

изображены палетки для интерпретации амплитудных и фазовых кривых МТЗ.

Внешне характер амплитудных кривых МТЗ сходен с кривыми зондирования методом сопротивлений, т. е. горизонтам с высоким сопротивлением отвечают максимумы на графиках р_т, горизонтам с низким сопротивлением — минимумы. Четкость, с которой отдельные горизонты проявляются на графиках, зависит от их мощности и контрастности по удельному сопротивлению. Характерная особенность графиков — наличие дополнительного экстремума в высокочастотной части кривой и переход величины р._г через нулевое значение при λ1/h1 = 8. В точке с абсциссой 8 и ординатой 1 располагается крест палетки.

Для кривых МТЗ, так же как и для кривых зондирований методом сопротивлений, справедлив принцип эквивалентности, т. е. при определенных соотношениях параметров разреза кривые МТЗ совпадают в пределах заданной погрешности.

Рис. 4.33. Палетки для интерпретации данных МТЗ: а, б — палетки соответственно амплитудных и фазовых кривых для двухслойного разреза, в, г — амплитудные кривые для трехслойных разрезов типов Ни К.

Для разрезов типов Н, А справедлива эквивалентность по S2 (так же как и для ВЭЗ). Для разрезов типов К и Q характерна эквивалентность по h2, т. е. при изменении в некоторых пределах h2 кривые МТЗ в разрезах таких типов практически совпадают. Из этого следует принципиальная возможность однозначного определения h 2 по кривым МТЗ К и Q с погрешностью, зависящей от области применимости принципа эквивалентности. Для установления этих границ Б. К. Матвеевым построены специальные номограммы, сходные по конструкции с аналогичными номограммами Пылаева.

Интерпретация двухслойных кривых МТЗ с помощью палеток заключается в том, что полевую кривую, построенную на прозрачном билогарифмическом бланке, накладывают на палетку и наилучшим образом совмещают с одной из теоретических кривых. При таком совмещении ордината креста палетки равна ρ1 а его абсцисса Т связана с h1 соотношением:

(4.60)

В настоящее время обработка данных измерений производится на компьютерах по специальным программам. В программах предусматриваются узкополосная фильтрация записей магнитотеллурического поля, выделение гармонических составляющих Е_x, Еy, Н_x, Н_y и Н_z для заданной последовательности периодов Т, вычисление импедансов Z_xy=£V Н_у\ Zyx=Еy/Н_x; Z_xx=Еx/Нx_; Zyy=Еy/Н_y, определение сдвига фаз между взаимно перпендикулярными

 

 

 

 

составляющими	Ех	и Ну,	Еу	И	Нх,	a Ti
формулам:
		Ртш		Т	Ех	/я/

также вычисление кажущихся сопротивлений по

Дополнительные импедансы Zxx, Zyy и вертикальную составляющую Н_z магнитного поля используют для анализа неоднородности среды и выбора соответствующей методики интерпретации результатов наблюдений.

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

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

Контрольные вопросы.

1. Физическая сущность метода естественного электрического поля (ЕП). Способы
измерений и условия их применения.

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

3. Физическая сущность и установки в методе электропрофилирования.

4. Чему равно значение потенциала в точке М от электродов А и В? (записать формулу).
Покажите токовые и эквипотенциальные линии от электродов А и В.

5. Нарисуйте вид кривой ВЭЗ, если геоэлектрический разрез имеет вид: pi=120 Ом*м,
Р2=50 Ом*м, рз=200 Ом*м, р4=80 Ом*м. Какой это тип кривой и разрез?

6. Физическая сущность метода и установки, применяемые в методе ВЭЗ.

7. Покажите, какой вид имеют аномалии естественного поля (ЕП) над скоплениями
сульфидов (съемка способом потенциала). Как определить, где находится рудный
объект?

8. В чем разница между истинным и кажущимся удельными электрическими
сопротивлениями?

9. Что понимается под удельным электрическим сопротивлением и удельной
проводимостью в электроразведке и каковы единицы измерения этих величин?

10. Запишите формулу для определения разности потенциалов между электродами MN от
электрода А (электрод В можно считать удаленным на бесконечность). Покажите
токовые и эквипотенциальные линии от электрода А.

11. В каких единицах измеряется значение коэффициента установки? Покажите
почему?

12. Для симметричной установки АМNB c коэффициентом k = 3.14 м на точке получено
значение рк = 52 Ом м. Чему будет равно значение рк, если АМ=BN=40 м, MN=20 м,
V=300mB, J=0.002 A?

13. Определите, в соответствии с вашим вариантом, можно ли считать указанную в табл.1
установку трехэлектродной, если среднее удельное сопротивление среды равно рср, и
ток равен J. Питающий электрод можно считать удаленным на бесконечность, если
его потенциал в точке N не превышает 3 % значения потенциала в этой точке от
ближнего питающего электрода.

Таблица 1

Вариант	Установка	Рср Ом м	J, А
	A10M5N120B		0.08
	A50M20N400B		0.1
	A80M15N100B		0.07
	A60M30N80B		0.09
	A40M20N300B		0.05
	A20M6N100B		0,06

 

Вариант	I	II	III	IV
Установ -ка	A10M5N120B	A50M20N400B	A80M15N100B	A60M30N80B
Рев.
J	0.08	0.1	0.07	0.09

14. Сущность электрического профилирования и зондирования (что общее и в чем
отличие?).

15. Какие разрезы называются двухслойными, трехслойными и какими параметрами они
характеризуются?

16. Почему кривые ВЭЗ изображаются в билогарифмическом масштабе?

17. В каких единицах измеряется истинное и кажущееся удельное электрическое
сопротивление?

18. Что такое коэффициент установки и в каких единицах он измеряется?

19. Что понимается под качественной и что - под количественной интерпретацией
кривых ВЭЗ?

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

21. Какие параметры измеряются и определяются в методе частотного зондирования
(ЧЗ)? Чем метод ЧЗ отличается от ВЭЗ?

22. Поясните сущность магнитотеллурического зондирования. Чем здесь достигается
изменение глубины исследования? Чем метод МТЗ отличается от ЕП?

23. За счет чего возможно применение методов электроразведки при поисках
нефтегазовых залежей?

Глава 5. РАДИОМЕТРИЯ И ЯДЕРНАЯ ГЕОФИЗИКА