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

книги / Сборник задач по подземной гидравлике

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

Решение. Формулу дебита скважины с учетом сжимаемости можно получить из формулы Дюпюи, заменяя объемный расход Q расходом Qm, а давление р функцией Лейбензона Р.

~ __ 2яkh Рк — Рс

~р \nRJrc

Для жидкости, подчиняющейся закону Гука с уравнением состояния р= р0еРж(75-73‘>)>функция Лейбензона

Р = Г pdp +

С = f

р0еР,к (Р_Ро) dp = —

еРж <Р_Ро) + С,

J

J

 

Рж

 

а

 

 

 

 

__Рс — _Р°—

(V Po) — 0^ж

в

 

С

Рж

 

 

Раскладывая ех в ряд и ограничиваясь тремя членами раз­ ложения

(e- ' = ‘ + ‘ + f + - £ +

)•

 

получим

 

 

 

е»ж К-р.) _ е0ж (Рс-р.) = рж(Рк _ Рс) [1 +

(р„ +

рс- 2р0)'|

И

 

 

 

_2яЫгро__Г j

 

Рс).

I

•у- (Рк + Рс — 2р„)] (Рк —

L

 

 

ц In------

 

 

 

ГС

Давления в последней формуле абсолютные. Если положить Ро=Рат, то можно записать формулу для Qm через избыточные давления р,: й рс

Qm = -nkhp*T{р* -Рс) [l

+

(р„ + РС)] •

 

\1 1п——

 

 

 

 

 

гс

 

 

 

 

Разность между

объемным

дебитом

с

учетом сжимаемости

и дебитом, определяемым по формуле Дюпюи, равна:

&Qm

_ ftkh (рк

рс) Рж (рк +

рс) __

Рат

 

.

RK

 

 

 

 

р In------

 

 

 

 

 

гс

 

 

= 3,14-0,4.1,02-10-12.15» (11,76—7,35). 106-4,64. IQ-ю .(П,76 + 7,35). 1QQ

 

 

 

104

 

 

1,02-Ю -з-2,3 -lg-------

 

 

 

 

0,1

 

 

= 0,642* 10“ 4 м3/с = 5,55 м3/сут,

что составляет от дебита,

определяемого по формуле

Дюпюи

п _

2nkh (рк — рс) _

 

6,28>0,4-1,02-10-12.15-(11,76-7,35)-10«.0,864.105 = ^

^

 

104

' у

1,02-Ю -з-2,3-lg-------

 

величину

-^ 2 - :Q = 0,00445 = 0,445%-

Рат

Следовательно, при установившемся режиме фильтрации дебит можно определить по формулам для несжимаемой жид­ кости.

IX. УСТАНОВИВШАЯСЯ ФИЛЬТРАЦИЯ ГАЗИРОВАННОЙ ЖИДКОСТИ

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

При фильтрации газированной жидкости рассматривают от­ дельно движение каждой из фаз, считая, что жидкая фаза дви­

жется в изменяющейся

среде,

состоящей

из частиц

породы и

газовых пузырьков, а

газовая

фаза — в

изменяющейся среде,

состоящей из породы

и жидкости. Полагая, что

фильтрация

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

каждой фазы, вводя коэффициенты фазовых

проницаемостей

&ж и &г, которые меняются в пласте от точки к точке:

<2ж =

 

 

Рж ds

 

Qr = — —

- у - ® (s).

(IX. 1)

рг

as

 

Здесь Q / — дебит свободного газа в пластовых условиях. Опытами Викова и Ботсета установлено, что фазовые про­

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

пор в данном элементе пористой среды. В результате опытов построены графики зависимостей относительных фазовых про­ ницаемостей k^==kmlk и k * = k r/k от насыщенности а для несце­

ментированных песков (рис. 60), для песчаников (рис. 61), известняков и доломитов (рис. 62); здесь k — абсолютная про­ ницаемость породы, определяемая из данных по фильтрации однородной жидкости.

Рис. 60

Рис. 61

О 20 W 60 80 в,°/о

Рис. 62

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

Г = -^-г-ат .

(IX.2)

Qw

юз

 

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

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

можно связать с давлением и построить график (рис. 63), где безразмерное давление

а

£ _ р Рг

Рж

Назовем функцией С. А. Христиановича выражение

H = ]k'x dp.

(IX.3)

о

Через функцию Христиановича дебит жидкой фазы записы­ вается по закону Дарси, в котором роль давления играет функ­ ция Я:

<Эж = - — ^

(IX . 4)

Рж

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

Qx = 2nkh(Як — Яс) t

(IX.5)

^ж1п —

Гс

а дебит жидкой фазы галереи шириной В в пласте длиной I равен

<Эж= — (Як~ Яг) Bh.

(IX.6)

Рж I

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

Н = НК— Як~ Яс In А . ,

(IX.7)

1

г

 

In-----

 

 

Гс

а при параллельно-струйной фильтрации — линейному закону

и = Ик- - к~ Нг X.

(IX.8)

При расчетах по методу Б. Б. Лапука значения функции Христиановича находят следующим образом. Путем графическо­ го интегрирования строят безразмерную функцию Христиано­ вича

Н* = Р\ k*Kdp*,

о

используя график &^(р*). Зависимость Я* от р* представлена на рис. 64 для трех значений а = 5 р г/|1жРат {1 — а = 0,020; 2

сс= 0,015; 3 — « = 0,0Ю). Определяют величину ^ = Г — ,затем

переходят от размерного давления к безразмерному при помощи формулы

р* = —V ;

ох*9)

Рать

 

по рис. 64 находят значение Я*, соответствующее подсчитанному

значению р*. Переходят к

размерной

функции Христиановича

Я

= Я*£рат.

(IX. 10)

Для нахождения давления в некоторой точке пласта снача­ ла определяют значение функции Я' по формуле (IX.7) или (IX.8), затем, используя график зависимости Я*(р*) (см. рис. 64), переходят к соответствующему значению давления.

Отметим, что функция Хри­ стиановича зависит, кроме дав- //* ления (величины переменной в пласте), от постоянного парамет­

ра a = S — рат, где

5 — объем-

 

Цж

 

растворимости

ный коэффициент

газа в жидкости.

 

 

 

И. А. Парным было отмечено,

что зависимость Н*(р*)

соглас­

но графику (см. рис. 64) в ши­

роком

диапазоне

значений

р*

изображается

почти прямой

ли­

нией

(при рс/рк^0,2),

поэтому

приближенно

можно

принять,

что

 

 

 

 

 

 

Н* = Ар* + В

(IX. И)

и, следовательно,

 

 

 

Нк- Н с =

А (рк- р

с), (IX.12)

где Аж 0,944— 21,43 а.

 

что даже если давление в пласте

Г

Б. Пыхачев

отмечает,

меняется в широких пределах, фазовая проницаемость А*, изме­

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

нкн с = Ъ*ж(ркрс).

(IX. 13)

З а д а ч а 87

 

В пласте имеет место фильтрация газированной нефти. Оп­ ределить, при каких насыщенностях жидкостью и газом фазо­ вая проницаемость для жидкости km равна фазовой проницае­ мости для газа kT. Найти величину этой фазовой проницаемо­ сти, если абсолютная проницаемость пористой среды £ = 0,8Д. Рассмотреть случаи, когда коллектор представлен несцементи­ рованным песком, песчаником, известняками и доломитами.

Указание. Воспользоваться графиками зависимостей фазо­ вых проницаемостей от насыщенности жидкостью порового про­ странства (см. рис. 60—62).

З а д а ч а 88

Через пористую среду, представленную несцементированным песком, фильтруется газированная жидкость. Абсолютная про­

ницаемость пористой среды k = 5 Д,

вязкость жидкости рж =

= 1 сП,

вязкость газа цг= 0,012 сП,

насыщенность жидкостью

порового пространства о = 65%.

 

Определить фазовые проницаемости &ж и kT\ сравнить сумму

фазовых

проницаемостей с абсолютной проницаемостью пористой

среды, найти отношения скоростей фильтрации жидкости и газа w jw r и скоростей движения vm/vr.

Ответ: km= 1,15 Д;

&г=0,75

Д; w jw r = 0,0184; vnJvr =

= 0,00991.

 

 

 

З а д а ч а

89

В полосообразном

пласте происходит установившаяся па­

раллельно-струйная фильтрация газированной жидкости по закону Дарси. Ширина пласта В = 600 м, длина пласта L = 3 км, мощность h = 10 м, абсолютная проницаемость пласта k = 150 мД, коэффициенты вязкости нефти и газа в пластовых условиях со­ ответственно равны рж=1,12 мПа-с, рг=0,014 мПа-с, коэффи­

циент растворимости

газа

в

нефти

S=l,22-10~5

м3/м3*Па,

газовый

фактор Г = 350 м3/м3. Давление

на контуре

питания

р,ч- = 14,7

МПа (150 кгс/см2), на забое

галереи поддерживается

давление рг=Ю ,8 МПа

(110 кгс/см2).

 

 

 

Определить дебит галереи и давление в точке, расположен­

ной на расстоянии х=2,5 км от контура питания.

функции

Указание. Воспользоваться

графиком

зависимости

Я* от безразмерного давления р*.

 

м3/сут, р = 11,5 МПа.

Ответ: Qm= 61 м3/сут,

(Qr)aT = 21300

В центре нефтяного пласта радиуса 7?и = 350 м находится эксплуатационная скважина радиуса гс= 0,1 м.

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

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

на

забое

скважины рс= 8,82 МПа

(90 кгс/см2),

давление

на

контуре пи­

тания

рк = 13,2 МПа

(135

кгс/см2), абсолютная

проницаемость

пласта Л = 0,1

Д,

мощность

пласта h= 10 м, коэффициенты вяз­

кости

нефти р„=1,2 мПа-с

и

газа

рг = 0,012

мПа-с,

коэффи­

циент

растворимости

газа в нефти £ = 1,53-10~5

м3/м3-Па,

газо­

вый фактор Г = 400 м3/м3, рат= 1 ,0 1 .105 Па.

 

 

ниже.

 

Зависимость Я* от р* для а = 0,015 приведена

 

р*

0

1

2

 

3

4

5

6

7

8

 

9

10

н*

0

0,1

0,3

 

0,6

0,95

1,32

1,72

2,15

2,61

3,08

3,56

р*

12

14

18

20

22

24

26

28

30

32

34

н*

4,56

5,65

7,85

9

10,18

11,36

12,56

13,76

15

16,25

17,50

Решение. Дебит нефти при установившейся

плоскорадиаль­

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

2nkh (Як — Нс)

Q H =

1п Лк

для чего найдем значения функции Христиановича Я к и Я с при

давлениях рк и рс. Подсчитаем коэффициент a = S — р&т, котоЬн

рый является параметром при определении функции Христиа­ новича Я:

а = 1,53-10- 5 -

' 1.01

105 = 0,0154.

Определим значение безразмерного

газового фактора

| = -Ии-Г =

-400 = 4

Рн

Ь *

 

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

Рк

Рк _

J135_

32,8;

4-1,03

£Рат

 

 

 

Pi

Pc

^_90_

21, 8.

4-1,03

£Рат

 

 

 

По таблице зависимости между безразмерными значениями давления р* и функции Христиановича Я* при а = 0,015 найдем Я* = 16,75 и Я* = 10,06 и перейдем к размерным значениям

Як = Я*|рат= 16,75-4-1,01-105 = 6,77 МПа,

Я с = Я ^ р ат = 10,06-4-1,01-10е = 4,06 МПа.

При этом дебит нефти

= 6,28-9,1-1,02-10—и -10-(6,77— 4,06)-10» = j 7g- JQ_ 3 мз/с =

Л

350

1,2-Ю-з.2,3-lg ——

=

154 м3/сут;

дебит газа

Qr.aT = QHг = 154400 = 61 600 м3/сут.

Распределение функции Христиановича в пласте определя­ ется по формуле

 

Н = НК

Я к -Я с

■In Як

 

 

 

 

 

 

 

 

In Як

 

 

 

 

 

 

 

 

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

получим, задаваясь

различными

значениями г, определяя соответствующие значения

Я

и

Я*

 

 

при заданных

RKy

гс,

Як

и

Я с,

 

 

и по

значениям

Я* — значения

 

 

р* и р. Результаты расчетов при­

 

 

ведены в табл.

10.

 

 

 

 

 

 

 

 

 

 

 

Таблица

10

 

 

Г, м

я.

Н*

 

р*

р,

МПа

 

 

 

МПа

 

(кгс/см2)

 

 

0,1

4,06

10,06

21,80

8,82(90)

 

 

1,0

4,83

12,00

25,07

10,1 (103)

 

 

10,0

5,60

13,85

28.00

11,3(115)

Рис.

65

100,0

6,35

15,70

31,12

12,5(128)

350,0

6,77

16,75

32,80

13,2(135)

Для построения индикаторной диаграммы задаемся различ­ ными значениями рс и для этих значений по формуле

2nfeft|paT ( Н'к - Hi)

Qn-

llHIn /?кЛс

6 ,28-0,1-1,02-10-и. 10-4-1,01-10* (16,75 — / / ’ )

350

1,2-10-3.2,3-lg

0,1

= 0,264 •10-3 (16,75 — Hi)

м3/с)

подсчитаем дебиты

QH (табл. 11, рис. 65).

 

 

 

 

Таблица 11

рс, МПа (кгс/см2)

*

н1

QH, ма/сут

Рс

 

 

12,25(125)

30,4

15,25

34,3

10,78(110)

26,7

12,98

86,4

8,82 (90)

21,8

10,06

154,0

4,90 (50)

12,1

4,61

276,0

0,98(10)

2,43

0,43

372,0

0,101 (1,03)

0,25

0,025

382,0

 

З а д а ч а 91

 

 

В пласте имеет

место

установившаяся

плоскорадиальная

фильтрация газированной нефти по закону Дарси.

Др =

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

депрессии

= 25 кгс/см2 = 2,45

МПа и

заданном газовом факторе

Г =

= 200м3/м3 будет более высокий дебит нефти, если пластовые

давления различны: 1)

рк = 9,8 МПа (100 кгс/см2); 2)

ри =

= 4,9 МПа

(50 кгс/см2).

Коэффициенты вязкости нефти

р,ж =

= 1 мПа-с

и газа \хг= 0,012 мПа-с, коэффициент растворимости

газа в нефти 5 = 1,73-10-5 м3/м3-Па.

 

Указание. Воспользоваться графиком зависимости Я* от р*.

Ответ:

 

 

 

Фж1

_ (Нк

*с)1 _ JJy

<3ж2

( я * - я ; ) 2

 

З а д а ч а

92

Сравнить дебиты при установившейся плоскорадиальной

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

газовых

факторах и одной и той

же

депрессии.

Отношение

Р ж / Ц г = Ю

0 ,

коэффициент растворимости газа

в

нефти

3’ =

= 1,02-10~5

м3/м3*Па, рат= 9,8-104

Па,

давление

на контуре

питания

ри=»11,76 МПа (120 кгс/см2),

давление

на

забое

сква-

,жины рс = 9,8 МПа (100 кгс/см2). Газовые факторы П = 300 м3/м3

и Г2 = 600 м3/м3. Пласт

представлен несцементированным пес­

ком.

(Qr)aT2/QraT 1 = 1,33.

Ответ: Qw \/Qn{ 2 = 1 >5;

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

Найти средневзвешенное по объему пористой среды значе­ ние функции Христиановича Й и соответствующее ему значение

давления

при установившейся

плоскорадиальной

фильтрации

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

радиусом У?к=1 км,

если

давление

на

контуре

питания

рн= 10,29

МПа

(105

кгс/см2),.

давление

на

забое

рс= 8,33

МПа

(85

кгс/см2),

отношений

Мт/м<ж = 0,01,

коэффициент растворимости

газа

в

нефти

S=*

= 1,02• 10-5 м3/м3-Па,

газовый

фактор Г = 400

м3/м3,

радиус

скважины

гс= 0,1

м. Пласт представлен

несцементированным

песком.

 

 

 

 

 

 

 

 

 

 

 

Решение. Средневзвешенное по объему пористой среды зна­

чение функции Христиановича определяется по формуле

 

Н = — f

HdQ . --------- !-------

 

HK DHC ln- ^ Л 2лhmrdrtt

 

 

nhm (R l-r^ - Ц

 

г

I

 

 

 

 

ln-^-

 

 

 

 

 

 

 

 

 

Нк

Яс

 

 

 

 

 

H„

2 In R K

Найдем значения коэффициентов Цх-Spa

|-1ж

Е = = J iL . Н-ж

и безразмерные давления

Р к , .= 26,2,

|Рат

Р с

= 21,2.

Рс = ЪРа

По графику зависимости Н* от р* при а = 0,01 (см. рис. 64) найдем

Нк = И и Hl = 7,

откуда

Як = Я*£рат= 11-4-9,8-104 = 4,31 МПа,

/ / с = Я ^ р ат = 7-4-9,8-104 = 2,74 МПа

и

Н = 4,31

4,31 — 2,74

= 4,23 МПа,

1000

2-2,3• lg

0,1

по