EA014282B1 - Геофизический способ разведки - Google Patents

Геофизический способ разведки Download PDF

Info

Publication number
EA014282B1
EA014282B1 EA200900300A EA200900300A EA014282B1 EA 014282 B1 EA014282 B1 EA 014282B1 EA 200900300 A EA200900300 A EA 200900300A EA 200900300 A EA200900300 A EA 200900300A EA 014282 B1 EA014282 B1 EA 014282B1
Authority
EA
Eurasian Patent Office
Prior art keywords
fourier transform
frequency
aliasing
wave numbers
seismic data
Prior art date
Application number
EA200900300A
Other languages
English (en)
Other versions
EA200900300A1 (ru
Inventor
Мичел Альберт Шоневайл
Original Assignee
Пгс Геофизикал Ас
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Пгс Геофизикал Ас filed Critical Пгс Геофизикал Ас
Publication of EA200900300A1 publication Critical patent/EA200900300A1/ru
Publication of EA014282B1 publication Critical patent/EA014282B1/ru

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Complex Calculations (AREA)

Abstract

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

Description

Область техники, к которой относится изобретение
Настоящее изобретение в целом относится к области геофизических способов разведки. Изобретение в особенности относится к области интерполяции или регуляризации трасс сейсмических данных.
Предшествующий уровень техники
В нефте- и газодобыче геофизические способы разведки обычно используются как одно из средств поиска и оценки подземных формаций. Геофизические способы разведки поставляют данные о геологическом строении земли, которые полезны для поиска и извлечения ценных минеральных ресурсов, в частности, для поиска залежей таких углеводородов, как нефть и природный газ. Хорошо известен такой геофизический способ разведки, как сейсморазведка. При наземной сейсморазведке сейсмический сигнал генерируется на поверхности или вблизи поверхности земли и затем распространяется вниз в геологическую среду. При морской сейсморазведке сейсмический сигнал также может распространяться вниз сквозь слой воды, покрывающий геологическую среду. Для генерации сейсмического сигнала применяются источники сейсмических волн; последние, распространяясь в землю, по меньшей мере, частично отражаются сейсмическими отражателями геологической среды. Как правило, такими сейсмическими отражателями являются границы подземных формаций, имеющих разные упругие свойства, в частности, разные скорости распространения звуковой волны и разные плотности породы, что ведет к изменениям акустического импеданса на их границах. Отраженные сейсмические волны улавливаются и записываются сейсмодатчиками (называемыми также сейсмоприемниками), расположенными на поверхности или вблизи поверхности земли, в слое воды, покрывающей геологическую среду, или на известной глубине в скважине.
Сейсмические источники, пригодные для генерирования сейсмического сигнала в наземной сейсморазведке, могут включать заряды взрывчатых веществ или вибраторы. В морской сейсморазведке обычно используются сейсмические источники, буксируемые судном и периодически активируемые для генерации акустических волн. Сейсмические источники для генерации акустических волн могут быть разных типов: это может быть небольшой заряд взрывчатого вещества, электрическая искра или дуга, морской вибратор и, чаще всего, пушка. В качестве таких сейсмических источников применяются водяная пушка, паровая пушка и, чаще всего, пневматическая пушка. Как правило, морской источник сейсмических колебаний включает не одиночный излучатель, а группу пространственно распределенных излучателей. Это в особенности справедливо для пневматических пушек, представляющих в настоящее время наиболее употребительный тип морских сейсмических источников.
К подходящим типам сейсмодатчиков обычно относят датчики скорости частиц, особенно в наземной сейсморазведке, и датчики давления воды, особенно в морской сейсморазведке. Иногда вместо или помимо датчиков скорости частиц используются датчики смещения частиц, датчики ускорения частиц или датчики градиента давления. Датчики скорости частиц и датчики давления воды в сейсморазведке обычно называют, соответственно, геофонами и гидрофонами. Сейсмодатчики могут развертываться по отдельности, но чаще они развертываются группами. Кроме того, в морской сейсморазведке датчики давления и датчики скорости частиц могут развертываться вместе; их располагают парами или парными группами.
При типовой организации морской сейсморазведки судно сейсморазведки идет по поверхности воды обычно со скоростью около 5 узлов и везет оборудование сбора сейсмических данных, в частности оборудование навигационного управления, управления сейсмическим источником, управления сейсмическими датчиками и записывающее оборудование. Оборудование управления сейсмическим источником в определенные моменты времени активирует сейсмический источник, буксируемый в подводном положении судном сейсморазведки. Судно сейсморазведки, буксирующее сейсмический источник, или другое судно тянет за собой находящиеся в подводном положении сейсмические косы, называемые также сейсмическими кабелями, которые представляют собой удлиненные напоминающие кабель конструкции. Как правило, судно сейсморазведки буксирует несколько сейсмических кос. Эти сейсмические косы содержат датчики для регистрации отраженных волн, вызванных сейсмическим источником и отраженных от отражающих границ. Обычно сейсмические косы содержат датчики давления, например гидрофоны, но предложены и сейсмические косы, содержащие помимо гидрофонов датчики скорости частиц воды, например геофоны, или датчики ускорения частиц, например акселерометры. Датчики давления и датчики движения частиц могут быть развернуты в непосредственной близости друг от друга и располагаться вдоль сейсмического кабеля парами или парными группами.
Итоговые сейсмические данные, собранные в процессе сейсморазведки, обрабатываются для выделения информации о геологической структуре и свойствах погребенных формаций в зоне проведения сейсморазведки. Эти обработанные сейсмические данные обрабатываются для визуального представления и анализа на потенциальное содержание углеводородов в составе погребенных формаций. Целью обработки сейсмических данных является извлечение из этих сейсмических данных возможно боле полной информации о погребенных формациях, чтобы дать адекватное изображение геологической среды. Для обнаружения мест геологической среды, где есть вероятность обнаружить залежи нефти, на сбор, обработку и интерпретацию сейсмических данных тратятся большие деньги. Результатом процесса построения на базе записанных сейсмоданных отражающих границ, которые определяют положение в не
- 1 014282 драх земли пластов, представляющих интерес, является временной или глубинный геологический разрез. Это изображение структуры геологической среды строится для того, чтобы дать возможность интерпретатору выделить места, в которых максимальна вероятность обнаружить залежи нефти. Для того чтобы убедиться в наличии месторождения нефти, надо пробурить скважину. Бурение скважин для определения наличия или отсутствия залежей нефти - чрезвычайно дорогое и затратное по времени предприятие. Поэтому существует непреходящая потребность в улучшении обработки и визуального представления сейсмических данных с тем, чтобы дать такое изображение структуры геологической среды, которое увеличит способность интерпретатора - независимо от того, проводит ли интерпретацию компьютер или человек - оценивать вероятность наличия месторождения нефти в конкретном месте геологической среды.
При сборе сейсмических данных, безразлично, наземном или морском, могут возникнуть две проблемы. Данные могут быть собраны при недостатке дискретных отсчетов (с элайсингом) или может быть неравномерна (нерегулярна) дискретизация собранных данных. Физические и экономические ограничения при проведении сейсморазведки часто приводят к тому, что сейсмические данные собираются с недостатком дискретных отсчетов или с неравномерной дискретизацией.
Данные, собранные при недостатке дискретных отсчетов, обычно называются данными с элайсингом (аПакеб ба!а). Теория дискретизации данных рекомендует не вводить в данные длину волны короче удвоенного интервала дискретизации. В противном случае для выделения особенности, соответствующей введенной длине волны, не хватит разрешения, и, следовательно, эта особенность будет искажена элайсингом.
Таким образом, временная частота неоднозначности, при которой начинается элайсинг, или частота
Найквиста, равна 2Д! в области частот и волновых чисел (области ί-к). Пространственная частота нек = — однозначности, при которой начинается элайсинг, или волновое число Найквиста, равно в пространственных координатах. Здесь Δΐ - временной интервал дискретизации в миллисекундах, а Дх - расстояние между точками наблюдения в предпочтительных единицах, например в метрах. Таким образом, большой интервал дискретизации в области время-пространство (ΐ-χ) соответствует низкой частоте Найквиста и малому волновому числу Найквиста в соответствующей £-к области.
Регистрация сейсмических данных с интервалами дискретизации, намного превышающими идеально подходящие, может отрицательно повлиять на дальнейшую обработку данных. Однако сбор данных сейсморазведки с более мелким интервалом дискретизации значительно увеличивает стоимость сбора сейсмических данных, особенно в случае трехмерной сейсморазведки. Поэтому вместо уменьшения интервала можно аппроксимировать отсутствующие данные, исходя из собранных. Таким образом, от интерполяции или экстраполяции сейсмических данных с пространственным элайсингом следует ожидать определенного экономического эффекта.
Интерполяция трасс данных без элайсинга, собранных с равномерной дискретизацией, проста. Интерполяция может быть выполнена, например, сверткой с 8шс-фильтром в пространственной области или продолжением посредством дополнения нулями волнового числа Найквиста сигнала с ограниченной полосой в область Фурье. Однако эта упрощенная интерполяция трасс предполагает интерполяцию с ортогональными базисными функциями. Если интерполяции трасс выполняются на сетке с нерегулярной дискретизацией, в данных происходит просачивание (1еакаде) энергии сигнала во все остальные частоты. Это просачивание энергии вызывается нерегулярностями дискретизации и граничными эффектами.
Таким образом, данные с неравномерной дискретизацией должны быть регуляризованы к ортогональной (регулярной) базисной сетке. Регуляризация трехмерных сейсмических данных требует генерации сейсмических трасс в местах, где при реальном расположении источника и приемника в процессе сейсморазведки регистрации данных не проводилось. Иными словами, сейсмические трассы, построенные по данным, собранным на нерегулярной сетке, интерполируются или экстраполируются на регулярную сетку.
В частности, морские сейсмические данные обычно дискретизированы нерегулярно и по пространственным направлениям скудны; причин тому много, и среди них снос кабеля, огибание препятствий, редактирование плохих трасс, а также соображения экономии. Но некоторые процедуры обработки сейсмических данных, включая такие как трехмерное подавление кратных волн, связанных с поверхностью, и трехмерная миграция волнового уравнения, требуют данных с регулярной дискретизацией. Наилучший способ получить трехмерные данные с регулярной дискретизацией - это собирать больше данных с большей избыточностью в направлении поперек профиля и в более широком азимутальном диапазоне, но это дорого и этого трудно добиться при существующих технологиях сбора морских данных. Поэтому регуляризация данных становится важным инструментом обработки.
Таким образом, имеется потребность в способе интерполяции трасс сейсмических данных, собранных одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией. В частности, имеется потребность в способе интерполяции трасс, уменьшающем просачивание энергии вследствие неравномерной дискретизации сейсмических данных с элайсингом.
- 2 014282
В качестве аналога предложенного изобретения можно рассмотреть способ интерполяции волн с сильным элайсингом, раскрытый в патентном документе США № 4922465, 01.05.1990, в котором суммарные сигналы участка трасс сейсмических данных разделяют посредством разложения волнового поля на частотные компоненты сигнала. Сигналы первого вида содержат частотные компоненты, идентифицируемые как компоненты с пространственным элайсингом, а сигналы второго вида содержат частотные компоненты, идентифицируемые как компоненты без элайсинга. Такие безэлайсинговые сигналы интерполируют известными методами, например, методом мпс интерполяции. Элайсинговые сигналы интерполируют методом Фр-дшбеб интерполяции для получения такого же количества трасс, что и в случае безэлайсинговых сигналов. Интерполированные сигналы первого и второго вида суммируют для получения участка сейсмических трасс, по существу, свободного от пространственного элайсинга.
Сущность изобретения
Настоящее изобретение предлагает способ интерполяции трасс сейсмических данных, которые могут быть собраны одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией. Оценочный спектр частот-волновых чисел формируют посредством применения способа антиликидж-преобразования Фурье к частотным составляющим без элайсинга сейсмических данных, преобразованных во временной области, и применения способа антиэлайсингового антиликидж-преобразования Фурье к частотным составляющим с элайсингом сейсмических данных, преобразованных во временной области. Согласно способу антиэлайсингового антиликидж-преобразования Фурье приписывают веса частотным составляющим с элайсингом в области частоты-волновые числа с использованием абсолютного спектра частот-волновых чисел, экстраполированного с частот без элайсинга на частоты с элайсингом. К оценочному спектру частот-волновых чисел применяют обратное временное и пространственное преобразование Фурье для формирования интерполяции трасс сейсмических данных.
Перечень чертежей
Настоящее изобретение и его преимущества будут прояснены в нижеследующем подробном описании и прилагаемых графических материалах, на которых представлены фиг. 1 - блок-схема, иллюстрирующая первый вариант осуществления настоящего изобретения для интерполяции трасс сейсмических данных, собранных одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией;
фиг. 2 - блок-схема, иллюстрирующая начальную часть второго варианта осуществления настоящего изобретения для интерполяции трасс сейсмических данных, собранных одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией;
фиг. 3 - блок-схема, иллюстрирующая промежуточную часть второго варианта осуществления настоящего изобретения для интерполяции трасс сейсмических данных; начальная часть показана на фиг. 2;
фиг. 4 - блок-схема, иллюстрирующая финальную часть второго варианта осуществления настоящего изобретения для интерполяции трасс сейсмических данных; начальная часть показана на фиг. 2;
фиг. 5 - блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для обработки частот без элайсинга в неравномерно дискретизированных сейсмических данных с фиг. 2;
фиг. 6 - блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для обработки частот с элайсингом в неравномерно дискретизированных сейсмических данных с фиг. 3;
фиг. 7 - спектр частот-волновых чисел на низких частотах без элайсинга;
фиг. 8 - экстраполированный спектр частот-волновых чисел, использованный в качестве весовой функции для высоких частот с элайсингом;
фиг. 9Α-9Ό - пример применения способа согласно изобретению к синтетическим сейсмоданным;
фиг. 10Α-10Ό - пример применения стандартного способа АЬРТ к синтетическим сейсмоданным с фиг. 9Α-9Ό;
фиг. 11Α-11Ό - пример применения способа согласно изобретению к полевым сейсмоданным и фиг. 12Α-12Ό - пример применения стандартного способа ΑΒΡΤ к полевым сейсмоданным с фиг. 11А-1Ш.
Хотя изобретение описывается далее в связи с предпочтительными вариантами осуществления, следует понимать, что они не ограничивают объема настоящего изобретения. Напротив, изобретение охватывает все альтернативные варианты, модификации и эквиваленты, которые могут быть включены в объем настоящего изобретения в соответствии с прилагаемой формулой изобретения.
Сведения, подтверждающие возможность осуществления изобретения
Самый простой способ интерполяции трасс - это одномерная интерполяция по определенному пространственному направлению для любой временной выборки посредством преобразования Фурье в определенном направлении области волновых чисел (к) с дополнением достаточным количеством нулевых выборок и обратное преобразование Фурье. Эквивалентная интерполяция трасс может быть выполнена в области частоты - волновые числа (ί-к) для каждого частотного среза данных. Указанную процедуру реализует интерполятор трасс - детерминистская, не зависящая от данных ыпс-функция (иначе - кардинальная синусная функция). Этот процесс хорошо работает, если в сейсмических данных нет недостатка дискретных отсчетов (элайсинга) или неравномерной дискретизации (на нерегулярной сетке).
- 3 014282
Интерполяция трасс как регуляризация сейсмических данных, имеющая целью оценить сейсмотрассы на пространственно регулярной сетке, используя собранные с нерегулярной дискретизацией данные, представляет собой проблему интерполяции/экстраполяции. Критическую роль при оценке частотных составляющих в области частоты - волновые числа играет преобразование Фурье, а соответствующее обратное преобразование Фурье воссоздает сейсмические данные на желаемой регулярной сетке вновь в пространственно-временной области. Основная проблема реконструкции данных на основе преобразования Фурье заключается в том, что глобальные базисные функции, такие как бае-функции. на нерегулярной сетке неортогональны. Такая неортогональность глобальных базисных функций преобразования Фурье приводит к тому, что энергия с одной частотной составляющей просачивается на другие, это явление называется просачиванием спектральных составляющих.
Настоящее изобретение предлагает способ интерполяции трасс сейсмических данных, которые могут быть собраны одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией. Изобретение базируется на известном способе интерполяции отсутствующих трасс - антиликиджпреобразовании Фурье (ЛЬЕТ: Аиб-1еакаде Еоштет ШцъГогт - антиликидж-преобразование Фурье). Способ АЕЕТ был предложен для решения проблемы регуляризации сейсмических данных путем ослабления эффекта частотного просачивания в данных, дискретизированных на нерегулярной сетке. Способ АЬЕТ реортогонализует глобальные базисные функции Фурье на нерегулярно дискретизированной сетке, что дает хорошую оценку спектра сигнала на нерегулярной сетке.
Для стандартных способов суммирования в преобразовании Фурье последовательность определения каждого из коэффициентов Фурье не влияет на конечный результат. Но в способах АЕЕТ порядок имеет критическое значение, так как коэффициенты Фурье с большей величиной (энергией) дадут больший вклад в энергию просачивания, чем коэффициенты с меньшей величиной. Таким образом, для того чтобы уменьшить просачивание, коэффициенты Фурье оцениваются итерационно, при этом в каждой точке оценивается коэффициент с максимальной энергией. После каждой оценки рассчитанная частотная составляющая (коэффициент Фурье) обнуляется корректировкой входных данных. Математически это эквивалентно удалению этой частотной составляющей из исходных входных сейсмических данных.
Эти новые остаточные входные данные используются для определения следующего коэффициента Фурье, вновь с применением того же критерия максимума энергии. Описанная итерационная процедура повторяется до тех пор, пока не будут определены все коэффициенты Фурье, то есть до тех пор, пока все откорректированные входные значения не будут стремиться к нулю (практически, не окажутся ниже некоторого порогового значения). Вообще говоря, глобальные базисные функции Фурье ортогональны только на регулярной сетке, то есть для регулярно дискретизированных данных. Описанное выше вычитание работает как механизм ортогонализации базисных функций Фурье на нерегулярной сетке. Иными словами, базисные функции Фурье реортогонализуются. Этот способ представляет практическое решение проблемы минимизации просачивания энергии с одной частоты на другую.
Если количество и диапазон составляющих Фурье достаточны, то окончательные откорректированные входные данные на нерегулярной сетке будут после всех операций вычитания стремиться к нулю. В этом случае реконструированные по полученным коэффициентам Фурье данные будут соответствовать результатам оригинальных замеров; это одно из требований, предъявляемых к желательному способу интерполяции.
Проблемы при использовании стандартного способа АЬЕТ возникают в связи с обработкой данных с элайсингом (в присутствии шума). Составляющая с элайсингом может быть такой же или большей величины, чем составляющая без элайсинга, и может быть ошибочно изъята из последовательности вычислений по способу АЬЕТ. Способ согласно изобретению развивает способ АЬЕТ путем использования информации от низких временных частот без элайсинга, помогающей провести деэлайсинг высоких временных частот. Конкретнее, в способе согласно изобретению используются низкие частоты без элайсинга для построения весовой функции, определяющей, какие спектральные составляющие в процедуре АЬЕТ рассчитываются (и удаляются) первыми. Эта весовая функция конструируется посредством экстраполяции спектра частот - волновых чисел для частот без элайсинга на частоты с элайсингом.
Фиг. 1-6 представляют собой блок-схемы, иллюстрирующие варианты осуществления настоящего изобретения для интерполяции трасс. Фиг. 1 и 2-4 показывают, соответственно, два варианта осуществления настоящего изобретения, в то время как на фиг. 5 и 6 представлены дальнейшие подробности способа, показанного на фиг. 2-4. Фиг. 7-12 иллюстрируют некоторые из методик, описанных в блок-схемах фиг. 1-6. Фиг. 6 и 7 изображают спектры, использованные при определении весов. Фиг. 9А-9Э и 10А-10Э представляют результаты интерполяции трасс на примере синтетических сейсмоданных в случаях применения способа согласно изобретению и стандартного способа АЬЕТ соответственно. Фиг. 11А-11Б и 12А-12Б представляют результаты интерполяции трасс на примере полевых сейсмоданных в случаях применения способа согласно изобретению и стандартного способа АЬЕТ соответственно.
Фиг. 1 - блок-схема, иллюстрирующая первый вариант осуществления настоящего изобретения для интерполяции трасс по входному набору сейсмических данных, собранных одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией.
В блоке 11 формируется оценочный спектр частот - волновых чисел. Первый, стандартный способ
- 4 014282 антиликидж-преобразования Фурье применяется к неэлайсинговым частотным составляющим сейсмических данных, преобразованных во временной области, а второй, нестандартный способ антиликиджпреобразования Фурье применяется к элайсинговым частотным составляющим сейсмических данных, преобразованных во временной области. Этот второй способ антиликидж-преобразования Фурье переносит абсолютный спектр в области частоты - волновые числа, экстраполированный с частот (обычно низких) без элайсинга на частоты (обычно высокие) с элайсингом, чтобы приписать веса частотным составляющим с элайсингом.
В блоке 12 обратное временное и пространственное преобразование Фурье применяется к оценочному спектру частот - волновых чисел из блока 11, при этом формируется оценочный набор данных. Оценочный спектр частот - волновых чисел получен обратным преобразованием из области частот - волновых чисел (области Г-к) вновь в область время - пространство (1-х), что дает в результате желаемую интерполяцию трасс входных сейсмических данных.
Фиг. 2 - блок-схема, иллюстрирующая начальную часть второго варианта осуществления настоящего изобретения для интерполяции трасс сейсмических данных, собранных одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией. Фиг. 2 развивает первый вариант осуществления настоящего изобретения, рассмотренный выше со ссылкой на фиг. 1.
В блоке 21 происходит получение входных сейсмических данных. Предполагается, что эти сейсмические данные получают в пространственно-временной области (1-х). Эти сейсмические данные могут быть одновременно и с недостатком дискретных отсчетов, и с неравномерной дискретизацией, поскольку способ согласно изобретению адаптирован к использованию при обоих этих условиях.
В блоке 22 временное преобразование Фурье применяется к входным сейсмическим данным, полученным в блоке 21. Для ускорения расчетов в качестве преобразования Фурье предпочтительно применять быстрое преобразование Фурье (ЕЕТ: Еа§1 Еоштет ТгапГогт - быстрое преобразование Фурье). Входные данные преобразуются из пространственно-временной области (1-х) в частотнопространственную (Г-х).
В блоке 23 определяется, какие частотные составляющие преобразованных в блоке 22 входных сейсмических данных являются составляющими без элайсинга, а какие - с элайсингом. Как правило, низкие частоты бывают без элайсинга, а высокие частоты - с элайсингом.
В блоке 24 выбирается первое число итераций, обозначаемое Ν1, для первого способа ЛЬЕТ, который должен быть применен к частотам без элайсинга, определенным в блоке 23. Число Ν1 для первого способа ЛЬЕТ определяется эмпирически.
В альтернативном варианте осуществления настоящего изобретения для первого способа ЛЬЕТ вместо первого числа итераций, Ν1, выбирается первый порог ε1. В этом случае первый способ ЛЬЕТ предпочтительно применяется с итерациями к каждому коэффициенту Фурье до тех пор, пока этот коэффициент не уменьшится до значения ниже первого порога ε1, а не к числу последовательных итераций. Способ согласно изобретению проиллюстрирован вариантом осуществления, использующим число итераций, исключительно в иллюстративных целях, и этот выбор не должен рассматриваться как ограничивающий объем изобретения.
В блоке 25 из тех частотных составляющих, которые, как определено в блоке 23, должны быть без элайсинга, выбирается частотная составляющая без элайсинга.
В блоке 26 осуществляется переход к блоку 51 на фиг. 5 для обработки частотной составляющей без элайсинга, выбранной в блоке 25. Там, в процессе обработки, показанном на фиг. 5, к частотам без элайсинга применяется первый (стандартный) способ ЛЬЕТ.
В блоке 27 получают из блока 59 с фиг. 5 оценочный спектр частот - волновых чисел для выбранной частотной составляющей без элайсинга.
В блоке 28 определяется, остались ли еще частотные составляющие без элайсинга, которые должны быть выбраны. Если такие частотные составляющие без элайсинга остались, процедура возвращается к блоку 25. Если же таких частотных составляющих без элайсинга не осталось, процедура продолжается переходом к следующему блоку 29.
В блоке 29 все оценочные спектры частот - волновых чисел для выбранных частотных составляющих без элайсинга из блока 27 объединяются в один безэлайсинговый оценочный спектр частот - волновых чисел.
В блоке 30 заканчивается процесс, изображенный на фиг. 2; продолжение - в блоке 31 на фиг. 3, где обрабатывается безэлайсинговый оценочный спектр частот - волновых чисел из блока 29.
Фиг. 3 - блок-схема, иллюстрирующая промежуточную часть второго варианта осуществления настоящего изобретения для процесса интерполяции трасс сейсмических данных, начало которого показано на фиг. 2.
В блоке 31 продолжается процедура блока 29 фиг. 2 с безэлайсинговым оценочным спектром частот волновых чисел.
В блоке 32 берутся абсолютные значения составляющих безэлайсингового оценочного спектра частот - волновых чисел из блока 31 и формируется абсолютный частотный спектр в области частот - вол
- 5 014282 новых чисел (области Г-к) для частот без элайсинга.
В блоке 33 этот абсолютный спектр частот - волновых чисел из блока 32 экстраполируется на частоты с элайсингом, в результате чего формируется экстраполированный абсолютный спектр частот волновых чисел в области частот - волновых чисел (Г-к) для частот с элайсингом. Как правило, к частотам с элайсингом относятся высокие частоты и волновые числа. Этот экстраполированный спектр содержит полезную информацию, которую несут низкие частоты без элайсинга, и эта информация улучшает выбор составляющих Фурье без элайсинга на высоких частотах с элайсингом.
На фиг. 7-8 даны примеры спектров, которые рассчитываются в блоках 32 и 33. На фиг. 7 показан спектр частот - волновых чисел, рассчитанный в блоке 32 при низких частотах, на которых нет элайсинга. На фиг. 8 показан экстраполированный спектр частот - волновых чисел, рассчитанный в блоке 33. Этот экстраполированный спектр используется в качестве весовой функции для высоких частот, часть которых - частоты с элайсингом.
В принципе, экстраполяция производится с низких частот на высокие и, следовательно, на большие полосы или волновые числа. Практически, низкочастотный спектр интерполируется по значениям и частот, и волновых чисел. Эта интерполяция может включать осреднение или сглаживание. На фиг. 8 показаны результаты интерполяции 2:1.
В блоке 34 выбирается второе число итераций, обозначаемое Ν2, для второго способа АЬРТ, который должен быть применен к частотам с элайсингом, определенным в блоке 23 на фиг. 2. Число Ν2 для второго способа АЬРТ определяется эмпирически. В одном из вариантов осуществления настоящего изобретения Ν1 = Ν2, хотя обычно Ν1 > Ν2. Однако в объем настоящего изобретения входит определение своего числа итераций Ν для каждой частотной составляющей, то есть частотная зависимость числа итераций.
В альтернативном варианте осуществления настоящего изобретения для второго способа АЬРТ вместо второго числа итераций, Ν2, выбирается второй порог ε2. В этом случае второй способ АЬРТ предпочтительно применяется с итерациями к каждому коэффициенту Фурье до тех пор, пока этот коэффициент не уменьшится до значения ниже второго порога ε2, а не к числу последовательных итераций. Способ согласно изобретению проиллюстрирован вариантом осуществления, использующим число итераций, исключительно в иллюстративных целях, и этот выбор не должен рассматриваться как ограничивающий объем изобретения.
В блоке 35 из тех частотных составляющих, которые, как определено в блоке 23 фиг. 2, должны быть с элайсингом, выбирается частотная составляющая с элайсингом.
В блоке 36 осуществляется переход к блоку 61 фиг. 6 для обработки частотной составляющей с элайсингом, выбранной в блоке 35. Там, в процессе обработки, показанном на фиг. 6, к частотам с элайсингом применяется второй (нестандартный) способ АЬРТ согласно изобретению.
В блоке 37 получают из блока 71 фиг. 6 оценочный спектр частот - волновых чисел для выбранной частотной составляющей с элайсингом.
В блоке 38 определяется, остались ли еще частотные составляющие с элайсингом, которые должны быть выбраны. Если такие частоты с элайсингом остались, процедура возвращается к блоку 35. Если же таких частот с элайсингом не осталось, процедура продолжается переходом к следующему блоку 39.
В блоке 39 заканчивается процесс, изображенный на фиг. 3; продолжение - в блоке 41 на фиг. 4, где обрабатываются безэлайсинговый оценочный спектр частот - волновых чисел из блока 31 и все оценочные спектры частот - волновых чисел для выбранных частотных составляющих с элайсингом, полученные из блока 37.
Фиг. 4 - блок-схема, иллюстрирующая финальную часть второго варианта осуществления настоящего изобретения для процесса интерполяции трасс сейсмических данных, начало которого показано на фиг. 2, а продолжение - на фиг. 3.
В блоке 41 продолжается процедура блока 39 фиг. 3: обрабатываются безэлайсинговый оценочный спектр частот - волновых чисел и все оценочные спектры частот - волновых чисел для выбранных частотных составляющих с элайсингом.
В блоке 42 все оценочные спектры частот - волновых чисел для выбранных частотных составляющих с элайсингом из блока 41 объединяются в один элайсинговый оценочный спектр частот - волновых чисел.
В блоке 43 безэлайсинговый и элайсинговый оценочные спектры частот - волновых чисел из блоков 41 и 42, соответственно, объединяются в общий оценочный спектр частот - волновых чисел.
В блоке 44 к общему оценочному спектру частот - волновых чисел из блока 43 применяется обратное пространственное преобразование Фурье, в результате чего формируется общий оценочный частотно-пространственный спектр. Обратное пространственное преобразование Фурье конструируется для преобразования трасс в желаемые положения, включая преобразование в положения отсутствующих трасс или в положения на регулярной (ортогональной) сетке. Для ускорения расчетов в качестве обратного преобразования Фурье предпочтительно применять обратное пространственное дискретное преобразование Фурье (ЭРТ: Э|5сгс1с Роштет ТгапГогт - дискретное преобразование Фурье) или неравномер
- 6 014282 ное быстрое преобразование Фурье (ΝΡΕΤ: Νοηιιηίίοηη Ра§1 Роштет ТтаикГотш - неравномерное быстрое преобразование Фурье). Общий оценочный спектр частот - волновых чисел обратно преобразуется из области частот - волновых чисел (области Г-к) в общий оценочный частотно-пространственный спектр в частотно-пространственной (Г-х) области.
В блоке 45 обратное временное преобразование Фурье применяется к преобразованному общему оценочному частотно-пространственному спектру из блока 44, в результате чего формируется общий оценочный набор данных. Для ускорения расчетов в качестве обратного преобразования Фурье предпочтительно применять обратное временное быстрое преобразование Фурье (РРТ). Общий оценочный частотно-пространственный спектр далее обратно преобразуется из частотно-пространственной (Г-х) области в общий оценочный набор данных в области время - пространство (1-х).
В альтернативном варианте обратное преобразование Фурье, применяемое в блоках 44 и 45, представляет собой двумерное (временное и пространственное) быстрое преобразование Фурье (РРТ). В обоих случаях в конечном результате оценочный спектр частот - волновых чисел обратно преобразуется из области частот - волновых чисел (области Г-к) вновь в область время - пространство (ΐ-х), в результате чего формируется желаемая интерполяция трасс входных сейсмических данных. Эта интерполяция трасс входных сейсмических данных может решить много задач, включая восполнение отсутствующих трасс и регуляризацию дискретизации трасс.
Для частот без элайсинга, определенных в блоке 23 фиг. 2, применяется стандартный способ ЛЬРТ. На фиг. 5 показана блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для обработки частот без элайсинга в неравномерно дискретизированных сейсмических данных с фиг. 2.
В блоке 51 получают выбранную частотную составляющую без элайсинга из блока 26 фиг. 2.
В блоке 52 формируется оценочный спектр частот - волновых чисел в области частот - волновых чисел (Г-к) для выбранной частотной составляющей без элайсинга, полученной в блоке 51. Этот оценочный спектр частот - волновых чисел исходно приравнивается нулю, а далее наращивается добавлением выбранных составляющих волновых чисел в блоке 55, описанном ниже.
В блоке 53 пространственное преобразование Фурье применяется к выбранной в блоке 51 частотной составляющей без элайсинга. Для ускорения расчетов в качестве преобразования Фурье предпочтительно применять дискретное преобразование Фурье (ЭРТ) или неравномерное быстрое преобразование Фурье (№РТ). Выбранная частотная составляющая преобразуется из частотно-пространственной области (Г-х) в область частот - волновых чисел (область Г-к).
В блоке 54 выбирается наибольшая составляющая волновых чисел преобразованной частотной составляющей из блока 53. Наибольшая составляющая волновых чисел - это составляющая Фурье с наибольшей величиной (энергией), возникшая из рассчитанного в блоке 53 пространственного преобразования Фурье.
В блоке 55 наибольшая составляющая волновых чисел, выбранная в блоке 54, добавляется к оценочному спектру частот - волновых чисел, который был сформирован с заданием начальных условий в блоке 52 для выбранной частотной составляющей без элайсинга.
В блоке 56 обратное пространственное преобразование Фурье применяется к выбранной в блоке 54 наибольшей составляющей Фурье. Для ускорения расчетов в качестве обратного преобразования Фурье предпочтительно применять обратное дискретное преобразование Фурье (ЭРТ) или обратное неравномерное быстрое преобразование Фурье (ЫЕРТ). Наибольшая составляющая Фурье обратно преобразуется из области (Г-к) частот - волновых чисел вновь в частотно-пространственную область (Г-х).
В блоке 57 обратно-преобразованная наибольшая составляющая, рассчитанная в блоке 56, вычитается из частотной составляющей без элайсинга, полученной в блоке 51. Эти вычитания итерационно формируют откорректированную частотную составляющую без элайсинга.
В блоке 58 определяется, достигнуто ли число Ν1 итераций в блоках 53-57 для частотной составляющей без элайсинга, полученной в блоке 51. Если число Ν1 итераций не достигнуто, процедура возвращается к блоку 53 для следующей итерации; если же число Ν1 итераций достигнуто, процедура Фиг. 5 продолжается переходом к блоку 59.
В альтернативном варианте осуществления настоящего изобретения, описанном в блоке 24 фиг. 2, определяется, стала ли уже откорректированная частотная составляющая без элайсинга из блока 57 меньше первого порогового значения ε1. Если она не стала меньше порогового значения ε1, процедура возвращается к блоку 53 для следующей итерации; если же она стала меньше порогового значения ε£, процедура фиг. 5 продолжается переходом к блоку 59.
В блоке 59 процедура фиг. 5 заканчивается и происходит возврат к блоку 27 фиг. 2, куда передается оценочный спектр частот - волновых чисел для выбранной частотной составляющей без элайсинга, итерационно сформированной в блоке 55.
Для частот с элайсингом, определенных в блоке 27 фиг. 2, применяется нестандартный способ ЛЬРТ согласно настоящему изобретению. На фиг. 6 показана блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для обработки частот с элайсингом в неравномерно дискретизированных сейсмических данных с фиг. 3.
- 7 014282
В блоке 61 получают выбранную частотную составляющую с элайсингом из блока 36 фиг. 3.
В блоке 62 формируется оценочный спектр частот - волновых чисел в области частот - волновых чисел (ί-к) для выбранной частотной составляющей с элайсингом, полученной в блоке 51. Этот оценочный спектр частот - волновых чисел исходно приравнивается нулю, а далее наращивается добавлением выбранных составляющих волновых чисел в блоке 67, описанном ниже.
В блоке 63 пространственное преобразование Фурье применяется к выбранной в блоке 61 частотной составляющей с элайсингом. Для ускорения расчетов в качестве преобразования Фурье предпочтительно применять дискретное преобразование Фурье (ΌΡΤ) или неравномерное быстрое преобразование Фурье (ΝΡΡΤ). Выбранная частотная составляющая преобразуется из частотно-пространственной области (ί-χ) в область частот - волновых чисел (область ί-к).
В блоке 64 экстраполированный абсолютный спектр частот - волновых чисел из блока 33 фиг. 3 применяется к преобразованной частотной составляющей, рассчитанной в блоке 63, чтобы определить весовые коэффициенты составляющих волновых чисел преобразованной частотной составляющей. Пример такого экстраполированного спектра частот - волновых чисел показан на фиг. 8.
В блоке 65 выбирается наибольшая составляющая волновых чисел преобразованной и взвешенной частотной составляющей из блока 64. Наибольшая составляющая волновых чисел - это составляющая Фурье с наибольшей величиной (энергией), возникшая из рассчитанного в блоке 53 пространственного преобразования Фурье.
В блоке 66 получают невзвешенную составляющую волновых чисел преобразованной частотной составляющей из блока 63, соответствующую наибольшей взвешенной составляющей волновых чисел, определенной в блоке 65.
В блоке 67 полученную в блоке 66 невзвешенную составляющую волновых чисел, соответствующую наибольшей составляющей волновых чисел, добавляют к оценочному спектру частот - волновых чисел, сформированному с заданием начальных условий в блоке 62 для выбранной частотной составляющей с элайсингом.
В блоке 68 обратное пространственное преобразование Фурье применяется к определенной в блоке 66 наибольшей невзвешенной составляющей. Для ускорения расчетов в качестве обратного преобразования Фурье предпочтительно применять обратное дискретное преобразование Фурье (ΌΡΤ) или обратное неравномерное быстрое преобразование Фурье (ΝΡΡΤ). Наибольшая невзвешенная составляющая обратно-преобразуется из области (ί-к) частот - волновых чисел вновь в частотно-пространственную область (ί-χ).
В блоке 69 обратно-преобразованная наибольшая невзвешенная составляющая, рассчитанная в блоке 68, вычитается из частотной составляющей с элайсингом, полученной в блоке 61. Эти вычитания итерационно формируют откорректированную частотную составляющую с элайсингом.
В блоке 70 определяется, достигнуто ли число Ν2 итераций в блоках 63-69 для частотной составляющей с элайсингом, полученной в блоке 51. Если число Ν2 итераций не достигнуто, процедура возвращается к блоку 63 для следующей итерации; если же число Ν2 итераций достигнуто, процедура продолжается переходом к блоку 71.
В альтернативном варианте осуществления настоящего изобретения, описанном в блоке 34 фиг. 3, определяется, стала ли уже откорректированная частотная составляющая без элайсинга из блока 69 меньше второго порогового значения ε2. Если она не стала меньше порогового значения ε2, процедура возвращается к блоку 63 для следующей итерации; если же она стала меньше порогового значения ε2, процедура фиг. 6 продолжается переходом к блоку 71.
В блоке 71 процедура фиг. 6 заканчивается и происходит возврат к блоку 37 фиг. 3, куда передается оценочный спектр частот - волновых чисел для выбранной частотной составляющей с элайсингом, итерационно сформированной в блоке 57.
На фиг. 9Ά-12Ό представлены для сравнения результаты интерполяции трасс в случаях применения способа согласно изобретению и стандартного способа ΆΕΡΤ. Фиг. 9Ά-9Ό и 10Ά-10Ό представляют результаты интерполяции трасс синтетических сейсмоданных, в то время как фиг. 11Ά-11Ό и 12Ά-12Ό представляют результаты интерполяции трасс полевых сейсмоданных.
На фиг. 9Ά-9Ό показан пример применения способа согласно изобретению к синтетическим сейсмоданным. На фиг. 9 А представлены оригинальные синтетические сейсмоданные. На фиг. 9В представлены входные данные с удалением трасс для моделирования данных с элайсингом. На фиг. 9С представлены данные, полученные интерполяцией с применением способа согласно изобретению. Фиг. 9Ό демонстрирует расхождения между оригинальными данными фиг. 9А и данными после интерполяции фиг. 9С; малые величины этих расхождений свидетельствуют о хорошем согласии.
На фиг. 10Ά-10Ό показан для сравнения пример применения стандартного способа ΑΕΡΤ к синтетическим сейсмоданным с фиг. 9Ά-9Ό. На фиг. 10А представлены оригинальные синтетические сейсмоданные, как и на фиг. 9А. На фиг. 10В представлены входные данные с удалением трасс для моделирования данных с элайсингом, как и на фиг. 9В. На фиг. 10С представлены данные, полученные интерполяцией с применением стандартного способа ΑΕΡΤ. Фиг. 9Ό демонстрирует расхождения между ориги
- 8 014282 нальными данными фиг. 10А и данными после интерполяции фиг. 10С; более значительные величины этих расхождений свидетельствуют о не таком хорошем согласии, как на вышеприведенной фиг. 9Ό для способа согласно изобретению.
На Фиг. 11Α-11Ό показан пример применения способа согласно изобретению к полевым сейсмоданным. На Фиг. 11А представлены оригинальные полевые сейсмоданные. На фиг. 11В представлены входные данные с удалением трасс для моделирования данных с элайсингом. На фиг. 11С представлены данные, полученные интерполяцией с применением способа согласно изобретению. Фиг. 11Ό демонстрирует расхождения между оригинальными данными фиг. 11А и данными после интерполяции фиг. 11 С.
На фиг. 12Α-12Ό показан для сравнения пример применения стандартного способа АЬЕТ к полевым сейсмоданным с фиг. 11Α-11Ό. На фиг. 12А представлены оригинальные полевые сейсмоданные, как и на фиг. 11А. На фиг. 12В представлены входные данные с удалением трасс для моделирования данных с элайсингом, как и на фиг. 11В. На фиг. 12С представлены данные, полученные интерполяцией с применением стандартного способа АЬЕТ. Фиг. 12Ό демонстрирует расхождения между оригинальными данными фиг. 12А и данными после интерполяции фиг. 12С.
Нормализованное среднеквадратичное отклонение (ΝΚΜ8: поппай/еб τοοΐ теап 5с.|иаге - нормализованное среднеквадратичное отклонение) данных после интерполяции, представленных на фиг. 11С (интерполяция согласно изобретению), от оригинальных данных, представленных на фиг. 11 А, составляет 70%, в то время как ΝΚΜ8 данных после интерполяции, представленных на фиг. 12С (интерполяция по стандартному способу АЬЕТ), от оригинальных данных, представленных на фиг. 12А, составляет 84%. Таким образом, превосходство способа согласно изобретению по сравнению со стандартным способом АЬЕТ показано и на примере полевых данных.
Применение способа антиэлайсингового антиликидж-преобразования Фурье согласно настоящему изобретению дает лучшую интерполяцию трасс сейсмоданных с элайсингом, чем применение только стандартного способа АЬЕТ. Стоимость дополнительных расчетов оценивается как весьма ограниченная. Раскрытый в настоящем описании способ согласно изобретению легко может быть распространен на многомерные варианты осуществления, включая трехмерный с двумя пространственными измерениями плюс время, четырехмерный с тремя пространственными измерениями плюс время и пятимерный с четырьмя пространственными измерениями плюс время. Эти пространственные измерения могут включать подмножество координат источника х, у и приемника х, у или, что равнозначно, подмножество координат средней точки вдоль профиля, средней точки поперек профиля, удаления и азимута. В альтернативных вариантах осуществления настоящего изобретения временная координата может быть частотной или глубинной.
Способ согласно изобретению может быть также распространен на способ восстановления изображений по одномерному и многомерному алгоритмам проекций на выпуклое множество (РОС8: рто_)ес1юп-оп1о-сопуех-5е15 - проекции на выпуклое множество). Возможны также и другие вариации, включая сглаживание спектра частот - волновых чисел для получения улучшенных значений весов и применение других схем взвешивания. Другие вариации включают - но не ограничительно - использование более высоких частот, или информации других выборок, или отличающихся способов оценки спектра частот волновых чисел при низких частотах, например преобразования Фурье по способу наименьших квадратов.
Следует учитывать, что все предшествующее было лишь подробным описанием специальных вариантов осуществления настоящего изобретения и что возможны многочисленные изменения, модификации и альтернативы раскрытых вариантов осуществления, не выходящие за объем настоящего изобретения. Таким образом, предшествующее описание не ограничивает объема настоящего изобретения. Напротив, объем настоящего изобретения определяется только пунктами прилагаемой формулы изобретения и эквивалентами.

Claims (21)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Геофизический способ разведки путем получения сейсмических данных, в котором формируют оценочный спектр частот - волновых чисел посредством применения способа антиликиджпреобразования Фурье к частотным составляющим без элайсинга сейсмических данных, преобразованных во временной области, и применения способа антиэлайсингового антиликидж-преобразования Фурье к частотным составляющим с элайсингом сейсмических данных, преобразованных во временной области, причем согласно способу антиэлайсингового антиликидж-преобразования Фурье приписывают веса частотным составляющим с элайсингом в области частоты-волновые числа с использованием абсолютного спектра частот - волновых чисел, экстраполированного с частот без элайсинга на частоты с элайсингом; и применяют обратное временное и пространственное преобразование Фурье к оценочному спектру частот - волновых чисел для формирования интерполяции трасс сейсмических данных.
  2. 2. Способ по п.1, в котором частоты без элайсинга представляют собой, по существу, низкие частоты сейсмических данных, а частоты с элайсингом представляют собой, по существу, высокие частоты сейсмических данных.
    - 9 014282
  3. 3. Способ по п.1, в котором сначала получают входные сейсмические данные в пространственновременной области; применяют временное преобразование Фурье к входным сейсмическим данным для формирования преобразованных сейсмических данных; определяют, какие частотные составляющие преобразованных сейсмических данных соответствуют частотам без элайсинга и какие частотные составляющие преобразованных сейсмических данных соответствуют частотам с элайсингом; и применяют способ антиликидж-преобразования Фурье к каждой частотной составляющей без элайсинга преобразованных сейсмических данных.
  4. 4. Способ по п.3, в котором временное преобразование Фурье представляет собой быстрое преобразование Фурье.
  5. 5. Способ по п.3, в котором определяют число итераций Ν1 для способа антиликиджпреобразования Фурье.
  6. 6. Способ по п.5, в котором определение числа итераций Ν1 включает определение порога ε1, для способа антиликидж-преобразования Фурье.
  7. 7. Способ по п.6, в котором осуществляют начальное приравнивание нулю оценочного спектра частот - волновых чисел для частотной составляющей без элайсинга; применяют пространственное преобразование Фурье к выбранной частотной составляющей без элайсинга; и выполняют следующие итерации для Ν1: выбирают наибольшую составляющую волновых чисел преобразованной частотной составляющей без элайсинга, добавляют наибольшую составляющую волновых чисел к оценочному спектру частот волновых чисел для частотной составляющей без элайсинга, применяют обратное пространственное преобразование Фурье к выбранной наибольшей составляющей волновых чисел, и вычитают обратнопреобразованную наибольшую составляющую из выбранной частотной составляющей без элайсинга для формирования откорректированной частотной составляющей.
  8. 8. Способ по п.7, в котором пространственное преобразование Фурье представляет собой дискретное преобразование Фурье.
  9. 9. Способ по п.7, в котором обратное пространственное преобразование Фурье представляет собой обратное дискретное преобразование Фурье.
  10. 10. Способ по п.7, в котором пространственное преобразование Фурье представляет собой неравномерное быстрое преобразование Фурье.
  11. 11. Способ по п.7, в котором обратное пространственное преобразование Фурье представляет собой обратное неравномерное быстрое преобразование Фурье.
  12. 12. Способ по п.7, в котором объединяют оценочные спектры частот - волновых чисел для частотных составляющих без элайсинга для формирования безэлайсингового оценочного спектра частотволновых чисел; определяют абсолютные значения составляющих безэлайсингового оценочного спектра частот - волновых чисел для формирования абсолютного спектра частот-волновых чисел; осуществляют экстраполяцию абсолютного спектра частот-волновых чисел на частоты без элайсинга для формирования экстраполированного абсолютного спектра частот-волновых чисел; и применяют способ антиэлайсингового антиликидж-преобразования Фурье к каждой из частотных составляющих с элайсингом преобразованных сейсмических данных.
  13. 13. Способ по п.12, в котором сглаживают абсолютный спектр частот - волновых чисел перед экстраполированием абсолютного спектра частот-волновых чисел на частоты без элайсинга.
  14. 14. Способ по п.12, в котором определяют число итераций Ν2 для способа антиэлайсингового антиликидж-преобразования Фурье.
  15. 15. Способ по п.14, в котором определение числа итераций Ν2 включает определение порога ε2 для способа антиэлайсингового антиликидж-преобразования Фурье.
  16. 16. Способ по п.7, в котором осуществляют начальное приравнивание нулю оценочного спектра частот - волновых чисел для частотной составляющей с элайсингом; применяют пространственное преобразование Фурье к выбранной частотной составляющей с элайсингом и выполняют следующие итерации для Ν2: приписывают веса преобразованной частотной составляющей с элайсингом с использованием экстраполированного абсолютного спектра частот-волновых чисел для формирования взвешенного элайсингового спектра частот - волновых чисел, выбирают наибольшую составляющую волновых чисел взвешенного элайсингового спектра частот-волновых чисел, получают невзвешенную составляющую волновых чисел, соответствующую наибольшей взвешенной составляющей волновых чисел, добавляют соответствующую невзвешенную составляющую волновых чисел к оценочному спектру частот - волновых чисел для частотной составляющей с элайсингом, рассчитывают обратное пространственное преобразование Фурье наибольшей невзвешенной составляющей волновых чисел и вычитают рассчитанную обратно-преобразованную наибольшую составляющую волновых чисел из частотной составляющей с элайсингом для формирования откорректированной частотной составляющей с элайсингом.
  17. 17. Способ по п.16, в котором применяют обратное пространственное преобразование Фурье к выбранной частотной составляющей с элайсингом для формирования интерполированных сейсмических данных вновь в частотно-пространственной области.
  18. 18. Способ по п.17, в котором обратное пространственное преобразование Фурье представляет со
    - 10 014282 бой обратное дискретное преобразование Фурье.
  19. 19. Способ по п.18, в котором обратное пространственное преобразование Фурье представляет собой обратное неравномерное быстрое преобразование Фурье.
  20. 20. Способ по п.17, в котором применяют обратное временное преобразование Фурье к выбранной частотной составляющей с элайсингом для формирования интерполированных сейсмических данных вновь в пространственно-временной области.
  21. 21. Способ по п.20, в котором обратное преобразование Фурье представляет собой обратное быстрое преобразование Фурье.
EA200900300A 2008-03-17 2009-03-11 Геофизический способ разведки EA014282B1 (ru)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US12/077,108 US7751277B2 (en) 2008-03-17 2008-03-17 Method for interpolating seismic data by anti-alias, anti-leakage Fourier transform

Publications (2)

Publication Number Publication Date
EA200900300A1 EA200900300A1 (ru) 2009-10-30
EA014282B1 true EA014282B1 (ru) 2010-10-29

Family

ID=40848502

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200900300A EA014282B1 (ru) 2008-03-17 2009-03-11 Геофизический способ разведки

Country Status (11)

Country Link
US (1) US7751277B2 (ru)
EP (1) EP2103959B1 (ru)
CN (1) CN101539634B (ru)
AU (1) AU2009200673B2 (ru)
BR (1) BRPI0900890B1 (ru)
CA (1) CA2658300C (ru)
EA (1) EA014282B1 (ru)
EG (1) EG26391A (ru)
MX (1) MX2009002932A (ru)
MY (1) MY150168A (ru)
SG (1) SG155833A1 (ru)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8014950B2 (en) * 2008-06-19 2011-09-06 Chevron U.S.A. Inc. System and method for seismic trace analysis
US8321134B2 (en) * 2008-10-31 2012-11-27 Saudi Arabia Oil Company Seismic image filtering machine to generate a filtered seismic image, program products, and related methods
US8239135B2 (en) * 2009-05-07 2012-08-07 Pgs Geophysical As Method for calculation of seismic attributes from seismic signals
US8619498B2 (en) * 2010-09-24 2013-12-31 CGGVeritas Services (U.S.) Inc. Device and method for calculating 3D angle gathers from reverse time migration
US9043155B2 (en) 2010-10-07 2015-05-26 Westerngeco L.L.C. Matching pursuit-based apparatus and technique to construct a seismic signal using a predicted energy distribution
US8862408B2 (en) * 2011-09-28 2014-10-14 Westerngeco L.L.C. Determining one or more target positions in an acquisition domain for processing survey data
US9541659B2 (en) * 2011-11-18 2017-01-10 Westerngeco L.L.C. Noise removal from 3D seismic representation
RU2488145C1 (ru) * 2012-01-10 2013-07-20 Министерство образования и науки РФ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Уральский государственный горный университет" Способ построения сейсмических изображений геологической среды
US9423518B2 (en) * 2012-02-09 2016-08-23 Pgs Geophysical As Method for processing dual-sensor streamer data with anti-alias protection
CN102636811B (zh) * 2012-04-10 2014-01-29 恒泰艾普石油天然气技术服务股份有限公司 一种海上二维地震资料中多次波的消除方法
US20140121977A1 (en) 2012-11-01 2014-05-01 Pgs Geophysical As Methods and systems for monitoring a petroleum reservoir
EP2784551A3 (en) * 2013-03-26 2015-10-28 CGG Services SA System and method for interpolating seismic data by matching pursuit in fourier transform
EP3004941A1 (en) * 2013-05-29 2016-04-13 CGG Services SA Processing of multi-sensor streamer data
CN104459770B (zh) * 2013-09-24 2017-06-16 中国石油化工股份有限公司 一种高维地震数据规则化方法
US9435903B2 (en) * 2013-10-11 2016-09-06 Chevron U.S.A. Inc. System and method for regularizing seismic data
US20150276955A1 (en) * 2013-11-06 2015-10-01 Robert H. Brune Method and System for Extending Spatial Wavenumber Spectrum Of Seismic Wavefields On Land Or Water Bottom Using Rotational Motion
US10261207B2 (en) 2014-12-18 2019-04-16 Pgs Geophysical As Seismic noise mitigation system and method
WO2016100797A1 (en) 2014-12-18 2016-06-23 Conocophillips Company Methods for simultaneous source separation
EP3109646A1 (de) * 2015-06-23 2016-12-28 Siemens Aktiengesellschaft Verfahren zur analyse eines signals sowie vorrichtung zur durchführung des verfahrens
EP3356974B1 (en) 2015-09-28 2021-12-15 ConocoPhillips Company 3d seismic acquisition
US10809402B2 (en) 2017-05-16 2020-10-20 Conocophillips Company Non-uniform optimal survey design principles
EP3714294B1 (en) * 2017-11-20 2024-01-03 Shearwater Geoservices Software Inc. Offshore application of non-uniform optimal sampling survey design
WO2019100068A1 (en) * 2017-11-20 2019-05-23 Conocophillips Company Offshore application of non-uniform optimal sampling survey design
CN108802820B (zh) * 2018-05-28 2019-10-11 中国石油天然气股份有限公司 一种深度域反假频方法、装置及系统
CN109001800B (zh) * 2018-07-20 2020-03-10 中国石油天然气股份有限公司 一种基于地震数据的时频分解与气藏检测方法及系统
CN109188535A (zh) * 2018-09-18 2019-01-11 中国科学院地质与地球物理研究所 地球物理数据处理的方法和装置
WO2020069143A1 (en) 2018-09-30 2020-04-02 Conocophillips Company Machine learning based signal recovery
US11346971B2 (en) * 2019-06-26 2022-05-31 Saudi Arabian Oil Company Imaging subterranean features using Fourier transform interpolation of seismic data
US11215725B2 (en) * 2019-07-17 2022-01-04 Saudi Arabian Oil Company Seismic processing workflow for orthogonal wide azimuth 3D surveys
CN113341220B (zh) * 2021-08-05 2021-11-02 中国空气动力研究与发展中心设备设计与测试技术研究所 含噪多频衰减实信号频率估计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4628492A (en) * 1984-01-11 1986-12-09 Mobil Oil Corporation Method of avoiding aliasing in slant stacking of seismic data
EP0551210A2 (en) * 1992-01-10 1993-07-14 Western Atlas International, Inc. Seismic surveying and interpolation of aliased seismic traces
US20020051003A1 (en) * 1998-09-21 2002-05-02 Michael Cosman Anti-aliased, textured, geocentric and layered fog graphics display method and apparatus

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4594693A (en) * 1983-11-04 1986-06-10 Mobil Oil Corporation Seismic trace interpolation using f-k filtering
US4922465A (en) * 1989-05-30 1990-05-01 Geco A/S Interpolation of severely aliased events
GB9320540D0 (en) 1993-10-06 1993-11-24 Ensign Geophysics Ltd Seismic data acquisition
US5677892A (en) * 1996-08-14 1997-10-14 Western Atlas International, Inc. Unaliased spatial trace interpolation in the f-k domain
US5617372A (en) * 1996-08-14 1997-04-01 Western Atlas International, Inc. Unaliased spatial trace interpolation in the f-k domain
US6115726A (en) * 1997-10-03 2000-09-05 Kromos Technology, Inc. Signal processor with local signal behavior
US7027929B2 (en) * 2003-11-21 2006-04-11 Geo-X Systems Ltd. Seismic data interpolation system
US7239578B2 (en) 2005-03-03 2007-07-03 John M. Robinson Removal of noise from seismic data using radon transformations

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4628492A (en) * 1984-01-11 1986-12-09 Mobil Oil Corporation Method of avoiding aliasing in slant stacking of seismic data
EP0551210A2 (en) * 1992-01-10 1993-07-14 Western Atlas International, Inc. Seismic surveying and interpolation of aliased seismic traces
US20020051003A1 (en) * 1998-09-21 2002-05-02 Michael Cosman Anti-aliased, textured, geocentric and layered fog graphics display method and apparatus

Also Published As

Publication number Publication date
EG26391A (en) 2013-09-26
CN101539634B (zh) 2014-02-05
US7751277B2 (en) 2010-07-06
AU2009200673A1 (en) 2009-10-01
CA2658300A1 (en) 2009-09-17
CN101539634A (zh) 2009-09-23
SG155833A1 (en) 2009-10-29
AU2009200673B2 (en) 2013-08-29
BRPI0900890B1 (pt) 2020-10-13
CA2658300C (en) 2014-10-28
BRPI0900890A2 (pt) 2010-04-06
EP2103959B1 (en) 2019-10-30
MY150168A (en) 2013-12-13
EP2103959A2 (en) 2009-09-23
EP2103959A3 (en) 2011-01-19
MX2009002932A (es) 2009-09-24
US20090231956A1 (en) 2009-09-17
EA200900300A1 (ru) 2009-10-30

Similar Documents

Publication Publication Date Title
EA014282B1 (ru) Геофизический способ разведки
US20210141114A1 (en) Method and system for separating seismic sources in marine simultaneous shooting acquisition
AU2016204073B2 (en) Method for separating seismic sources in marine seismic surveys
EP2372400B1 (en) Method of imaging the earth's subsurface during marine seismic data acquisition
US9151856B2 (en) Separating interfering signals in seismic data
EA017293B1 (ru) Способ очистки данных морской сейсмоприемной косы с нерегулярным расположением приемников от ложных отражений
EP2322956B1 (en) Method for full-bandwidth deghosting of marine seismic streamer data
US9448317B2 (en) Method for swell noise detection and attenuation in marine seismic surveys
EP2336809A2 (en) Method for Attenuating Interference Noise in Dual-Sensor Seismic Data
EA022531B1 (ru) Способ определения сейсмического атрибута по сейсмическим сигналам
EP2330443A2 (en) Method for full-bandwidth source deghosting of marine seismic streamer data
US10386518B2 (en) Device and method for deblending simultaneous shooting data using an iterative approach
US20220342103A1 (en) Noise Attenuation Methods Applied During Simultaneous Source Deblending and Separation

Legal Events

Date Code Title Description
MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): AM BY KG MD TJ TM

QB4A Registration of a licence in a contracting state
QZ4A Registered corrections and amendments in a licence
MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): AZ KZ RU