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

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

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

поля. С учетом этих замечаний может проводиться интерпретация результатов измерения элементов эллипса поляризации магнитного поля.

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

 

 

 

 

 

 

 

Т а б л и ц а 1

 

 

 

 

 

 

Частотные

характеристики

 

Графики

распределения

Подобные

Im Н1

 

 

 

 

 

Подобные Re Н'

Подобные

Н

хх

\тН

, <р ,

<?

нх-т,

на-но,

 

 

х

 

х

а

 

 

 

 

 

 

 

 

 

Re Нх

Подобные

Н1у

Im Ну,

 

Hj,

Re Ну,

1а

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

Г Л А В А III

ПЕРЕМЕННЫЕ ЭЛЕКТРОМАГНИТНЫЕ ПОЛЯ В ПРИСУТСТВИИ ГОРИЗОНТАЛЬНОСЛОИСТОЙ СРЕДЫ

ПРОСТЕЙШИЕ ВИДЫ ПОЛЕЙ В ОДНОРОДНОЙ СРЕДЕ

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

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

В соответствии с уравнением Максвелла (1.2) условием суще­ ствования поля, близкого к однородному магнитному, является весьма низкая проводимость среды (о - > - 0) .

 

Если # = const, то на основании формулы (1.1) rot Е = const = С,

а

циркуляция электрического поля по любому замкнутому контуру

L

равна:

где 5 — площадь, ограниченная контуром L . Задаваясь различной

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

L

31

откуда

т. е. электрическое поле линейно растет с удалением от оси симмет­ рии. Такое поле аппроксимирует поле круглой петли с током в ее средней части (в окрестности оси симметрии).

Условием существования поля, близкого к однородному элек­ трическому, является в соответствии с уравнением (1.1) его доста­ точно низкая частота. Закон изменения магнитного поля в этом слу­ чае также в значительной мере произволен.

Плоские волны

Наиболее общим решением уравнения Гельмгольца (1.14) или (1.17) в однородном пространстве в декартовой прямоугольной си­ стеме координат (х, у, z) является линейная комбинация функций вида

~

^

i Ckxx + lkvy +

lkzz)

 

F =

F°e

у

,

(III.

V

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

V V V

торных потенциалов. Величины kx, ky, kz должны удовлетворять

соотношению k2 + k2 + k2z = k2,

где & —волновое число среды. Вы-

брав ориентацию

осей

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

образом,

чтобы kx = 0,

можно записать выражение ( I I I . 1 ) в виде

 

 

 

 

 

 

 

(III.2)

В общем случае, поскольку ky

= ksinQ и

& z = & c o s 0 ,

удовлетворяя

V

V

V

 

 

V

соотношению k2y + k2z = k2, могут иметь различные аргументы, угол 0

следует считать комплексным. Мгновенные значения векторов поля, соответствующие выражению (II1-2), записываются следующим об­ разом:

v

 

v

+ik". Это выражение определяет волно-

где ky = k' +ik"

 

и kz = k'

У

У

Z

z

вое поле, в котором поверхности равных фаз и амплитуд представ­

ляют собой плоскости k'yy + k'zt=Ci

и kyy + k"z = C2

(в общем

слу­

чае непараллельные). Такое поле называют плоской

(в общем

слу­

чае неоднородной) волной [6]. Если

угол 6—действительный

(ар-

32

гументы ky и kz одинаковы), то формула (111.3) характеризует пло­ скую (однородную) волну, распространяющуюся под углом б к оси Ог.

Если под

V

 

 

 

 

то

в

соответствии

F понимать компоненту поля Ех,

с уравнением

(1.1)

напряженность

магнитного

поля

будет

иметь

две компоненты по осям у

кг:

 

 

 

 

 

 

 

В= (ёу -

i - ег

EV

< V +

V ) = X

( i

cos

0 -

i , s i n 9 )

X

\ У tofj.

z

COM. /

 

COJA

4 У

 

 

г

/ '

ч

 

 

x £ o e * G - n a + . c o . T 0

 

 

 

 

 

(III.4)

где e y и ez — орты соответствующих координатных осей. В неодно-

родных плоских волнах {ky и kz имеют разные аргументы) компо­ ненты Ну и Hz будут сдвинуты по фазе и магнитное поле эллипти­ чески поляризовано. В силу этого оно всегда имеет продольную компоненту, совпадающую с направлением распространения волны (перпендикулярную к поверхности равных фаз). Такая волна назы­ вается магнитной или поперечной электрической [42]. Если под F понимать компоненту магнитного поля Нх, то волна будет электри­ ческой, или поперечной магнитной (электрическое поле эллиптиче­ ски поляризовано и имеет продольную компоненту). Всякое плоское поле можно представить как сумму поперечных электрических и магнитных волн. В плоской однородной волне отсутствуют продоль­ ные компоненты поля, плоскости равных фаз и амплитуд совпадают между собой.

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

z = j L = a = j £ . = е -

г

т

( I 1 L 5 )

Это отношение называется также волновым

(или

характеристи­

ческим) сопротивлением среды. Скорость перемещения плоскостей с некоторой заданной фазой у однородных плоских волн (фазовая скорость) определяется действительной частью волнового числа k':

а ослабление амплитуды

(затухание)

—экспоненциаль-

v=—гг,

k

 

 

 

 

 

ным множителем, зависящим от его мнимой части k":

t~h"z.

Мнимая часть волнового числа k"

называется

коэффициентом

поглощения. При изменении z на величину

б =

l/k",

называемую

толщиной скин-слоя, напряженность

поля

однородной плоской

волны убывает в е раз.

 

 

 

 

3 Заказ № 271

33

В квазистационарном случае действительная и мнимая части волнового числа равны и растут пропорционально корню квадрат­ ному из проводимости и частоты:

соответственно

 

 

 

 

 

* = 1 / ^ и

8 = 1 / Х .

 

Г

ар.

 

г <J[X(o

 

 

Поля реальных возбудителей могут

быть

представлены

в виде

суперпозиции бесконечного

числа

неоднородных плоских

волн.

В дальней зоне на участках, размеры

которых

малы по сравнению

с расстоянием до возбудителя, их поля аппроксимируются однород­ ными плоскими волнами.

Цилиндрические волны

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

FZ=F(0)H^(^P),

(III

v

фаз и ампли­

где F° = const. Для этой волны поверхности равных

туд представляют собой коаксиальные круговые цилиндры. В зави-

V

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

сутствует компонента Hz, во второй — компонента

Ег.

В поперечно-магнитном

поле

формула ( I I I . 6 )

выражает элек­

трическое поле линейного тока /

(бесконечно длинного кабеля), со­

впадающего с осью Oz, в однородной среде:

 

E,=

-^-H$\kv).

(III

Напряженность магнитного поля в этом случае имеет вид:

Н9=-!%-

Н[Х) (kP).

(III.8

v

В ближней зоне (р = | kp | <С 1)

Е * = -

[ln р - ° > 1 1 5 9 - 1 - f ] ;

( 1 П - 9 )

Я * = 1 ^ Г [ ! -

^T—IT^P

- 0 , 6 1 5 9 ) ] .

(ШЛО)

34

В дальней зоне (р = | kp | ^> 1)

(111.11)

(III . 12)

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

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

 

1

 

ваются, как

сдвиг фазы между электрическим и

магнитным

 

TP

 

 

л

 

полями стремится к — .

 

Импеданс

цилиндрической волны в ближней зоне

зависит от

расстояния р до источника:

 

 

Z = ^ ^ - / o ) j j L P ( l n / ? - 0 , 1 1 5 9 - / ^ - ) ,

(111.13)

а в дальней зоне совпадает с импедансом плоской волны:

 

 

Z = ^ = - 1 ^ K ^ .

№ 1 4 )

Можно считать, что в этой зоне на небольших участках цилиндри­ ческая волна аппроксимируется плоской волной. В поперечном электрическом поле формула (111.6) при FZ = EZ выражает элек­ тромагнитное поле бесконечного соленоида.

Сферические волны

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

3*

35

В случае магнитного диполя магнитный векторный потенциал

А* выражается формулой

А* = Ж-^-,

(111.15)

где M = IS — момент диполя, эквивалентный моменту витка с пло­ щадью S и током /, а г — расстояние от диполя до точки наблюде­ ния. Векторный потенциал совпадает по направлению с моментом диполя.

Электрический векторный потенциал А имеет в этом случае также одну компоненту, но направленную по координате ср:

A=AV=M-^-(\~ikr)slnB.

(III.16

Компоненты электромагнитного поля и его импеданс равны:

НГ=М-^-(1-ikr)cosB,

 

HB=M-^-(l-ikr-k^slaB,

(111.17)

E,=M-^^-(\-ikr)sinB;

 

 

 

E

 

 

 

1 _

ikr

 

 

 

В ближней зоне (p—\kr\-*-0)

 

с точностью

до

членов 0 (р3 ) бу­

дем иметь:

 

 

 

 

 

 

 

 

 

 

 

^ = = ^ r ( 1

+ ^ ) c

o s 6 >

 

 

 

 

 

Я е = - 4 ^ ( 1 - ' 4 - ) 8 1 п 9 '

 

 

( Ш - 1 9 )

 

 

Мш*

(л

,

, р*

 

 

 

 

 

 

4кг'2

 

 

 

 

 

 

 

 

 

Z = - ^

- =

* ( 0 { i r ( l + / / > 2 ) .

 

 

(111.20)

 

 

 

в

 

 

 

 

 

 

 

В

первом приближении

магнитное

поле

магнитного

диполя

в этой зоне совпадает с магнитным

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

диполя. Однако

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

поле

переменного

магнитного ди­

поля

вследствие

явления

индукции и в ближней

зоне

отлично

от нуля и сдвинуто по фазе относительно магнитного поля

на угол

я_

—. С увеличением расстояния точки наблюдения до диполя г магнитное поле убывает как 1/г3, а электрическое — как 1/г2. Обе

36

компоненты магнитного поля сравнимы по величине, поле имеет су­ щественную продольную компоненту Нг. На одинаковом расстоя­ нии от диполя она вдвое превышает поперечную компоненту.

В дальней зоне

( р =

\kr\-*-oo)

с точностью до

членов о ( — j

имеем:

 

 

 

 

 

 

tfe

= -

M

- ^ s i n e ,

(Ш.21)

 

Е9

 

JZg^-slnB;

 

 

 

 

 

 

(111.22)

В этой зоне продольная компонента магнитного

поля Нг прене­

брежимо мала по

сравнению

с

# е — поле становится поперечным.

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

быть получено

из выражений (III . 15) — (111.22) путем следующих

замен: М-^Р/о,

шц~+о, Н —~Е, А~ А*

ОПРЕДЕЛЕНИЕ ПОЛЕЙ НЕКОТОРЫХ ВОЗБУДИТЕЛЕЙ В ПРИСУТСТВИИ ГОРИЗОНТАЛЬНО-СЛОИСТОЙ СРЕДЫ

Чтобы не усложнять решение поставленной задачи непринципи­ альными деталями, ограничимся изучением пространства, состоя­ щего из трех областей, характеризующихся электромагнитными свойствами | оо, M-ol. I °и M-i и 102, р-г| и разделенных плоско-па­ раллельными поверхностями. В соответствии с формой поверхно­ стей раздела введем декартову ортогональную или цилиндрическую систему координат (s, t, z) и направим ее ось Oz вертикально вверх — перпендикулярно к слоистости среды. Начало координат расположим на поверхности раздела сред с индексами 0 и 1.

Пусть электрический векторный потенциал Л° первичного поля некоторого возбудителя, расположенного в плоскости z = A, в одно­

родной и изотропной среде с волновым числом &о = У шцо00

может

быть описан в этой системе координат только одной

ортогональной

компонентой А0, которая

может быть представлена в виде

 

 

 

Л ° = Г { а ° ) = Г { е - ' л ° | г

- л 1 } ,

 

(111.23)

где

т о = уя2 k20,

а Г{а°)

представляет

собой линейный

опера­

тор

интегрального

преобразования над

функцией

a0 =e _ m »| z _ f t l

37

с действительным параметром интегрирования X, имеющим размер­ ность волнового числа. Ядро оператора Г{} может содержать функции координат s a t параметра X, но не должно зависеть от координаты z и изменяться при переходе через границы сред. Знак

у корня /По = уЯ2 k\ выбирается так, чтобы его действительная

часть была положительна. В частности, при Х = 0 это означает, что mn = —iko. Такой выбор знака у то обеспечивает необходимое убы­ вание первичного поля в бесконечности.

Решение поставленной задачи, выраженное через векторный по­ тенциал А, должно удовлетворять внутри областей с индексами О, 1 и 2 соответствующим уравнениям Гельмгольца, на поверхностях

раздела этих

областей-—граничным условиям

вида (1.18), а в

ок­

рестности

возбудителя — приближаться

к его

первичному

полю

( I I I . 2 3 ) . В соответствии с последним требованием векторный

потен­

циал А должен содержать компоненту

As,

записанную в

форме,

аналогичной

( I I I . 2 3 ) . Поскольку в общем

случае й1\А°фО,

то,

для

того чтобы удовлетворить второму граничному

условию (1.18), не­

обходимо ввести еще одну компоненту векторного потенциала

А,

появляющуюся лишь в случае неоднородной среды — Az.

 

 

Тогда

граничные условия (1.18) могут быть

переписаны

в виде

двух последовательно рассматриваемых систем уравнений: системы

относительно

As

 

 

 

 

 

 

dAs

 

dAs

 

 

^

= ц ж Л , ж ; -3J-=—^±

 

(Ш.24)

и системы, объединяющей As и Az, —

 

 

 

 

AZ=AZ

J _ d i v ^ =

- j

i - d i v J w

(111.25)

при z = 0 и z = -—hi, где hi — мощность первого слоя.

 

Так как в

ортогональной декартовой

и

цилиндрической

систе­

мах координат векторное уравнение Гельмгольца относительно век­ тора А, содержащего две компоненты — As и Аг, разделяется на два

независимых

скалярных

уравнения

относительно этих

компонент,

то с учетом

выражения

(III . 24) задача

определения

As

становится

вполне самостоятельной.

 

 

 

 

 

 

 

Будем искать ее решение в виде:

 

 

 

 

 

 

в среде с индексом О

 

 

 

 

 

 

 

 

Л , о = Г { а , 0

) = Г { а ° + а 1 0 } = Г { е - т ° 1

г - л |

+ 4 2 ) е - т " (

г + Л )

К

(111.26)

в средах с индексами

1 и 2

 

 

 

 

 

 

 

 

Л , 1 1 2 = Г { а , ь 2 } = Г { а 1 ( 2 2 е - т . , 2 г + ^ 2 е т 1 , 2 г } ;

 

(Ш.27)

где m.ii2 =

i X 2 — &2 4

а а<2> и Р( 2 ) — неизвестные

функции

электро­

магнитных

свойств и геометрии среды (верхний

индекс

обозначает

38

количество слоев в структуре; в общих

случаях

он

будет

опус­

каться). Функцию

а<2) уместно назвать

функцией

отражения. По­

скольку выражение

(III . 23) удовлетворяет уравнению

Гельмгольца

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

0, а оператор Г{ } не зависит от свойств

среды,

то и выражения ( I I I . 2 6 ) , (III . 27) удовлетворяют

тому же

уравне­

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

конечности,

необходимо положить

а<^ = 0.

Подставляя

выражения

(111.26), (III . 27) в

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

(III . 24) и символически со­

кращая не зависящий от констант

среды

оператор Г{

} , получим

граничные

условия

относительно

as,

эквивалентные

равенству

(111.24):

при

2

= 0

 

 

 

 

 

 

 

 

м 4

1ма* = — м £ = — д а т

° й ;

 

 

 

да\

да.

да®

 

_

h

 

 

 

~дГ -

-ЯГ = -

- а Г =

-

« о е - " ' * ;

(Ш-24')

при

2

= — h i .

 

 

 

 

 

 

 

 

 

да.

да.

 

 

 

 

 

 

- д Г ~ П ) Г = 0 -

 

 

 

( J " - 2 4 ' )

Эти

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

собой

систему

линейных

алгебраических уравнений относительно искомых функций а<2> и |3(2А Поскольку в нашем рассмотрении основной интерес пред­ ставляет поле в верхнем полупространстве, выпишем выражение для функции отражения а<®, введя удобное для дальнейшего па­ раллельное обозначение а<® (р,):

(2)

— а 0

(2) , v

Яр] +

w] 2 e~

2 , "] f t '

, ш 9 R .

а0

 

(Р-)—

 

2 я ,1 Й 1

,

(Ш.28)

где

 

 

 

 

 

 

 

 

 

Пц=пМ=

m i V 4 ~ m

P l

т

 

(111.29)

В случае однородного полупространства

( « 1 2 = 0)

 

 

 

 

4 1 ) ( ^ ) = " o i .

 

 

 

(Ш.28')

Переходя к определению вертикальной компоненты векторного потенциала Аъ, будем, избегая в нашем рассмотрении ненужной общности, полагать найденную компоненту A s декартовой, т. е.

As — Ах.

39

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