книги / Сборник задач по подземной гидравлике
..pdfРешение. Формулу дебита скважины с учетом сжимаемости можно получить из формулы Дюпюи, заменяя объемный расход 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
по