EA004857B1 - Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности - Google Patents
Сейсмические способ и устройство для формирования панели подобия и вычисления наклона отражающей поверхности Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000012545 processing Methods 0.000 claims abstract description 18
- 230000006870 function Effects 0.000 claims description 34
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000007423 decrease Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000010845 search algorithm Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 241000234282 Allium Species 0.000 description 1
- 235000002732 Allium cepa var. cepa Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis 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.
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)
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 | 中国石油天然气集团公司 | 用于低信噪比地震资料的速度分析方法 |
EP2748740A4 (en) * | 2011-11-07 | 2015-12-16 | Landmark Graphics Corp | SEISMIC IMAGING SYSTEMS AND METHODS EMPLOYING A STACK BASED ON CORRELATION |
NO345771B1 (no) * | 2012-06-13 | 2021-07-26 | Schlumberger Technology Bv | Seismisk sporingsattributt |
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)
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 |
-
2001
- 2001-01-05 GB GBGB0100207.0A patent/GB0100207D0/en not_active Ceased
-
2002
- 2002-01-07 US US10/250,849 patent/US6952379B2/en not_active Expired - Fee Related
- 2002-01-07 AU AU2002216298A patent/AU2002216298B2/en not_active Ceased
- 2002-01-07 EP EP02734864A patent/EP1348138B1/en not_active Expired - Lifetime
- 2002-01-07 WO PCT/GB2002/000044 patent/WO2002059648A1/en not_active Application Discontinuation
- 2002-01-07 EA EA200300769A patent/EA004857B1/ru not_active IP Right Cessation
- 2002-01-07 CA CA002433544A patent/CA2433544A1/en not_active Abandoned
- 2002-01-07 DE DE60235606T patent/DE60235606D1/de not_active Expired - Lifetime
- 2002-01-07 CN CNB028034732A patent/CN1308703C/zh not_active Expired - Fee Related
-
2003
- 2003-07-01 NO NO20033028A patent/NO20033028L/no not_active Application Discontinuation
Also Published As
Publication number | Publication date |
---|---|
EP1348138B1 (en) | 2010-03-10 |
US6952379B2 (en) | 2005-10-04 |
DE60235606D1 (de) | 2010-04-22 |
US20040196739A1 (en) | 2004-10-07 |
WO2002059648A1 (en) | 2002-08-01 |
CN1484770A (zh) | 2004-03-24 |
AU2002216298B2 (en) | 2005-09-22 |
CN1308703C (zh) | 2007-04-04 |
GB0100207D0 (en) | 2001-02-14 |
NO20033028D0 (no) | 2003-07-01 |
CA2433544A1 (en) | 2002-08-01 |
EA200300769A1 (ru) | 2003-10-30 |
EP1348138A1 (en) | 2003-10-01 |
NO20033028L (no) | 2003-09-05 |
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 | |
EP1141750A1 (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 | |
CN116256801B (zh) | 基于图像融合的深地油气精准导航断层表征方法与系统 | |
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 | |
US4835745A (en) | Method of radial seismic data collection for improved salt structure imaging | |
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) | Способ определения минимально достаточного размера бина для размещения источников и приемников при проведении сейсмической съемки | |
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 |