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

книги из ГПНТБ / Якубовский, Ю. В. Электроразведка учебник

.pdf
Скачиваний:
30
Добавлен:
22.10.2023
Размер:
17.65 Mб
Скачать

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

Рис. 13. Пояснение к расчету поля точечного источника тока при наличии в разрезе плоской поверхности раздела.

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

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

31

Если бы нѳ было второй среды с удельным сопротивлением р2, заполняющей правое полупространство, то потенциал U х в точке М, ■связанный с источником тока I , был бы равен

Ui /рі/4ягуім. (II.9)

Наличие поверхности раздела двух сред с различным удельным сопротивлением обусловливает изменение величины потенциала в точке М. Это можно объяснить тем, что среда с удельным сопроти­ влением р2 влияет на плотность тока в области точки М. Если при этом р2 > рх, то плотность тока возрастет за счет того, что в объем, занятый средой с удельным сопротивлением р2, пойдет тока меньше, чем в том случае, если бы этот объем был занят средой с удельным сопротивлением рг Если же р2 <; рх, плотность тока уменьшится, так как среда с удельным сопротивлением р2 отвлечет на себя больше тока вследствие своей лучшей проводимости. Аналогичные рассужде­ ния можно привести и в отношении потенциала U2 в точке М , находящейся в среде с удельным сопротивлением р2.

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

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

раздела сред с разным удельным

сопротивлением, а также вблизи

и на большом удалении от источников.

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

выражаются двумя требованиями.

Во-первых, потенциалы в точках, бесконечно близко расположен­

ных к поверхности раздела двух сред, но по обе стороны от этой поверхности, должны быть равны, т. е.

и г = и л.

(11.10)

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

Во-вторых, должна соблюдаться непрерывность нормальной составляющей вектора плотности тока ]л при переходе из одной

среды

в другую:

 

Ома

 

 

 

В

соответствии с законом

 

 

 

 

In'-

 

 

L

ди

 

 

 

р

дп

где dU/dn — производная по направлению

нормали.

Тогда второе граничное условие можно записать следующим

образом:

dU1

 

дЦа

 

 

 

1

1

 

(И -11)

 

Р і

дп

р 2

дп

 

 

 

 

:32

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

источника потенциал вблизи него равен

т. е- при г

О U

-V оо как 1/г.

 

 

 

У с л о в и е

н а б е с к о н е ч н о с т и

означает, что

потен­

циал источника с возрастанием г стремится к нулю.

и з о ­

На этих положениях основан м е т о д з е р к а л ь н ы х

б р а ж е н и й .

Сущность его заключается

в том, что при

расчете

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

Таким образом, поле точечного источника в неоднородной среде рассчитывают как поле в однородной среде, но от совокупности

действительного

и

фиктивных источников. Последние

помещают

в

точках зеркального изображения действительного

источника

в

поверхностях

раздела.

 

 

 

 

Рассмотрим несколько примеров расчета полей в неоднородных

средах

при помощи

метода

зеркальных

изображений.

 

 

Поле

точечного источника

при наличии

в разрезе плоской по­

верхности раздела сред. Пусть источник тока I помещен в точке А среды с удельным сопротивлением р2 (рис. 13). Для учета влияния плоской поверхности раздела сред с удельным сопротивлением рх и р2 поместим в точку А', являющуюся зеркальным изображением точки А, фиктивный источник тока Г . Поле в среде с удельным сопро­ тивлением рх будем рассчитывать как поле обоих источников в без­

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

рх:

£7і = /рі/4яг-і-/,р1/4яг\

(11.12)

где г и г' — расстояния от точки М, в которой определяется поле, до точек А и А'.

Потенциал U2 в точках среды с удельным сопротивлением р2 будем рассчитывать как поле источника тока /", расположенного в точке А, причем вся среда считается однородной и обладает удель­ ным сопротивлением р2:

£/2 = /"р2/4яг"; (НДЗ)

здесь г" — расстояние от точки М ' до точки А.

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

3 Заказ 512

33

условия, а именно — при расположении точек М и М ' на поверхности раздела должны соблюдаться равенства

и

1

dU1

1

 

Рі

дп

р2 дп *

Подставим в эти уравнения выражения (11.12) и (11.13). При этом

примем во внимание,

что при совмещении точек М и М '

х = d и

г = г = г" = УX2 +

у 2 + z2. Кроме того, направление

нормали

к поверхности раздела совпадет с направлением оси х , поэтому для определения нормальной составляющей плотности тока нужно про­

дифференцировать потенциальную функцию по х,

т. е.

dU

dU

д

I

Ір \

 

Ір

X

 

Ip

d

дп

дх

дх \4лг )

г3

 

г3 *

Тогда получим

 

 

 

 

 

 

 

 

 

 

 

ірі/4лг +

F p1IAnr‘ = /"р2/4яг",

 

1

/рі d

.

1

/'pi

d ____1

7"p2

d

Pi

r3

T" p!

 

4л (r')3

p2

 

(r")3

или с учетом изложенного выше

 

 

 

 

 

 

/рі/4яг +

FpJAnr = /"р2/4лг,

 

___1

/Рі

d

I

1

/'pi

d

 

1

/"рз

r f .

 

Pi

r3

■ pi

г3 —

p2

r3 ’

здесь знак минус у второго члена меняется на плюс, так как г' имеет направление, противоположное направлению оси х.

После упрощения получаем уравнения

/рі-Ь /'Р і = /"р2.

Решая эти уравнения относительно неизвестных F и/", находим:

1" = (Ра— Рі)/(Ра + Рі) 1 = knl>

(II. 14)

Г =

2рі//(ря +

Pi) = (1 - klt) /,

(II. 15)

где

*i2 = (Pa—Pi)/(P2 + Pi)-

 

 

 

Безразмерный коэффициент

к называется

к о э ф ф и ц и е н ­

т о м о т р а ж е н и я .

Это название не совсем удачно, так как не

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

Цифровой индекс у коэффициента к 12 соответствует индексам контактирующих сред. Величина к 12 при любых значениях рх и р2 меняется в пределах от + 1 до —1, т. е. по абсолютной величине не превышает единицы.

34

Теперь мы можем написать выражения для потенциалов U1 и

U 2 и полупить искомое решение задачи:

 

Uг

 

 

(Н.16)

и а

Тр2

(1 —кц)

(И.17)

г"

 

 

Смысл полученного решения уясним на примере поля точечного источника, расположенного в среде с удельным сопротивлением рх вблизи контакта со средой с бесконечно высоким сопротивлением р2 =

=оо. В этом случае к 12 =

=1, т. е. фиктивный ис­ точник имеет ту же силу, что и действительный.

Если

точка

наблюдения

 

отстоит от действительного

 

источника

на расстоянии,

 

значительно

большем, чем

 

расстояние

 

между

этим

 

источником

и контактом,

 

или

если

источник

поля

 

расположен

на контакте,

 

то г -> г' и поле и г в точке

Рис. 14. Учет влияния границы земля — воздух.

наблюдения

согласно вы­

ражению

(11.16)

удваи­

 

вается. Во второй среде потенциал U2 будет равен нулю, т. е. ток в среду бесконечно высокого сопротивления не потечет.

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

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

характерную особенность, что у дневной

поверхности нормальная

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

равна нулю. Физически

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

3*

35

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

В том частном случае, когда источник расположен на поверх­ ности земли, его отражение совпадает с самим источником, и для учета влияния поверхности раздела земля — воздух силу тока, питающего источник, удваивают.

Поле точечного источника при наличии в разрезе вертикальной поверхности раздела сред. Пусть источник тока I расположен на по­ верхности земли в точке А на расстоянии d от вертикального кон­ такта сред с pj и р2. Данная задача отличается от предыдущей нали­ чием второй поверхности — границы земля — воздух. Отразим ниж­ нее полупространство в верхнее в этой поверхности и таким образом учтем ее влияние. Тогда условие задачи сведется к условиям преды­ дущей задачи и мы сможем для подсчета потенциала воспользоваться выражениями (11.16) и (II.17), подставив в них вместо тока / вели­

чину 21:

 

и * = І Н - г + - т Д .

(ІШ >

 

<I U 9 >

Расчетные формулы упрощаются, если профиль наблюдений провести через источник поля А перпендикулярно к линии контакта сред (рис. 15). Поместив начало координат в точку А и направив ось X по линии профиля, получим:

г т

І р 1

(

1 I

k l 2 ^

Ul

2 п

\

X _t_

(2d— x) )

TT

 

(1 — fcl2)

Uz -

 

X

Найдем также выражения для напряженности поля:

р _

d U i

_; /рі

,_1_______ ^12

~1

 

дх

2 л

Lx2

(2d— x)2_J*

F

djy2 _

Д>2

(1—^іг)

■^2 —

дх

2п

ж2

( 11. 20)

(11.21)

( 11. 22)

(11.23)

Формулы (11.20) и (11.22) позволяют рассчитывать поле в среде с р1( а формулы (11.21) и (11.23) — в среде с р2.

На основе общего решения задачи можно графически изобразить электрическое поле точечного источника, расположенного у поверх­

ности раздела двух

сред. На рис. 16 представлено поле для двух

случаев: р2< Р і и

ра > р г

Токовые линии на границе преломляются, причем для каждой линии соблюдается следующее условие:

Pi/Pa = tK «a/tg <*i»

где а 1 и а 2 — углы между направлением токовой линии и нормалью к границе в точке их преломления.

Р>. Рг

Рис. 15. Пояснение к расчету поля точечного источника тока, расположенного у поверхности Земли при наличии в разрезе вертикальной поверхности раздела.

В' E E h

Рис. 16. Электрическое поле точечного источника, расположенного у по­ верхности раздела сред с различным удельным сопротивлением.

а — р2 в pj/5; б — р2 = 5рі; линии: 1 — токовые, 2 — эквипотенциаль­ ные.

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

37

в верхней среде плотность тока возрастает (ас ней возрастает и раз­ ность потенциалов между точками М и N), во втором случае — умень­ шается (соответственно уменьшается и разность потенциалов между точками М и N).

Следует обратить внимание на то, что ток вследствие преломле­ ния токовых линий в плохо проводящем пласте стремится течь поперек пласта; переходя из плохо проводящего пласта в пласт высокой проводимости, он стремится течь вдоль пласта.

Поле точечного источника

при наличии

в

разрезе гори­

зонтальной поверхности раз­

дела сред. Пусть мы

имеем

двуслойный

с

геоэлектриче­

ский разрез

одной

гори­

зонтальной

 

поверхностью

раздела

 

сред

с сопротивле­

ниями р2

и р2 (рис. 17). Мощ­

ность пёрвого

слоя

примем

равной

h v

Рассчитаем

на

поверхности

земли

поле

то­

чечного

 

источника

тока

I,

расположенного

также

на

дневной

поверхности

в точ­

ке А.

 

 

всего

отразим

Прежде

нижнее

 

проводящее

полу­

пространство

в

 

поверхности

земля — воздух.

 

Получим

среду

с

двумя

поверхно­

стями раздела — истинной I

и фиктивной I I

— и с источ­

ником А, питающим полное пространство током силой 21. Будем определять потен­ циал в произвольной точке

М, находящейся от источника на расстоянии г, методом зеркаль­ ных изображений. Потенциал U1 в точке М будет создаваться точечным источником А и фиктивным источником А ѵ являющимся зеркальным отражением источника А в границе I.

Согласно формуле (11.14) источник А г должен

отдавать в среду

ток

 

 

 

Ij = 2//с12,

 

 

 

 

 

 

 

 

поэтому в соответствии с формулами (11.12) и (11.16)

Ur

2ірі

, Ірі

Ірі

IkiiPi

^ Ipx

I

IkiiPi____

4яг

4Jtr'

2яг

2nr'

~~ 2ixr

2n/r2 + (2M)2 *

 

83

Этим самым мы соблюли условия на границе /. Однако необ­ ходимо выполнить также условия на границе II. Для этого нужно отразить в ней не только источник А, но и источник А х (на рис. 17 эта операция отражения — вторая по счету — показана дугообраз­

ными стрелками

2).

Источник

А

своим зеркальным изображением

создаст фиктивный

источник

А и ,

расположенный от точки М

на

расстоянии гп и

от поверхности

наблюдения

на расстоянии 2hx\

в точке зеркального изображения фиктивного

источника A j

по­

явится еще один

фиктивный источник А П1 на

расстояниях гиІ

от

точки М и 4/zj от точки А.

 

 

 

 

Ток источника А п должен быть равен

 

 

а ток источника

 

І ц — ш іаі

 

 

А ІІ:

 

 

 

 

І ц і = І і к 12 — (2 //е12) к 12 = 2Ік і2 -

Но появившиеся источники А и и А ш нарушают условия на гра­ нице I. Для их соблюдения придется отразить источники А и и А іи в границе I (на рис. 17 — третья операция отражения обозначена стрелками 5); получим еще два фиктивных источника А ІѴ и А ѵ на расстояниях от точки М и поверхности наблюдения соответственно гІѴ, 4h x и rv , бh x. Токи источников А 1Ѵ и А ѵ будут

 

IIV —

~ (2 1 2 ) ^12 =

2Ік \ 2,

 

 

 

І у = І ц і к 12 — (2 Ік \ 2) к 12 =

2Ік \ 2.

 

 

Однако

теперь нарушены

условия на

границе

II.

Операцией

отражения

4 (рис. 17) они

восстанавливаются (источники A VJ и

Лг а ), но

одновременно снова нарушаются условия на границе /.

Проведем

пятую операцию отражения и т. д. Эти операции, оче­

видно, следует продолжать до бесконечности.

М

будет пред­

Таким образом, искомый потенциал в точке

ставлять бесконечную сумму потенциалов, создаваемых действи­

тельным источником А и

фиктивными

источниками Aj, А п ,

А ш ,

А п:

 

 

j j

__

7рі

j _______ ^ ia P i_______ I______Ik\2Pi______ I_________7fcj2Pi_________ I

1

 

2 л r

> 2л Y r2 + (2Ai)2

2я / г2 + (2/гх)2

/ r 2 + (2 • 2AX)2

I_________ ^ 1 2 !-11________г_______ IkiaPi________ I_______ ^ ' ' 1 2 Pi__________ |_

^

2 л У гг2 -і-(2 -2 к 1р ~ '~ 2 л Ѵ г 2 ^ ( 3 - 2 .Ь 1)2

]Л -2+ ( 3 • 2AX)2

_7pi_

&12

(11.24)

2 л

n=1 /r2+(2nÄX)2

 

 

 

Получим выражение для

напряженности поля,

для чего про­

дифференцируем почленно выражение (11.24) по г:

 

дѴ1

 

(

 

 

 

00

 

Е1

7рі

j

1

I

2

 

_____ /;,|І2Г_____ . (11.25)

дг

 

 

 

I

г®

"г*

 

£ [r 2.y(2refei)2]s /l

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

Аналогичным путем можно решить задачу для разреза с двумя

горизонтальными поверхностями

раздела. Получаемые при этом

 

с

 

выражения

для потенциала и на­

 

 

пряженности

поля

также имеют

 

 

 

 

 

 

вид сходящихся рядов, в целом

 

 

 

аналогичных

формулам (11.24) н

 

 

 

(11.25),

но

имеющих под знаком

 

 

 

суммирования более сложные по­

 

 

 

стоянные

множители,

зависящие

 

 

 

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

отражения к 12

 

 

 

и

к 23

и

мощностей слоев.

 

 

 

 

 

 

Поле

точечного

источника

 

 

 

в однородной анизотропной среде.

 

 

 

Однородные

 

и

неоднородные

сре­

 

 

 

ды, в которых

мы рассматривали

 

 

 

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

поля,

подразуме­

 

 

 

вались изотропными. Однако, как

 

 

 

указывалось

 

в

гл. I,

горные

по­

 

 

 

роды часто

 

обладают

микроани­

 

 

 

зотропией.

 

 

 

 

 

 

 

 

 

 

 

 

 

Влияние микроанизотропии по­

 

 

 

род на

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

поля

оце­

 

 

 

ним на примере

поля

точечного

 

 

 

источника. Выше было показано,

 

 

 

что в

однородной

 

среде это поле

 

 

 

характеризуется

 

центральной

 

 

 

(сферической)

симметрией,

т.

е.

 

 

 

токовые линии радиально

и

рав­

Рис. 18. Влияние анизотропии

на токо­

номерно

расходятся

во

все

 

сто­

 

вые линии.

 

роны

из

 

точки

 

расположения

верхности

 

 

источника, эквипотенциальные по­

являются семейством

концентрических

окружностей

(см. рис. 5). В прямоугольной системе координат х, у,

z с началом

в точке расположения источника эти особенности поля могут

быть

выражены

следующими

соотношениями:

 

 

 

 

 

 

 

 

 

 

 

1x1х = ІуІУ= 7z/z,

E J x = Ey/у = EJz.

 

 

 

 

 

В однородной безграничной анизотропной среде поле точечного источника по сравнению с полем в изотропной среде деформируется. Характер этой деформации можно видеть на рис. 18, на котором анизотропная среда представлена в виде параллельно переслаива­ ющихся слоев с сопротивлениями рх и р2. Токовые линии 1 и 2 точечного источника А, находящегося в слое с сопротивлением р2

40

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