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

Повышение продуктивности и реанимация скважин

.pdf
Скачиваний:
43
Добавлен:
12.08.2019
Размер:
10.67 Mб
Скачать

vk.com/club152685050 | vk.com/id446425943

ГЛАВА

11

ПРОГНОЗИРОВАНИЕ ЭФФЕКТИВНОСТИ СКВАЖИННЫХ ОБРАБОТОК

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

Технологический эффект виброволнового воздействия определяется совокупным влиянием множества разнородных геолого-физических и технологических параметров, таких, например, как проницаемость, пористость, глинистость и карбонатность пласта, вязкость пластового флюида, толщина продуктивных интервалов пласта, забойное и пластовое давление, начальная и текущая продуктивность, режимы виброволнового воздействия и др. Поэтому для исследования данного процесса и адекватного выбора прогнозирования его показателей необходимо использовать компьютерные модели, построенные с использованием методов многофакторного статистического анализа.

В гл. 11 описывается прогнозное моделирование, осуществленное в рамках системного подхода к выбору объектов, назначению технологических операций, реагентов, общей орга-

334

vk.com/club152685050 | vk.com/id446425943

низации внедренческих работ, реализованное при ОПР и внедрении технологий на Ново-Елховском месторождении АО "Татнефть".

На основе накопленной информации по исходным геологофизическим и технологическим параметрам проведенных вышеописанных опытно-промышленных работ осуществляли прогнозное моделирование процесса обработки ПЗП с использованием метода главных компонент [58].

Преимущества данного подхода особенно проявляются на начальных этапах пробного моделирования при необходимости учета в модели большого количества факторов и недостатке числа "наблюдений". Применение метода главных компонент позволяет заменять совокупность большого числа взаимосвязанных признаков на некоторую совокупность относительно малого числа некоррелированных параметров – обобщенных факторов процесса главных компонент с сохранением всей информации об изменчивости процесса.

Для моделирования процесса обработки нагнетательных скважин, вскрывающих девонские и бобриковские пласты Но- во-Елховского месторождения, было определено восемь переменных – факторов процесса, оказывающих влияние на показатель процесса – приращение коэффициента удельной приемистости. Это – удельная максимальная приемистость скважины по истории работы, удельная приемистость на момент перед обработкой, отношение значения пластового давления к гидростатическому, общая толщина продуктивных интервалов, коэффициент вариации толщин продуктивных интервалов, средневзвешенный коэффициент пористости, коэффициент депрессии режима виброволновой обработки, число комплексирования с реагентными операциями. Коэффициент удельной приемистости – коэффициент приемистости, отнесенный к 1 м толщины продуктивного интервала пласта, наиболее адекватно учитывает эффект осуществляемого воздействия непосредственно на среду коллектора призабойной зоны скважины. Коэффициент депрессии режима обработки, отношение забойного давления режима виброволнового воздействия к гидростатическому давлению глубины интервала перфорации определяют характер процесса очистки призабойной зоны (количественное выражение режима депрессии или репрессии на пласт). Коэффициент вариации толщин продуктивных интервалов скважины и другие подобные параметры ка-

335

vk.com/club152685050 | vk.com/id446425943

чественно оценивают неоднородность пластов и изменчивость геолого-промысловых параметров объекта.

В табл. 11.1 представлены исходные данные по проведенным на Ново-Елховском месторождении виброволновым обработкам нагнетательных скважин с различными геологофизическими и технологическими условиями.

На первом этапе осуществляется корреляционный анализ

Таблица 11.1

Исходные данные модели технологической эффективности процесса обработки нагнетательных скважин, вскрывающих пласты девонского и бобриковского горизонтов Ново-Елховского нефтяного месторождения

Номер

 

 

 

Параметры процесса

 

 

 

сква-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

жины

y

x1

x2

x3

x4

x5

x6

x7

x8

523

0,10

1,11

0,11

0,96

4,00

0,00

19,30

1,23

1,00

3189

0,00

1,15

0,26

1,10

3,00

0,33

14,50

1,00

2,00

3086

0,21

2,10

0,37

1,06

5,00

0,32

17,70

1,05

2,00

3093

0,10

0,15

0,00

0,97

8,00

0,69

17,00

0,98

1,00

2518

0,71

0,47

0,14

1,13

3,00

0,00

18,00

1,80

2,00

745

1,22

0,41

0,00

0,98

3,00

0,00

23,00

0,40

4,00

2122

0,22

0,22

0,00

1,25

4,00

0,00

17,30

1,10

1,00

4428

0,03

0,15

0,03

1,09

3,20

0,25

19,00

0,40

3,00

4437

0,11

1,79

0,04

1,12

2,40

0,50

19,90

0,40

3,00

134

0,08

0,63

0,02

1,07

14,0

0,33

15,00

0,20

3,00

689

0,41

0,86

0,10

1,09

5,60

0,40

19,80

0,25

2,00

4468

0,60

0,43

0,15

1,06

6,40

0,23

18,00

1,10

3,00

8432

0,19

0,67

0,23

0,91

3,90

0,23

17,50

1,10

3,00

7214

0,09

0,79

0,19

0,99

2,80

0,00

22,00

0,24

4,00

7267

0,07

0,43

0,14

0,80

3,60

0,00

19,50

0,20

4,00

7180

0,16

0,68

0,15

0,81

3,00

0,00

21,00

1,20

1,00

7184

0,47

2,50

0,14

0,81

1,40

0,00

22,00

1,23

2,00

7203

0,07

0,43

0,09

0,81

6,00

0,00

23,00

1,24

1,00

7277

0,24

0,72

0,10

0,86

2,00

0,00

21,00

1,16

2,00

7279

0,11

2,10

0,21

0,96

4,20

0,00

20,20

1,03

1,00

7281

0,39

1,40

0,10

0,87

2,00

0,00

20,50

1,15

2,00

7301

0,42

1,20

0,16

0,89

2,00

0,00

19,90

1,12

2,00

7375

0,10

0,80

0,10

0,72

3,00

0,00

22,20

1,40

2,00

П р и м е ч а н и е. y – приращение коэффициента удельной приемистости

после обработки, 10–1м2/(сут МПа); x1 – удельная максимальная приемистость

скважины по истории работы, 10–1м2/(сут МПа); x2 – удельная приемистостьна

момент перед вибровоздействием, м3/(сутатмм); x3 – коэффициент пластового

давления (отношение пластового давления к гидростатическому); x4 – общая

толщина продуктивных интервалов, м; x5 – коэффициент вариации толщин

продуктивных интервалов; x6 – средневзвешенный коэффициент пористости,

%; x7 – коэффициент депрессии режима виброобработки; x8 – число комплек-

сирования с физико-химическими мероприятиями.

 

 

 

 

 

 

 

 

 

 

 

 

 

336

 

 

 

 

 

 

 

 

 

vk.com/club152685050 | vk.com/id446425943

исходного материала. Результаты расчетов выборочных

337

vk.com/club152685050 | vk.com/id446425943

Таблица 11.2

Корреляционная матрица

Показатель

 

 

 

 

 

 

 

 

Коэффициенты корреляции

 

 

 

 

 

 

 

Средние

Среднеквад-

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значения

ратичные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

факторы

 

y

 

x1

 

x2

х3

x4

 

x5

 

x6

 

x7

x8

 

отклонения

y

 

1,000

 

–0,055

 

–0,170

0,096

–0,190

 

–0,238

0,281

 

0,098

0,272

0,27

0,28

x1

 

–0,055

 

1,000

 

0,512

–0,127

–0,294

 

–0,043

0,138

 

0,137

–0,164

0,92

0,66

x2

 

–0,170

 

0,512

 

1,000

–0,111

–0,232

 

–0,093

 

–0,152

 

0,262

–0,020

0,12

0,09

x3

 

0,096

 

–0,127

 

–0,111

1,000

0,245

0,434

 

–0,604

 

–0,227

0,061

0,97

0,13

x4

 

–0,190

 

–0,294

 

–0,232

0,245

1,000

0,478

 

–0,505

 

–0,295

–0,021

4,15

2,67

x5

 

–0,238

 

–0,043

 

–0,093

0,434

0,478

1,000

 

–0,565

 

–0,337

0,028

0,14

0,20

x6

 

0,281

 

0,138

 

–0,152

–0,604

–0,505

 

–0,565

1,000

 

0,032

0,066

19,45

2,33

x7

 

0,098

 

0,137

 

0,262

–0,227

–0,295

 

–0,337

0,032

 

1,000

–0,656

0,91

0,45

x8

 

0,272

 

–0,164

 

–0,020

0,061

–0,021

0,028

0,066

 

–0,656

1,000

2,22

1,00

 

 

 

 

 

 

 

 

 

 

 

Таблица 11.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Показатель

 

 

 

 

 

 

 

Критерии надежности коэффициентов корреляции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и факторы

y

 

x1

 

 

x2

 

x3

 

x4

 

x5

 

 

x6

 

x7

x8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1,00

 

0,27

 

 

0,84

 

0,46

 

0,95

 

1,21

 

1,46

 

0,47

1,41

x1

1,00

 

1,00

 

 

3,33

 

0,62

 

1,54

 

0,21

 

0,68

 

0,67

0,81

x2

1,00

 

1,00

 

 

1,00

 

0,54

 

1,17

 

0,45

 

0,74

 

1,35

0,10

x3

1,00

 

1,00

 

 

1,00

 

1,00

 

1,25

 

2,57

 

4,57

 

1,15

0,30

x4

1,00

 

1,00

 

 

1,00

 

1,00

 

1,00

 

2,97

 

3,25

 

1,55

0,10

x5

1,00

 

1,00

 

 

1,00

 

1,00

 

1,00

 

1,00

 

3,98

 

1,82

0,14

x6

1,00

 

1,00

 

 

1,00

 

1,00

 

1,00

 

1,00

 

1,00

 

0,15

0,32

x7

1,00

 

1,00

 

 

1,00

 

1,00

 

1,00

 

1,00

 

1,00

 

1,00

5,53

x8

1,00

 

1,00

 

 

1,00

 

1,00

 

1,00

 

1,00

 

1,00

 

1,00

1,00

338

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk.com/club152685050 | vk.com/id446425943

коэффициентов корреляции и их критерии надежности, оценивающие достоверность полученных коэффициентов корреляции, приведены в виде корреляционной матрицы и матрицы критериев надежности (табл. 11.2 и 11.3).

Анализ полученных результатов показывает, что существует зависимость между выбранными показателем эффективности виброволнового воздействия и геолого-технологическими факторами. При этом также наблюдается существенная корреляция между факторами процесса, т.е. данные обладают свойствами мультиколлинеарности. Это затрудняет проведение анализа и определение коэффициентов корреляционных моделей. Трудно установить, какие факторы влияют наиболее существенно и определить степень этого влияния. В результате регрессионного анализа получена линейная зависимость, прогнозирующая показатель при заданных значениях факторов. Критериями надежности данной модели являются мера идентичности и коэффициент множественной корреляции, характеризующие часть дисперсии показателя, которую описывает данная регрессия, F – отношение, характеризующее меру надежности данного выражения (в сравнении с процентилями распределения Фишера), оценка несмещенности и доверительные интервалы (вычисленные с помощью распределения Стьюдента), характеризующие ошибки предсказания. Числовые значения этих критериев приведены в табл. 11.4.

Наличие большого числа факторов (восемь) при небольшом числе наблюдений (23 скважины) делает необходимым для увеличения надежности прогнозной модели отбрасывание некоторого их числа, однако затруднительно определить влияние этого отсечения методом регрессионного анализа. Это вызвало необходимость следующего этапа моделирования – применения метода главных компонент.

На этом этапе появляется задача замены исходных взаимосвязанных признаков на некоторую совокупность не-

Таблица 11.4

Результаты регрессионного анализа по исходным факторам процесса

Искомое уравнение регрессии для показателя y: y = –0,22+ –0,1612х6+ +0,03х7

Мера идентичности: 0.09 Коэффициент множественной корреляции: 0,30 F-критерий: 0,97

Критерий несмещенности: 0,269 Доверительный интервал (90 %): 0,394 Попадание: 91,304

339

vk.com/club152685050 | vk.com/id446425943

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

– главных компонент. Эти новые переменные отражают внутренние, объективно существующие закономерности, не поддающиеся непосредственному наблюдению.

Вкачестве исходной информации для метода используется полученная ранее корреляционная матрица (см. табл. 11.2), применяемая как ступень для дальнейшего анализа исходных признаков. Задача заключается в линейном преобразовании m

признаков X1, X2, ..., Xm в новый набор случайных величин K1, K2, ..., Km, что делает их некоррелированными и располагает в порядке убывания их дисперсий. При этом подавляющую долю дисперсии, которая характеризует изменчивость процесса, обычно содержат первые несколько компонент. Выделенные компоненты рассчитываются для исходных объектов, при этом по группируемости объектов – скважин в координатах главных компонент извлекается ценная дополнительная информация, позволяющая классифицировать объекты по заданному признаку. Регрессионная модель, в которой в качестве переменных факторов используются главные компоненты, позволяет устранять отмеченные выше трудности.

Врезультате проведенного компьютерного анализа было извлечено восемь компонент процесса. Полученная матрица весов переменных – исходных факторов, входящих в главные компоненты, собственные значения главных компонент и вклады компонент в общую дисперсию приведены в табл. 11.5. Первая компонента объясняет 34 % всей дисперсии, вторая – 22 %. Основную долю дисперсии процесса (72 %) описывают три первые главные компоненты. Отметим, что относительно малое число первых главных компонент, накапливающих основную дисперсию процесса, еще раз свидетельствует о сильной коррелированности исходных факторов.

340

vk.com/club152685050 | vk.com/id446425943

Таблица 11.5

Результаты анализа процесса обработок скважин методом главных компонент

Показа-

Главные компоненты процесса обработки ВДХВ нагнетательных скважин

тель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

процесса

к1

 

к2

 

к3

 

к4

к5

 

к6

 

к7

 

к8

x1

–0,419

 

0,427

 

0,584

0,292

 

0,379

0,207

0,173

 

 

–0,007

x2

–0,333

 

0,539

 

0,610

 

–0,034

–0,414

 

–0,035

 

–0,222

 

0,055

x3

0,679

 

0,268

 

0,098

 

–0,540

0,308

0,210

 

–0,091

 

0,135

x4

0,725

 

0,098

 

–0,223

0,468

 

–0,275

0,318

0,006

 

 

0,139

x5

0,750

 

0,271

 

0,146

0,283

 

0,226

 

–0,449

 

–0,021

 

0,097

x6

–0,706

 

–0,566

 

–0,048

0,174

 

0,241

0,012

 

–0,181

 

0,240

x7

–0,558

 

0,590

 

–0,447

 

–0,169

–0,131

 

–0,109

0,202

 

 

0,203

x8

0,258

 

–0,660

 

0,584

 

–0,174

–0,240

 

–0,054

0,222

 

 

0,131

 

 

 

Собственные числа главных компонент

 

 

2,712

 

1,732

 

1,337

 

0,766

 

0,669

 

0,406

 

0,211

 

 

0,166

 

 

 

 

 

 

 

 

 

 

 

 

Вклады компонент в общую дисперсию

 

 

0,339

 

0,217

 

0,167

 

0,096

 

0,084

 

0,051

 

0,026

 

 

0,021

 

 

 

 

 

 

 

 

 

 

Накопленные доли дисперсий по главным компонентам

 

 

0,339

 

0,555

 

0,723

 

0,818

 

0,902

 

0,953

 

0,979

 

 

1,000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Анализ матрицы весов главных компонент позволил провести классификацию обобщенных факторов процесса. Так, в первую главную компоненту с наибольшими весами входят коэффициент пластового давления, общая толщина продуктивных интервалов, коэффициент вариации толщин продуктивных интервалов и средневзвешенный коэффициент пористости. Ее можно классифицировать как обобщенную гидропроводность скважины. Вторая главная компонента с наибольшими весами включает коэффициент депрессии режима виброволновой обработки, число комплексирования с реагентными операциями. Ее следует классифицировать как обобщенный фактор виброволнового режима обработки. Третья главная компонента процесса с наибольшими весами включает число комплексирования с физико-химическими мероприятиями и факторы, определяющие "восприимчивость" скважины к проведению физико-химического воздействия, – максимальную приемистость и удельную приемистость скважины на момент перед воздействием. Ее можно классифицировать как обобщенную характеристику физико-химического компонента воздействия.

Каждая главная компонента выражается через исходные факторы процесса в виде

341

vk.com/club152685050 | vk.com/id446425943

kr = 1 (a1rx1+ a2rx2 + ....+ amr xm) Vr

где air – элементы матрицы весов главных компонент (см. табл. 11.5); xi – переменные факторы процесса; Vr – собственные значения главных компонент.

В табл. 11.6 представлены главные компоненты, рассчитанные для набора исходных факторов промысловых объектов.

На основе полученных результатов извлекается ценная информация по "распознаванию" промысловых объектов по признаку эффективности для проведения виброволновых обработок скважин. На рис. 11.1 представлено распределение обработанных нагнетательных скважин девонских и бобриковских пластов в координатах первых главных компонент. Видно, что все скважины весьма отчетливо разделяются на три группы: группу скважин с высоким эффектом (приращение коэффици-

Таблица 11.6

Главные компоненты, рассчитанные по исходным факторам процесса обработки нагнетательных скважин

Номер

 

Главные компоненты по исходным факторам процесса

 

сква-

 

 

 

 

 

 

 

 

 

 

жины

к1

к2

к3

к4

к5

к6

к7

к8

523

–4,06

–5,58

–0,73

6,08

5,71

4,31

–13,80

34,33

3189

–2,83

–4,40

0,25

4,23

4,18

3,02

–9,11

27,41

3086

–3,31

–5,06

0,25

6,55

4,95

5,14

–11,06

33,73

3093

–1,98

–4,85

–1,63

7,94

3,05

6,20

–12,69

34,45

2518

–3,87

–5,52

–0,52

4,45

4,88

2,96

–11,75

33,27

745

–4,71

–8,48

0,54

5,54

6,22

3,12

–15,11

40,21

2122

–3,29

–5,18

–1,03

5,11

4,71

3,99

–12,95

31,51

4428

–3,57

–6,78

0,15

4,89

5,02

2,94

–13,04

34,05

4437

–4,20

–6,67

1,00

5,30

6,70

2,91

–12,54

34,80

134

0,34

–4,81

–1,24

10,85

–0,53

11,48

–9,08

37,04

689

–3,28

–6,36

–0,31

7,12

5,12

5,21

–14,35

36,32

4468

–2,67

–5,93

–0,42

6,57

3,27

5,37

–11,31

36,12

8432

–3,30

–5,85

0,15

5,13

4,14

3,43

–10,78

33,21

7214

–4,55

–8,06

0,91

5,35

6,07

3,16

–14,31

38,46

7267

–3,66

–7,35

0,69

5,28

4,59

3,44

–12,34

35,30

7180

–4,74

–6,32

–0,80

5,80

6,40

3,28

–15,65

35,83

7184

–5,62

–6,67

0,65

5,51

8,09

2,84

–13,97

36,69

7203

–4,42

–6,87

–1,52

7,99

5,77

5,56

–17,37

41,26

7277

–4,89

–6,77

–0,18

4,96

6,53

2,43

–14,60

35,76

7279

–4,36

–5,66

–0,25

6,82

6,49

5,04

–14,05

35,55

7281

–4,86

–6,44

0,14

5,10

6,75

2,77

–13,63

35,00

7301

–4,67

–6,28

0,11

4,87

6,39

2,66

–13,38

34,14

7357

–5,03

–7,03

–0,45

5,92

6,49

3,15

–15,24

38,50

342

vk.com/club152685050 | vk.com/id446425943

ента удельной приемистости ∆Kпр > 0,3), группу скважин со средним эффектом (приращение ∆Kпр от 0,09 до 0,24) и группу скважин с низким эффектом воздействия (∆Kпр < 0,09).

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

Следующим этапом моделирования является регрессионный анализ на основе определенных главных компонент. Полученные результаты представлены в табл. 11.7. Сравнение данных таблицы с соответствующими данными, установленными в результате регрессионного анализа по исходным факторам, показывает увеличение надежности по F-критерию и уменьшение доверительных интервалов, отвечающим одина-

Рис. 11.1. Распределение скважин в координатах главных компонент модели эффективности виброволновых обработок нагнетательных скважин с ис-

пользованием технологий ВДХВ и ВПВ:

1 – группа скважин сильного эффекта обработки (∆Kпр > 0,3); 2 – группа скважин среднего эффекта обработки (∆Kпр от 0,09 до 0,24); 3 – группа скважин слабого

эффекта обработки (∆Kпр < 0,09)

343