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

MISCELLANEOUS / Geophysics / Geophysics Methods Костицын В. И

..pdf
Скачиваний:
6
Добавлен:
03.05.2023
Размер:
6.07 Mб
Скачать

Электроразведка

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

Способность пород поляризоваться, т.е. накапливать за-

ряд при пропускании тока, а затем разряжаться после отключения этого тока, оценивается коэффициентом поляризуемости . Величина вычисляется в процентах как отношение напряжения, которое остается в измерительной линии MN по истечении определенного времени (обычно 0,5–1 с) после размывания токовой цепи (ΔUВП), к напряжению в той же линии при пропускании тока ( U), т.е.

UВП 100% / U .

(3.

 

4)

Поляризация – это сложный электрохимический процесс, протекающий при пропускании через породу постоянного или низкочастотного переменного (до 10 Гц) тока. Наибольшей поляризуемостью ( = 10–40%) отличаются руды с электронной проводимостью (сульфиды, сульфосоли, некоторые самородные металлы, отдельные окислы, графит, антрацит). Природа этих потенциалов ВП связана с так называемой концентрационной и электродной поляризацией рудных минералов. Поляризуемость до 2–6% имеют обводненные рыхлые осадочные породы с глинистыми частицами. Поляризуемость их обусловлена деформациями внешних обкладок двойных электрических слоев, возникающих на контакте твердой и жидкой фаз. Большинство изверженных, метаморфических и осадочных пород, насыщенных минеральной водой, слабо поляризуются ( < 2 %).

3.8.3. Диэлектрическая и магнитная проницаемость

Диэлектрическая ( ) и магнитная ( ) проницаемости

играют значительную роль в образовании электрических полей лишь при высоких частотах. Относительная диэлектрическая

169

В.К. Хмелевской, В.И. Костицын

проницаемость = п / 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

Электроразведка

ловиями сопряжения, решаются прямые задачи электроразведки, т.е. получаются аналитические или численные значения E и H, которые соответствуют заданному геоэлектрическому разрезу. В теории электроразведки прямые задачи решаются для разных физико-геологических моделей (ФГМ) сред.

Под ФГМ понимаются абстрактные геоэлектрические разрезы простой геометрической формы, которыми аппроксимируются реальные геолого-геофизические разрезы. Сложность решения прямых задач заключается в выборе моделей, близких к реальным, но таких, чтобы для избранного типа первичного поля удалось получить хотя бы приближенное решение для E или H. Для этого применяется математическое моделирование с использованием современных компьютеров. В недалеком прошлом основным способом решения прямых задач для сложных ФГМ и разных по структуре типов полей являлось физическое моделирование на объемных или плоскостных моделях сред.

Наиболее простые модели геологических сред:

однородное изотропное пространство или полупространство с одинаковыми электромагнитными свойствами (решения над ними называются соответственно первичным или нормальным полем источника);

анизотропное пространство или полупространство с электромагнитными свойствами, различающимися в направлении и вкрест слоистости пород;

одномерные неоднородные среды, в которых свойства меняются в одном направлении (например, вертикальные контакты двух сред, ряд вертикальных пластов или горизонтально слоистая среда с разными );

двухмерные неоднородные среды, в которых электромагнитные свойства меняются в двух направлениях (например, наклонные пласты или цилиндры, простирающиеся вдоль одного направления и отличающиеся по от вмещающих горных пород);

трехмерные неоднородные среды, в которых свойства меняются по трем направлениям (самой простой из подобных моделей является шар с разными , или в однородном полупространстве).

171

В.К. Хмелевской, В.И. Костицын

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

Существуют различные подходы к решению прямых задач с помощью уравнения (3.5). Любое правильное решение, удовлетворяющее всем физическим требованиям, единственно и корректно. Под корректностью понимается такое решение, в котором малым изменениям исходных данных соответствуют малые приращения расчетных параметров.

3.10. Нормальные поля в электроразведке

Под нормальным полeм понимается электромагнитное поле того или иного источника над однородным изотропным полупространством с неизменными электромагнитными свойствами. Из прямой задачи для поля точечного источника постоянного тока на земной поверхности можно получить нормальные поля постоянных электрических токов для разных установок или комбинаций питающих (АВ) и приемных (МN) электродов.

В электроразведке чаще всего применяются четырехэлектродные установки АМNВ (рис. 3.2). К одному питающему электроду (например А) подключается положительный полюс источника тока, к другому (В) – отрицательный. Разность потенциалов на приемных электродах (МN) от электрода А, определенная по полученной выше формуле (3.1), равна

U A

J

1

1

.

 

 

 

 

 

2

 

AM

 

AN

 

 

 

 

Аналогичным образом можно получить разность потенциалов от отрицательного полюса В, но величину тока следует принять равной (- J).

172

AN BM

Электроразведка

Рис. 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

В.К. Хмелевской, В.И. Костицын

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

Электроразведка

По этим же формулам можно рассчитать некоторое сопротивление над реальным, неизвестным и практически всегда неоднородным полупространством. Тогда оно называется ка-

жущимся сопротивлением (КС или к).

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

3.11. Электрическое поле точечного источника постоянного тока над двухслойной средой

Простейшей прямой задачей для электроразведки ме-

тодом сопротивлений является задача об электрическом поле и кажущемся сопротивлении на поверхности полупространства, представляющего собой верхнее полупространство – воздух, а нижнее – двухслойную горизонтально слоистую среду с мощностью верхнего слоя h1, нижнего h2 = , УЭС слоев 1, 2 и 3 =

(воздух) (рис. 3.3).

Поставленная задача может быть решена с помощью уравнения (3.2), которое при f = 0 превращается в уравнение Лапласа U = 0, где U – потенциал в любой точке М с напряженностью электрического поля E U / r .

Рис. 3.3. Решение прямой задачи о поле точечного источника постоянного тока над двухслойной средой методом зеркальных отражений

175

В.К. Хмелевской, В.И. Костицын

Однако ее можно быстро решить методом зеркальных отражений. Согласно правилам метода зеркальных отражений уравнение Лапласа и физические требования, в том числе граничные условия, выполняются, если потенциал в одномерной среде, где расположен точечный источник, принять равным сумме потенциалов этого источника (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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Электроразведка

 

Подставив в эту формулу производную

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

В.К. Хмелевской, В.И. Костицын

Рис. 3.4. Двухслойная палетка ВЭЗ: 1 и 2 – теоретические и полевая кривые

Одномерные прямые задачи электроразведки для много-

слойных горизонтально слоистых сред при любых первичных полях все-таки сводятся к аналитическим формулам при расчете КС. В результате принято строить кривые КС, аналогичные приведенным на рис. 3.4.

Двухмерные и трехмерные прямые задачи электрораз-

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

3.12. Принципы решения обратных задач электроразведки

Накопленный материал по физическому и математическому моделированию прямых задач электроразведки позволил создать методы решения обратных задач, т.е. определять те или иные параметры геоэлектрического разреза по наблюденным графикам E, H или, например, кривым КС. Решение обратных задач неоднозначно в силу некорректности, как и всех обратных задач математической физики. Некорректность проявляется в

178