Файл: Даев Д.С. Высокочастотные электромагнитные методы исследования скважин.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.06.2024

Просмотров: 85

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

div Д = 0;

 

 

(2.4)

div H ~ 0,

 

 

(2.5)

где H — напряженность магнитного поля; E — напряженность

электрического поля; у — электропроводность среды;

г — диэлек­

трическая проницаемость' среды;

р — магнитная проницаемость.

Решение системы уравнений

(2.2) —(2.5)

приводит

к волново­

му уравнению

 

 

 

 

Ѵ2Л +

£М =

0,

 

(2.6)

где Л — вектор-потенциал магнитного

типа

(2?= rotiT)

k — волно­

вое число.

 

 

решение

волнового

Для определения поля достаточно найти

уравнения, удовлетворяющее краевым условиям задачи.

 

Не воспроизводя здесь хорошо известного решения, выпишем выражения для компонент поля магнитного диполя в однородной среде. В сферической системе координат с началом, совпадающим с диполем, момент которого ориентирован по оси Ѳ= 0, они имеют

вид [53]

 

 

Дф =

(корМ/4лД2)еад (1 — ikR) sin Ѳ;

(2.7)

HR

= (2M/4nR3) eikR(1 — ikR) cos Ѳ;

(2.8)

Яо = (A4/4nR3) eikR (1 - ikR kW ) sin Ѳ,

(2.9)

где R — расстояние от диполя до точки наблюдения;

 

 

k2 = ш2ер 4- гурсо. .

(2.10)

Поскольку нас интересует влияние среды только на поведение поля и не важны его абсолютные значения, целесообразно упро­ стить выражения (2.7) —(2.9), выразив компоненты £ф. HR и Не в относительных'единицах—долях квазистационарного поля в воздухе, т. е. поля, вычисленного без учета токов смещения.

Компоненты поля в воздухе имеют вид

Дф о =

,(ісорМ/4лД2) sin Ѳ;

(2.11)

HRо-

(2M/4nRs) cos Ѳ;

(2.12)

Hqо — (M/4nR3) sin Ѳ.

(2.13)

Компоненты поля магнитного диполя в относительных едини­ цах запишутся в виде

е9 = еІІ* ( 1 — іЩ ;

(2. 14)

hR = elkR( l - i k R ) ;

(2.15)

he == eikR (1 ikR kW ).

(2.16)

26


Зависимость поля от параметров среды определяется волновым числом k. Представляя k в виде комплексного числа k = a+ib и выделяя действительную и мнимую части, имеем

(2.17)

(2.18)

Как видно из выражений (2.17) — (2.18), влияние среды на по­ ле в сильной мере зависит от отношения сое/у, т. е. от соотноше­ ния между токами смещения и токами проводимости'. При ие/у^> ^>1 (преобладание токов смещения) формулы (2.17) и (2.18)

упрощаются:

 

!

а = со)/ер.

 

(2.19)

Полагая (,і= ро= 4л-ІО7 Г/м и бо=10_9/36я Ф/м, имеем

 

а — со

 

(2.20)

где с= 3 - ІО3 м/ с, е*= е/ео.

 

 

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

 

 

6 = (V/2)/р7~е

 

(2.21)

или

 

 

6 = 1 8 8 у //і* .

'

• (2.22)

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

При сое/ѵ<§Сі1, т. е. при преобладании токов проводимости, фа­ зовая постоянная и коэффициент поглощения численно равны:

а = b = у уцо)/2 .

(2.23)

В этом случае >е среды не влияет на поле.

Представляет интерес установить более четко пределы приме­ нимости приближенных выражений (2.19) —(2.22). Разлагая вряд

(2.17) и (2.18) по степеням ше/у и у/сое,

можно убедиться, что при­

ближенные зависимости (2.19) — (2.22)

позволяют определить а и

b с ошибкой не более 10%, если сое/у

1,1, а формула (2.23) обес­

печивает такую же точность при сое/у^0,2.

Фазовая постоянная и коэффициент поглощения определяют

длину волны X и глубину проникновения поля б. Длина волны со­

ответствует расстоянию, на котором фаза плоской волны меняется на 2л: Х = 2л/а.1

1 Здесь и далее имеется в виду плотность токов смещения и проводимости.

27


Глубина проникновения соответствует 'расстоянию, на котором по­ ле плоской волны убывает в е раз,

б = 1/6.

Из выражений (2.14) — (2.16) следует, что поле переменного магнитного диполя является комплексным и характеризуется в каждой точке двумя составляющими — реальной и мнимой, или амплитудным значением и фазой. В дальнейшем нас будет интере­ совать главным образом поле на оси магнитного диполя, которое в цилиндрических координатах в соответствии с выражением (2.15) записывается в виде

h2 = Hz/HZo =--eikz ( l - ik z ) ,

(2.24)

где HZo — вертикальная компонента поля в воздухе. Момент диполя направлен здесь по оси 2.

Идентичность выражений (2.14) и (2.15) позволяет утверждать, что выводы о поведении компоненты поля hz будут справедливы и в отношении компоненты электрического поля еч, если вместо z иметь в виду расстояние R.

Произведя в уравнении (2.24) разделение па реальную и мни­

мую части, имеем

 

 

 

Rehz =

e_iz[(l ~г bz) cosaz -\-azsinaz\\

(2.25)

Im hz =

e~bz [(1 -f- bz) sin az az cos az].

(2.26)

Амплитуда hz запишется в виде

 

 

\ h z \ = |

Hz/HZo I - e~bzY (1

+ b z f + (azr.

(2.27)

Фаза равна

 

 

 

Ф — az — arc tg [az/(1 +

6z)].

(2.28)

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

dl,

(2.29)

где dl — элемент'длины рамки.

Если поперечные размеры приемной катушки малы по сравне­

нию с расстоянием до генератора, т. е. z « R , то,

подставляя в

(2.29) выражение для £ф

(2.7),

получаем

 

і -

і оци^

- ë kz (1 — ikz),

(2.30)

 

2jiz3

 

где Л4ппЯп— момент приемной катушки; па — количество витков приемной катушки; sn — площадь витка.

28


Электродвижущая сила в воздухе при пренебрежении, токами смещения будет, очевидно, равна

 

с?о —

Mp/WrMn

 

(2.31)

 

2jiz3

 

 

 

 

 

Наведенная э. д. с. в приемной катушке зонда,

находящегося

в среде, выраженная

в относительных единицах,

будет

 

. <§/<?0 =

е^(1 -і7 гг).

 

(2.32)

Сравнивая уравнения (2.24)

н (2.32), видим,

что поле диполя

и сигнал в приемной

катушке

зонда, выраженные

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

в единицах поля или э. д. с. в воздухе, описываются идентичными выражениями. Аналитические выражения и расчетные графики для амплитуды поля |/iz| можно одновременно рассматривать как со­ ответствующие зависимости для амплитуды э. д. с. в приемной рамке зонда <§/<§о- Аналогично связаны реальная компонента по­ ля Re Іи и реактивная компонента э. д. с. <gреакт/ ё о; мнимая ком­ понента поля Im hz и активная компонента э. д. с. <§акт/(§о. Иден­ тичность (2.24) и (2.32) показывает, что фаза э. д. с. и фаза поля также описываются одинаковыми зависимостями.

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

Двухкатушечный зонд в однородной среде

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

Токи проводимости значительно превышают токи смещения (сое/уСІ). Диэлектрическая проницаемость среды не влияет на поведение поля. Действительная и мнимая части волнового числа равны и не зависят от е [см. (2.23)]. При обычно встречающихся значениях е и р горных пород этот случай охватывает широкий диа­ пазон частот, начиная от самых низких и кончая единицами мега­ герц. Выражения для амплитуды и фазы поля принимают сле­ дующий вид:

I К I

= е - ^ ] / (1 -Г рг)М -(рг)2;

(

(2-33)

ф =

pz — arc tg [pz/(l + pz)],

 

(2.34)

где p = a = b=У уцй)/2.

 

 

 

 

Амплитуда и фаза поля зависят здесь от параметра pz.

 

 

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

Рас­

А. Ближняя зона (область

малых параметров) |Az|->-0.

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

(длина зонда) значительно

меньше

длины волны в среде.

Разлагая в выражении (2.24) показатель­

29