Расчет динамики средневзвешенного по объему порового пространства пластового давления в газовой залежи

СОДЕРЖАНИЕ

ВВЕДЕНИЕ………………………………………………………........................              6

ИСХОДНЫЕ ДАННЫЕ……………………………………………………........              7

1 Расчет динамики средневзвешенного по объему порового

   пространства пластового давления в газовой залежи. ……………...………   8

2 Расчет динамики средневзвешенного по объему порового

   пространства пластового давления в газоконденсатной залежи. .…………              12

ЗАКЛЮЧЕНИЕ………………………………………………………………….              16

БИБЛИОГРАФИЧЕСКИЙ СПИСОК………………………………..…………              17

ВВЕДЕНИЕ

Курсовая работа по дисциплине «Численные методы решения задач нефтегазопромысловой механики» содержит:

1 Расчет динамики средневзвешенного по объему порового пространства пластового       давления в газовой залежи.

2 Расчет динамики средневзвешенного по объему порового пространства пластового       давления в газоконденсатной залежи.

Динамика пластового давления определяется запасами газа, его свойствами и темпами отбора.

Для каждой задачи расчеты выполняются с временным шагом 1 год для 5-го периода.

ИСХОДНЫЕ ДАННЫЕ

Таблица 1 – Исходные данные

1 Расчет динамики средневзвешенного по объему

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

                            газовой залежи.

Пусть залежь разрабатывается с постоянным во времени отбором . При известном на момент времени среднем пластовом давлении текущие запасы в залежи определяются по формуле:

, где

- объем порового пространства;

- коэффициент сверхсжимаемости.

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

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

,

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

Функциональные зависимости и от давления заданы в виде полиномов от безразмерного давления . Это означает что, .

Итерационный процесс на каждом временном шаге следует выполнять до достижения точности 0,1 атм, т.е. когда .

Решение

Считаем для 1-го года разработки.

атм.

атм.

Продолжаем считать до точности 0,1. Остальные значения считаются по аналогии и представлены в таблицах 1.1 – 1.5

Таблица 1.1 – Результаты вычислений для 1-го года разработки

Таблица 1.2 – Результаты вычислений для 2-го года разработки

Таблица 1.3 – Результаты вычислений для 3-го года разработки

Таблица 1.4 – Результаты вычислений для 4-го года разработки