RU2144683C1

RU2144683C1 - Способ обработки сейсмического сигнала и разведки месторождений

Info

Publication number: RU2144683C1
Application number: RU96115275/28A
Authority: RU
Inventors: Майкл С. Бахорич (US); Майкл С. Бахорич; Стивен Л. Фамер (US); Стивен Л. Фамер
Original assignee: Амоко Корпорейшн
Priority date: 1994-12-12
Filing date: 1995-10-05
Publication date: 2000-01-20
Also published as: WO1996018915A1; US5838564A; CN1121619C; CN1138902A; AU4133396A; CA2179901A1; NO962731D0; EG20609A; NO962731L; EP0736185B1; EP0736185A1; MX9603026A; US5563949A; AU696742B2; NO311316B1; CA2179901C

Abstract

Использование в сейсмической разведке месторождений для идентификации структурных и стратиграфических характеристик в трех измерениях. Способ разведки углеводородов предусматривает следующие операции: получение комплекта сигналов сейсмических трасс, распределенных по заранее определенному трехмерному объему толщи земли, разделение трехмерного объема на множество вертикально смещенных и главным образом удаленных друг от друга горизонтальных срезов и разделение по меньшей мере одного из указанных срезов на множество ячеек, причем каждая из указанных ячеек имеет участки по меньшей мере трех локализованных в ней сейсмических трасс, измерение поперек каждой из указанных ячеек кросс-корреляции между одной парой трасс, лежащих в первой вертикальной плоскости, для получения значения инлайн (in-line), и измерение кросс-корреляции между другой парой трасс, лежащих в другой вертикальной плоскости, для получения значений кросслайн (cross-line), комбинирование указанного инлайн значения с указанным кросслайн значением для получения значения когерентности для каждой из ячеек и визуализацию значений когерентности. Изобретение позволяет значительно ускорить составление карт структурных деталей и выявить детали с ложными соотношениями. 5 с. с. и 26 з.п. ф-лы, 6 ил.

Description

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

Обычные 2-D сейсмические данные получают вдоль линий (см. линии 10 и 11 на фиг. 1), которые образованы решетками сейсмоприемников на участке побережья до прибойной зоны или решетками гидрофонных сейсмоприемных кос (стримеров), которые пересекают прибрежную водную зону. Сейсмоприемники и гидрофоны работают в качестве датчиков и принимают энергию, которая передана вглубь земли и отражена назад к поверхности земли от поверхностей раздела 12 скальной породы нижнего горизонта. Энергия обычно создается на поверхности земли при помощи вибросейсмических аппаратов, которые передают с поверхности импульсы встряхивания грунта с заранее определенными интервалами и частотами. При работе с водной поверхности для этой цели обычно используют пневматические пушки. Тонкие изменения в распределении энергии, возвращенной на поверхность, часто отражают вариации стратиграфического, структурного и жидкостного содержания подземных резервуаров.

При использовании 3-D сейсмических данных принцип остается таким же, однако линии и решетки размещают более тесно (см. фиг. 1 и 2), что позволяет получать более детальное перекрытие нижнего горизонта. При таком перекрытии с очень высокой плотностью требуется записывать, хранить и обрабатывать чрезвычайно большие объемы цифровых данных, перед тем, как может быть произведена их окончательная интерпретация. Обработка данных требует применения огромных компьютерных ресурсов и сложного программного обеспечения, чтобы усилить сигнал, отраженный от нижнего горизонта, и выделить его из сопровождающих шумов, которые маскируют сигнал.

После того, как произведена обработка данных, бригада геофизиков производит компиляцию и интерпретацию сейсмической 3-D информации в форме 3-D куба (см. фиг. 4), который эффективно представляет характеристики нижнего горизонта. При использовании куба данных информация может выводиться на индикацию в самом разном виде. Могут быть сделаны карты горизонтальных временных срезов на выбранной глубине (см. фиг. 5). При использовании компьютерных рабочих станций интерпретатор может делать срезы по полю для исследования выходов резервуара на различных горизонтах. Могут быть также произведены вертикальные срезы или сечения в любом направлении, с использованием сейсмических или скважных данных. Временные карты могут быть конвертированы в глубину для получения структурной интерпретации на конкретном уровне.

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

Так как 3-D сейсмические данные используются для разведки и интерпретации данных, то возникает необходимость в их усовершенствовании.

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

В соответствии с настоящим изобретением предлагается способ разведки углеводородов. Этот способ содержит следующие операции: получение комплекта сигналов сейсмических трасс, распределенных по заранее определенному трехмерному объему толщи земли; разделение трехмерного объема на множество вертикально смещенных и главным образом удаленных друг от друга горизонтальных срезов; разделение каждого из срезов на множество ячеек, которые подразделены на идущие в боковом направлении ряды и колонки, таким образом, что участки по меньшей мере трех главным образом вертикально идущих сейсмических трасс локализованы в них; измерение поперек каждой ячейки кросс-корреляции между одной парой трасс, лежащей в одной вертикальной плоскости, для получения значения инлайн (in-line), и измерение кросс-корреляции между другой парой трасс, лежащей в другой вертикальной плоскости, для получения значения кросслайн (cross-line), которые служат оценкой временного (по времени) падения (спада) в инлайн и кросслайн направлениях; комбинирование инлайн значения с кросслайн значением для получения значения когерентности для каждой ячейки; и визуализация значений когерентности поперек по меньшей мере одного горизонтального среза.

Предложенная техника особенно хорошо подходит для интерпретации ложных планов (плоскостей дефектов) в пределах 3-D сейсмического объема и для обнаружения тонких стратиграфических характеристик в 3-D (информации). Это происходит потому, что сейсмические трассы, перерезанные ложной линией (линией дефекта), обычно имеют другой сейсмический характер, чем трассы на любой стороне от ложной. Измерение похожести трасс (то есть когерентности или непрерывности 3-D) вдоль временного среза позволяет выявить контуры низкой когерентности вдоль этих ложных линий. Такие значения когерентности позволяют выявить критические детали нижнего горизонта, которые не всегда ясны на традиционных сейсмических сечениях. За счет вычисления когерентности вдоль серии временных срезов можно при помощи ложных контуров идентифицировать ложные планы или поверхности.

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

На фиг. 1 показано расположение сейсмоприемников для получения 3-D сейсмических данных из нижнего горизонта земной толщи, необходимое для проведения обработки в соответствии с настоящим изобретением. На фиг. 2 показан вид сверху расположения фиг. 1. На фиг. 3 приведено расположение сейсмических трасс, лежащих в плоскости, проходящей через ряд сейсмоприемников, показанных на фиг. 2. На фиг. 4 приведено графическое представление информации, полученной при обработке 3-D данных. На фиг. 5 приведено графическое представление горизонтального временного среза 3-D сейсмических данных, обработанных в соответствии с известным состоянием техники. На фиг. 6 приведено графическое представление горизонтального временного среза 3-D сейсмических данных, обработанных в соответствии с настоящим изобретением.

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

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

Следующей операцией является выработка "куба нарушения сплошности". Это производится приложением алгоритма когерентности к 3-D сейсмическим данным. Этот алгоритм может приобретать различные формы. Вне зависимости от выбранной формы его функционирование заключается в сравнении аналогичности прилегающих друг к другу районов сейсмических данных внутри 3-D сейсмического объема. Если сегмент трассы аналогичен соседнему (например, в инлайн или кросслайн направлениях), то ему присваивается значение низкого нарушения сплошности; если сегмент трассы не аналогичен соседнему, то ему присваивается значение высокого нарушения сплошности.

На фиг. 2 приведен вид в плане участка 3-D сейсмического объема. Для того, чтобы измерить нарушение сплошности, участок трассы A сравнивают со смежными участками трассы B и C. Один из путей вычисления аналогичности трасс описан ниже.

Нулевая средняя задержанная кросс-корреляция в направлении инлайн (в x-направлении) между трассой u(t, x, y) и u(t, x+dx, y) при временной задержке "tlag" мс определяется в соответствии с выражением:
(1)

в котором

и

являются автокорреляциями, использованными для нормализации кросс-корреляции, и где W+W является длиной в мс корреляционного окна. Важно выбрать W достаточно широким для того, чтобы оценка среднего нуля была действенной. Могут быть использованы и другие способы нормализации (например, произведение энергий трасс и т.п.). Приемлемыми являются значения порядка сейсмических импульсоидов. В частности, кросс-корреляция является одним из способов комбинирования двух колебаний для измерения аналогичности колебаний. Автокорреляция является способом комбинирования колебания с самим собой. См. публикацию Шерифа. Энциклопедический словарь геофизической разведки. Выпущена обществом геофизиков-разведчиков в г.Тулза, штат Оклахома, США.

Нулевая средняя задержанная кросс-корреляция в направлении кросслайн (в y-направлении) между трассой u(t, x, y) и u(t, x, y+dy) при временной задержке "tlag" мс определяется в соответствии с выражением:

где

Направление видимого временного падения в направлениях x и y считают той задержкой (то есть tlagx и tlagy), которая имеет самую большую (то есть наиболее положительную) кросс-корреляцию. Эти значения равны p_x(t, tlagx) и p_y(t, tlagy).

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

Это значение (или аттрибутив) служит в качестве достаточно надежной оценки нарушения сплошности сигнала в геологической формации, а также нарушения сплошностей сигнала через ложные и ошибочные несоответствия.

Компьютерная программа
Далее приводится упрощенная программа ФОРТРАН 77 для осуществления указанных вычислений.

Пусть имеется трасса "x" из 3-D сейсмического амплитудного объема и две соседних трассы "y" (в направлении инлайн) и "z" (в кросслайн направлении). В этом случае стандартная подпрограмма COH позволяет вычислить выходную трассу "rho", содержащую коэффициенты когерентности, с использованием алгоритма прогонки окна кросс-корреляции, где
"mins" и "maxs" представляют собой минимальные и максимальные показатели образца для всех четырех трасс;
"inwinl" является длиной окна в образцах;
"nlags" задает число задержек (относительных временных сдвигов), которые нужно сделать с каждой стороны "0" в кросс-корреляции; и
"sr" является интервалом образца (выборки) в мс.

При каждой выборке стандартная подпрограмма CROSS производит вычисление серии нормализованных коэффициентов кросс-корреляции, возвращая самые большие коэффициенты для каждого направления в "rho1" и "rho2". Временной сдвиг, при котором случаются максимальные коэффициенты, возвращается в "tshf1" и "tshf2"; эти времена не используются. Стандартная подпрограмма COH вызывается повторно, один раз для каждой трассы во входном сейсмическом амплитудном объеме, для получения объема 3-D данных или "куба когерентности", содержащего коэффициенты когерентности.

subroutine coh (x, y, z, rho, mins, maxs, iwinl, nlags, sr)
real x(mins:maxs), y(mins:maxs), x(mins:maxs)
real rho(mins:maxs)
ihwin = iwin1/2
doj = mins+ihwin, maxs-ihwin
k = j - ihwin
call cross (x(k), iwinl, y(k), iwinl, nlags, sr, rho1, tshf1)
call cross (x(k), iwinl, z(k), iwinl, nlags, sr, rho2, tshf2)
rho(J) = sqrt (rho1^*rho2)
enddo
return
end
subroutine cross (x, nx, y, ny, lags, sr, peak, tshift)
real x(0:nx-1), y(0:ny-1), sr, peak, tshift
parameter (maxlags=128)
real g(-maxlags:+maxlags)
double precision xx,yy
nlags = max(0, min(lags, maxlags))
tshift = 0.0
peak = 0.0
xx = 0.0
yy = 0.0
ks = 0
do ix = 0, nx-1
xx = x(ix)^**2+xx
enddo
if (xx.eq. 0.0) return
do iy = 0, ny-1
yy = y(iy)^**2+yy
enddo
if(yy.eq.0-0)return
do is = -nlags, +nlags
g(is) = 0.0
do it = 0, nx-1
if (it-is.ge.0) then
if (it-is.le.ny-1)then
g(is) = g(is) + x(it)^*y(it-is)
endif
endif
enddo
if (abs(peak).lt.abs(g(is))) then
peak = g(is)
ks = is
endif
enddo
tshift=ks^*sr
Могут быть использованы, например, рабочие станции интерпретации результатов типа Landmark или GeoQuest, позволяющие интерпретировать ложные значения и стратиграфические характеристики путем загрузки куба нарушения сплошности в качестве сейсмического объема. Может быть использовано программное обеспечение визуализации (например, Landmark's SeisCube) для быстрого выполнения срезов через объем нарушения сплошности с целью установления комплексных ложных соотношений. Визуализация нарушения сплошности позволяет сократить время цикла интерпретации при ее использовании избирательно относительно подлежащих интерпретации сейсмических линий, что позволяет исследователю работать в зонах вблизи ложных и слабых данных. Кроме того, при этом могут быть просто идентифицированы и интерпретированы тонкие стратиграфические характеристики и комплексные ложные сигналы, которые трудно обнаружить при традиционной сейсмической визуализации. На фиг. 5 и 6 рядом произведено сравнение одной и той же сейсмической информации, которая выводится на индикацию и обрабатывается при помощи как известных ранее методов, так и в соответствии с настоящим изобретением. Ложные линии четко видны на фиг. 6.

Карты когерентности были получены для множества 3-D разведок. При глубинах приемлемого качества данных, ориентировочно 90% ошибок (ложных сигналов) легко могут быть обнаружены. На картах когерентности были обнаружены ложные сигналы, которые плохо видны на сейсмических сечениях, но четко присутствуют на картах когерентности благодаря надежности предложенного способа и перспективе карт для ложных картин (картин распределения дефектов). Так как карты когерентности могут быть получены для не интерпретированных временных срезов, то настоящее изобретение предлагает средство для значительного ускорения составления карт структурных деталей и для выявления деталей с ложными соотношениями, которые в противном случае могли бы быть выявлены только в результате сложной утомительной работы.

Специфические примеры
2-D карты когерентности были разработаны вдоль захваченных горизонтов и позволили четко идентифицировать сланцевые диапиры в прибрежных водах Нигерии.

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

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

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

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

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

Дальнейший анализ
Куб нарушения сплошности позволяет четко выявить ложные планы как зоны высокого нарушения сплошности. Однако эти зоны не обязательно являются зонами низкого отношения сигнал/шум. Способ усиления этих ложных зон предусматривает применение так называемого "среднего планарного оператора". Дефекты в нижнем горизонте земли обычно выражаются сами по себе как плоскости или поверхности. В случае изогнутой поверхности дефекта ряд небольших плоскостей может быть использован для аппроксимации поверхности дефекта. В соответствии с этим аспектом настоящего изобретения, небольшой планарный оператор используется для усиления (то есть в качестве "фильтра") идентификации тонких стратиграфических характеристик. Прежде всего, выбирают небольшой участок сейсмических данных вокруг центрального значения.

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

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

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

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

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

Claims

1. Способ локализации подземных характеристик, дефектов и контуров, включающий в себя: а) получение трехмерных сейсмических данных, перекрывающих заранее определенный объем толщи земли, b) разделение указанного объема на решетку относительно малых трехмерных ячеек, отличающийся тем, что каждая трехмерная ячейка характеризуется наличием по меньшей мере трех разделенных в боковом направлении и главным образом вертикальных сейсмических трасс, локализованных в ней, и тем, что проводят с) измерение в каждой указанной ячейке подобия указанных по меньшей мере трех трасс относительно двух заранее определенных направлений и d) визуализацию указанного подобия трасс указанных ячеек для образования двухмерной карты подземных характеристик.

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

3. Способ по п. 2, отличающийся тем, что операция (1) включает в себя операцию нормализации каждого значения инлайн и каждого значения кросслайн.

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

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

6. Способ по п. 2, отличающийся тем, что операция (1) включает в себя операции вычисления средней нулевой задержанной кросс-корреляции в указанном направлении инлайн и вычисления средней нулевой задержанной кросс-корреляции в указанном направлении кросслайн.

7. Способ по п. 6, отличающийся тем, что операция (2) включает в себя операции идентификации наиболее положительного значения указанной средней нулевой задержанной кросс-корреляции в указанном направлении инлайн и идентификации наиболее положительного значения указанной средней нулевой задержанной кросс-корреляции в указанном направлении кросслайн.

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

9. Способ по п.1, отличающийся тем, что в операции (с) указанные заранее определенные направления являются взаимно перпендикулярными, причем указанное подобие указанных ячеек измерено как функция кросс-корреляции между двумя трассами в первом направлении и как функция кросс-корреляции между двумя трассами в направлении, которое перпендикулярно указанному первому направлению.

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

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

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

13. Изделие, используемое при сейсмической разведке, при которой трехмерные сейсмические данные, содержащие отраженную сейсмическую энергию, регистрируются как функция времени для создания серии сейсмических трасс, и когда используется компьютер, адаптированный к обработке таких сейсмических трасс, отличающееся тем, что оно включает в себя среду, которая может быть считана компьютером и которая содержит команды для указанного компьютера для осуществления способа, содержащего следующие операции: а) определение трехмерных сейсмических данных для заранее заданного объема, b) сравнение подобия прилегающих районов указанных трехмерных сейсмических данных для указанного объема путем: 1) разделения указанного объема на по меньшей мере один горизонтальный срез и разделения указанного по меньшей мере одного среза на множество ячеек, которые подразделены на идущие в боковом направлении ряды и колонки, причем каждая из указанных ячеек имеет участки по меньшей мере трех сейсмических трасс, простирающихся через них, включающих в себя первую трассу и вторую трассу, которые лежат в одной плоскости, и третью трассу, которая вместе с указанной первой трассой лежит в другой плоскости, которая находится главным образом под прямым углом к указанной первой плоскости, 2) измерения поперек каждой указанной ячейки подобия указанных трасс, лежащих в указанной первой плоскости, для получения первого значения подобия, и измерения подобия указанных трасс, лежащих в указанной второй плоскости, для получения второго значения подобия, 3) комбинирования указанного первого значения подобия и указанного второго значения подобия в каждой ячейке для получения одного значения подобия, которое является репрезентативным подобием для указанных сейсмических трасс в указанной ячейке, и 4) запоминания указанного одного значения подобия каждой ячейки указанного по меньшей мере одного горизонтального среза в форме, позволяющей визуализировать указанные комбинированные значения подобия как сейсмическую атрибутивную карту.

14. Изделие по п.13, отличающееся тем, что указанная среда содержит команды для указанного компьютера для осуществления операции (2) путем определения кросс-корреляции между указанной первой трассой и указанной второй трассой, и путем определения кросс-корреляции между указанной первой трассой и указанной третьей трассой.

15. Изделие по п.14, отличающееся тем, что указанная среда содержит команды для компьютера для осуществления операции (2) путем определения средней нулевой задержанной кросс-корреляции между указанной первой трассой и указанной второй трассой, и путем определения средней нулевой задержанной кросс-корреляции между указанной первой трассой и указанной третьей трассой.

16. Изделие по п. 15, отличающееся тем, что указанная среда содержит команды для компьютера для осуществления операции (3) путем идентификации наиболее положительной из каждых из указанных средних нулевых задержанных кросс-корреляций.

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

18. Устройство для сейсмической разведки, которое включает в себя: а) среду с записью, которая может быть считана компьютером и которая содержит команды для обработки информации, отличающееся тем, что оно предусматривает проведение следующих операций: (1) обеспечение доступа к набору данных, который включает в себя сейсмические сигнальные трассы, распределенные по заданному трехмерному объему толщи земли, (2) разделение указанного трехмерного объема на множество распределенных вертикально и главным образом смещенных друг от друга горизонтальных срезов и разделение по меньшей мере одного из указанных срезов на множество ячеек, которые подразделены на идущие в боковом направлении ряды и колонки, причем каждая из указанных ячеек имеет участки по меньшей мере трех сейсмических трасс, локализованных в ней, при этом каждый из указанных участков указанных трасс проходит главным образом через указанные ячейки, причем первая трасса и вторая трасса в указанной ячейке лежат в первой плоскости, а третья трасса и указанная первая трасса лежат в другой плоскости, которая главным образом расположена под прямым углом, к указанной первой плоскости, (3) вычисление в каждой из указанных ячеек кросс-корреляции между указанными трассами, лежащими в указанной первой плоскости, для получения значения инлайн, и вычисление кросс-корреляции между указанными трассами, лежащими в указанной другой плоскости, для получения значения кросслайн, причем эти значения представляют собой оценку падения по времени (временного падения) в направлении инлайн и в направлении кросслайн, и (4) комбинирование указанных значений инлайн и кросслайн для получения значения подобия для каждой из указанных ячеек.

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

20. Устройство по п.18, отличающееся тем, что операция (4) дополнительно предусматривает операцию вычисления средней нулевой задержанной кросс-корреляции в указанном направлении инлайн и средней нулевой задержанной кросс-корреляции в указанном направлении кросслайн.

21. Устройство по п.20, отличающееся тем, что операция (4) дополнительно предусматривает операции идентификации наиболее положительного значения указанной средней нулевой задержанной кросс-корреляции в указанном направлении инлайн, и идентификации наиболее положительного значения указанной средней нулевой задержанной кросс-корреляции в указанном направлении кросслайн.

22. Устройство по п.21, отличающееся тем, что операция (4) дополнительно включает в себя операцию вычисления среднего между указанными двумя наиболее положительными значениями.

23. Устройство для сейсмической разведки, которое включает в себя: а) среду с предварительной записью, которая может быть считана компьютером и которая содержит команды для обработки информации, отличающееся тем, что оно предусматривает проведение следующих операций: (1) выборку из памяти трехмерных сейсмических данных, которые перекрывают заданный объем толщи земли, (2) цифровую сортировку указанных данных с получением решетки относительно малых трехмерных ячеек, причем каждая из указанных ячеек имеет по меньшей мере три смещенные в боковом направлении и главным образом вертикальные сейсмические трассы, локализованные в ней, (3) вычисление в каждой из указанных ячеек значения подобия, исходя из данных указанных по меньшей мере трех трасс, относительно двух заданных направлений, (4) запоминание указанных значений подобия указанных ячеек в компьютере для вывода на индикацию двухмерной карты подземных характеристик, представленных указанными значениями подобия.

24. Устройство по п.23, отличающееся тем, что в операции (2) указанные два заранее заданных направления являются взаимно перпендикулярными, причем каждое значение подобия вычислено как функция кросс-корреляции между двумя трассами, лежащими в одном из указанных взаимно перпендикулярных направлений, и как функция кросс-корреляции между двумя трассами, лежащими в другом из указанных взаимно перпендикулярных направлений.

25. Устройство по п.24, отличающееся тем, что указанное значение подобия вычислено как функция наиболее кросс-корреляции в указанном первом направлении и наибольшей кросс-корреляции в указанном другом направлении.

26. Устройство по п.25, отличающееся тем, что указанное значение подобия в операции (3) является функцией геометрически среднего указанных двух наибольших кросс-корреляций.

27. Сейсмическая карта, подготовленная за счет обеспечения доступа компьютера к комплекту данных, который включает в себя сейсмические сигнальные трассы, распределенные по заданному трехмерному объему толщи земли, при использовании способа обработки указанных данных, который отличается тем, что включает в себя следующие операции: (1) разделение указанного трехмерного объема на множество распределенных вертикально и главным образом смещенных друг от друга горизонтальных срезов и разделение по меньшей мере одного из указанных срезов на множество ячеек, которые подразделены на идущие в боковом направлении ряды и колонки, причем каждая из указанных ячеек имеет участки по меньшей мере трех сейсмических трасс, локализованных в ней, при этом каждый из указанных участков указанных трасс проходит главным образом через указанные ячейки, причем первая трасса и вторая трасса в указанной ячейке лежат в первой плоскости, а третья трасса и указанная первая трасса лежат в другой плоскости, которая главным образом расположена под углом к указанной первой плоскости, (2) вычисление в каждой из указанных ячеек кросс-корреляции между указанными трассами, лежащими в указанной первой плоскости, для получения значения инлайн, и вычисление кросс-корреляции между указанными трассами, лежащими в указанной другой плоскости, для получения значения кросслайн, причем эти значения представляют собой оценку падения по времени в направлении инлайн и в направлении кросслайн, (3) комбинирование указанных значений инлайн и кросслайн для получения одного значения подобия для каждой ячейки, и (4) вывод на печать указанных значений подобия с получением распечатки значений подобия указанных ячеек по меньшей мере для одного из указанных горизонтальных срезов.

28. Сейсмическая карта по п. 27, отличающаяся тем, что до проведения операции (4) указанные значения подобия заносят в память в цифровом виде.

29. Сейсмическая карта по п. 27, отличающаяся тем, что операция (2) предусматривает проведение следующих операций: вычисление нулевой средней задержанной кросс-корреляции в указанном направлении инлайн и вычисление нулевой средней задержанной кросс-корреляции в указанном направлении кросслайн.

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

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