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

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

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

2.7. Установившееся движение сжимаемой (упругой) жидкости и газа

2.7.1.Установившееся движение сжимаемой жидкости

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

Сжимаемость жидкости характеризуют коэффициентом объемного сжатия (упругости):

βж = −

1

 

dV .

(2.91)

 

 

V0

dP

 

Коэффициент объемного сжатия жидкости показывает, на какую часть от первоначального значения изменится объем жидкости при изменении давления на единицу. Знак «минус» означает, что с увеличением давления объем уменьшается.

Выразим объем флюида следующим образом:

V =

M

,

(2.92)

 

ρ

 

 

где М – масса рассматриваемого объема жидкости (считается постоянной).

Продифференцируем данное выражение:

dV = −M

dρ

.

(2.93)

 

 

ρ2

 

Подставляя эти выражения в формулу (2.91), получим

β

 

=

ρ

M

dρ

=

dρ

,

(2.94)

ж

M

ρ2 dP

ρ dP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

71

разделяя переменные в (2.94), имеем

dρ

= βж dP.

(2.95)

ρ

 

 

Полагая, что коэффициент объемного сжатия жидкости есть величина постоянная, и интегрируя полученное уравнение в пределах по давлению от Р0 (атмосферное давление) до Р (любое давление) и по плотности от ρ0 (плотность при атмосферном давлении) до ρ (плотность при давлении Р), находим

ρ

dρ

P

ρ

 

= βж (P P0 ).

 

 

ρ

= βж dP = ln

 

 

 

(2.96)

ρ

0

 

ρ0

 

P0

 

 

 

 

 

Из уравнения (2.96) имеем

 

 

 

 

 

 

 

 

eβж (PP0 )

=

 

ρ

,

(2.97)

 

 

 

 

или

 

 

 

 

ρ0

 

 

 

 

 

 

 

 

 

 

 

ρ = ρ0 eβж (PP0 ).

(2.98)

Уравнение (2.98) выражает зависимость плотности жидкости от давления. Такие формулы принято называть уравнениями состояния. Выражение (2.98) есть уравнение состояния упругой (сжимаемой) капельной жидкости.

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

Эффекты, связанные с изменением свойств жидкости и пористой среды в зависимости от давления, можно учесть, заменив давление на некоторую функцию:

72

 

k (P) ρ (P)

 

 

P =

μ (P)

dP + C,

(2.99)

 

 

 

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

Пренебрегая изменением проницаемости и вязкости при изменении давления (ввиду малости таких изменений), функцию Лей-

бензона можно записать в следующем виде:

 

P = ρ (P)dP + C.

(2.100)

Между движением сжимаемой и несжимаемой жидкости суще-

ствуют определенные аналогии:

 

 

 

 

 

 

 

 

 

 

 

Несжимаемая жидкость

 

Сжимаемая жидкость

 

Объемный расход Q

 

Массовый расход QМ

 

Давление Р

 

Функция Лейбензона P

Объемная скорость фильтрации w

 

Массовая скорость фильтрации ρw

С учетом уравнения состояния (2.98) функцию Лейбензона

можно определить по формуле

 

 

 

 

 

P = ρ (P)dP + C = ρ0

eβж (PP0 )dP + C =

ρ0

+ C.

(2.101)

 

 

 

 

βж

 

Для установившегося плоскорадиального потока сжимаемой (упругой) капельной жидкости массовый расход определяется по формуле

Q =

2πkh

 

P

P

(2.102)

μ

к

c

,

 

 

м

 

 

rк

 

ln rc

где Pк , Pc – значения функции Лейбензона соответственно на контуре питания и в скважине.

73

С учетом (2.98) формулу (2.102) можно записать в виде

Q =

2πkh

 

ρк − ρc

,

(2.103)

 

 

м

μ βж

rк

 

ln rc

где ρк, ρс – плотность жидкости соответственно на контуре питания и в скважине.

Изменение функции Лейбензона вдоль линий тока

P = P

Pк Pс

 

ln

rк

 

 

 

 

 

 

 

к

 

 

ln

rк

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

c

 

 

 

 

 

Рк Рс

 

 

 

 

 

Р = Р

+

ln

 

r

.

 

 

с

 

 

 

 

rк

 

r

 

 

 

 

 

 

 

 

 

 

c

 

 

 

ln r

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

Изменение плотности жидкости вдоль линий тока

ρ = ρк

 

ρк − ρc

ln

rк

 

 

 

 

 

r

 

 

 

 

ln

 

rк

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ = ρ

 

+

ρк − ρc

ln

 

r

.

 

 

 

 

 

c

 

 

 

rк

r

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

ln r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

(2.104)

(2.105)

(2.106)

(2.107)

2.7.2. Установившееся движение идеального газа

Газ отличается от капельной жидкости в первую очередь большей сжимаемостью. Согласно закону Бойля – Мариотта газ занимает объем, обратно пропорциональный давлению, под действием которого он находится (при условии сохранения постоянной температуры). Расстояние между отдельными молекулами газа велико, поэто-

74

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

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

P

=

Pатм

= RT,

(2.108)

ρ

ρатм

 

 

 

где ρ, ρатм – плотность газа соответственно при давлении Р и при атмосферном давлении Ратм; R – газовая постоянная; Т – температура.

Из формулы (2.108) можно получить уравнение состояния идеального газа:

ρ = ρ

 

 

P

.

(2.109)

атм

 

 

 

P

 

 

 

 

атм

 

По аналогии с движением сжимаемой жидкости при фильтрации газа используют функцию Лейбензона:

Р = ρdP + C = ρ

 

P

dP + C =

ρ

P2

 

(2.110)

 

 

атм

 

.

P

2

атм

 

 

P

 

 

 

 

атм

 

атм

 

 

 

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

При плоскорадиальной фильтрации газа по линейному закону его массовый расход определится по формуле

75

Q =

2πkh

 

Р

Р

πkh ρ

Р2

Р2

(2.111)

μ

 

к

 

 

с

=

атм

 

к

 

 

с

.

 

 

 

 

 

 

 

 

 

м

 

 

 

rк

μ P

 

 

rк

 

 

 

 

ln r

атм

ln r

 

 

 

 

 

 

c

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

Объемный расход, приведенный к атмосферному давлению:

Q

=

Q

=

πkh

 

P2

P2

 

м

 

к

 

 

c

.

(2.112)

ρ

μ P

 

 

 

 

атм

 

 

 

 

 

rк

 

 

 

атм

 

атм

 

ln r

 

 

 

 

 

 

 

 

 

 

 

 

 

c

2.7.3.Установившееся движение реального газа

впростейших фильтрационных потоках при линейном

инелинейном законах фильтрации

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

газа

в

недеформируемом пласте

(k = const) и при переменной вязкости газа (μ ( p) = ναr)

функция

Лейбензона может быть представлена в следующем виде:

 

р = k

ρ ( p)

dp + C =

kpатм

 

pdp

+ C,

(2.113)

 

pатм

μ ( p) z ( p)

 

μ ( p)

 

 

 

где z ( p) – коэффициент сверхсжимаемости реального газа (опреде-

ляется по графику Д. Брауни).

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

Q =

2πkh ( рк рc )

 

2πkhpатм

 

pк

 

p

dp.

 

=

 

 

 

(2.114)

 

 

 

 

 

 

м

ln

rк

 

 

 

Ратм ln

rк

 

p

μ ( p) z ( p)

 

 

 

 

 

 

 

 

 

 

rc

 

 

rc

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ =

μr + μc

и z =

zr

+ zc

,

 

(2.115)

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

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

76

Q

=

Q

=

πkh ( pк2 pc2 )

.

(2.116)

м

 

 

p

 

rк

атм

 

Р μ z ln

 

 

 

 

атм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

атм

rc

 

 

 

 

 

 

 

 

 

При нелинейном законе фильтрации на основе двучленной формулы

p2

p2

 

p

μ z

Q

r

 

p

Р

z β Q2

=

атм

 

атм

ln

к

+

атм

атм

 

атм

. (2.117)

 

 

 

 

 

 

 

 

к

c

 

 

πkh

 

 

r

 

2π2 h2 r

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

c

 

 

Решая последнее уравнение относительно расхода, можно выразить объемный дебит газовой скважины при нормальных условиях.

Распределение функции Лейбензона вдоль линий тока описывается формулой

P = P

Pк Pс

ln

rк

.

(2.118)

rк

 

к

 

 

r

 

ln rc

Распределение давления

P =

P2

P2

P2

ln

r

 

к

c

к

.

(2.119)

 

rк

 

 

к

 

 

r

 

ln rc

Рис. 2.24. Сравнение распределения давления для несжимаемой жидкости и идеального газа

77

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

2.7.4. Исследование газовых скважин методом установившихся отборов (при установившихся режимах)

Индикаторная диаграмма газовой скважины в координатах ( P, Qатм ) является параболой и обработке не подлежит (рис. 2.25). Для установления закона фильтрации газа индикаторные диаграммы принято строить в координатах (Pк2 Pc2 ,Qатм ). Линейность графика

в таких координатах свидетельствует о линейности закона фильтрации газа.

Рис. 2.25. Индикаторные диаграммы газовой скважины в условиях линейного закона фильтрации

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

P2

P2

= A Q

+ B Q 2

,

(2.120)

к

c

атм

атм

 

 

где А – коэффициент фильтрационного сопротивления, учитывающий потери давления на трение; В – коэффициент, учитывающий инерционную составляющую фильтрационного сопротивления:

78

A

=

μ Pатм ln

rк

 

,

 

(2.121)

r

 

 

 

 

πkh

 

 

 

 

 

 

 

 

c

 

 

 

 

B =

 

 

ρатм Pатм β

 

 

,

(2.122)

 

2 h2 r

 

k

 

 

 

c

 

 

 

 

 

где β – экспериментальная константа пористой среды.

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

a =

11,6μ z p

 

Rпр

 

+ C1

+ C2

 

 

T

 

 

 

 

атм ln

 

 

 

пл

,

(2.123)

rc

 

 

π kh

 

 

 

 

 

 

 

Tст

 

 

b =

 

ρатм z pат Tпл2

 

1−

 

rc

 

+ C + C

 

,

(2.124)

 

2 h2 l r T 2

R

 

 

 

 

 

3

 

 

4

 

 

 

 

c

ст

 

 

пр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где l – коэффициент макрошероховатости породы; С1С4 коэффициенты, учитывающие несовершенство по характеру и степени вскрытия в линейной и квадратичной частях уравнения притока; Rпр – приведенный радиус влияния скважины:

R =

Rσ

 

Qc

,

(2.125)

 

 

пр

2

 

Qc + 0,5Qσ

 

 

 

 

где Rσ – среднеарифметическое от расстояний до соседних скважин; Qc – дебит скважины; Qσ – суммарный дебит соседних скважин.

Индикаторная диаграмма газовой скважины, в случае фильтра-

 

 

P2

P2

 

ции по нелинейному закону, строится в координатах

 

к

c

,Qатм

Qатм

 

 

 

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

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

79

Рис. 2.26. Обработка индикаторной диаграммы по двухчленной формуле притока

Внекоторых случаях индикаторная диаграмма, перестроенная

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

Рис. 2.27. Качественный вид индикаторных диаграмм зависимости от факторов влияния: 1 – стандартная; 2 – очищение призабойной зоны; 3 – реальные свойства газа, скапливание жидкости или породы на забое, неполное восстановление пластового и забойного давлений, изменение фильтрационно-емкостных свойств; 4 – многопластовая залежь

Влияние реальных свойств газа. В ранее представленных фор-

мулах используются осредненные величины z и μ. В некоторых случаях изменение температуры и давления на разных режимах приво-

80