УДК 550.3 : 622.02 (075.8)
К. Л. Дудко, А. И. Шиканов
ВЛИЯНИЕ ГЕОЛОГИЧЕСКОГО СТРОЕНИЯ ЭЛЕКТРОПРОФИЛЯ НА КАЖУЩЕЕСЯ ЭЛЕКТРОСОПРОТИВЛЕНИЕ ПОРОД
В геологическом отношении железорудные месторождения Сибири имеют много общего [1]. Углы падения рудных тел составляют, как правило, 70-900. Вмещающая толща сложена крутопадающими массивами скарна, сланцев, известняков, порфирита, сиенита и др. При этом массивы характеризуются дизъюнктивными тектоническими нарушениями, высокой трещиноватостью. Трещинная тектоника ориентирована в основном в меридиональном направлении.
Применяемое на рудниках для прогноза горных ударов электропрофилирование проводится в выработках, ориентированных горизонтально. При проведении электропрофилирования с контура горных выработок линия профиля может пересекать границы различных по составу слоев пород, плоскости нарушений и крупных трещин. Поэтому источник поля, линия наблюдений и поверхности раздела могут занимать различные положения в пространстве.
В работе [2] рассмотрена задача влияния вертикального плоского контакта двух сред с удель-
РМ А
1
42 44 46 48 50 52 54
Рис. 1 Распределение КУЭС над контактом
сланцев и сиенита 1 - модель удельных сопротивлений (р^/й )' 2 - теоретический график КУЭС; 3 - измеренные значения КУЭС; 42-54 - номер пикета
ными электросопротивлениями (УЭС) р и Р2 (р > р2) на параметры поля точечного источника. Получены формулы для расчета кажущегося удельного электросопротивления (КУЭС) вдоль профиля.
Результаты расчетов по этим формулам приведены на рис.1, 2.
Результаты расчетов показывают, что контакт двух сред существенно влияет на распределение КУЭС по профилю измерений. Отражение тока от границ вносит существенный вклад в величины КУЭС, отклоняя их значения от истинного УЭС тем больше, чем больше коэффициент отражения тока.
В работе [2] также рассмотрена модель неоднородной среды, представляющей полупространство из нескольких крутопадающих пластов с плоскими границами.
Если расстояние между границами превышает размеры измерительной установки, то задача сводится к выделению одиночных контактов. Иначе
32 34 36 38
Рис. 2 Распределение КУЭС над контактом сланцев и диорита
1 - модель удельных сопротивлений (Р2 /р ) ;
2 - теоретический график КУЭС; 3 - измеренные значения КУЭС; 32-38 - номер пикета
12
К. Л. Дудко, А. И. Шиканов
необходимо учитывать экранирующее влияние соседних пластов.
Для случая, когда крайние пласты имеют неограниченную мощность, получены следующие выражения для потенциальных функций и) и
градиентов потенциала ди& в каждом слое:
иг (г, г) = } (аге~тг + Ьгет) 10 (тг )с1т
(1)
= -¡(ате~т + Ь,тет) 10(тг )йт д 0
(2)
где ¡о(тг) - функция Бесселя нулевого порядка
первого рода; т - переменная интегрирования, имеющая смысл частоты повторяемости цилин-
дрической функции на единицу длины; а|, Ь -постоянные коэффициенты (I = 1,2,...5) .
В виду неограниченности первого и пятого слоев коэффициенты а1 = 1р/ 4п и Ь5 = 0.
Остальные коэффициенты можно найти с помощью граничных условий, предусматривающих равенство потенциалов и нормальных составляющих плотности тока на границах раздела сред.
Решив систему уравнений для выбранного варианта и определив соответствующие коэффициенты аI, Ь , получаем выражения для вычисления потенциалов и их градиентов во всех пластах. Полученные выражения представляют собой несобственные интегралы, содержащие функции Бесселя: Применительно к геоэлектрике существует несколько алгоритмов вычисления таких интегралов.
При их вычислении получается функциональ-
Р/Р1
Р/Р1 А
/
;
•
1
с к № _ с. и: __
ЕГ и
ч
!
2
V
V
1
И § О
0 25 50 75
XX «к «ч» ЖС XX XX
XX хх XX XX XX
2, М
58 60 62 64
— «
Я" о я
0> — И - — - <и
-3ч и § Ее И -а ц
о — к а — - о
0 25 50 75 100 125 150 175
ъ, м
А А АЛЛ •V» «V»
*ч» *ч# *ч» •ч. >» А А А А *Ч» ААА
Рис.3. Распределение КУЭС над пластом из сланца
1 - модель удельных сопротивлений (р2 /р );
2 - теоретический график КУЭС; 3 - измеренные значения КУЭС; 58-64 - номер пикета
12 14 16 18 20 22 24 26
Рис.4. Распределение КУЭС над сближенными пластами из сиенита и полевого шпата
1 - модель удельных сопротивлений (р> /р );
2 - теоретический график КУЭС; 3 - измеренные значения КУЭС; 12-26 - номер пикета
r =V(х - хс)2 + (У - Ус)2 ; х, y, z - координаты расчетной точки; Хо, yo, zo - координаты источника тока. Для расчета коэффициентов аппроксимации As и Bs определяют корни разностного уравнения через обратную матрицу тёплицевой структуры, формируют полином из полученных коэффициентов, определяют его экспоненциальные корни
p = e Bs и решают систему линейных алгебраических уравнений, сформированную корнями полинома в виде матрицы Вандермонда относительно коэффициентов As. При этом получаемые коэффициенты Bs не знакопеременны (основное условие использования тождества Вебера-Липшица).
Рассмотренный выше алгоритм вычисления значений функции Ф для n - слойной модели среды был реализован на языке JavaScript с использованием программ, приведенных в работе [4]. Результаты расчетов приведены на рис. 3, 4.
Результаты расчетов и сравнение их с экспериментальными данными (рис.1-4) показывает, что геологическое строение электроразведочного профиля существенно отражается на характере распределения КУЭС вдоль профиля. Это необходимо учитывать при обработке результатов измерений для оценки напряженного состояния массива.
СПИСОК ЛИТЕРАТУРЫ
1. Курленя М. В. Геомеханические проблемы разработки железорудных месторождений Сибири / М. В. Курленя, А. А. Еременко, Б. В. Шрепп // Новосибирск, Наука, 2001. - 184 с.
2. Матвеев Б. К. Электроразведка.- Москва, Недра, 1990. - 368 с.
3. Зубов К. Н. Метод расчета заземляющих устройств произвольной конфигурации в неоднородных грунтах / К. Н. Зубков, А. Е. Немировский // Липецк, Вести высш. учеб. заведений Черноземья № 2, 2010. - 21-26 с.
4. Кафтанова Ю.В. Специальные функции математической физики. Часть I. - Харьков, Новое слово, 2009. - 178 с.
Авторы статьи
Дудко Константин Львович,
ассистент каф. теоретической и геотехнической механики КузГТУ e-mail: constantin.dudko@yandex.ru
Шиканов Алексей Иванович
канд.техн.наук, доцент каф. теоретической и геотехнической механики КузГТУ e-mail: shai.tgm@mail.ru
Поступило в редакцию 09.12.2014
ный ряд, коэффициенты которого зависят от геометрических и электромагнитных параметров модели. Эти коэффициенты умножаются на собственные функции, зависящие от величины разноса. Применительно к нашей задаче (небольшие по сравнению с электроразведочными работами разносы электродов) целесообразно использовать метод расчета потенциала точечного источника тока в многослойной плоскопараллельной модели среды, предложенный в работе [3]. Аппроксимацию подынтегральной функции
Ф = ае~тг + Ьет
проводят экспоненциальным многочленом по методу Прони, что позволяет произвести взятие несобственного интеграла, используя тождество Ве-бера-Липшица (аналогичный подход аппроксимации экспонентами - алгоритм Андерсона - используется в электроразведке):
ф = j ф(т, z0, z )10 (mr )dm =
A,
о
<» N N
= í lo(mr) Z Ле- mBsdm = ; ; o s=1 s=1д/ r + B;
(3)
где N - старшая степень аппроксимирующего многочлена, соответствующая количеству слоев исследуемой модели среды;
As, Bs - коэффициенты аппроксимации;