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

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

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

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

Добавлен: 11.04.2024

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

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

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

Из рассмотрения описанных алгоритмов можно рекомендовать

следующую

схему

редактирования цифровых

каротажных данных

на ЭВМ.

 

 

 

1.

Ввод

данных

в ЭВМ.

 

2.

Обработка контрольных записей нуль-

и стандарт-сигналов

с целью уточнения масштаба записи (цены разряда цифровых дан­

ных) и положения нуля.

 

3.

Увязка меток глубины и числа слов (циклов) между ними.

4.

Распаковка слов, исправление ошибок в записи отдельных чисел

(сбоев

преобразования, записи и ввода) и увязка данных по глубине

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

записи геофизических величин.

По

окончании редактирования

числовых массивов M j геофизиче­

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

для

печати

(табл.

5).

 

 

 

 

Т а б л и ц а 5

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценка

качества отредактированных геофизических данных

 

 

Характер погрешности (штрафа)

 

Число

Процент

Оценка (балл)

 

 

 

ошибок

ошибок

Ошибка

установки

масштаба

 

 

0

3

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

2

 

1

Расхождение

числа

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

меток,

0

 

3

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

Sa 20

1

Сбои

сннхронизащтп

при записи слова

(цик­

0

3=0,2

3

ла)

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

Ошибка

измерения

 

(записи)

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

0

 

3

величины

 

 

 

 

 

 

 

3=0,2

2

 

 

 

 

 

 

 

 

 

 

> 1

1

Несогласование

измерений с

контрольными

0

 

3

данными

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

> 2 0

1

Общая

оценка

качества

 

 

Средний

 

 

 

 

 

 

 

 

 

 

 

балл

14. ТРАНСФОРМАЦИЯ ГЕОФИЗИЧЕСКИХ Д А Н Н Ы Х

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

102


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

Поэтому

трансформация

заключается

в следующем:

 

 

а)

 

выделение всех прослоев по экстремумам или точкам перегиба;

б) отсчет существенных значений против

прослоев;

 

 

в)

 

исправление

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

значений;

 

 

 

 

г) присвоение исправленных значений всем точкам в интервале

прослоев и запись трансформированных

данных по точкам в форме,

принятой для исходных геофизических дан­

 

 

 

 

ных:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l j 11' a

j

12'

 

 

а/1

a

211

u Ci;

 

., a

12m2i

 

 

 

 

 

 

j'lmn

j

j

22i

 

 

 

 

 

 

 

 

 

 

«/Vil

 

 

 

 

 

 

 

 

 

(33)

 

 

 

 

где

ö

/ ll

Я/12'

 

.,

=

a/i m ,

— исправлен-

 

 

 

 

ное существенное значение для первого

 

 

 

 

прослоя;

а/ Ѵ і

= ауѵ

 

 

=

аjvmv

 

 

ТО ЖѲ,

 

 

 

 

для ѵ-го прослоя;

 

j

— номер

параметра;

 

 

 

 

ту — число

точек

в ѵ-м прослое.

 

 

 

 

 

 

 

 

Возможна

более

плотная

 

(экономичная)

 

 

 

 

форма

записи

трансформированной

 

кривой:

 

 

 

 

 

Й Д> тѵ> а /2і т

 

 

;

ßjy тпѵі • • •

(34)

 

 

 

 

Для выделения прослоев необходимо при­

 

 

 

 

менять способы отбивки

границ,

автомати­

 

 

 

 

чески

обеспечивающие

фильтрацию

помех

 

 

 

 

измерения

и

регистрации

(например,

из-

 

 

 

 

зубренность кривой). Таким способом яв­

 

 

 

 

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

гра­

Рис. 40.

Примеры

отсчета

ниц.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

средних

значений о

кривых

 

 

 

прослоев

в

зависимости

от

бокового

(а) и

радиоактив­

Против

ного (б) каротажа.

вида каротажа и мощности прослоя

(пла­

 

 

 

 

ста)

отсчитывают

средние

значения

а с р ,

экстремальные

(макси­

мальные

а т а х ,

минимальные

а т і

п ) .

Средние значения

отсчиты­

вают в прослое, если мощность его больше граничного значения НТ.

Величину

Нг устанавливают для каждого вида каротажа.

Среднее

показание

обычно

вычисляют

не во всем

интервале

прослоя

Z K z n

( Z K '

z n глубины соответственно кровли

и подошвы), а в той

части

его

z'K—z'n,

где влияние

вмещающих

пород незначительно

(с учетом зоны экрана). Примеры отсчета средних значений для кри­ вых Б К и НГК показаны на рис. 40. Экстремальные значения от­ считывают при мощности прослоя меньше граничного значения.

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

103


прослоях (# - < Нт). Исправление существенных значений целесо­ образно производить умножением их на поправочные коэффициенты Kj, определяемые по кривым зависимости Kj — ф х, С2), где Сг и С2 — величины, влияющие на результаты геофизических измерений.

По исправленным существенным значениям (а'ср, a'max и a'min) определяют приведенные' значения геофизических величин по фор­ муле

 

I а'іч a / l

I

 

°/v = п. -a-

 

'

( 3 5 )

 

— "vi

 

 

где

a;-v — исправленное существенное

значение

против ѵ-го пласта;

ajx

и ay, — исправленные существенные значения против опорных

пластов, например против глин (ау 1 ) и плотных (или чистых) пород (а/г)-

Приведенные существенные значения присваивают всем точкам в интервале пласта и записывают эти значения по форме (33) или (34).

15. У В Я З К А Д А Н Н Ы Х КАРОТАЖА ПО Г Л У Б И Н Е

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

1. Увязка каротажных кривых взаимной корреляцией. Способ

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

[46] заключается

в вычислении

функ­

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

k+n

 

 

 

 

 

 

 

 

 

R (ô)0 , -

(

3

6

)

4 ' 0 l

f h+n

k+n

 

 

4 '

'i=/t i=k

для двух последовательностей чисел: аоі (основной кривой) и а(согласуемой кривой). Одна из них сдвигается относительно другой каждый раз на величину шага дискретизации А по глубине до полу­ чения максимума функции взаимной корреляции. Число сдвигов ô = ±1 одной последовательности чисел относительно другой до положения максимума функции взаимной корреляции соответствует величине смещения Az(- = Al кривых по глубине. Эта величина учи­ тывается при интерпретации. Интервал корреляции (количество чисел ІѴК) обычно выбирается в зависимости от характера коррели­ руемых кривых. Для Туймазинского месторождения этот интервал достаточно взять равным 20—25 м (7ѴК = 100—125).

Этот способ увязки промыслово-геофизических данных был опро­

бован на фактическом материале. Результаты опробования в скв. 718

104


ПІкапово (см. рис. 33) при А = 0,2 м приведены в табл. 6. Как видно, величины смещения кривых разнородных геофизических величин определяются методом взаимной корреляции со значительными ошибками (до 2 м). Из этого следует, что способ взаимной корреля­ ции пригоден для увязки по глубине кривых только однородных геофизических величин (ПС, КС, записанных градиент-зондами с АО < Я і др.) и не годится для увязки по глубине в общем случае.

 

 

 

 

 

Т а б л и ц а 6

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

 

относительно кривой КС А0,8М0,Ш

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

 

 

(Д = 0,2

ы)

 

 

 

 

 

Смещение в величинах Д

Ошибка в величинах

 

 

 

определен!тое на ЭВМ

 

 

 

Вид измерения

фактиче­

без транс-

с транс­

без транс­

с транс­

 

 

ское

формации

форма­

 

 

 

формацші

форма­

 

цией

 

 

 

цией

 

 

 

 

 

 

 

 

К С ,

з о н д A4M0.5N

6

5

4

1

2

К С ,

з о н д В7,5А0,75М

6

10

3

4

3

КС,

з о н д АБК - 3

—4

- 1

- 5

- 3

—1

К С, м н к р о з о и д

2

8

3

6

1

ПС

 

0

10

2

10

2

Н Г К

 

0

7

0

7

0

КС,

з о н д А0,4М0,Ш

2

2

 

0

 

2. Увязка по глубине трансформированных кривых. Результаты определения величины смещения кривых при помощи

взаимной корреляции улучшаются, если кривые предварительно трансформировать в графики прямоугольной формы (см. рис. 33). Погрешность определения смещения (увязки кривых по глубине) в этом случае значительно меньше (см. табл. 6). Это показывает, что увязку геофизических данных по глубине на ЭВМ можно осуще­ ствлять с небольшой погрешностью (не более 2А) путем корреляции трансформированных кривых по формуле (36) или сопоставлением последовательности глубин z0 ѵ и zjv границ пластов, выделенных соответственно по основной и привязываемым к ней кривым. Для не­ которых каротажных кривых можно сопоставлять глубины харак­ терных точек (экстремумов, перегибов и др.).

3. Корреляция границ пластов.

глубин z o v и zjv

 

Алгоритм корреляции последовательностей

гра­

ниц пластов заключается в следующем:

 

 

 

 

а) определяют положения границ пластов

(характерных

точек)

на коррелируемых кривых и

запоминают их

глубины

z0 ѵ

и

Zjv;

б) для каждой границы z0 ѵ ,

выделенной ка основной

кривой, на­

ходят идентичную ей границу z;-v- на другой кривой, расположенной

105


к ней ближе других в

интервале z o V

± Az. Величина Az задается

в зависимости от предполагаемого смещения кривой;

 

 

в) вычисляют

абсолютные

разности

глубин | z o v

Z/v'| для

всех

идентичных пар

границ

и их

сумму

 

 

 

 

к

 

 

 

 

 

•s(ô) = 2l2ov - (z/ v + A-o)|,

ô = o , i , . . . ,

I;

(37)

здесь к — число пар

границ.

 

 

 

Сумму вычисляют после каждого сдвига границ z;-ѵ< на величину шага квантования до получения минимума суммы S (ô) при б = ± Z. Величина сдвига I последовательности глубин границ zjv относи­ тельно последовательности глубин границ z0 ѵ основной кривой определяет величину смещения кривых по глубине.

4. Увязка каротажных данных по опорной кривой.

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

16. ОЦЕНКА КАЧЕСТВА Д А Н Н Ы Х Б К З И ДРУГИХ ЗОНДОВ КАРОТАЖА СОПРОТИВЛЕНИЙ

Точность результатов интерпретации материалов бокового каро­ тажного зондирования (БКЗ) в значительной степени определяется качеством исходных данных. Перед обработкой материалов Б К З необходимо проверить их качество: установить погрешность опреде­ ления удельного сопротивления бурового раствора и в случае боль­ шого ее значения уточнить это сопротивление, выявить кривые со­ противления с низким качеством измерений (неудовлетворительные замеры) и, если невозможно их исправить, то исключить из обработки.

При определении на ЭВМ удельного сопротивления пластов по данным Б К З и комплекса измерений разнотипными зондами проверка качества исходных данных имеет еще большее значение, чем при ручной обработке. Способы контроля качества в обоих случаях могут быть одними и теми же.

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

на