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

(169)

для зависимости коэффициента макрошероховатости / от дав­ления имеем

7/ \       1 л—5,15 i   i    \аь 11,27                                                                                                                                            /17ГП

1{р) = Ю     lko(p!pH) ko\    .                                                                                        (1'0)

Получим формулу притока газа к скважине при начальном снижении давления от ри до текущего р. Зависимость пара­метров пласта и газа от давления считаем степенной. Уравне­ние (165) нельзя решить аналитически. Как показано в работах [48, 59, 74J, нарушение закона Дарси существенно в призабой-ной зоне, т. е. второй член в (165) в наибольшей степени акти­вен в призабойной зоне скважины, поэтому можно принять, что

k(P)        _ А_(Рс)______ ,h/r,\._________ (171)

1 (Р) ^ (Р)      t (Рс) I

После интегрирования уравнения (165) с учетом (169) — (171) в пределах по г от Rc до /?„ и по р от рстши до р имеем

,     Rk


ч


(172)


1S7


где

/0= Ю-515Й'27.                                                                                       (173)

Уравнение (172) справедливо для упругопластичной пори­стой среды при рс^Рстт- Уравнение (172) впервые получено в диссертации Н. М. Кульпиной и описывает фильтрацию газа при нарушении закона Дарси в упругодеформируемой пористой среде.

При значениях забойного давления рс выше ранее достигну­того минимального рс min необходимо учитывать возникшие упругопластические деформации коллектора. В этом случае закон изменения проницаемости по радиусу имеет вид (150). Коэффициент проницаемости при учете упругопластических свойств пласта зависит от распределения минимального дав­ления по координате г.

Распределение минимального давления pmin{r) вокруг сква­жины по радиусу определяется из выражения

... —ч

где

Л ___     ЦогнРат  1„   Rk .

2nk0hpl

2нРатРат

Распределение коэффициента проницаемости по координате г для упругопластичного пласта имеет следующий вид:

k (р, г) = k0 [pmin (r)lpHfk о-°ь [P (r)/pH)aK                                                                                       (175)

Здесь

Уравнение, описывающее установившую фильтрацию ре­ального газа при нарушении закона Дарси в упругопластичной пористой среде, запишется в виде:

/ Р (Г) \ak~a\x~az _ ( Pmin (Г)\ak o~ak   ф =    ^о^нРат       <*Г
\   Рп   )                            
\      Рн     /                                         2л^0Л         г

,      ^нРатРаЖпл (Р)    ф  dr                                                                                                                                                                          (176)

2А2                                  г3

188


где

, .       \р

^Упл (Р) = ~

1о |

При т]й = 0 или ak = ahQ из (176) можно получить уравнение установившейся фильтрации реального газа в упругодеформи-руемой среде, т. е. уравнение (172).

Также частным случаем из (176) получается уравнение, описывающее фильтрацию реального газа в коллекторах, под­верженных пластическим деформациям. При пластических де­формациях т|7г = оо или а/| = 0, тогда (176) перепишется в виде

(Pmin (Q/РнЛ ° «fa        ЕУнРат     j£_   ,    2нРатРаЖл (р)2 _^L    Л 77)

г                4я^                   г*


Здесь

^Пл (Р) = (Р/Р„)^ (P


min


Решение (177) получим ниже. Уравнение (176) можно ре­шить только приближенно, так же, как задачу об установив­шейся фильтрации газа в упругопластичной пористой среде при сохранении закона Дарси. Примем коэффициент изменения проницаемости ak постоянным для каждого цикла, но меняю­щимся от цикла к циклу. Коэффициент изменения проницае­мости аь является функцией минимального давления Ртт{г), которое различно по координате г. Мы же допускаем, что а>ь не зависит от г и его значение зависит от минимального забой­ного давления Pcmin, т. е. выражение (175) примет вид:

k (р, г) = k0 [pc т1пи)а* *~ak [p {г)1ря\\                                                                                       (178)

где

Решение (176) с учетом приближенного подхода запишется в виде:

Ра /               \ Ри У                               I    \     Рн

(179)

Здесь

(Рс/Рн)

ГУпл\Рс)

(о    ■ /о )l'27(ako~ak) Lo (Pc mm/Рн)

е = 1 -Yak — a^~az.

189