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

книги из ГПНТБ / Светов Б.С. Теория, методика и интерпретация материалов низкочастотной индуктивной электроразведки

.pdf
Скачиваний:
3
Добавлен:
24.10.2023
Размер:
11.24 Mб
Скачать

Вторичное поле магнитного типа в данном случае очень мало и представляет собой поле бесконечного соленоида. Вторичное поле электрического типа, как и в случае шара, слабо зависит от ча­ стоты, а при 0i^>ao становится независимым и от электропровод­ ности цилиндра. При со' = 0 вторичное поле электрического типа не исчезает и совпадает с соответствующим возмущением постоянного

поля, когда вектор Е перпендикулярен оси цилиндра.

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

проанализированного

случая — ^-поляризованной

плоской

волны.

Соответствующие выражения могут быть получены из

формулы

(IV.45) путем замены

<р—>-ф — ф0 и подстановки

вместо

Е°

и Я 0

значений первичного поля линейного тока на оси

цилиндра.

 

СФЕРОИД

Определение вторичного поля

Известны решения задач об идеально проводящем сфероиде, находящемся в поле плоской волны или коаксиального электриче­ ского либо магнитного диполя [28]. При некоторых упрощающих предположениях можно найти решения тех же задач и для неиде­ ально проводящего сфероидального проводника. Они обычно нахо­ дятся через потенциалы Абрагама электрического и магнитного ти­ пов [36], которые удовлетворяют скалярным уравнениям Гельм­ гольца в сфероидальных координатах. Потенциалу Абрагама электрического (магнитного) типа соответствует одна компонента магнитного (электрического) поля, направленная по координате ц> относительно оси вращения сфероида. Полное поле определяется, как и в случае шара, суммой полей электрического и магнитного типов.

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

Для полей электрического типа примем за

основу

следующую

постановку

задачи. Пусть

сфероид с

конечной

проводимостью <Ti

находится

в однородной

вмещающей

среде с

проводимостью со

и пусть он возбуждается

коаксиальным электрическим

диполем.

Устремим расстояние от точки наблюдения до диполя к бесконечно­ сти, а частоту поля к нулю с тем, чтобы перейти к первично-одно­

родному электрическому полю и заменить уравнение

Гельмгольца

более простым

уравнением Лапласа. Вместе с тем избавимся

от рассмотрения

сфероидальных функций и сможем

выразить ре-

110

Рис. 26. Сфероид в переменном электромаг нитном поле

шение через присоединенные функ­ ции Лежандра от сфероидальных координат.

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

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

Чтобы избавиться от рассмотрения уравнения Гельмгольца и внутри проводника, устремим его проводимость к бесконечности и заменим условия сопряжения (1.8) краевыми условиями равенства

нулю

тангенциальной составляющей

электрического

поля

(1.10).

В обеих задачах будем рассматривать вытянутый сфероид.

 

 

 

Введем вытянутую сфероидальную систему координат (рис. 26),

в которой координаты (£, т), <р) связаны с декартовыми

прямоуголь­

ными координатами (х, у,

г)

соотношениями

 

 

 

 

 

Х=сУ

(S2— 1)(1 — tf) COS <р,

}> =

cV(¥—

1)(1 — fj2) Sin ср,

Z={%q

 

 

 

 

 

 

 

 

 

 

(IV.49)

и изменяются в областях

1 ^ g ^ o o ,

— l

^ r j ^ l

, 0 ^ < р ^ 2 я .

Коорди­

натными поверхностями

являются:

g = const — семейство

вытяну­

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

2с;

TJ =

= const — семейство софокусных двуполостных

гиперболоидов

вра­

щения с тем же междуфокусным

расстоянием;

ц> = const — семей­

ство

полуплоскостей, проходящих

через ось 0z под углом <р к

пло­

скости Oxz. В этих координатах

сформулированные

выше

задачи

можно решать непосредственно через компоненты магнитного поля

Я ф (для полей электрического типа) и электрического

поля Ev (для

полей магнитного типа).

 

 

При отсутствии зависимости поля

от координаты

ср эти компо­

ненты удовлетворяют

уравнению

 

 

Уё^-§т[(У¥=Т)(Н,,

Ef)}+VT^~rX

X

[(VT=W)(H9,

£ , ) ] = ( )

(IV.50)

111

и граничным условиям на поверхности сфероида при £ = £oi

Н ъ = н Ч ~£ ' - Ж { У ¥ = ^ И ^ ^ --Ж^Г^НЛ

(IV.51)

 

i - ^

W ^ Z

j

— ^

i l

^ E

j .

(IV.51'

В случае

идеально

проводящего

 

сфероида

граничные условия

(IV.51') заменяются на одно:

 

 

 

 

 

 

 

 

 

£

= 0 .

 

 

 

 

(IV.51")

Частные

решения

уравнения

(IV.

50)

выражаются

через

присо­

единенные функции Лежандра первого и второго рода с азимуталь­ ным числом т = 1:

Pn{'l)\P\® . (IV.52)

Первое решение конечно в начале координат, а второе — регулярно в бесконечности.

Первичное поле возбудителя запишем, как обычно, на основе представления функции Грина уравнения Лапласа для неограни­

ченного пространства

в сфероидальной системе координат [44]:

со

п

 

л = 0

т — 0

 

X

Рп Wo Рп (ri) Q? Go) Рп (6),

(IV.53)

где (go, т)о, <ро) —координаты точки расположения возбудителя

поля

i<| |o» a R— расстояние между диполем и точкой наблюдения. По­

лагая T](j=l (диполь

на

оси сфероида), т = 0 (отсутствует

зависи­

мость поля от координаты ср)

и g 0 - > ° °

(диполь удален

в бесконеч­

ность), с учетом соотношений

Я = гохЛ, £ = го1:Л* найдем

выраже­

ния для первичного

поля: соосного

электрического диполя

(при

H ° = l B

u

Р ' ( ч )

P l ^ ' « ^

" Г Р'

W Р ! «>

 

(IV.S4)

и соосного магнитного диполя

(при оо-»-0

А*а=-^—^-\

 

 

 

 

\

 

 

 

4nR

I

 

 

El

 

Р\ (ч) Р\ а ) =

to^-f

 

Я{ (v)) Я!

ft).

(IV.55)

Здесь учтено, что с|о = го — расстояние

от

точки расположения

ди­

поля до центра сфероида, и обозначено

Е°

и Н° — напряженности

первичного электрического и магнитного полей в начале

координат.

112

Будем искать суммарное поле в случае возбуждающего элек­ трического диполя в среде с индексом 0 —

Я 9 = а 0 £ о

р\ { r i ) [р[ (g) + a (а) Q\ (?)]

(IV.56)

и в среде с индексом 1 —

Нщ = ЪЕ° \ Р\ (т)) [р (а) Я} (£)] •

(IV.56')

Подставляя эти выражения в граничные условия (IV . 51), полу­ чим

а ( а ) =

 

 

?J—

_ - .

(IV.57)

I

° J 2

6 о , - 1

g„i ^

«1 g g , - l j

 

Поле магнитного диполя, возмущенного идеально проводящим сфероидом, находящимся в непроводящей среде, записывается в виде

E^i^H^PlwlPlft+airiQl®],

 

(IV.58)

где, На основе (IV.51")

 

 

 

2

£ 0 , - 1

£ 2 , - 1

 

Обсуждение полученных результатов

 

 

В экваториальной плоскости сфероида, где

г ) = 0 , г = с| — рас­

стояние точки наблюдения от центра

сфероида, c'=c|oi радиус

сечения сфероида, 6 = с У | 2

э 1 1 его

большая

полуось, вторичное

магнитное поле электрического типа может быть записано в сле­ дующем обобщенном виде [52]:

где

 

 

Hv=E°auaSR,

 

 

 

 

 

(IV.60)

 

 

 

1

 

fl2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

с3

 

\^

3

а2

+ с2

 

 

 

 

 

аа2

+ с2 /

0 \

/ g j - ^ g

+ c

2

с

Л

J?o.a2

+

c2\ '

 

I

'1 /

/ а 2 + с2

с

 

Уа? + с2

\

а 1

а 2

/

 

 

 

 

 

 

 

 

 

 

(IV.61)

Множитель S отражает связь величины вторичного магнитного поля со свойствами проводника — степенью его вытянутое™

8 Заказ № 271

ИЗ

и относительной электропроводностью. Множитель R характеризует

распределение поля в пространстве и зависит

только

от

степени

вытдаутости

проводника.

 

 

 

 

 

 

 

 

 

При 6/а->-1 (сфероид вырождается в шар)

получим

выражение

вторичного поля для сферического проводника —

 

 

 

 

 

 

1 — °0

R-

d2_

 

 

 

 

 

(IV.61')

 

1 + 2

/-2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н9=Е°ап

 

r2

 

 

 

 

 

(IV.60')

 

1 + 2 -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при 6/а—>-оо

(сфероид вырождается

 

в цилиндр)—соответствую­

щие выражения для цилиндрического проводника •

 

 

 

 

 

 

 

 

 

1—*!-

 

 

R-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(IV.61")

 

 

 

HID

= £ 0 ( ° i - * o ) -Д22 r

,

 

 

 

 

 

 

 

 

 

(IV.60")

 

 

 

совпадающие

с

полем

элек­

 

 

 

трического типа

в

формуле

 

 

 

(IV.58)

при

\i\ = \io и

 

 

 

 

Множитель 5 может быть

 

 

 

приближенно

записан

в сле­

 

 

 

дующем

обобщенном

виде:

 

 

 

 

 

1

 

 

 

(IV.62)

 

 

 

 

 

 

 

 

 

 

 

 

где коэффициент

т изменя­

 

 

 

ется в зависимости от вытя-

 

 

 

нутости

тела

в

незначитель­

 

 

 

ных пределах — от

1,8

до 2

 

 

 

и может

полагаться

равным

 

 

 

2, а

значение

/

однозначно

 

J вЬ/а определяется

 

 

отношением

 

 

 

полуосей

сфероида

и

нахо-

Рис 27. График функции / для сферой-

 

Дится

с

помощью

графика

дальнего проводника

 

(рис.

27).

 

При

b/a ~> 1

114

а при 6 / а - > оо / - ^ 0 . Для проводников, у которых / > 2

0 i

(сферического типа), характерно незначительное увеличение вто­ ричного поля при повышении электропроводности и сильная его зависимость от степени вытянутости проводника. Проводники, у ко­

со

торых / < 2

(цилиндрического типа),

отличаются незначитель-

0 i

 

 

 

их вытянутости

ным возрастанием вторичного поля при увеличении

и сильной зависимостью

поля от удельной электропроводности.

Рассмотрим поля магнитного типа.

 

 

В плоскости экваториального сечения Я| = 0, а

 

 

 

1

с У~а2 +~с2

*

H

r i = = H ° ы

Уа2 + с2

 

 

У а2 + с2 — с

а 2

 

х ( v2'

- i n £ S 4 + C V

<iv-63>

\ У г2 -\- с2

у г2 + с2 — с J

 

Если, как и в предыдущем рассмотрении, зафиксировать

а и устре­

мить с-*-0, то

 

 

Я „ = / / ° - ^ - .

(IV.63')

При этом вторичное поле совпадает с полем сферического про­

водника.

 

 

При с-У- оо вторичное поле

 

 

Я , = Я о ^ 1 п ^

(IV.63")

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

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

Вторичные поля электрического и магнитного типов

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

8*

115

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

Взяв это строгое математическое понятие за основу, будем тол­ ковать его более широко с физической точки зрения. Назовем вто­ ричным полем магнитного типа такое поле, которое связано с то­ ками, замыкающимися внутри локального проводящего объекта, а вторичным полем электрического типа — такое поле, которое свя­ зано с токами, проходящими через локальный проводник, но замы­ кающимися вне негоi . В области достаточно низких частот, где воз­

буждающее их электромагнитное поле может

рассматриваться как

невзаимодействующая сумма

электрического

и магнитного

полей,

первое из них связано с индуцированием

вихревых токов

внутри

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

интер­

претироваться как результат

воздействия

на проводник

первич­

ного электрического поля, которое в первом приближении

может

считаться постоянным.

 

 

 

 

Вторичные поля магнитного типа для достаточно широкого класса локальных проводящих объектов могут быть записаны в обобщенном виде:

 

hl=-m==*(p2)~R0)>

(IV-64)

V

) — комплексная функция

параметра

проводника, опреде­

где Ф ( р 2

ляющая

частотную характеристику

Я 1 , a R

(г)—действительная

функция координат возбудителя поля и точки наблюдения, опреде­ ляющая распределение Я 1 в пространстве. Из такой записи следует, что с точностью до произвольного множителя частотная характери­ стика вторичного поля слабо зависит от расположения возбудителя и приемника, а распределение поля на съемочном планшете слабо изменяется с частотой. Частотные характеристики вторичных полей магнитного типа в первом приближении не зависят от способа за­ дания и вида измеряемой компоненты поля, проводимости вмещаю­ щей среды, а в случае немагнитных тел—-также от формы и рас­ положения проводника. Это позволяет по частотной характеристике сферического проводника, возбуждаемого однородным магнитным полем (см. рис. 22), определить параметр аномального объекта р2 , записываемый в виде р а = тсо, где т.— некоторая функция электро­ магнитных свойств и размеров проводящего объекта. Она изме­ ряется в секундах и может быть названа постоянной времени про­ водника. Конкретная форма записи т для тел различной формы находится из решений электродинамических задач или на основе физического моделирования.

1 В геофизической литературе вместо термина вторичное поле «магнитного типа» употребляются идентичные по содержанию термины «вихревые» или «ин­ дукционные» поля, а вместо термина вторичное поле «электрического типа» — «дивергентные» (или «кондуктивные») поля.

116

Связь частотных характеристик вторичных полей магнитного типа с постоянной времени проводящих объектов лежит в основе применения методов низкочастотной индуктивной электроразведки при поисках месторождений проводящих руд. Вид частотной зави­ симости этих полей от параметра р2 =та> (см. рис. 22) позволяет рациональным выбором частоты улучшить соотношение между ано­ малиями от объектов с различной электропроводностью и повы­ сить глубинность поиска хорошо проводящих рудных тел. Знаниепостоянной времени аномальных объектов дает возможность в не­ которых случаях определить их геологическую природу.

Частотные характеристики вторичных полей магнитного типа су­ щественно изменяются в зависимости от магнитной проницаемости аномальных объектов (см. рис. 21), но в отличие от магнитораз­ ведки на наблюдаемые аномалии не влияет вектор остаточной на­

магниченности. Это

составляет основу

применения низкочастотной

электроразведки при поисках и разведке месторождений

магнитных

руд в условиях сложного естественного магнитного поля.

Вторичные поля

электрического

типа могут быть

записаны

в следующем обобщенном виде:

 

 

где функции 5 и R для достаточно общего случая сфероидального проводника определены (IV . 61), а величина а представляет собой эквивалентный радиус сечения тела, перпендикулярного напряжен­ ности электрического поля. Величина аномалий электрического типа существенно зависит от импеданса первичного поля, проводи­ мости вмещающей среды и вытянутости проводника в направлении возбуждающего электрического поля. Частотные характеристики этих полей определяются частотной характеристикой импеданса

поля и функцией Ф, на которую также влияет вытянутость провод­

ника и его относительная электропроводность. Для

изометрических

проводников она близка 1, а для вытянутых совпадает с функцией V

кругового цилиндра (IV . 46) .

 

Сложная

зависимость

вторичных полей электрического типа

от большого

количества

факторов, не связанных

непосредственно'

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

Всилу большого различия свойств вторичных полей магнитного

иэлектрического типов методику низкочастотной электроразведки надо строить так, чтобы аномальные объекты обусловливались пре­ имущественно полями одного типа. Применительно к определенным

117

геологическим задачам этого можно в некоторой мере добиться, выбирая вид возбудителя поля и его ориентировку относительно простирания искомых объектов. Свойства вторичных полей элек­ трического и магнитного типов описаны в ряде работ [30, 52, 58].

ПРИБЛИЖЕННЫЕ РАСЧЕТЫ ВТОРИЧНЫХ ПОЛЕЙ МАГНИТНОГО ТИПА В ОБЛАСТИ МАЛЫХ ПАРАМЕТРОВ

Определение вторичного поля

Как уже было показано в гл.

I , для вторичных полей

магнитного

типа

от

локальных

немагнитных

проводников всегда

существует

первое

приближение к решению,

пропорциональное (k\—k2Q),

ко­

торое

получается

путем подстановки в интегральное

уравнение

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

Как уже указывалось выше, такое разделение возможно лишь в некоторых системах координат и при определенном виде первич­ ного поля, поэтому несколько упростим постановку задачи, что по­ зволит получить ее приближенное решение для значительно более

широкого класса граничных поверхностей. При

этом будем опи­

раться на более общее физическое определение

полей

магнитного

типа, данное в предыдущем

разделе.

 

 

 

Рассмотрим одиночный проводник с проводимостью

01,

поме­

щенный в непроводящую

вмещающую среду

(ао = 0),

считая

I kit \ <Cl, где / характеризует размеры локального проводника. При

= О на поверхности проводника

нормальная компонента плотно­

сти тока должна равняться нулю

(/п = 0) и вторичное поле, возни­

кающее за счет токов, индуцированных внутри

проводящего объ­

екта, будет относиться к магнитному типу. В этих

предположениях

сформулированная задача может

быть сведена

к

последователь­

ному решению двух более простых задач.

Первая из них заключается в определении токов, индуцирован­ ных в проводящем теле первичным полем #°, которые должны удо­

влетворять внутри проводника уравнению

(IV.66)

rotj=k2H°,

а на его поверхности — условию / п = 0. Поскольку d i v / = 0,

сформу­

лированная задача имеет единственное решение [35].

 

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

ков и при условии / и = 0 решается достаточно просто путем

вычис­

ления интегралов вида

 

118

по объему V проводящего тела в трехмерном случае, или

J l ( Р ) = = ^ Г } j ( Q ) l n -W7Q) dS<*

(IV'67'>

по сечению цилиндрического тела 5 в двумерном случае.

Плоская задача

Будем считать, что цилиндрический проводник с сечением об­ щего вида вытянут в направлении оси Oz и находится в однородном магнитном поле, направленном по оси Оу.

Тогда вектор плотности тока имеет только одну компоненту j Z r удовлетворяющую уравнениям

 

- § г = 0 '

т = - ^ я ° -

< I V - 6 8 >

Если совместить начало координат с центром тяжести сечения

тела, то

 

 

 

 

 

jz=-k\H\.

 

(IV.69)

Подставляя формулу (IV.69)

в ( I V . 6 7 ' ) ,

получим следующие

выражения для вектор-потенциала

и компонент вторичного магнит­

ного поля:

 

 

 

 

К (х,

k\H°

,

1

dxx dyx,

у ) = - - ^ г

J

х, I n ( Л Q )

s

где p (P, Q) =

i(x-Xi)2+(y-yif.

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

сводится к квадратурам

(IV . 70) . Интегралы вида (IV.70) для тел

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

ментарных или табулированных функциях. В более общем

случае

произвольного сечения цилиндрического

проводника

эти

инте­

гралы могут быть найдены

приближенными

методами,

например,

графически. Необходимые

сведения по этому

поводу

содержатся,

в работах О. М. Косенкова

и Б. С. Светова

[34, 83]. В этих работах

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

119/

Соседние файлы в папке книги из ГПНТБ