Автор работы: Пользователь скрыл имя, 10 Декабря 2011 в 18:55, курсовая работа
В данной работе рассчитывается динамика средневзвешенного по объему порового пространства пластового давления в газовой и газоконденсатной залежи методом последовательного приближения (метод итерации).
Он заключается в следующем:
1. В качестве первого приближения, уместно взять приближение с предыдущего временного шага, т.е. принимаем Рн = Рν.
2. По формулам приведённым ниже рассчитываем Р*, z(Рн*), z (Рν*).
3. Рассчитываем Рν+1 и находим разницу двух давлений |Рν+1-Pν|
ВВЕДЕНИЕ 4
1 ТЕОРИЯ РАСЧЕТА ДИНАМИКИ СРЕДНЕВЗВЕШЕННОГО ДАВЛЕНИЯ 5
1.1 Газовая залежь 5
2.2 Газоконденсатная залежь 5
2 РАСЧЕТ ДИНАМИКИ СРЕДНЕВЗВЕШЕННОГО ДАВЛЕНИЯ 7
2.1 Исходные данные 7
2.2 Газовая залежь 7
2.3 Газоконденсатная залежь 10
2.4 Распределение средневзвешенного давления 12
ЗАКЛЮЧЕНИЕ 13
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 14
СОДЕРЖАНИЕ
ВВЕДЕНИЕ 4
1 ТЕОРИЯ РАСЧЕТА ДИНАМИКИ СРЕДНЕВЗВЕШЕННОГО ДАВЛЕНИЯ 5
1.1 Газовая залежь 5
2.2 Газоконденсатная залежь 5
2 РАСЧЕТ ДИНАМИКИ СРЕДНЕВЗВЕШЕННОГО ДАВЛЕНИЯ 7
2.1 Исходные данные 7
2.2 Газовая залежь 7
2.3 Газоконденсатная залежь 10
2.4 Распределение средневзвешенного давления 12
ЗАКЛЮЧЕНИЕ 13
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 14
В данной работе рассчитывается динамика средневзвешенного по объему порового пространства пластового давления в газовой и газоконденсатной залежи методом последовательного приближения (метод итерации).
Он заключается в следующем:
1. В качестве первого приближения, уместно взять приближение с предыдущего временного шага, т.е. принимаем Рн = Рν.
2. По
формулам приведённым ниже
3. Рассчитываем Рν+1 и находим разницу двух давлений |Рν+1-Pν|
4. Итерационный процесс на каждом временном шаге следует выполнять до достижения точности 0,1 атм., то есть когда |Рν+1-Pν|<0,1.
5. В качестве второго приближения, принимаем приближение с предыдущего временного шага, т.е. Рν=Рν+1. Дальнейшие расчёты аналогичны. Расчёт считается законченным, если достигается точность 0,1 атм.
Приведенная
расчетная схема в силу своего
учебного характера не может полностью
корректно применяться для
Пусть
залежь разрабатывается с постоянным
во времени отбором
. При известном на момент времени
среднем пластовом давлении
текущие запасы газа в залежи определяются
по формуле:
где - объем порового пространства;
- коэффициент сверхсжимаемости.
Давление
в данной формуле берется в
атмосферах. Начальные запасы газа
находятся по аналогичной формуле
от начального давления
. Поскольку текущие и отобранные запасы
в каждый момент времени в сумме равны
начальным, то справедливо соотношение
являющееся для любого уравнением относительно искомого давления.
Указанное
уравнение для фиксированного
можно решить с помощью следующей
итерационной процедуры:
В
целом задача аналогична предыдущей
за исключением того, что по причине
выпадения в пласте конденсата поровый
объем, занятый газовой фазой, становится
переменной величиной. Данный объем
корректируется величиной пластовых
потерь конденсата и определяется долей
выпавшей углеводородной жидкости, являющейся
функцией давления. В рассматриваемом
случае текущие запасы газа в залежи определяются
по формуле:
где - объемная доля жидкой фазы.
Соотношение
между начальными, текущими и отобранными
запасами принимает вид
,
(при
начальном давлении жидкая
Залежь
разрабатывается с постоянным во
времени дебитом Q. Начальное пластовое
давление равно
. Объём порового пространства
. Функциональные зависимости
и
от давления заданы в виде полиномов
от безразмерного давления
и имеют следующий вид
и
. Исходные
данные возьмем из таблицы 2.1
Таблица 2.1 – исходные данные
Вариант | Q, 106 м3/год | ||||
79 |
220 |
1,5 |
29,7 |
a0=1
a1= -1,1 a2=1,1 |
b0=0,07
b1=0,28 b2= -0,35 |
В
качестве первого приближения возьмем
давление с предыдущего временного шага,
т.е.
. Рассчитаем безразмерное давление
:
Найдем значение
полинома
:
Вычислим значение
давления
, t = 1 год:
Проверим условие нахождения искомого давления :
условие не выполнено.
Тогда вычислим
Проверим
условие
не выполнено. Итерационный процесс на
каждом временном шаге выполнять до достижения
точности 0,1 атм, то есть когда
. Дальнейшие вычисления давления
на 1 год сведем в таблицу 2.2.
Таблица 2.2 – Результат расчета
220 | 1,00000 | 1,00000 | 200,2000 | 19,8000 |
200,2 | 0,91000 | 0,90991 | 182,1640 | 18,0360 |
182,164 | 0,82802 | 0,84336 | 168,8398 | 13,3242 |
168,8398 | 0,76745 | 0,80368 | 160,8976 | 7,9421 |
160,8976 | 0,73135 | 0,78388 | 156,9321 | 3,9655 |
156,9321 | 0,71333 | 0,77506 | 155,1669 | 1,7652 |
155,1669 | 0,70530 | 0,77136 | 154,4272 | 0,7397 |
154,4272 | 0,70194 | 0,76986 | 154,1257 | 0,3015 |
154,1257 | 0,70057 | 0,76925 | 154,0042 | 0,1215 |
154,0042 | 0,70002 | 0,76901 | 153,9555 | 0,0487 |
Имеем давление на конец 1-го года равное 154,0042 атм.
Аналогично
рассчитаем давления на пятилетний период
и сведем в таблицу 2.3.
Таблица 2.3 – Расчет средневзвешенного давления в газовой залежи
1 | 2 | 3 | 4 | 5 |
220 | 1,00000 | 1,00000 | 200,2000 | 19,8000 |
200,2 | 0,91000 | 0,90991 | 182,1640 | 18,0360 |
182,164 | 0,82802 | 0,84336 | 168,8398 | 13,3242 |
168,8398 | 0,76745 | 0,80368 | 160,8976 | 7,9421 |
160,8976 | 0,73135 | 0,78388 | 156,9321 | 3,9655 |
156,9321 | 0,71333 | 0,77506 | 155,1669 | 1,7652 |
155,1669 | 0,70530 | 0,77136 | 154,4272 | 0,7397 |
154,4272 | 0,70194 | 0,76986 | 154,1257 | 0,3015 |
154,1257 | 0,70057 | 0,76925 | 154,0042 | 0,1215 |
154,0042 | 0,70002 | 0,76901 | 153,9555 | 0,0487 |
| ||||
154,1257 | 0,70057 | 0,76925 | 138,7730 | 15,3527 |
138,773 | 0,63079 | 0,74382 | 134,1843 | 4,5887 |
134,1843 | 0,60993 | 0,73829 | 133,1880 | 0,9963 |
133,188 | 0,60540 | 0,73722 | 132,9945 | 0,1935 |
132,9945 | 0,60452 | 0,73702 | 132,9579 | 0,0366 |
132,9579 | 0,60435 | 0,73698 | 132,9510 | 0,0069 |
132,951 | 0,60432 | 0,73697 | 132,9497 | 0,0013 |
132,9497 | 0,60432 | 0,73697 | 132,9494 | 0,0002 |
132,9494 | 0,60432 | 0,73697 | 132,9494 | 0,0000 |
132,9494 | 0,60432 | 0,73697 | 132,9494 | 0,0000 |
| ||||
133,1880 | 0,60540 | 0,73722 | 118,3975 | 14,7905 |
118,3975 | 0,53817 | 0,72660 | 116,6924 | 1,7052 |
116,6924 | 0,53042 | 0,72602 | 116,5985 | 0,0939 |
116,5985 | 0,52999 | 0,72599 | 116,5939 | 0,0046 |
116,5939 | 0,52997 | 0,72599 | 116,5937 | 0,0002 |
116,5937 | 0,52997 | 0,72599 | 116,5937 | 0,0000 |
116,5937 | 0,52997 | 0,72599 | 116,5937 | 0,0000 |
116,5937 | 0,52997 | 0,72599 | 116,5937 | 0,0000 |
116,5937 | 0,52997 | 0,72599 | 116,5937 | 0,0000 |
116,5937 | 0,52997 | 0,72599 | 116,5937 | 0,0000 |
| ||||
116,5939 | 0,52997 | 0,72599 | 102,2191 | 14,3748 |
102,2191 | 0,46463 | 0,72638 | 102,2737 | 0,0546 |
102,2737 | 0,46488 | 0,72636 | 102,2710 | 0,0027 |
102,271 | 0,46487 | 0,72636 | 102,2712 | 0,0001 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
102,2712 | 0,46487 | 0,72636 | 102,2712 | 0,0000 |
| ||||
102,2712 | 0,46487 | 0,72636 | 87,8893 | 14,3819 |
87,88927 | 0,39950 | 0,73611 | 89,0694 | 1,1802 |
89,06943 | 0,40486 | 0,73496 | 88,9297 | 0,1397 |
88,92974 | 0,40423 | 0,73509 | 88,9459 | 0,0161 |
88,94588 | 0,40430 | 0,73507 | 88,9440 | 0,0019 |
88,94401 | 0,40429 | 0,73508 | 88,9442 | 0,0002 |
88,94423 | 0,40429 | 0,73508 | 88,9442 | 0,0000 |
88,9442 | 0,40429 | 0,73508 | 88,9442 | 0,0000 |
88,9442 | 0,40429 | 0,73508 | 88,9442 | 0,0000 |
88,9442 | 0,40429 | 0,73508 | 88,9442 | 0,0000 |
Информация о работе Расчет динамики средневзвешенного давления