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