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

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

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

вие переходной зоны. В зависимости от выполнения этого

условия

глубину zA

точки Атах

записывают в ячейку R

(z0) или в ячейку

 

 

 

 

R (z m a x

m a x ) хранения глубины мак-

 

 

Вход

 

симум-максиморума.

После

этого

Поиск

Q,max

max1 гmax max

повторяют

действия,

начиная

с

пункта

 

3,

описанного

алгоритма.

и аmin mtri'•

min

min

 

 

 

 

 

10.

Первое

условие

существо­

Поиск

точна

 

вания

переходной

зоны

(63)

не

 

выполнено — переходной зоны нет.

 

 

 

 

11.

Засылка

глубины

z m a x m a x B

Нет

 

 

 

ячейку R (zm a x )

хранения

первого

 

 

 

 

существенного максимума на кри­

 

 

 

 

вой

КС большого

зонда в интер­

Определение

Q m n ^

моде­

вале z m i n

V

2m invn, если z m a x m a x ^>

 

ли ПЗ • 3

 

l> zmax-

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

12.

 

Определение

 

глубины

 

Вычисление

 

zn 2 = z m a x

подошвы

пласта

вы-

 

 

 

 

сокого

 

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

и

запись

 

 

 

 

zn

и

величины

^ 4 т а х

в

табли­

 

 

 

 

цу ТГ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С о п о с т а в л е н и е

 

г р а ­

 

 

 

 

н и ц

п л а с т о в .

Влияние

зоны

 

 

 

 

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

на результаты от­

 

 

 

 

бивки

границ

пластов

по

кри­

 

 

 

 

вой КС, записанной малым зон­

 

 

 

 

дом,

можно

установить

сопоста­

 

 

 

 

влением ее с кривой КС, заре­

 

 

 

 

гистрированной

большим

зондом.

 

 

 

 

Кривые (и зонды), которые ис­

 

 

 

 

пользуются

для

отбивки

границ

 

 

 

 

пластов,

называют

основными

в

 

 

 

 

комплексе БКЗ . Алгоритм сопо­

 

 

 

 

ставления

основных

кривых КС

 

 

 

 

и отбивки

границ

пластов

по ним

 

 

 

 

показан на рис. 50. Сопоставление

 

 

 

 

производят в интервале минимум—

'-тахтах

 

 

минимум

[37]. Оно уточняет по­

 

 

ложение

границ

пластов,

опре­

 

 

 

 

Запоминание

 

деляет характер пласта и способ

 

отсчета

 

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

значений

 

 

 

 

 

 

 

 

 

кажущегося

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

про­

Рис. і 9 . Уточнение границы пласта (BHK,

тив

него.

В

результате

сопоста­

ГВК) в переходной зоне.

вления

 

границ,

 

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

 

 

 

 

малому

 

и

большому

основным

зондам, получают семь следующих вариантов отбивки границ пла­

стов между минимумами z m i n V

и z m i n V + i :

1) пачка из двух пластов

низкого и высокого сопротивлений;

границы, отбитые по малому и большому основным зондам, совпадают.

130

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

3) в подошве пласта высокого сопротивления отмечается пере­ ходная зона. В этом случае граница проводится по водо-нефтяному контакту, т. е. переходная зона включается в пласт высокого со­ противления;

I

Вычисление

R(8)

в интервале min-min

Отказ

Выаор (2.

вариант I

Отказ

Вариант 2

Рис. 50. Сопоставление границ пластов, выделенных по кри­ вым КС, заппсаппым малым п большим зондами.

4) пачка из двух пластов низкого и высокого сопротивлений по показаниям большого зонда в интервале пачки нет границы пласта — возможно повышающее проникновение раствора в верхней части пласта;

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

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

9* 131

7) пачка состоит из трех пластов; границы, определенные по основ­ ным зондам, не совпадают.

Для выделения пластов по кривой КС большого зонда исполь­ зуют тот же алгоритм (см. рис. 48), что и для малого зонда. Отличие состоит лишь в том, что интервал обработки ( z m j n ѵ zm\n v+]) задается в результате интерпретации кривой КС малого зонда. В про­ цессе реализации этого алгоритма так же, как и при обработке кри­ вой КС, записанноймалым зондом, проверяются условия (66) и (67). Если максимумы кривой КС в заданном интервале не удовлет­ воряют этим условиям, то для отбивки границ и определения мощ­ ности пластов используют только данные малого зонда. Варианты

отбивки границ устанавливают по числу максимумов

77imax

н а кривой

малого зонда. В случае,

когда »іпШ Х = 0, отбивают

пласт

«ступень»

(вариант 5). Е с л и / » т а х

== 0, то отбивают два пласта, в одном из кото­

рых предполагают повышающее проникновение раствора (вариант 4). При выполнении условий (66) и (67) проверяют наличие суще­

ственных максимумов на кривой малого зонда. При те^ах =

0 пласт

не интерпретируют. Если mäa x

0, сравнивают глубины

выделен­

ных границ. Предварительно из нескольких существенных макси­ мумов в рассматриваемом интервале па кривой КС большого зонда (прп тЯХ >> 1) выбирают один, первый (нижний, если градиентзонд подошвенный), удовлетворяющий условию (67), и его глубина засылается в ячейку z^a x .

Далее проверяют совпадение глубин границ пласта высокого сопротивления г™пах и z^àx, определенных соответственно по кри­ вым КС малого и большого основных зондов. При небольшом рас­ хождении границ I Sm'iJx Zmäx | ^ 0,8 в качестве границы пласта высокого сопротивления принимают показания малого основного зонда, пачка разбивается на два пласта — пласт низкого и пласт высокого сопротивлений (вариант 1). Для отсчета существенных зна­ чений КС против пластов используют блоки I и I I (см. с. 141—144) соответственно для пласта низкого и высокого сопротивлений. В слу­ чае расхождения между границами более чем на 0,8 м анализируется число существенных максимумов на кривой КС малого основного зонда.

При невыполнении условий z^ax — г ^ х

< 0

и m^lx ;> 1

макси­

мум на кривой КС большого основного

зонда

расположен

выше,

чем на кривой КС малого зонда, фиксируют расхождение в отбивке

подошвы пласта высокого сопротивления (вариант

2). В случае

пітіх >> 1 пачка разделяется

на три пласта путем

сопоставления

максимумов кривых между

собой.

 

При совпадении максимумов на кривых малого и большого основ­

ных зондов, т. е.

I Zmàx — zmàx I <С 1,2 м, пачка разбивается на три

пласта — низкого

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

мощностью Нъ промежуточного

сопротивления — Н2 и высокого

сопротивления — Н3, границы ко­

торых совпадают по обоим основным зондам (вариант 6). Снятие су­

щественных значений

КС против выделенных пластов производят

по соответствующим

алгоритмам (см. рис. 56—58).

132

ш

г

н

3

X

 

 

 

 

Рис. 51. Результаты автоматической интерпретации данных БКЗ в скв. 50 Туймазы (выделение пластов, отбивка их границ, отсчет существенных значений, определение удельного сопро­

тивления).

1 — пласт высокого сопротивления; 2 — пласт еще более высокого сопротивления; з — пласт низкого сопротивления; 4 —неинтерпретируемый на ЭВМ пласт; I—результаты

ручной интерпретации; XI — то же, на ЭВМ.

При несовпадении максимума, отбитого на кривой КС большого основного зонда, ни с одним из максимумов КС малого зонда в рас­ сматриваемом интервале (вариант 7) предусмотрена дальнейшая интерпретация по варианту 6 или 2 по выбору интерпретатора.

Если разность глубин между максимумом на кривой КС большого основного зонда и подошвой рассматриваемой пачки

zmax^ т и і ѵ ^ З м ,

(69)

то возможно существование переходной зоны. Необходимо провести анализ кривой КС большого зонда между Zmàx ^ -^тѴп V' Его произ­ водят по алгоритму уточнения границ в переходной зоне (см. рис. 49).

При ^тах •> Zmax сопоставление границ не представляется воз­ можным и разбивка пачки на три пласта проводится только по дан­ ным малого основного зонда.

Программы автоматического выделения пластов на ЭВМ для определения их удельного сопротивления опробованы в скважинах различных районов. Отбивка границ производилась по данным КС двух зондов: А0,8М0,1 и А4М0,5. Для обоих зондов граничное зна­ чение коэффициента кажущейся анизотропии л.к.г было принято рав­ ным 1,05, а модель границ пласта состояла из десяти точек (т — 10).

Обработка и интерпретация данных зондов Б К З на ЭВМ по ука­ занной программе дали удовлетворительные результаты [24]. Ошибка определения глубин границ пластов в большинстве случаев (90—95%) не превышала одного шага квантования А (0,2 мм). Существенные расхождения в отбивке границ на ЭВМ и ручным способом наблю­ даются в разрезе с неоднородными пластами, где определение гра­ ниц пластов представляет трудную задачу, и результаты, получен­ ные различными интерпретаторами, обычно не совпадают. Пример автоматической отбивки границ пластов на ЭВМ по кривым КС,

зарегистрированным

 

зондами

Б К З , дан

на рис. 51.

 

Полуавтоматическая

отбивка

границ пластов при интерпретации

 

 

 

данных БКЗ

 

 

Рассмотрим алгоритм полуавтоматической отбивки границ

пластов [38]

(рис.

52).

 

 

 

1. Вызов

чисел

аь

представляющих

основную кривую

КС, и

присвоение знака «минус» числам с признаком характерной

точки.

2. В зависимости от знака числа я(- и наличия признака границы

в ячейке R (т) выполнение следующих

операций:

 

а) знак числа а,

отрицательный — засылка знака разности Ад,-

и числа at соответственно в ячейки R ( ± )

и R (иг). При знаке «плюс»

в ячейке признака кровли, и подошвы R (к/п) счетчик пластов ѵ

увеличивается на единицу;

 

 

 

б) знак числа aL положительный и в ячейке R (m) хранится число,

отличное от нуля, — вычисление разности Да, — а — а,_г,

опреде­

ление знака Ад,- и сравнение его с хранимым в ячейке R ( ± ) знаком. Если знаки Ааг и й ( ± ) одинаковые, то управление передается счет-

134

чику i + 1 г> N. Когда знаки разные, отбиваются границы пласта: подошва при знаке «плюс» и кровля при знаке «минус».

3. После отбивки подошвы пласта вычисляется глубина по­ дошвы za V и определяется р т а х .

Рис. 52. Блок-схема про­ граммы полуавтоматиче­ ской отбивки границ пластов п отсчета суще­ ственных значений кажу­ щегося сопротивления

по данным БКЗ.

1

Восстановление R(m) и перемена знака

Я(к/л)

Да

Коней.

Пет

Ma

О/— R(m)

Знак A Q

Нет

Крааля

Да Т Подошва Сч ($)*•/

Вычисление z,

Уточнение

Начало

I

Присвоение знака „ ми­ нус"числам с харак­

терной точкой

Нет

Кровля

Подошва

Вычисление

Определение HgM

I

Определение

Поправка на

Печать 3

4. После выделения кровли пласта вычисляются глубина кровли

zK. V и мощность Нѵ

пласта, определяются мощность Нвм

и кажу­

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

р в м вмещающих пород, отсчитывается

среднее

значение кажущегося сопротивления пласта р о р и вносится поправка на его мощность Ну.

135

р т і п )
р э к с
р т а к

 

5. Если

отсчеты p m a x

и р С | ,

сделаны не для всех зондов (/ +

1 =S

где

I — число зондов),

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

уточнение

подошвы

пласта

по

кривой

следующего

зонда

и берутся

отсчеты р с р

и p D M .

 

Глава VIII

ОПРЕДЕЛЕНИЕ УДЕЛЬНОГО СОПРОТИВЛЕНИЯ ПЛАСТОВ

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

Определение удельного сопротивления пластов состоит из сле­ дующих этапов:

1)

выделение пластов п отбивка их границ;

2)

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

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

3)

уточнение удельного сопротивления бурового раствора;

4)

определение удельных сопротивлений пласта и зоны проникно­

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

Выделение пластов и отбивка их границ рассмотрены в преды­ дущей главе.

19. ОПРЕДЕЛЕНИЕ СУЩЕСТВЕННЫХ З Н А Ч Е Н И Й КАЖУЩЕГОСЯ УДЕЛЬНОГО СОПРОТИВЛЕНИЯ ПЛАСТОВ

В зависимости от типа и длины зонда для определения удельного сопротивления пластов используют экстремальные и средние зна­ чения кажущегося удельного сопротивления [14, 20]. Экстремальные значения (максимальные и минимальные снимают с кривых КС, записанных градиент-зондами длиной, близкой или превышающей мощность пласта, и зондами бокового и индукционного каротажа. Средние значения р с р кажущегося удельного сопротивле­ ния отсчитывают с кривых КС, полученных градиент-зондами длиной, меньшей мощности пласта. По В. Н. Дахнову, при подсчете среднего сопротивления исключают зону экранирования на участке, равном длине зонда:

h+n

136

где

п — число точек в

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

экранирования;

к — порядковый номер

числа, соответствующего

максимуму кри­

вой

КС в подошве пласта высокого сопротивления.

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

Из анализа кривых КС, полученных на каротажном интеграторе ЭКСМ, установлено, что оптимальные значения р о п т для пластов

недостаточно большой мощности (H

3—4L,

где L — длина зонда)

не могут интерпретироваться

при помощи палеток Б К З . Чтобы рас­

ширить область применения

палеток

Б К З , по

результатам измере­

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

Поправочные коэффициенты на

мощность

 

 

#опт =

Р™/РсР;

#экс = Рк./Рэкс

(71)

сведены в

палетки

[5] и

таблицы [50], используемые при

машинной

обработке

данных

Б К З .

Несмотря на то, что палетки и таблицы по­

правочных коэффициентов построены для двуслойных кривых зон­ дирования, они могут быть применены для пластов с проникновением фильтрата бурового раствора. Получаемая при этом поправка оптимального значения на мощность не превышает в большинстве случаев 10% и максимального значения 15%. Это обстоятельство позволило использовать таблицы поправочных коэффициентов на мощность при определении удельного сопротивления пластов на ЭВМ. Для этого показания зондов длиной меньше 0,95/? приводят к показа­ ниям против бесконечно мощных пластов.

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

зондов Б К З

эти

пределы

следующие.

 

1. Пласт высокого сопротивления. Оптимальные сопротивления

отсчитывают

при

длине

зондов Lj

0,75ff;

показания зондов при

Lj >> 0,3/7 исправляют

на мощность. Максимальные сопротивления

отсчитывают

при 0,7511 <! L,- <J 0,95/7 и L t

• > # ; показания зондов

исправляют

при

Ь} << 0,95/7.

 

 

2. Пласт низкого сопротивления. Оптимальные сопротивления

отсчитывают при

Lj < 0,3/7

и L; - < 0,5/7 (рис. 53), если

 

 

 

 

Рвм/Рп ^

15.

(72)

137

Минимальные

 

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

 

отсчитывают при

длине

зондов

Lj £ > 0 , 5 #

и Lj

!>0,ЗЯ,

если условие (72) не выполняется.

При Lj

\> H

 

на показания градиент-зондов

большое

влияние

оказывает мощность Нвы

 

вмещающей среды (см. табл. 10).

 

Как видно из табл. 10, при малой мощности пласта и вмещающей

среды

(Lj

г > H +

0,5JyBJ

минимальные

значения

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

сильно искажены и использовать

их для интерпретации Б К З неце­

 

 

 

 

 

 

 

 

 

лесообразно.

Поэтому

для от­

 

 

 

 

 

 

 

 

 

счета

минимальных

значений

 

 

 

 

 

 

 

 

 

принято

следующее

ограниче­

 

 

 

 

 

 

 

 

 

ние:

Lj<H

+ 0,5HBU.

(73)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В тех случаях, когда

удель­

 

 

 

 

 

 

 

 

 

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

пласта близ­

 

 

 

 

 

 

 

Рвн/Рп

ко к удельному сопротивлению

 

 

 

 

 

 

 

вмещающих пород

п

= (0,8-f-

 

 

 

 

 

 

 

 

 

Рпс. »3. Изменение

положения

(L/H) точек

 

— 1,3) р в м 1 , вместо

экстремаль­

отхода

кривых оптимального

зондирования

 

ных

значений

можно

 

исполь­

и БКЗ с увеличением

отношения

р в м / р п -

 

зовать значения

р о п

 

против

Шифр

кривых — относительная

мощность

 

т

 

 

пластов

H/d.

 

 

 

пласта

без

исправления на

 

 

 

 

 

 

 

 

 

мощность. Это правило

можно

распространить

на

пласт

с неглубоким

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

раствора

(рис.

54) для отсчета существенных значений

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

с кри­

вых КС, записанных большими зондами (вариант 4 отбивки границ1 ).

 

 

 

 

 

 

 

 

 

Т а б л и ц а

10

Влияние конечной мощности вмещающей среды р п м / Р с = 100

 

на минимальные

значения

кажущегося удельного

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

 

в

пласте

низкого

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

п /Рс = 10; H/d =

8)

 

H

=

2'id

 

 

 

 

 

 

Я в м = с о

 

DM

1

 

 

 

 

 

 

 

 

Lj/d

Pmirr 0 m

- m

6, %

 

Lj/d

Pmin- 0 m ' m

6, %

Pmin- ° M

l M

H

 

18,20

 

1

 

H

19.80

10

18

 

# + 0 , 3 3 # в м

 

30,80

 

3

#

+ 0,5tfB M

25,00

17

30

 

# - г - # в м

 

1,44

 

2000

 

 

1,05

2000

28

 

П р и м е ч а й

н e. ö =

( P „ , l n / P m i n )

100.

 

 

 

 

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

1 Варианты отбивки границ см. главу V I I , раздел 18.

138

- 4. Несовпадение границ пласта высокого сопротивления, опре­ деленных по малым и большим градиент-зондам (вариант 2 отбивки границ х ) . Такое несовпадение границ может быть встречено в про­ дуктивном пласте с водо-нефтяным контактом с большим влиянием зоны проникновения на показания малого зонда. В этом случае гра­ ница пласта, полученная по данным большого зонда, обычно совпа­ дает с ВНК. Граница, отбитая по показаниям малого зонда, может наблюдаться ниже. Применяя правила отсчета существенных зна­ чений КС против пластов низкого и высокого сопротивлений, при несовпадении границ существенные значения р о п т и р э к с для зондов

AôdMOJdN

Al6dMÜ,5N A3ldM0,5N

і і з

ч о і о і

ом-м

Рис. 54. Кривые КС (а) и БКЗ (б), полученные

на

каротажном

интеграторе

ЭКСМ для пласта с удельным сопротивлением р п =

р п зоной

проникнове­

ния р з п = 25 р с н D =

2d.

 

 

 

2 — пласт с H =

16d; 2 — вмещающая среда; з — крест скважины.

с Lj ^

L M необходимо

определить с учетом границ, полученных по

показаниям

малого основного

зонда длиной

L M , а для зондов

с Lj > L 6

— с

учетом границ,

полученных

большим основным

зондом

длиной

L 6 .

 

 

 

5. Ступенчатое изменение удельного сопротивления пласта, когда прослой высокого сопротивления расположен со стороны сближенных электродов зонда (вариант 5 отбивки границ х ) или со стороны уда­ ленного электрода (вариант 6 отбивки границ 1 ) . В первом случае против промежуточного прослоя следует отсчитывать оптимальные

значения р 0 П т , если длина зондов Lj «s 0,5Н3

3 мощность про­

межуточного слоя). Кривые КС, полученные

зондами с L; ->>0,5#„

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

Pmin при

Lj<$0,9H3.

 

 

В л и я н и е в м е щ а ю щ е й с р е д ы н а к а ж у щ е е с я

у д е л ь н о е

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

Точность

результатов

 

X

 

 

1 Варианты

отбивки границ см. главу V I I , раздел 18.

 

13S

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