Материал: Математическое моделирование в процессах разработки и нефте-газодобычи

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

СПБГУАП| Институт 4 группа 4736

где

На решения.

 

 

25r

2

 

 

 

t

 

,

 

 

 

w

(1.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- коэффициент пьезопроводности.

 

рисунке изображено сравнение

 

логарифмической аппроксимации и точного

PD

Сравнение точного решения и логарифмической аппроксимации для модели постоянного дебита

10

9

8

7

6

5

4

3

 

 

 

Точное решение

 

 

 

 

 

 

 

 

2

 

 

 

Логарифмическая аппроксимация

 

1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

0.02

0.04

0.06

0.08

0.1

0.12

0.14

 

 

 

 

tDA

 

 

 

Рисунок 1.9. Сравнение точного решения и логарифмической аппроксимации для модели постоянного дебита

Для динамики забойного давления при вышеизложенном условии имеет место формула, которая получается из (1.32) путем подстановки r rw :

p

wf

t

 

 

p0

Q

 

 

4 t

 

ln

r 2

4 kh

 

 

 

 

w

 

 

 

 

 

 

 

 

,

(1.34)

1.4.2. Аппроксимация решения на временах неустановившегося режима для модели постоянного давления

Рассмотрим модель постоянного давления. В работе [21] говорится, что при анализе дебита скважины, работающей на неустановившемся режиме с постоянным забойным давлением с точностью около 1% можно пользоваться приближенной формулой:

18

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

 

2 kh

P P

 

q

0

wf

 

 

 

 

t

 

 

 

 

ln 1

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

w

Формула (1.35) имеет место при следующем условии на время:

 

r

r

 

 

t

 

 

2

,

e

w

 

 

 

 

 

 

 

 

 

(1.35)

(1.36)

где re - радиус дренирования, rw - радиус скважины.

Многие специалисты при анализе неустановившегося режима скважины,

работающей с постоянным давлением, пользуются соответствующей формулой (1.34) для неустановившегося режима скважины, работающей с постоянным дебитом, предполагая,

что давление в формуле (1.34) постоянно, а дебит при этом меняется. Иногда, построенное таким образом решение называется «перевертышем», отражая его суть. На Рисунок 1.10

показано сравнение точного решения, с аппроксимацией Чекалюка [21], а также с

«перевертышем». График построен в безразмерных координатах, поэтому выводы верны для всей области изменения параметров пласта. Из рисунка видно, что аппроксимация Чекалюка имеет более точное приближение. Заметим, что с математической точки зрения

«перевертышем» в модели постоянного давления пользоваться некорректно.

Сравнение решений

QD

0.17

0.16

Точное решение

0.15

Аппроксимация Чекалюка

"Перевертыш"

0.14

0.13

0.12

0

0.05

0.1

0.15

0.2

tDA

Рисунок 1.10. Сравнение точного решения, аппроксимации Чекалюка и «перевертыша»

19

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

Как уже было замечено ранее, за неустановившимся режимом следует поздний неустановившийся – переходный режим между неустановившимся и псевдоустановившимся или установившимся режимом (если псевдоустановившийся режим отсутствует). Пока отсутствуют аппроксимации этого режима, при его описании можно пользоваться лишь точным или численным решением.

1.4.3.Время протекания неустановившегося режима

Вразных источниках встречается разные данные о продолжительности

неустановившегося режима. Очевидно, что мерой характерного времени протекания неустановившегося режима в ограниченных пластах есть:

 

 

ttr C

A

,

(1.37)

 

 

 

 

 

 

 

 

где

A

- площадь пласта или

области дренирования,

- коэффициент

пьезопроводности, C - безразмерный коэффициент пропорциональности.

Проблема возникает в отыскании численного значения коэффициента C в формуле

(1.37). Рассмотрим скважину, работающую с условием постоянного дебита в центре кругового однородного изотропного пласта. Для общности построим на одном графике в безразмерных координатах решения зависимости забойного давления от времени для скважины, находящейся в бесконечном пласте, скважины, находящейся в пласте с поддержанием постоянного давления на границе, а также для скважины, находящейся в пласте с отсутствием перетока на границе. По оси абсцисс отложим безразмерное время,

нормированное на площадь пласта:

tDA

t

.

A

 

 

По оси ординат отложим безразмерный перепад давления:

(1.38)

p

D

 

2 kh pq

,

(1.39)

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

Смысл введения безразмерных координат заключается в общности построенных графиков в этих координатах для произвольного разброса параметров пласта, так как зависимости «нормируются» на различные характеристики пласта для получения безразмерного комплекса.

20

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

Время протекания неустановившегося режима в модели постоянного дебита

pD

12

 

 

 

 

 

10

 

 

 

 

 

8

 

 

 

 

 

6

 

 

 

 

 

4

 

 

Постоянное давление на границе

 

 

 

 

 

2

 

 

Отсутствие перетока через границу

 

 

 

 

 

 

 

 

 

Бесконечный пласт

 

 

0

 

 

 

 

 

0

0.1

0.2

0.3

0.4

0.5

 

 

 

tDA

 

 

Рисунок 1.11. Время протекания неустановившегося режима в модели постоянного дебита

Из Рисунок 1.11 видно, что решения совпадают на неустановившемся режиме.

Расхождение начинается при времени tDA 0.1, при этом времени решения для разных граничных условий отличаются не более 1%. Это определяет коэффициент пропорциональности в формуле (1.37). Поэтому для скважины в центре кругового пласта имеет место оценка для времени неустановившегося режима:

 

 

 

A

 

r

2

t

 

0.1

 

 

 

 

e

 

 

 

 

 

 

tr

 

 

 

 

 

 

 

 

 

r

2

 

 

0.3

e

 

 

,

(1.40)

. В книге [57] произведена оценка времени неустановившегося режима для скважин,

находящихся в пластах различной формы, а также в различных частях этого пласта,

показано, что формула (1.40) является верхней оценкой продолжительности неустановившегося режима в модели постоянного дебита.

Рассмотрим модель постоянного давления. Аналогично предыдущему пункту построим зависимости дебита от времени в безразмерных координатах для трех случаев:

бесконечного пласта, пласта с поддержанием постоянного давления на границе и пласта с отсутствием перетока на границе (Рисунок 1.12).

Безразмерный дебит определяется по формуле:

21

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

q

 

 

q B

 

D

2 kh p

,

 

 

 

 

 

 

(1.41)

Время протекания неустановившегося режима в модели постоянного давления

QD

0.16

0.14

0.12

0.1

0.08

0.06

0.04

Постоянное давление на границе

 

0.02

Отсутствие перетока через границу

 

 

Бесконечный пласт

0

 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

tDA

Рисунок 1.12. Время протекания неустановившегося режима в модели постоянного давления

Из графика видно, что расхождение начинается при tDA 0.1.

Расчеты показывают, что и в модели постоянного давления при безразмерном времени, меньшем чем 0,1: tDA 0.1, решения различаются также не более чем 1%.

Анализируя результат [20], [21] и наши расчеты, можно без ограничения общности дать достаточно точную оценку протекания неустановившегося режима следующей формулой:

t

tr

 

0.1

A

.

(1.42)

22

Контакты | https://new.guap.ru/i03/contacts