Расчеты проведены путем поиска аппроксимирующего полинома с помощью метода наименьших квадратов. В общем виде метод может быть описан следующим образом. Пусть в точках XI, Х2, . . . , хп заданы (точно или приближенно) значения уь У2,..... уп некоторой функции у •= f (x). Функция f (х) аппроксимируется
т полиномом Рщ(х) ш 2 at• х/ степени ш, коэффициенты которого
1=0 определяются путем минимизации функционала
п
F(a0, ai, . . . , am) = 2 (yi - Pm(x0)
i=\
по переменным a0, ai, . . . , am.
Прогноз падения пластового давления составлялся отдельно по каждому профилю для обеих частей газогидродинамической системы
52
Рис. 1. Карта изобар водонапорного бассейна на уровне начального ГВК (-1190 м):
1 - контур газоносности; 2 - изобары; 3 - в числителе - номер скважины,
в знаменателе - величина пластового давления
53
р,МПа |
a a юм |
8 9 Ю 11 |
Рис.2. Эпюры пластового давления в газовой залежи и подстилающем водоносном бассейне (продольный профиль I)
в каждой конкретной точке (скважине). При этом в качестве аргумента функции f берется номер года, а в качестве ее значения - соответствующее давление, т.е. xi = i (i = О, 1, ...,4), a yi - это давление в скважине в (1984+i) году. Степень искомого полинома m полагалась равной 1, т.е. предполагается, что давление в отдельной скважине аппроксимируется зависимостью вида Зо + ai х i, а минимизируемый функционал принимает вид
4 F(ao, ai) - 2 (Pi -ao-ai x i) ,
i=0
где Pi - измеренное давление газа в скважине в соответствующем году.
54
S 6 7 8 9 |
р.МПа
Профиль 11
Профиль Ж
Рис. 3. Эпюры пластового давления в газовой залежи и подстилающем водоносном бассейне (поперечные профили П, Ш): А - западная часть водоносного бассейна; В - границы газовой залежи; С - восточная
часть водоносного бассейна;
линия давления: 1 - в газовой залежи; 2 - в водоносной части; 3' и 3" - прогнозные значения пластового давления в газовой и водоносной частях
Минимизация проводится с помощью методов покоординатного спуска и золотого сечения. Полученные значения коэффициентов ао и ai используются для прогноза давления в скважине на год вперед. Прогнозируемое давление Рпр в скважине в 1989 г. вычислялось как значение полинома Pi(x) в точке Xi = 5, т.е. РПр - а0 + 5 х ai, где ао, ai - вычисленные ранее коэффициенты.
В табл. 1 даны значения коэффициентов а0, ai аппроксимирующего полинома Pi(x), прогнозируемое давление Рпр, а также значение функционала F.
55
Таблица 1
Прогноз падения пластового давления в системе: газовая залежь - водоносный бассейн (Уренгойское месторождение)
Коэффициенты |
Прогнози- |
||||
Значение |
полинома |
руемое |
|||
Профиль |
Скважина |
функцио- |
давление |
||
нала F |
а0 |
al |
(1989 г.) |
||
Рпр, МПа |
|||||
1 |
0.194 |
12.01 |
-0.547 |
9.28 |
|
2 |
0.0457 |
11.39 |
-0.533 |
8.72 |
|
3 |
0.011 |
10.52 |
-0.49 |
8.07 |
|
4 |
0.007 |
10.34 |
-0.45 |
8.09 |
|
1 |
5 |
0.004 |
10.21 |
-0.417 |
8.13 |
газовая |
6 |
0.008 |
10.73 |
-0.503 |
8.21 |
залежь |
7 |
0.05 |
11.15 |
-0.507 |
8.52 |
8 |
0.0139 |
11.91 |
-0.627 |
8.73 |
|
9 |
0.009 |
12.10 |
-0.61 |
9.05 |
|
10 |
0.003 |
12.33 |
-0.553 |
9.56 |
|
11 |
0.010 |
12.50 |
-0.403 |
10.44 |
|
1 |
0 |
13.1 |
-0.6 |
10.01 |
|
2 |
0.0082 |
11.92 |
-0.4 |
9.92 |
|
3 |
0.003 |
10.62 |
-0.33 |
8.97 |
|
4 |
0.014 |
10.37 |
-0.35 |
8.64 |
|
1 |
5 |
0.016 |
10.67 |
-0.38 |
8.75 |
водонос- |
6 |
0.0045 |
11.51 |
-0.45 |
9.24 |
ный пласт |
7 |
0.023 |
11.88 |
-0.39 |
9.93 |
8 |
0.012 |
12.15 |
-0.36 |
10.33 |
|
9 |
0.043 |
12.21 |
-0.32 |
10.62 |
|
10 |
0 |
12.7 |
-0.3 |
11.2 |
|
11 |
0 |
13.4 |
-0.3 |
11.9 |
|
5 |
0 |
10.3 |
-0.4 |
8.3 |
|
6 |
0 |
10.2 |
-0.4 |
8.2 |
|
газовая |
7 |
0.017 |
10.5 |
-0.34 |
8.8 |
залежь |
8 |
0.051 |
10,81 |
-0.38 |
8.89 |
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.