книги из ГПНТБ / Якубовский, Ю. В. Электроразведка учебник
.pdfГлава XII
ИНДУКТИВНЫЕ МЕТОДЫ ЭЛЕКТРОРАЗВЕДКИ
§1. ОБЩИЕ СВЕДЕНИЯ ОБ ИНДУКТИВНЫХ МЕТОДАХ
Киндуктивным относятся такие методы электроразведки, в кото рых электромагнитное поле создается при помощи незаземленных контуров, обтекаемых переменным током. В качестве таких контуров могут служить рамки малых размеров, большие прямоугольные петли, очень длинные кабели и др. Электромагнитное поле этих источников (обычно его магнитная компонента) исследуется на днев ной поверхности, в скважинах или над поверхностью земли (назем ные, скважинные и аэроэлектроразведочные варианты индуктив ных методов).
Другой особенностью индуктивных методов является то, что расстояние между источником поля и точками его наблюдения выби рают таким, чтобы оно было меньше длины волны, т. е. поле иссле дуется в ближней (индукционной) зоне источника поля.
Рис. 193 иллюстрирует сущность тех явлений, которые опре деляют характер переменных электромагнитных полей, наблюдаемых в индуктивных методах электроразведки.
Первичное переменное поле источника индуцирует в проводящем теоэлектрическом разрезе в и х р е в ы е т о к и , интенсивность которых определяется характером геоэлектрического разреза, типом и расположением источника, а также зависимостью поля от времени. Схематично можно выделить токи, наведенные в покровных отло жениях, в рудовмещающих породах и в локальных рудных залежах. Магнитные поля вихревых токов накладываются на первичное поле незаземленного контура и являются одним из источников наблюда емых аномалий. При рудопоисковых исследованиях особый интерес представляют те аномалии, которые связаны с вихревыми токами в рудных залежах. При геологическом картировании большую роль играют аномалии от вихревых токов во вмещающих покровных отложениях.
В том случае, когда вихревые токи, наведенные в проводящем разрезе, пересекают поверхности раздела сред с различным удельным
311
сопротивлением, на этих поверхностях накапливаются электриче ские заряды, поле которых обусловливает перераспределение токов в окружающей среде. В частности, токи в окрестностях хорошо проводящих рудных тел концентрируются в этих телах подобно тому, как это происходит в постоянных полях (см. гл. Ill), токи в окрестностях плохо проводящих тел обтекают последние и т. д. Аномалии, связанные с подобным перераспределением токов, при нято называть к о н ц е н т р а ц и о н н ы м и .
Если горные породы и руды, слагающие геоэлектрический разрез, обладают повышенной магнитной восприимчивостью, то, намагни-
Рис. 193. Проводящий магнитный разрез в пере менном индуктивно-возбу жденном поле.
1 — покровные |
отложения; |
|
2 — рудовмещающие |
поро |
|
ды; 3 — проводящая |
маг |
|
нитная рудная |
залежь; 4 — |
|
источник первичного |
поля; |
|
силовые линии: а — первич ного поля, б— вихревых то
ков в рудном теле; в — вих
ревых токов |
в покровных |
отложениях; |
г — вихревых |
токов в рудовмещающих по родах; д — индуцированных магнитных зарядов.
чиваясь в первичном магнитном поле, они создают аномалии, которые мы в дальнейшем будем называть м а г н и т о с т а т и ч е с к и м и .
Аномалии вихревого, концентрационного и магнитостатического типов сложным образом накладываются друг на друга, создавая суммарное аномальное поле. Специальным выбором типа источника поля, его зависимости от времени и измеряемого параметра поля можно усилить одни аномальные эффекты, подавить другие и таким образом специализировать модификации индуктивных методов для решения тех или иных геологических задач.
Теория индуктивных методов электроразведки строится на основе расчетов и моделирования электромагнитных полей при наличии проводящих геоэлектрических разрезов. Вследствие математических затруднений расчеты удается выполнить лишь для простейших разрезов. В более сложных случаях приходится прибегать к модели рованию.
Для теории рудопоисковых модификаций индуктивных методов большую роль играют расчеты переменных электромагнитных полей при наличии локальных проводящих тел, для картировочных моди фикаций — расчеты полей в присутствии однородного или горизонтально-слоистого разреза.
312
Пусть в первично-однородном магнитном поле Н0, меняющемся с частотой со, находится проводящее тело сферической формы, обла дающее проводимостью у и радиусом а (рис. 194). Проводимость вмещающей среды будем полагать равной нулю. Изменение магнит ного потока через любой контур внутри проводящей сферы приводит к появлению э. д. с. в этом контуре и соответственно к току, сила которого прямо пропорциональна э. д. с. в контуре и обратно про порциональна его сопротивлению.
Вихревые токи в сфере циркулируют по окружностям, располо женным в плоскостях, перпендикулярных к направлению первичного поля и имеющих центры на оси симмет рии, совпадающей с направлением пер вичного поля.
Величина суммарного поля, наблю даемого в точке с координатами г и Ѳ, может быть вычислена по следующей формуле, которая здесь приводится без вывода:
Hr= Hx0 (1 —-О/т]3) cos Ѳ,
#ѳ = -tf,o(l+ £ >/2if)sin0;
здесь |
Hx0 — первичное |
магнитное |
|
|
|||
поле; |
г и Ѳ — сферические |
координаты |
|
|
|||
точки, |
для |
которой |
рассчитывается |
Рис. 194. Проводящая сфера в пер |
|||
поле; ц = г]а — расстояние |
от |
точки |
вично-однородном |
переменном маг |
|||
наблюдения |
до центра |
сферы, |
выра |
нитном |
нож*. |
||
женное в радиусах сферы; D — коэффи |
|
|
|||||
циент, |
зависящий от частоты |
поля, проводимости, магнитной про |
|||||
ницаемости сферы и ее радиуса.
Из приведенных выражений следует, что аномальные компоненты поля при наличии сферы
Нга = НГ—Нх0cos Q= HXо cos Ѳ,
Я Ѳа = # ѳ + н *оsin Ѳ= Нхо 2 ^- Sin Ѳ
совпадают с полем магнитного диполя, помещенного в центр сферы и имеющего момент, прямо пропорциональный первичному полю и зависящий от коэффициента D. Физический смысл этого коэффи циента можно выяснить, если принять Ѳ = 0 и г = а, т. е. располо жить точку наблюдения на поверхности сферы. В этом случае
Яѳ = 0, Нга/Нх0= D,
т.е. функция D равна аномальному полю на поверхности сферы. График зависимости этой функции от безразмерного параметра
ра2 = усоца2 изображен на рис. 195. На графике приведена также зависимость фазы суммарного поля от того же параметра.
313
При постоянных величинах сферы и ее проводимости графики на рис. 195 характеризуют зависимость аномального поля и фазы суммарного поля от его частоты. При достаточно низких частотах (в области ра" <( 3) амплитуда аномального поля и фаза суммарного поля растут приближенно пропорционально частоте. Это объяс няется тем, что э. д. с., наведенная в любом контуре внутри сферы, растет пропорционально частоте, а скин-эффект внутри проводя щего тела при столь низких частотах еще не сказывается. G повыше нием частоты рост амплитуды аномального поля замедляется
и в области достаточно высоких частот (при ра2 > 1 2 0 ) она практи чески остается постоянной и равной амплитуде первичного поля (D = 1) на поверхности сферы. Эта весьма характерная особенность в частотной зависимости аномалии объясняется тем, что с ростом частоты скин-эффект приводит к выжиманию вихревых токов к по верхности сферы, т. е. к уменьшению сечения проводника и таким образом к увеличению его сопротивления. В результате рост э. д. с. с частотой компенсируется увеличением сопротивления и плотность вихревых токов, а следовательно, и аномальное поле перестают зависеть от частоты.
Фаза суммарного поля достигает максимума в области ра2 «=* 10 и затем по мере роста частоты уменьшается. Это объясняется тем, что в области достаточно высоких частот индуктивное сопротивление сферы значительно больше ее активного сопротивления.
Расчет и моделирование переменных электромагнитных полей для тел более сложной формы (цилиндра, пласта и др.) свидетель ствуют о том, что частотная зависимость аномалии от этих тел имеет общий характер, подобный описанному выше для сферы. Оказалось,
314
что частотные характеристики аномалий от проводящих тел различ ной формы, построенные в полулогарифмическом масштабе, при близительно совпадают, если по оси абсцисс откладывать пара метр pQ2, где Q — некоторый характерный размер, равный для сферы
радиусу, деленному на ]/2, для цилиндра — радиусу, а для пласта — произведению мощности на длину по падению.
|%л|%
Рис. 196. Обобщенная параметрическая характеристика,
о — амплитуда вторичного поля для тел: I — с изометрическим сечением, I I — вытянутых по падению; б — тангенс фазы вто ричного поля tg срвт: I — для пластины, II — для цилиндра,
шара-
График зависимости амплитуды и фазы аномального поля от
параметра pQ2 |
полупил название |
о б о б щ е н н о й |
п а р а м е |
т р и ч е с к о й |
х а р а к т е р и с т и к и а н о м а л и и |
(рис. 196). |
|
Описанные |
выше особенности |
в частотных характеристиках |
|
аномалий позволяют выбрать оптимальную частоту при работе ин дуктивными методами электроразведки. Очевидно, что при поисках хорошо проводящих рудных тел частота должна быть такой, чтобы амплитудные и фазовые аномалии от этих тел были максимальны, а аномалии от неоднородностей во вмещающей сфере (зон дробления>
315
•покровных отложений и др.) — минимальны. Исходя из первого требования, рабочую частоту желательно выбирать такой, чтобы параметр pQ2 приблизительно равнялся 10—12. Так, например, если на исследуемой площади рудные залежи имеют изомерную форму, близкую к сфере с радиусом а = 1 0 0 м, а удельное сопро тивление руды равно 0,1 Ом-м, то оптимальную рабочую частоту можно определить из следующего соотношения:
|
pQ2^ 12; |
2_усора2^ 1 2 , |
||
т. е. |
4-12 |
|
||
СО; |
' 400; |
|||
10 • 1002 • 4я • ІО-7 |
||||
отсюда |
|
|||
|
|
|
||
|
/ = £ ~ |
| § ~ 60 Г* |
||
Повышение частоты сверх оптимальной практически не обусло вливает увеличение амплитуды аномалии, а ведет к уменьшению фазовых аномалий и, что весьма нежелательно, к возрастанию интен
сивности аномалий |
от |
хорошо проводящих покровных |
отложений |
и неоднородностей |
в |
рудовмещающих породах — зон |
дробления, |
графитизированных разностей пород и др. |
|
||
Ниже показано, что обобщенная характеристика может быть использована при интерпретации результатов полевых наблюдений для определения параметров источников наблюдаемых аномалий.
Предположим теперь, что первичное магнитное поле Нх0 в момент
времени |
t = 0 |
меняется |
по ступенчатому закону, т. е. при t За 0 |
HXQ > 0 , |
при t |
0 Hxo = |
0. |
Всоответствии с законом индукции в момент резкого изменения первичного поля в проводящей среде индуцируются вторичные токи, интенсивность и распределение которых таково, что они стремятся сохранить неизменным первичное поле внутри проводника.
Вначальный момент времени индуцированные токи циркулируют на поверхности сферы вдоль окружностей, перпендикулярных к оси х
и соосных ей. С течением времени токи проникают в глубь сферы и одновременно затухают вследствие расхода энергии на выделение тепла.
Математически описанный выше переходный процесс в сфери ческом проводнике может быть описан следующими выражениями:
Her (t) = Нхо |
cos QLC(at), |
(XII.2) |
|
Нео(t) = Y н х0 |
sin 0Lc (at), |
||
|
|||
где |
р-0'-(л£)гt |
|
|
|
(XII.3) |
||
^ <“') = « 2 |
% лк), ; |
||
к=1 |
|
|
|
а = 1 /ура2. |
(XII.4) |
||
316
Таким образом, и в данном случае поле вихревых токов, наведен ных в проводящей сфере, совпадает с полем диполя, помещенного в центре сферы с осью, направленной против первичного поля. Момент этого диполя убывает со временем по закону бесконечной суммы экспонент, входящей в выражение для переходной функ ции Lc. Скорость убывания характеризуется параметром а , име ющим размерность, обратную времени. Из выражения (XII.4) сле дует, что в хорошо проводящих телах вихревые токи затухают медленнее, чем в плохо проводящих.
График зависимости Lc от t изображен на рис. 197.
B/So
Рис. 197. График зависимости I, от t. |
Рис. 198. Обобщенная переходная характери- |
с |
стика. |
|
1 — для цилиндра; 2 — для сферы. |
Расчеты, выполненные для тел цилиндрической и сфероидальной формы, а также моделирование — для тел более сложной формы, свидетельствуют о сходстве зависимостей нестационарного поля от времени для локальных проводящих тел различной формы. Таким
образом, оказалось |
возможным построить |
о б о б щ е н н у ю |
п е р е х о д н у ю |
х а р а к т е р и с т и к у , |
описывающую не |
стационарное поле для тел различной формы, так, как это было сделано выше для гармонически меняющихся полей.
Обобщенная переходная характеристика (рис. 198) представляет собой зависимость переходной функции L от некоторого обобщенного параметра
^обобщ 1 /YH C?!
где Q — квадрат характерного размера тела, аналогичного введен ному выше размеру локального проводника для частотных характеристик.
Отметим, что для рудных залежей, сложенных хорошо проводя щими массивными рудами, параметр а обобщ в зависимости от раз мера залежи и ее проводимости колеблется от первых единиц до нескольких десятков с-1.
317
Выбор времени регистрации переходного процесса, оптимального для хорошо проводящих руд, определяется соотношением между неустановившимся сигналом от искомых рудных залежей и сигналом от вихревых токов, наводимых в момент ступенчатого изменения первичного поля во вмещающих и покровных отложениях. Расчеты, подтвержденные опытом полевых работ, свидетельствуют о том, что при поисках хорошо проводящих руд неустановившийся сигнал следует регистрировать в интервале от одной-двух до 20—30 мс.
§2. НИЗКОЧАСТОТНЫЕ ИНДУКТИВНЫЕ МЕТОДЫ
Вэту группу входят индуктивные методы, в которых приме няются поля, меняющиеся по гармоническому закону. Классифика цию низкочастотных методов принято осуществлять по типу применяемого в них источника поля.
Метод незаземленной петли (НП). При исследованиях методом незаземленной петли первичное поле создается витком провода в виде
прямоугольника со сто
Вронами от нескольких со тен метров до 1 — 2 км, расположенным на поверх ности земли. Магнитное поле петли измеряют на профилях, расположенных либо внутри петли, в ее центральной части, либо на профилях, располо женных вне петли перпен
Рис. 199. Основные |
модификации |
индуктивных ме |
дикулярно |
к |
ее |
длин |
|||
ной стороне (рис. 199, а). |
|||||||||
|
|
тодов. |
|
||||||
а — метод незаземленной |
петли; |
б — метод беско |
Обычно измеряют |
ампли |
|||||
нечно длинного |
кабеля; |
в — дипольно-индуктивное |
туду и фазу или только |
||||||
профилирование; |
г — дипольно-индуктивное зонди |
||||||||
рование; Г — генератор; |
П — приемник; К — ка |
амплитуду |
вертикальной |
||||||
|
бель; Н т |
— петля. |
составляющей поля, хотя |
||||||
не исключена |
возможность измерения |
||||||||
горизонтальней |
составля |
||||||||
ющей.
Метод бесконечно длинного кабеля (БДК). В этом методе источ ником поля служит длинный прямолинейный кабель, заземленный на концах (рис. 199, б). Магнитное поле кабеля изучают вдоль про филей, ориентированных перпендикулярно к кабелю или парал лельно ему. Длину профилей и их расстояние от кабеля выбирают таким образом, чтобы во всех точках наблюдения магнитное поле токов, введенных в землю через питающие электроды, было пре небрежимо мало по сравнению с магнитным полем токов, текущих в кабеле. Таким образом, несмотря на наличие заземлений в пита ющей линии, метод бесконечно длинного кабеля по существу индуктивный.
При работе методом БДК, так же как и при работе методом не заземленной петли, измеряют амплитуду и фазу или только ампли-
318
туду вертикальной и реже горизонтальной составляющей магнитного поля. В обоих методах применяют сходную генераторную и измери тельную аппаратуру.
Дипольное индуктивное профилирование (ДИП). Эта модификация индуктивных методов характеризуется тем, что источником первич ного поля служит магнитный диполь — многовитковая рамка не большого диаметра (около 1 м), обтекаемая переменным током. Магнитное поле, так же как и в других индуктивных методах, изме ряют при помощи приемной рамки (измерительного диполя). Оба диполя в процессе полевых работ перемещают вдоль профиля без изменения взаимного положения (рис. 199, в).
Дипольное индуктивное зондирование (ДИЗ). Отличается от про филирования тем, что изучается зависимость поля неподвижного
магнитного диполя |
от расстояния |
0 |
|
|
||||||
между ним и измерительным ди |
|
|
||||||||
полем (дистанционное |
зондирова |
|
|
|||||||
ние) |
или зависимость |
поля от ча |
|
|
|
|||||
стоты при |
неподвижных |
генера |
|
|
|
|||||
торном и измерительном |
диполях |
|
|
|
||||||
(частотное зондирование). |
магнит |
|
|
|
||||||
Способы |
измерения |
200. |
Блок-схема |
микровольтметра |
||||||
ного |
поля. |
При |
полевых |
рабо |
||||||
Рис.-П> |
ИМА. |
|||||||||
тах |
индуктивными методами при |
ОТНОСИТ е л ь н ы е способы |
||||||||
меняются |
а б с о л ю т н ы е |
и |
||||||||
измерения магнитного поля. |
|
|
|
|
||||||
Абсолютные измерения. При абсолютных измерениях определяют амплитуду и фазу, либо только амплитуду или фазу различных пространственных составляющих магнитного поля в точке измере ния. К абсолютным следует отнести измерение элементов эллипса поляризации вектора магнитного поля.
Для абсолютных измерений наиболее широко применяют а м п л и т у д н ы й в о л ь т м е т р . На рис. 200 приведена блок-схема вольтметра ИМА (измеритель магнитной амплитуды), собранного на унифицированных блоках электроразведочной аппаратуры, изго товляемых заводами геофизического приборостроения.
В качестве датчика 1 в вольтметре ИМА применена секциониро ванная катушка с ферритовым сердечником; при помощи надлежащей коммутации секций и подстроечных конденсаторов датчик может быть настроен на одну из четырех частот — 78, 312, 1250 и 2500 Гц. Датчик устанавливается на штативе и при помощи двухосного угло мерного устройства может быть ориентирован по любому выбранному заранее направлению, компоненту поля вдоль которого желательно измерить. Сигнал с датчика поступает на выносной усилитель 2, снабженный калиброванным делителем напряжения. Усилитель укреплен на штативе датчика. С выхода выносного усилителя сигнал по соединительному шлангу поступает на семиступенчатый делитель напряжений 3, а затем на усилитель 4, снабженный активным фильтром 5, настроенным на рабочие частоты. Фильтр позволяет
319
подавлять внешние помехи и внутренние шумы схемы. Амплитудный детектор 6 выпрямляет усиленный и отфильтрованный сигнал. Микроамперметр 7 предназначен для измерения этого сигнала, а также для контроля напряжения питания, которое осуществляется бата реей 8 через стабилизатор напряжения 9.
Втом случае, когда необходимо использовать поля, меняющиеся
счастотой менее 78 Гц, для абсолютных измерений амплитуды маг нитного поля может быть применен комплект аппаратуры ИКС-50, в котором измерительный прибор укомплектован магнитоиндук ционным датчиком.
Относительные измерения. При относительных измерениях
определяют отношение амплитуд и разность |
фаз магнитного поля |
в двух точках исследуемой площади или в |
одной и той же точке, |
И зм ерит ельны й п ульт |
|
Рис. 20І. Блок-схема афиметра АФИ-4.
но для двух пространственных составляющих. Для относительных измерений применяют а ф и м е т р ы (амплитудно-фазовые измери тели). Блок-схема афиметра АФИ-4 изображена на рис. 201.
Датчиками магнитного поля в афиметре служат две многовитковые рамки с ферритовыми сердечниками.
Измеряемые напряжения с датчиков поля 1 поступают на вынос ные усилители 2, расположенные около рамок или измерительных электродов, и затем на градуированные в децибеллах делители напряжения 3, при помощи которых амплитуда измеряемых сигналов делится без изменения их фазы. С делителей напряжения сигналы поступают на градуированные в градусах фазовращатели 4, которые изменяют фазы этих сигналов без изменения их амплитуды. Посту пающие с выхода фазовращателей сигналы смешиваются на сумма торе 5 и затем подаются на вход индикатора нуля, представляющего собой высокочувствительный избирательный усилитель 6 со стре лочным прибором 7 на выходе.
Процесс измерения с афиметром заключается в том, что при помощи делителей напряжения уравнивают амплитуды сигналов в обеих рамках; используя фазовращатели 4, сдвиг по фазе этих сигналов доводят до 180°. Поскольку на сумматор поступают два равных по амплитуде и противоположных по фазе сигнала, индика тор нуля отмечает момент уравнивания нулевым показанием стре лочного прибора (или отсутствием звука в телефоне). Отсчет
320
