книги / Мониторинг гидравлического разрыва пласта на основе математической обработки геолого-промысловых данных
..pdfОкончание табл. 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,553− 0,00142t |
-0,82 |
|
|
|
|
39 |
a
б
Рис. 3.10. Зависимость Z-1 от t: а – в интервале 08.11.2009–11.04.2012; б – в интервале 08.11.2009–11.01.2012
40