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

книги из ГПНТБ / Сохранов Н.Н. Машинные методы обработки и интерпретации результатов геофизических исследований скважин

.pdf
Скачиваний:
7
Добавлен:
25.10.2023
Размер:
14.73 Mб
Скачать

ВхоЗ

\

 

Подготовка программы

 

к работе

 

И: =\Ч- I

Нет

" „ > 5 „ Л « , ==: 1,05rf„Ap > 1 5 р С 1

 

 

Да

 

Подсчет числа п.

R:=s: = 0; р < « » : =р<°>; р'*': = „ „ „ „ .

мальное число

Вычисление ^15"'

Нет

 

Да

 

Да

 

 

РГ^РГ

 

РГ--=РГ%

с

]

с

 

 

J : =5+ I

 

 

Да

«ей

1—J Я йеО VS> 1

0 3 '

 

 

 

Выход

Вычисление рт6,

 

6* sse Нет АВОСТ

 

 

Да f-

 

 

 

 

Wem

 

Занесение рс(0)

 

 

 

В таблицу р^1

m —0

Нет

АВОСТ

 

 

" < " о п

Да

 

лг

 

 

АВОСТ

Вычисление m

 

 

 

Присвоение меток L.

 

 

тф 1

Да

 

АВОСТ

 

 

Настройка блока счета п.\

: отключением помеченного \

L,,

отключение блока

J

2 Б „

Вычисление К:

р^: = помеченное p-f

Lj Ф ІИКА^Ф max i r j

нет

Kj„:-0

 

 

Поиск Lj — (LA или Lg)

(л, =",- IAnfi = \)V(пА = 1Лпв=п• -1)

Вычисление Б Ѵ ,

Да se.- s^?

увеличение ппп на f

Рис. 41. Блок-схема алгоритма контроля качества.

X знак присваивания; величина, записанная слева от знака присваивания, полагается равной

значению арифметического выражения, записанного справа от знака.

р\Ѵ

Б л о к

П.

Производят

проверку

знака

величины

p\t'.

Если

^

0,

вычисляют

балл

Б

— оценка

качества

измерений

для

данного пласта. Для этого определяют величину p^o.s

в

предположе­

нии, что допустимая

погрешность

измерений составляет

0,5

о. Если

полученное

при

этом

значение

pj^o.s ^

0, то считается, что качество

измерений

хорошее,

Б д

=

2;

при

plsl,b

<L 0

качество

измерений

удовлетворительное,

Б^

=

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные

значения

Бц и р с

( 1

поступают

на

счетчики,

сумми­

рующие их; показание счетчика опорных пластов поп

 

увеличивается

на единицу. Далее управление передается в блок I (вызов информа­

ции о следующем пласте с номером

ц. +

1 и т. п.).

 

 

 

 

 

 

 

 

Если исходная величина р(^<С0,

то выявляют зонды с плохим

качеством измерений. Для этого подсчитывают число

ПА

ЗОНДОВ,

интервалы

[А/,

В А

которых

не

пересекаются

(max

А

—В{

£>0)

с интервалом [max А, В Л зонда

ЬА

с максимальным

 

значением

А,

и

число

пв

зондов,

интервалы

[AjB/]

 

которых

 

не

 

пересекаются

(min В — Aj<ü0)c

интервалом ,

min В] зонда L

B с минимальным

значением В. Если состояния счетчиков чисел пАж

пв

удовлетворяют

условиям

1 <

пл

•< п3 1 Е 1 <

 

пв <

и3 — 1 или

условию

ПА

+

+

пв

ф

п3,

то данный пласт исключается из числа пластов,

по кото­

рым можно уверенно выявить зонд с плохим качеством измерений.

При ПА =

1 И ПВ = — 1 делается вывод, что показания

зонда с дли­

ной L B не согласуются с показаниями других зондов. Если пА = п3

— 1 и пв

1, то делается вывод, что показания зонда

с длиной L

A

не согласуются с показаниями других зондов. При этом показания счетчиков непересечений (несогласий) для соответствующих зондов п} увеличиваются на единицу (если это не самый большой зонд для рассматриваемого пласта); кроме того, увеличивается на единицу показание поп счетчика числа опорных пластов, против которых вы­ явлено по одному некачественному замеру КС.

Для исправления показаний выявленных таким образом зондов с ошибочными значениями кажущихся удельных сопротивлений

вычисляют

отношение

 

 

 

 

 

 

Я/и =

Р/Ѵ/Р/іч»

(40)

где pjllp

— кажущееся

удельное

сопротивление, рассчитанное для

данного

зонда по р п / р с ,

установленному по показаниям остальных

зондов

комплекса.

 

 

 

Вычисленное значение K J ß вместе с величинами Я ц

и L S запоми­

нается в специальной таблице СТ.

 

Б л о к

I I I . После того как перебраны все нласты,

производится

анализ полученных результатов. В случае, когда число опорных пла­

стов

о п меньше

некоторого наперед заданного

предельного

числа

поп

g опорных

пластов, при которых можно проводить уточнение

удельного

сопротивления бурового раствора, сравниваются

пока­

зания 7гоп и п соответствующих

счетчиков. Если

поп

>> п,

возмож­

ности для

уточнения удельного

сопротивления pt

нет.

При

поц

<j п

111

путем анализа таблицы CT выявляют зонды с числом ïij непересече­ ний, превосходящим значение

 

С2 (Сj ~f

мо п ),

 

 

 

 

где Сj — показания счетчика

числа

замеров

КС

/-ым

зондом

(по

таблице CT);

С 2 — постоянная

выбираемая в

пределах

0,5—1.

 

Зонды с

iij s= С 2 (Сj ~г поп)

отмечают (помечают L-), число

m

таких зондов считают счетчиком т. При m >> 1 анализ

прекращают.

При m = 1

повторяют процесс уточнения величины удельного со­

противления

бурового раствора (блок I) без учета

показаний зонда

с числом непересеченин, превосходящих допустимое. При повторном контроле значения удельного сопротивления бурового раствора опе­ ратор суммирования баллов качества отключается (сумма Б ѵ не изменяется). Если после повторного определения опять наблю­

дается

неравенство

поп

>< ?гоп

3 ,

то

анализ

прекращается;

если

о п

будет

больше и о п - 3 ,

то операции

проводятся

применительно к

дан­

ному

случаю.

 

 

 

 

 

 

 

 

 

 

При

о п

;> поп

g

вычисляется

дисперсия

а 2

арифметического

среднего

из

значений

р ^ ,

определенных

для

опорных

пластов.

Если

число

опорных

пластов не

очень велико,

например

меньше

наперед заданного числа поп 3 1 , то в качестве дополнительного ис­ пользуется также значение ре01. Из имеющихся значений р"^ удель­ ного сопротивления бурового раствора отбраковываются величины, наиболее отличающиеся от среднего, например по критерию За.

Затем по оставшимся значениям определяется отношение 2 Ро8 ) /^0 п и вычисляется для него дисперсия а\. Если полученная таким образом величина дисперсии меньше допустимой погрешности е определения

удельного сопротивления бурового раствора, то 2 Pc цМоп считают истинным удельным сопротивлением бурового раствора. При стх •> е возникает аварийная ситуация, выход из которой представляется сделать интерпретатору (принять рс за истинное р с п и т. п.).

Общая оценка качества измерений дается средним баллом Бср, вычисляемым как сиеднее из баллов по отдельным опорным пластам:

Если показания m счетчика помеченных длин зондов L; - не равны нулю, на основе анализа CT вычисляют корректирующие коэффи­ циенты Kj. Их величина определяет необходимость (К,- =5 Amin) и возможность (Kj ^ Ктп.л) исправлений показаний /-го зонда.

Глава VII

РАСЧЛЕНЕНИЕ РАЗРЕЗА СКВАЖИН

17. Р А С Ч Л Е Н Е Н И Е РАЗРЕЗА СКВАЖИН ПО ГЕОФИЗИЧЕСКИМ Д А Н Н Ы М

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

Рис. 42. Форма кривых против различных пластов.

а — пачка

прослоев

высокого

и низкого

сопротивления;

б •—

однородный

пласт

высокого

сопротивления;

в — пласт с

плавным

изменением

сопроти­

вления

по глубине (переходная

зона);

г

— пласт

с

водо-неф-

тяным

контактом;

J — пласт

высокого сопротивления,

выде­

ленный на

ЭВМ

по

экстрему­

мам кривой КС малого гра- дпсит-зопда; 2— пласт высокого сопротивления, выделспный по кривым КС малого и большого

зондов путем объединения

близ­

ких по

сопротивлению

fnpo-

слоев; 3 — переходная

зона;

4 — зона

проникновения

(сту­

пень); 5 — пласт низкого

сопро­

 

тивления.

 

Ё

Р

по разрезу; рис. 42). Для однородных одиночных пластов этот ана­ лиз сводится к нахождению на кривой характерных точек (макси­ мумы, минимумы, точки перегиба и др.), положение которых соот­ ветствует границам пластов (прослоев).

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

8 Заказ 111

113

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

Отбивка границ пластов по всем характерным точкам привела бы к излишне дробному расчленению разреза, усложняющему интер­ претацию данных каротажа. Например, при интерпретации данных

Б К З выделение всех прослоев неоднородного пласта

(рис. 42, а)

значительно увеличило бы объем работ и не дало бы

возможность

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

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

Таким образом, расчленение разреза скважин (выделение пластов и отбивка их границ) по промыслово-геофизическим данным можно свести к следующим этапам:

1)нахождение характерных точек (экстремумов, точек перегиба) кривых, соответствующих однородным прослоям и пластам;

2)выделение среди них существенных путем объединения в пласт отдельных прослоев с близкими значениями исследуемого физи­ ческого свойства;

 

3) уточнение границ пластов с

плавным изменением

свойств

к

границам.

 

 

 

Рассмотрим различные способы

расчленения разреза

скваяшн

по

геофизическим данным.

 

 

Аналитический способ

В аналитическом способе [26, 27] отбивки границ характерные точки выделяются путем анализа разности геофизических величин àdj; в соседних точках (градиента кривой при постоянном шаге квантования):

Aß/i = ац -

(42)

114

где / — порядковый номер (шифр) интерпретируемого геофизи­ ческого параметра.

Положение характерных точек определяют по смене знака этой разности (точки экстремумов кривых кажущегося сопротивления, снятых градиент-зондами) и по максимуму абсолютного его значения (точки перегиба кривых ПС, гамма-каротажа, кажущегося сопро­

тивления, снятых зондами индукционного и бокового

каротажа).

За характерные точки могут быть взяты также точки со

значениями

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

аг = ^(атіп-гК0ата:{),

(43)

где Кг и Кг — коэффициенты, зависящие от вида каротажа и гео­ метрических параметров пласта и вмещающих пород.

Критериями выделения существенных характерных точек — объединение прослоев, мало отличающихся по величине исследу­ емого параметра, — в аналитическом способе могут служить следу­ ющие выражения [26, 27, 37]:

 

 

а тах

^тіп^*^!

 

 

(44)

 

 

|zm J x z m i n l > T ) ;

 

(45)

 

 

 

I Датах I > С;

 

 

(46)

 

 

"max -

° т In >

 

ff.

 

( 4 7 )

 

 

 

flmax max

- > # i ,

 

(48)

 

 

 

 

 

 

 

где а т а х и а т І п

— соседние

максимум и

минимум

кривой;

zmax и

z m m глубины

максимума

и

минимума; Д а т а х — максимальная

разность

значения в интервале

между максимумом

и минимумом;

и, и, С, К

и К1

— некоторые постоянные,

устанавливаемые

на осно­

вании опыта обработки данных каротажа по району. В некоторых случаях вместо максимумов и минимумов используют максимум-

макеиморум a m a x m a x

и минимум-миниморум

а т і а т і п .

использован

Аналитический

способ отбивки границ

пластов

в ряде программ [27 , 37 , 41], а алгоритм нахождения

характерных

точек в последовательности чисел а/ 7 применяется во всех програм­ мах отбивки границ пластов. Программы отбивки границ пластов аналитическим способом опробованы на ряде месторождений с пес- чано-глинистым разрезом. Выделение пластов этим способом по кри­ вой ПС обычно дает хорошие результаты (рис. 43, II). Однако при отбивке границ по кривой КС не всегда получают удовлетворитель­ ные результаты. Это связано с тем, что в сложном разрезе трудно

8*

115

подобрать значения х, л, С и К в критериях определения суще­ ственных характерных точек. Указанное обстоятельство является. иедостатком аналитического способа.

Рис. 43. Отбивка границ пластов по кривой ПС способом моделирования границ

(/) и аналитическим способом с выделением ЦІ) и без выделения (111) сущест­ венных характерных точек.

Математическое моделирование

В способе математического моделирования [46] границ пластов характерные точки выделяют взаимной корреляцией последователь­ ности чисел а/ 7 , представляющих интерпретируемую кривую, с по­ следовательностью чисел ап!]-, изображающих модель кривой на границе пласта. В качестве модели ат]- обычно используют последо­ вательность чисел

 

amj

= ï при

7 =

1,2, ... , т / 2 ;

 

 

 

 

 

(49)

amj

= — 1

при

7 =

— - j - 1 , — + 2, . . ., m,

где m — длина

(количество

чисел)

модели.

Экстремумы функции взаимной корреляции, определяемой для этой модели по формуле (36), располагаются между точкой перегиба кривой и следующим за ней экстремумом. Поэтому выделение пла­ стов этим способом сводится к следующему:

1)вычисляют функцию взаимной корреляции для последователь­ ности чисел а,-,- интерпретируемой кривой и последовательности чисел a„,j модели границы пластов, перемещая каждый раз модель на один шаг квантования по глубине;

2)вблизи экстремума функции взаимной корреляции на интер-

116

претируемой кривой (последовательности а; 7 ) находят характерные точки. Поиск характерных точек осуществляется аналитическим способом.

Математическое моделирование границ обладает следующими существенными преимуществами перед аналитическим способом:

1) помехи измерения, иззубренности кривой, погрешности пре­ образования и флуктуации кривой практически не оказывают влия­

ния на результаты

выделения

характерных точек;

 

 

2) детальность расчленения разреза на пласты можно регулиро­

вать изменением

 

длины

(числа

т)

модели

границ; однако

число т,

в модели

должно

быть

ие

 

меньше

двух,

 

 

20 Он-м

так как при m = 2 функция

взаимной

 

 

корреляции для

модели

(49) будет

равна

 

 

 

• разности

Да;(-

и

результаты

выделения

 

 

 

характерных

точек

способами

моделиро­

 

 

 

вания границ и аналитическим будут

 

 

 

идентичны.

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно из рис. 43, на котором

 

 

 

приведены

результаты

расчленения

раз­

 

 

 

реза по кривой ПС разными способами,

 

 

 

иззубренность кривой не влияет на ре­

 

 

 

зультаты отбивки

границ

способом

мате­

 

 

 

матического моделирования (см. рис. 43,1).

 

 

 

При отбивке же границ аналитическим

 

 

 

способом

без

выделения

 

существенных

 

 

 

характерных

точек иззубренность

кривой

 

 

 

приводит к излишне дробному расчлене­

 

 

 

нию разреза (см. рис. 43,

777).

Хорошее

 

 

 

расчленение

разреза

оказалось

возмож­

Рис. 44. Модели кривой для

ным получить при отбивке границ по

уточнения

границ

пластов;

взаимной корреляцией.

существенным

экстремумам,

выделяемым

1 — кривая

КС; г

— модели;

на кривой

по

подобранным

для

разреза

з — точки

моделей,

использу­

емые для корреляции.

критериям

(44—46).

 

 

 

 

 

 

 

 

 

 

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

целесообразно использовать

следующие

модели:

 

атак

1) переходной

зоны — прямая,

соединяющая

максимум

с минимумом

amin

(рис. 44,

/ ) , или

с

точкой

с

ординатой

патІѢ

 

 

•атіі = пат,п

+ (і-1)

" m a * ~ " a m i "

,

 

(50>

где 1 < п <

2;

 

 

 

 

 

 

 

И Т

2) пласта

низкого

сопротивления — вогнутая

гипербола,

опи­

рающаяся на точки максимума и минимума

(рис. 44, ТУ)

 

 

П •

П •

! " п м х ~ "min

 

 

 

I N

 

"m ит • • ч " ш ш т

 

{-(-1) ' * ' " *

 

\OLJ

3) пласта

высокого

сопротивления — выпуклая

гипербола,

опи­

рающаяся на точки максимума и минимума (рис. 44, /77)

 

_

(т ~\~2) Ащах

n min

i

"max

g min

, • • • .

l^r>\

«m!n. • - M

J^Çl

 

I

—y

(04)

где i — номер точки модели по ходу обработки

геофизических

дан­

ных. Длина моделей m всегда равна

длине

участка анализируемой

кривой.

 

 

 

 

 

 

 

 

 

Результаты корреляции интерпретируемой кривой о;-(- с различ­

ными моделями amj показаны в табл.

 

7. Как видно, из табл. 7, по

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

 

 

 

 

Т а б л и ц а 7

Значения

фупкцпп

взаимной корреляции кривой

КС (см. рис. 44)

 

 

с различными

моделями

 

Модель

Вогпутая

Прямая линия

Выпуклая

гипербола

гипербола

Значение

функции

0,96

0,88

0,40

1)коррелируемость кривой с вогнутой гиперболой (R (б) >>0,95) означает, что в исследуемом интервале находится пласт низкого сопротивления и максимум <зт а х является существенным;

2)коррелируемость кривой с выпуклой гиперболой указывает на пласт высокого сопротивления и необходимость перенесения гра­ ницы пласта из максимума в точку А, удаленную на наибольшее расстояние от прямой, проведенной через максимум и минимум иссле­

дуемой кривой;

3) коррелируемость кривой с прямой линией определяет нали­ чие переходной зоны.

Указанные выше преимущества математического моделирования при отбивке границ пластов позволили широко использовать его для расчленения разреза скважин по каротажным кривым. На его основе во ВНИИГеофизике составлены алгоритм и программа автомати­ ческой отбивки границ пластов по кривым КС (см. следующий раздел).

118

Полуавтоматический способ

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

Рис. і5. Результаты обработки данных БКЗ па ЭВМ.

V7' п I — 7 — соответственно кровли и подошвы выделенных пластов.

Примерные границы пластов отмечаются интерпретатором на диаграмме (см. рис. 12 и 45). При преобразовании диаграмм в цифро­ вую форму преобразователем Ф001 эти отметки наносятся на перфо­ ленту (перфокарту). Глубины границ могут задаваться также табли­ цей. Отмеченные числа находятся ЭВМ и вблизи них определяются характерные точки, которые принимаются за существенные (границы пластов). По этому способу составлены алгоритм и программа полу­ автоматической отбивки границ пластов при интерпретации Б К З (см. ниже).

119>

Соседние файлы в папке книги из ГПНТБ