Фильтрация газов в трещиноватых коллекторах, страница 70

Здесь kB[pB(t)] — коэффициент фазовой проницаемости для воды в обводненной збне пласта при давлении в обводненной зоне /?в(0; Л — средняя эффективная мощность пласта; уъуУ)— противодавление поступившей в залежь воды на высоту у _удельным весом ув.

Значение kB(pB(t)) определяется из соотношения

К Й» (01 - К ехр {- акв - р, (t)]}.                                                                                       (66)

Здесь kQ — коэффициент проницаемости для воды при рп; €Lk — коэффициент изменения проницаемости*

Определить ah возможно двумя путями. В первом случае для коллекторов конкретного месторождения пк определяется по результатам газодинамических исследований скважин по мето­дике, изложенной в [74Q, во втором случае используются резуль­таты исследований кернов   рассматриваемого   месторождения.

Среднее давление а обводненной зоне пласта определяется л о формуле [58]


p(R8x)-p{t) In R3/R(t)


(67)


Здесь времени t соответствует условное время т.

Из (67) следует, что для определения давления в обводнен­ной зоне пласта принимается логарифмический закон (как и в •случае фильтрации однородной несжимаемой жидкости). Давле­ние на текущей границе раздела газ — вода приравнивается к средневзвешенному давлению в газонасыщенной части залежи.

Подынтегральную функцию k*(t) формулы (60) определяем, используя зависимость

k* (t) = ехр {- ak [pa - рвод (/)]}.                                                                                      (68)

Здесь />Вод(0—средневзвешенное давление в возмущенной части водоносного пласта (за пределами радиуса залежи R3). Одним из известных способов, например согласно [96], находим

(1


вод


 д (0] h \ (N (t) - 1) (


3   (N(t)-l)(VN(t)-\)


122


Здесь N


УШ_____ }.                                            (69)


Совместное решение уравнений (53), (54), (63) — (67) дает выражение для среднего дебита поступающей в залежь воды за интервал времени (tAt, t):

 (70)


Здесь А =


 Пп R3/R(t)


p (fp - fp^,) I

 К       У




o - fora_t) qB (т - At) -


(72)


1


(а^н ехр [— атн — р (t))]



- QB (* - Л/)) - z


2H


(t)


 (73)


Алгоритм расчета поступления в залежь воды базируется на методе последовательных приближений. Пусть на момент t—А^ решение задачи известно. Тогда методика расчетов для момента t следующая.

1. В первом приближении считаем, что

p(t)~~p(t-At), j3~pB(t-

Люд ^ Люд(0,   К(0 ^K{t — H{t)~R{t-At),  y{t)~y{t—Af),

z [p (t)) ~z(p(t- АО),   k* (0 c*k*(t- At).

2. По формуле (64) для времени т определяем безразмерное
время fo. Затем определяются значения функции p(fo). В пер­
вом приближении принимаем, что r(t)=x(tAt)+At.

3.  По формулам (70) — (73) определяем средний дебит по­
ступающей в залежь воды за промежуток времени [t—A^, t].

123


4.  По формуле (54) вычисляем суммарное количество воды,
(поступающей в залежь на момент t

5.  Уточняем aQ(t), используя следующую зависимость:

aQ (t) = aQH exp [- am (pH - p (t))] - QB (t).

6.  Используя (53), находим среднее давление в газовой зале­
жи на момент t

7.  Находим текущее положение границы раздела газ — вода
и высоту подъема y(t):

exp [— am (pn p {t))\ h(a — a0CT)

8. Определяем давление на начальном положении газоводя-юго контакта к моменту t:

9. По формуле (67) определяем pn{t).

10.  Используя рвО), по уравнению (66) уточняем значение

ЫО)

11.  По (69) определяем среднее давление в возмущенной зо­
не водоносного пласта рводО) и уточняем   значение   k*(t) по
формуле (68).

12.  Уточняем значение т, соответствующее t. Для этого (60)
запишем следующим образом:

t~ At                           t                         t—м

т=  Г k*{t)dt+ f k*(t)dt=--
o                      ^д^

Исходные параметры для гидродинамических расчетов основ­ных параметров разработки газовой залежи в условиях водона­порного режима приведены ниже.