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

книги / Нефтегазовая гидромеханика

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

В случае, если поверхность имеет произвольную форму и ха-

рактеризуется двумя главными радиусами R1 иR2 ,

то

P′ − P′′ = σ

1

+

1

.

(2.193)

 

R2

 

R1

 

 

Выражение (2.193) является обобщением формулы Лапласа (2.186) на случай произвольной формы поверхности жидкости. При R1 = R2 = R выражение (2.186) переходит в (2.193).

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

P′ − P′′ =

σ

,

(2.194)

 

R

 

 

где R – радиус цилиндрической поверхности жидкости. Термодинамический потенциал Гельмгольца позволяет также

определить зависимость поверхностного натяжения σ от температу-

ры T. Будем считать, что давления жидкости

Pи газа Р′′

поддер-

живаются одинаковыми:

P′ = P′′.

 

Тогда,

 

с учетом того, что

dV ′′ = −dV , имеем выражение для свободной энергии в виде

 

 

 

dΨ = −SdT + σdSповерхн.

 

 

 

(2.195)

Из этого выражения следует

 

 

 

 

 

 

 

 

∂σ

 

S

 

 

1

 

 

δQ

 

 

 

 

= −

 

 

= −

 

 

 

 

.

(2.196)

 

 

 

 

 

Sповерхн

 

 

T

 

 

 

 

 

 

T Sповерхн

 

T

 

 

Sповерхн T

 

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

 

δQ

 

 

 

qповерх =

 

 

 

,

(2.197)

S

 

 

 

 

 

 

 

 

поверхн T

 

 

то формула (2.196) примет вид

101

 

∂σ

= −

qповерхн

.

(2.198)

 

 

 

 

 

T Sповерхн

 

Т

 

 

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

Поэтому поверхностное натяжение σ уменьшается с повышением температуры, причем быстрота этого уменьшения обратно пропорциональна абсолютной температуре.

2.8.3. Взаимное вытеснение жидкостей и газов

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

102

Горизонтальное вытеснение нефти водой

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

 

= 1− σСВ − σОН ,

(2.199)

σн

где σСВ – насыщенность порового пространства связанной (реликтовой) водой; σОН – насыщенность порового пространства остаточ-

ной (невытесняемой) нефтью.

Уравнения, описывающие процесс движения жидкости при поршневом вытеснении, те же, что при фильтрации однородной жидкости. Задача решается методом смены стационарных состояний. Необходимо зафиксировать возможные положения фронта вытеснения. Обычно лучше всего задать некоторый шаг перемещения фронта вытеснения:

Lф = L.

В системе уравнений Ю.П. Борисова меняется лишь внешнее фильтрационное сопротивление перед первым действующим рядом (батареей) скважин:

1

=

 

μв

Lф +

 

 

μн

 

(Lк Lф ),

(2.200)

a

k h

 

a k h

 

 

 

 

 

 

 

 

 

 

 

 

=

 

μ

в

 

ln

r

+

 

μ

н

ln

rф

.

(2.201)

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2π k h

rф

 

2π k h

 

r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

103

Задача сводится к определению зависимости отбора жидкости из всех рядов от положения фронта вытеснения (текущего ВНК)

Qi = f (Lф ), Qi = f (rф ). Продолжительность интервала времени, в течение которого фронт вытеснения перемещается на один шаг,

 

 

 

 

m a h

L

 

 

 

 

 

 

t =

σ

н

,

 

 

 

 

 

 

 

Qi (Lф )

 

 

 

 

 

 

 

 

 

 

r )

2

2

 

 

σн

m π h (rф +

 

rф

 

t =

 

 

 

 

 

 

 

 

 

.

 

 

 

 

Qi (rф )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вертикальное вытеснение нефти водой

(2.202)

(2.203)

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

 

 

kг

f (η),

 

Wн = σн

m hн

 

(2.204)

kв

 

 

 

 

104

где

η =

b

,

b

вскрытая толщина пласта; hн – нефтенасыщенная

h

 

 

 

 

 

 

 

н

 

 

 

толщина пласта;

kг , kв – проницаемость горизонтально и вертикаль-

но напластованию.

Вытеснение газа водой

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

Положение контура вытеснения можно определить из уравнения материального баланса.

Для одномерного потока идеального газа

x

=

x2

+

 

k pатм Qг

 

t2 .

(2.205)

m2 μ

 

) a h

в

 

0

 

в

(L

x

 

 

 

 

 

 

 

x

0

 

 

 

Время вытеснения газа из пласта

T =

m a h (Lx x0 ) p0

.

(2.206)

 

 

pатм Qг

 

В круговой залежи зависимость между положением фронта вытеснения и временем выглядит следующим образом:

t = m r0 h

μ

в

r2

 

r

 

r2

 

r2

 

r

 

r2

 

 

 

 

0

ln

в

в

 

0

ln

0

0

.

(2.207)

k pатм Qг

2

rk

4

2

rk

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вытеснение нефти с помощью активных примесей

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

105

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

Полимеры – это высокомолекулярные соединения, органического и неорганического происхождения, состоящие из так называемых мономерных звеньев.

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

Адсорбция – концентрирование какого-либо вещества на поверхности жидкого или твердого тела в результате самопроизвольного перехода его из объема фазы.

Адсорбент – вещество, на поверхности которого происходит адсорбция.

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

Адсорбция на одной плоской поверхности определяется термодинамическим уравнением адсорбции Гиббса – отношением избытка вещества в поверхностном слое к площади, разделяющей поверхности, т.е. адсорбция

Г= −

∂σ

 

c

,

(2.208)

c

RT

 

 

 

 

где σ – поверхностное натяжение; с – концентрация растворенного в объеме вещества; R – газовая постоянная; Т – абсолютная температура.

Уравнение линейного вытеснения нефти водой, учитывающее сорбционные процессы:

(v c) = −m

sc +

b(s) g (s)

 

,

(2.209)

 

 

 

x

в

 

 

m

 

 

 

 

t

 

 

 

106

где vв – скорость фильтрации воды; vвc – объемная скорость фильтрации ПАВ; g (s) – количество вещества, адсорбированного объемом пористой среды при насыщенности раствором, равной s; b(s) – безразмерная, экспериментально определяемая функция, ха-

рактеризующая степень уменьшения адсорбции вследствие неполного контакта удельной поверхности породы с раствором.

2.8.4.Неизотермическая фильтрация жидкости

впористых средах

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

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

Полная энергия единицы массы системы может быть представлена в виде суммы трех членов:

E = Ew + Ez + U

,

(2.210)

A

 

 

где Еw – кинетическая энергия; Еz – потенциальная энергия; U – внутренняя энергия (термодинамическая функция внешних параметров Р, V и T); А – тепловой эквивалент работы.

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

107

При изучении нагревания пористого тела потоком фильтрующейся жидкости И.А. Чарным предложено использовать закон теплопередачи Ньютона в виде

q = α (T − θ),

(2.211)

где q – мощность теплового потока, нагревающего единицу объема пористого тела; α – постоянный коэффициент теплопередачи через единицу площади поверхности смачивания; θ – температура пористого тела; Т – температура жидкости; Ω – поверхность смачивания

вединице объема пористой среды.

Вдействительности процесс выравнивания температур намного сложнее и зависит также от размеров зерен и пор, коэффициентов теплопроводности и скоростей фильтрации.

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

Воснову вывода полного уравнения энергии для потока упругой жидкости в упругой пористой среде положено общее уравнение энергии потока сжимаемой жидкости, данное И.А. Чарным и записанное для конечного объема V:

 

P

+ U

 

v

2

 

 

z +

+

 

 

ρgdq +

ρg

 

 

 

A

 

2g

 

+

 

z + U

+

v2

 

ρgdV = Qвн + W

,

(2.212)

 

 

 

t

 

A

 

 

 

A

вн

 

 

 

 

2g

 

 

 

где z – высотная отметка; Р – давление; U – внутренняя энергия упругой жидкости; v – скорость потока; ρg – удельный вес жидкости; g – ускорение свободного падения; t – время; А – механический эквивалент тепловой энергии; Qвн и Wвн – тепловая и механическая мощности, подводимые извне к наблюдаемому объему.

108

В единице объема пористого тела объем жидкости соответствует пористости m.

Учитывая это и переходя для удобства дальнейших преобразований к пределу V 0, формулу (2.212) запишем в виде

 

P

 

U

 

v2

 

 

U

 

v2

 

 

m div z +

 

+

 

+

 

 

ρv

+ m

 

z +

 

+

 

 

ρg

=

ρg

A

 

 

A

 

 

 

 

2g

 

 

t

 

2g

 

 

 

 

 

 

=

 

1

 

Qвн

 

+ Wвн .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.213)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A V

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Раскрывая квадратные скобки в уравнении (2.213), получаем

 

 

P

+ U

 

 

v

2

 

 

 

 

 

 

 

 

 

 

 

P

 

+ U

 

 

 

v

2

 

 

mρg v grad z +

 

+

 

 

 

 

 

+ m z +

 

 

+

 

 

 

×

 

ρg

2g

 

 

 

ρg

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

2g

 

 

 

 

 

 

 

 

 

 

 

 

 

v

2

 

 

 

 

 

 

 

 

 

v

2

 

∂ρg

 

× div(ρgv) + mρg

 

z

 

+ U +

 

 

 

+ m z + U

+

 

 

=

 

t

 

 

 

 

 

2g

 

t

 

 

 

 

 

 

A

 

 

 

2g

 

 

 

 

A

 

 

 

 

 

 

 

 

=

 

1

 

Qвн

+ Wвн .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.214)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A V

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div(ρgv) = −

 

∂ρg

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.215)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то четвертое слагаемое в уравнении (2.214) сокращается, и тогда

 

 

P

+ U

 

v

2

 

 

 

 

 

 

 

 

 

 

+ U

 

 

v

2

 

 

mρg v grad z +

 

+

 

 

 

 

+ mρg

z

+

 

 

 

ρg

 

 

 

 

 

t

 

 

 

 

 

 

 

A

 

 

2g

 

 

 

 

 

 

 

 

 

A

 

 

2g

 

m

 

P

 

 

∂ρg

=

 

1

 

Qвн

+

 

Wвн

.

 

 

 

 

 

 

 

 

 

(2.216)

ρg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

A V

 

V

 

 

 

 

 

 

 

 

 

 

Используя подстановки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

P

 

∂ρg

 

= mρgP V

,

 

 

 

 

 

 

 

 

 

 

 

 

(2.217)

ρg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где – удельный объем вещества,

109

а также

 

ω = mv,

(2.218)

где ω – скорость фильтрации, с учетом основных гидродинамических соотношений:

I = U + A

P

,

(2.219)

ρg

 

 

 

 

 

Tds = dU + APdV ,

 

(2.220)

где s – энтропия системы, уравнение (2.213) приводится к виду

 

 

 

 

 

 

v

2

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρg w grad z +

 

+

+ m ρg

×

 

2g

 

 

 

 

 

 

 

 

 

A

 

 

 

 

v2

 

T s

 

1

 

 

(2.221)

×

 

z +

 

 

+

 

 

 

 

=

 

 

Qвн + Wвн.

 

 

A

 

 

A

 

t

2g

 

t

 

 

 

 

 

 

2.9. Интерференция скважин

 

 

Явление интерференции

(взаимодействия)

скважин

состоит

в том, что под влиянием изменения режима работы (в том числе пуска, остановки) одной скважины, называемой возмущающей, изменяются дебиты и (или) забойные давления группы скважин (реагирующих), эксплуатирующих тот же пласт (рис. 2.34). Допустим, что скважина 1 эксплуатировалась таким насосом, что динамический уровень в ней стоял у приема насоса Е. Горизонтальная прямая D1D2 указывает положение начального статического уровня. При работе скважины 1 в скважину 2 устанавливается «возмущенный» статический уровень С2С2. Кривая В1С2 является кривой распределения давления при работе скважины 1. Допустим, что скважину 2 также пустили в эксплуатацию таким же насосом на ту же глубину. Динамический уровень – В2В2. До пуска скважины 2 в эксплуатацию депрессия, а следовательно, и дебит, в скважине 1 определялась величиной S (разницей между статическим и динамическим уровнями). Если при работающей скважине 2 остановить скважину 1, то

110