книги из ГПНТБ / Светов Б.С. Теория, методика и интерпретация материалов низкочастотной индуктивной электроразведки
.pdfполя. С учетом этих замечаний может проводиться интерпретация результатов измерения элементов эллипса поляризации магнитного поля.
На основе вышесказанного можно утверждать, что для линейно поляризованных и достаточно малых вторичных полей графики рас пределения и частотные характеристики различных элементов поля весьма сходны и по их виду может быть проведена некоторая клас сификация величин, измеряемых в различных модификациях индук тивной электроразведки (табл. 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