![](/user_photo/_userpic.png)
MISCELLANEOUS / Geophysics / Geophysics Methods Костицын В. И
..pdf![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9171x1.jpg)
Электроразведка
милливольт у руд с электронопроводящими минералами (сульфиды, графит, антрацит). В целом зависит от многих природных факторов (минерального состава, глинистости, пористости, проницаемости, влажности, минерализации подземных вод и др.).
Способность пород поляризоваться, т.е. накапливать за-
ряд при пропускании тока, а затем разряжаться после отключения этого тока, оценивается коэффициентом поляризуемости . Величина вычисляется в процентах как отношение напряжения, которое остается в измерительной линии MN по истечении определенного времени (обычно 0,5–1 с) после размывания токовой цепи (ΔUВП), к напряжению в той же линии при пропускании тока ( U), т.е.
UВП 100% / U . |
(3. |
|
4) |
Поляризация – это сложный электрохимический процесс, протекающий при пропускании через породу постоянного или низкочастотного переменного (до 10 Гц) тока. Наибольшей поляризуемостью ( = 10–40%) отличаются руды с электронной проводимостью (сульфиды, сульфосоли, некоторые самородные металлы, отдельные окислы, графит, антрацит). Природа этих потенциалов ВП связана с так называемой концентрационной и электродной поляризацией рудных минералов. Поляризуемость до 2–6% имеют обводненные рыхлые осадочные породы с глинистыми частицами. Поляризуемость их обусловлена деформациями внешних обкладок двойных электрических слоев, возникающих на контакте твердой и жидкой фаз. Большинство изверженных, метаморфических и осадочных пород, насыщенных минеральной водой, слабо поляризуются ( < 2 %).
3.8.3. Диэлектрическая и магнитная проницаемость
Диэлектрическая ( ) и магнитная ( ) проницаемости
играют значительную роль в образовании электрических полей лишь при высоких частотах. Относительная диэлектрическая
169
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9172x1.jpg)
В.К. Хмелевской, В.И. Костицын
проницаемость = п / 0 (где п и 0 – диэлектрические проницаемости породы и воздуха) показывает, во сколько раз увеличивается емкость конденсатора, если вместо воздуха в него поместить данную породу. Величина
меняется от нескольких единиц (у сухих осадочных пород) до 80 (у воды) и зависит, в основном, от процентного содержания воды и от минералогического состава породы. У изверженных пород
меняется от 5 до 12 единиц, у осадочных – от 2–3 (у сухой) до 16–40 (у полностью насыщенной водой породы). Как отмечалось выше, магнитная проницаемость подавляющего большинства пород равна магнитной проницаемости воздуха. Лишь у ферромагнетиков относительная магнитная проницаемость может возрастать до 10 единиц (см. 2.7).
3.9. Общие принципы решения прямых задач электроразведки
В основе теории электроразведки лежат уравнения Максвелла, используемые в макроскопической электродинамике. Они включают все основные законы электромагнетизма (законы Ома, Ампера, Кирхгофа и др.) и описывают поля в разных средах. Из уравнений Максвелла получается дифференциальное уравнение, названное телеграфным. Решая его, можно получить электрическую (E) компоненту поля в средах вдали от источника с электромагнитными параметрами , , :
E |
|
|
|
E |
|
|
|
|
2 E |
, |
|
|
||||
|
|
|
|
t |
|
|
|
|
t 2 |
|
|
(3. |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
2 |
E |
|
2 |
E |
|
|
2 |
E |
|
5) |
||||
|
|
|
|
|||||||||||||
где E |
|
|
|
|
|
. |
|
|||||||||
|
x2 |
|
y 2 |
|
|
|
||||||||||
|
|
|
z 2 |
|
Дифференцирование ведется по декартовым координатам (х, у, z) и времени (t). Уравнение для магнитной (H) компоненты поля аналогично.
Если геоэлектрический разрез известен, то с помощью уравнения (3.5) и физических условий задачи, называемых ус-
170
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9173x1.jpg)
Электроразведка
ловиями сопряжения, решаются прямые задачи электроразведки, т.е. получаются аналитические или численные значения E и H, которые соответствуют заданному геоэлектрическому разрезу. В теории электроразведки прямые задачи решаются для разных физико-геологических моделей (ФГМ) сред.
Под ФГМ понимаются абстрактные геоэлектрические разрезы простой геометрической формы, которыми аппроксимируются реальные геолого-геофизические разрезы. Сложность решения прямых задач заключается в выборе моделей, близких к реальным, но таких, чтобы для избранного типа первичного поля удалось получить хотя бы приближенное решение для E или H. Для этого применяется математическое моделирование с использованием современных компьютеров. В недалеком прошлом основным способом решения прямых задач для сложных ФГМ и разных по структуре типов полей являлось физическое моделирование на объемных или плоскостных моделях сред.
Наиболее простые модели геологических сред:
однородное изотропное пространство или полупространство с одинаковыми электромагнитными свойствами (решения над ними называются соответственно первичным или нормальным полем источника);
анизотропное пространство или полупространство с электромагнитными свойствами, различающимися в направлении и вкрест слоистости пород;
одномерные неоднородные среды, в которых свойства меняются в одном направлении (например, вертикальные контакты двух сред, ряд вертикальных пластов или горизонтально слоистая среда с разными );
двухмерные неоднородные среды, в которых электромагнитные свойства меняются в двух направлениях (например, наклонные пласты или цилиндры, простирающиеся вдоль одного направления и отличающиеся по от вмещающих горных пород);
трехмерные неоднородные среды, в которых свойства меняются по трем направлениям (самой простой из подобных моделей является шар с разными , или в однородном полупространстве).
171
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9174x1.jpg)
В.К. Хмелевской, В.И. Костицын
В порядке увеличения сложности структуры первичных полей, а значит, возрастания сложности решения прямых задач, используемые для электроразведки поля можно расположить в следующей последовательности: точечные и дипольные источники постоянного тока, плоские гармонические электромагнитные волны, сферические волны дипольных гармонических или импульсных источников, цилиндрические волны длинного кабеля и т.п.
Существуют различные подходы к решению прямых задач с помощью уравнения (3.5). Любое правильное решение, удовлетворяющее всем физическим требованиям, единственно и корректно. Под корректностью понимается такое решение, в котором малым изменениям исходных данных соответствуют малые приращения расчетных параметров.
3.10. Нормальные поля в электроразведке
Под нормальным полeм понимается электромагнитное поле того или иного источника над однородным изотропным полупространством с неизменными электромагнитными свойствами. Из прямой задачи для поля точечного источника постоянного тока на земной поверхности можно получить нормальные поля постоянных электрических токов для разных установок или комбинаций питающих (АВ) и приемных (МN) электродов.
В электроразведке чаще всего применяются четырехэлектродные установки АМNВ (рис. 3.2). К одному питающему электроду (например А) подключается положительный полюс источника тока, к другому (В) – отрицательный. Разность потенциалов на приемных электродах (МN) от электрода А, определенная по полученной выше формуле (3.1), равна
U A |
J |
1 |
1 |
. |
||
|
|
|
|
|
||
2 |
|
AM |
|
AN |
||
|
|
|
|
Аналогичным образом можно получить разность потенциалов от отрицательного полюса В, но величину тока следует принять равной (- J).
172
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9175x1.jpg)
Электроразведка
Рис. 3.2. Расположение питающих (А и В) и приемных (М и N) электродов в разных установках метода сопротивлений: а – четырехэлектродной, б – срединного градиента, в – симметричной четырехэлектродной, г – трехэлектродной, д – двухэлектродной, е – дипольной радиальной, ж – дипольной азимутальной
Разность потенциалов от обоих электродов АВ равна су-
перпозиции |
|
U A и |
U В : |
|
|
|
|
|
|
|||
U |
J |
1 |
1 |
1 |
1 |
. |
( |
|||||
2 |
|
|
AM |
|
AN |
|
BM |
|
BN |
3.6) |
Если МN установить в центре АВ так, чтобы АМ = BN, АN = ВМ, то получим формулу для расчета сопротивления при работе с симметричной четырехэлектродной установкой:
AM AN |
|
U |
(3. |
|
|
|
|
. |
7) |
MN |
J |
Потенциал двухэлектродной установки АМ (А и N удалены в бесконечность) можно получить из формулы (3.6), приняв
BN , т.е.
173
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9176x1.jpg)
В.К. Хмелевской, В.И. Костицын
U Jr / 2AM .
В методах сопротивлений применяется и ряд других установок. Для глубинных исследований используются различные дипольные установки (рис. 3.2). Если приемный диполь (МN) перпендикулярен радиусу (r) между его центром и центром питающего диполя (АВ), а угол между радиусом и питающей линией АВ (θ) находится в пределах 70º θ 110º, то такая установка называется азимутальной. Частным случаем азимуталь-
ной является экваториальная установка (θ = 90º). Если прием-
ный диполь (МN) направлен вдоль r, а -30º θ +30º, то такая установка называется радиальной. Частным случаем радиальной установки является осевая (θ = 0º).
Для каждой установки можно получить формулы, по которым рассчитывается коэффициент установки. Так, для азимутальной установки
K 2 r 3 / AB MN
q ,
для радиальной –
K r 3 / AB
MN
p ,
где p и q – коэффициенты, мало отличающиеся от единицы и определяемые по специальным номограммам.
При работе с любой установкой рассчитывается по формуле для нормального поля:
K |
U |
, |
(3. |
J |
8) |
где U – разность потенциалов на МN; J – ток в АВ; а K – коэффициент устaновки, зависящий лишь от расстояний между электродами.
174
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9177x1.jpg)
Электроразведка
По этим же формулам можно рассчитать некоторое сопротивление над реальным, неизвестным и практически всегда неоднородным полупространством. Тогда оно называется ка-
жущимся сопротивлением (КС или к).
Расчет нормальных полей для других источников (гармонических, импульсных) более сложен, но в любом случае принято изучать кажущееся сопротивление.
3.11. Электрическое поле точечного источника постоянного тока над двухслойной средой
Простейшей прямой задачей для электроразведки ме-
тодом сопротивлений является задача об электрическом поле и кажущемся сопротивлении на поверхности полупространства, представляющего собой верхнее полупространство – воздух, а нижнее – двухслойную горизонтально слоистую среду с мощностью верхнего слоя h1, нижнего h2 = , УЭС слоев 1, 2 и 3 =
(воздух) (рис. 3.3).
Поставленная задача может быть решена с помощью уравнения (3.2), которое при f = 0 превращается в уравнение Лапласа U = 0, где U – потенциал в любой точке М с напряженностью электрического поля E U / r .
Рис. 3.3. Решение прямой задачи о поле точечного источника постоянного тока над двухслойной средой методом зеркальных отражений
175
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9178x1.jpg)
В.К. Хмелевской, В.И. Костицын
Однако ее можно быстро решить методом зеркальных отражений. Согласно правилам метода зеркальных отражений уравнение Лапласа и физические требования, в том числе граничные условия, выполняются, если потенциал в одномерной среде, где расположен точечный источник, принять равным сумме потенциалов этого источника (A) и всех его многократных отражений от границ раздела (A1, A2, A3,…) с коэффициен-
тами отражений, равными на границе I K12 |
2 |
1 / 2 1 , |
|
а на границе II K13 |
3 1 / 3 1 1 (так как 3 |
). |
На рис. 3.3 показано расположение источников при следующих обозначениях:
AM r , R R |
2 |
|
r 2 2h |
2 |
, |
|
|
||||||
|
|
1 |
|
|
|
|
1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
R |
4 |
r 2 |
|
4h |
2 ,... R |
n |
R |
n 1 |
r 2 2nh 2 |
, |
||
3 |
|
|
|
1 |
|
|
|
|
1 |
|
где |
n = 1, 2, 3,…, . |
|
|
|
|
|
|
|||||||||
|
Таким образом, искомое выражение для потенциала имеет |
|||||||||||||||
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U |
|
J |
1 1 |
2 |
|
|
|
K12n |
|
. |
(3 |
||||
|
|
2 |
|
|
r |
|
|
|
|
|
.9) |
|||||
|
|
|
|
|
n 1 r 2 (2nh1 )2 |
|
|
|||||||||
|
Выражение для кажущегося сопротивления – формуле |
|||||||||||||||
(3.1) – можно записать в виде |
|
|
|
|||||||||||||
|
|
|
2 |
|
r 2 |
|
U |
, |
|
|
|
|
||||
|
|
K |
J |
|
MN |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где |
U / MN |
|
E – напряженность электрического поля. |
|
||||||||||||
|
Но E |
|
|
U / r , поэтому |
|
|
|
|||||||||
|
|
2 |
|
r 2 |
|
U |
. |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
K |
J |
|
r |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
176 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9179x1.jpg)
Электроразведка
|
Подставив в эту формулу производную |
U / |
r из выра- |
|||||||||||||||||||||
жения (3.9), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
2 |
r 2 |
|
J |
1 |
|
1 |
|
|
|
|
|
|
K n r |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
12 |
|
|
. |
|
|
|
|
||
|
K |
|
J |
2 |
|
|
|
r 2 |
|
|
|
|
|
(2nh1 )2 3 / 2 |
|
|
|
|
||||||
|
|
|
|
|
|
n 1 r 2 |
|
|
|
|
|
|||||||||||||
|
Откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 n 1 |
|
|
K12n r 3 |
|
|
|
|
. |
|
|
|
|
|
|
|
(3. |
||||
|
K |
1 |
r 2 |
(2nh1 )2 |
3 / 2 |
|
|
|
|
|
|
|
|
|
10) |
|||||||||
|
Анализируя эту формулу, можно найти асимптотические |
|||||||||||||||||||||||
выражения |
|
K , равные |
1 и |
|
2 . В самом деле, |
при r |
0 будет |
|||||||||||||||||
K |
1 , при r |
|
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
1 |
2 |
|
K n |
|
|
|
1 |
2 |
|
|
K12 |
|
, |
|
|
|
|
|
|
||
|
K |
1 |
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|||||||||
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
n 1 |
|
|
|
|
|
|
1 |
K12 |
|
|
|
|
|
|
|
|||||
так как K |
|
2 |
1 |
/ |
2 |
|
1 |
, |
а |
K n равна |
K |
12 |
/ 1 |
K |
12 |
как |
||||||||
|
12 |
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
сумма членов геометрической прогрессии.
С помощью формулы (3.10), справедливой для трехэлектродной и симметричной четырехэлектродной градиентустановок, принято строить теоретические двухслойные кривые
– графики зависимости lg( K / 1 ) от lg r / h1 . Они называются
двухслойными теоретическими кривыми ВЭЗ (вертикальное электрическое зондирование), или двухслойной палеткой ВЭЗ
(рис. 3.4).
Более громоздкое решение получается в задаче о поле точечного источника над многослойной горизонтально-слоистой средой, а еще сложнее решение для такой же среды, но при возбуждении поля дипольными гармоническими или импульсными источниками.
177
![](/html/55006/802/html_5GD4XIfrCx.Xtj6/htmlconvd-6ez8I9180x1.jpg)
В.К. Хмелевской, В.И. Костицын
Рис. 3.4. Двухслойная палетка ВЭЗ: 1 и 2 – теоретические и полевая кривые
Одномерные прямые задачи электроразведки для много-
слойных горизонтально слоистых сред при любых первичных полях все-таки сводятся к аналитическим формулам при расчете КС. В результате принято строить кривые КС, аналогичные приведенным на рис. 3.4.
Двухмерные и трехмерные прямые задачи электрораз-
ведки решаются с помощью аналитических формул лишь для тел простой формы (шар, пласт, цилиндр) в однородной среде. В более общих случаях получаются лишь приближенные численные решения, найденные с помощью компьютеров.
3.12. Принципы решения обратных задач электроразведки
Накопленный материал по физическому и математическому моделированию прямых задач электроразведки позволил создать методы решения обратных задач, т.е. определять те или иные параметры геоэлектрического разреза по наблюденным графикам E, H или, например, кривым КС. Решение обратных задач неоднозначно в силу некорректности, как и всех обратных задач математической физики. Некорректность проявляется в
178