Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Физика и гидравлика нефтяного пласта

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
6.42 Mб
Скачать

Сравнивая (11.41) с (11.40), видим, что k\i/E= т есть время релаксации системы, т. е. время, в течение которого напряжение снижается в 2,7 раза от своей первоначальной величины.

Определим значения производных в точках / = 0 и t = to. Дифференцируя (11.39), имеем

 

 

-g - =

£ | _ e

 

**

 

 

e'(0)d0 +

e'(Oj.

(Н.42)

Подставляя

в

(11.42)

t =

0 и

t =

t0,

получим

значение

произ­

водных в этих точках:

 

 

 

 

 

 

 

 

.

 

da

 

 

 

 

 

 

 

 

 

 

 

 

^ “1 = чг u = - E z ' ( 0 ) = E ^

 

 

 

(11.43)

tga2 =

 

 

 

I

п

(/fi

 

~

 

de

 

 

 

 

 

 

 

 

ky. ■°° +

Е Ч Г

t=t,

гг~ + Е -^т-

t=ta

 

 

 

 

 

dt

 

По значениям

производ­

 

 

 

 

 

ных в точках t =

0 и t =

to,

d t

 

 

 

 

а также

по

отрезку

т

(см.

 

 

 

 

 

рис.

11.10)

не

представит

 

 

 

 

 

труда найти

искомые пара­

 

 

 

 

 

метры модели.

 

 

 

 

рас­

 

 

 

 

 

Последовательность

 

 

 

 

 

чета

такова:

 

 

второй

вет­

 

 

 

 

 

1. Ординате

 

 

 

 

 

ви

кривой

a =

F{t)

 

 

(на

 

 

 

 

 

рис.

11.10 равной

о0/е)

соот­

 

 

 

 

 

ветствует

аргумент — время

 

 

 

 

 

релаксации т.

 

 

 

 

 

 

 

 

 

 

 

2. Модуль упругости Е =

 

 

 

 

 

= tga.il (dz/dt) t= o-

 

 

 

Рис. НЛО. Кривая изменения напряжения при

 

 

 

заданной

скорости деформации (по модели

3. Если модель

выбрана

 

 

Максвелла)

 

правильно,

то

должно

по

 

 

 

 

 

крайней мере соблюдаться приближенное равенство

 

 

 

 

 

 

 

сто

|

\u-b/uijt=st

 

 

 

 

 

 

 

 

 

 

(de/rff)/ <

 

 

 

 

 

 

 

 

 

 

 

+

 

c w o , : , W - * * *

 

 

 

 

 

 

 

 

 

 

 

 

 

4. По известным значениям модуля £ и г находим величину

k\i =

тЕ.

точно

можно

 

определить

параметры при обработке

Более

 

всей кривой до точки А по способу наименьших квадратов по

формуле

(П.39)1.

Пусть

теперь задано напряжение a = a(t). Тогда, согласно

(П.38), имеем

e = -g -o (/) +

1

$ a(0)d0 =

^ [ a ( O +

S - ^ - d e l ,

(11.44)

 

 

о

L

о

J

 

Для использования способа необходимо в (11.39) задать —e'(t).

41

где о (< )— напряжение, действующее в момент времени t; 0 — «текущее» время.

Для конкретизации задачи примем следующую закономер­ ность изменения a (t):

 

a (t) =

a0t/t0

(0 < / <

/ 0)1

(11.45)

 

o(t) =

o0

 

(t> to )

У

 

 

 

Согласно (11.38),

будем

иметь

следующие значения

дефор­

мации:

 

 

 

 

 

 

 

 

e(C

=

 

 

1

 

(О < / < / „ ) ,

(11.46)

е(/,

to) =

-jh

1 +

t____ (о

(t > t 0)■

(11.47)

т

При стремлении

to-^-0

получим

значение деформации

е(0)

в начале времени при внезапном приложении постоянного на­ пряжения:

 

е(0) = о0/Е.

 

 

(11.48)

В дальнейшем постоянное растягивающее напряжение обус­

ловливает согласно (11.47) неограниченное

удлинение системы

с постоянной скоростью

 

ста

_ ст0

 

£(<) =

de

(11.49)

dt

тЕ

 

 

тем большей, чем больше приложенное напряжение и меньше вязкость системы.

Это явление называется ползучестью. Как следует из при­

веденного анализа [см. уравнение

(11.46)], в начальный момент

времени до достижения постоянного

напряжения

скорость

dtldt удлинения может быть неограниченно большой.

 

Пусть в момент времени t = t\ снято напряжение. Тогда при

t > t\ имеем

 

 

 

 

«(0 = - ^

= - ^ - = const.

 

(11.50)

Иными словами, после снятия напряжения деформация из­

меняется по сравнению с значением ее

(1 -[—

в конце пе­

риода /1 (рис. 11.11) на величину Оо/Е, а далее эта деформация не изменяется.

Параллельное соединение обоих элементов (упругого и вяз­ кого) дает модель Кельвина—Фойгта (рис. 11.12). Уравнение для этого случая выводится следующим образом.

Суммарное напряжение упругого и вязкого элементов

о = Оу„р + овязк.

(11.51)

42

Но Зупр—

авязк =

k^de/dt, а

потому

 

 

 

о =

Ее + ftp. dzjdt =

Е (е +

т dzjdt).

(11.52)

При заданной

скорости деформации

dz/dt =

z'(t) получим

искомое напряжение

j е' (0) db j=

 

 

 

о =

Е [те' (0

+

Е [те' (/) + z(t) -

z (0)] =

= E[Tz'(t) + z

Рис. 11.11. График развития деформации при за-

Рис. 11.12.

Модель вязкоупругого

данном ступенчатом изменении напряжения

тела

Кельвина—Фойгта

Если скорость деформации с самого начала не равна нулю,

то напряжение в

начале деформации

также отлично от нуля:

 

 

 

а(0) =

£те'(0).

 

 

 

(11.54)

В дальнейшем

монотонный

рост

напряжения

при

z'(t) =

= е'(0 ) =

const носит линейный характер:

 

 

 

 

 

а =

Е [те' (0) + е' (0) t\.

 

 

(11.55)

А между тем для модели Максвелла при тех

же

условиях

в начальный момент времени о(0) = 0 ,

а рост напряжения, со­

гласно уравнению

(11.41), подчиняется закономерности

 

 

о

ETZ' (0) (1 — е

т)-

 

 

(11.56)

Пусть

в некоторый

момент времени to

d z/ d t = 0

(что до­

стигается

закреплением

траверсы

модели

Кельвина— Фойгта

43

(см. рис. 11.12). Тогда, согласно

 

(11.52),

напряжение

о0=

= E e (t0) = E e '(0 )to (предполагается, что

до

этого

скорость

деформации была постоянна).

 

 

 

 

 

 

 

Для модели же Максвелла напряжение будет снижаться по

экспоненциальному закону (11.41).

 

a(t) решение (11.52)

имеет

При

заданном напряжении

a =

 

вид

 

/

/

в

 

 

 

 

 

 

I

 

 

 

 

 

 

X Г

X

 

 

 

(11.57)

 

Е = ~ЁГе

 

) е

 

° (0) dQ.

 

 

 

 

 

0

 

 

 

 

 

 

При t =

to удлинение, согласно

(11.57), будет

 

 

 

 

1

~

Г

 

о (0) dQ.

 

 

(11.58)

 

е° — Ех е

 

\ е

 

 

 

 

 

 

 

 

 

В этот момент напряжение снимается, и далее удлинение

будет изменяться по экспоненциальному закону

 

 

 

е =

е0е

 

 

-

 

 

 

(11.59)

При о =

Сто имеем из (11.57)

и (11.58)

 

 

 

 

 

е = _^ ~ (е '

-

l),

t < t o ,

 

 

(11.60)

 

Ь =

 

l),

t =

t0.

 

 

(11.61)

Математические выражения

моделей

Максвелла

и

Кель-

вина— Фойгта не всегда пригодны для описания реальных неф­ тей и особенно полимерных растворов.

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

ФИЛЬТРАЦИЯ вязкопластичной жидкости

Обобщим задачу фильтрации вязкопластичной жидкости. Рассмотрим стационарное движение вязко-пластичной жидкости в цилиндрической трубе (капилляре).

Распределение скоростей определяется из следующего урав­ нения [42]:

* - T z ! h - r !+ - ^ ' ' + c ' l" ' + c ’ - <IL62>

где то — предельное напряжение сдвига.

Для определения постоянных интегрирования Ci и С2 ис­

пользуем следующие граничные условия:

на стенке трубы

Г = г ь v = 0;

(11.63)

на некотором расстоянии Го от оси трубы г = го касательное напряжение т = то, что дает согласно уравнению Бингама— Шведова

 

т= То +

цп-g -,

 

 

 

(Н.64)

 

-Д^-= 0 при г =

г0.

 

 

(11.65)

Искомое расстояние г0 определим из условия равновесия

цилиндрического жидкого слоя

 

 

 

 

 

 

 

 

 

2яг0/.т0 =

яГо Ар,

 

 

 

(11.66)

откуда

 

 

 

 

 

 

 

 

 

 

 

 

Го=

-

^ .

 

 

 

 

 

(11.67)

Используя граничные условия

(11.63)

и

(11.65),

получим

 

 

 

 

 

 

 

 

 

 

 

(П'681

Выражение

(11.68) справедливо

для

Г\ ^ г

>

г0. При

г = г0

имеем из (11.67)

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

(п -69)

Для Г о ^ г ^ О

скорость

постоянна

 

и

равна

»о

(рис.

11.13).

Расход жидкости через трубу

 

 

 

 

 

 

 

 

Qi =

J 2яг© dr +

ягоТ)о-

 

 

(11.70)

Подставляя значение v и Vo соответственно из

(П.68) и

(11.69)

в уравнение (11.70), получим

 

 

 

 

 

 

 

 

 

я г 4. А р

 

 

 

 

 

 

 

 

(11.71)

Q' =

~ i k r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

45

Предположим, что в уравнении (11.67) ro = r\, тогда найдем наибольший перепад давления Д/?о, при котором вязкопластичная жидкость находится в предельном равновесии:

отсюда

 

 

Ар0 =

2£л0/г1,

 

 

 

 

 

(11.72)

 

г01г\ =

ApoAV

 

 

 

 

 

(11.73)

 

 

 

 

 

 

 

Заменяя в уравнении (П.71) г0/г\

через

Ар0/Ар,

получим

формулу Букингема

 

 

 

 

 

 

 

 

 

 

Qi

пг\Ар

1 - X

APo

1 ^ аро у

 

 

(11.74)

 

 

-

Ар

 

Ар /

 

 

 

 

Обычно третьим членом

этого

уравнения

 

 

можно

пренебречь,

поэтому

окончательно

 

 

имеем

 

 

 

 

 

 

 

 

 

 

 

Qi =

яН Др

1

 

-*£°Л

I

Др

«

Л

 

 

 

8|xnZ.

 

 

Ар

J ’

^

V *

 

 

 

 

 

 

 

 

 

 

 

(11.75)

 

 

В уравнении

(11.75)

отброшен

индекс 1

 

 

при г.

Используем

теперь

эту

зависимость

Рис, 11,13. К расчету

для получения

уравнения

фильтрации вязко­

фильтрации

вязко­

пластичной жидкости аналогичного уравнению

пластической

жидко­

сти в пористой среде Козени для вязкой жидкости. При этом заме­ нив в (11.75) величину L через \L (фактиче­ скую длину пути фильтрации) и, умножив на п (общее число

капилляров), получим

Q-

ппг*

Ар

1 —

_4_

Ар0

(11.76)

а ц

3

Др J

Далее и для этого случая остаются справедливыми следую­ щие соотношения:

п л гЦ

__

гп%

(11.77)

т = ■ Р

 

F

 

Из (11.75) и (11.77) получим

 

 

 

 

Q=='F (^S2p‘ ) l [ 5 ' [ 1

4

Др0

(11.78)

3

Др

 

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

 

* =

-d | r -

 

 

(Н-79)

Таким образом,

k Ар

 

 

 

 

Q = F

I -

4

Дро

(11.80)

(Хп7

3

Др

46

Уравнению (11.80) можно придать другой вид, используя (11.77), (11.78) и (11.79):

2gZ.To

ILXQS _ nfiLig

(11.81)

Лр0 =

2S*£

Окончательно имеем следующую формулу для определения расхода вязкопластичной жидкости через пористую среду:

Q = F

* Ар 2

Ftrfixо

(11.82)

 

Un*

 

 

Эта зависимость представляется в виде прямой линии, пере­ секающей ось абсцисс в точке с координатой

Дро= 4 '

Lm2т0

(11.83)

kSi

 

 

ФИЛЬТРАЦИЯ ж идкостей, ОПИСЫВАЕМАЯ СТЕПЕННЫМ ЗАКОНОМ

Приведенную выше методику используем для построения реологической модели, соответствующей степенному закону фильтрации (11.23) — (11.28).

(11.23) перепишем так:

 

Ч - 5 г ) л ,= т -

 

 

 

 

(и-84)

Из условия равновесия цилиндрического слоя имеем

 

 

 

т

Др

г-

 

 

 

 

 

 

 

т=

ж

 

 

 

 

 

С учетом (11.84) получим

Др

\‘/Я| 1/я,

 

 

 

 

 

dv

(

 

 

 

(11.85)

 

dr

Л 2*оL )

г

'

 

 

 

 

 

 

 

 

Разделяя переменные и учитывая, что при r =

r\

v =

0, най­

дем распределение скоростей

 

 

 

 

 

 

 

 

V 2*оL )

 

1

 

(rJ+ 1/"< _

r' + 1/«i)

 

(11.86)

 

1+ 1/Я1

 

 

 

 

 

 

 

Далее определим расход

\l/«i

г3-1" 17" 1

 

 

 

 

п

/

Д

 

 

 

 

Ар

 

 

3+1/Я| *

 

 

(11.87)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В уравнении (11.87)

отброшен

индекс

1

при

г.

При

п\ = 1

из (11.87) получим формулу Пуазейля. Умножая величину Qi на число капиллярных трубок, заменяя L через |L, определим общий поток:

пя.

 

 

Q

г Г •

(” -8«

3 + —

п\

 

 

47

Если в уравнении

Д

+ i - )

(11.89)

®ср —

 

(3+ v )

 

для ньютоновских

жидкостей « 1 = 1 , /го = Ц,

то оно превра­

щается в обычный закон Дарси.

 

Дадим общий вывод для произвольного

закона Y = M т).

Заменяя у через — dv/dr, получим

 

 

— а г - М т ) .

(П-90)

Напряжения трения на цилиндрическом элементе жидкости радиуса г в цилиндрическом капилляре и на стенке его соот­ ветственно равны

г Ар

_ г \ д р

т

г

 

(11.91)

т — 2L > T l — 2L

Т1

Г1

 

где г\ — радиус капилляра. Из (11.90) и (11.91) имеем

dv

(11.92)

dr

Интегрируя, получим

О (0 = Sfi (т| ~гу ) dr.

Расход жидкости в капилляре

Г\

Q\ = $2xrv(r)dr.

о

В результате интегрирования по частям получим

гi

Qi = * \ r2fi

У ) аг-

Используя (11.91), имеем из (11.95)

(11.93)

(11.94)

(11.95)

 

 

о Ар

 

 

 

2Ц.

(11.96)

Qi =

АрЗ

J т2fi (т) dj.

 

 

о

 

При этом, так же

как в

предыдущих

случаях, L заме­

нено \L.

 

 

 

48

Расход через п трубок Q = Q\n.

Используя соотношения теории фильтрации Кармана— Козени, получим после ряда преобразований следующую наибо­ лее общую зависимость для фильтрации:

УТ к'!*АР

т'1Ч

Q = Tj!s£

S

(11.97)

Проблема аномальных

жидкостей вызывает интерес

перс­

пективностью применения водорастворимых полимеров для за­ воднения нефтяных пластов. По некоторым данным «полимер­ ным» заводнением можно увеличить конечный коэффициент нефтеотдачи на 10%. Поэтому рассмотрение начнем с харак­ теристики полимерных растворов.

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

При сильно разбавленном растворе полимера описанная сет­ чатая (пространственная) структура может быть существенно упрощена, если исключить взаимное влияние полимерных цепей и рассматривать каждую такую цепь изолировано.

Сформулированные выше признаки полимера предопреде­ ляют такие его физические свойства, как гибкость цепи, т. е. способность ее изменять форму под влиянием теплового движе­ ния звеньев или внешнего поля; высокоэластичность — способ­ ность к очень большим обратимым деформациям; переход из высокоэластичного состояния в вязкотекучее и др.

Представляется возможным использование моделей, описы­ вающих состояние жидкостей для полимеров и их растворов [72] (см. рис. II.9 и 12).

Рассмотрим теперь полимеры, применяемые для загущения закачиваемой в пласт воды.

Наилучшие результаты показал полиакриламид (ПАА). Подвижность раствора этого полимера в пористой среде оказа­ лась меньше рассчитанной стандартными методами. Отношение первого значения подвижности к второму называется фактором сопротивления. Согласно [63], изменение подвижности объяс­ няется неньютоновским характером течения, а также адсорб­ цией и механическим улавливанием полимера в пористой среде.

Однако отсутствие необходимой информации об эксперимен­ тах исключает возможность оценки роли каждого фактора.

4 Заказ K s 283

49

По-видимому, для такой оценки необходимо получение по адекватной методике реологических параметров полимерного раствора.

В работе [63] справедливо отмечается, что для изучения реологической характеристики полимерных растворов необхо­ димы высокоточные вискозиметры, обеспечивающие поддержа­ ние постоянного расхода в процессе опыта.

Особенно важно точное значение параметров в области ма­ лых градиентов скорости, которые характерны для пластовых условий.

На высокоточном вискозиметре были определены характе­ ристики полимерных растворов. Ниже приводятся некоторые результаты этих исследований.

Анализ механизма течения разбавленных растворов полиме­

ров приводит

к выводу

об уменьшении эффективной вязкости

с увеличением

градиента

скорости (см. на рис. 11.4, кривую 2).

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

Такие клубки подвергаются воздействию пары сил (вслед­ ствие разности скоростей слоев) и перемещаются, вращаясь, вдоль слоев потока. Пока градиент скорости мал, это движение не сопровождается деформацией частиц, при этом в принципе асимметричные клубки ориентируются длинными осями в на­ правлении потока. С ростом градиента скорости происходит де­ формация полимерной частицы, возрастает ее асимметрия и, как следствие, увеличивается вязкость. Но, с другой стороны, ориентационные эффекты полимерных клубков должны способ­ ствовать уменьшению вязкости. В результате влияния этих двух факторов (деформации и ориентации) вязкость может изме­ няться с ростом градиента скорости в ту или другую сторону. Однако для реальных макромолекул влияние второго фактора сильнее, чем первого, поэтому в большинстве случаев вязкость сильно разбавленного раствора полимера должна уменьшаться при увеличении градиента скорости.

Указанное подтверждается специальными экспериментами, показывающими, что при сдвиговом течении водного раствора полиакриламида (ПАА) течение носило псевдопластический характер (см. на рис. II.4 кривую 2).

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

БО