Проблемы эффективности производства на северных нефтегазодобывающих предприятиях (доклады и сообщения). Том 2, страница 25

Расчеты проведены путем поиска аппроксимирующего пол­инома с помощью метода наименьших квадратов. В общем виде метод может быть описан следующим образом. Пусть в точках 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