
книги из ГПНТБ / Ширковский, А. И. Добыча и подземное хранение газа учебное пособие
.pdfПроизводя выкладки, аналогичные для случая фильтрации га за по закону Дарси, получим выражение для дебита газа, пола гая k, р и г независимыми от р и S
|
pI - pI = a q + bq \ |
|
(62) |
|
где при Rk^ R c и p* = const |
|
|
|
|
д _ П6р0г0Го Jn |
Rk_ . |
в _ |
Р*РаТо |
(63) |
я khpaTc |
Rc |
2я2Л2раТ^/?с0,746-104 |
||
|
р* |
63-10° |
|
(64) |
|
|
|
||
О п р е д е л е н и е г е о л о г о - ф и з и ч е с к и х |
п а р а м е т р о в |
п о р и с т о й с р е д ы по к о э ф ф и ц и е н т а м ф и л ь т р а ц и о н н ы х с о п р о т и в л е н и й / ! и В.
Определим коэффициент проницаемости из выражения для А
1 16[1q ZqT Q
(65)
nAhpnTc
где Si — коэффициент учета несовершенства скважины по степени
вскрытия |
(определяется из таблиц [1] по параметрам |
^ = ~ и |
|
_ |
D |
— коэффициент учета несовершенства скважин |
по ха |
Лс= у~); |
рактеру вскрытия пласта — определяется по приближенной форму ле Е. М. Минского
Здесь п — число работающих перфорационных отверстий, прихо дящихся на 1 м вскрытой мощности пласта; Ro— радиус полусфе рической каверны в пласте, прилегающей к перфорационному отверстию в металлической обсадной колонне и цементном камне за колонной (Ro зависит от типа перфоратора и крепости горных пород).
При применении перфоратора ПК-ЮЗ радиус полусферической каверны (в мм) можно определить из выражений:
для гранулярных газоносных коллекторов (песков и песчани ков)
R0 = 31,7m'и при 0,15 < т < 0,30; |
(67) |
* В Международной системе за единицу динамической вязкости газа р<> принят 1 Н -с/м2=1 Па-с. 1 пауз (Пз) = 102 сантипуаз (сПз) = 1 дин-с/см2=-- = 10-‘ Н • с/м2=10” 1 Па-с.
40
для карбонатных коллекторов |
|
|
|
||
Ro= \ 5 0 m u |
при |
0,01 < |
m < |
0,10. |
(68) |
Найдем отношение — из |
выражения |
для коэффициента |
В |
||
т |
|
|
|
|
|
63-io°pa7jj(i + j; + |
i;) |
2/; |
(69) |
||
2пг№Яст\ 0,746104раВ |
|
||||
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
(70) |
52 » |
— \ . |
|
|
(71) |
|
|
|
зпЩ30 |
|
|
|
Подсчитаем коэффициент пористости m путем деления значения к, определенного по выражению для А на величину klm, получен ную выше
m = |
(72) |
Определим структурный коэффициент т по формуле А. И. Ширковского
1,29-10з
(73)
(k — в Д, m — в долях единицы).
Вычислим коэффициент извилистости пористой среды
— = 0,645т0' 545, |
(74) |
где L и L0— длина соответственно криволинейного и прямолиней ного порового канала среднего радиуса. Рассчитаем коэффициент формы газонасыщенного коллектора
к = 1,55т0’455. |
(75) |
Определим удельную поверхность пористой среды (в см2/см3)
(76)
1,02-10-8
J V .
Найдем объем связанной воды SB путем подбора из уравнения
h |
|
0 , 1 8 2 - S B |
|
|
— = |
(1 _ |
SB). 10 571 |
. |
(77) |
m |
' |
в/ |
|
|
41
Определим коэффициенты газонасыщенности рг, |
абсолютной |
|
проницаемости k 0 и пористости т 0 |
|
|
Рг 1 *^в> k 0 |
тп |
(78) |
(1 ~ S B)2 |
(1 — S B) |
|
Нестационарная изотермическая плоскорадиальная фильтрация реального газа в однородном пласте при справедливости закона Дарси описывается нелинейным дифференциальным уравнением в частных производных параболического типа
1 |
d ^ |
k (p )p |
dp |
__ d |
m (p )p |
(79) |
|
r |
dr |
p, (p) z (p ) |
dr |
dt |
Z(P) J |
||
|
Интегрирование его возможно при применении методов ли неаризации или численных методов решения с использованием ЭВМ.
Обработка кривых нарастания забойного давления в скважинах без учета притока газа в ствол скважины при ее закрытии на устье по схеме бесконечного пласта проводится по формуле
Pi = P i 0 + |
PQo in |
|
+ 5 Qo + P' Jg |
t |
(80) |
|||
|
|
|
P in |
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
P\ = a |
+ |
P' Ig t, |
|
(81) |
||
где |
|
|
|
|
|
|
|
|
|
a — P3 о + |
PQo 1° |
„ |
+ B Q o , |
|
(82) |
||
|
|
|
|
|
P in |
|
|
|
|
p - |
t W |
a r ° |
, |
P' = |
2,3BQ0, |
|
(83) |
|
r |
2nkahTc |
|
K |
^ 0 |
|
|
|
Qo— установившийся дебит, с которым |
скважина |
работала |
до ее |
|||||
остановки (при |
^ = 0); |
р3.о — забойное давление, |
соответствующее |
|||||
дебиту Qo; Rc.n— приведенный радиус скважины |
|
|
||||||
|
|
kh |
42,4Q0paT0Zo |
|
|
/0<ч |
||
|
|
~ |
= — |
Р |
----- • |
|
(°4) |
|
|
|
Р |
|
*С |
|
|
|
|
При известном значении В определим коэффициент пьезопро |
||||||||
водности х и mh |
|
|
|
|
|
|
|
|
—— = 0,445 exp 2 3 |
|
|
x — kpn |
(85) |
||||
Rlr, |
|
|
|
|
|
|
m \i |
|
m h = |
2,25 — |
PH |
exp |
- 2 ,3 |
|
2\ T |
(86) |
|
|
|
P'
42
В формулах приняты следующие размерности: — |
в |
Qo |
|
|
(х |
сПз |
|
в тыс. м3/сут; Т0 и Тс в °К; k в Д; |
в 1/с; х в см2/с; |
Rc в см; |
mil |
в м; В в ( - W - )2- ( ^ ) 2.
тыс.м3 см2
Величины а и р определяются по кривой нарастания забойного давления после остановки скважины, построенной в координатах p i — lg^ (рис. 5) : а — отрезок, отсекаемый прямолинейным участ
ком этой зависимости на оси |
ординат |
(при |
t = \ |
с); |
Р’’— тангенс |
|||||
угла наклона этого прямоли |
|
|
|
|
|
|
|
|
||
нейного отрезка к оси времени. |
|
|
|
|
|
|
|
|
||
По кривой восстановления |
|
|
|
|
|
|
|
|
||
забойного |
давления в беско |
500 |
|
|
=fi‘=39,5 |
|
|
|||
нечном по протяженности пла |
|
|
|
|
|
|||||
000- |
/яiO3 |
|
|
|
|
|||||
сте без учета притока газа в |
|
|
|
|
||||||
скважину |
после |
ее закрытия |
300 |
|
|
|
|
|
|
|
можно определить запасы газа |
|
0,5 |
1,0 |
1,5 |
г,0 |
2,5 |
3,0 tyt |
|||
(в м3) в объеме дренирования |
О |
|||||||||
скважины |
при |
стандартных |
|
|
|
|
|
|
|
|
условиях |
|
|
Рис. |
5. |
График |
зависимости |
р3—lg i |
|||
2,3.2,25-10 |
|
для |
скв. 50 |
Султангуловского |
газового |
|||||
(87) |
|
|
месторождения |
|
||||||
Q3 = |
|
|
|
|
||||||
|
2р'0,864.105 |
|
|
|
|
|
|
|
|
где Qo— установившийся дебит, с которым скважина работала до остановки, в тыс. м3/сут; рп— статическое пластовое давление в объеме дренирования в кгс/см2; t0—-время полного восстановле
ния забойного давления до |
пластового в с; р '— тангенс |
угла |
на |
|
клона прямолинейного участка зависимости p \ —p\i}gt) |
к оси |
lg/ |
||
в (кгс/см2)2/цикл |
|
|
|
|
о |
/ ______2 , 3/lQ () |
|
(88) |
|
Р |
- |
2 In R J R C'B |
|
|
|
|
Здесь RK= 0,472 R — радиус контура области питания при уста новившихся условиях работы скважины; Rc,n— приведенный ра диус скважины, т. е. радиус такой гидродинамически совершенной по степени и характеру вскрытия пласта скважины, которая имеет одинаковый дебит с несовершенной скважиной радиуса Rc при одинаковой разности давлений рк — р3. Принимаем, что контурное давление рк приближенно равно начальному давлению рн
Ясп = |
Rcec. |
(89) |
R = Дс + |
/2 ,2 5 х* . |
(90) |
43
Запасы газа в удельном объеме дренирования скважины мож но определить по коэффициенту фильтрационного сопротивления А
|
2,25104р | t0 In -§^~ |
(91) |
||
Q3 = |
|
|
|
|
Из выражения (91) |
следует, что |
при одинаковых |
значениях |
|
In RK/Rc.n |
|
|
|
|
|
-=£i- = /lQoa _ |
|
(92) |
|
|
Р2 |
-^2^02 |
|
|
Коэффициент качества вскрытия пласта скважиной С можно |
||||
определить так: |
|
|
|
|
1 п - ^ |
=1п - ^ + С; |
(93) |
||
|
Rr |
Rc |
|
|
С = 2,3/lQ0 — In |
Rk |
(94) |
||
|
|
2Р' |
|
|
Пример 6. Определить коэффициенты абсолютной проницаемости и открытой |
||||
пористости, объем связанной |
воды, |
начальную |
газонасыщенность, |
структурный |
коэффициент и удельную поверхность пористой среды по данным исследования гидродинамически совершенной по степени и характеру вскрытия пласта скв. 50 Султангуловского газового месторождения по данным исследований на устано вившихся режимах.
Исходные данные для расчетов: |
рн= 24,4 кгс/см2; Д= 0,68; ра = 1,03 кгс/см2; |
|||||
h = 6,5 м; /?с=0,06 м; pic = 0,012 сПз; |
7'0=286°К; |
7'С= 293°К; |
Д = 1,04 сут/тыс. м3; |
|||
В = 5,68 • 10~4 (сут/тыс. м3)2; |
г„=1; |
Rк |
ct=472; |
|3' = 38,5; L = 260 м; |
||
= 1000; |
||||||
Qc= 203,5 тыс. м3/сут; р3.о=19 |
кгс/см2. |
|
|
|||
Определим коэффициент проницаемости из выражения для А |
||||||
|
116(г0Г0 In Rk_ |
116-0,012-286-2,3-3 |
|
|||
k = |
Rc |
|
0.428Д. |
|||
AnTchpa |
1,04-3,14-6,5-293-1,03 |
|||||
|
|
Найдем отношение kjm из уравнения для В
63-102рат2
Г- 2я2й2 Т* Rc0 ,746В
63-102-0,68-1.205-2862
2,28.
2-9,9-6,52-9,81 -2932-0,06-1,03-104-5,68-10~4
k
— = 1,73.
Определим пористость т, разделив k на к[т
0,428
0,247.
1,73
44
Вычислим структурный коэффициент
т = |
1,29-Юз |
10,7 — 3,22 lg |
(тГ
1,29-Юз
—— гу7 [Ю,7 ~ 3-22 (* -73)]- * = 9 -9- (1.73)
Найдем коэффициент извилистости
I = 0 ,645т0 ,545 = 0,645-9,90> 545 = 2,24.
Определим коэффициент формы газонасыщенного коллектора
Я = 1,55т0,455 = —- = — • = 4,42. £ 2,24
Вычислим удельную поверхность пористой среды
F =
к\ ,02-10~8 V/, У»
т
0,247
= 593 см2/см8.
0,428-1,02-10-« у /, 1/2
'9,9
0,247
Зависимости фазовой проницаемости к и пористости т от объема связанной воды 5 В примем в таком виде:
* •= M l — SB)2, |
т = |
та ( \— SB) , |
SB= |
0,182 — 0,1 lg — , |
||
Тогда |
|
|
|
|
|
т0 |
0,182—SB |
|
|
0,182—SB |
|||
|
|
|
||||
— (1 — SB)- 10 |
оЛ |
|
или 1,73 = (1 — SB)-10 |
ол |
||
Отсюда SB= 0,15. |
|
|
|
|
абсолютной |
проницаемости ко |
Определим коэффициент газонасыщенности рг, |
||||||
и открытой пористости т 0 |
|
|
|
|
|
|
pr = 1 — SD= 1 — 0,15 = |
0,85; |
|
||||
|
ко |
|
0,428 |
=0.592Д; |
|
|
|
(1 — 0,15)2 |
|
||||
|
|
|
|
|
||
|
т 0 = |
0,247 |
0,291. |
|
||
|
1 — 0,15 |
|
||||
|
|
|
|
|
|
Найдем геолого-физические параметры пласта kh/\i, к, х, mh, klm поданным нестационарных исследований, используя кривую нарастания забойного давле ния, показанную на рис. 5.
kh 42,4-203,5-1,03-286
=225Д-м/сПз;
I1 — |
38,5-293 |
|
225р0 |
225-0,012 |
|
Л |
|
0.416Д; |
|
6,5 |
45
х = 0,445/q exp |
2,3 ( a - p l o - B Q l ) |
|
|
|
|
|
P' |
|
|
: 0,445-36-exp 2,3 |
[472 — 361 — 5,68- IQ—4 (203,5)2] |
n |
3510; |
|
|
38,5 |
= |
||
|
|
|
|
|
kh |
pH„„ |
24,4 |
|
|
mh = |
^ - |
= 22 500— 7 - = 156,5 cm; |
|
|
(X |
x |
3510 |
|
|
|
к |
0,012 |
|
|
— = |
3510—----= 1,73. |
|
|
|
|
m |
24,4 |
|
|
Из данных расчетов следует, |
что значение к/т, полученное из выражений |
|||
для коэффициентов фильтрационного сопротивления А и В |
при |
исследовании |
скв. 50 на установившихся отборах, совпадает е этим же значением, получен ным при обработке кривой нарастания забойного давления.
Пример 7. Определить начальные запасы газа Султангуловского газового месторождения по кривой нарастания забойного давления и по коэффициенту фильтрационного сопротивления А в скв. 50.
Исходные данные |
для расчета: |
Q0 = 203,5 |
тыс. |
м3/сут; рн = 24,4 кгс/см2; |
Р' = 38,5 (кгс/см2)2/цикл; |
Л = 1,04 сут/тыс. м3 (кгс/см2)2; |
Rk |
||
— = 103; |
||||
По графику рис. 5 определяем 70= 1780 с. |
|
Rc |
||
|
|
|||
Запасы газа, подсчитанные после |
полной выработки месторождения |
|||
|
<2ф= 1,79-10° мз. |
|
|
|
Определим запасы газа по кривой нарастания забойного давления |
||||
2,3-2,25-107-203,5-24,42- 1780 |
1,67-10° м». |
|||
|
2-38,5-0,864-105 |
|||
|
|
|
||
Погрешность расчета |
|
|
|
|
|
( 1 ,7 9 - |
1,67)-100 |
|
|
6 н = ---------- --------------- = 6>7%- |
||||
Запасы газа, подсчитанные по коэффициенту фильтрационного сопротив |
||||
ления А |
|
|
|
|
2,25-107-24,42- 1780-2,3-3 |
1,825-10° ма. |
|||
Qa = —1 |
1 |
------ = |
Л1,04-0,864-105
Погрешность расчета
(1,79— 1,825)-100 ЬА = — ---------------------- = - 1 ,9 6 % .
Пример 8. Определить коэффициент качества вскрытия пласта С по следую
щим исходным данным: |
рн = 243 |
кгс/см2; 5 = 0,0148 |
(сут/тыс. |
м3)2 (кгс/см2)2; |
Qo= 1300 тыс. м3/сут; р |0 |
=4844 |
(кгс/см2)2; /?с = 6,5 м; |
|3'= 2200 |
(кгс/см2) 2/'цикл; |
время работы скважины до остановки ^р=45 ч; х=1500 см2/с. |
|
|||
Найдем Л<Эо |
|
|
|
|
AQo = р\ — Рз.о— BQl = 59 000 — 4844 — 0,0148-13002 = |
||||
= |
29 156 (сут./тыс.м3) (кгс/см2)2. |
|
46
Вычислим радиус влияния скважины |
|
|
||
|
R = - / 2,25-1500-45-3600 |
= 1160 м. |
|
|
|
RK= |
0,472-/? = 547 м. |
|
|
Определим коэффициент С |
|
547 |
|
|
|
2,3-29-156 |
|
|
|
~ |
2-2200 |
2,3 lg. 6,5-10-2 = 6 |
,2. |
|
§ 10. Определение |
числа работающих |
перфорационных каналов |
и радиуса полусферических каверн по изменению забойной температуры
При фильтрации газа в пористой среде температура его на забое скважины изменяется вследствие дросселирования газа и изменения его кинетической энергии
Д/0 = Д + Дt2, |
(95) |
где Д^о — общее понижение температуры газа |
на забое скважины |
при установившемся режиме в °С; Д^ — понижение температуры газа в процессе дросселирования в °С (оно определяется по эн тальпийным диаграммам или по среднему значению коэффициента
Джоуля — Томсона е); At2— изменение температуры, изменением скорости газового потока при движении питания до забоя скважины.
|
Д/1 = |
е(рк — Р3)- |
AmW32 |
|
Д/а = |
Ат К |
— Э |
||
2gCf |
21С~Р ’ |
|||
|
при wK<£w3.
вызванное от контура
(96)
(97)
Здесь Ат — механический эквивалент тепла |
в ккал/кгс-м; |
U7K, W3— средние истинные скорости движения газа |
на контуре |
питания и забое скважины в м/с; Ср — удельная теплоемкость газа при постоянном (забойном) давлении и температуре в ккал/кг-°С; g — ускорение земного притяжения в м/с2.
Измерив Д^0 при установившемся движении газа, можно опре
делить истинную скорость движения газа на забое |
|
||
Щ = j / " [Д/„ — е (рк — Рз)1 |
• |
(98) |
|
Просветность забоя скважины |
|
|
|
П = -^_ = |
■° 2зРа— , |
|
■(99) |
а»3 |
Рзг и Р 3w 3 |
|
|
где v3— скорость фильтрации газа на забое скважины; Q — дебит скважины при ра= 1 кгс/'см2 и / 3;' /^ —-площадь фильтрации на забое скважины.
47
Просветность забоя скважины выражается также зависимостью (поверхностная пористость)
U = F J F S, |
(100) |
где Fi — площадь всех перфорационных отверстий
= |
F3= 2nRch, |
(101) |
Приравнивая выражения (99) и (100), найдем, число работаюющих отверстий п
4Qz3pa |
( 102) |
|
p3wazaTid2 |
||
|
Определив число работающих перфорационных каналов, зная пористость т, величину и структуру коэффициентов фильтрацион ного сопротивления А и В, можно вычислить средний радиус полу сферической каверны Ro путем подбора из равенства
[ |
о |
lnt + E‘ |
1 |
лАНраТс |
nR0 |
||
|
63106p3Tq |
Rc |
|
|
3n*R3 |
||
= m |
|
|
|
|
|
(103) |
|
|
2.42'i2/?c T'0,746104paB |
где I' и £' — коэффициенты учета несовершенства скважины по степени вскрытия пласта.
Пример 9. Определить число работающих перфорационных каналов и сред ний радиус полусферической каверны по данным исследований скв. 7 Мессоях-
ского газового месторождения. |
/?с = 0,1 м; #„=100 |
м; |
Л = 2,15 |
(сут/тыс. м3) |
|
Исходные данные: /0=14°С; |
|||||
(кгс/см2)2; |
5 = 0,00375 (сут/тыс. м3)2 (кгс/см2); h = 10 |
м; |
ро=0,014 сПз; zo=0,83; |
||
m = 0,2; |
ср = 0,75 ккал/кг° С; |
е= 0,467° С/кгс/см2; |
|
Q= 522,4 |
тыс. м3/сут; |
р3 = 68,3 кгс/см2; р„ = 78,8 кгс/см2; ^3=9°С.
Перфорация скважины осуществлена перфоратором ПКС-80 с линейной плот
ностью «1 = 10 |
выстрелов на 1 м, диаметр |
отверстия d = 8 |
мм. |
Скважина совер |
||
шенна по степени вскрытия пласта (gi и |
=0). Истинная |
скорость фильтрации |
||||
на забое скважины |
|
|
|
|
|
|
ш3 = /[ ( 1 4 — 9) — 0,467 (78,8 — 68,3)] 2-9,8-0,6-427 |
= |
22,4 м/с. |
||||
Находим число работающих перфорационных отверстий |
|
|
|
|||
|
_________522,4-Юз-0,83-1,03(273+ 14) |
|
|
|
||
П ~ 68,3-22,5-0,864-105(273 + 20) 3,14-0,0082 |
= |
|
||||
65 |
65 |
|
|
|
|
|
или п0 = — |
== — = 6,5 |
|
|
|
|
|
Определяем средний радиус полусферической каверны |
|
|
|
|||
|
Rо |
116-0,014-287-0.83 |
_ 1 ______■ |
|
|
|
|
3,14.2,15-10-293 |
2,3 lg |
|
|
|
|
|
|
6,5#„ |
|
|
48
63’106-0,56-1.205-2872 |
( 1 + |
О-1 |
\ |
|
|
3 ,4 2 -2 ^ |
) |
= 0,2 _2-9,9-9,81-0,1-2932.3,75-10-3.1,03-102-0,746-104 |
|||
или |
|
|
|
0,596 |
|
0,1 |
*/» |
|
126,6/?§ ] |
|
|
|
|
|
|
Путем подбора находим /?0 = 0,037 м = 3,7 |
см. |
|
|
§ 11. Оборудование устья газовой скважины
Оборудование устья газовой скважины предназначено для об вязки верхних концов обсадных колонн и фонтанных труб, герме тизации межтрубного пространства и соединений, осуществления контроля за технологическим режимом эксплуатации скважин. Оно состоит из колонной и трубной головок и фонтанной «елки».
К о л о н н а я головка обвязывает верхние концы кондуктора и эксплуатационной колонны, герметизирует межтрубное пространст во, является опорой трубной головки с фонтанной елкой. Т р у б н а я головка служит для подвески фонтанных труб и герметизации межтрубного пространства между эксплуатационной колонной и фонтанными трубами. На нее непосредственно устанавливают фонтанную елку крестового или тройникового типа.
Основным элементом фонтанной елки крестового тина является крестовина, тройникового — тройник.
На рис. 6 изображено оборудование устья газовой скважины с елками крестового и тройникового типов.
Фонтанная елка крестового типа имеет небольшую высоту, удобна в обслуживании, уравновешена. Применяется она в том случае, если в потоке газа отсутствуют твердые взвеси, газообраз ные или жидкие коррозионные агенты, в результате действия ко торых крестовина может разрушиться и скважина выйдет из эксплуатации.
Фонтанная арматура тройникового типа имеет два тройника: верхний является рабочим, нижний — резервным, который исполь зуется только во время ремонта или замены верхнего. Арматура этого типа имеет большую высоту (до 5 м от поверхности земли), неудобна в обслуживании, неуравновешена. Применяется она в особо сложных условиях эксплуатации скважины — при наличии твердых взвесей в потоке газа, вызывающих абразивный износ оборудования, газообразных или жидких коррозионных агентов (углекислота, сероводород, пропионовая, масляная или другие кис лоты жирного ряда), при резких колебаниях давления и темпера туры. В верхней части фонтанной елки устанавливается буферный патрубок с манометром, показывающим давление на головке сква жины (буферное давление).
49