Еще лет 50 назад были начаты исследования с целью обоснования возможности применения бесконтактных измерений в электроразведке методом сопротивлений. Необходимость применения такой модификации возникает при проведении измерений в зонах многолетнемерзлых пород, скальных грунтов, при наличии снежного покрова или такого непроводящего искусственного покрытия, как асфальт или бетон. Кроме того, отсутствие необходимости заземления токовых и измерительных электродов позволяет значительно ускорить процесс измерений. Некоторые результаты таких исследований были опубликованы, например, в работах [4][7-9].

Теоретическое обоснование применения бесконтактных измерений и интерпретации получаемых результатов основывалось не на «строгом» решении соответствующей прямой задачи электродинамики, а на приближенных подходах. И, судя по недавним публикациям наших коллег, положение с теоретическим обоснованием метода до сих пор существенно не изменилось. Например, часть работы [2] посвящена обоснованию тех ограничений, которые имеет применяемая приближенная методика расчетов.

Ниже представлены некоторые, полученные на основе «строгого» решения соответствующей прямой задачи электродинамики, результаты математического моделирования для расположенной на небольшой высоте h над проводящим полупространством предельной дипольно осевой установки BAMN.

Понятно, что проведение бесконтактных измерений в методах сопротивления возможно лишь при переменном токе I в питающей линии AB. Некоторые результаты математического моделирования при переменном токе I = I0 × cos (2π · f · t) в электроразведке методом сопротивлений были приведены в статьях [5][6]. Основой для получения решения задачи послужила методика, описанная в книге [3]. Эта же методика была применена при получении приведенных ниже результатов моделирования.

Начнем с относительно простой модели. На рисунке 1 показана модель, для которой были получены результаты расчетов, приведенные в статье [6] и на рисунке 2 в этой работе. Горизонтальная плоскость S разделяет полупространства V1 и V2. Генераторный (AB) и измерительный (MN) диполи расположены в полупространстве V1 (в воздухе) на оси X, параллельной S, на малой высоте h1 над границей S. Компоненту Ex напряженности электрического поля E, которой пропорционально напряжение ΕMN поля E в предельно короткой измерительной линии MN, определяет несобственный интеграл в смысле главного значения ([1], с. 96—101). Подынтегральная функция содержит функции Бесселя первого рода нулевого и первого порядка. Для вычисления таких интегралов было применено преобразование Эйлера.

Представленные в этой работе результаты моделирования получены при частоте f = 16 кГц. Это «рабочая» частота, которая была выбрана для нескольких видов аппаратуры при бесконтактных измерениях в методах сопротивлений [9]. При наших расчетах для проводящего полупространства была задана диэлектрическая проницаемость ε2 = 10. Хотя при указанном значении f выбор величины ε2 несуществен, так как при любых, реальных для горных пород, значениях ε можно пренебречь влиянием токов смещения в нижнем полупространстве. То есть в проводящем полупространстве электромагнитное поле — квазистационарное:

2π · f · ε0 · ε2 · ρ2 << 1, (1)

где ε0 — электрическая постоянная.

Напротив, при тех же частотах f и при очень высоких значениях ρ1 в верхнем полупространстве (в воздухе) справедливо неравенство:

2π · f · ε0 · ε1 · ρ1 >> 1. (2)

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

На рисунке 2 представлены зависимости кажущихся удельных электрических сопротивлений ρк(Ex) и ρк(|Re Ex|) от разноса x при h = 5 см и при трех значениях удельного электрического сопротивления ρ2. Приведенные на этом рисунке значения ρк определены двумя способами. Величины ρк(Ex) определены «стандартным» способом — по амплитуде |Ex| = √(Re Ex)2 + (Im Ex)2 компоненты Ex, где Re Ex — реактивная составляющая, меняющаяся синфазно либо в противофазе со сторонним током I в отрезке AB, а Im Ex — активная составляющая, отличающаяся по фазе от тока I на π/2. То есть ρк(Ex) = K × |Ex| / I0, где K = π · x3 — коэффициент предельной дипольно-осевой установки, расположенной на поверхности проводящего полупространства. Понятно, что для показанной на рисунке 1 модели на величину |Ex| может оказывать существенное влияние поле зарядов ±e, величины которых не зависят от ρ2. Значения ρк(|Re Ex|) определены по модулю реактивной компоненты Re Ex. Понятно, что на эту компоненту не оказывают влияния заряды ±e на концах отрезка AB.

На рисунке 2 видим, что при малых разносах на величину ρк(Ex) практически не влияют значения ρ2, так как величина |Ex| зависит в основном от зарядов ±e. При бóльших разносах и очень высоких значениях ρ2 (1000 Ом·м) величины ρк(Ex) соответствуют значениям ρ2, но при более низких ρ2 (100, 10 Ом·м) определить ρ2 по ρк(Ex) вряд ли возможно. При дальнейшем увеличении разноса значения ρк(Ex) выходят на известную в методах сопротивлений «индукционную асимптоту».

Совсем по-другому «выглядят» приведенные на рисунке 2 зависимости ρк(|Re Ex|) от разноса x. У них нет не зависящих от ρ2 аномально высоких значений ρк при малых разносах, так как составляющая Re Ex не зависит от поля зарядов ±e на концах отрезка AB. Но при этом зависимости ρк(|Re Ex|) от разноса «осложнены» тем, что с увеличением разноса составляющая Re Ex меняет знак. Поэтому для определения ρ2 по ρк(|Re Ex|) в зависимости от f и ожидаемых пределов изменения ρ2 на основе результатов математического моделирования следует выбрать «оптимальный» разнос. Из показанных на рисунке 2 результатов моделирования следует, что при частоте f = 16 кГц и 10 Ом·м ≤ ρ2 ≤ 1000 Ом·м таким оптимальным является разнос x ≈ 8—10 м.

Перейдем к результатам моделирования, полученным для более сложной и приближающейся к возможным реальным условиям при проведении бесконтактных измерений модели. На рисунке 3 показана трехслойная модель горизонтально-слоистой среды. Верхнее полупространство V1 — воздух, нижнее полупространство V3 — однородные по ρ горные породы, а слой V2 при высоких значениях ρ2 может отвечать асфальту или бетону, а при низких ρ2 — например, влажному почвенному слою. При разработке алгоритма расчетов для этой модели потребовалось получить решение системы 8 уравнений с 8 неизвестными.

На рисунке 4 приведены зависимости ρк(Ex) и ρк(|Re Ex|) от разноса x при высоком значении удельного сопротивления ρ2. Эти зависимости не имеют коренных отличий от тех, которые были приведены на рисунке 2. Присутствие в проводящем полупространстве слоя V2 с высоким ρ2 сказывается на этих зависимостях примерно так, как если бы в модели на рисунке 2 высота h1 была заменена на большее значение h2 (см. рис. 3).

Зависимости ρк(Ex) и ρк(|Re Ex|) от разноса x на рисунке 5 получены при низком значении удельного сопротивления ρ2. В этом случае, особенно при высоком значении ρ3, определение истинного удельного сопротивления ρ3 как по ρк(Ex), так и по ρк(|Re Ex|) может быть затруднено. Это, по-видимому, обусловлено особенностями распределения поверхностных зарядов, индуцированных на границах S1, S1 при ρ1>>ρ2<<ρ3.

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

Установлено, что использование реактивной составляющей поля позволяет исключить это влияние и получить более устойчивую связь между кажущимся и истинным удельным электрическим сопротивлением. Для двухслойной модели определен оптимальный диапазон разносов (около 8—10 м для частоты 16 кГц).

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

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