Автор: Пользователь скрыл имя, 28 Марта 2012 в 09:33, курсовая работа
Курсовая работа по дисциплине «Численные методы решения задач нефтегазопромысловой механики» содержит:
1 Расчет динамики средневзвешенного по объему порового пространства пластового давления в газовой залежи.
2 Расчет динамики средневзвешенного по объему порового пространства пластового давления в газоконденсатной залежи.
Динамика пластового давления определяется запасами газа, его свойствами и темпами отбора.
Для каждой задачи расчеты выполняются с временным шагом 1 год для 5-го периода.
ВВЕДЕНИЕ………………………………………………………........................ 6
ИСХОДНЫЕ ДАННЫЕ……………………………………………………........ 7
1 Расчет динамики средневзвешенного по объему порового
пространства пластового давления в газовой залежи. ……………...……… 8
2 Расчет динамики средневзвешенного по объему порового
пространства пластового давления в газоконденсатной залежи. .………… 12
ЗАКЛЮЧЕНИЕ…………………………………………………………………. 16
БИБЛИОГРАФИЧЕСКИЙ СПИСОК………………………………..………… 17
Таблица 1.5 – Результаты вычислений для 5-го года разработки
5 год | ||||
Точность | ||||
160,1867 | 0,533956 | 0,800922 | 140,1614 | 20,02528 |
140,1614 | 0,467205 | 0,80086 | 140,1506 | 0,010825 |
Окончательный результат представлен в таблице 1.6
Таблица 1.6 – Динамика пластового давления за 5 лет разработки
Год разработки | Давление, атм | Давление, МПа |
1 | 239,7129 | 23,97129 |
2 | 207,28 | 20,728 |
3 | 182,0653 | 18,20653 |
4 | 160,1867 | 16,01867 |
5 | 140,1506 | 14,01506 |
На рисунке 1.1 представлен график зависимости пластового давления во времени.
Рисунок 1.1 – Динамика пластового давления во времени
2 Расчет динамики средневзвешенного по объему
порового пространства пластового давления в
газоконденсатной залежи.
В целом задача аналогична предыдущей за исключением того, что по причине выпадения в пласте конденсата поровый объем, занятый газовой фазой, становится переменной величиной. Данный объем корректируется величиной пластовых потерь конденсата и определяется долей выпавшей углеводородной жидкости, являющейся функцией давления. В рассматриваемом случае текущие запасы газа в залежи определяются по формуле:
,
где – объемная доля жидкой фазы.
Соотношение между начальными, текущими и отобранными запасами принимает вид
,
(при начальном давлении жидкая фаза отсутствует), откуда следует итерационная схема для расчета динамики пластового давления:
Решение
Считаем для 1-го года разработки.
атм.
атм.
Продолжаем считать до точности 0,1. Остальные значения считаются по аналогии и представлены в таблицах 2.1 – 2.5
Таблица 2.1 – Результаты вычислений для 1-го года разработки
1 год | ||||||
φ | z/1-φ | Точность | ||||
300 | 1 | 1 | 0 | 1 | 275 | 25 |
275 | 0,916667 | 0,938889 | 0,037222 | 0,975188 | 268,1766 | 6,823428 |
268,1766 | 0,893922 | 0,92414 | 0,046416 | 0,969123 | 266,5088 | 1,667816 |
266,5088 | 0,888363 | 0,92066 | 0,048601 | 0,967691 | 266,115 | 0,393784 |
266,115 | 0,88705 | 0,919846 | 0,049113 | 0,967356 | 266,0229 | 0,092147 |
Таблица 2.2 – Результаты вычислений для 2-го года разработки
2 год | ||||||
φ | z/1-φ | Точность | ||||
266,0229 | 0,886743 | 0,919656 | 0,049233 | 0,967278 | 241,8194 | 24,2035 |
241,8194 | 0,806065 | 0,87494 | 0,078045 | 0,949005 | 237,2513 | 4,568103 |
237,2513 | 0,790838 | 0,867669 | 0,082898 | 0,946099 | 236,5248 | 0,726468 |
236,5248 | 0,788416 | 0,866547 | 0,083653 | 0,945654 | 236,4135 | 0,111295 |
236,4135 | 0,788045 | 0,866376 | 0,083768 | 0,945586 | 236,3966 | 0,016947 |
Таблица 2.3 – Результаты вычислений для 3-го года разработки
3 год | ||||||
φ | z/1-φ | Точность | ||||
236,3966 | 0,787989 | 0,86635 | 0,083786 | 0,945576 | 212,7546 | 23,64202 |
212,7546 | 0,709182 | 0,835006 | 0,105763 | 0,933763 | 210,0966 | 2,657954 |
210,0966 | 0,700322 | 0,832103 | 0,107923 | 0,93277 | 209,8733 | 0,22331 |
209,8733 | 0,699578 | 0,831865 | 0,108101 | 0,93269 | 209,8552 | 0,018053 |