Здесь 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) дает выражение для среднего дебита поступающей в залежь воды за интервал времени (t—At, 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(t—At)+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 ^д^
Исходные параметры для гидродинамических расчетов основных параметров разработки газовой залежи в условиях водонапорного режима приведены ниже.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.