О влиянии поверхностной проводимости на электропроводность горных пород

Титов К.В.(Шо@кт.4096.8рЬ.е(ш)

Всероссийский НИИ разведочной геофизики им. А.А. Логачева (ВИРГ-Рудгеофизика)

ВВЕДЕНИЕ

В соответствии с эмпирическим законом Арчи отношение удельного электрического сопротивления породы р и равновесного с ней раствора р0 является постоянной величиной, не зависящей от концентрации ионов в растворе. Указанное отношение, названное впоследствии параметром пористости F, связано с пористостью породы n и извилистостью путей тока через породу $ (электрической извилистостью) [1]:

F = - = - (1)

р0 n

Такая интерпретация закона Арчи основана на представлении о породе как о совокупности капилляров сколь угодно сложной формы, которая учитывается с помощью электрической извилистости. Поскольку определение последней на практике затруднительно, используют другую форму записи закона Арчи [1]:

F = a n-m, (2)

где а, m - эмпирические коэффициенты. Как правило, a=1, m~2.

Оценка пористости породы на основе закона Арчи необходима для определения гидравлической проницаемости. Многочисленные

экспериментальные данные показали, что закон Арчи справедлив для пород, не содержащих заметного количества тонкодисперсной (глинистой) фракции [2, 3]. Наличие в породе тонкодисперсного материала ведет к аномальному поведению параметра пористости, выраженному в его зависимости от концентрации ионов (электропроводности) равновесного с породой раствора. Такая аномалия связана с влиянием поверхностной проводимости [4]. Поверхностная проводимость вызвана наличием избытков ионов на границе раздела твердой и жидкой фазы в двойном электрическом слое (ДЭС). Если концентрация ионов в равновесном растворе высока, электропроводность пористой среды обеспечивается главным образом объемной электропроводностью раствора, и влияние поверхностной

проводимости на общую электропроводность - невелико. При низкой концентрации ионов объемная проводимость раствора невелика и роль поверхностной проводимости становится существенной. Эта концепция привела к представлению о двойной воде , имеющейся в горных породах. Одним видом воды считается та, которая находится в тонких порах глины, а другим - вода в относительно широких порах между зернами песка [1].

Иной подход к интерпретации электропроводности породы основан на так называемых законах смеси, когда осуществляется теоретическая оценка электропроводности системы, состоящей, по крайней мере, из двух компонентов с известной электропроводностью [5, 6, 7]. Такое представление согласуется с гранулярной моделью, в соответствии с которой, порода представляется матрицей, в которую помещены включения. В качестве матрицы рассматривается вода с конечной электропроводностью, в роли включений - частицы твердой фазы, электропроводность которой, как правило, принимается бесконечно малой.

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

1. ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ

В первую очередь, наше внимание привлек факт высокой электропроводности кимберлитов - пород, содержащих промышленные скопления алмазов, отмеченный практически повсеместно при наземных и скважинных электроразведочных наблюдениях [8]. Кимберлит практически не содержит минералов с электронной проводимостью и, поэтому, проводит электрический ток по ионному механизму. Его удельное электрическое сопротивление составляет около 10 Ом.м, что существенно меньше, чем для большинства других пород, насыщенных пресными водами. На рис.1 показано распределение удельного электрического сопротивления в верхней части кимберлитовой трубки Архангельская месторождения им. М.В. Ломоносова (Архангельская область) и в перекрывающих ее слабосцементированных песчаниках. Минимальное значение удельного сопротивления кимберлитов составляет 7.6 Омм, максимальное - 23 Омм, среднее геометрическое - 11 Омм. Как видно из рис.1, удельное электрическое сопротивление кимберлита существенно ниже, чем у перекрывающих его песчаников.

глинистый песчаник

260 песчаник 15 кимберлит

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

Для изучения природы аномалии были выполнены лабораторные измерения зависимости параметра пористости кимберлита от сопротивления равновесного с ним раствора. Порошок кимберлита привели в равновесие с растворами KCl различной концентрации. Удельное электрическое сопротивление измеряли с помощью стандартной четырехэлектродной установки, откалиброванной на специально приготовленных растворах KCl с известной электропроводностью. Результаты измерений показывают существенную зависимость параметра пористости от удельного сопротивления раствора (Рис.2). Были выполнены два цикла измерений, когда концентрация раствора увеличивалась и уменьшалась. Оба цикла дали хорошее совпадение точек на графике. В области низкого сопротивления равновесного раствора (доминирует объемная проводимость) параметр пористости составил 3.5, что отвечает пористости 53% (при значении показателя степени в законе Арчи m=2, формула (2)) и соответствует пористости порошка кимберлита. В области высокого сопротивления равновесного раствора (доминирует поверхностная проводимость) параметр пористости достигает величины 0.18, т.е. электрическое сопротивление пористой системы оказывается много ниже сопротивления равновесного с ней раствора. Этот эффект известен в коллоидной химии как коллоидная сверхпроводимость [4].

5.00 -t

4.00

з.оо -4

2.00

1.00

0.00

100.00

Область коллоидной сверхпроводимости

Рис. 2. Зависимость параметра пористости кимберлитового порошка от удельного сопротивления равновесного раствора. Кресты обозначают результаты измерения при уменьшении концентрации раствора, черные треугольники - с увеличением концентрации. При значении сопротивления раствора более 2 Омм отмечается коллоидная сверхпроводимость.

2. ТЕОРЕТИЧЕСКИЕ ПРЕДСТАВЛЕНИЯ

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

2.1 ИСПОЛЬЗОВАНИЕ КАПИЛЛЯРНОЙ МОДЕЛИ ПОРОДЫ И ОСРЕДНЕНИЯ ЭЛЕКТРОПРОВОДНОСТИ РАСТВОРА ПО СЕЧЕНИЮ ПОР

Для исследования пород, содержащих тонкодисперсный материал, используют так называемый истинный параметр пористости, свободный от влияния поверхностной проводимости [1]. Истинный параметр пористости определяется при высокой электропроводности равновесного с породой раствора,

когда влияние поверхностной проводимости пренебрежимо мало по сравнению с объёмной:

F0 = lim -, (3)

Ро ->0 р0

где р, р0 - удельное электрическое сопротивление породы и равновесного с ней раствора соответственно.

В. В. Кормильцев [9] предложил для теоретических оценок использовать относительную электропроводность раствора в поре: (7

Х=- (4)

где 7 - электропроводность раствора, усредненная по сечению поры и по всем порам, (70=1/р0 - электропроводность равновесного раствора.

Параметр пористости, истинный параметр пористости и относительную электропроводность можно связать простым соотношением:

F (70) = - 7 . (5)

В формуле (5) разделено влияние пористости и извилистости (характеризуемое истинным параметром пористости) и поверхностной проводимости (характеризуемой относительной электропроводностью). В соответствии с ней может быть определена зависимость относительной электропроводности от электропроводности равновесного раствора х(ао) по экспериментальным данным (по отношению истинного параметра пористости и параметра пористости).

Для теоретической оценки относительной электропроводности необходимо сделать предположение о форме пор. В. В. Кормильцев [9] предположил, что порода пронизана цилиндрическими порами, и сделал соответствующие оценки, решив численно сопряженную задачу о потенциале и концентрации в ДЭС. Мы предлагаем использовать более простую оценку для плоского ДЭС поскольку, во-первых, для относительно грубодисперсных сред кривизна зерен существенно менее чем дебаевский радиус, характеризующий толщину ДЭС; во-вторых, частицы глинистых минералов имеют форму пластин [10] и для описания пор между ними скорее подходит геометрия щели, чем цилиндра.

Для случая плоской межфазной границы относительная электропроводность получается путем последовательного определения микроскопического электрического потенциала в ДЭС и концентрации ионов, а затем, - усреднения концентрации по ширине щелевой поры. Решение получается аналитически:

о sh(-)

Х = 2 ~---+1, (6)

sh( y - -) sh( У) о

где y = ln

дебаевский радиус,

2d - ширина поры, zF(p1

Г1 =-- безразмерный потенциал на внутренней границе диффузного слоя,

ф1 (В) - электрический потенциал на той же границе, z - валентность ионов,

F=9.65 104 (Кл/моль) - постоянная Фарадея,

s=80 . 8.85 . 10-12 (Ф/м) - абсолютная диэлектрическая проницаемость воды, R=8.31 (Дж . (моль . K)-1) - универсальная газовая постоянная, T (K) - абсолютная температура,

С0 (моль/м3) - равновесная концентрация ионов в свободном растворе (т.е., например, концентрация на осевой плоскости широкой поры, где не сказывается влияние ДЭС).

По экспериментальным данным в соответствии с формулой (5) мы определили относительную электропроводность кимберлита и сопоставили с теоретическими оценками по формуле (6) (Рис. 3). Здесь же показаны экспериментальные данные, полученные на гидрослюдяной синей глине кембрийского возраста (Ленинградская область), и, полученные американскими исследователями на порошке гранита Вестерли [11].

102 100 102 d/8

Рис. 3 Сопоставление экспериментальных и теоретических данных о поверхностной проводимости пород. Шифры теоретических кривых отвечают значениям модуля межфазного потенциала ф1 (мВ).

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

Для объяснения этого явления мы обратились к экспериментальным данным Д.А. Фридрихсберга и Ю.П. Шишкина [12]. Безразмерный межфазный потенциал Г1 можно представить как разность полного межфазного потенциала Г0 и приращения потенциала на слое Гельмгольца Аф01:

Г1 =Г0 - Аф01. (7)

В эксперименте на керамической диафрагме (ф0=-200 мВ), с ростом концентрации KCl от 510-3 N до 50 10-3 N, значения -ф1 закономерно

уменьшались от 195 до 170 мВ. Опытные данные (Рис. 4) можно аппроксимировать зависимостью:

ф1 =ф0

1 VC0 2

где С0 выражена в моль/литр.

Указанная зависимость может быть обобщена:

ф1 =(0(1 AC 0), (9)

где A, b - постоянные, определяемые из опыта, A=0.5 ... 1 (моль ), b=0.5...1. Для безразмерных потенциалов формула (8) принимает вид:

1 Ж

1 Ж

(10)

В соответствии с формулой (10) потенциалы ф1 и Г1 состоят из постоянной части, определяющейся ф0 (или Г0) и переменной (поправки), зависящей от концентрации раствора.

Предложенная зависимость (10) является аппроксимацией данных в относительно узком диапазоне концентраций, она не применима для всего диапазона минерализации природных вод. Так, при С0> 4 моль/литр, в соответствии с формулой (10) Г1 меняет знак, что противоречит исходному положению о том, что С0 - концентрация противоионов. Можно считать, что формула (10) справедлива только для разбавленных растворов.

На основе формулы (10) можно связать потенциал Г1 с электропроводностью свободного раствора:

1 2

(11)

где k=Fu;3z, u - электролитическая подвижность ионов.

Для разбавленных растворов KCl k=1.3 10-2Смм2 (моль)-1, поэтому, в частности, для KCl:

Ц (1 4.4).

(12)

4 -3.5

3 2.5

1п(фо-ф1) 2 -1.5

1 -0.5 0

Рис. 4. Зависимость разности потенциалов на границе слоя Гельмгольца от концентрации в миллинормальном диапазоне. Точки отвечают экспериментальным данным Д.А.Фридрихсберга и Ю.П.Шишкина, сплошная линия - аппроксимация по формуле (10).

С учетом предложенной поправки мы вновь вычислили зависимость относительной электропроводности от d/8 и сопоставили теоретические кривые с экспериментальными зависимостями для кимберлита (Рис. 5). Экспериментальные данные хорошо соответствуют теоретической зависимости, отвечающей ф0 =300 мВ. Полученный результат позволяет предположить, что кимберлит характеризуется более высоким межфазным потенциалом по сравнению с другими породами, что объясняет его высокую электропроводность.

По соответствию экспериментальных и теоретических зависимостей можно оценить значения межфазного потенциала ф1 и раскрытие пор 2d (см. формулу (6)). Оценка характерных размеров пор для глины, кимберлита и гранита

8 7 6

приводит к значениям соответственно 10 , 10 , 3 10 м, что не противоречит имеющимся экспериментальным данным для глин (10-8 [1, 10]), и почти в десять раз превышает оценку для гранитов Вестерли (210-7 [13]). Таким образом, данные о зависимости относительной электропроводности породы от электропроводности равновесного раствора позволяют оценить межфазный потенциал и характерную ширину пор, т. е. степень дисперсности породы. Оценка дисперсности, вероятно, менее надежна, чем межфазного потенциала.

100-

0.01

-гтттттт,-i .......

0.1 1 10

р0, Омм

50 100

Рис. 5. Теоретические зависимости х(р0) для среды со средним размером пор 0.1 мкм и экспериментальная зависимость х(р0) для кимберлитового порошка. Цифрами даны значения модуля потенциала ф0, треугольниками показаны данные, отвечающие увеличению концентрации раствора, крестами -уменьшению.

2.2 ИСПОЛЬЗОВАНИЕ ГРАНУЛЯРНОЙ МОДЕЛИ ПОРОДЫ И ФОРМУЛЫ СМЕСИ

Для вычисления эффективной электропроводности (эффективной диэлектрической проницаемости) двухкомпонентной смеси в разное время был предложен ряд формул. Как отмечается в работе [5], формула Бруггемана, модифицированная Ханаи, наилучшим образом удовлетворяет экспериментальным данным в широком диапазоне концентрации включений. В петрофизической практике эта формула использована впервые Сеном и др.[6] и носит название формулы Бруггемана-Ханаи-Сена (БХС):

т

сс0 2 с /

(13)

с с т

где n - доля воды в двухкомпонентной системе (пористость), d - коэффициент деполяризации частиц твердой фазы, принимающий значения от 0 до 1. В частности, для длинных цилиндрических частиц, параллельных направлению электрического поля d=0, для сферических частиц d=1/3, для тонких пластин, перпендикулярных полю d=1. Коэффициент деполяризации связан с показателем экспоненты в законе Арчи: