EA004857B1 - Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности - Google Patents

Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности Download PDF

Info

Publication number
EA004857B1
EA004857B1 EA200300769A EA200300769A EA004857B1 EA 004857 B1 EA004857 B1 EA 004857B1 EA 200300769 A EA200300769 A EA 200300769A EA 200300769 A EA200300769 A EA 200300769A EA 004857 B1 EA004857 B1 EA 004857B1
Authority
EA
Eurasian Patent Office
Prior art keywords
traces
seismic
sample
slope
summation
Prior art date
Application number
EA200300769A
Other languages
English (en)
Other versions
EA200300769A1 (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 EA200300769A1 publication Critical patent/EA200300769A1/ru
Publication of EA004857B1 publication Critical patent/EA004857B1/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times

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)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

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

Description

Настоящее изобретение относится к способу обработки сейсмических данных и в частности - к способу формирования панели подобия. Оно также относится к устройству для обработки сейсмических данных.
Предшествующий уровень техники
После сбора сейсмических данных их обрабатывают для получения информации о геологической структуре внутренней области грунта. Для обработки сейсмических данных необходимо знать или оценить, каким образом изменяется скорость сейсмических сигналов в грунте (этот этап также известен как определение «поля скоростей»). Точность информации, полученной относительно внутренней области грунта, будет зависеть от точности поля скоростей, использованного при обработке сейсмических данных.
Этап определения или оценки скорости сейсмической энергии по сейсмической трассе обычно известен как «выделение скорости». По существу он заключается в идентификации конкретного признака в сейсмической трассе и в приписывании этому признаку конкретной траектории волны сейсмической энергии. После этого скорость сейсмической энергии определяют на основании длины траектории волны сейсмической энергии и времени пробега волны сейсмической энергии, имеющей эту траекторию.
В настоящее время выделение скорости осуществляют вручную, что требует чрезвычайно интенсивного труда, а затраты, необходимые для обработки сейсмических данных весьма значительны.
Предлагалось много способов автоматического выделения скорости с использованием компьютера. Однако надежность таких способов низкая, и требуется выполняемый вручную контроль качества результатов автоматического выделения скорости. Контроль качества осуществляется столь же долго, как и выделение скорости вручную.
Поэтому существует необходимость в способе выделения скорости, который требовал бы меньших затрат труда, чем современные способы, но при этом был бы более надежным, чем современные способы автоматического выделения скорости.
Во многих современных способах используют «панели подобия». Панели подобия широко используются в сейсмической разведке в качестве средства для определения полей скоростей, предназначенных для дальнейшей обработки сейсмических данных.
Панель подобия состоит из карты изолиний, в осях времени и скорости в плоскости, а по вертикальной оси, то есть по оси, перпендикулярной к плоскости время-скорость, отсчитывается подобие. Цель выделения скорости заключается в нахождении функции в пределах пространства время-скорость, которая проходит через самые высшие точки карты в изолиниях, наряду с тем, что удовлетворяет ограничивающим условиям, таким, как непрерывное возрастание времени и ограничения, касающиеся инверсий скорости. В известном способе автоматического выделения скорости, основанном на процедуре итерации, сделана попытка обнаружить наивысшие точки карты в изолиниях путем возмущения «порождающей» функции так, чтобы склоны максимумов переместились на карте кверху, а затем максимизировать интеграл при функции.
Одна из проблем, связанных с этим известным способом, заключается в том, что этим способом можно находить локальный максимум, эффективно обнаруживая предгорье максимума, а не сам максимум. Такие «предгорья» часто вносятся в карту в изолиниях шумом в сейсмических данных. Как только способ автоматического выделения скорости ошибочно сориентируется на предгорье, а не на главный максимум, он будет оставаться верным предгорью и не будет обнаруживать главный максимум. Этим способом автоматического выделения не будет выделяться правильная скорость, и следовательно, необходимо вручную контролировать качество результатов.
Была сделана попытка повышения точности известных способов автоматического выделения скорости путем снижения шума в сейсмических данных. При этом уменьшается высота предгорий и следовательно, уменьшается вероятность того, что способ автоматического выделения скорости будет обнаруживать предгорье, а не главный максимум. Использование данных из других точек, находящихся в непосредственной близости от точек анализа, представляет собой способ увеличения объема данных и следовательно, повышения отношения сигнала к шуму. Проблема, связанная с этим способом, заключается в том, что, если обрабатываемые сейсмические данные были получены от наклонной отражающей поверхности, то есть от поверхности, которая не является горизонтальной, то максимумы подобия от соседних точек не будут наблюдаться в одно и то же время. Это означает, что при использовании данных более чем от одной точки анализа происходит размывание максимумов подобия, и это размывание будет сводить на нет любое повышение точности, обусловленное повышенным отношением сигнала к шуму. В результате этот способ используют нечасто. Кроме того, если используют этот способ, то только в очень ограниченных пределах, в которых данные от весьма небольшого количества соседних точек могут быть просуммированы до появления значительного размывания максимумов подобия.
Краткое изложение существа изобретения
В основу настоящего изобретения поставлена задача создания способа обработки сейс3 мических данных, сформированных при отражении сейсмической энергии от отражающей поверхности.
Способ, согласно изобретению включает в себя этапы: выбор η-ой выборки трасс и по меньшей мере еще одной ближайшей выборки трасс; и формирование панели подобия, связанной с местом, соответствующим η-ой выборке из выбранной выборки, с учетом наклона отражающей поверхности.
В принципе, в указанном способе можно использовать любую выборку из сейсмических трасс, содержащих когерентную энергию, которая пригодна для формирования подобия. Примеры удовлетворяющих требованиям выборок трасс включают в себя выборки трасс средней точки, выборки трасс общей точки изображения или выборки трасс общей глубинной точки.
Согласно изобретению во время формирования панели подобия учитывают наклон отражающей поверхности поблизости от места анализа. Это предотвращает или по меньшей мере уменьшает размывание максимумов подобия при присоединении данных от соседних выборок трасс для повышения отношения сигнала к шуму. Уменьшение или исключение размывания обеспечивает возможность объединения данных большего количества выборок трасс. Обычно суммируют данные более чем 20 соседних выборок трасс без возникновения значительного размывания. Это, в свою очередь, обеспечивает повышение отношения сигнала к шуму в сейсмических данных с уменьшением вероятности того, что алгоритм автоматического выделения скорости будет обнаруживать дополнительный максимум, а не основной максимум.
Поэтому достигаемый технический эффект от использования настоящего изобретения заключается в повышении точности существующих способов автоматического выделения скорости и тем самым уменьшении необходимости в ручном контроле результатов автоматического выделения скорости.
Настоящее изобретение можно использовать не только при автоматическом выделении скорости, но также и при ручном выделении, поскольку максимумы подобия более четко воспринимаются глазом человека.
Преимущества, обеспечиваемые настоящим изобретением, можно пояснить следующим образом. При выделении скорости используют не только панель подобия, но и также панели результатов суммирования с последовательностью многих скоростных функций. Набор панелей результатов суммирования с последовательностью многих скоростных функций состоит из локализованных просуммированных изображений, полученных из выборок трасс в месте, для которого осуществляют анализ скорости. Панели результатов суммирования с последовательностью многих скоростных функ ций формируют путем многократного суммирования трасс в выборках, используя иную скоростную функцию для каждого процесса суммирования. Диапазон скоростей, охватываемых скоростными функциями, используемыми для формирования панелей результатов суммирования с последовательностью многих скоростных функций, должен быть рассчитан на полное включение диапазона погрешностей вблизи начальной скоростной функции, то есть он должен включать диапазон скоростей, в пределах которого предполагается нахождение реальных скоростей. Поэтому характеристику суммирования можно оценить количественно для диапазона скоростных функций. При скоростной функции суммирования, которая близка к реальной сейсмической скоростной функции, вступления волн в просуммированных сейсмических данных будут резко очерченными и хорошо определенными. По мере того, как скоростная функция суммирования отдаляется от истинной скоростной функции сейсмической энергии, вступление волн в просуммированных данных становится более диффузным и в конечном счете исчезает. Чтобы повысить вероятность правильного выделения скорости, оператор, выделяющий скорость, в дополнение к панели подобия может использовать панели результатов суммирования с последовательностью многих скоростных функций.
В известном способе автоматического выделения скорости нельзя использовать панели результатов суммирования с последовательностью многих скоростных функций, а можно использовать только панель подобия. Настоящее изобретение предлагает способ включения информации панелей результатов суммирования с последовательностью многих скоростных функций в панель подобия. Тем самым обеспечивается способ автоматического выделения скорости на основе панели подобия. В предпочтительном варианте осуществления настоящего изобретения способ осуществляется путем идентификации вступлений на панелях результатов суммирования с последовательностью многих скоростных функций при использовании поиска наклона и последующего суммирования с целью формирования панели подобия вдоль определенного наклона, а не при постоянном времени, как известно в технике.
Изобретение не ограничено обработкой сейсмических данных, содержащих волны, отраженные только от одной отражающей поверхности, а может быть применено к сейсмическим данным, содержащим отраженные волны, связанные с более чем одной отражающей поверхностью.
Наклон отражающей поверхности в районе сейсмической разведки может быть уже известен из предсуществующих данных. Например, наклон может быть известен на основании более ранних сейсмических разведок, проведенных в этом месте, или наклон можно получить из предсуществующих сейсмических данных, собранных в месте сейсмической разведки. В качестве альтернативы наклон можно оценить из других данных, например из геологических данных. Если наклон известен или может быть определен одним из этих способов, то информация о наклоне может быть использована в настоящем изобретении. Однако наклон отражающей поверхности в месте сейсмической разведки известен не всегда, вследствие чего в предпочтительный вариант осуществления изобретения включено определение наклона отражающей поверхности на основании сейсмических данных.
В основу настоящего изобретения поставлена также задача создания устройства для обработки сейсмических данных, сформированных при отражении сейсмической энергии от отражающей поверхности.
Устройство согласно изобретению содержит средство для выбора η-ой выборки трасс из принятой сейсмической энергии и по меньшей мере другой ближайшей выборки трасс и средство для формирования панели подобия, связанной с местом, соответствующим η-ой выборке из выбранных выборок, с учетом наклона отражающей поверхности.
В предпочтительном варианте осуществления устройство содержит программируемый процессор данных.
В изобретении также предусмотрено запоминающее средство, содержащее программу для процессора данных из устройства, определенного выше.
Краткое описание чертежей
В дальнейшем изобретение поясняется описанием предпочтительного варианта воплощения со ссылками на сопровождающие чертежи, на которых:
фиг. 1 изображает схему отражения волн от пластов в разрезе, иллюстрирующую наземную сейсмическую разведку;
фиг. 2 - схему формирования панели подобия, согласно изобретению;
фиг. 3 - схему формирования панели подобия, согласно другому варианту осуществления настоящего изобретения;
фиг. 4 - схему определения угла наклона из панели результата суммирования с последовательностью многих скоростных функций, согласно изобретению;
фиг. 5 - блок-схему устройства для обработки сейсмических данных, согласно изобретению.
Описание предпочтительного примера воплощения изобретения, в котором панель подобия формируют из выборок трасс общей средней точки.
На фиг. 1 показана схема сбора сейсмических данных при сейсмической разведке. Сейсмические данные собирают, используя группу сейсмических источников и сейсмических приемников. В случае сейсмической разведки, осуществляемой на суше, данные можно собирать, используя, например, взрывные заряды в качестве сейсмических источников и геофоны в качестве сейсмических приемников. Схема расположения сейсмических источников и приемников известна в области наземной сейсмической разведки, поскольку все источники и приемники находятся на поверхности грунта.
Первый, второй и третий сейсмические источники 1, 2 и 3 согласованно действуют с первым, вторым и третьим приемниками 4, 5 и 6, соответственно. Источники и приемники расположены относительно средней точки 7. Для упрощения предполагается, что формация или порода 8 ниже источников и приемников является изотропной и имеет локальные первую и вторую горизонтальные отражающие поверхности 9 и 10. Сейсмическая энергия, формируемая при возбуждении первого источника 1, отражается от локальных отражающих поверхностей 9 и 10 и принимается каждым из приемников 4, 5 и 6.
Для упрощения будут рассматриваться только пути прохождения сейсмической энергии, включающие в себя отражения, которые происходят непосредственно ниже средней точки 7. Например, будет рассматриваться энергия, принимаемая на первом приемнике 4 только в результате возбуждения первого источника 1, которая отражается ниже средней точки 7, энергия, принимаемая на втором приемнике 5 в результате возбуждения второго сейсмического источника 2, и энергия, принимаемая на третьем приемнике 6 в результате возбуждения третьего сейсмического источника 3. Точка 7 представляет собой общую среднюю точку (ОСТ) для всех путей прохождения энергии (фиг. 1). Для других путей прохождения сейсмической энергии (не показаны) точка 7 не будет являться средней точкой. При типовой сейсморазведочной расстановке будет иметься большое число возможных путей прохождения сейсмической энергии, и эти пути прохождения будут иметь большое количество соответствующих общих средних точек.
После сбора необработанных сейсмических данных отраженные сигналы (известные как трассы), принятые каждым из приемников в результате возбуждения источников сейсмической энергии, обрабатывают для получения изображения внутренней области грунта. Во время обработки сейсмических данных, полученных при типовой сейсмической разведке, трассы сначала сортируют так, чтобы трассы, имеющие одну и ту же общую среднюю точку, были сгруппированы друг с другом. Группа трасс, относящихся к общей средней точке, известна как «выборка трасс общей средней точки». Она
Ί позволяет в ряде мест исследовать геологию под линией источников и приемников.
Три выборки 11 (фиг. 2) трасс общей средней точки показаны справа, при этом предполагается, что существуют выборки трасс, соответствующие общим средним точкам с номерами (п-1), п, (п+1). Каждая выборка состоит из набора сейсмических трасс, которые характеризуют амплитуду сейсмической энергии, принятой на приемнике, в зависимости от времени. Каждая трасса в выборке трасс общей средней точки соответствует одному приемнику в группе приемников сейсморазведочной расстановки. Трассы в выборках трасс общей средней точки являются вертикальными, то есть вертикальная ось характеризует время, а горизонтальная ось характеризует амплитуду сейсмической энергии, принятой на приемнике.
Ниже описан известный способ формирования панели подобия из выборок 11 трасс общей средней точки. Панель подобия с номером п, то есть панель подобия для п-ой общей средней точки, обычно формируют из п-ой выборки трасс общей средней точки. Данные из соседних выборок трасс общей средней точки также используют при формировании п-ой панели подобия, что пояснено выше, дополнительные данные можно брать только из нескольких дополнительных выборок трасс общей средней точки. Дополнительные данные берут из (п-1)-ой и (п+1)-ой выборок трасс общей средней точки. Панель подобия формируют, используя соотношение
где: Лук - представляет собой индивидуальные трассы в выборках трасс общей средней точки, _) - идентифицирует индивидуальные трассы в выборке трасс общей средней точки, суммирование по _) осуществляется над трассами в выборке трасс общей средней точки от _)=1 до _)=М, при этом в каждой выборке 11 трасс общей средней точки имеются М трасс.
Индексом к идентифицируется выборка трасс общей средней точки. В этом примере из уровня техники выборки трасс с номерами (п-1), п и (п+1) общей средней точки используют для формирования п-ой панели подобия, так что суммирование по к осуществляют от п-1 до п+1. В этом способе суммирование по к осуществляют при постоянном значении времени, и это эквивалентно предположению о том, что трассы были получены при отражении от горизонтальной отражающей поверхности, то есть угол наклона отражающей поверхности равен нулю.
Конечное суммирование осуществляют по индексу 1. Эта операция представляет собой суммирование отсчетов в пределах окна в сейс мических данных. Сейсмический сигнал является непрерывным во времени, но в процессе регистрации сигнала в цифровой системе в последовательные дискретные моменты времени из сигнала берутся отсчеты (обычно на регулярном временном интервале, известном как интервал отсчетов). Поэтому сейсмическая трасса Лук описывается рядом чисел, которые характеризуют амплитуду сейсмического сигнала в момент каждого отсчета. Отсчеты с 1=1 до N характеризуют временное «окно» в каждой сейсмической трассе.
Как отмечено выше, поскольку в известном способе суммирование по индексу к осуществляют в предположении, что отражающая поверхность имеет нулевой угол наклона, то при суммировании обычно можно использовать только несколько выборок трасс общей средней точки до того момента, как максимумы в панели подобия станут размытыми. В заявленном изобретении при формировании панели подобия учитывают угол наклона отражающей поверхности.
В заявленном способе для формирования панели подобия используют две или более выборок трасс общей средней точки. Панель подобия для места, соответствующего п-ой выборке трасс, образуют из п-ой выборки трасс и по меньшей мере одной ближайшей выборки трасс. В предпочтительном варианте осуществления используют р соседних выборок трасс по каждую сторону от п-ой выборки трасс общей средней точки. Поскольку при формировании панелей подобия в соответствии с настоящим изобретением учитывают угол наклона, то для формирования панели подобия можно использовать большее число выборок трасс общей средней точки, чем в известном способе. При формировании панели подобия можно использовать много выборок трасс общей средней точки, например можно использовать 10 соседних выборок трасс по каждую сторону от п-ой выборки трасс общей средней точки.
Угол наклона отражающей поверхности можно оценить на основании сейсмических данных, например путем анализа вступлений отражений вблизи п-го местоположения общей средней точки, а оценку угла наклона можно использовать при формировании панели подобия. В качестве альтернативы, если имеется предшествующая информация относительно угла наклона в месте разведки или угол наклона можно определить из предшествующих сейсмических данных, то такую предшествующую информацию можно использовать при формировании панели подобия.
На фиг. 3 справа показаны выбранные выборки 11 трасс общей средней точки, и при этом предполагается, что они представляют собой выборки трасс, соответствующие общим средним точкам с номерами (п-р) , ...(п-1), п, (п+1) ,...
(п+р) . Если по каждую сторону от п-ой выборки трасс выбрать 10 соседних выборок трасс, то р=10.
Выборки 11 трасс общей средней точки показаны после того, как трассы в каждой выборке скорректированы с учетом переменного горизонтального расстояния (известного как «удаление») между сейсмическим источником и приемником, использованными для получения этой трассы. Трассы в выборке общей средней точки получены при различных удалениях, так что в различных трассах отдельно взятое сейсмическое событие будет происходить в разные моменты времени. Если этап корректировки необработанных сейсмических данных к нулевому удалению выполнить правильно, то в каждой скорректированной трассе отдельно взятое отражение будет наблюдаться в один и тот же момент времени. Трассы из выборок 11 общей средней точки скорректированы по удалению с использованием коррекции нормального (или гиперболического) приращения времени (НПВ), при осуществлении которой предполагается, что скорость сейсмической энергии во внутренней области грунта является постоянной и изотропной, но изобретение не ограничено этой конкретной коррекцией по удалению.
Скорость, используемая при коррекции нормального приращения времени, может быть получена любым известным способом. Например, ее можно определить путем анализа скорости на предшествующем этапе обработки сейсмических данных, или ее можно получить в результате анализа скорости при предшествующей сейсмической разведке, выполненной в месте разведки. В качестве альтернативы скорость, используемую при коррекции нормального приращения времени, можно наилучшим образом оценить на основе любой другой доступной геофизической и/или географической информации о месте разведки.
В соответствии с изобретением суммирование по индексу к (то есть суммирование выборок трасс с общей средней точкой) не осуществляют при постоянном значении времени. Вместо этого суммирование по к осуществляют вдоль наклона отражающей поверхности, определяемого в результате анализа трасс или оцениваемого другими способами.
В предпочтительном варианте осуществления панель подобия формируют, используя соотношение
N / η+Ρ Μ V δ(.Σ 5*·] .
(2Р + 1)(М-1) Ν η+Ρ М * Σ Σ Σ4
(2)
Суммирование по _), к и ί в уравнении (2) в основном соответствует суммированию по ф к и ί в уравнении (1). Однако суммирование по к в уравнении (2) осуществляют вдоль наклона, тогда как суммирование по к в уравнении (1) примера из уровня техники осуществляют при постоянном значении времени. Поэтому в настоящем изобретении при формировании сейсморазведочной панели учитывается наклон.
Поскольку в настоящем изобретении учитывают наклон при формировании панели подобия, то большее число выборок трасс общей средней точки можно включить в суммирование до начала размывания максимумов подобия. Суммирование по к в уравнении (2) осуществляют от к=(п-р) до к=(п+р), что соответствует суммированию р соседних выборок трасс общей средней точки по каждую сторону от п-ой выборки трасс общей средней точки. В итоге это означает суммирование (2р+1) выборок трасс. В противоположность этому суммирование по к в уравнении (1) осуществляют от к=(п-1) до к=(п+1). Поскольку в уравнении (2) для формирования панели подобия используют (2р+1) выборок трасс общей средней точки, нормирующий множитель за пределами скобки в уравнении (2) равен 1/(2р+1). В противоположность этому нормирующий множитель в уравнении (1) равен 1/3, поскольку используют только 3 выборки трасс общей средней точки.
После того как панель подобия вычислена для п-ой выборки трасс, она может использоваться под управлением любого известного алгоритма автоматического выделения скорости, основанного на подобии. В качестве альтернативы она может использоваться оператором, занимающимся выделением скорости вручную.
Как отмечено выше, в предпочтительном варианте осуществления изобретения угол наклона определяют на основании сейсмических данных. Пример выполнения такой операции описан ниже со ссылкой на фиг. 4. В этом способе на основании сейсмических данных формируют набор панелей результатов суммирования с последовательностью многих скоростных функций, а угол наклона отражающей поверхности определяют из панелей результатов суммирования с последовательностью многих скоростных функций.
В способе использован алгоритм поиска наклона. Цель использования алгоритма поиска наклона заключается в определении поля наклонов в каждой точке анализа скорости. Поле наклонов представлено таблицей значений наклонов по одному значению наклона для каждого момента отсчета из панели результата суммирования с последовательностью многих скоростных функций. Чтобы получить поле одного наклона, сначала образуют отдельное поле наклонов для каждой панели результата суммирования с последовательностью многих скоростных функций. Затем эти поля наклонов объединяют с образованием поля одного наклона путем выбора со всех панелей в каждый момент времени вступлений с максимальной когерентностью.
На фиг. 4 показан способ определения наклона на примере отдельной панели в конкретный момент Т отсчета, то есть определения угла наклона отражающей поверхности, которая создает в сейсмических данных отраженную волну в момент Т времени. Показана одна панель 12 результата суммирования с последовательностью скоростных функций, по которой в момент Т времени выполняют алгоритм поиска. Процесс повторяют для всех моментов выборок на всех панелях результатов суммирования с последовательностью многих скоростных функций.
Поиск наклона осуществляют между границами «максимальный угол наклона» и «минимальный угол наклона». Эти границы характеризуют верхний и нижний пределы вероятного значения угла наклона и могут быть получены, например, из предсуществующей геологической информации о месте разведки. Сначала вычисляют подобие выборок трасс вдоль наклона при нижнем пределе угла наклона, при минимальном угле наклона. Затем увеличивают угол наклона и повторно вычисляют подобие, процесс повторяют для последующих увеличенных значений угла наклона до тех пор, пока не будет достигнут верхний предел угла наклона, максимальный угол наклона. Угол наклона, при котором получается наибольшее значение подобия, выбирают в качестве наклона в момент Т времени. Чтобы выбрать наклон для поля общего наклона, значение подобия также используют для сравнения с когерентностью вступлений, идентифицированных в то же самое время на других панелях результатов суммирования с последовательностью многих скоростных функций.
Для поиска наклона требуется значительный объем вычислительных работ. Поэтому при практической реализации способа описанный выше процесс можно не проводить по всем панелям результатов суммирования с последовательностью многих скоростных функций или при всех моментах взятия отсчетов из панелей результатов суммирования с последовательностью многих скоростных функций, по которым осуществляют определение. Можно уменьшить объем требуемых вычислений путем выбора панелей результатов суммирования с последовательностью многих скоростных функций и определения поля наклонов на некотором временном интервале взятия отсчетов из первых отсчетов. После этого промежуточные значения поля наклонов могут быть вычислены путем интерполяции.
Для удобства в описанных выше вариантах осуществления предполагалось, что имеется единственная отражающая поверхность, хотя на практике сейсмические данные содержат вступления волн от более чем одной отражающей поверхности. Изобретение не ограничено сейсмическими данными, содержащими вступления волн только от одной отражающей поверхности, а может быть применено к сейсмическим данным, содержащим отраженные волны от более чем одной отражающей поверхности.
В способе, описанном со ссылкой на фиг. 4, наклон отражающей поверхности определяют в момент взятия каждого отсчета. Поскольку данные от различных отражающих поверхностей будут приходить в разные моменты взятия отсчетов, то, если этот способ применить к сейсмическим данным, содержащим данные от более чем одной отражающей поверхности, будет определяться наклон каждой отражающей поверхности, от которой в сейсмических данных имеется отраженная волна (если в конкретный момент времени взятия отсчета волны не обнаружено или если в конкретный момент времени взятия отсчета наклон не обнаружен, то угол наклона в этот момент взятия отсчета принимается равным нулю). После этого значение угла наклона, определенное для каждой отражающей поверхности, можно использовать при формировании панели подобия.
В варианте осуществления изобретения, в котором наклон отражающей поверхности определяют на основании предшествующей информации о месте разведки или на основании предшествующих данных, полученных в месте разведки, желательно получать отдельные значения угла наклона для каждой отражающей поверхности, что даст возможность усилить вступления волны в сейсмических данных.
На фиг. 5 приведена блок-схема устройства 13 для обработки данных согласно настоящему изобретению. В соответствии со способом в устройстве могут обрабатываться выборки трасс из принятой сейсмической энергии для получения панели подобия.
Устройство 13 содержит программируемый процессор 14 данных с программным запоминающим устройством 15, например в виде постоянного запоминающего устройства, сохраняющего программу для управления процессором 14 данных для обработки трасс сейсмических данных заявленным способом. Устройство также содержит долговременную память 16 с чтением/записью для хранения, например, данных, которые должны сохраняться в отсутствие питания. «Рабочая» или «сверхоперативная память» для процессора данных образована запоминающим устройством 17 с произвольным доступом. Устройство 18 ввода данных предусмотрено для приема данных и команд пользователя. Устройство вывода 19 предусмотрено для отображения информации, относящейся к ходу и результату процесса обработки. Устройством вывода может быть, например, принтер, визуальный индикатор или выходное запоминающее устройство.
Сейсмические данные, предназначенные для обработки, могут поступать через устройство 18 ввода данных или по выбору могут выда13 ваться машинно-считываемым запоминающим устройством 20.
Программа, предназначенная для работы системы и осуществления способа, описанного выше, хранится в программном запоминающем устройстве 15, которое может быть реализовано в виде полупроводникового запоминающего устройства, например в виде хорошо известного постоянного запоминающего устройства. Однако программа может успешно храниться в любом другом подходящем запоминающем средстве, например, на магнитном носителе 15а данных («дискете») или на компакт-диске 15Ь постоянной памяти.

Claims (15)

1. Способ обработки сейсмических данных, сформированных при отражении сейсмической энергии от отражающей поверхности, заключающийся в том, что осуществляют выбор η-ой выборки трасс и по меньшей одной другой, ближайшей выборки трасс, формируют панель подобия, связанную с местом, соответствующим η-ой выборке из выбранных выборок, с учетом наклона отражающей поверхности.
2. Способ по п.1, отличающийся тем, что при формировании панели подобия осуществляют суммирование выбранных выборок вдоль наклона отражающей поверхности.
3. Способ по п.1 или 2, отличающийся тем, что на этапе выбора осуществляют выбор η-ой выборки и р (р>1) соседних выборок по каждую сторону от η-ой выборки.
4. Способ по любому из пп.1-3, отличающийся тем, что при формировании панели подобия осуществляют суммирование
N ( η+Ρ Μ V
Ν η+Ρ М *
Σ Σ ΣΑ» 1=1 к=п-р 1=1 ч 7 причем суммирование по к как числителя, так и знаменателя дроби осуществляют вдоль наклона отражающий поверхности, где Аут представляет собой отсчет трассы, М - число трасс в каждой выборке, N - число временных отсчетов.
5. Способ по п.3 или 4, отличающийся тем, что используют р=10.
6. Способ по любому из пп.1-5, отличающийся тем, что дополнительно определяют на
Подобие клон отражающей поверхности на основании сейсмических данных.
7. Способ по п.6, отличающийся тем, что при определении наклона определяют набор панелей результатов суммирования с последовательностью многих скоростных функций на основании выбранных выборок и определяют наклон на основании панелей результатов суммирования с последовательностью многих скоростных функций.
8. Способ по любому из пп.1-5, отличающийся тем, что дополнительно определяют наклон отражающей поверхности на основании предшествующих данных.
9. Способ по любому из пп.1-8, отличающийся тем, что выборки сейсмических трасс представляют собой выборки общей средней точки.
10. Устройство для обработки сейсмических данных, сформированных при отражении сейсмической энергии от отражающей поверхности, содержащее средство для выбора η-ой выборки трасс из принятой сейсмической энергии и по меньшей мере одной другой ближайшей выборки трасс, средство для формирования панели подобия, связанной с местом, соответствующим η-ой выборке из выбранных выборок, с учетом наклона отражающей поверхности.
11. Устройство по п.10, отличающееся тем, что дополнительно содержит средство для суммирования выбранных выборок трасс вдоль наклона отражающей поверхности.
12. Устройство по любому из пп.10, 11, отличающееся тем, что дополнительно содержит средство для определения наклона отражающей поверхности на основании сейсмических данных.
13. Устройство по п.12, отличающееся тем, что содержит средство для определения набора панелей результатов суммирования с последовательностью многих скоростных функций на основании выбранных выборок трасс и для определения наклона на основании панелей результатов суммирования с последовательностью многих скоростных функций.
14. Устройство по любому из пп.10-13, отличающееся тем, что содержит программируемый процессор данных.
15. Запоминающее средство, содержащее программу для процессора данных для устройства по п.14.
EA200300769A 2001-01-05 2002-01-07 Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности EA004857B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GBGB0100207.0A GB0100207D0 (en) 2001-01-05 2001-01-05 A method of processing seismic data
PCT/GB2002/000044 WO2002059648A1 (en) 2001-01-05 2002-01-07 A seismic method and apparatus for generating a semblance panel and computing the reflector dip

Publications (2)

Publication Number Publication Date
EA200300769A1 EA200300769A1 (ru) 2003-10-30
EA004857B1 true EA004857B1 (ru) 2004-08-26

Family

ID=9906261

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200300769A EA004857B1 (ru) 2001-01-05 2002-01-07 Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности

Country Status (10)

Country Link
US (1) US6952379B2 (ru)
EP (1) EP1348138B1 (ru)
CN (1) CN1308703C (ru)
AU (1) AU2002216298B2 (ru)
CA (1) CA2433544A1 (ru)
DE (1) DE60235606D1 (ru)
EA (1) EA004857B1 (ru)
GB (1) GB0100207D0 (ru)
NO (1) NO20033028L (ru)
WO (1) WO2002059648A1 (ru)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102006003535A1 (de) * 2006-01-24 2007-08-02 Schott Ag Verfahren zur Temperaturbeeinflussung einer Schmelze
CN101419293B (zh) * 2007-10-25 2011-05-25 中国石油天然气集团公司 一种提高地震数据信噪比的方法
FR2946153B1 (fr) * 2009-05-27 2011-06-10 Cggveritas Services Sa Procede de surveillance d'une zone du sous-sol, notamment lors d'operations de fracturation provoquee
CN102096103B (zh) * 2010-12-03 2012-08-15 中国石油天然气集团公司 用于低信噪比地震资料的速度分析方法
MX2014004128A (es) * 2011-11-07 2015-04-14 Landmark Graphics Corp Sistemas y metodos de generacion de imagenes sismicas que emplean apilamiento basado en correlacion.
GB2505042B (en) * 2012-06-13 2014-10-08 Schlumberger Holdings Seismic trace attribute
US11531129B2 (en) 2019-05-30 2022-12-20 Saudi Arabian Oil Company Picking seismic stacking velocity based on structures in a subterranean formation
US11899147B2 (en) * 2021-08-12 2024-02-13 Saudi Arabian Oil Company Method and system for seismic denoising using omnifocal reformation

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1011351B (zh) * 1985-06-19 1991-01-23 斯托拉尔公司 煤层结构成象的勘测方法
US4994747A (en) * 1988-01-14 1991-02-19 Stolar, Inc. Method and apparatus for detecting underground electrically conductive objects
US5153858A (en) * 1991-07-09 1992-10-06 Landmark Graphics Corporation Method for finding horizons in 3D seismic data
US5930730A (en) * 1994-12-12 1999-07-27 Amoco Corporation Method and apparatus for seismic signal processing and exploration
GB9508525D0 (en) * 1995-04-27 1995-06-14 Geco As Method of processing seismic data
FR2751757B1 (fr) * 1996-07-25 1998-08-28 Elf Aquitaine Methode d'acquisition et de traitement de donnees sismiques reflexion pour l'exploration d'un milieu a tectonique complexe
US5933789A (en) * 1996-11-12 1999-08-03 Atlantic Richfield Company Method and system for applying dispersive normal moveout corrections to seismic survey signals
US5978314A (en) 1997-03-21 1999-11-02 Exxon Production Research Company Method for determining seismic velocities
GB9904101D0 (en) 1998-06-09 1999-04-14 Geco As Subsurface structure identification method

Also Published As

Publication number Publication date
NO20033028D0 (no) 2003-07-01
NO20033028L (no) 2003-09-05
CN1308703C (zh) 2007-04-04
DE60235606D1 (de) 2010-04-22
CA2433544A1 (en) 2002-08-01
EA200300769A1 (ru) 2003-10-30
EP1348138A1 (en) 2003-10-01
US20040196739A1 (en) 2004-10-07
EP1348138B1 (en) 2010-03-10
AU2002216298B2 (en) 2005-09-22
WO2002059648A1 (en) 2002-08-01
CN1484770A (zh) 2004-03-24
US6952379B2 (en) 2005-10-04
GB0100207D0 (en) 2001-02-14

Similar Documents

Publication Publication Date Title
US6058074A (en) Method and system for detecting hydrocarbon reservoirs using amplitude-versus-offset analysis with improved measurement of background statistics
CA2558000C (en) Method, media, and signals for processing seismic data to obtain a velocity field
Schneider Developments in seismic data processing and analysis (1968–1970)
RU2694621C1 (ru) Способ и устройство для обработки сейсмических данных
US4964103A (en) Three dimensional before stack depth migration of two dimensional or three dimensional seismic data
US8687463B2 (en) Method for processing at least two sets of seismic data
WO2000031567A1 (en) Hydrocarbon edge detection using seismic amplitude
US5978314A (en) Method for determining seismic velocities
US5073876A (en) Geophysical exploration using near surface structure corrections developed from common endpoint gather stacked traces
AU2004318849B2 (en) Generalized 3D surface multiple prediction
EA004857B1 (ru) Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности
US4403313A (en) Process and an apparatus for seismic geophysics with processing by focuses
US6249746B1 (en) Automated seismic isochron analysis
US6980482B2 (en) Seismic survey system
US20060253257A1 (en) Method of processing seismic data for AVO or AVOA characterization
AU2002216298A1 (en) A seismic method and apparatus for generating a semblance panel and computing the reflector dip
CN104570115A (zh) 一种面波衰减方法及装置
CA2533011C (en) Method for bispectral picking of anelliptical nmo correction parameters
US5629905A (en) Method of datuming seismic data and method of processing seismic data
RU2779518C1 (ru) Способ определения минимально достаточного размера бина для размещения источников и приемников при проведении сейсмической съемки
CN116256801B (zh) 基于图像融合的深地油气精准导航断层表征方法与系统
Lee et al. Pre-stack sea bottom detection and swell correction within an expected hyperbolic range for noisy high-resolution 8-channel airgun seismic data
CN114460636A (zh) 基于炮检对信噪比地震照明方法、装置、电子设备及介质
CN115774288A (zh) 横波矢量反演建模方法及装置
WO1998020368A1 (en) Aligning seismic traces

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 AZ BY KZ KG MD 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): RU