Влияние концентрации солей на вязкость воды при течении в пористой среде

УДК 532.546

Влияние концентрации солей на вязкость воды при течении в пористой среде

В. В. Кадет, А. С. Корюзлов

Кафедра подземной гидромеханики РГУ нефти и газа им. И. М. Губкина Ленинский пр-т, 65, Москва, Россия, 119991

В данной работе представлена перколяционная модель электрокинетического течения раствора электролита в пористой среде. Проанализировано влияние ионной концентрации раствора, дзета-потенциала поверхности порового пространства и вида порометри-ческой кривой на скорость течения. Спланирован и поставлен эксперимент, проведено сравнение экспериментальных данных с результатами численного расчёта.

Ключевые слова: пористая среда, течение электролита в пористой среде, теория пер-коляции, двойной электрический слой, эффект электровязкости.

1. Введение

В ряде экспериментов по исследованию течения флюидов в пористых средах отмечалось превышение эффективной, т.е. наблюдаемой в эксперименте, вязкости по сравнению с её классическим значением для капельной жидкости. В экспериментальных работах [1-3] убедительно показано, что в тонких каналах наблюдается превышение гидравлического сопротивления по сравнению с величинами гидравлического сопротивления, предсказываемыми классической теорией, причем течение в микроканалах зависит от природы жидкости и геометрических параметров канала.

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

2. Модель процесса

Воспользуемся подходом, развитым в [4], который позволяет получать макроскопические законы течения в пористой среде на базе детального анализа течения флюида в отдельном элементе такой среды — поровом канале. В рамках этого подхода в качестве модели порового пространства используется регулярная пространственная решётка, узлами которой являются поры, а связями — капиллярные каналы между ними. Предполагается, что капиллярные каналы имеют круговую цилиндрическую форму. Характеристиками такой решётки служат длина ребра (связи) I, функция плотности распределения поровых каналов по радиусам (порометрическая кривая) f (г) и число рёбер, выходящих из одного узла (координационное число решётки) в.

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

Статья поступила в редакцию 7 декабря 2007 г.

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

В качестве первого шага рассмотрим течение по отдельному каналу радиуса а. Когда жидкость течёт за счёт градиента давления через микроканал, ионы подвижной части ДЭС увлекаются жидкостью, вызывая электрический ток, называемый током протекания. Накопление ионов в области, определяемой направлением течения, является причиной возникновения электрического потенциала, который принято называть потенциалом протекания. В свою очередь потенциал протекания способствует появлению тока ионов в направлении, обратном направлению течения жидкости (ток проводимости).

Массовая сила, действующая в такой системе, будет определяться взаимодействием потенциала протекания ф с зарядом ионного облака объёмной плотностью ре. Таким образом, установившееся течение в капилляре будет описываться уравнением Навье-Стокса следующего вида

1 ё / ёи\ „ „ ,

- Реуф = 0'

где и — скорость течения, р — динамическая вязкость капельной жидкости, Ур — градиент давления, Уф — градиент потенциала протекания вдоль оси канала. Объёмная плотность заряда ре определяется при этом уравнением Пуассона [5,6]

_ _

Ре ££° г ёг \ ёг ) '

где £о — диэлектрическая проницаемость вакуума, а Ф — потенциал ДЭС. Подставив (2) в (1) и проинтегрировав со следующими граничными условиями: и(а) = 0, ф(а) = £ (условие прилипания и условие равенства дзета-потенциала на стенке канала величине £), и и'(0) = ^'(0) = 0 (условия симметричности решения относительно оси канала), получим выражение для скорости течения в канале с учётом влияния потенциала ДЭС

и =4^

Отсутствие адсорбционного слоя в модели (потенциал на стенке равен £) связано с особенностями межфазных процессов на границе кварц-раствор электролита, при фильтрации воды в песчаниках. Наиболее часто в пластовых условиях кислотность лежит в пределах рН=3-8. Согласно исследованиям, проведённым в работах [7, 8], в этом диапазоне кислотности на поверхности контакта кварца (составляющего большую часть поверхности пор песчаника) и водного раствора соли образуется достаточно тонкий адсорбционный слой Штерна, поэтому в данной работе используется модель ДЭС Гуи-Чэпмена.

Как видно из (3), скорость зависит от потенциала протекания, для определения которого запишем закон сохранения заряда. Объёмный ток имеет аксиальное направление и складывается из двух частей: конвекционного тока (тока протекания) и тока проводимости, при этом его плотность есть

] = ге

где т — подвижность ионов, г — их валентность, е — заряд электрона, п+ и п— — объёмные плотности катионов и анионов соответственно.

Когда ток проводимости становится равным току протекания, устанавливается равновесное состояние (установившийся режим течения). Локальная концентрация ионов подчиняется закону распределения Больцмана [5] п± = по ехр^ге^/квТ), где по — суммарная объёмная концентрация ионов. Тогда (4) может быть записано в виде

] = Реи - ауф,

где проводимость жидкости а = 2шгепо сЬ(геф/квТ), а

ре = ге

Проинтегрировав уравнения (3) и (5) по области поперечного сечения канала, получим выражения для объёмного расхода д и тока I

д = Ар + ¿12Дф,

I = ¿12 Ар + ¿22 Аф.

Здесь Ар = Ур/, Аф = Уф/, / — длина канала. Значения коэффициентов ¿ц, ¿12, ¿21, ¿22 зависят от свойств жидкости и размеров канала

4 2 ^ 2 2 ^ 2

_ пееоа% пе2е0^ па2а

¿11 = , ¿12 = ¿21 = -;-ЬЪ ¿22 =--;-^2--;-Ьз

8^/ /

С, = 1 - 2J -líd,, a = 2j £

a o a

d

d(r/a) _

G = 2/ ch(01 dr.

0 0

Для их определения необходимо знать распределение потенциала в канале. Подставляя выражение (6) в (2), получаем уравнение Пуассона-Больцмана

1 ё / ёф\ 2геп0 ^ / геф \ г ёг V ёг / ее0 \кВТу

В области потенциалов ф, которая нас интересует (характерные значения потенциала от ~ 20 до ~ 70 мВ), функция гиперболического синуса может быть линеаризована

1А (г ёф^ = к2ф

г ёг V ёг ) '

где 1/к — толщина ДЭС, выраженная через параметр Дебая-Хюккеля к = (2г2е2П0/ее0квТ)1/2 [9].

Решив его с граничными условиями ф(а) = £, ф'(0) = 0 (условие равенства дзета-потенциала на стенке канала величине £ и условие симметричности решения относительно оси канала), получим распределение потенциала в сечении канала

ф = «Ш •

10

Здесь 10 — модифицированная функция Бесселя первого рода нулевого порядка, ка обычно называют безразмерным электрокинетическим радиусом (отношение радиуса канала к толщине ДЭС).

Теперь 64, Ь2, Ьз приобретут следующий вид

Gi =1 - ™, G = (ka)2

r/i2

fca/0

a

G =2/ — ch í zeC Jo(kr)A d, J a2 UbT/o(fca)/

При установившемся электрокинетическом течении в канале ток проводимости равен току протекания, следовательно, общий ток I в (7) равен нулю. Учитывая также, что Ьц = Ьц, найдём перепад давлений на концах канала

Арг = q/

Перепад давления на проводящем пути, составленном из чередующихся узлов и капилляров, выражается следующим образом [4,10]

N

-P = ^ДК,

i= 1

где i — номер капилляра, а N — общее число капилляров, составляющих рассматриваемую проводящую цепочку. Учитывая, что перепад давления —P имеет место на характерном макроскопическом расстоянии L = N1 (т.е. внешний градиент давления VP = ДР/L), и, поделив обе части (10) на L, получим

VP=^

i=i

Если минимальный радиус капилляра, принадлежащего данному проводящему пути, равен r, а радиусы капилляров распределены с плотностью f (r) (далее будет считаться, чтс^(r) — нормирована), то в (11) можно перейти от суммирования к интегрированию

vp = Т Др

Подставляя (9) в (12), получаем выражение для градиента давления вдоль проводящей цепочки

VP = 1 I (Li

1 - L12/L22)

r

Выразив из (13) расход по проводящей цепочке, имеем

1

J -^f

q = VP1

f (r) .dr

Вычислим количество путей с минимальным радиусом Г1, приходящихся на единицу площади поперечного сечения потока в пористой среде. Вероятность того, что капилляр имеет радиус г ^ Г1

ИЪ

Г1

Критический радиус гс связан с порогом протекания по решётке соотношением [4]

где Б — размерность пространства, в — координационное число решётки.

Для того чтобы посчитать общую проводимость решётки капилляров, необходимо ввести классификацию проводящих путей, которая позволит просуммировать их проводимости. В работе [4] была введена иерархия проводящих путей по минимальному радиусу входящих в них капилляров, согласно которой в процесс фильтрации сначала вовлекаются наиболее «толстые» цепочки, а затем постепенно «подключаются» все более «тонкие», то есть с меньшим минимальным радиусом капилляров. Следовательно, при любом Г1 < гс совокупность капилляров с радиусами больше г1 образуют бесконечный кластер — связную систему капилляров, по которой может осуществляться перенос флюида в системе.

Для того чтобы получить общую проводимость решётки, зафиксируем произвольный Г1 < гс. Соответствующий ему бесконечный кластер имеет нерегулярную ячеистую структуру с характерным размером ячейки [4,10]

Я = I

1 - ^ъ

(Г1) - жь(п)]

где V — индекс радиуса корреляции, равный 0, 85 для Б = 3. Тогда концентрация п(г 1) ориентированных коллинеарно действующему внешнему градиенту давления проводящих путей, составленных из пор с г ^ Г1, с учётом (15), (16) будет выражаться следующей формулой

/ Мёш .

п(г1) = Я =

1 - ^ъ

Следовательно, распределение числа проводящих путей (приходящихся на единицу площади поперечного сечения) по минимальному радиусу г1 есть

Р (п) = -

ёп(г1) ёг1

Просуммировав расходы по всем проводящим путям с использованием (14) и (18), получим закон фильтрации

т = V Р1

о

/

ёг

-1

Р(г) ёг.

Для того чтобы провести анализ эффекта электровязкости, удобно получить представление о зависимости скорости течения от градиента давления, в которое вязкость жидкости и проницаемость среды входят в явном виде. Если сравнить (19) с законом Дарси |го| = (К/р) |УР|, то видно, что

К =,

о

/

откуда

Рэ = К

' /

ёг

ёг

1

Р (г) ёг,

1

Р (г) ёг

1

В работе [10] получено выражение для проницаемости пористой среды, моделируемой пространственной решёткой капилляров

К = 27^

гЪ\ — 2v

с

/ (г)ёг

0 Г1

'/

оо оо

где I(г1) = П8 / /(г)/г4ёг(/ /(г)ёг)-1. Таким образом, подставив (21) в (20), по-

Г1

лучаем выражение

Г1

27vг-2

Ъл-2^

/ (г)ёг

Рэ =

0 Г1

'/Ы I

/

ёг

-1

Р(г) ёг

V

с

г

С ' с

с

г

определяющее эффективную вязкость электролита в процессе его течения в пористой среде.

3. Эксперимент

С целью проверки полученной зависимости для эффективной вязкости авторами был спланирован и проведён эксперимент по фильтрации дистиллированной воды и растворов соли ^С1 различной концентрации (10-3 Н, 10-2Н и 0,1 Н) через керновые образцы песчаника из скважин одного из месторождений Западной Сибири. Для каждого из образцов была получена порометрическая кривая, а (-потенциал определялся методом потенциала протекания, в котором учитывается вклад поверхностной проводимости в тонких капиллярах [11]. Характеристики образцов керна приведены в табл. 1.

Таблица 1

Характеристики образцов керна

Образец Пористость, д.ед. Плотность, кг/м3 Проницаемость по газу, мД

1 0,0974 2340 0, 99

2 0,12 2310 1, 29

Рис. 1. Порометрические кривые образцов. Пунктирной линией обозначен первый образец, сплошной линией — второй

Из приведённых данных (табл. 1, рис. 1) видно, что второй образец обладает большей проницаемостью, и в нём преобладают более крупные капилляры (2-4 мкм).

Установка, на которой проводился эксперимент (см. рис. 2), состоит из керно-держателя (1), электродов из нержавеющей стали (2), патрубков (3), микронасоса '^ЫМР-304" , ёмкости с раствором, измерительного капилляра (4) со средним диаметром 3, 5 мм и длинной 250 мм, со шкалой с ценой деления 0, 5 мм и абсолютной погрешностью измерения 0,05 мм. Измерение удельной электропроводности солевого раствора осуществлялось с помощью кондуктометра типа РИЬ-Т, концентрация определялась весовым методом на аналитических весах ET-600W. Относительная погрешность измерения объёмной скорости движущегося электролита не превышала 6%, относительная погрешность измерения проводимости и потенциала протекания не превышала 5%, а относительная погрешность измерения перепада давления — 4%.

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

4

Рис. 2. Схема экспериментальной установки: 1 — образец керна, 2 — электроды из нержавеющей стали, 3 — впускной и выпускной патрубки, 4 — измерительный капилляр

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

Рис. 3. Отношение эффективной вязкости к вязкости капельной жидкости, данные

эксперимента и теоретические кривые. Пунктирная линия — образец №1, сплошная — образец №2; pH=7,1,

Из приведённых графиков (рис. 3) видно, что эффект электровязкости более значителен при течении через керн с меньшей проницаемостью. Это обусловлено тем, что в этом образце присутствует большее число капилляров, радиус которых сопоставим с толщиной ДЭС. На этом же графике показаны теоретические кривые, полученные для следующих параметров: длина ребра решётки I =110 мкм (образец 1), I = 180 мкм (образец 2), координационное число в = 6, диэлектрическая проницаемость раствора е = 80, корректировочный коэффициент 7, учитывающий отсутствие перетоков между проводящими цепочками равен 1, 45.

4. Заключение

Построена перколяционная модель течения электролита в пористой среде. В рамках этой модели получена аналитическая зависимость скорости фильтрации раствора электролита от вида порометрической кривой, дзета-потенциала и концентрации ионов электролита. Результаты расчётов согласуются с данными проведённого авторами эксперимента. Установлено, что при наличии достаточного числа каналов с радиусом, сравнимым с толщиной ДЭС, скорость фильтрации

электролита значительно (в разы) замедляется по сравнению со скоростью фильтрации при отсутствии влияния ДЭС, что интерпретируется как возрастание вязкости (эффект электровязкости).

Литература

1. Tuckermann D. B., Pease R. F. W. High Performance Heat Sinking for VLSI // IEEE Electron Device Lett. — Vol. 5, No 2. — 1981.

2. Pfahler J. N. Liquid Transport in Micron and Submicron Size Channels: Ph.D. thesis / Department of Mechanical Engeneering and Applied Mechanics, Univ. of Pennsylvania. — 1992.

3. Peng X. F., Peterson G. P., Wang B. X. Frictional Flow Characteristics of Water Flowing Through Rectangular Microchannels // Exp. Heat Transfer. — No 7. — 1994.

4. Селяков В. И., Кадет В. В. Перколяционные модели процессов переноса в микронеоднородных средах. — М.: Недра, 1995. — 224 с.

5. Фролов Ю. Г. Курс коллоидной химии (Поверхностные явления и дисперсные системы): учебник для вузов. — М.: Химия, 1982. — 400 с.

6. Pride S. Governing Equations for the Coupled Electromagnetics and Acoustics of Porous Media // Phys. Rev. B. — Vol. 50. — 1994. — Pp. 15678-15696.

7. Her R. K. The Chemistry of Silica. — NY: John Wiley, 1979.

8. Revil A., Leroy P. Governing Equations for Ionic Transport in Porous Shale // Gepphys. Res. — Vol. 109. — 2004.

9. Синайский Э. Г. Гидродинамика физико-химических процессов. — М.: Недра, 1997. — 340 с.

10. Кадет В. В., Максименко А. А. Принципы аналитического описания течения жидкости в решёточных моделях пористых сред // Изв. РАН. — Т. 1. — 2000. — С. 79-83.

11. Lua F., Howa T. Y., Kwok. D. Y. An Improved Method for Determining Zeta Potential and Pore Conductivity of Porous Materials // Colloid Interface Sci. — Vol. 299. — 2006. — Pp. 972-976.

UDC 532.546

Influence of Salinity on Water Flow in Porous Media V.V. Kadet, A. S. Koryuzlov

Oil and Gas Underground Hydromechanics Department Gubkin Russian State University of Oil and Gas Leninskiy prospekt, 65, Moscow, Russia, 119991

The percolation model for electrokinetic flow of electrolyte in porous media is developed. In the context of this model taking into account ionic concentration of solution, zeta-potential, pore-metric curve and surface properties of porous media analytic dependence for filtration rate is obtained. The model results are in good agreement with experimental data obtained by authors (filtration of saline water in sandstone samples). Determined that in the presence of sufficient amount of capillaries with radius comparable with the thickness of EDL, the rate of filtration noticeably decreases as against the velocity of filtration without influence of EDL, that is interpreted as electro-viscous effect in porous media.