![](/user_photo/_userpic.png)
книги из ГПНТБ / Якубовский, Ю. В. Электроразведка учебник
.pdfс которыми можно встретиться при изучении геоэлектрических раз резов различного типа. А это, в свою очередь, требует от нас знания расчетов полей различных источников в электрически неоднородных средах. Однако эти расчеты очень сложны и их удается выполнить лишь при некоторых упрощающих предположениях. Для этого реальный геологический разрез, каждая порода в котором в той или иной мере электрически неоднородна, заменяют идеализированным разрезом, состоящим из отдельных областей с электрически однород ными средами. Изменение свойств этих сред на их границах про исходит скачкообразно. Сложные поверхности, по которым контак тируют породы в разрезе, заменяют простыми — плоскими, сфериче-
Рис. 13. Пояснение к расчету поля точечного источника тока при наличии в разрезе плоской поверхности раздела.
сними, цилиндрическими и др. В одних случаях эти упрощения невызывают значительных отклонений этих разрезов от реального гео логического (например, при замене разреза с горизонтальными или слабо наклоненными пластами осадочных пород горизонтально слоистым геоэлектрическим разрезом) и проводимые расчеты поз воляют количественно определить мощность и глубину залегания отдельных пластов. В других случаях при сложной геологической обстановке расчеты полей, проводимые на основе указанных выше упрощений, дают возможность лишь качественно судить о место положении, размере или характере того геологического объекта, который обусловил наблюденную аномалию электрического поля. Однако уже и этого качественного суждения достаточно для понима ния результатов электрической разведки.
Характер влияния неоднородностей разреза на электрические поля источников, а также особенности расчета полей в неоднород ных средах рассмотрим на следующем примере. Пусть в среде с удельным сопротивлением рх, заполняющей левое полупространство, в точке А расположен источник тока силой / (рис. 13).
31
Если бы нѳ было второй среды с удельным сопротивлением р2, заполняющей правое полупространство, то потенциал U х в точке М, ■связанный с источником тока I , был бы равен
Ui —/рі/4ягуім. (II.9)
Наличие поверхности раздела двух сред с различным удельным сопротивлением обусловливает изменение величины потенциала Uг в точке М. Это можно объяснить тем, что среда с удельным сопроти влением р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лг ) |
4л |
г3 |
|
4л |
г3 * |
||||
Тогда получим |
|
|
|
|
|
|
|
|
|
|
|
|
ірі/4лг + |
F p1IAnr‘ = /"р2/4яг", |
|
||||||||
1 |
/рі d |
. |
1 |
/'pi |
d ____1 |
7"p2 |
d |
||||
Pi |
4л |
r3 |
T" p! |
|
4л (r')3 |
p2 |
|
4л |
(r")3 |
||
или с учетом изложенного выше |
|
|
|
|
|
||||||
|
/рі/4яг + |
FpJAnr = /"р2/4лг, |
|
||||||||
___1 |
/Рі |
d |
I |
1 |
/'pi |
d |
|
1 |
/"рз |
r f . |
|
|
Pi |
4л |
r3 |
■ pi |
4л |
г3 — |
p2 |
4л |
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) |
|
4л |
г" |
|||
|
|
Смысл полученного решения уясним на примере поля точечного источника, расположенного в среде с удельным сопротивлением рх вблизи контакта со средой с бесконечно высоким сопротивлением р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 - |
|
2Я |
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
![](/html/65386/283/html_d9uWUGbaPj.EeE3/htmlconvd-cT5R9x38x1.jpg)
в верхней среде плотность тока возрастает (ас ней возрастает и раз ность потенциалов между точками М и 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 и |
|||
Лг а ), но |
одновременно снова нарушаются условия на границе /. |
||||
Проведем |
пятую операцию отражения и т. д. Эти операции, оче |
||||
видно, следует продолжать до бесконечности. |
М |
будет пред |
|||
Таким образом, искомый потенциал Uг в точке |
ставлять бесконечную сумму потенциалов, создаваемых действи
тельным источником А и |
фиктивными |
источниками 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 |
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) |
|
дг |
|
|
||||||
|
2л |
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