Указания к выполнению лабораторных работ по курсу «Комплексирование геофизических методов»
Гравитационный и магнитный потенциалы одного и того же источника связаны функционально. В случае постоянной плотности и однородной намагниченности эта связь описывается уравнением Пуассона, которое в случае двумерных тел, можно записать в виде:
(1)
или
(2)
где U – потенциал магнитного поля, Vz, и Vx – производные гравитационного потенциала (Vz – сила тяжести), Z – вертикальная составляющая вектора индукции, I – модуль вектора намагниченности, σ – избыточная плотность, G – гравитационная постоянная, φ – уголмежду осью z и вектором намагниченности, Vzz – вертикальный градиент, Vxz - горизонтальный градиент силы тяжести.
Количественные соотношения между гравитационными и магнитными аномалиями можно рассматривать в двух аспектах: 1)- определение различных характеристик аномалеобразующего тела по изолированным гравитационным и магнитным аномалиям и 2)- разделение полей различных особенностей геологического строения.
Без доказательств запишем некоторые формулы и определения.
Нормой функции f(x) в функциональном пространстве L2 называется величина
(3)
скалярным произведением функций f1 и f2 называется интеграл
(4)
Исходя из этих формул и определений потенциалов и их свойств, можно получить следующие соотношения
(5)
(6)
Используя предыдущие соотношения и уравнения (1) и (2) получаем
(7)
(8)
Из (7) и (8) получаем:
(9)
и
(10)
В вычислительном отношении предпочтение нужно отдать формулам (10) (более точно пересчитывать измеренное поле Vz в поле Vzz, чем измеренное поле Z в U). Чтобы еще более снизить влияние ошибок исходных данных и вычислений на значения I/σ и φ, магнитное поле следует пересчитать на некоторую высоту или осреднить с некоторым радиусом, а вертикальный градиент гравитационного поля вычислять на той же высоте или после осреднения с тем же радиусом.
Можно использовать и следующие формулы
(11)
Разделение полей комплексированием гравитационного и магнитного полей рассмотрим на простейшем примере. Пусть одно аномалеобразующее тело порождает гравитационную и магнитную аномалию («рудное» тело Р), другое – только гравитационную («безрудное» тело S), и пусть известно отношение I/σ и угол φ. Тогда магнитное поле «рудного» тела можно пересчитать в гравитационное. Наиболее просто определить вертикальный () или горизонтальный () градиент силы тяжести. Для этого достаточно выделить из поля Z симметричную и антисимметричную составляющие. Из (2) следует, что
и (12)
. (13)
Далее необходимо трансформировать наблюдаемое поле Vz в поле Vzz или Vxz. Потом вычесть из этих полей соответственно или . Получим или , тем самым, разделив поля.
Очевидно, что таким случаем область применения комплексирования полей для разделения аномалий далеко не исчерпывается.
Лабораторная работа № 1
Определение петрофизических характеристик по данным грави- и магниторазведки
Дано: измеренные значения вертикальной составляющей индукции аномального магнитного поля (Z, в нТл), вертикальный (Vzz) и горизонтальный (Vxz) (в Э) градиенты силы тяжести, полученные пересчетом измеренной аномалии силы тяжести (смотри таблицу вариантов). Измерения проводились вдоль профиля, проходящего вкрест простирания структуры, уподобляемой горизонтальному цилиндру. Расстояние между пикетами 100 м.
Определить отношение I/σ и угол наклонения вектора намагниченности.
Для этого необходимо
1. По выражению (3) определить , , .
2. По определению (4) определить скалярные произведения и .
3. По выражениям (10) и (11) определить искомые величины.