Учет наноявлений в расчетах многофазной фильтрации в пористых средах Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

2558

Фундаментальные проблемы теоретической и прикладной механики Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2558-2560

УДК 532.5:532.6

УЧЕТ НАНОЯВЛЕНИЙ В РАСЧЕТАХ МНОГОФАЗНОЙ ФИЛЬТРАЦИИ

В ПОРИСТЫХ СРЕДАХ

© 2011 г. А.Я. Хавкин

Институт проблем нефти и газа РАН, Москва

[email protected]

Поступила в редакцию 24.08.2011

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

Ключевые слова: расчеты технологических показателей, добыча нефти, наноразмерные явления, капиллярный гистерезис, ионообмен.

Наноявления в системе «нефть-вода -порода»

К числу основных наноявлений при вытеснении нефти относятся капиллярный гистерезис и ионообмен между закачиваемой водой и пластовой системой «нефть -вода-порода» [1-8].

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

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

заводнения.

Анализ структуры распределений насыщенностей в пористых средах с учетом наноявлений выявил значительное влияние наноявлений на распределение фаз в пористых средах.

Гидродинамическая модель с учетом капиллярного гистерезиса

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

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

«нефть-вода-порода».

Предложенная автором запись закона фильтрации с указанием на различные размеры кластеров нефти (введением индекса у длины кластера нефти) имеет (без учета нелинейных эффектов) вид:

(1)

где Ц - скорости движения фаз (/ = 1 относится к водной фазе), к - проницаемость,/ - относительные фазовые проницаемости, цг- - вязкости фаз, 1Я - длина ганглии.

Скорость движения ганглии нефти Ц может быть найдена из соотношения

Ц2 = Ц1(В1 - В2 /1 X 1 = 1Я /1Ъ (2)

Bi =

hv-1

/1^2 Yi

Y 5 +

k/G Y /А

* Л

Y 4-

B2 =

^iUi ^iUi

/2m,iY2

2Y1Y3

lo(G2 ) =

k/1P12

MlU1

1 -

MaU1

-1

I1 = lo(0),

где О1 - начальные градиенты для фильтрации фаз, у;- - параметры, показывающие, в какое количество раз увеличивается от воздействия активной примесью вязкость нефти (у = 1), гистерезис капиллярного давления Р12 (у = 2), длина 1Я целика нефти (у = 3) и начальный градиент для фильтрации нефти (у = 4), уменьшается подвижность закачиваемого агента (у = 5); 10 и /1 - характерные линейные размеры неподвижных блобов неньютоновской и ньютоновской нефти, О3 - начальный градиент для фильтрации водной фазы.

Математическое ожидание объема закачанной воды в момент достижения всеми элементами длины больше Я добывающей скважины, выраженного в долях порового объема V, и математическое ожидание доли добытой нефти из целиков длины больше Я из всего пласта (нефтеотдача) П в момент V в принятых обозначениях (ф(х, I) -распределение ганглий по координате и длинам) равны:

П( Я) =

ад,

Конкретный вид ф(х, l) зависит от свойств флюидов и условий вытеснения.

Гидродинамическая модель с учетом ионообмена в системе «нефть—вода-порода»

Пусть ю = ю(с, 5) — количество адсорбированной воды, W — полное водосодержание единицы объема пористой среды, т0 — начальная пористость, Am — изменение пористости, U. — скорости фильтрации фаз (индекс i = 1 относится к воде); a = a(C, S) — количество сорбированных солей в единице объема пористой среды; D — ко -эффициент диффузии;/— относительные фазовые проницаемости; ^ — вязкости фаз; Pi — давление в фазах; P — капиллярное давление; а — эффективное поверхностное натяжение на разделе фаз; J — функция Леверетта, р* — отношение плотностей воды в поровом пространстве к плотности воды в адсорбированном состоянии.

Тогда система уравнений с учетом ионообмена в системе «нефть—вода—порода» имеет вид: dW

-----+ divU1 = 0, W = mS + ю, (5)

dt

m = m0 -Am, Am = р*ю,

d

— [m(1 - S)] + divU2 = 0,

dt

d

— (mSC + a) + div(CU1) = div( Dgrad С),

dt

Ui = -k (C, S) / (S, C)/ ц i (C )grad p., i = 1,2, p2 -P1 = P(S,C) = Gyfm/kJ(S,C).

[(1 -S0)HL0R1Y1 jjlW)Nq*x,l)dldx, (3) №2. С. 36—43'

0R

Список литературы

1. Хавкин А.Я. Государственный фонд алгоритмов и программ СССР, ПОО3958, 1979, 12 с. Алгоритмы и программы. Информ. бюлл. ВНТИЦентра. М., 1980. №1. С. 51.

2. Хавкин А.Я. // Нефтяное хозяйство. 1991. №5. С. 23-27.

3. Хавкин А.Я. // Нефть и битумы. Казань, 1994. Т. II. С. 458-463.

4. Хавкин А.Я. // Докл РАН. 1998. Т. 358, №2. С. 193-195.

5. Хавкин А.Я. Гидродинамические основы разработки залежей нефти с низкопроницаемыми коллекторами. М.: МО МАНПО, 2000. 525 с.

6. Хавкин А.Я. // Отечественная геология. 2007.

V (R) = -

Lp f1(R1 - У) г1 Ф(У, l)dl mL(1 - S0) J B1 R 1 - pl-1 W

7. Хавкин А .Я. // Вестник РАН. 2009. № 6. С. 519—522.

8. Хавкин А.Я. Наноявления и нанотехнологии в добыче нефти и газа. М.: ИИКИ, 2010. 692 с.

*

ACCOUNTING FOR NANOPHENOMENA IN THE ANALYSIS OF MULTIPHASE FILTRATION

IN POROUS MEDIA

A.Ya. Khavkin

As a result of the working experience in the oil industry, the author has come to believe that the discrepancy between the results of calculations of oil recovery technological parameters and the actual parameters is explained by not accounting for the physical and chemical nano-dimensional phenomena of multiphase filtration in porous media: a capillary hysteresis and ion-exchange between the injected water and the «oil-water-rock» system.

Keywords: technological parameters calculations, oil recovery, nano-dimensional phenomena, capillary hysteresis, ion-exchange.