Файл: Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 171
Скачиваний: 0
и при водонапорном режимах) реальные свойства газа учитывались лишь в уравнении материального баланса. При учете продвижения воды в газовую залежь получается экономия в 28 скважин на конец 15-го года разработки. Однако здесь не учтена возможность обводне ния и выбытия скважин из эксплуатации.
Из табл. 17 следует, что потребное число скважин в первые годы разработки как в условиях газового, так и водонапорного режима (для рассмотренных исходных геолого-промысловых данных) практи чески совпадает. На ранних этапах проектирования разработки га зовых месторождений отсутствуют достоверные данные о параметрах водоносного пласта и самой залежи. Поэтому, например, при соста влении проекта опытно-промышленной эксплуатации месторождения исходят из допущения о проявлении газового режима. Как видно из табл. 17, такое допущение в ряде случаев не вносит значительных погрешностей в расчеты.
§ 5. Определение показателей разработки газовой залежи в период падающей добычи газа при водонапорном режиме
Расчеты продвижения воды в газовую залежь при падающей до быче газа следует осуществлять на электрических моделях или ЭВМ исходя из рассмотрения двумерных задач с подвижной границей раз дела газ— вода. Это относится к крупным и средним по запасам газа месторождениям, так как к началу падающей добычи накапливаются значительная информация о залежи и данные о ее разработки. Для небольших месторождений с достаточно однородными коллектор скими свойствами при определении перспектив разработки или доразработки могут использоваться приближенные решения.
Для падающей добычи в случае водонапорного режима справед ливо следующее дифференциальное уравнение истощения газовой залежи:
Щ { 0= |
J _____ d_ |
р (t) ай (г) |
( 1) |
Par dt |
2 Ср ) |
Здесь ай (t) — текущий газонасыщенный объем залежи.
Разделим переменные и проинтегрируем уравнение (1) следу ющим образом:
t |
|
|
ай p |
|
Par J |
n ( t ) q { t ) d t |
|
(2 ) |
|
Vt-At |
2 Cp) |
|||
t-&t |
|
|
|
|
где |
|
|
|
|
ай (t) p (t) |
__ ай (t — At) p (t — At) |
|
||
z[p{t)] |
’ Jt A‘ |
*lP {t — At)] |
|
Примем, что число скважин в общем случае в период падающей добычи изменяется в результате обводнения и некоторого добури-
191
вания новых скважин с целью регулирования движением границы раздела газ— вода, приобщения недренируемых участков пласта и т. д. При решении задачи по шагам приближенно можно оценивать количество обводняющихся скважин за рассматриваемый период времени A t * . Зависимость п — п (t) считается известной, в частном случае п = const.
Пользуясь ранее применявшимся приемом, результат интегриро вания (2) записываем в виде:
Ра\ М - [n (t — At) q ( t — At) + |
n (t) q (t)] = |
__ aQ(t —At ) p (t — At) |
aQ (t) p (t) |
zip {t — At)] |
z[p(*)] |
Для простоты примем, что в период падающей добычи (вследствие невысоких пластовых давлений) можно пренебречь влиянием реаль ных свойств газов на форму записи уравнения притока газа к сква жине. Тогда из уравнения притока газа к средней скважине, запи санного с учетом, например, уравнения технологического режима эксплуатации скважин при заданной допустимой депрессии на пласт, находим
(4)
Значение текущего газонасыщенного объема залежи в момент времени t будет
|
aQ (t) = aQH— [(?в (t — At) + (q^(t — At) + |
A qB (t)) At]. |
(5) |
|||
Здесь использованы обозначения предыдущих параграфов. |
||||||
Подставляя (4) и (5) в (3), получаем |
|
|
|
|||
рат At |
n ( t — A t ) q ( t — At) + |
n ( t ) ^ — - ^ - + |
| / ^ |
р - + Й(2* д |
" б) |
|
2 |
||||||
|
|
|
|
|
||
|
ай (t — At) р (t — At) |
(aQH— [<?в (t — At) - f (gB(t - |
At) + |
|||
|
z [p («— Д*)1 |
|||||
|
z [p («)] |
|
|
|
||
|
|
+ Aq3 (t)) AH). |
|
|
(6) |
|
Давление на стенке укрупненной скважины (на расстоянии R 3) |
||||||
Р (# э , *) = Рн 2nkh А ?в(0 |
^ ( f o — fo„_x) |
2nkh |
2 Aq* iP ( f o - f o / . i ) . |
* Например, по данным разработки иногда удается получить четкую за висимость процента обводняющих скважин от процента обводнения газовой залежи (или количества поступившей в залежь воды) [21].
192
Здесь используются и значения AqBj(j — 1,2, . . .), найденные в результате расчетов для периодов нарастающей и постоянной до бычи газа;
fo — *0/1-! = - ^ ; A?bi = ?bi; 4 = 0; t = nAt.
Связь между дебитом воды в момент времени t и давлениями на расстояниях R3 и R устанавливаем, как и ранее, с использованием формулы Дюпюи:
?В (t — At) + AqB(t) = — 2лквк |
[p(R3, t) —(p (t) + PBy («))]. |
(8) |
|||||||||||
|
|
|
7 |
-П3 |
|
|
|
|
|
|
|
|
|
Уравнение (8) с учетом (7) |
записывается в виде: |
|
|
|
|||||||||
qB(t —At) + AqB(t)-. |
|
|
Rs |
Рв- |
|
М-в |
P(fo — fo„_i)AgB(t)- |
|
|||||
|
P b I n - |
|
2nkh |
|
|||||||||
|
|
Л (О |
|
|
|
|
|
|
|
|
|
|
|
/г-1 |
|
|
|
|
|
|
|
|
|
|
|
||
М-в |
2 |
А?в ,-Р(fo — foM) - |
(p (t) - |
pBy (*)) |
|
|
|||||||
2nkh |
/'=i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Отсюда имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь |
|
AqB(t) = |
- ^ - - ^ - p { t ) . |
|
|
|
|
(9> |
|||||
|
|
/г-1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сг = Са Рв — |
2 |
Л?в |
(fo _ |
f° м ) — рвУW |
|
|
|||||||
|
|
|
/=i |
|
|
|
|
|
|
|
|
|
|
|
|
—qB(t-At) In |
|
Лз |
> |
|
|
|
|
|
|||
|
|
|
|
|
|
л (i) |
|
2яРв*вЛ |
|
|
|||
' 2 |
R( t) “Г" |
к*” Р(1о — fo^j); |
С3~ |
|
|
||||||||
Со = 1п |
Лз |
|
|
|
|
|
|
|
|
|
|
|
|
Подставляем (9) |
в (6): |
|
|
|
|
l/ — |
|
|
|
||||
Рат Д* и (/ — A/) q (t — At ) ~ - n (t) ( 2S |
|
6(2p (Q - fl) |
|
||||||||||
|
У |
ш |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
В |
|
|||
|
|
aQ (t — A t)p (t — At) |
|
|
p (t) |
■ X |
|
|
|
||||
|
|
z[p(« — A<)] |
|
|
|
|
|
|
|
||||
|
|
|
|
z[p(«)] |
|
|
|
|
|||||
X (a£JH- [<?B(f - |
At) + |
qB( t - At) At + -g- At - |
-g_ p (t) A<] ) . |
(10) |
Уравнение (10) удобно решать путем подбора. Значение р (t) подбирается таким, чтобы тождественно выполнялось соотноше ние (10). При этом в первом приближении члены, зависящие от ре шения на момент времени t, принимаются согласно имеющемуся
13 Заназ 1013 |
193 |
решению на момент времени t — At. Иначе говоря, принимается, что
z [p (t)]^ z [p (t —At)]; |
R(t)f**R(t —At); у (t)^ y (t —At). |
(И) |
|
После определения |
р (t) уточняется величина коэффициента |
||
сверхсжимаемости газа |
z |
[р (£)]. С использованием формул |
(5) |
и (9) определяем текущее значение газонасыщенного объема порового пространства ccfi (t). Найденная величина ссй (t) позволяет по гра
фику зависимости у — f (aQ) уточнить величину подъема воды в залежь у (t). Суммарное количество воды, которое поступит в га зовую залежь ко времени t, будет
(?в (О = Q B(t— At) + [qB(t— At) - f AqB(i)] At.
Найденное значение QB(t) используется для уточнения величины текущего радиуса газоносности R (t).
Уточненные значения z [р (£)], R (t) и у (t) используются в рас четах второго приближения и т. д. Как обычно, процесс последо вательных приближений продолжается до тех пор, пока результаты расчетов в последнем и предпоследнем приближениях будут разли чаться на величину не более заданной погрешности е. При малом шаге At можно ограничиться первым приближением (11).
По найденной искомой величине р (t) и формуле (4) опреде ляем q (t). После этого переходим к расчетам для следующего вре менного слоя. В результате последовательного использования уравне
ний (10) и (4) определяем искомые зависимости р = р (t) и q = q (t). Тогда из уравнения технологического режима эксплуатации сква
жин при допустимой депрессии на пласт р (t) — рс (t) = б (() нахо дим зависимость рс = рс (t). Изменение отбора газа из месторожде ния в период падающей добычи определяется уравнением
Q{t) = n (t) q (t).
При учете разнодебитности скважин результат интегрирования уравнения (1) записываем в виде:
Q доб (0 Q доб |
— ДО |
а£3 (t — At) р (t — At) |
aQ (t) p (t) |
paTz [p (t — At)] |
• (12) |
||
|
|
paTz [ p ( f) l |
|
Тогда вместо |
уравнения (10) имеем |
|
a£2 (t— At) p (t — At)
< ? д о б (0 - < ? ц о б ( t - A t ) =
paTz [p (t — At)]
; - г | й г ( а о , - [ е . « - М ) + , , « - Д 0 Д < + - £ м - - ^ « ) м ] ) .
( 1 3 )
Решая квадратное уравнение (13) относительно р (t), получаем
? ( * ) = - i + V i |
+ т * |
(14) |
|
194
Со дг
С г Р а т 2 [ р (< )]
Рат2 [p |
(aQH— |
— M) + qB(t —&t)At + -^- A * ]); |
(0] |
|
|
__ |
aQ(t — Дt) p (t — Дt) |
+ <?доб № — A<) — <?доб ( 0 -
p aTz [p (f — AOJ
Перед тем как определить р (t) из уравнения (14), задаемся добы
тым количеством газа на момент времени t — фдоб (£). |
Кроме того, |
в первом приближении принимаем условия (11). |
|
Определяем р (t). С использованием уравнений |
притока газа |
к каждой скважине и соответствующих уравнений технологического
режима находим значения q{ (t) |
(i = 1, 2, . . |
п). По |
найден |
ным qt (t) уточняется значение добытого количества газа |
ко вре |
||
мени t, например, по формуле |
|
|
|
|
|
п |
|
<?дов(*) = едоб(*-Л *) + -4г |
£(* — Д*) + |
2 ? / ( 0 . |
|
|
|
1=1 |
|
Уточняются значения ъ [р (01» К (0 и У(0 и расчеты повторяются во втором приближении и т. д. в описанном уже порядке.
§ 6. Использование метода электроаналогий для решения двумерных задач
сподвижной границей раздела газ—вода
Впредыдущих параграфах рассмотрены достаточно простые инженерные методы расчета продвижения воды в газовые залежи
иопределения основных показателей их разработки. Эти методы достаточно общи, учитывают наиболее существенные факторы и более предпочтительны по сравнению с попытками отыскания аналити ческих решений соответствующих краевых задач. При отыскании аналитических решений возникают значительные трудности с учетом реальных свойств газа, противодавления поступающей в залежь воды, наличия защемленного газа и фазовой проницаемости в обвод ненной зоне пласта и т. д.
Остановимся на характерных допущениях и особенностях рас
смотренных решений задач по продвижению воды в газовые залежи. 1. Газовая залежь (практически вне зависимости от конфигура
ции) представляется в виде укрупненной скважины.
2. В пределах залежи фильтрация газа не рассматривается, т. е. не учитываются потери давления в области газоносности. (Достаточно схематично потери давления в периферийной зоне учи тываются при рассмотрении разработки залежи в условиях водо напорного режима и размещения скважин в центральной зоне [20]).
13* |
195 |