Файл: Министерство науки и высшего образования рф федеральное государственное бюджетное учреждение высшего образования Российский государственный геологоразведочный университет имени Серго Орджоникидзе мгри.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.11.2023

Просмотров: 56

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Лабораторная работа №7 «Построение математической модели продуктивного пласта на основе метода множественной регрессии»

Пусть имеется два поля, два вида каротажа опорной скважины (F1; F2).

Введем функцию Y, отражающую гипотезу Н0 о продуктивном пласте вдоль эталонного объекта. Это значит, что:



Матрицы F и Y будет иметь следующий вид:

F1

F2

Y

10,373

6,76

0

4,373

9,373

0

13,373

7,373

0

15,373

-13,627

0

14,373

-10,627

0

13,373

-12,627

0

6,76

-7,627

0

6,76

-6,627

0

5,76

-6,627

0

11,76

-13,627

0

9,76

-6,627

0

7,76

0,76

0

5,76

-1,24

0

11,28

5,72

0

-6,28

6,72

0

-8,28

22,72

10

-16,28

11,72

10

-15,28

12,72

10

-14,28

14,72

10

-11,72

12,72

10

-14,28

10,72

10

-9,4

7,4

0

12,373

-0,24

0

18,373

3,76

0

10,373

-5,627

0

8,373

-10,627

0

9,373

-9,627

0

7,373

-11,627

0

15,373

-7,627

0

4,76

1,76

0

7,76

-0,24

0

9,76

2,76

0

0,76

2,76

0

-0,24

2,76

0

14,373

1,76

0

14,373

-10,62

0

13,373

-12,27

0

6,76

-7,67

0

2,73

-6,62

0

5,76

-6,67

0


Будем полагать, что между полем Y и полями F1 и F2 существует линейная корреляционная зависимость:



Отыщем уравнение, определяющее связь значений функции Y со значениями геологических полей F1 и F2 в каждой точке скважины.

В матричном виде эта зависимость примет вид:



где – заданные значения функции, а - вычисленные значения по уравнению регрессии у=FA.

Задача состоит в нахождении коэффициентов А,



Для решения системы уравнений наберем следующее выражение:



Для начала к матрице с переменными Fj добавляем единичный столбец. Получим матрицу :

1

1

1

1

1



1

1

1

1

10,373

4,373

13,373

15,373

14,373



13,373

6,76

2,73

5,76

6,76

9,373

7,373

-13,627

-10,627



-12,27

-7,67

-6,62

-6,67

Умножим матрицы (FT·F):

40

202,818

-23,735

202,818

4650,147

-2554,3

-23,735

-2554,3

3349,515

Умножим матрицы FT и Y:


60

-801,2

853,2





???????? ???? =

Находим обратную матрицу [???????? ????]−1:


0,036758

-0,00251

-0,00166

-0,00251

0,000542

0,000395

-0,00166

0,000395

0,000588




[???????? ????]−1 =


Далее и находим коэффициенты А по ранее известной формуле:

???? [???????? ????]−1 ???????? ????


0,036758

-0,00251

-0,00166

-0,00251

0,000542

0,000395

-0,00166

0,000395

0,000588


























































60

-801,2

853,2

2,806

-0,24753

0,085847



???? =


В результате получим вектор-столбец искомых коэффициентов А.

Полученные на эталонном объекте значения коэффициентов позволяют записать уравнение
регрессии,

???? = 2,806 - 0,24753 · ????1 + 0,085847 · ????2

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

x

y

1

0,818749

2

2,528219

3

0,128797

4

-2,16904

5

-1,66397

6

-1,58814

7

0,477979

8

0,563826

9

0,811351

10

-1,27473

11

-0,17875

12

0,950451

13

1,273808

14

0,504963

15

4,937354

16

6,805955

17

7,841842

18

7,680164

19

7,604332

20

6,798973

21

7,260945

22

5,768009

23

-0,27723

24

-1,41899

25

-0,24464

26

-0,17882

27

-0,3405

28

-0,01714

29

-1,65396

30

1,778874

31

0,864604

32

0,627094

33

2,854822

34

3,102348

35

-0,60059

36

-1,66337

37

-1,55749

38

0,474287

39

1,561954

40

0,80766