Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MISCELLANEOUS / Geophysics / Geophysics Methods Костицын В. И

..pdf
Скачиваний:
6
Добавлен:
03.05.2023
Размер:
6.07 Mб
Скачать

Магниторазведка

ме того, при расчетах можно считать, что

0 , где

0

4 10 7

– магнитная проницаемость воздуха.

 

 

 

 

2.12. Прямая и обратная задачи для намагниченного вертикального бесконечно длинного столба (стержня)

1. Прямая задача. Пусть на глубине h залегает вершина бесконечно длинного столба (вертикального цилиндра или стержня) сечением s (рис. 2.4). Его можно представить как тело

одного полюса (m) с интенсивностью намагничения

(J ) , на-

правленного вдоль оси z , и «магнитной массой» m

Js . Ниж-

ний полюс столба расположен очень глубоко, поэтому его влиянием можно пренебречь и считать, что вся «масса» сосредоточена на вершине столба.

Необходимо найти напряженность поля вдоль профиля x над телом. Потенциал от верхнего полюса столба в точке P будет равен потенциалу точечной массы:

U

m

 

 

m

 

 

.

(2.14)

 

 

 

 

 

 

r

 

 

 

 

 

 

 

x2

h2

 

 

 

 

 

 

 

 

 

 

 

Составляющие поля выражаются производными потенциала по соответствующим осям координат:

Z a

dU

 

 

J sh

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dh

 

(x

2

 

h

2

)

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H a

 

dU

 

 

 

J sx

 

 

 

 

 

 

 

,

(2.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

(x

2

h

2

)

3

2

 

 

 

 

 

5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J s

 

 

T

Z 2

H 2

 

.

 

 

 

 

(x 2 h 2 ) 5

 

a

a

a

 

2

 

Используя полученные формулы, можно построить графики напряженности поля (рис. 2.4). Легко видеть, что над столбом будут максимумы Ta и Z a , а значения их одного знака,

119

В.К. Хмелевской, В.И. Костицын

положительные при вертикальной J a . Горизонтальная состав-

ляющая (H) слева будет иметь максимум, а справа – минимум. Вдали от столба аномалии исчезают. В плане над таким столбом изолинии Ta и Z a будут иметь вид концентрических окружно-

стей одного знака.

Рис. 2.4. Магнитное поле вертикального бесконечно длинного столба

2. Обратная задача. Решение уравнений (2.15) дает возможность по характерным точкам на графиках определить глубину залегания верхней кромки вертикального бесконечно длинного столба (h). Так, центр столба находится в точке, где x=0, а

Z

m ax

T

m ax

Js / mh 2 .

 

 

 

Для точек, удаленных на расстояния x Z 1 2 от начала координат, в которых Z равно половине максимального,

Z

1/ 2

Js h / 2h 2

Js h / (x2

1/ 2

h 2 ) 3 / 2 .

 

 

Z

 

Решив это уравнение, получим

x Z 1 2 0,7 h . Аналогич-

ным образом находятся связи и между другими характерными

120

Магниторазведка

точками x T 1/ 2 , xH Э (экстремумы на составляющей H), xZH

(абсциссы точек пересечения Z и H). В результате для расчета h по абсолютным значениям этих параметров получаются следующие формулы:

h 1,4

xH Э

1,3

xZ 1/ 2

 

xT 1/ 2

 

xZH

.

(2.16)

 

 

 

 

 

 

 

 

 

 

Зная h, можно оценить величину магнитной массы:

m Js Z

m ax

h2 T

m ax

h2

3,67 H

h2 .

(2.17)

 

 

 

m ax

Так как J Tср æ, где Tср – среднее значение полного вектора напряженности поля в изучаемом районе, а æ – магнитная восприимчивость столба, то

æs = m / Tср .

Отсюда, если известна æ по измерениям на образцах, то можно определить площадь поперечного сечения столба (s).

2.13. Прямая и обратная задачи для вертикально намагниченного шара

1. Прямая задача. Пусть вертикально намагниченный шар с центром на глубине H залегает под началом координат (рис. 2.5). Необходимо определить напряженность поля вдоль профиля x. Потенциал шара можно представить как потенциал диполя, помещенного в его центре. Поэтому, согласно (2.14), потенциал шара с магнитным моментом M = JV (или магнитной массой m = M) равен

U

M cos

 

MH

 

JVH

(2.1

 

 

 

 

 

.

r 2

 

r 3

 

(x2 H 2 )3 2

8)

121

В.К. Хмелевской, В.И. Костицын

Рис. 2.5. Магнитное поле шара

Отсюда, взяв производные, найдем элементы магнитного поля шара

Za

dU

 

JV (2H 2 x2 )

,

dH

 

(x2 H 2 )5 2

H a

dU

 

3 J VHx

(2.1

 

 

 

,

dx

 

(x2 H 2 )5 2

9)

 

 

 

 

 

 

 

 

 

 

 

JV 4H 2

x2

.

 

 

 

 

 

 

 

 

T

 

Z 2

H 2

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

(x2 H 2 )5 2

 

 

 

 

 

 

 

 

 

Анализ этих формул и построенных по ним графиков по-

казывает,

 

 

что

 

над

центром

шара

(x

0)

будут

Z

max

T

 

 

2 JV / mH 3 , а H

0 . При x

 

аномалии ис-

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чезают. При x

 

2 H получим

Za

0 , при x

 

2 H

Za

0,

а при x

 

 

 

 

 

 

 

 

 

 

 

 

2 H

Za 0 .

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в плане над шаром изолинии Za и

Ta

бу-

дут иметь вид концентрических окружностей. При этом изолинии Za будут двух знаков, а Ta – одного.

122

Магниторазведка

2. Обратная задача. Решение уравнений (2.19) теми же приемами, что и для столба, дает возможность по характерным точкам на графиках найти глубину центра вертикально намагниченного шара:

H a

1,8

xZ1 2

1,8

xZH

 

1,5

xT1 2

0,7

xZ 0

0,5

xZ min

 

,

(2.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где xZ1 2

и xT1 2 – абсциссы точек половины Za , Ta ,

xZ 0 – точки

с Za 0;

xZ min – точки с Za

Zm in .

 

 

 

Зная H, можно оценить магнитную массу шара (m):

m JV Z

max

h 3 2 T

h 3 2 .

 

max

 

Отсюда, так как J Tср æ, то получим æV m / Tср. Если известны Tср и æ, то можно определить объем шара.

2.14. Прямая и обратная задачи для вертикально намагниченного тонкого пласта бесконечного простирания и глубины

Пусть на глубине h параллельно оси y расположен беско-

нечно длинный вертикальный пласт (с толщиной l, меньшей глубины залегания), намагниченный вертикально (рис. 2.6). Определим для простоты лишь Z a вдоль оси x.

Поскольку нижняя часть пласта расположена глубоко, то влияние магнитного полюса глубоких частей пласта будет мало, и можно считать, что магнитные массы сосредоточены вдоль поверхности в виде линейных полюсов. Магнитная масса единицы длины пласта равна dm/ dy Jl .

Разобьем пласт на множество тонких «столбов». Тогда притяжение пласта будет складываться из притяжения всех элементарных столбов, а вертикальная составляющая его магнитного притяжения будет равна интегралу в пределах от до

123

В.К. Хмелевской, В.И. Костицын

(по оси y) выражения для притяжения элементарного столба. Потенциал элементарного тонкого столба равен

dU dm / R J l dу / x2 y2 h 2 ,

а вертикальная составляющая –

dZ

a

d (dU ) / dh J l h d у / x2

y 2 h2 3 / 2 ,

 

 

 

откуда Z равно

Z a

Jl hdy

 

2J lh

(2

 

 

 

.

 

(x2 y 2 h2 )3 2

(x2 h2 )

.21)

Рис. 2.6. Магнитное поле тонкого пласта бесконечного простирания

График Za будет иметь максимум над центром пласта и асимптотически стремиться к нулю при удалении от пласта. В плане над пластом будут вытянутые аномалии Za одного знака. Анализируя формулу (2.21), можно найти связи между глубиной залегания пласта (h) и x1/2, т.е. абсциссой графика, где

Za

Z m ax / 2 ; h x 1 2 .

 

 

 

 

Магнитная

масса

единицы

длины

равна

m

Jl Z m ax h / 2 .Заменив J ≈ æTср , получим læ = m/Tср. Зная

Tср и æ, можно рассчитать ширину пласта.

 

 

124

 

 

 

 

 

Магниторазведка

2.15. Прямая и обратная задачи для вертикально намагниченного горизонтального цилиндра бесконечного простирания

Пусть на глубине H параллельно оси y расположен бес-

конечно длинный цилиндр с магнитным моментом единицы длины, равным M Js , где J – интенсивность намагничива-

ния, постоянная для всего цилиндра и направленная вертикально, s – поперечное сечение цилиндра (рис. 2.7). Требуется определить напряженность поля вдоль оси x . Поле такого цилиндра можно считать эквивалентным полю бесконечного числа вертикальных магнитных диполей, центры которых расположены по оси цилиндра.

Рис. 2.7. Магнитное поле горизонтального цилиндра бесконечного простирания

Потенциал в точке P от элементарного диполя определяется согласно уравнению (2.8):

dU

dM cos

 

JH dv

 

JH ds dy

,

R2

 

R2

 

R2

 

 

 

 

125

В.К. Хмелевской, В.И. Костицын

где R

x2 y2 H 2 .

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

U

JH

 

ds

dy

 

 

 

.

 

S

(x2 y 2 H 2 )3 / 2

Tак как

ds s

 

, а

 

 

 

 

 

 

dy

 

 

 

 

 

2

,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2

 

y 2

H 2 )3 2

 

 

 

 

x2

H 2

 

U

 

2 Js H

 

 

 

 

 

2 M H

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2

H 2 )

 

 

 

 

(x2 H

2 )

 

 

 

 

 

 

 

 

(2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Za

 

dU

 

2M

 

 

H 2

x2

 

 

.22)

 

 

 

 

 

 

 

 

 

.

 

 

 

 

dH

 

 

 

 

 

(x2

H 2 )

 

 

 

Можно видеть, что при x 0 будет максимум Za, а при x H Za 0 . При x0 значения Za будут отрицательны, при xH – положительны. В плане над горизонтальным

цилиндром будут вытянутые аномалии двух знаков.

При решении обратной задачи глубину залегания цилинд-

ра можно определить по формуле: H

x0

0,7

xm in

,

где

x0 и xmin – абсциссы точек, в которых Za

0 и Za

 

 

Zm in . Зная

H, можно найти погонную массу цилиндра

M

Z

max

H 2

/ 2 .

 

 

 

 

 

 

 

 

Заменив J ≈ æTср, получим æs = M/Tср. Зная Tср и æ, можно рассчитывать площадь сечения цилиндра.

126

Магниторазведка

2.16. Численные методы решения прямых и обратных задач магниторазведки

Пользуясь формулой (2.13), можно решать прямые задачи для тел других форм и невертикальной намагниченности. Практически эти расчеты реализуются с помощью ПЭВМ. Обратные задачи можно решать методом сравнения наблюденных графиков или карт аномальных магнитных полей с теоретически рассчитанными для меняющихся геометрических параметров и магнитных свойств. Получив наименьшие расхождения между ними, теоретические параметры совпавшей модели переносят на реальные объекты. Они играют роль одного из эквивалентных решений.

2.17. Измеряемые параметры геомагнитного поля

Измерения магнитного поля Земли и его вариаций проводят на стационарных пунктах – в магнитных oбсерваториях, которых насчитывается на Земле около 150. Магнитное поле измеряется и при магниторазведочных работах. Абсолютные определения полного вектора напряженности геомагнитного поля сводятся к измерению, как правило, трех его элементов (например Z, D, H). Для этого применяют сложные трехкомпонентные магнитные приборы, которые называются магнитными теодолитами и вариационными станциями.

При геологической разведке измеряют абсолютные Z, T и относительные (по отношению к какой-нибудь исходной, опорной точке Z, T) элементы.

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

2.18. Оптико-механические магнитометры

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

127

В.К. Хмелевской, В.И. Костицын

компасе) для измерений приращений горизонтальной составляющей в двух точках (ΔH), либо вокруг горизонтальной оси для измерения приращений вертикальной составляющей (ΔZ). Углы отклонения , пропорциональные H или Z, определяются с помощью специальной оптической системы. Сняв отсчеты по магнитометру в двух точках (n) и (n0), можно определить приращение

Z c (n n0 ) ,

где с – цена деления магнитометра, которую определяют путем градуировки с помощью эталонировочных магнитов. На этом принципе был построен магнитометр, названный весами Шмидта, применявшийся в магниторазведке для измерения Z свыше 50 лет. Среди отечественных магнитометров к этому типу относились полевые приборы М-2, М-18, М-27, а также приборы для измерения магнитных cвойств образцов М-14 и астатические магнитометры. Погрешности в определениях Z с помощью таких магнитометров составляют ± 2–5 нТл.

2.19. Феррозондовые магнитометры

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

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

128