Файл: Сохранов Н.Н. Машинные методы обработки и интерпретации результатов геофизических исследований скважин.pdf

Скачать файл
Заказать решение

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 132

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ш

г

н

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

Смотрите также файлы