EA026615B1 - Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке - Google Patents

Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке Download PDF

Info

Publication number
EA026615B1
EA026615B1 EA201490431A EA201490431A EA026615B1 EA 026615 B1 EA026615 B1 EA 026615B1 EA 201490431 A EA201490431 A EA 201490431A EA 201490431 A EA201490431 A EA 201490431A EA 026615 B1 EA026615 B1 EA 026615B1
Authority
EA
Eurasian Patent Office
Prior art keywords
duration
time delay
time
sources
activation
Prior art date
Application number
EA201490431A
Other languages
English (en)
Other versions
EA201490431A1 (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 EA201490431A1 publication Critical patent/EA201490431A1/ru
Publication of EA026615B1 publication Critical patent/EA026615B1/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3861Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas control of source arrays, e.g. for far field control
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02Generating seismic energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Unknown Time Intervals (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Farming Of Fish And Shellfish (AREA)

Abstract

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

Description

В течение нескольких последних десятилетий в нефтегазовой отрасли значительные средства инвестируются в разработку методов морской сейсморазведки, позволяющих получать сведения о подземных формациях, расположенных ниже водного слоя, в целях обнаружения и извлечения ценных полезных ископаемых, таких как нефть. Высокоразрешающие сейсмические изображения подземной формации имеют существенное значение для количественной интерпретации сейсморазведочных данных и мониторинга продуктивных пластов. При обычной морской сейсморазведке сейсморазведочное судно буксирует сейсмический источник, при этом то же или другое судно буксирует одну или несколько сейсмических кос, образующих поверхность сбора сейсмических данных ниже водной поверхности и выше подземной формации, исследуемой на наличие месторождений полезных ископаемых. Судно везет оборудование для сбора сейсмических данных, такое как оборудование навигационного управления, оборудование управления сейсмическими источниками, оборудование управления сейсмическими приемниками и записывающее оборудование. Оборудование управления сейсмическими источниками активирует сейсмоисточник, обычно представляющий собой группу элементов-источников, таких как пневмопушки, инициируя генерирование акустических импульсов в выбранные моменты времени. Сейсмический источник часто описывают как возбуждаемый или активируемый с целью генерации акустического импульса. Каждый акустический импульс представляет собой звуковую волну, распространяющуюся в воде в нижнем направлении и проходящую в подземную формацию. На каждой границе раздела между различными типами пород часть звуковой волны передается, а другая часть отражается обратно к водному слою, распространяясь в направлении поверхности. Сейсмические косы, буксируемые позади судна, представляют собой удлиненные кабельные конструкции. Каждая сейсмическая коса содержит некоторое количество сейсмических приемников или многокомпонентных датчиков, регистрирующих волновые поля давления и(или) движения частиц, связанные со звуковыми волнами, отражающимися обратно к водному слою от подземной формации. При традиционном сборе сейсмических данных сейсмические данные, формируемые приемниками, записываются в течение конечного периода времени и рассматриваются в качестве исходящих из единого источника. При сборе данных от источника одновременного действия несколько источников возбуждают в пределах небольшого временного окна, называемого интервалом активации, при этом записанные сейсмические данные нельзя считать исходящими из единого источника: их рассматривают в качестве исходящих из нескольких источников и записываемых с перекрыванием по времени. При реализации способов обработки следует принимать во внимание перекрывание этих записей. Один из таких способов состоит в активном разделении записываемых данных по источникам, при этом для всех сейсмограмм можно предположить, что акустическая волна исходит из единого источника. В нефтегазовой отрасли ведется поиск систем и способов, снижающих сейсмические помехи в сейсмических данных, формируемых приемниками при сборе данных от источника одновременного действия.
Перечень фигур
На фиг. 1А-1В показаны вид сверху и вертикальный вид сбоку одного из вариантов морской системы сбора сейсмических данных.
На фиг. 2 показаны путь судна и линия времени морской съемки для морской системы сбора сейсмических данных.
На фиг. 3 показана гистограмма примера временных задержек со случайными начальными значениями для восьми источников.
На фиг. 4 показано представление временных задержек с помощью круговых временных шкал.
На фиг. 5А показан интервал активации, разделенный на группы по длительности.
На фиг. 5В показан пример списка временных интервалов групп длительности.
На фиг. 6 показано случайное упорядочение групп длительности.
На фиг. 7 показан круговой график случайно упорядоченных групп длительности, представленных на фиг. 6.
На фиг. 8 показана линия времени морской съемки с первым интервалом активации, увеличенным, чтобы продемонстрировать последовательность возбуждения источников.
На фиг. 9 показан пример сдвига временных задержек источников.
На фиг. 10 показана линия времени морской съемки с первым и вторым интервалами активации, увеличенными, чтобы продемонстрировать последовательность возбуждения источников.
На фиг. 11 показан пример сдвига временных задержек источников.
На фиг. 12 показана линия времени морской съемки с тремя последовательными интервалами возбуждения, увеличенными, чтобы продемонстрировать последовательность возбуждения источников.
На фиг. 13 показана блок-схема способа рандомизации времени возбуждения источника одновременного действия.
На фиг. 14 показана блок-схема способа сдвига временных задержек, вызываемых в блоке 1308 блок-схемы, представленной на фиг. 13.
На фиг. 15 показан иллюстративный пример обобщенной компьютерной системы, реализующей способ рандомизации времени возбуждения источников одновременного действия.
- 1 026615
Сведения, подтверждающие возможность осуществления изобретения.
В настоящем описании изобретения представлены вычислительные системы и способы оптимизации рандомизации порядка возбуждения нескольких источников при сборе данных от источника одновременного действия СДИОД (§§А). Несколько источников возбуждаются в интервалах активации, разделенных периодами ожидания, во время которых источники не возбуждаются. Раскрываемые ниже системы и способы относятся к псевдорандомизации порядка возбуждения источников в пределах каждого интервала активации. Порядок возбуждения нескольких источников в пределах интервала активации псевдорандомизирован так, чтобы способы обработки данных приемника позволяли оптимальным образом разделять записанные данные из сейсмограмм отдельных источников. Иными словами, эти способы и системы псевдорандомизируют время возбуждения нескольких источников в интервалах активации таким образом, чтобы каждый источник возбуждался в другой момент времени, и чтобы разность времен между возбуждением двух любых последовательных источников была больше, чем заданный параметр.
При традиционном сборе сейсмических данных сейсмические данные, формируемые приемниками и записываемые в течение конечного периода времени, рассматриваются в качестве поступающих от единого источника. При СДИОД (§§А) сейсмические волновые поля, генерируемые несколькими источниками, участвуют в формировании данных за конечный период времени одним и тем же набором приемников. Альтернативно, данные приемников могут формироваться и записываться непрерывно, при этом записанные волны от различных источников интерферируют друг с другом. Цель СДИОД (§§А) заключается в уменьшении времени, затрачиваемого на сбор сейсмических данных, и(или) увеличении разнообразия данных по таким показателям, как кратность, азимут и удаления.
Один подход к обработке данных приемников заключается в разделении источников, при котором сейсмограммы по данным единственного источника получают так, как если бы между источниками не существовало никакой интерференции. Такой подход, называемый активным разделением, обеспечивает преимущества за счет произвольного возбуждения нескольких источников. При активном разделении один из нескольких источников выбирают в качестве первичного, при этом данные приемников согласуют с временем возбуждения первичного источника и упорядочивают в соответствующей области, такой как область общей точки приема или область равных удалений. В идеальном случае рандомизация времени возбуждения для нескольких источников проявляет себя в том, что волны от оставшихся вторичных источников будут интерферировать и выглядеть некогерентными, тогда как данные приемников от первичного источника будут представляться в качестве когерентной (регулярной) волны. Вследствие этого можно использовать способы, разделяющие источники во времени на основе когерентной фильтрации. Однако стратегии рандомизации временных задержек для источников при СДИОД (§§А) могут привести к появлению волны от вторичного источника, представляющейся когерентной в составе сейсмических данных, формируемых приемниками, что может произойти, когда среди некоторого количества интервалов активации два или несколько времен задержки для любых источников одновременного действия идентичны или близки друг к другу. В результате при использовании методов разделения на основе когерентности фильтр когерентности, действующий в пределах некоторого пространственновременного окна, оказывается не в состоянии определить, что волны от вторичных источников являются некогерентными, вследствие чего полученные результаты будут неоптимальными.
В отличие от этого в описанных ниже вычислительных системах и способах применяется методика рандомизации нескольких источников, согласно которой порядок или последовательность возбуждения нескольких источников рандомизируется от одного интервала активации до другого, чтобы избежать когерентности волн, создаваемых вторичными источниками. Рандомизированный порядок возбуждения для нескольких источников гарантирует, что в пределах длины оператора методики разделения разность времен возбуждения между последовательными трассами в пространственно-временном окне достаточно велика, чтобы волна от вторичного источника стала некогерентной.
На фиг. 1А-1В представлены соответственно вид сверху и вертикальный вид сбоку одного из вариантов морской системы сбора сейсмических данных, состоящей из четырех сейсморазведочных судов 101-104. Как показано на чертеже, каждое из сейсморазведочных судов буксирует по меньшей мере один из восьми отдельных источников, обозначенных §1, §2, §3, §4, §5, §6, §7 и §8. Источники §1-88 представляют собой пример источника одновременного действия, работа которого более подробно раскрыта ниже. Конкретно, сейсморазведочное судно 101 буксирует источник §1, сейсморазведочное судно 102 буксирует источники §2 и §3, сейсморазведочное судно 103 буксирует источники §4 и §5, а сейсморазведочное судно 104 буксирует источники §6, §7 и §8. Другие схемы проведения сейсморазведки могут предусматривать различные другие количества и варианты распределения источников между сейсморазведочными судами, количество которых также может быть различным. На фиг. 1А сейсморазведочное судно 104 также буксирует шесть отдельных сейсмических кос 106-111, размещенных в водоеме. В этом примере один конец каждой сейсмической косы присоединен к сейсморазведочному судну 104 при помощи кабеля для передачи данных сейсмической косы. Сейсмические косы 106-111 образуют плоскую горизонтальную поверхность сбора данных приемников по отношению к свободной поверхности водоема. Однако на практике поверхность сбора данных приемников может плавно изменяться под воздействием активных морских течений и погодных условий. Иными словами, хотя сейсмические косы 106-111
- 2 026615 показаны прямыми, на практике буксируемые сейсмические косы могут совершать волнообразные колебания вследствие динамических условий водоема, в который погружены сейсмические косы. Следует отметить, что поверхность сбора данных приемников не ограничена горизонтальной ориентацией по отношению к свободной поверхности 112. Сейсмические косы могут буксироваться на глубинах, ориентирующих поверхность сбора данных приемников по некоторым углом относительно свободной поверхности 112, или таким образом, чтобы одна или несколько сейсмических кос находились на различных глубинах. Следует также отметить, что поверхность сбора данных приемников не ограничивается шестью сейсмическими косами. На практике поверхности сбора данных приемников могут состоять из различного количества сейсмических кос, от одной до 20 и более. Некоторые варианты осуществления содержат одну или несколько сейсмических кос, буксируемых одним или несколькими сейсморазведочными судами 101-104.
Фиг. 1А содержит плоскость ху 112, а фиг. 1В - плоскость χζ 114 декартовой системы координат, имеющей три ортогональные пространственные оси, обозначенные х, у и ζ. Система координат используется для указания ориентации и расположения координат в водоеме. Направление оси х указывает положение точки в направлении, параллельном длине сейсмических кос, и называется продольным направлением. Направление оси у указывает положение точки в направлении, перпендикулярном оси х и практически параллельном свободной поверхности, и называется поперечным направлением. Направление оси ζ указывает положение точки, перпендикулярное плоскости ху (т. е. перпендикулярное свободной поверхности), причем положительная ось ζ направлена вниз относительно свободной поверхности. Глубину сейсмических кос ниже свободной поверхности можно оценивать в различных точках вдоль сейсмических кос при помощи глубиномеров, прикрепленных к сейсмическим косам. Например, глубиномеры могут измерять гидростатическое давление или использовать акустические измерения расстояний. Глубиномеры могут объединяться с регуляторами глубины погружения, такими как параваны или отводители, управляющими глубиной и положением сейсмических кос во время буксирования сейсмических кос в водоеме. Глубиномеры, как правило, располагаются с интервалами около 300 м вдоль каждой сейсмической косы. Отметим, что в некоторых вариантах осуществления буи, прикрепленные к сейсмическим косам 106-111, могут использоваться для сохранения ориентации и глубины сейсмических кос ниже свободной поверхности.
На фиг. 1А заштрихованные прямоугольники 116 представляют приемники или датчики, разнесенные по длине каждой сейсмической косы. Сейсмические косы 106-111 представляют собой длинные кабели, содержащие линии питания и линии передачи данных, присоединяющие приемники 116 к оборудованию для сбора сейсмических данных, размещенному на борту сейсморазведочного судна 104. При одном типе сейсмической разведки каждый приемник представляет собой многокомпонентный датчик, содержащий один или несколько геофонов, каждый из которых регистрирует смещение в водоеме во времени путем регистрации движения частиц, скоростей или ускорений, и гидрофон, регистрирующий изменения давления воды во времени. Сейсмические косы 106-111 и сейсморазведочные суда 101-104 содержат специальные электронные измерительные устройства и оборудование обработки данных, позволяющие соотносить показания приемников с абсолютным положением на свободной поверхности и абсолютным пространственным положением относительно произвольной трехмерной системы координат.
Каждый из источников §1-88 состоит из одного или нескольких элементов источника, таких как гидропушки, пневмопушки или морские вибраторы, подвешенные на поплавках ниже свободной поверхности. Пневмопушки и гидропушки выбирают таким образом, чтобы они имели различные размеры камер, и специально располагают в пределах группы так, чтобы генерировать результирующий акустический импульс. Морские вибраторы могут использоваться в режиме ступенчатой функции для генерации результирующего акустического импульса. В другие вариантах осуществления источники 81-88 могут представлять собой источники электромагнитного (ЭМ) поля. Источник ЭМ поля содержит два возбуждающих электрода, расположенных на противоположных концах кабеля для образования электродипольной передающей антенны, подключаемой к сейсморазведочному судну при помощи возбуждающего вводного кабеля. Возбуждающие электроды ЭМ поля и кабель могут располагаться горизонтально или вертикально. На фиг. 1В показан вертикальный вид сбоку сейсморазведочных судов 101-104, буксирующих источники 81-88 ниже свободной поверхности 118 и выше подземной формации 120. Сейсмические косы 106-111, буксируемые сейсморазведочными судами 102-104, не показаны для удобства. Кривая 122 представляет твердую поверхность на дне водоема, расположенного выше подземной формации 120. Подземная формация 120 представляет собой некоторое количество подземных слоев, сложенных из осадков и пород. Кривые 124, 126 и 128 представляют собой границы раздела между подземными слоями различного состава. Заштрихованная область 130, ограниченная сверху кривой 132, а снизу кривой 134, представляет подземную залежь углеводородов, глубина и позиционные координаты которой могут определяться путем анализа сейсмических данных, собранных в ходе морской сейсморазведки. По мере продвижения сейсморазведочных судов 101-104 над подземной формацией 120 источники 81-88 формируют кратковременные волны давления, называемые акустическими импульсами, с пространственными и временными интервалами, подробнее описанными ниже. Для удобства на фиг. 1В
- 3 026615 изображен акустический импульс, распространяющийся наружу от источника 81 в виде волнового поля 136 давления, представленного полуокружностями возрастающего радиуса с центром в источнике 81. Волновые фронты, по существу, показаны на фиг. 1В в поперечном разрезе вертикальной плоскости. Часть волнового поля 136 давления, распространяющаяся наружу в нижнем направлении, называется первичным волновым полем, в конечном счете, достигающим твердой поверхности 122 подземной формации 120, причем в этой точке первичное волновое поле частично отражается от твердой поверхности 122, а частично преломляется вниз в твердую подземную формацию 120, становясь упругими волнами внутри подземной формации 120. Иными словами, в водоеме акустический импульс состоит из продольных волн давления или Р-волн, тогда как в твердой подземной формации 120 волны содержат как Рволны, так и сдвиговые волны, или 8-волны. В пределах подземной формации 120 на каждой границе раздела между материалами различных типов или при скачкообразном изменении плотности или одной или нескольких различных других физических характеристик или параметров, нисходящие волны частично отражаются и частично преломляются. Вследствие этого каждая точка поверхности 122 и каждая точка границ раздела 124, 126 и 128 становится потенциальным вторичным точечным источником (пунктом возбуждения), из которого акустические и упругие волны соответственно могут исходить вверх в направлении сейсмических кос (не показаны), буксируемых сейсморазведочным судном 104, в ответ на акустический импульс, возбуждаемый источником 81, и нисходящие упругие волны, возбуждаемые импульсом давления. Как показано на фиг. 1В, вторичные волны значительной амплитуды обычно распространяются из точек, расположенных на твердой поверхности 122 или вблизи от нее, таких как точка 138, а также из точек, расположенных на границах раздела в подземной формации 120 или в непосредственной близости от них, таких как точки 140 и 142. Третичные волны, называемые волнами-спутниками со стороны приемника, порождаются вторичными волнами, отражающимися от свободной поверхности 118 назад в направлении сейсмических кос и подземной формации 120.
Вторичные волны обычно испускаются в различные моменты времени в диапазоне времени, следующем за первоначальным акустическим импульсом. Точка на твердой поверхности 120, такая как точка 138, воспринимает возмущение давления, соответствующее первоначальному акустическому импульсу, быстрее, чем точка внутри подземной формации 120, такая как точки 140 и 142. Аналогичным образом, точка на твердой поверхности 122, расположенная непосредственно под источником 81, воспринимает акустический импульс давления быстрее, чем более отдаленная точка на твердой поверхности 122. Таким образом, моменты времени, в которые вторичные волны и волны более высокого порядка испускаются из различных точек внутри подземной формации 120, связаны с расстоянием в трехмерном пространстве, на которое эти точки отстоят от источника 81. Следует отметить, что акустические импульсы, излучаемые другими источниками 82-88, распространяются наружу в качестве волновых полей давления, взаимодействующих с подземной формацией 120, таким же образом, как акустический импульс, излучаемый описанным выше источником 81.
Однако акустические и упругие волны распространяются с различными скоростями в различных материалах, а также в одном материале при разных давлениях. Поэтому времена пробега первичного волнового поля и вторичного волнового поля, испускаемого в ответ на первичное волновое поле, представляют собой комплексные функции расстояния от источника 104, а также материалов и физических характеристик материалов, через которые проходит первичная волна. Кроме того, вторичные расширяющиеся волновые фронты могут изменяться при пересечении волновыми фронтами границ раздела и при изменении скорости звука в средах, через которые проходит волна. Суперпозиция волн, испускаемых из подземной формации 120 в ответ на первичное волновое поле, обычно представляет собой очень сложное волновое поле, содержащее информацию о формах, размерах и характеристиках материалов подземной формации 120, в том числе информацию о формах, размерах и местах расположения различных отражающих объектов в подземной формации 120, представляющей интерес для сейсморазведчиков.
Подземные формации, расположенные под водоемом, можно также исследовать с использованием методик морской донной съемки. Один из вариантов этих методик реализуют при помощи морских донных кос МДК (ОВС). МДК (ОВС) аналогичны буксируемым сейсмическим косам, описанным выше, в том отношении, что МДК (ОВС) содержат некоторое количество пространственно разнесенных приемников, например приемников, размещаемых через каждые 25-50 м, однако МДК (ОВС) проложены по подошве 122 водного слоя или близи от нее. МДК (ОВС) могут быть подключены при помощи электронного соединения к стоящему на якоре регистрирующему судну, обеспечивающему питание, подачу команд и управление измерительной аппаратурой, а также телеметрическую передачу данных датчиков на записывающее оборудование, расположенное на борту судна. Альтернативно, методики морской донной съемки могут быть реализованы при помощи автономных систем, состоящих из приемных узлов, размещаемых и извлекаемых при помощи дистанционно управляемых аппаратов. Приемные узлы могут располагаться на подошве 122 водного слоя или вблизи от нее, например, на расстоянии 400 м друг от друга, образуя довольно грубую сетку. Автономные системы приемников, как правило, реализуют при помощи одного из двух типов систем приемных узлов. Первая система приемных узлов представляет собой кабельную систему, в которой приемные узлы соединены друг с другом кабелями и присоединены к стоящему на якоре регистрирующему судну. Кабельные системы снабжены питанием, поступающим на каж- 4 026615 дый приемный узел по кабелю, а сейсмические данные возвращаются на регистрирующее судно по кабелю или при помощи радиотелеметрии. Во второй системе приемных узлов применяются автономные приемные узлы, имеющие источник питания ограниченной мощности, однако для загрузки записанных сейсмических данных приемные узлы обычно приходится поднимать со дна. Используя как МДК (ОВС), так и автономные приемные узлы, сейсмические суда, оснащенные двумя или несколькими сейсмоисточниками, как описано выше со ссылкой на фиг. 1А и 1В, могут применяться для возбуждения акустических импульсов с пространственными и временными интервалами при перемещении судов возбуждения по свободной поверхности.
Количество источников одновременного действия не ограничивается восемью источниками §1-88, буксируемыми четырьмя сейсморазведочными судами 101-104, как показано на фиг. 1А и 1В. На практике количество источников, составляющих источник одновременного действия, может варьировать в диапазоне от всего двух до 10 или большего количества источников, при этом источники могут распределяться между любым подходящим количеством сейсморазведочных судов. В целях удобства и краткости для раскрытия способов и систем управления источником одновременного действия в псевдоодновременном режиме с целью достижения временного перекрытия при записи волновых полей, генерируемых источниками, используется пример комплекта из восьми источников 81-88. В целом, раскрываемые ниже способы могут распространяться на η источников одновременного действия, где η - целое число, большее 1.
На фиг. 2 показана морская система сбора сейсмических данных, изображенная на фиг. 1А-1В, проходящая путь судна, представленный стрелкой-указателем 202. Стрелка-указатель 204 представляет линию времени морской съемки, связанную с работой источников 81-88 по мере прохождения системой сбора сейсмических данных пути 202 судна. Прямоугольники представляют интервалы активации, во время которых возбуждаются источники 81-88, а двусторонние стрелки-указатели 210 представляют периоды ожидания между интервалами активации. Например, первые три интервала 206-208 активации могут иметь длительность приблизительно 1 с, при этом периоды 210 ожидания могут изменяться в диапазоне 6-8 с. По мере прохождения системой сбора сейсмических данных пути 202 судна источники 8188 возбуждаются в псевдоодновременном режиме, раскрываемом ниже, в пределах каждого интервала активации, при этом во время периодов ожидания источники 81-88 не возбуждаются. Раскрываемый ниже псевдоодновременный режим возбуждения источников 81-88 приводит к временному перекрытию при записи волновых полей, генерируемых источниками. Псевдоодновременное возбуждение источников 81-88 означает, что каждый из источников 81-88 возбуждается с уникальной рандомизированной временной задержкой. Рандомизированные временные задержки возбуждения позволяют избежать усиливающей интерференции излучаемых волновых полей. Способы генерации случайного распределения временных задержек возбуждения для источников 81-88 раскрыты ниже.
На фиг. 3 показана гистограмма примера временных задержек со случайными начальными значениями для восьми источников 81-88. Вертикальная ось 302 идентифицирует восемь источников 81-88, а горизонтальная ось 304 представляет пример интервала активации, изменяющегося в диапазоне от 0 до 1000 мс (т. е. 1 с). Все восемь источников 81-88 возбуждаются в соответствии с временными задержками, представленными заштрихованными полосами, идущими от оси 302 источников. Случайные начальные временные задержки задаются произвольно с ограничением, согласно которому случайно присваиваемые начальные значения предотвращают усиливающую интерференцию. Отметим, что временные задержки в интервалах активации неодинаковы и не могут быть ближе друг к другу, чем определенный пользователем параметр. Длительность каждой временной задержки соответствует длине каждой полосы и отмечена вдоль временной оси 304. Например, источнику 81 присвоена временная задержка 125 мс, представленная длиной полосы 306, а источнику 85 присвоена временная задержка 726 мс, представленная длиной полосы 308. Вследствие этого набор временных задержек, показанных на фиг. 3, называют случайными начальными значениями, а каждую временную задержку называют временной задержкой со случайным начальным значением. Случайные начальные значения используются для определения случайных временных задержек при возбуждении источников 81-88 в каждом интервале активации.
На фиг. 4 показан альтернативный способ представления временных задержек со случайными начальными значениями относительно интервала активации с помощью внутренней круговой временной шкалы 402 и внешней круговой временной шкалы 404. Круговые шкалы 402 и 404 используются ниже, чтобы продемонстрировать концепцию сдвига временных задержек со случайными начальными значениями и перенос назад сдвинутых временных задержек, больших, чем длительность интервала активации, раскрытую ниже. Внутренняя круговая временная шкала 402 представляет интервал активации и образована окружностью с 20 штрихами шкалы, расположенными с равными интервалами 50 мс. Штрихи шкалы обозначены моментами времени, разделенными интервалами 50 мс в пределах интервала активации 1000 мс. Внешняя круговая временная шкала 404 образована окружностью с 8 штрихами шкалы, представляющими временные задержки со случайными начальными значениями, первоначально присвоенными источникам 81-88. В следующем описании псевдослучайного упорядочения времен возбуждения источников 81-88 в пределах интервала активации внешняя круговая шкала 404 времени повернута относительно фиксированного положения внутренней круговой временной шкалы 402, чтобы проиллю- 5 026615 стрировать, как изменяются временные задержки при различных интервалах активации.
Для каждого интервала активации линии времени морской съемки порядок возбуждения источников 81-88 является псевдорандомизированным. Псевдорандоминизированный порядок возбуждения источников 81-88 в пределах интервала активации определяется добавлением псевдослучайно выбранного временного сдвига к временным задержкам со случайными начальными значениями источников 81-88. Способ псевдослучайного выбора порядка возбуждения источников 81-88 в пределах интервала активации будет раскрыт теперь со ссылкой на фиг. 5-12. Длительность интервала активации делится на длину оператора, например на длину локального окна данных, в котором данные, формируемые приемниками, фильтруются по степени когерентности. Поскольку сейсмические данные генерируются в виде трасс, длина оператора может определяться по количеству трасс в некоторой области. Например, чтобы определить результаты фильтрования для трассы х, один из примеров локального окна данных определяют от трассы х-4 до трассы х+4. По этим 9 трассам, лежащим в пределах локального окна данных, определяют результат фильтрования для трассы х. Так, для каждой трассы в полной сейсмограмме результат работы фильтра определяют по фильтру когерентности с длиной оператора, составляющей 9 трасс. В описанном ниже примере длительность интервалов активации, составляющая 0-1000 мс, делится на длину оператора, составляющую 20 трасс. На практике длина оператора может состоять из любого требуемого количества трасс. Вследствие этого интервал активации делят на 20 субвременных интервалов, называемых группами длительности, где каждая группа длительности имеет длительность около 50 мс. На фиг. 5А показан интервал активации 0-1000 мс, разделенный на 20 групп длительности, при этом каждая группа длительности обозначена дт, где т=1,2,3, ...,20 - индекс группы длительности. На фиг. 5В показан список из 20 групп длительности и временные интервалы, связанные с каждой из 20 групп длительности. Отметим, что в данном примере отношение количества групп длительности к количеству трасс на длине оператора окна данных составляет один к одному, что может подходить для области источниковприемников или области равных удалений.
После образования групп длительности их размещают в произвольном порядке с использованием псевдорандомизации. На фиг. 6 показан случайный порядок расположения 20 групп длительности, изображенных на фиг. 5А, в направлении, представленном стрелкой-указателем 602. Например, группа д7 длительности произвольно выбрана в качестве первой группы длительности, группа д12 длительности произвольно выбрана в качестве второй группы длительности и так далее, при этом группа д17 длительности произвольно выбрана в качестве последней группы длительности. Затем группы длительности располагают в виде кругового графика, представленного на фиг. 7. На фиг. 7 показана окружность 702 с 20 точками, расположенными с равными интервалами и обозначенными по часовой стрелке произвольно упорядоченными номерами групп длительности, изображенных на фиг. 6. Например, точки 706-708 идентифицированы с помощью первых трех групп д7, д12 и д5 длительности из числа псевдослучайно упорядоченных групп длительности, изображенных на фиг. 6, а точка 710 идентифицирована с помощью последней группы д17 длительности из числа псевдослучайно упорядоченных групп длительности, изображенных на фиг. 6.
Круговой график псевдослучайно упорядоченных групп длительности, показанный на фиг. 7, и случайные начальные значения, показанные на фиг. 4, используются в сочетании для определения произвольной последовательности возбуждения источников 81-88 в каждом интервале активации на линии времени морской съемки в соответствии со следующим. Для первого интервала активации на линии времени морской съемки идентифицируют наименьшее время задержки из числа случайных начальных значений и группу длительности, в которую попадает наименьшее время задержки. Например, случайные начальные значения, представленные на фиг. 3 и 4, демонстрируют, что источнику 81 присвоено наименьшее время задержки 125 мс, попадающее в группу д3 длительности, как показывает список субвременных интервалов групп длительности, представленный на фиг. 5В. Вследствие этого, возбуждение источников 81-88 в первом интервале активации на круговом графике, показанном на фиг. 7, начинается в группе д3 длительности. На фиг. 7 группы длительности, присвоенные первым трем интервалам активации, идентифицированы обведенными номерами 1, 2 и 3, при этом группа д3 обозначена обведенным номером 1, показывающим, что группа д3 длительности содержит источник с наименьшим временем задержки.
В начале первого интервала активации источники 81-88 возбуждаются в соответствии с временными задержками со случайными начальными значениями. На фиг. 8 показана линия 204 времени морской съемки с первым интервалом 206 активации, увеличенным, чтобы продемонстрировать последовательность возбуждения источников 81-88. В этом примере источники 81-88 возбуждаются в соответствии с временными задержками со случайными начальными значениями после начала интервала 206 активации. Например, источник 81 возбуждается спустя 125 мс после начала интервала активации, затем, спустя 212 мс после начала интервала активации, происходит возбуждение источника 87, вслед за этим в соответствии с временными задержками со случайными начальными значениями возбуждаются источники 84, 82, 86, 85 и 88 пока, спустя 920 мс после начала интервала активации, не произойдет возбуждение источника 83. После того как произойдет возбуждение всех источников 81-88, источники сбрасываются и по истечении интервала активации длительностью 1000 мс до начала второго интервала 207 активации про- 6 026615 должается период ожидания.
Для второго интервала 207 активации временная задержка для источника с наименьшей случайной начальной временной задержкой выбирается из соседней группы длительности на круговом графике, показанном на фиг. 7. Например, источник 81 имеет наименьшую случайную начальную временную задержку, при этом ему присваивается произвольно выбранная временная задержка 115 в группе §15 длительности. Иными словами, источнику 81 может быть присвоено любое время 115 во временном интервале 701-750 мс. Временные задержки для оставшихся источников 82-88 сдвигают, добавляя разность между временной задержкой 115 и случайной начальной временной задержкой 125 мс для источника 81 (т. е. ΐι5-125 мс) к каждой из случайных начальных временных задержек оставшихся источников 82-88. Если сдвинутое время задержки больше, чем длительность интервала активации, сдвинутое время задержки регулируют, вычитая длительность интервала активации так, чтобы сдвинутое время задержки попадало внутрь интервала активации. Предположим, что временная задержка ®5, выбранная для источника 81 в группе §15 длительности, составляет 725 мс. Разность между выбранным временем задержки ®5-725 мс для источника 81 и случайной начальной временной задержкой для источника 81 равна 600 мс (т. е. 725125 мс). Каждая из случайных начальных временных задержек для оставшихся источников 82-88 также сдвигается на 600 мс. Например, сдвинутые временные задержки для источников 87 и 84 составляют 812 мс (т. е. 212+600 мс) и 902 мс (302+600 мс) соответственно. Однако добавление 600 мс к случайным начальным временным задержкам для источников 82, 86, 85, 88 и 83 дает сдвинутые временные задержки, превышающие 1000 мс и, таким образом, выходящие за пределы длительности интервала активации. Сдвинутые временные задержки, превышающие интервал активации, корректируют путем вычитания длительности интервала активации, при этом сдвинутые временные задержки помещают обратно в интервал активации, сохраняя значения разности между моментами времени активации. Иными словами, сдвинутые временные задержки, превышающие длительность интервала активации, возвращают назад в интервал активации, а также называют возвращаемыми временными задержками. Например, источник 82 имеет случайную начальную временную задержку 420 мс, сдвинутую к 1020 мс (т. е. 420+600 мс). Поскольку сдвинутое время задержки 1020 мс превышает 1000 мс, сдвинутое время задержки возвращают назад в интервал активации, вычитая 1000 мс, чтобы получить возвращаемую временную задержку 20 мс.
На фиг. 9 показан пример сдвига временных задержек источников 81-88 для второго интервала 207 активации с использованием круговых временных шкал 402 и 404, описанных выше со ссылкой на фиг.
4. Присвоение сдвинутых временных задержек источникам 81-88 эквивалентно повороту внешней временной шкалы 404 в направлении по часовой стрелке на 600 мс, как показано стрелкой-указателем 902. Места расположения штрихов шкалы, связанных с каждым из источников 81-88, относительно внутренней временной шкалы 402, можно использовать для присвоения сдвинутых временных задержек во втором интервале 207 активации. Сдвинутые временные задержки, попадающие в область внешней временной шкалы 404, идентифицируемую штриховой стрелкой-указателем 904, также представляют собой возвращаемые временные задержки. В начале второго интервала активации источники 81-88 возбуждаются в соответствии с присвоенными им сдвинутыми временными задержками.
На фиг. 10 показана линия 204 времени морской съемки с первым и вторым интервалами 206 и 207 активации, увеличенными, чтобы продемонстрировать последовательность возбуждения источников 8188. Во втором интервале 207 активации источники 81-88 возбуждаются в соответствии с присвоенными им сдвинутыми временами задержки после начала интервала 207 активации, как показано на круговой временной шкале 404 на фиг. 9. Например, источник 82 представляет собой первый источник, возбуждаемый через 20 мс после начала интервала возбуждения, затем, спустя 148 мс после начала интервала активации, происходит возбуждение источника 86, а источники 85, 88, 83, 81 и 84 возбуждаются в соответствии со связанными с ними сдвинутыми временными задержками, представленными на фиг. 9. После того как произойдет возбуждение всех источников 81-88, источники сбрасываются и по истечении интервала активации длительностью 1000 мс до начала второго интервала 208 активации продолжается период ожидания.
Для третьего интервала 208 активации временную задержку для источника с наименьшей присваиваемой временной задержкой выбирают из следующей группы длительности на круговом графике, показанном на фиг. 7. Как описано выше со ссылкой на фиг. 3, источник 81 является источником с наименьшей случайной начальной временной задержкой, при этом ему присваивается сдвинутое время задержки 11 в группе §1 длительности, представляющей собой временной интервал 0-50 мс, как показано на фиг. 5В. Сдвинутые временные задержки для оставшихся источников 82-88 вычисляют, вычитая сначала разность между временной задержкой 11 и случайной начальной временной задержкой 125 мс для источника 81 (т. е. Д-125 мс), а затем добавляя разность к каждой из случайных начальных временных задержек, связанных с оставшимися источниками 82-88. Например, предположим, что произвольно выбранная временная задержка 11 для источника 81 составляет 25 мс, что соответствует группе §1 длительности. Разность между выбранной временной задержкой Д=25 мс для источника 81 и случайной начальной временной задержкой для источника 81 равна 100 мс (т. е. 25-125 мс). Сдвинутые временные задержки для оставшихся источников 82-88 вычисляют, добавляя разность 100 мс к каждой случайной начальной вре- 7 026615 менной задержке для источников 82-88, как описано выше со ссылкой на фиг. 3. Например, сдвинутые временные задержки для источников 87 и 84 составляют 112 мс (т. е. 212-100 мс) и 202 мс (302-100 мс) соответственно.
На фиг. 11 показан пример сдвига временных задержек источников 81-88 для третьего интервала 208 активации с использованием круговых временных шкал 402 и 404, описанных выше со ссылкой на фиг. 4. Присвоение сдвинутых временных задержек источникам 81-88 эквивалентно повороту внешней временной шкалы 208 в направлении против часовой стрелке на 100 мс, как показано стрелкойуказателем 1102. Места расположения штрихов шкалы, связанных с каждым из источников 81-88, относительно внутренней временной шкалы 402 можно использовать для присвоения сдвинутых временных задержек источникам 81-88 в третьем интервале 207 активации. В начале третьего интервала 208 активации источники 81-88 возбуждаются в соответствии с присвоенными им сдвинутыми временными задержками, представленными на фиг. 11.
На фиг. 12 показана линия 204 времени морской съемки с первым, вторым и третьим интервалами 206-208 активации, увеличенными, чтобы продемонстрировать последовательность возбуждения источников 81-88 в каждом из интервалов активации. В третьем интервале 208 активации источники 81-88 возбуждаются в соответствии с присвоенными им временными задержками после начала интервала 208 активации. Сдвинутые временные задержки представлены на внешней временной шкале 404 на фиг. 9. Например, источник 81 представляет собой первый источник, возбуждаемый через 25 мс после начала интервала активации, затем, спустя 112 мс после начала интервала активации, происходит возбуждение источника 87, а источники 84, 82, 86, 85, 88 и 83 возбуждаются в соответствии со связанными с ними сдвинутыми временными задержками, представленными на фиг. 11.
Как правило, по истечении интервала активации источники сбрасываются, после чего до начала следующего интервала активации продолжается период ожидания. Способ предусматривает выбор следующей группы дп длительности из кругового графика, при этом временные задержки для следующего интервала активации сдвигаются на основе временной задержки, произвольно выбираемой для источника с наименьшей случайной начальной временной задержкой следующим образом. Пусть Ό представляет длительность интервала возбуждения, 1т представляет сдвинутую временную задержку из группы дт длительности для источника с наименьшей случайной начальной временной задержкой, ζ представляет случайную начальную временную задержку для источника с наименьшей случайной начальной временной задержкой, а 11 представляет случайную начальную временную задержку для ί-го источника одновременного действия. Вычисление сдвинутых временных задержек может осуществляться при помощи следующего псевдокода:
для (ί=1; ί<η; ί++) {Т1=11+(1т-13); если Τ1>Ό то Т11 -Ό; /* возвращаемая временная задержка */возвращаемое Т1}.
Когда начинается интервал активации, источники возбуждаются в соответствии со сдвинутыми временными задержками. Процесс выбора группы длительности из кругового графика и сдвига случайных начальных временных задержек может повторяться до завершения сейсмосъемки.
На фиг. 13 показана блок-схема способа рандомизации времени возбуждения источника одновременного действия. В блоке 1301 длительность Ό интервалов активации на линии времени морской съемки делят на длину оператора, выбранного для фильтра когерентности, согласно описанию, приведенному выше со ссылкой на фиг. 5А-5В. В блоке 1302 группы длительности помещают на круговой график на основе произвольного упорядочивания с использованием псевдорандомизации, согласно описанию приведенному выше со ссылкой на фиг. 6 и 7. В блоке 1303 генерируется набор случайных начальных значений для п источников одновременного действия, согласно описанию приведенному выше со ссылкой на фиг. 3. В блоке 1304 идентифицируют источник с наименьшей случайной начальной временной задержкой и группу длительности на круговом графике, в которую попадает наименьшая случайная начальная временная задержка. В блоке 1305 источники одновременного действия возбуждаются в соответствии с временными задержками. В блоке 1306, когда сейсмосъемка завершена, исполнение способа приостанавливается; в противном случае управление переходит к блоку 1307 и исполняются операции, указанные в блоках 1307-1309, 1305 и 1306. В блоке 1307 выбирают следующую группу длительности на круговом графике. В блоке 1308 процедура сдвинуть временные задержки, описанная выше со ссылкой на фиг. 9, вызывается для выполнения сдвига временных задержек источников в соответствии с группой длительности, выбранной в блоке 1307. В блоке 1309 период ожидания между интервалами активации продолжается до его истечения согласно описанию, приведенному выше со ссылкой на фиг. 2.
На фиг. 14 показана блок-схема способа сдвига временных задержек, вызываемых в блоке 1308 блок-схемы, представленной на фиг. 13. В блоке 1401 временная задержка 1т представляет собой временную задержку, произвольно выбранную из группы длительности, выбранной на шаге 1304 или шаге 1307 на фиг. 13 для источника с наименьшей случайной начальной временной задержкой. В блоке 1402 индекс ί источника устанавливается в начальное состояние 1. В блоке 1403, если индекс ί источника меньше, чем количество п источников одновременного действия, управление передается на блок 1404, в котором вычисляется сдвинутая временная задержка Т1; в противном случае набор {Т1} сдвинутых временных задержек для интервала активации возвращается в блок 1308 способа, раскрытого на фиг. 13. В
- 8 026615 блоке 1405, если сдвинутая временная задержка Т1; больше, чем длительность И интервала активации, управление передается на блок 1406, в котором длительность И вычитается из временной задержки Т1 для получения возвращаемой временной задержки, находящейся в пределах данного интервала активации; в противном случае управление передается на блок 1407, в котором индекс ί источника увеличивается на 1, и операции, исполняемые в блоках 1403-1407, повторяются.
В раскрытом выше варианте осуществления количество групп длительности равно количеству трасс на длине оператора окна данных. Вследствие этого для каждого интервала активации группу длительности выбирают псевдослучайным образом для определения рандомизированных временных задержек в пределах каждого интервала активации, что удобно в случае области источников-приемников или области равных удалений, в которых отношение количества групп длительности к количеству трасс в операторе длины составляет 1:1. Однако следует отметить, что оператор разделения может применяться и к другим областям, таким как область общей глубинной точки. Шаг между трассами в области общей глубинной точки больше, чем в области равных удалений, что зависит от соотношения расстояния между точками активации и расстояния между различными трассами в пределах одной активации (т. е. часто 2 или 4). Вследствие этого количество трасс на длине оператора необязательно равно количеству групп длительности, используемых для реализации раскрытого выше способа псевдорандомизации. Варианты осуществления содержат выбор количества групп длительности на основе длины оператора в области, где выполняется разделение. Иными словами, в области общей глубинной точки отношение количества групп длительности к количеству трасс или длине оператора может составлять 2:1 или 4:1. Например, когда длина оператора в области общей глубинной точки (т. е. количество трасс в области общей глубинной точки) равна пяти при соотношении шага источников/трасс, равном 2:1, количество групп длительности на круговом графике равно 9. Конкретно, рассмотрим девять активаций (выстрелов) 1, 2, 3, 4, 5, 6, 7, 8 и 9 двойных источников (с двумя состояниями) в области общей глубинной точки. Для линии общей глубинной точки, принадлежащей выстрелам 2, 4, 6 и 8 с первым устойчивым состоянием, только выстрелы 1, 3, 5, 7 и 9 со вторым устойчивым состоянием вносят вклад в длину оператора, состоящую из 5 трасс. Выстрел 10 вносил бы вклад в линию общей глубинной точки со вторым устойчивым состоянием. Вследствие этого рандомизировать пришлось бы только 9 последовательных выстрелов. Отметим, что когда выстрелы каждого из двух устойчивых состояний обрабатываются по отдельности, соотношение возвращается к 1:1. Альтернативно, если длина оператора в области общей глубинной точки равна пяти при соотношении шага источников/трасс, равном 4:1, количество групп длительности на круговом графике равно 20.
На фиг. 15 показан один иллюстрирующий пример обобщенной компьютерной системы, реализующей результативный способ рандомизации времени возбуждения источников сейсмосигнала и, таким образом, представляющей собой систему обработки данных. Внутренние компоненты многих малых, средних и больших компьютерных систем, а также специализированных систем хранения данных на базе процессора могут быть описаны применительно к этой обобщенной архитектуре, хотя каждая конкретная система может содержать компоненты, подсистемы и параллельные системы с архитектурой, аналогичной этой обобщенной архитектуре. Компьютерная система содержит один или несколько центральных процессоров (ЦП) 1502-1505, одно или несколько электронных запоминающих устройств 1508, соединенных с ЦП при помощи шины 1510 ЦП/подсистемы памяти или нескольких шин, первого моста 1512, соединяющего шину 1510 ЦП/подсистемы памяти с дополнительными шинами 1514 и 1516, или другие типы средств высокоскоростного соединения, в том числе многоканальные высокоскоростные последовательные линии связи. Эти шины или последовательные линии связи, в свою очередь, соединяют ЦП и память со специализированными процессорами, такими как графический процессор 1518 и один или несколько дополнительных мостов 1520, соединенных с высокоскоростными последовательными каналами или с несколькими контроллерами 1522-1527, такими как контроллер 1527, обеспечивающий доступ к различным типам машиночитаемых носителей, таких как машиночитаемый носитель 1528, электронным дисплеям, устройствам ввода и другим подобным компонентам, подкомпонентам и вычислительным ресурсам. Электронные дисплеи, включая экран визуального отображения, громкоговорители и другие выходные интерфейсы и устройства ввода, в том числе мыши, клавиатуры, сенсорные панели и другие подобные входные интерфейсы, вместе образуют входные и выходные интерфейсы, которые позволяют компьютерной системе взаимодействовать с людьми-пользователями. Машиночитаемый носитель 1528 представляет собой устройство хранения данных, включая электронное запоминающее устройство, накопитель на оптических или магнитных дисках, И8В-накопитель, флэш-память и другие подобные устройства хранения данных. Машиночитаемый носитель 1528 может использоваться для хранения машиночитаемых команд, кодирующих описанные выше способы вычислений, а также для сохранения кодированных данных при выполнении операций сохранения, причем кодированные данные могут считываться из него при выполнении операций считывания компьютерными системами, системами хранения данных и периферийными устройствами.
После возбуждения всех источников в пределах интервала активации, каждый источник согласуют с набором записанных данных. Рассмотрим, например, источник одновременного действия с η различными источниками. Набор данных, записанных в интервале активации, согласуют с источником §ш,
- 9 026615 сдвигая время, связанное с набором записанных данных, на такую же временную задержку, присвоенную источнику §т в интервале активации. При упорядочивании в конкретной области, такой как область общей точки приема или область общей глубинной точки, волна, формируемая источником §т18, является когерентной, тогда как волны, формируемые другими источниками, некогерентны. То же самое относится к согласованию записанных данных с каждым из других источников. Записанные данные могут быть согласованы с любым из источников при условии, что волна когерентна, а волна от других источников некогерентна.
Различные описанные здесь варианты осуществления не являются исчерпывающими или ограничивающими настоящее изобретение точными формами, раскрытыми выше. Например, любое количество различных вариантов реализации способа компьютерной обработки, позволяющих осуществить способы рандомизации времен возбуждения источника одновременного действия могут проектироваться и разрабатываться с помощью различных языков программирования и компьютерных платформ путем изменения различных параметров реализации, включая управляющие структуры, переменные, структуры данных, модульную организацию и другие подобные параметры. Системы и способы рандомизации времен возбуждения источника одновременного действия могут быть реализованы почти в реальном времени при проведении морской сейсморазведки подземной формации. Термин почти в реальном времени относится к временной задержке, вызванной передачей и обработкой данных, достаточно кратковременными, чтобы обеспечить возможность своевременного использования задержек времени возбуждения, вычисленных во время сбора сейсмических данных в ходе сейсморазведки. Например, почти в реальном времени может относиться к ситуации, когда длительность формирования временных задержек и передачи соответствующих сигналов возбуждения источнику одновременного действия незначительна. В других вариантах осуществления временные задержки для каждого интервала активации и периоды ожидания сбора сейсмических данных в ходе сейсморазведки могут вычисляться заранее и храниться на машиночитаемом носителе. После начала съемки источники одновременного действия могут возбуждаться в соответствии с заданной последовательностью временных задержек и периодов ожидания. Варианты осуществления включают в себя также другие методы рандомизации при условии, что рандомизированные временные задержки для вторичных источников находятся в пределах длины одного оператора и не могут быть идентичными или находиться в пределах задаваемого пользователем минимального периода времени. Окно данных может также содержать поперечный компонент, например компонент, расположенный перпендикулярно или под углом к продольному направлению. Вследствие этого, когда окно данных содержит продольный и поперечный компоненты, оно не ограничивается длиной одномерного оператора, а имеет вместо этого площадь оператора. Раскрытые выше способы могут использоваться при любом типе морской сейсморазведки, например при разведке с использованием прямой или спиральной траектории съемки.
Следует понимать, что вышеприведенное описание раскрытых вариантов осуществления представлено для того, чтобы позволить любому специалисту воспроизвести или использовать настоящее изобретение. Различные модификации этих вариантов осуществления будут очевидны специалистам, при этом общие принципы, определенные в настоящем описании, могут применяться к другим вариантам осуществления без выхода за пределы существа и объема раскрытия настоящего изобретения. Таким образом, настоящее раскрытие изобретения не ограничивается представленными в нем вариантами осуществления, но предполагает самую широкую сферу применения, согласующуюся с принципами и новыми признаками, раскрытыми в настоящем описании.

Claims (50)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Способ активации нескольких источников, буксируемых одним или более сейсморазведочными судами, в котором генерируют набор временных задержек с начальными значениями; присваивают каждую временную задержку с начальным значением одному из источников;
    и для каждого интервала активации линии времени морской съемки случайным образом выбирают временную задержку, которая меньше длительности интервала активации;
    псевдослучайным образом сдвигают временные задержки с начальными значениями на основе указанной временной задержки для генерации сдвинутых временных задержек в интервале активации, причем каждая сдвинутая временная задержка связана с одним из указанных нескольких источников, и возбуждают несколько источников в интервале активации на основе присвоенных сдвинутых временных задержек.
  2. 2. Способ по п.1, отличающийся тем, что генерация набора временных задержек с начальными значениями включает в себя произвольную генерацию временных задержек с начальными значениями, находящихся в пределах длительности интервалов активации.
  3. 3. Способ по п.1, в котором делят длительность интервалов активации линии времени морской съемки на набор групп длитель- 10 026615 ности;
    размещают в псевдопроизвольном порядке указанные группы длительности в виде графика групп длительности;
    и для каждого интервала активации линии времени морской съемки псевдослучайным образом выбирают группу длительности из указанного набора групп длительности;
    случайным образом выбирают временную задержку из выбранной группы длительности; и псевдослучайным образом сдвигают временные задержки с начальными значениями, генерируя сдвинутые временные задержки путем присвоения временной задержки источнику с наименьшей начальной временной задержкой из указанных нескольких источников, вычисления разности между временной задержкой и наименьшей начальной временной задержкой и добавления указанной разности к начальным временным задержкам оставшихся источников с получением сдвинутых временных задержек оставшихся источников.
  4. 4. Способ по п.3, отличающийся тем, что деление длительности интервалов активации на набор групп длительности содержит определение количества групп длительности на основе отношения количества групп длительности к количеству трасс на длине оператора окна данных, равного или большего единице.
  5. 5. Способ по п.3, отличающийся тем, что псевдослучайный выбор группы длительности содержит следующие шаги: псевдослучайным образом располагают группы длительности в наборе в виде последовательности и для каждого интервала активации выбирают группу в соответствии с этой последовательностью.
  6. 6. Способ по п.3, отличающийся тем, что вычисление сдвинутых временных задержек для нескольких источников содержит следующие шаги: выбирают временную задержку в произвольно выбранной группе длительности; вычисляют разность между временной задержкой и выбранной временной задержкой и добавляют эту разность к временной задержке каждого из нескольких источников для генерации сдвинутой временной задержки для каждого источника.
  7. 7. Способ по п.6, отличающийся тем, что, если сдвинутое время задержки, связанной с источником, больше, чем длительность интервалов активации, то он содержит также шаг вычитания длительности интервала активации из сдвинутой временной задержки для генерации возвращаемой сдвинутой временной задержки, связанной с указанным источником, причем эта возвращаемая сдвинутая временная задержка находится в пределах интервала активации.
  8. 8. Способ по п.1, отличающийся тем, что сбрасывают несколько источников для следующего интервала активации и позволяют истечь периоду ожидания между интервалами активации.
  9. 9. Способ по п.1, отличающийся тем, что возбуждение нескольких источников в интервале активации на основе сдвинутых временных задержек содержит возбуждение источников в порядке от наименьшей до наибольшей сдвинутой временной задержки по отношению к началу интервала активации.
  10. 10. Способ по п.1, отличающийся тем, что каждый источник содержит пневмопушки, гидропушки, морские вибраторы или источник электромагнитного поля.
  11. 11. Компьютерная система, содержащая один или более процессоров;
    одно или более устройств хранения данных;
    и программу, хранящуюся в одном или более устройствах хранения данных, которая при исполнении одним или более процессорами обеспечивает выполнение операций способа по п.1.
  12. 12. Система по п.11, отличающаяся тем, что генерация набора временных задержек с начальными значениями включает в себя произвольную генерацию временных задержек с начальными значениями, находящихся в пределах длительности интервалов активации.
  13. 13. Система по п.11, отличающаяся тем, что указанная программа управляет компьютерной системой на предмет выполнения ею следующих операций:
    деления длительности интервалов активации линии времени морской съемки на набор групп длительности;
    размещения в псевдопроизвольном порядке указанных групп длительностей в виде графика групп длительности;
    и для каждого интервала активации линии времени морской съемки выбора псевдослучайным образом группы длительности из указанного набора групп длительности; выбора случайным образом временной задержки из выбранной группы длительности; и сдвигания псевдослучайным образом временных задержек с начальными значениями для генерирования сдвинутых временных задержек путем присвоения временной задержки источнику с наименьшей начальной временной задержкой из указанных нескольких источников, вычисления разности между временной задержкой и наименьшей начальной временной задержкой и добавления указанной разности к начальным временным задержкам оставшихся источников с получением сдвинутых временных задержек оставшихся источников.
  14. 14. Система по п.13. отличающаяся тем, что деление длительности интервалов активации на группы
    - 11 026615 длительности содержит определение количества групп длительности на основе отношения количества групп длительности к количеству трасс на длине оператора окна данных, равного или большего единице.
  15. 15. Система по п.13, отличающаяся тем, что псевдослучайный выбор группы длительности содержит расположение псевдослучайным образом групп длительности в наборе в виде последовательности и для каждого интервала активации выбор группы в соответствии с этой последовательностью.
  16. 16. Система по п.13, отличающаяся тем, что вычисление сдвинутых временных задержек для нескольких источников содержит выбор временной задержки в произвольно выбранной группе длительности; вычисление разности между временной задержкой и выбранной временной задержкой и добавление этой разности к временной задержке каждого из источников для генерации сдвинутой временной задержки для каждого источника.
  17. 17. Система по п.16, отличающаяся тем, что, если сдвинутое время задержки, связанной с источником, больше, чем длительность интервалов активации, программа предназначена также для вычитания длительности интервала активации из сдвинутой временной задержки для генерации возвращаемой сдвинутой временной задержки, связанной с указанным источником, причем эта возвращаемая сдвинутая временная задержка находится в пределах интервала активации.
  18. 18. Система по п.11, отличающаяся тем, что программа предназначена для сбрасывания нескольких источников для следующего интервала активации и позволения истечь периоду ожидания между интервалами активации.
  19. 19. Система по п.11, отличающаяся тем, что возбуждение нескольких источников в интервале активации на основе сдвинутых временных задержек содержит возбуждение источников в порядке от наименьшей до наибольшей сдвинутой временной задержки по отношению к началу интервала активации.
  20. 20. Система по п.11, отличающаяся тем, что программа предназначена для сохранения сдвинутых временных задержек в одном или более устройствах хранения данных.
  21. 21. Система по п.11, отличающаяся тем, что каждый источник содержит пневмопушки, гидропушки, морские вибраторы или источник электромагнитного поля.
  22. 22. Машиночитаемый носитель данных, содержащий машиночитаемые команды, закодированные на нем, которые обеспечивают выполнение операций способа по п.1.
  23. 23. Носитель по п.22, отличающийся тем, что машиночитаемые команды обеспечивают деление длительности интервалов активации линии времени морской съемки на набор групп длительности;
    размещение в псевдопроизвольном порядке указанных групп длительностей в виде графика групп длительности;
    и для каждого интервала активации линии времени морской съемки выбор псевдослучайным образом группы длительности из набора групп длительности;
    выбор случайным образом временной задержки из выбранной группы длительности;
    и сдвигание псевдослучайным образом временных задержек с начальными значениями для генерирования сдвинутых временных задержек путем присвоения временной задержки источнику с наименьшей начальной временной задержкой из указанных нескольких источников, вычисления разности между временной задержкой и наименьшей начальной временной задержкой и добавления указанной разности к начальным временным задержкам оставшихся источников с получением сдвинутых временных задержек оставшихся источников.
  24. 24. Носитель по п.23, отличающийся тем, что деление длительности интервалов активации на группы длительности содержит определение количества групп длительности на основе отношения количества групп длительности к количеству трасс на длине оператора окна данных, равного или большего единице.
  25. 25. Носитель по п.23, отличающийся тем, что псевдослучайный выбор группы длительности содержит расположение псевдослучайным образом групп длительности в наборе в виде последовательности и для каждого интервала активации выбор группы в соответствии с этой последовательностью.
  26. 26. Носитель по п.23, отличающийся тем, что вычисление сдвинутых временных задержек для нескольких источников содержит выбор временной задержки в произвольно выбранной группе длительности; вычисление разности между временной задержкой и выбранной временной задержкой и добавление этой разности к временной задержке каждого источника для генерации сдвинутой временной задержки для каждого источника.
  27. 27. Носитель по п.26, отличающийся тем, что, если сдвинутое время задержки, связанной с источником, больше, чем длительность интервалов активации, машиночитаемые команды обеспечивают также вычитание длительности интервала активации из сдвинутой временной задержки для генерации возвращаемой сдвинутой временной задержки, связанной с указанным источником, причем эта возвращаемая сдвинутая временная задержка находится в пределах интервала активации.
  28. 28. Носитель по п.23, отличающийся тем, что машиночитаемые команды обеспечивают сбрасывание двух или более источников для следующего интервала активации и позволение истечь периоду ожидания между интервалами активации.
  29. 29. Носитель по п.23, отличающийся тем, что возбуждение источников в интервале активации на
    - 12 026615 основе сдвинутых временных задержек содержит возбуждение источников в порядке от наименьшей до наибольшей сдвинутой временной задержки по отношению к началу интервала активации.
  30. 30. Носитель по п.22, отличающийся тем, что каждый источник содержит пневмопушки, гидропушки, морские вибраторы или источник электромагнитного поля.
  31. 31. Морская сейсморазведочная система, содержащая несколько источников и компьютерную систему по п.11.
  32. 32. Система по п.31, отличающаяся тем, что генерация набора временных задержек с начальными значениями включает в себя произвольную генерацию временных задержек с начальными значениями, находящихся в пределах длительности интервалов активации.
  33. 33. Система по п.31, отличающаяся тем, что указанная программа управляет компьютерной системой на предмет выполнения ею следующих операций:
    деления длительности интервалов активации линии времени морской съемки на набор групп длительности;
    размещения в псевдопроизвольном порядке указанных групп длительностей в виде графика групп длительности;и для каждого интервала активации линии времени морской съемки выбора псевдослучайным образом группы длительности из набора групп длительности; выбора случайным образом временной задержки из выбранной группы длительности; и сдвигания псевдослучайным образом временных задержек с начальными значениями для генерирования сдвинутых временных задержек путем присвоения временной задержки источнику с наименьшей начальной временной задержкой из указанных нескольких источников, вычисления разности между временной задержкой и наименьшей начальной временной задержкой и добавления указанной разности к начальным временным задержкам оставшихся источников с получением сдвинутых временных задержек оставшихся источников.
  34. 34. Система по п.31. отличающаяся тем, что деление длительности интервалов активации на набор групп длительности содержит определение количества групп длительности на основе отношения количества групп длительности к количеству трасс на длине оператора окна данных, равного или большего единице.
  35. 35. Система по п.31, отличающаяся тем, что псевдослучайный выбор группы длительности содержит расположение псевдослучайным образом групп длительности в наборе в виде последовательности и для каждого интервала активации выбор группы в соответствии с этой последовательностью.
  36. 36. Система по п.31, отличающаяся тем, что вычисление сдвинутых временных задержек для нескольких источников содержит выбор временной задержки в произвольно выбранной группе длительности; вычисление разности между временной задержкой и выбранной временной задержкой и добавление этой разности к временной задержке каждого из нескольких источников для генерации сдвинутой временной задержки для каждого источника.
  37. 37. Система по п.36, отличающаяся тем, что, если сдвинутое время задержки, связанной с источником, больше, чем длительность интервалов активации, программа предназначена для вычитания длительности интервала активации из сдвинутой временной задержки для генерации возвращаемой сдвинутой временной задержки, связанной с указанным источником, причем эта возвращаемая сдвинутая временная задержка находится в пределах интервала активации.
  38. 38. Система по п.36, отличающаяся тем, что программа предназначена для сбрасывания источников для следующего интервала активации и позволения истечь периоду ожидания между интервалами активации.
  39. 39. Система по п.31, отличающаяся тем, что возбуждение источников в интервале активации на основе сдвинутых временных задержек содержит возбуждение источников в порядке от наименьшей до наибольшей сдвинутой временной задержки по отношению к началу интервала активации.
  40. 40. Система по п.31, отличающаяся тем, что программа предназначена для сохранения сдвинутых временных задержек в одном или более устройствах хранения данных.
  41. 41. Система по п.31, отличающаяся тем, что каждый источник содержит пневмопушки, гидропушки, морские вибраторы или источник электромагнитного поля.
  42. 42. Способ генерации моментов активации источников, используемый для выполнения способа по п.1, содержащий следующие этапы:
    генерацию случайной последовательности временных задержек с начальными значениями; присвоение одной временной задержки с начальным значением из указанных временных задержек с начальными значениями каждому источнику из нескольких источников; и для каждого интервала активации линии времени морской съемки случайный выбор временной задержки, которая меньше длительности интервала активации; случайную генерацию временного сдвига на основе указанной временной задержки; и применение этого временного сдвига к каждой временной задержке с начальным значением.
  43. 43. Способ по п.42, в котором этап генерации случайной последовательности временных задержек с начальными значениями включает в себя генерацию временных задержек с начальными значениями в
    - 13 026615 пределах длительности интервалов активации.
  44. 44. Способ по п.42, в котором этапы случайной генерации временного сдвига и применения этого временного сдвига к каждой временной задержке с начальным значением содержат деление длительности интервалов активации на группы длительности;
    произвольное упорядочивание указанных групп длительности по круговому графику;
    при этом для первого интервала активации из указанных интервалов активации линии времени морской съемки присвоение указанным источникам случайной последовательности временных задержек с начальными значениями в качестве времени возбуждения источников и идентификацию группы длительности в круговом графике, содержащей наименьшую временную задержку из последовательности временных задержек;
    а для любого другого интервала активации после первого случайный выбор временной задержки из соседней группы длительности на круговом графике, расчет временного сдвига как разность между наименьшей временной задержкой в последовательности временных задержек и произвольно выбранной временной задержкой, добавление этого временного сдвига к временным задержкам с начальными значениями, присвоенным каждому источнику, и присвоение указанным источникам сдвинутых временных задержек в качестве времени возбуждения источников.
  45. 45. Способ по п.44, дополнительно содержащий, если сдвинутая временная задержка, связанная с источником, больше, чем длительность интервала активации, шаг вычитания длительности интервала активации из сдвинутой временной задержки для генерации возвращаемой сдвинутой временной задержки, связанной с указанным источником, в пределах интервала активации.
  46. 46. Способ по п.44, в котором деление длительности интервалов активации на группы длительности содержит определение количества групп длительности на основе отношения количества групп длительности к количеству трасс на длине оператора окна данных, равного или большего единице.
  47. 47. Способ по п.42, дополнительно содержащий этапы, на которых сбрасывают указанные несколько источников для следующего интервала активации и позволяют истечь периоду ожидания между интервалами активации.
  48. 48. Способ по п.42, дополнительно содержащий возбуждение источников в порядке от наименьшей до наибольшей сдвинутой временной задержки по отношению к началу интервала активации.
  49. 49. Способ по п.42, дополнительно содержащий возбуждение источников в каждом интервале активации с присвоенными временными задержками, сдвинутыми на временной сдвиг, сгенерированный для этого интервала активации.
  50. 50. Способ по п.42, дополнительно содержащий сохранение сдвинутых временных задержек для каждого интервала активации в устройстве хранения данных.
EA201490431A 2013-03-15 2014-03-05 Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке EA026615B1 (ru)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US13/835,711 US10359528B2 (en) 2013-03-15 2013-03-15 Systems and methods for randomizing firing times of simultaneous sources in marine surveys

Publications (2)

Publication Number Publication Date
EA201490431A1 EA201490431A1 (ru) 2014-09-30
EA026615B1 true EA026615B1 (ru) 2017-04-28

Family

ID=50238270

Family Applications (2)

Application Number Title Priority Date Filing Date
EA201490431A EA026615B1 (ru) 2013-03-15 2014-03-05 Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке
EA201790023A EA033617B1 (ru) 2013-03-15 2014-03-05 Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке

Family Applications After (1)

Application Number Title Priority Date Filing Date
EA201790023A EA033617B1 (ru) 2013-03-15 2014-03-05 Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке

Country Status (9)

Country Link
US (2) US10359528B2 (ru)
EP (1) EP2778721B1 (ru)
CN (1) CN104049276B (ru)
AU (3) AU2014201584B2 (ru)
BR (1) BR102014005983B1 (ru)
CA (2) CA2843979C (ru)
EA (2) EA026615B1 (ru)
MX (1) MX355338B (ru)
MY (1) MY184805A (ru)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9103943B2 (en) * 2011-11-28 2015-08-11 Fugro-Geoteam As Acquisition and processing of multi-source broadband marine seismic data
US9188693B2 (en) * 2012-03-26 2015-11-17 Apache Corporation Method for acquiring marine seismic data
US10551515B2 (en) * 2013-07-30 2020-02-04 Westerngeco, L.L.C. Determining an interval between activations of at least one survey source
US11391857B2 (en) * 2013-12-30 2022-07-19 Pgs Geophysical As Methods and systems for attenuating residual acoustic energy in seismic data
US10670757B2 (en) * 2014-02-26 2020-06-02 Pgs Geophysical As Methods and systems for quantifying coherency and constraining coherency-based separation in simultaneous shooting acquisition
US20170363760A1 (en) * 2014-12-05 2017-12-21 Cgg Services Sas Multi-vessel seismic data acquisition system
US10878142B2 (en) * 2015-10-26 2020-12-29 Pgs Geophysical As Bipole source modeling
GB2553744B (en) * 2016-04-29 2018-09-05 Advanced Risc Mach Ltd Graphics processing systems
US10571589B2 (en) * 2016-08-17 2020-02-25 Pgs Geophysical As Constraint of dithering of source actuations
US10871588B2 (en) 2016-12-14 2020-12-22 Pgs Geophysical As Seismic surveys with increased shot point intervals for far offsets
US11280924B2 (en) 2017-05-17 2022-03-22 Saudi Arabian Oil Company Suppressing noises in seismic data
US9910175B1 (en) * 2017-07-12 2018-03-06 Edward Majzlik Marine seismic survey system for generating and collecting data and forming a seismic image
US10649108B2 (en) 2017-10-19 2020-05-12 Cgg Services Sas System and method for generating dithering sequences for seismic exploration
US11175425B2 (en) * 2017-12-18 2021-11-16 Pgs Geophysical As Survey design for data acquisition using marine non-impulsive sources
US10310124B1 (en) * 2018-02-28 2019-06-04 Lawrence Scott Floating vessel based system for generating a multidimensional seismic data set for a target area
AU2019275488B2 (en) 2018-05-23 2022-03-03 Blue Ocean Seismic Services Limited An autonomous data acquisition system and method
US11009618B2 (en) 2018-09-04 2021-05-18 Sercel System and method for generating dithering sequences with minimum value for seismic exploration
CN109240317B (zh) * 2018-11-21 2021-06-04 哈尔滨工程大学 考虑螺旋桨故障的海底地震检波飞行节点有限时间构型包含控制方法
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2390902A (en) * 2002-07-17 2004-01-21 Pgs Americas Inc Firing sequence of multiple sources in a marine seismic survey
WO2005062075A1 (en) * 2003-12-15 2005-07-07 Chevron, U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US20110158042A1 (en) * 2009-12-30 2011-06-30 Nicolae Moldoveanu Randomization of Data Acquisition in Marine Seismic and Electromagnetic Acquisition
US20110188343A1 (en) * 2010-02-04 2011-08-04 Guigne Jacques Y Random Transmitter Placement Method For Stationary Seismic Imaging

Family Cites Families (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3697938A (en) 1969-11-07 1972-10-10 Seismic Computing Corp Seismic prospecting with random injected signal
US4739858A (en) * 1987-03-02 1988-04-26 Western Atlas International, Inc. Spectrally-shaped air gun arrays
US4953657A (en) 1987-11-30 1990-09-04 Halliburton Geophysical Services, Inc. Time delay source coding
SU1543357A1 (ru) 1988-01-07 1990-02-15 Киевское геофизическое отделение Украинского научно-исследовательского геологоразведочного института Способ пространственной сейсморазведки по методу многократных перекрытий
US4908801A (en) * 1989-05-30 1990-03-13 Teledyne Exploration Real-time simulation of the far-field signature of a seismic sound source array
US5198979A (en) * 1989-09-26 1993-03-30 Shell Oil Company Seismic migration of multiprocessor computer
US5596349A (en) * 1992-09-30 1997-01-21 Sanyo Electric Co., Inc. Image information processor
US5924049A (en) 1995-04-18 1999-07-13 Western Atlas International, Inc. Methods for acquiring and processing seismic data
US5764516A (en) * 1995-12-29 1998-06-09 Atlantic Richfield Company Method and system for surface-consistent phase and time lag correction of seismic data
US5761152A (en) * 1996-10-29 1998-06-02 Pgs Exploration (Us), Inc. Method and system for increasing fold to streamer length ratio
US6574567B2 (en) 2001-01-23 2003-06-03 Pgs Americas, Inc. Weighted slant stack for attenuating seismic noise
CA2348355A1 (en) * 2001-05-24 2002-11-24 Cloakware Corporation General scheme of using encodings in computations
US6545944B2 (en) 2001-05-30 2003-04-08 Westerngeco L.L.C. Method for acquiring and processing of data from two or more simultaneously fired sources
US6751559B2 (en) * 2002-09-10 2004-06-15 Pgs Exploration (Uk) Limited Method for suppressing noise from seismic signals by source position determination
FR2890177B1 (fr) 2005-09-01 2007-11-02 Cie Generale De Geophysique Sa Procede et systeme d'acquisition de donnees sismiques
US7935545B2 (en) * 2007-03-30 2011-05-03 Tokyo Electron Limited Method and apparatus for performing a site-dependent dual patterning procedure
US8559265B2 (en) * 2007-05-17 2013-10-15 Westerngeco L.L.C. Methods for efficiently acquiring wide-azimuth towed streamer seismic data
US7869303B2 (en) * 2007-08-14 2011-01-11 Pgs Geophysical As Method for noise suppression in seismic signals using spatial transforms
US8345510B2 (en) 2008-06-02 2013-01-01 Pgs Geophysical As Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s)
WO2010019957A1 (en) * 2008-08-15 2010-02-18 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US8395966B2 (en) 2009-04-24 2013-03-12 Westerngeco L.L.C. Separating seismic signals produced by interfering seismic sources
US8553496B2 (en) 2010-02-09 2013-10-08 Ion Geophysical Corporation Seismic source separation
US8579545B2 (en) * 2010-03-02 2013-11-12 Fairfield Industries Incorporated Apparatus and methods for an ocean bottom seismic sensor deployment vehicle
US9395460B2 (en) * 2010-07-19 2016-07-19 Conocophillips Company Continuous composite relatively adjusted pulse
US8339896B2 (en) 2010-08-16 2012-12-25 Pgs Geophysical As Method for separating seismic sources in marine seismic surveys
EP2649472B1 (en) * 2010-12-09 2016-05-11 BP Corporation North America Inc. Seismic acquisition method and system
US9513389B2 (en) * 2011-05-13 2016-12-06 Saudi Arabian Oil Company Coupled time-distance dependent swept frequency source acquisition design and data de-noising
US9075162B2 (en) 2011-11-10 2015-07-07 Pgs Geophysical As Method and system for separating seismic sources in marine simultaneous shooting acquisition
US9261619B2 (en) * 2012-01-03 2016-02-16 Cggveritas Services Sa Method and device for marine seismic acquisition
US8724428B1 (en) * 2012-11-15 2014-05-13 Cggveritas Services Sa Process for separating data recorded during a continuous data acquisition seismic survey
US9360575B2 (en) * 2013-01-11 2016-06-07 Fairfield Industries Incorporated Simultaneous shooting nodal acquisition seismic survey methods
US10520623B2 (en) 2013-05-31 2019-12-31 Westerngeco L.L.C. Methods and systems for marine survey acquisition

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2390902A (en) * 2002-07-17 2004-01-21 Pgs Americas Inc Firing sequence of multiple sources in a marine seismic survey
WO2005062075A1 (en) * 2003-12-15 2005-07-07 Chevron, U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US20110158042A1 (en) * 2009-12-30 2011-06-30 Nicolae Moldoveanu Randomization of Data Acquisition in Marine Seismic and Electromagnetic Acquisition
US20110188343A1 (en) * 2010-02-04 2011-08-04 Guigne Jacques Y Random Transmitter Placement Method For Stationary Seismic Imaging

Also Published As

Publication number Publication date
EP2778721A3 (en) 2016-02-10
CA2843979A1 (en) 2014-09-15
BR102014005983A2 (pt) 2015-10-06
AU2014201584A1 (en) 2014-10-02
CA3092480A1 (en) 2014-09-15
CA2843979C (en) 2022-11-01
EA201490431A1 (ru) 2014-09-30
MX355338B (es) 2018-04-13
AU2018202056B2 (en) 2019-11-07
US10359528B2 (en) 2019-07-23
EP2778721B1 (en) 2023-01-25
US11422278B2 (en) 2022-08-23
AU2020200770B2 (en) 2021-08-05
CN104049276B (zh) 2019-08-20
AU2018202056A1 (en) 2018-04-12
BR102014005983B1 (pt) 2021-05-25
MX2014003199A (es) 2014-11-03
US20190086570A1 (en) 2019-03-21
EA201790023A3 (ru) 2017-10-31
EA201790023A2 (ru) 2017-09-29
EP2778721A2 (en) 2014-09-17
EA033617B1 (ru) 2019-11-11
US20140269169A1 (en) 2014-09-18
AU2020200770A1 (en) 2020-02-20
MY184805A (en) 2021-04-23
AU2014201584B2 (en) 2017-12-14
CN104049276A (zh) 2014-09-17

Similar Documents

Publication Publication Date Title
EA026615B1 (ru) Система и способ рандомизации времени возбуждения источников сейсмосигнала в морской сейсморазведке
US11099287B2 (en) Seismic data processing
AU2018220021B2 (en) Systems and methods for removing acquisition related effects from seismic data
AU2015268879B2 (en) System and method of building complex earth models
EP3710867B1 (en) Noise attenuation of multiple source seismic data
US20090326895A1 (en) Technique and system for seismic source separation
EA037851B1 (ru) Способ сейсмической съемки
AU2014218351A1 (en) Attenuating noise by shot repetition
US11391857B2 (en) Methods and systems for attenuating residual acoustic energy in seismic data
EP3948360A1 (en) Low-frequency seismic survey design
AU2015205963A1 (en) Methods and systems to operate marine sources to avoid air in water effects
GB2520124A (en) Methods and systems for attenuating noise in seismic data
CA3106943A1 (en) Fast plane-wave reverse time migration
GB2530399A (en) Methods and systems to operate marine sources to avoid air in water effects

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 TJ TM

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