Все научные статьи
Титов К.В. О влиянии поверхностной проводимости на электропроводность горных пород
Научная статья
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1013а Титов K.B.(tito@kt4096.spb.edu)
Всероссийский НИИ разведочной геофизики им. А.А. Логачева (ВИРГ-Рудгеофизика)
ВВЕДЕНИЕ
В соответствии с эмпирическим законом Арчи отношение удельного электрического сопротивления породы р и равновесного с ней раствора ро является постоянной величиной, не зависящей от концентрации ионов в растворе. Указанное отношение, названное впоследствии параметром пористости F, связано с пористостью породы п и извилистостью путей тока через породу $ (электрической извилистостью) [1]:
г,Р&1
F = Ч = Ч(1)
Роаа п
Такая интерпретация закона Арчи основана на представлении о породе как
о совокупности капилляров сколь угодно сложной формы, которая учитывается с
помощью электрической извилистости.а Посколькуа определениеа последней на
практике затруднительно, используют другую форму записи закона Арчи [1]:
F = a-n"m,аа (2)
где a, m - эмпирические коэффициенты. Как правило, а=1, ш~2.
Оценка пористости породы на основе закона Арчи необходима для
определенияаа гидравлическойаа проницаемости.а Многочисленные
экспериментальные данные показали, что закон Арчи справедлив для пород, не содержащих заметного количества тонкодисперсной (глинистой) фракции [2, 3]. Наличие в породе тонкодисперсного материала ведет к аномальному поведению параметра пористости, выраженному в его зависимости от концентрации ионов (электропроводности) равновесного с породой раствора. Такая аномалия связана с влиянием поверхностной проводимости [4]. Поверхностная проводимость вызвана наличием избытков ионов на границе раздела твердой и жидкой фазы в двойном электрическом слое (ДЭС). Если концентрация ионов в равновесном растворе высока, электропроводность пористой среды обеспечивается главным образома объемной электропроводностьюа раствора,аа иа влияниеа поверхностной
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1014 Иной подход к интерпретации электропроводности породы основан на так называемых законах смеси, когда осуществляется теоретическая оценка электропроводности системы, состоящей, по крайней мере, из двух компонентов с известной электропроводностью [5, 6, 7]. Такое представление согласуется с гранулярной моделью, в соответствии с которой, порода представляется матрицей, в которую помещены включения. В качестве матрицы рассматривается вода с конечной электропроводностью, в роли включений - частицы твердой фазы, электропроводность которой, как правило, принимается бесконечно малой.
Цель настоящей работы состоит в теоретическом исследовании влияния поверхностной проводимости на электропроводность пород на основе капиллярной и гранулярной моделей.
1. ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ
В первую очередь, наше внимание привлек факт высокой электропроводности кимберлитов - пород, содержащих промышленные скопления алмазов, отмеченный практически повсеместно при наземных и скважинных электроразведочных наблюдениях [8]. Кимберлит практически не содержит минералов с электронной проводимостью и, поэтому, проводит электрический ток по ионному механизму. Его удельное электрическое сопротивление составляет около 10 Ом'м, что существенно меньше, чем для большинства других пород, насыщенных пресными водами. На рис.1 показано распределение удельного электрического сопротивления в верхней части кимберлитовой трубки "Архангельская" месторождения им. М.В. Ломоносова (Архангельская область) и в перекрывающих ее слабосцементированных песчаниках. Минимальное значение удельного сопротивления кимберлитов составляет 7.6 Ом'м, максимальное - 23 Ом'м, среднее геометрическое - 11 Ом'м. Как видно из рис.1, удельное электрическое сопротивление кимберлита существенно ниже, чем у перекрывающих его песчаников.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 1015а
0 20 40 60 80 100 120аа 140аа 160 100 100 140 73 150 150 150 240 86 190 86 350 23 490 240 200 240 240 10- 260 20- 15 15 30- 8.0 7.6 40- 7.6 7.6 11 11
глинистый песчаник
Рис. 1. Геоэлектрический разрез верхней части кимберлитовой трубки "Архангельская" и перекрывающих ее пород. Цифрами показаны значения удельного электрического сопротивления пород по данным вертикального электрического зондирования, Ом'м.
Для изучения природы аномалии были выполнены лабораторные измерения зависимости параметра пористости кимберлита от сопротивления равновесного с ним раствора. Порошок кимберлита привели в равновесие с растворами КС1 различной концентрации. Удельное электрическое сопротивление измеряли с помощью стандартной четырехэлектродной установки, откалиброванной на специально приготовленных растворах КС1 с известной электропроводностью. Результаты измерений показывают существенную зависимость параметра пористости от удельного сопротивления раствора (Рис.2). Были выполнены два цикла измерений, когда концентрация раствора увеличивалась и уменьшалась. Оба цикла дали хорошее совпадение точек на графике. В области низкого сопротивления равновесного раствора (доминирует объемная проводимость) параметр пористости составил 3.5, что отвечает пористости 53% (при значении показателя степени в законе Арчи ш=2, формула (2)) и соответствует пористости порошка кимберлита. В области высокого сопротивления равновесного раствора (доминирует поверхностная проводимость) параметр пористости достигает величины 0.18, т.е. электрическое сопротивление пористой системы оказывается много ниже сопротивления равновесного с ней раствора. Этот эффект известен в коллоидной химии как коллоидная "сверхпроводимость" [4].
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1016а
100.00
Область коллоидной "сверхпроводимости"
Рис.а 2.а Зависимость параметра пористости кимберлитового порошка от
удельного сопротивления равновесного раствора. Кресты обозначают результаты
измерения при уменьшении концентрации раствора, черные треугольники - с
увеличением концентрации. При значении сопротивления раствора более 2 Ом м
отмечается коллоидная "сверхпроводимость."
2. ТЕОРЕТИЧЕСКИЕ ПРЕДСТАВЛЕНИЯ
Поскольку поверхностная проводимость является свойством межфазной границы, при теоретическом рассмотрении она может быть отнесена как к жидкой, так и к твердой фазе. Рассмотрим последовательно обе возможности.
2.1 ИСПОЛЬЗОВАНИЕ КАПИЛЛЯРНОЙ МОДЕЛИ ПОРОДЫ И ОСРЕДНЕНИЯ ЭЛЕКТРОПРОВОДНОСТИ РАСТВОРА ПО СЕЧЕНИЮ ПОР
Для исследования пород, содержащих тонкодисперсный материал, используют так называемый "истинный" параметр пористости, свободный от влияния поверхностной проводимости [1]. "Истинный" параметр пористости определяется при высокой электропроводности равновесного с породой раствора,
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1017 F0= Hm Ч, (3)
Ро^ Ро где р, ро - удельное электрическое сопротивление породы и равновесного с ней раствора соответственно.
В.В. Кормильцев [9] предложил для теоретических оценок использовать относительную электропроводность раствора в поре:
<т
Х=Ч (4)
^о
где <7 - электропроводность раствора, усредненная по сечению поры и по всем порам, Оо=1/ро - электропроводность равновесного раствора.
Параметр пористости, "истинный" параметр пористости и относительную электропроводность можно связать простым соотношением:
F(<70) = -^-.аа (5)
В формуле (5) разделено влияние пористости и извилистости (характеризуемое "истинным" параметром пористости) и поверхностной проводимости (характеризуемой относительной электропроводностью). В соответствии с ней может быть определена зависимость относительной электропроводности от электропроводности равновесного раствора %(оо) по экспериментальным данным (по отношению "истинного" параметра пористости и параметра пористости).
Для теоретической оценки относительной электропроводности необходимо сделать предположение о форме пор. В.В. Кормильцев [9] предположил, что порода пронизана цилиндрическими порами, и сделал соответствующие оценки, решив численно сопряженную задачу о потенциале и концентрации в ДЭС. Мы предлагаем использовать более простую оценку для плоского ДЭС поскольку, во-первых, для относительно грубодисперсных сред кривизна зерен существенно менее чем дебаевский радиус, характеризующий толщину ДЭС; во-вторых, частицы глинистых минералов имеют форму пластин [10] и для описания пор между ними скорее подходит геометрия щели, чем цилиндра.
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1018а %=dЧiЧ+1,а (6)
о
где у = In
'аа (Г ^ Vаа v 2 J)
яа lа qRTаа *
о=Чаа -------- - дебаевскии радиус,
zFy2C0 2d - ширина поры,
zFcpi
Г-у =--------- безразмерный потенциал на внутренней границе диффузного слоя,
2RT
ф1 (В) - электрический потенциал на той же границе,
z - валентность ионов,
Р=9.65 104(Кл/моль) - постоянная Фарадея,
s=80 ' 8.85 ' 10" (Ф/м) - абсолютная диэлектрическая проницаемость воды,
R=8.31 (Джа (мольа К)" ) - универсальная газовая постоянная,
Т (К) - абсолютная температура,
Соа (моль/м ) - равновесная концентрация ионов ва свободном растворе (т.е.,
например, концентрация на осевой плоскости широкой поры, где не сказывается
влияние ДЭС).
Поаа экспериментальнымаа даннымаа ваа соответствииаа саа формулойаа (5)аа мы
определилиа относительнуюа электропроводностьа кимберлитаа иа сопоставилиа с
теоретическимиаа оценкамиаа поаа формулеаа (6)аа (Рис.аа 3).аа Здесьаа жеаа показаны
экспериментальные данные, полученные на гидрослюдяной синей глине
кембрийского возраста (Ленинградская область), и, полученные американскими
исследователями на порошке гранита Вестерли [11].
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1019а
1000 100 Ч
X Е
ж Кимберлит
^ Синяя глина
? Граниты Вестерли
10"2 10а 102а d/8
Рис. 3 Сопоставление экспериментальных и теоретических данных о поверхностной проводимости пород. Шифры теоретических кривых отвечают значениям модуля межфазного потенциала (pi (мВ).
Экспериментальные данные для гранитов и глин хорошо соответствуют теоретическим зависимостям, отвечающим невысоким значениям межфазного потенциала, что согласуется с известными значениями ^-потенциала. Для кимберлитов не удается добиться хорошего соответствия экспериментальных и теоретических зависимостей.
Для объяснения этого явления мы обратились к экспериментальным данным Д.А. Фридрихсберга и Ю.П. Шишкина [12]. Безразмерный межфазный потенциал Tiаа можноаа представитьаа кака разностьаа полногоаа межфазногоаа потенциалаа Гоаа и приращения потенциала на слое Гельмгольца Афоь zF
A#V |
(7) |
Г =Г
2RT |
L1аа Lо
Ваа экспериментеаа н керамическойаа диафрагмеаа (фо=-200аа мВ),аа саа ростом концентрацииаа КС1аа отаа 510"3аа Nаа доаа 50 10"3аа N,аа значенияаа -cpiаа закономерно
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1020а аппроксимировать зависимостью:
^1=^0
/
1
(8)
где Со выражена в моль/литр.
Указанная зависимость может быть обобщена:
(9)
-1/2ч
где A, b - постоянные, определяемые из опыта, А=0.5 ... 1 (моль"аа ), Ь=0.5...1. Для безразмерных потенциалов формула (8) принимает вид:
Г,=
z<p0
50
1-
= ГА
1-
(10)
В соответствии с формулой (10) потенциалы cpi и Ti состоят из постоянной части, определяющейся фо (или Го) и переменной (поправки), зависящей от концентрации раствора.
Предложенная зависимость (10) является аппроксимацией данных в относительно узком диапазоне концентраций, она не применима для всего диапазона минерализации природных вод. Так, при Со> 4 моль/литр, в соответствии с формулой (10) Г\ меняет знак, что противоречит исходному положению о том, что Со - концентрация противоионов. Можно считать, что формула (10) справедлива только для разбавленных растворов.
На основе формулы (10) можно связать потенциал Г\ с электропроводностью свободного раствора:
|
Ti=r0
(11)
где k=Fu3z, иэ - электролитическая подвижность ионов.
Дляа разбавленныха растворова КС1а к= 1.3 10" См'м (моль)" ,аа поэтому,аа в частности, для КС1:
(12) |
r1=r0(l-4.4V^).
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1021а 3 2.5
ln(<Po-<Pi)аа 2
1.5
1
0.5
0
1аа 10а 100
С0, mN
Рис. 4. Зависимость разности потенциалов на границе слоя Гельмгольца от концентрации в миллинормальном диапазоне. Точки отвечают экспериментальным данным Д.А.Фридрихсберга и Ю.П.Шишкина, сплошная линия - аппроксимация по формуле (10).
С учетом предложенной поправки мы вновь вычислили зависимость относительной электропроводности от d/8 и сопоставили теоретические кривые с экспериментальными зависимостями для кимберлита (Рис. 5). Экспериментальные данные хорошо соответствуют теоретической зависимости, отвечающей | фо | =300 мВ. Полученный результат позволяет предположить, что кимберлит характеризуется более высоким межфазным потенциалом по сравнению с другими породами, что объясняет его высокую электропроводность.
По соответствию экспериментальных и теоретических зависимостей можно оценить значения межфазного потенциала (pi и раскрытие пор 2d (см. формулу (6)). Оценка характерных размеров пор для глины, кимберлита и гранита приводит к значениям соответственно 10" , 10" , 3 10" м, что не противоречит имеющимся экспериментальным данным для глин (10" [1, 10]), и почти в десять раз превышает оценку для гранитов Вестерли (210" [13]). Таким образом, данные о зависимости относительной электропроводности породы от электропроводности равновесного раствора позволяют оценить межфазный потенциал и характерную ширину пор, т.е. степень дисперсности породы. Оценка дисперсности, вероятно, менее надежна, чем межфазного потенциала.
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1022а 1 10
Ро, Омм
Рис. 5. Теоретические зависимости х(Ро) для среды со средним размером пор 0.1 мкм и экспериментальная зависимость х(Ро) для кимберлитового порошка. Цифрами даны значения модуля потенциала фо, треугольниками показаны данные, отвечающие увеличению концентрации раствора, крестами -уменьшению.
2.2 ИСПОЛЬЗОВАНИЕ ГРАНУЛЯРНОЙ МОДЕЛИ ПОРОДЫ И ФОРМУЛЫ СМЕСИ
Для вычисления эффективной электропроводности (эффективной
диэлектрической проницаемости) двухкомпонентной смеси в разное время был
предложен ряд формул. Как отмечается в работе [5], формула Бруггемана,
модифицированнаяаа Ханаи,аа наилучшим образом удовлетворяет
экспериментальным данным в широком диапазоне концентрации включений. В петрофизической практике эта формула использована впервые Сеном и др. [6] и носит название формулы Бруггемана-Ханаи-Сена (БХС):
<у - а*.
<j - <тп
^0_
,d
= п
(13)
где п - доля воды в двухкомпонентной системе (пористость), d - коэффициент деполяризации частиц твердой фазы, принимающий значения от 0 до 1. В частности, для длинных цилиндрических частиц, параллельных направлению электрического поля d=0, для сферических частиц d=l/3, для тонких пластин, перпендикулярных полю d=l. Коэффициент деполяризации связан с показателем экспоненты в законе Арчи:
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1023а \- а
Для учета поверхностной проводимости воспользуемся тем фактом, что
поле двухслойной сферической частицы эквивалентно полю однородной частицы
саа некоторойаа эффективнойаа электропроводностью.аа Еслиаа электропроводность
внутренней сферы равна нулю, а второй слой - тонкий с проводимостью S, то
эффективная электропроводность двухслойной сферической частицы с внешним
радиусом а составляет:
,аа S
атэ=Ч. (15)
la
Этот подход при моделировании электропроводности был использован в работах [3, 7]. Поверхностную проводимость можно вычислить по формуле Бикермана [14] в пренебрежении электроосмосом (эффект второго порядка малости по сравнению с ошибками, связанными с заменой природных полиионных электролитов модельным электролитом, содержащим по одному аниону и катиону валентности равной единице):
S = 4SzFu3C0(chTl-\)(16)
Подстановка результатов вычислений по формуле (16) в выражение (15) позволяет вычислить эффективную электропроводность частиц твердой фазы с учетом поверхностной проводимости. Дальнейшая подстановка эквивалентной электропроводности твердой фазы в формулу БХС (13) позволяет вычислить зависимости параметра пористости от пористости породы (Рис. 6).
В гидрогеологической и нефтепромысловой практике удельное сопротивление пород используют для оценки пористости и коэффициента фильтрации. Графики 1,2 на рис. 6 иллюстрируют возможность оценки коэффициента фильтрации неглинистых пород по электрическому сопротивлению. Уменьшение сопротивления с ростом пористости при слабом влиянии поверхностной проводимости сопровождается ростом коэффициента фильтрации. Эта зависимость характерна, например, для скальных пород. Напротив, увеличение сопротивления с ростом пористости характерно при доминирующем влиянии поверхностной проводимости (кривые 5, 6). Эта зависимость характерна для суглинков и глин. Кривая с индексом 4, совпадающая
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1024
1.Е+03 Рп
1.Е+02 1.Е+01
1.Е+00
1.Е-01
1.Е-02
0.001аа 0.01аа 0.1аа 1
п
Рис. 6. Зависимости параметра пористости Рп=р/ро от пористости для двухкомпонентной смеси при разных значениях концентрации свободного раствора Со (по формуле БХС) с учетом поверхностной проводимости. Концентрация раствора С0: 1 - 1.7 N, 2 - l.TlO^N, 3 - 1.710"2N, 4 - 8.4 Ю-3 N, 5 -1.7'10"3N, 6 - 1.7'10"4N. Потенциал ф0=150 мВ, диаметр частиц - 10"6м, электролитическая подвижность ионов иэ= 4 10" м (В с)" , показатель степени d=l/2 ( что соответствует показателю цементации в законе Арчи ш=2).
Изопроводимость возникает при равенствеаа эффективной
электропроводности включений <тт и электропроводности свободного раствора.
Для принятых при вычислениях значений размера частиц, межфазного потенциала и подвижности ионов состояние изопроводимости соответствует концентрации Со=8.4 mN (480 мг/л NaCl), достаточно характерной для природных вод. Приведенные результаты согласуются с результатами вычислений, приведенными в работе [4].
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1025а 1.а При теоретическом описании и интерпретации электропроводности пород
поверхностная проводимость может быть отнесена как к жидкой, так и к твердой
фазе.
2. Высокая электропроводность кимберлитов может быть объяснена
влиянием поверхностной проводимости, связанной с аномально высоким
межфазным потенциалом.
3.аа Влияние поверхностной проводимости пород приводит к следующим
петрофизическим следствиям:
- зависимость параметра пористости от пористости и проницаемости имеет прямой или обратный характер при разном соотношении объемной проводимости раствора и поверхностной проводимости;
- при некотором соотношении значений концентрации порового раствора, межфазного потенциала, дисперсности параметр пористости оказывается независимым от пористости.
БЛАГОДАРНОСТИ
Автор признателен профессору Санкт-Петербургского университета В.А. Комарову за критику и обсуждение положений настоящей работы, А.Н. Слюсареву - за участие в лабораторных экспериментах, Ф.Г. Атрощенко - за предоставленную возможность провести полевые наблюдения на месторождении им. Ломоносова.
ЛИТЕРАТУРА
- Добрынин В.М., Вендельштейн Б.Ю. и Кожевников Д.А. Петрофизика. - М.: Недра,-1991.-368 С.
- Keehm Y., Mukerji Т., Nur A. Computational rock physics at the pore scale: Transport properties and diagenesis in realistic pore geometries // The leading edge. -2001.-V. 20. -C. 180- 183.
- Revil A., Cathles L.M., Losh S. Electrical conductivity in shaly sands with geophysical applications // Journ. Geoph. Res. -1998. -V.103B. -C. 23925-23936.
- Фридрихсберг Д.А. Курс коллоидной химии. - Л.:Химия, -1974. - 352 С.
- Челидзе Т.Л., Деревянко А.И., Курилленко О.Д. Электрическая спектроскопия гетерогенных систем. -Киев: Наукова думка, 1977. -231 С.
Электронный журнал ИССЛЕДОВАНО В РОССИИа 1026а
Все научные статьи