ИНТЕРПРЕТАЦИЯ ДАННЫХ
МАГНИТОТЕЛЛУРИЧЕСКОГО ЗОНДИРОВАНИЯ И ПРОФИЛИРОВАНИЯ
Общие замечания
Интерпретация кривых МТЗ - обзор подходов
Интерпретация кривых МТЗ по асимптотам и характерным
точкам
Обработка данных МТП
Пример интерпретации данных МТЗ/МТП
Порядок выполнения работы
Отчетные материалы
Общие
замечания
Методы
магнитотеллурического зондирования (МТЗ) и профилирования (МТП) основаны
на изучении естественного переменного электромагнитного поля Земли, имеющего
космическое происхождение. Методы МТЗ и МТП применяются при исследовании
рельефа поверхности высокоомных и низкоомных геоэлектрических горизонтов
на глубинах от 1 до 10 км и более.
Магнитотеллурические
зондирования представляют собой разновидность электромагнитных зондирований,
основанную на регистрации колебаний магнитотеллурического поля в широком
диапазоне периодов (от долей сек до нескольких часов и суток). Глубина
зондирования благодаря скин-эффекту определяется периодом регистрируемых
колебаний.
Информацию
о строении разреза несут горизонтальные компоненты электрического и магнитного
полей, измеряемые на поверхности Земли. По их отношению вычисляется входной
импеданс среды
![](image002.gif)
Зависимость
импеданса от периода колебаний изображают в виде амплитудных и фазовых
кривых МТЗ, являющихся частотными характеристиками изучаемого разреза.
В случае горизонтально-слоистой среды (ГСС) входной импеданс связан с параметрами
среды следующим образом
![](image004.gif)
где h
i - мощность i-слоя,
ri
- его сопротивление, а k
i - волновое число.
Амплитудные
и фазовые кривые зондирования в случае ГСС отражают изменение сопротивления
горных пород с глубиной. При построении амплитудных кривых используют величину
кажущегося (или эффективного) сопротивления. В практических единицах:
![](image006.gif)
Кривые
rT МТЗ,
как и
rK ВЭЗ,
строят на бланках с двойным логарифмическим масштабом. Для экспериментальных
кривых по оси ординат откладывают
rT,
а по оси абсцисс
![](image008.gif)
.
Для теоретических кривых по оси ординат откладывают
rT/
r1,
а по оси абсцисс
l1/h
1,
где
![](image010.gif)
-
длина волны в первом слое.
Магнитотеллурическое
профилирование (МТП), в отличие от МТЗ изучает вариации естественных
полей в небольшом диапазоне периодов (обычно от 10 до 80 сек). Этот диапазон
относится к наиболее регулярной во времени и пространстве части спектра.
Если МТП в узком диапазоне периодов соответствует асимптотическим значениям
кривой МТЗ, то из МТП можно получить информацию о таких параметрах разреза,
как суммарная продольная проводимость S (при высокоомном опорном горизонте)
или суммарная мощность Н (при низкоомном опорном горизонте). Изучая поле
в узком диапазоне периодов можно резко сократить время наблюдения вариаций,
сгустить сеть полевых точек. При комплексировании МТП с МТЗ в благоприятных
условиях удается получить важные количественные характеристики разреза,
свойственные МТЗ, и высокую детальность исследований по профилю, свойственную
МТП.
Методика
интерпретации кривых МТЗ сначала была разработана для горизонтально-слоистых
сред. Практические кривые МТЗ, полученные в условиях горизонтально-неоднородных
сред, отличаются от кривых для ГСС. Эти отличия называются искажениями
кривых МТЗ. При интерпретации полевых данных, неизбежно приходится иметь
дело с искаженными кривыми, и процент таких кривых может быть довольно
высок. Однако вопросы выявления и анализа искажений, качественной и количественной
интерпретации искаженных кривых выходят за пределы данной задачи.
Количественную
интерпретацию кривых МТЗ в рамках модели ГСС выполняют с помощью следующих
способов: 1) асимптотических; 2) алгебраических или дифференциальных трансформаций;
3) палеточных (см. рисунок); 4) по координатам экстремальных точек; 5)
методом подбора на ЭВМ.
Идея
первого способа состоит в том, что в случае опорного горизонта высокого
удельного сопротивления кривая rT
имеет восходящую ветвь, наклоненную под 63њ. Проводя касательную к этой
асимптоте по пересечению этой асимптоты с горизонтальными линиями rT=10
или rT=1, можно
определить величину суммарной продольной проводимости S до высокоомного
основания разреза.
Идея
второго способа состоит в преобразовании кривой
![](image012.gif)
в кривую, как можно более близкую к
rИСТИН=f
(H
ИСТИН) или S
ИСТИН = f (H
ИСТИН).
В
основе третьего способа лежит графическое совмещение практических кривых
МТЗ с теоретическими кривыми.. Многослойные кривые интерпретируют по частям,
используя принцип эквивалентных замен. Практические результаты интерпретации
могут быть получены в пределах действия принципа эквивалентности. Эти пределы
практической неоднозначности в оценке параметров промежуточных слоев могут
быть определены с помощью номограмм эквивалентности Б.К.Матвеева. Установлено,
что на переменном токе принцип эквивалентности проявляется в более узких
пределах, чем на постоянном токе. Кроме того, для разрезов типа K и Q вместо
принципа эквивалентности по Т2 на переменном токе действует
эквивалентность по h2, то есть значение сопротивления промежуточного
высокоомного слоя почти не влияет на оценку мощности этого слоя.
Четвертый
способ интерпретации основан на теоретических и экспериментальных связях
координат характерных точек кривых МТЗ с параметрами разреза. Подробнее
он будет рассмотрен ниже.
Интерпретация
на ЭВМ наиболее распространена в настоящее время. С помощью ЭВМ параметры
геоэлектрического разреза находятся путем минимизации функционала невязки:
![](image016.gif)
характеризующего среднеквадратичное отклонение экспериментальной кривой
rT
от модельной кривой
rTq.
Минимизацию выполняют, корректируя параметры модели q. Начальные значения
q (нулевое приближение) выбирают с использованием имеющейся геолого-геофизической
информации.
Интерпретация
кривых МТЗ по асимптотам и характерным точкам
1. Определение суммарной
продольной проводимости S надопорной толщи.
2. Определение удельного
сопротивления опорного горизонта и периодов ТМИН соответствующих
минимуму кривых rТ.
3. Уточнение границ главного
и дополнительного интервалов МТП.
4. В благоприятных условиях
по минимуму кривых rТ
определяется среднее продольное сопротивление rL
надопорной толщи.
5. Определение суммарной
мощности надопорной толщи.
Наиболее просто величина
S определяется по восходящей ветви кривых
rТ,
если она наклонена к оси абсцисс под углом 63њ (
rN=
?).
В этом случае:
![](image018.gif)
или
![](image020.gif)
,
где
![](image022.gif)
и
![](image024.gif)
- абсциссы точки пересечения линии S с линиями
rТ=1
Омћм и
rТ=10
Омћм. Если интерпретируемая кривая МТЗ имеет четкий минимум, то величина
S может быть приближенно оценена как
![](image026.gif)
,
(
Tмин,
rмин)
- координаты точки минимума кривой МТЗ
Если
восходящая ветвь кривой rТ
наклонена к оси абсцисс под углом 63њ, то полагают, что rN=?.
При меньших углах наклона восходящей ветви величину rN
оценивают путем количественной интерпретации.
[В
начало пункта]
Определение среднего продольного сопротивления rL
надопорной толщи.
В
благоприятных условиях rL
может быть определено по ординате минимума кривой МТЗ rL=Prмин,
где P - коэффициент, зависящий от соотношения параметров разреза (таблица
ниже).
h2/h1
|
r2/r1
|
P
|
1-2
|
1/9-1/39
|
0.825
|
2-5
|
1/4-1/39
|
1
|
5-10
|
2/3-1/39
|
1.15
|
>10
|
2/3-1/39
|
1.3
|
В
случае трехслойного разреза типа H может быть предложена более точная формула
![](image028.gif)
,
которая применима при известном
r2.
Значения q=1.2 при
n2=1-2
и q=1.15 при
n2>2.
Если
на части кривых минимум не выражен, величина
rL
может быть определена по корреляционной связи между нею и S (по графику,
представляющему зависимость
rL(S)).
Формула Гуммеля говорит о линейной связи между этим величинами. Вместе
с тем из-за логнормального распределения, которое, как известно, характерно
для величин, применяемых в электроразведке, вполне возможен какой-то вариант
линейной связи между логарифмами этих величин. Такую связь следует опробовать,
если предположение о линейной связи
rL(S)
приводит к противоречию с геологическими данными и здравым смыслом
[Возврат
к МТП].
В
случае опорного горизонта (основания разреза) высокого сопротивления, суммарная
мощность надопорной толщи определяется по формуле Гуммеля
![](image030.gif)
(размерность
- метры). Для определения H (размерность - километры) используются также
приближенные формулы Т.Н.Завадской
![](image032.gif)
или
Г.Д.Цекова
![](image034.gif)
,
где
r2 - сопротивление
второго слоя в разрезе типа H, а С - коэффициент, зависящий от
n2
и
m2: при
n?4,
С=1, при
n2<4,
1<C<2.2.
В
разрезах типа H (и А) эту величину находят по ниспадающей ветви кривой
МТЗ. Если сопротивление подстилающего (2-го) проводящего слоя равно 0,
то:
![](image036.gif)
,
где
![](image038.gif)
- абсцисса и ордината любой точки на ниспадающей под 63њ амплитудной кривой
МТЗ. В случае более пологой нисходящей ветви (
r2
не равно 0), мощность 1-го слоя можно оценить по двухточечной формуле Б.
К. Матвеева:
[В начало пункта]
Наибольшее
применение МТП находит при изучении разрезов с высокоомным опорным горизонтом
(трехслойные разрезы типа Н и A, или многослойные, сводящиеся к разрезам
типа Н и A), когда сопротивление основания значительно превосходит среднее
продольное сопротивление надопорной толщи.
В
этом случае основным параметром, получаемым в результате интерпретации
МТП является суммарная продольная проводимость S, определяемая по
главной
формуле МТП ![](image042.gif)
или
![](image044.gif)
(последнее - при непроводящем основании). Значение
rN
при необходимости определяется по опорным МТЗ. Главная формула МТП обеспечивает
точность определения S не хуже 10% в интервале периодов, существенно превышающих
Т
мин, а именно T/T
MIN?2.3,
то есть если период, на котором выполняются измерения МТП, лежит в области
правой асимптоты кривой МТЗ.
Кроме
главной формулы иногда используется
дополнительная формула МТП ![](image046.gif)
(при непроводящем основании), по которой величину S можно определить при
Т
МТП близком к Т
мин. Формула применима (см. рисунок)
если 1<T/T
мин<4.
По
значениям SМТП строят графики и карты, дающие представление
о тектонике исследуемого района. В случае высокоомного основания максимумы
S отвечают прогибам фундамента, а минимумы S - поднятиям.
Чтобы
определить участки профиля, на которых следует использовать ту или иную
формулу, можно построить по данным МТЗ на том же профиле график изменения
отношения T/Tмин по профилю.
Зная
среднее продольное сопротивление надопорной толщи и закономерности его
изменения на изучаемой площади, можно определить
глубины залегания опорного
горизонта (для разрезов типа Н и А) в каждой точке МТП по формуле Гуммеля.
Информацию о
rL
и его изменениях по площади получают из опорных МТЗ, или из структурных
скважин или данных сейсморазведки (в последних случаях определяется
rL
по известному Н и измеренному S) с использованием корреляционной связи
между
rL и S
(см. раздел
'Интерпретация данных МТЗ'). Кроме того,
для определения Н можно по данным МТЗ непосредственно установить корреляционную
связь между Н и S. Эта связь будет, скорее всего, линейной (см., впрочем,
замечания о связи
rL
и S в разделе
'Интерпретация данных МТЗ').
В
платформенном районе проведены рекогносцировочные исследования методами
МТЗ и МТП с целью прослеживания кровли палеозойского фундамента. Получено
5 точек МТЗ, отстоящих друг от друга на 50 км по профилю, ориентированному
с запада на восток. В промежутках между МТЗ с шагом 10 км выполнено 16
точек МТП на периоде Т=80 с (
![](image050.gif)
).
На
пикете 50 профиля МТЗ/МТП расположена опорная скважина, разрез по данным
бурения и каротажа - трехслойный, первый слой мощностью 1 км представлен
известняками и песчаниками с удельным сопротивлением 32 Омћм, второй -
толщей переслаивания песков и глин, насыщенной минерализованными водами,
мощностью 2 км, удельным сопротивлением 2 Омћм, а третий - высокоомными
(более 1000 Омћм) метаморфизованными породами палеозойского фундамента.
На основе геологических сведений о строении разреза можно предполагать,
что мощность первого слоя и сопротивления слоев в разрезе выдержаны в пределах
профиля.
Кривые
МТЗ и график импеданса МТП представлены на рисунке. Интерпретация кривых
МТЗ проводится по асимптотам и характерным точкам. (заметим, что эти кривые
НЕ могут быть подобраны в рамках горизонтально-слоистой модели).
Все
кривые имеют тип Н с выраженной правой асимптотой, составляющей угол 63њ
с осью периодов. Поэтому сопротивление основания разреза можно считать
бесконечно большим.
Суммарная
продольная проводимость определяется по периоду пересечения правой асимптоты
с уровнем
rТ=10
Омћм. Для определения среднего продольного сопротивления на точках 100,
150, 200, имеющих выраженный минимум, используются координаты этого минимума.
На точке 50 для этого используются данные по скважине. Расчет сделан по
формуле для известного по данным бурения
r2=2
Омћм. На точке 0 минимума нет, поэтому для нее среднее продольное сопротивление
устанавливается по линейной регрессии
rL(S)=2.404ћ10
-3ћS+0.721,
полученной по соответствующему
графику.
Мощность
надопорной толщи определена по формуле Гуммеля.
Мощность
первого слоя может быть установлена по кривым на точках 150 и 200, имеющим
наиболее выраженные минимумы, с применением двухточечной формулы Матвеева.
На пикете 50 мощность первого слоя известна по данным бурения. На пикетах
0 и 100 мощность первого слоя принята равной среднему по известным значениям.
Для
более детального отслеживания кровли основания разреза используются данные
МТП. Для расчета суммарной продольной проводимости по данным МТП использованы
главная и дополнительная формулы МТП. Соотношение периодов
МТП и минимума кривой МТЗ меняется по профилю. На фрагменте 0-80 ТМТП/Тмин>4,
здесь обязательно применение главной формулы МТП. На фрагменте 140-200
ТМТП/Тмин<2.3, здесь обязательно применение дополнительной
формулы МТП. На остальной части профиля применимы обе формулы, при расчете
использована главная формула.
Оценка
глубины кровли основания разреза (мощности надопорной толщи) выполнена
по линейной регрессии Н(S)=7.635ћS-3901, полученной по
графику
Н(S) по результатам интерпретации данных МТЗ.
В
основании разреза залегают метаморфические породы удельным сопротивлением
свыше 1000 Омћм. Кровля основания падает с запада на восток от глубины
около 1 км до глубины 15 км. На фрагменте 40-100 имеется сводообразное
поднятие кровли от глубины 6-7 км на краях фрагмента до глубины 3-3.5 км
на пикетах 60-70. Незначительное поднятие отмечается также на пикетах 140-150
(только по данным МТП).
На
основании залегает песчано-глинистый слой удельным сопротивлением 2 Омћм.
Низкое сопротивление объясняется насыщенностью пород слоя минерализованной
водой. Мощность слоя изменяется от сотен метров на западе до 13-14 км на
востоке. На фрагменте 40-100 (поднятие кровли основания) мощность слоя
сокращается от 4-6 км по краям до 2-2.5 км в середине фрагмента над самым
поднятым участком кровли.
С
поверхности разрез перекрыт слоем известняков и песчаников сопротивлением
32 Омћм мощностью 1-1.5 км.
Порядок
выполнения работы
1.
Построить кривые МТЗ и график импеданса МТП по профилю.
2.
Определить S надопорной толщи.
3.
Определить rL
по кривой у скважины.
4.
ОпределитьrL
по кривым с выраженным минимумом. Значение r2
взять по скважине.
5.
Определить зависимость rL
от S , определить по ней rL
для кривой МТЗ без минимума.
6.
Определить суммарную мощность надопорной толщи по формуле Гуммеля.
7.
Оценить h1 по кривым с выраженной левой ниспадающей ветвью h1
по формуле Матвеева
8.
Построить график ТМТП/Тмин> по профилю, установить
интервалы применения формул МТП, оценить S по данным МТП, построить график
S по профилю.
9.
По рассчитанным S определить одним из способов мощность надопорной толщи.
10.
Построить геоэлектрический разрез по результатам интерпретации.
Отчетный
материал
Графика - кривые МТЗ,
графики импеданса и S по МТП, графики, использованные для установления
связей между величинами, геоэлектрический разрез.
Пояснительная записка.
Пример оформления графики
и образец текста пояснительной записки приведены в разделе 'Пример интерпретации
данных МТЗ/МТП'.
[В начало]