К решению задачи об изолированной неоднородности в пористой среде с электрокинетическими свойствами

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УЧЕНЫЕ ЗАПИСКИ ПЕТРОЗАВОДСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Декабрь, № 8. Т. 1 Технические науки 2014
УДК 53. 09
ВАЛЕРИЙ МИХАЙЛОВИЧ ЛЕВИН
доктор физико-математических наук, профессор, научный сотрудник, Мексиканский институт нефти (Мехико, Мексика)
vlevine @imp. mx
К РЕШЕНИЮ ЗАДАЧИ ОБ ИЗОЛИРОВАННОЙ НЕОДНОРОДНОСТИ В ПОРИСТОЙ СРЕДЕ С ЭЛЕКТРОКИНЕТИЧЕСКИМИ СВОЙСТВАМИ*
Решается задача об изолированной неоднородности в однородной изотропной среде с электрокинетическими свойствами. Выводятся интегральные уравнения для электромеханических полей внутри и в окрестности неоднородности. Для эллипсоидальной неоднородности эти уравнения решаются в явной аналитической форме.
Ключевые слова: электрокинетическая среда, изолированная неоднородность, электромеханические поля, эллипсоидальная неоднородность
Одной из важнейших задач фильтрации жидкости в случайно неоднородных пористых средах является замена реальной неоднородной среды на однородную с некоторыми эффективными (макроскопическими) свойствами и эквивалентной реакцией на внешние воздействия (задача гомогенизации). Решение задачи об изолированной неоднородности в неограниченной однородной среде является базовой для многих гомогенизационных схем [5]. В настоящей работе рассматривается задача об изолированной неоднородности в среде, обладающей электрокинетическими свойствами. Пусть в такой среде имеется замкнутая область V (включение) с другими электрокинетическими характеристиками. Векторы скорости фильтрации жидкости U t (x) и плотности электрического тока It (x) связаны в среде с неоднородностью с векторами градиента давления T (x) и электрического поля Et (x) линейными соотношениями [2] (х — произвольная точка в трехмерной среде)
Ui (x) = -kij (x)Tj (x) — av (x) Ej (xX (1
1 (x) = -atJ (x)Tj (x) — atj (x) Ej (x)
и удовлетворяют следующим дифференциальным уравнениям
Э kUk (x) = 0, Э kIk (x) = 0. (2)
В уравнениях (1) и (2) обозначено: kj (x) — отношение тензора проницаемости к коэффициенту вязкости жидкости, o. (x) — тензор электропроводности и atj (x) — электроосмотический тензор, Эk = Э / Эxk. Эти тензоры представляют собой кусочно-постоянные функции координат, принимающие значения kj, aH, o., если x е V, и значения
О о о УУУ
kj, ay, Oj, если x? V.
Обозначим через V (x) характеристическую функцию области V, занимаемой включением:
Г1 (x eV)
V (x) = (3)
10 (x? V). (3)
Эта функция позволяет представить тензоры kj (x), aj (x), Oj (x) в виде следующих сумм
k. (x) = k0 + klV (x), k1 = k. — k0,
ij ij ij ij ij ij
Oj (x) = a0 +a]V (x), al =au -al, (4)
Oj (x) = ° + 0v (x), Ol = Oj -O0.
Пусть p (x) и y/(x) представляют собой скалярные потенциалы полей T (x) и Et (x): T (x) = Эp (x), Et (x) = Э/(x). Подставляя соотношения (1) в уравнения (2) с учетом формул (4), можем записать
kj ЭtЭ jP (x) + a° Э (x) = f (x),
a° Э г Э jP (x) + 00 Э. Э/(x) = ф (x), где обозначено
f (x) = Э .V (x) [ kj Э. р (x) + a Э/(x)], Ф (x) = Э .V (x) Ю Э. р (x) + Oj Э/(x)].
(5)
(6)
Система дифференциальных уравнений (5) эквивалентна следующей системе интегральных уравнений
p (x) = p0 (x) + J G (x — x'-) f (x'-)dx'- + Jr (x
/(x) = /о (x) + Jr (x — x0 f (x'-)dx'- + Jg (x
x'-)m (x'-)dx'-,
(7)
x'-)p (x '-^dx.
В этих уравнениях p0(x) и /0(x) — поро-вое давление и потенциал электрического поля, которые были бы в среде без неоднородности при заданных условиях на бесконечности, G (x), Г (x), g (x) — компоненты функции Грина системы (7), удовлетворяющие системе уравнений
k0 Э гЭ jG (x) + 0° Э гЭ j Г (x) = -S (x), a° ЭгЭ jG (x) + 0 ЭгЭ j Г (x) = 0
kj ЭгЭ j Г (x) + a° Э гЭ jg (x) = 0,
a0 ЭгЭ j Г (x) + 0 ЭгЭ jg (x) = -S (x),
где S (x) — функция Дирака.
(8)
(9)
© Левин В. М., 2014
К решению задачи об изолированной неоднородности в пористой среде с электрокинетическими свойствами 109
Если основная среда изотропна, то есть
к0 = к, 8., а0 =а, 8., & lt-г° =ап8., (10)
то компоненты функции Грина имеют вид
G (х) = р-, Г (х) = -а 1
Д0 4яг'-
Д0 4лт'-
g (х) = К
1
(11)
Д" 4лг
, Д0 = Кр — а0, Г = |х|.
Пусть материал включения также изотропен:
К = К8., а =а8у, р = рЛ, (12)
к1 = к — к0, а1=а — а0, а1=а — а0.
Дифференцируя обе стороны уравнений (7) по координатам, получим
T (х) = T0(х) + JКу (х — X)[diTy (X) + d2Ey (х'-)]dx, (13)
V
E. (х) = E0 (х) + J К. у (х — х) [d3Ty (х'-) + d4Ey (У) ] dE. (14)
V
Здесь обозначено:
T0 (х) = д,-Р0 (x), Ei (х) =д& amp-0 (x),
(1 ^
К. (х) = дд
i у V '- I у
к4лг j
dx = Д- (р К — а0а), d2 = Д- (ра — ра), Д0 Д0
d3 =-(к0а1 — а0 к1), d4 = - (к0р — а1а0).
(15)
Если х е V, то уравнения (13) и (14) представляют собой систему уравнений для определения полей T (х) и Et (х) внутри V. Если эти поля внутри V известны, то поля T (х) и Et (х) вне этой области восстанавливаются из уравнений (13) и (14) однозначно. Таким образом, поля T (х) и Et (х) внутри V являются основными неизвестными задачи. Заметим, что ядро в этих интегральных уравнениях — формально реинтегрируемая функция с особенностью |х| в нуле. Для того чтобы придать смысл этому интегралу, будем рассматривать К. (х) как обобщенную функцию. Регуляризация интегралов, связанных с действием интегрального оператора с ядром К у (х) на гладкие финитные функции, приведена, в частности, в [1].
Введем следующие символические векторы и матрицы:
F (х)
T (х)
E. (х)],
F0(х)
T 0(х)
E,°(х) ],
M0 = -
дп


(16)
L1 =
к, а [ К у (х) 0 & quot-
1 1 а р, K (х) = у 4 у 0 Ку (х) _
Тогда два уравнения (13) и (14) можно записать в виде одного уравнения
F (х) = F0 (х) + J K (х — х) М 0L1F (х/)8х/. (17)
V
Уравнение, аналогичное уравнению (17), можно получить и для пары функций J (х) = = [(х), /(х)] в однородной среде с изолированной неоднородностью. Введем для этой цели соотношения, обратные (1), в той же символической краткой форме
F (х) = М (х^ (х), М (х) = [L (х)] 1,
L (х)
к (х) а (х)
а (х) р (х)
(18)
Умножая теперь обе стороны уравнения (17) на L0 и используя (18), получим
L0M (x)J (х) = J0(х) +
Г (19)
+J L0K (х — х,)M°L1M (х) J (х^х.
V
С учетом соотношений
L0M = L0 (M0 + M1) = I0 + L0M1, M1 = M — M0,
L1M = (L — L0) M = I0 + L0 M1 = -L0M1, (20)
где I0 — единичная 2×2-матрица, уравнение (19) можно переписать следующим образом
J (х) = J 0(х) — J S (х — х%° M1 J (Xd'-,
V
S (х)
S у (х) 0
0
Sv (х)
, S у (х) = 8у8(х) + К у (х).
(21)
В общем случае уравнения (19) и (21) могут быть решены лишь численно (эффективный метод численного решения уравнений такого типа описан в [4]). Для включений эллипсоидальной формы и полиномиальных внешних полей эти уравнения имеют явное аналитическое решение. В этом случае поля внутри включения также полиномы той же степени, что и внешние (теорема о полиномиальной консервативности, см. [3]). В частности, для однородных внешних полей поля F+=[T+, E.+ ] и J+ = [U. +, I+ ] внутри эллипсоидального включения также однородны и определяются выражениями
F +=Д^ F0, Д^ =(I + AM0 L1)-1, (22)
J + = Д-7 J0, Д-7 =(i + BL0M1)-1, (23)
A =
[ А 0 '- [ B. 0 '- [ 8. 0 & quot-
у, в = у, I = у 8у _
0 Ау _ 0 B у _ 0
. (24)
Здесь Ау и B у — тензоры с постоянными компонентами. В системе координат, связанной с главными осями эллипсоида, эти тензоры представляются в виде (по индексу i не суммировать!)
110
В. М. Левин
Aik _ Ai8ik, Bik _ 8ik Aik
A _-Ьт J
^ л
d?
(Or + ?)& lt-J (aL + %){(?, + %){al + ?,)
(25)
(L6)
T+ = D [{8km + d4Akm) E — drAkE ] ,
E+ _ D. 71 [-d2A, T° + (8,+ dLA,)E0 ]
i ik [ 2 km m km 1 km f m ]
U+_ C- [{ + c4Bm Em — ClBkml
0
m I ?
i: _с- [-C3 Bk"um+{+Cl Bkm e: ]. Здесь обозначено
Dik _ (8iP + dL Aip)(8pk + d4 Apk) — d2d3 AipApk, Ck _ (8p + C1Bip)(8pk + C4Bpk) — C2C3B pBpk,
CL _ ko^i — (X0P1, Cl = ^01 — KPl,
C3 = a04 — ®ofil, C4 _®0K1 — ®oPl,
°0 A, а а0
A 0' _a — a0
k0 A b 1 b II
A o,
Полученные решения покрывают широкий спектр форм неоднородностей: сфера, цилиндр, эллиптическая игла, диск. В частности, для сферической неоднородности Aj = 8 / 3, Bij _ 28ij / 3 и формулы (27) и (28) преобразуются в следующие
где aL, а2, а3 — полуоси эллипсоида. В развернутой форме формулы (22) и (23) имеют вид
(27)
(28)
(29)
(30)
T +__ 1
'- _ D

1+4
3
Т& lt-0 — Е e 0
i 3 i
-^ T0 + (d ^ 1 + 2 1 0
_ 3 i 1 3 7
e: _-
i D
D _ L + d1 + d4 + k1°1 -а1
(31)
u+ _-
i C
LV
1 C
2c.
1 + -43
/
2C
-23 U0 + 3i
9A"
Lc
U 0 -f2_ 10
i 3 i
1+lCl 3
(32)
C _ 1 + & quot-L (C1 + C4) + 9A0(51*L — ^1L).
Таким образом, при равенстве «совместной» постоянной, а нулю полученные формулы распадаются на известные решения о неоднородности в среде с гидравлической и электрической проводимостями.
* Работа выполнена при поддержке Программы стратегического развития ПетрГУ на 2012−2016 гг.
СПИСОК ЛИТЕРАТУРЫ
1. Михлин С. Г. Многомерные сингулярные интегралы и интегральные уравнения. М., 1962. 254 с.
2. AdlerP., Mityushev V. Effective medium approximation and exact formulae for electrokinetic phenomena in porous media // J. Phys. A: Math Gen. 2003. № 36. P. 391−404.
3. Eshelby J. The determination of the elastic fields of an elliptical inclusion, and related problems // Proc. of the Royal Soc. of London. 1957. A241. P. 376−391.
4. K a n a u n S. K., Levin V. M. Self-Consistent Methods for Composites. Vol. 1. Static Problems. Springer, 2008. 376 p.
5. Kanaun S. K., Levin V. M. Effective field method in the theory of heterogeneous media // Kachanov M., Sevostianov I. (Eds.). Effective Properties of Heterogeneous Materials. Springer, 2013. P. 199−283.
Levin V. M., Mexican Oil Institute (Mexico City, Mexico)
ON PROBLEM SOLUTION OF ISOLATED INHOMOGENEITY IN POROUS MEDIUM UNDERGOING
ELECTROKINETIC PHENOMENA
The problem of isolated inhomogeneity in porous medium with electro-kinetic phenomena is solved. The integral equations for the electro-mechanical fields in such medium are developed. For the ellipsoidal inhomogeneity these equations are solved in an explicit analytical form.
Keywords: electrokinetic medium- isolated heterogeneity- electromechanical field- ellipsoidal heterogeneity.
REFERENCES
1. M i k h l i n S. G. Multi-dimensional singular integrals and integral equations. Мoscow, 1962. 254 p.
2. A d l e r P., M i t y u s h e v V. Effective medium approximation and exact formulae for electrokinetic phenomena in porous media // J. Phys. A: Math Gen. 2003. № 36. P. 391−404.
3. Eshelby J. The determination of the elastic fields of an elliptical inclusion, and related problems / / Proc. of the Royal Soc. of London. 1957. A241. P. 376−391.
4. Kanaun S. K., Levin V. M. Self-Consistent Methods for Composites. Vol. 1. Static Problems. Springer, 2008. 376 p.
5. Kanaun S. K., Levin V. M. Effective field method in the theory of heterogeneous media // Kachanov M., Sevostianov I. (Eds.) Effective Properties of Heterogeneous Materials. Springer, 2013. P. 199−283.
Поступила в редакцию 11. 07. 2014

ПоказатьСвернуть
Заполнить форму текущей работой