Методы повышения эффективности процессов добычи и транспорта газа, страница 62

о

В случае если значение функции f(t) известно в дискрет­ные моменты U, равностоящие с шагом At, то /г-й момент <Xk вычисляется по формуле

aft=y] tkif(tt)M,                                                                                       (IV.9>

где N— число точек наблюдения.

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

Методика определения параметров однородного пласта при условии линейного радиального притока жидкости к скважине-заключается в следующем.

Вычисляются моментные характеристики:

а0 = J [Ар - р (01 dt,    аг = j [Ар - р (/)] Шг
о                                     о

(IV.10)

114


Q(t) tdt.

Здесь Ар — установившаяся депрессия на пласт до останов­ки скважины.

Параметры пласта вычисляются по формулам [30]:



 - PjAp - 2PO (a, - -^

&__    16   L     ______      V         V    / J___________________ (IV 11)


Здесь Q — дебит скважины до остановки.

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




х        25                      «oQo

где Qo — стационарный дебит скважины, приведенный к нор-нальным условиям; Ар2 — установившаяся разность квадрата пластового рпл и забойного р3 давлений;

«о = f [Pi - р2dt,     a[= $t [pi - ра (0] dt,
о                                               о

flo = J Q(t)dt,  p't = jtQ(t)dt;

p(t)—текущее забойное давление,  определяемое по кривой восстановления давления; Q(t)—дебит,   приведенный к нор-

kh мальным   условиям;   — —капилляропроводность; % — пьезо-

г*

проводность; R радиус залежи.

Пример IV.1. Ниже приведены данные о восстановлении дав­ления в газовой скв. 104 Вуктыльского газоконденсатного ме-

8*    П5


сторождения. Установившийся дебит скважины Qo = 525 тыс. м3/сут. Установившаяся разность квадратов забойного и пла­стового давлений равнялась примерно 15 600 (кгс/см2)2; А = = 28 (кгс/см2)2/сут/тыс. м3, Д=0,12   (кгс/см2) 2/сут2/(тыс.  м3)2.

t, сР% (кгс/см2)*    /, орс(кгс/см*)»

10

42 277

20

59 061

30

62678

40

62 910

100

63 137

160

63 400

220

63 664

340

64 392

520    64 591

820    65124

1 300   65 392

3100     65 861

6 700    66 331

19 300   66 872

40 900    67 075

—       67 080

Для подсчета параметров фильтрации в данном случае мож­но  пользоваться   следующими   формулами:   ао' = 22-1О6 (кгс/

/см2)-с, ai/=14«1010 (кгс/см2)-с. Отсюда находим #2/х = =40 000 с, kh/ii = 900 Д • см/сП.

ОПРЕДЕЛЕНИЕ  ФИЛЬТРАЦИОННЫХ  СОПРОТИВЛЕНИЙ

Уравнение, описывающее работу газового пласта, записывается в виде

где Т—характерное время переходного процесса; а иЬ — ко­эффициенты фильтрационных сопротивлений; Q(t)—суммар­ный дебит газа, приведенный к нормальным условиям; Др2(О = =Р2пл—P2a(t)t причем рпп — пластовое давление, p3(t) —забой­ное давление.

Коэффициенты а, Ь, Т определяются в следующей последо­вательности (И. М. Аметов и др.).

Вычисляются функции

t                                   t

ф (А = 1------------- ;      F(t) =

L(t) =—------------ .---------------------- (IV.20)

Так как давление и расход известны в дискретные моменты, то, записав уравнение (IV. 19) в момент 4 и вычтя из него ана­логичное уравнение, записанное в момент t\, получим

 - Ф (к) = a [F (tk) - F &)] + Ь [L (tk) - L &)].

1


Разделив обе части этого уравнения на F(th)—F(t{), полу-

чим

(IV.21)

F{tk)~-F(tt)            '      i

Таким образом, перестроив опытные данные в координатах Ф(0-Ф(^)       I До-


получим прямую линию.