Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Мониторинг гидравлического разрыва пласта на основе математической обработки геолого-промысловых данных

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
19.49 Mб
Скачать

Окончание табл. 3.3

Сква-

Статистические характеристики показателей

Статистические

жина №

Класс 1,

Класс 2,

характеристики

n = 69

n = 81

показателей

 

 

 

 

3,1859

206

17,9 ± 6,5

13,7 ± 9,2

0,0017

12,3335

 

 

 

 

 

 

0,0020

 

 

 

-6,5112

207

29,6 ± 3,4

33,4 ± 3,6

0,0000

37,8521

 

 

 

 

 

 

0,0000

 

 

 

1,5660

212

17,5 ± 13,9

14,0 ± 13,0

0,1194

2,4564

 

 

 

 

 

 

0,2928

 

 

 

-8,5039

216

15,0 ± 10,9

25,6 ± 2,4

0,0000

61,5640

 

 

 

 

 

 

0,0000

 

 

 

-0,7117

217

10,3 ± 10,8

11,2 ± 3,1

0,4777

0,5313

 

 

 

 

 

 

0,7666

 

 

 

6,6508

231

18,1 ± 11,7

9,6 ± 3,2

0,0000

45,1903

 

 

 

 

 

 

0,0000

 

 

 

-6,5101

234

46,1 ± 8,9

53,9 ± 5,9

0,0000

39,2385

 

 

 

 

 

 

0,0000

По скважинам № 214, 215, 68, 222, 230, 229, 228, 64, 213 сред-

ние значения QН до и после ГРП статистически различаются. Также статистические различия получены по скважинам № 207, 65, 203, 216, 102, но статистического различия в средних значениях не получено по скважине № 217. По скважинам № 204, 206, 212, 227, 235, 236, 231 средние значения QН также отличаются, но только в обрат-

31

ную сторону: после проведения ГРП значения QН весьма уменьшились. Для скважины № 57 статистическое различие в средних значениях также зафиксировано. Здесь необходимо отметить, что различия наблюдается при невысоких дебитах нефти.

Для более полного статистического обоснования полученных данных выполнено сравнение распределений по QН до и после ГРП по скважинам. Примеры распределений QН по скважинам № 65 и № 217, которые имеют максимальные и минимальные значения критерия t, представлены в табл. 3.4.

Таблица 3.4

Распределение значений QН

 

 

 

Интервалы варьирования QН, т/сут

 

 

Класс

0–5

5–10

10–15

15–20

20–25

25–30

30–35

35–40

40–45

45–50

 

 

 

 

 

Скважи

на № 6

5

 

 

 

1

0,073

0,319

0,376

0,132

0,028

0,072

2

0,197

0,345

0,444

0,014

 

 

 

 

 

Скважи

на № 217

 

 

 

1

0,289

0,377

0,145

0,016

0,173

2

 

0,444

0,456

0,086

0,014

Анализ показывает, что наблюдаются различия в распределениях значений QН для изучаемых классов для скважины № 65 и значительно меньшие различия для скважины № 217.

Для оценки статистического различия плотностей распределений QН применялась статистика (критерий) Пирсона χ2 , которая рассчитывается по следующей формуле:

e

 

 

 

 

2

1

 

M1

M 2 ,

χ2 = N1N2

 

 

i=1

M1 + M 2

N1

N2

 

где N1, N2 – соответственно количество значений QН в классах 1, 2; M1 ,M 2 – количество значений, попавших в заданный интервал, соответственнодля двух изучаемыхклассов; е– количество интервалов.

Значения критерия χ 2 приведены в табл. 3.4, из которой видно, что по данному параметру распределения для классов 1 и 2 по

32

QН статистически различаются для скважины № 65 и не различаются для скважины № 217. Значения χ 2 вычислены по всем изучаемым скважинам и подтвердили выводы, полученные по критерию t. Поле корреляции между t и χ 2 приведено на рис. 3.6, из которого видно,

что значения t и χ 2 хорошо коррелируют между собой. Для более

полного представления о том, как изменяются значения QН во времени, построим графики зависимости с учетом времени проведения ГРП в скважине № 221. Для шести скважин, наиболее близко расположенных от скважины № 221, зависимости приведены на рис. 3.6.

χ2

Рис. 3.5. Поле корреляции

Данные рис. 3.6 свидетельствуют, что зависимости между QН и t до и после ГРП отличаются. Изменения значений QН во времени до проведения ГРП в основном достаточно сложные, а после ГРП – повышаются. Уравнения регрессии QН от t по скважинам до и после проведения ГРП приведены в табл. 3.5.

Значения r для зависимостей между QН и t до проведения ГРП в скважине № 221 изменяются от «-0,14» до «0,65», после – от «0,56» до «0,84». Отметим, что значения r во втором случае для всех скважин больше, чем в первом. Среднее значение r до проведения ГРП равно « +0,176 ± 0,330 », после проведения ГРП « +0,708 ± 0,095 »,

по критерию t они являются статистически различными

33

(t = −3,79048; p = 0,003541). Совместный график изменения значений QН и t по всем шести скважинам представлен на рис. 3.7.

Скв. № 64 до ГРП

 

Скв. № 64 после ГРП

 

Скв. № 214 до ГРП

 

Скв. № 214 после ГРП

 

 

 

 

 

 

 

Скв. № 215 до ГРП

 

Скв. № 215 после ГРП

 

Скв. № 222 до ГРП

 

Скв. № 222 после ГРП

 

 

 

 

 

 

 

 

Скв. № 228 до ГРП

 

Скв. № 222 после ГРП

 

Скв. № 229 до ГРП

 

Скв. № 229 после ГРП

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.6. Зависимости изменения QН во времени до и после ГРП

 

 

 

 

 

 

 

 

 

 

Таблица 3. 5

 

 

 

Зависимость QН от t

 

 

Сква-

Уравнение регрессии

r

Уравнение регрессии

 

r

жина

до ГРП

после ГРП

 

 

 

 

 

 

 

 

 

 

 

QН = −409,189 + 0,01066t

 

 

 

64

QН = −591,293 + 0,0153t

0,65

 

 

0,75

214

QН = 57,227 0,00109t

-0,05

 

QН = −126,925 + 0,00354t

 

0,75

215

QН = −269,029 + 0,0074t

0,32

 

QН = −137,202 + 0,00400t

 

0,69

222

QН = −368,903 + 0,0101t

0,41

 

QН = −312,882 + 0,00863t

 

0,84

228

QН = −89,318 0,00148t

-0,13

 

QН = −80,133 + 0,00291t

 

0,56

229

QН = 168,605 0,00306t

-0,14

 

QН = −158,307 + 0,00505t

 

0,66

Все

QН = −152,345 + 0,0045t

0,16

 

QН = −204,106 + 0,0058t

 

0,33

скв.

 

 

 

 

 

 

 

 

 

 

 

34

 

 

 

 

 

 

 

 

 

 

 

До ГРП После ГРП

Рис. 3.7. Зависимости изменения QН во времени до и после ГРП для скважин, минимально удаленных от скважины № 221

Можно отметить, что зависимости между QН и t до и после ГРП отличаются друг от друга. Изменение данных QН во времени до проведения ГРП характеризуются значительным диапазоном значений – 0 до 65 т/сут, и распределение значений во времени имеет очень сложный характер. Значения QН во времени после проведения ГРП в основном повышаются. Отметим, что во временном интервале от 08.10.2009 г. до 11.01.2012 г., в который входит время проведения ГРП, наблюдается особенность, заключающаяся в том, что соотношения между QН и t значительно отличаются от таковых до этого временного диапазона и после него. Сравнение уравнений (см. табл. 3.5) показывает, что во втором случае наблюдается не только повышение значений r практически в два раза, но и свободный член уравнения, и коэффициент при t также больше. Это свидетельствует, что действительно после проведения ГРП в скважине № 221 произошло изменение работы системы добычи нефти в центральной части объекта разработки, при этом необходимо отметить, что во времени эти изменения продолжаются.

35

Докажем, что данные изменения не являются случайными. Для этого будет использован пошаговый линейный дискриминантный анализ. Возможность построения линейной дискриминантной

функции (ЛДФ) сводится к следующему: если обозначить через Xij

значения дебитов нефти во времени по всем изучаемым скважинам до проведения ГРП в скважине № 221, то в результате может быть получена матрица W1 порядка результатов m и n1 наблюдений над этой выборкой:

 

X11 X12 ….X1n1

 

 

 

X21

X22 ….X2n1

 

W

=

.

1

 

 

 

 

 

 

 

 

Xm2 ….Xmn1

 

 

Xm1

 

Обозначим через Xij1 значения дебитов нефти во времени по

всем изучаемым скважинам после проведения ГРП, получим матрицу W2 порядка m и n2:

 

 

X111 X121 .X1n1

 

 

 

 

1

1

1

 

 

W =

X 21

X 22

.X 2n2

 

,

2

 

 

 

 

 

 

 

 

 

 

 

 

 

X 1

.X 1

 

 

 

X 1

 

 

 

 

m1

m2

mn2

 

 

где m – число скважин; n1, n2 – количество значений дебитов нефти до проведения ГРП и после такового.

Для построения линейной дискриминантной функции составляются матрицы центрированных сумм квадратов и смешанных произведений, по ним вычисляется выборочная матрица. Далее для определения коэффициентов линейной дискриминантной функции находится обратная выборочная ковариационная матрица – матрица С. Затем вычисляются граничные значения дискриминантных функций

(R0 ), что делит выборку на три множества. Надежность классификации вычисляются с помощью критерия Пирсона χ2 . В результате реализации данного метода для шести скважин, находящихся наи-

36

более близко к скважине № 221 (условной первый ряд), получена следующая линейная дискриминантная функция:

Z-1 = −0,119622Qскв.228 0,083035Qскв.222

0,024023Qскв.214

+

н

н

 

 

н

 

+ 0,055011Qскв.215 0,016267Qскв.64

+ 7,243166

 

н

н

 

 

 

 

при R = 0,844; χ2 = 181,98; p = 0,000,

где Qскв.228

– индивидуальные

 

 

н

 

 

 

значения QН по каждой скважине.

Отметим, что в построение данной функции не включали значения QН по скважине № 229, которая наиболее удалена от скважины № 221. В построении линейной дискриминантной функции приняли участие 83,3 % скважин. По данной функции вычислены значения Z-1 и определены вероятности к классу после ГРП – Р(Z-1). Зависимость Р(Z-1) от Z приведена на рис. 3.8.

P(Z-1)

Z-1

Рис. 3.8. Зависимость Р(Z-1) от Z-1

При изменении значений Z-1 от отрицательных к положительным величина Р(Z-1) закономерно снижается. Среднее значение Z-1 для первого класса равно «2,153», для второго «-2,894». Доля верного распознавания составила 97,3 %. Изменение значений Z-1 во времени приведено на рис. 3.9.

37

Анализ графика показывает, что во времени до проведения ГРП наблюдается значительное изменение значений Z-1, особенно в левой его части. Здесь отметим, что в диапазоне времени от 07.06.2009 г. до даты проведения ГРП изменение значений Z-1 от t имеет более выраженный направленный вид, чем до этого времени.

Z-1

До ГРП После ГРП

Рис. 3.9. Зависимость Z-1 от t

Зависимость Z-1 от t после проведения ГРП характеризуется закономерным повышением значений Z-1 во времени, но в пределах данного графика также имеется два вида соотношений между Z-1 от t. Уравнения регрессии, описывающие данные соотношения только из условия до и после ГРП (без учета вышеотмеченных особенностей), приведены в табл. 3.6.

Сравнение этих уравнений показывает, что зависимость, построенная по значениям QН до ГРП, характеризуется сравнительно более низким значением коэффициента r, а также значениями свободного члена и коэффициента при t в уравнении регрессии по сравнению с этими же значениями в уравнении, построенным по данным QН после ГРП.

38

 

 

 

 

 

 

Таблица 3. 6

 

 

Зависимость Z-1 от t

 

 

 

 

 

 

Уравнение регрессии

 

Пара-

Уравнение регрессии

r

 

r

метр

до ГРП

 

 

 

после ГРП

 

 

 

 

 

Первый

Z-1 = 20,634 0,000476t

 

 

 

Z-1 = 44,966 0,00110t

 

ряд

-0,27

 

 

-0,83

скважин

 

 

 

 

 

 

 

Это, вероятно, может свидетельствовать о том, что выполненное ГРП в скважине №221 практически полностью перестроило процесс формирования значений QН как по отдельным скважинам, так и по их совокупности, наиболее близко расположенных к скважине №221. Отметим, что данные по скважине №221 при построении ЛДФ не использовались. Одновременно учитываем, что во временном интервале от 08.01.2009 до 11.01.2012 г., в который входит осуществление ГРП, наблюдается особенность, заключающаяся в очень сильной корреляции, которуюможно описатьследующим видом:

Z -1 = 141,667 0,00347 t, r = −0,95.

Графическое изображение соотношений до и после данного временного диапазона приведено на рис. 3.9.

В пределах корреляционного поля в данном временном диапазоне выделяется группа точек, имеющих другое соотношение между Z-1 и t, чем на других частях корреляционного поля.

Данное подполе разделено на две части – до ГРП и после него, построены уравнения регрессии и вычислены коэффициенты r. Уравнениярегрессии приведены втабл. 3.7, полякорреляции– нарис. 3.10.

 

 

 

Таблица 3. 7

Зависимость Z-1 от t в интервале времени до и после ГРП

 

 

 

 

 

 

Группа

Уравнение регрессии

r

Уравнение регрессии

r

значе-

до ГРП

после ГРП

ний

 

 

 

 

 

 

Левая

Z-1 = 0,531+ 0,0000613t

0,02

Z-1 = 104,362 0,00256t

-0,53

Правая

Z-1 = −124,242 0,00304t

-0,94

Z-1 = 58,5530,00142t

-0,82

 

 

 

 

39

a

б

Рис. 3.10. Зависимость Z-1 от t: а – в интервале 08.11.2009–11.04.2012; б – в интервале 08.11.2009–11.01.2012

40