RU2694621C1 - Способ и устройство для обработки сейсмических данных - Google Patents
Способ и устройство для обработки сейсмических данных Download PDFInfo
- Publication number
- RU2694621C1 RU2694621C1 RU2018105481A RU2018105481A RU2694621C1 RU 2694621 C1 RU2694621 C1 RU 2694621C1 RU 2018105481 A RU2018105481 A RU 2018105481A RU 2018105481 A RU2018105481 A RU 2018105481A RU 2694621 C1 RU2694621 C1 RU 2694621C1
- Authority
- RU
- Russia
- Prior art keywords
- seismic
- selection
- correlation coefficient
- traces
- time
- Prior art date
Links
Images
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/366—Seismic filtering by correlation of seismic signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
- G01V2210/532—Dynamic changes in statics, e.g. sea waves or tidal influences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Abstract
Изобретение относится к области геофизики и может быть использовано для обработки сейсмических данных. В настоящем раскрытии представлены способ и устройство для обработки сейсмических данных, относящиеся к области геологических исследований. Способ содержит суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-й подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-й подборки сейсмических трасс в обе стороны по порядку; и выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума. Технический результат – повышение точности получаемых данных. 3 н. и 14 з.п. ф-лы, 4 ил.
Description
ОБЛАСТЬ ТЕХНИКИ
Для настоящего раскрытия испрашивается приоритет заявки CN201610204476.6 на патент Китая, поданной 1 апреля 2016 года под названием ʺMethod and apparatus for processing seismic dataʺ, все содержание которой включено в эту заявку путем ссылки.
УРОВЕНЬ ТЕХНИКИ
В способах сейсмической разведки исследуют подземные геологические условия при распространении в пластах искусственно возбуждаемой сейсмической волны. Когда сейсмическая волна, возбуждаемая в определенном месте на грунте, распространяется под землей, отраженные волны или преломленные волны образуются на границах пласта вследствие различной упругости и возвращаются на земную поверхность. Специальные приборы можно использовать для регистрации этих волн и анализа характеристик регистрируемых волн (например, времени распространения и формы пульсации волн), чтобы точно измерять глубины и формы этих границ и оценивать литологию, физические и жидкостные свойства пластов.
Поскольку сейсмическая запись является зашумленной, сейсмические волны обычно возбуждают и принимают на различных местах и получают сейсмические волны, отраженные от одной и той же подземной точки, и набор этих сейсмических волн, отраженных от одной и той же точки, называют подборкой сейсмических трасс. Путем суммирования сигналов шум может быть исключен и отношение сигнала к шуму может быть повышено.
В сейсмических данных, собираемых в полевых условиях, для одной и той же точки отражения будет иметься время отражения сейсмической волны, которое зависит от удаления и скорости распространения сейсмической волны, и поэтому для решения задачи суммирования сейсмических сигналов необходимо корректировать время регистрации сейсмических волн для различных удалений в соответствии с временной задержкой между автономным возбуждением и автономным приемом отраженных сейсмических волн, и этот процесс называют коррекцией за нормальное приращение времени (НПВ).
Однако скорость распространения сейсмической волны является неизвестной, поэтому коррекцию за нормальное приращение времени обычно реализуют, изначально предусматривая ряд скоростей распространения сейсмических волн для выполнения процесса коррекции за нормальное приращение времени подборки сейсмических трасс.
Поскольку объем сейсмических данных, собираемых в полевых условиях, очень большой, человеку невозможно выполнить скоростной анализ для каждой точки отражения сейсмической волны. Обычно скоростной анализ выполняют для одной точки отражения сейсмической волны на интервале нескольких десятков трасс, при этом скорость из подборки трасс получают без скоростного анализа путем интерполяции скорости из подборки трасс для, которой скоростной анализ был выполнен. Поскольку скорости сейсмических волн не являются одинаковыми, интерполированные скорости в большей или меньшей степени имеют погрешности, и результат коррекции за нормальное приращение времени не может приводить к достижению идеального эффекта, поэтому трудно выправлять события во всех подборках сейсмических трасс в горизонтальное состояние.
СУЩНОСТЬ РАСКРЫТИЯ
Согласно вариантам осуществления настоящего раскрытия для разрешения проблем из предшествующего уровня техники предложены способ и устройство для обработки сейсмических данных.
Согласно вариантам осуществления настоящего раскрытия предложен способ обработки сейсмических данных, содержащий суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; и выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Кроме того, согласно вариантам осуществления настоящего раскрытия предложено устройство для обработки сейсмических данных, содержащее запоминающее устройство, процессор и компьютерную программу, сохраняемую в запоминающем устройстве и выполняемую процессором, при этом процессор выполняет следующие операции при выполнении компьютерной программы: суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Кроме того, согласно вариантам осуществления настоящего раскрытия предложен машиночитаемый носитель данных, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять следующие операции: суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Применением технических решений согласно вариантам осуществления достигаются следующие положительные эффекты: при суммировании множества подборок сейсмических трасс; вычислении коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выборе K-той подборки сейсмических трасс; вычислении точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; и выполнении остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума это решение объединяют с решением из предшествующего уровня техники, в котором реализуется процесс коррекции за нормальное приращение времени при анализе скоростей распространения сейсмических волн, так что события всех подборок модельных трасс выправляются в горизонтальное состояние, в результате чего исключаются погрешности, обусловленные интерполяцией, в решении из предшествующего уровня техники, и повышается точность коррекции за нормальное приращение времени до суммирования.
Следует понимать, что приведенное выше общее описание и последующие конкретные варианты осуществления являются только примерными и иллюстративными, а не ограничивающими объем настоящего раскрытия.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Чтобы более ясно описать технические решения согласно вариантам осуществления настоящего раскрытия, при описании вариантов осуществления будут использоваться чертежи, которые кратко представлены ниже. Очевидно, что в нижеследующем описании чертежами только иллюстрируются некоторые варианты осуществления настоящего раскрытия, и специалист в данной области техники может получить другие чертежи без затрат творческих усилий.
На чертежах:
фиг. 1 - блок-схема последовательности действий способа обработки сейсмических данных согласно варианту осуществления настоящего раскрытия;
фиг. 2 - блок-схема последовательности действий еще одного способа обработки сейсмических данных согласно варианту осуществления настоящего раскрытия;
фиг. 3а - схематичное представление подборок сейсмических трасс при обычной обработке сейсмических данных согласно варианту осуществления настоящего раскрытия;
фиг. 3b - схематичное представление подборок сейсмических трасс после процесса выравнивания энергий амплитуд согласно варианту осуществления настоящего раскрытия;
фиг. 3с - схематичное представление подборок сейсмических трасс после процесса остаточной коррекции за нормальное приращение времени подборок сейсмических трасс из фиг. 3b согласно варианту осуществления настоящего раскрытия; и
фиг. 4 - структурная схема устройства для обработки сейсмических данных согласно варианту осуществления настоящего раскрытия.
ПОДРОБНОЕ ОПИСАНИЕ ВАРИАНТОВ ОСУЩЕСТВЛЕНИЯ
Чтобы задача, технические решения и преимущества вариантов осуществления настоящего раскрытия стали понятней, ниже с обращением к чертежам и подробному описанию будет доходчиво пояснена сущность содержания настоящего раскрытия. После усвоения вариантов осуществления настоящего раскрытия любой специалист в данной области техники сможет изменять и модифицировать технические средства и способы, усвоенные из настоящего раскрытия, без отступления от сущности и объема настоящего раскрытия.
На фиг. 1 показана блок-схема последовательности действий способа обработки сейсмических данных согласно варианту осуществления настоящего раскрытия. Как показано на фиг. 1, способ содержит этап 101 суммирования подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом.
При этом все S подборок сейсмических трасс представляют собой подборки трасс общей точки отражения (ОТО) и подборки трасс общей средней точки (ОСТ) и регистрируются при различных удалениях для получения подборок трасс с различными порядковыми номерами.
На фиг. 3а показано схематичное представление подборок сейсмических трасс, где по горизонтальной координате показано удаление (единица измерения: метр) и каждый сегмент обозначает идентичное удаление, тогда как по вертикальной координате показано время (единица измерения; миллисекунда) и каждый сегмент обозначает подборку сейсмической трассы в зависимости от удаления, и эти подборки являются набором сейсмических трасс в зависимости от удаления.
При этом перед этапом 101 подборки сейсмических трасс обрабатывают способом коррекции за нормальное приращение времени, при выполнении которого анализируют скорости распространения сейсмических волн. В этом варианте осуществления оценивают набор скоростей суммирования сейсмических волн, чтобы вычислить время отражения сейсмической волны; значение разности между временем отражения сейсмической волны и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении выбирают в качестве величины поправки за нормальное приращение времени; выполняют коррекцию за нормальное приращение времени исходных подборок сейсмических трасс; и ведут наблюдение, чтобы определить, что текущая скорость суммирования сейсмических волн является надлежащей, когда достигается горизонтальное состояние подборок трасс при этой скорости суммирования сейсмических волн. Процесс наблюдения оператор обычно завершает на глазок.
Время отражения сейсмической волны может быть вычислено по следующей формуле:
где tx обозначает время отражения сейсмической волны, x обозначает удаление, t0 обозначает время отражения сейсмической волны, при котором удаление равно 0, и va обозначает скорость суммирования сейсмической волны.
При этом заданный интервал может означать все S подборок сейсмических трасс или часть их.
Суммирование подборок сейсмических трасс означает добавление амплитуд подборок сейсмических трасс различных порядковых номеров на основании выборочных точек.
Этап 102: вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции.
Этап 103: вычисление по порядку величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
При этом величина остаточной поправки за нормальное приращение времени означает величину поправки, используемую в процессе дальнейшей коррекции сейсмических данных после выполнения процесса коррекции за нормальное приращение времени, при котором анализируют скорости распространения сейсмических волн.
В варианте осуществления этап 103 реализуют на следующих этапах, на которых выполняют:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции M-той подборки сейсмических трасс с (M+1)-ой подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для M-той подборки сейсмических трасс, где K-2≥M≥1 и M является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс (например, в данном случае (N+1)-ой подборки сейсмических трасс) во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
При этом время отражения сейсмической волны из подборки сейсмических трасс во временном окне может быть средним значением времен отражения сейсмических волн в соответствующих сейсмических трассах из подборки сейсмических трасс в точке оптимума или временем отражения сейсмической волны в любой одной из сейсмических трасс.
Этап 104: выполнение остаточной коррекции за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
При этом остаточная коррекция за нормальное приращение времени означает процесс дальнейшей коррекции сейсмических данных после выполнения способа коррекции за нормальное приращение времени, при выполнении которого анализируют скорости распространения сейсмических волн.
В варианте осуществления настоящего раскрытия при суммировании множества подборок сейсмических трасс; вычислении коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выборе K-той подборки сейсмических трасс; вычислении точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; и выполнении остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума это решение объединяют с решением из предшествующего уровня техники, в котором реализуется процесс коррекции за нормальное приращение времени при анализе скоростей распространения сейсмических волн, так что события во всех подборках модельных трасс выправляются в горизонтальное состояние, в результате чего исключаются погрешности, обусловленные интерполяцией, в решении из предшествующего уровня техники, и повышается точность коррекции за нормальное приращение времени до суммирования.
На фиг. 2 показана блок-схема последовательности действий другого способа обработки сейсмических данных согласно варианту осуществления настоящего раскрытия. Как показано на фиг. 2, способ содержит этап 201 выполнения процесса выравнивания энергий амплитуд каждой подборки сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ).
При этом процесс выравнивания энергий амплитуд каждой подборки сейсмических трасс выполняют на следующих этапах.
На первом этапе подсчитывают энергию Е0 амплитуд в заданном интервале временного окна (например, от 1000 до 3000 мс):
где Q обозначает количество подсчитанных трасс (суммарное количество подборок трасс), Р обозначает количество выборочных точек и aj обозначает значение амплитуды.
На втором этапе вычисляют энергию Е1 амплитуд для каждой подборки сейсмических трасс:
На третьем этапе вычисляют коэффициент k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1.
На четвертом этапе выполняют процесс выравнивания амплитуд для каждой подборки сейсмических трасс, используя коэффициент выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс (значений амплитуд множества выборочных точек), и обозначает значение амплитуды после процесса выравнивания амплитуд.
В варианте осуществления до этапа 201 способ также содержит 1) возбуждение сейсмических волн управляемым сейсмическим источником и регистрацию данных о сейсмических волнах; 2) регистрацию при использовании управляемого сейсмического источника и выполнение анализа поглощения в выбранном временном окне для получения амплитудного спектра в пунктах взрыва или пунктах приема в выбранном временном окне: 3) реализацию процесса коррекции за нормальное приращение времени при анализе скоростей распространения сейсмических волн.
Этап 202: суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс для получения модельной трассы, при этом S является целым числом.
В варианте осуществления нижеследующую формулу используют для суммирования подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ), чтобы получить модельную трассу:
где m1 и m2 обозначают порядковые номера подборок сейсмических трасс и aij обозначает значение амплитуды j-той выборочной точки из i-той трассы; m1, m2, i и j являются положительными целыми числами, 1≤m1<m2≤S, i∈[m1, m2].
При этом различные порядковые номера подборок сейсмических трасс показывают, что эти подборки сейсмических трасс соответствуют различным удалениям; а выборочные точки выбирают из амплитудного спектра в заданном временном интервале.
При этом заданный интервал [m1, m2] может быть интервалом всех или части S подборок сейсмических трасс.
При этом значение амплитуды выборочной точки может быть средним значением положений во времени подборки сейсмических трасс в выборочной точке или значением любой сейсмической трассы из подборки сейсмических трасс.
Этап 203: вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции.
Нижеследующую формулу используют для вычисления коэффициента корреляции модельной трассы с подборкой сейсмических трасс:
где обозначает коэффициент корреляции x-вой трассы с модельной трассой, aj обозначает значение амплитуды j-той выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-той выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
Этап 204: выполняемое по порядку вычисление величин остаточных поправок за нормальное приращение времени (НПВ) для S подборок сейсмических трасс.
В варианте осуществления этап 204 реализуют на следующих этапах, на которых выполняют:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции M-той подборки сейсмических трасс с (M+1)-ой подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для M-той подборки сейсмических трасс, где K-2≥M≥1 и M является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс (например, в данном случае (N+1)-ой подборки сейсмических трасс) во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
Оптимальное положение каждой подборки сейсмических трасс определяют по порядку от K-той трассы, коэффициент корреляции каждой подборки сейсмических трасс с модельной трассой вычисляют при последовательном скольжении по подборке сейсмических трасс, положение выбирают в качестве точки оптимума в подборке сейсмических трасс, когда коэффициент корреляции является максимальным, и величину временного сдвига, соответствующую точке оптимума обозначают Δt. При упомянутом выше вычислении первоначально выбранную K-тую трассу используют в качестве модельной трассы, в то время как при последующих вычислениях предшествующую трассу используют, в свою очередь, в качестве модельной трассы для текущей трассы, например, (М+1)-ую трассу используют в качестве модельной трассы для М-той трассы и N-ную трассу используют в качестве модельной трассы для (N+1)-ой трассы.
При этом коэффициент корреляции подборки сейсмических трасс с модельной трассой во временном окне вычисляют по формуле:
где обозначает коэффициент корреляции подборки сейсмических трасс с модельной трассой во временном окне, в котором расположена j-тая выборочная точка, aj обозначает значение амплитуды j-той выборочной точки в модельной трассе , Aa обозначает среднее значение амплитуд модельной трассы во временном окне, bj обозначает значение амплитуды j-той выборочной точки в подборке сейсмических трасс и Ba обозначает среднее значение амплитуд подборки сейсмических трасс во временном окне, где j является положительным целым числом.
При этом заданный временной интервал может быть от T0-dT до T1-dT, протяженность временного окна может быть Т (от Т0 до Т1), а конкретные значения могут быть выбраны на основании текущих требований. При упомянутом выше вычислении два скольжения временного окна перекрываются друг с другом, так что выбранная точка оптимума имеет достаточную точность.
Этап 205: выполнение остаточной коррекции за нормальное приращение времени (НПВ) S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок для S подборок сейсмических трасс.
В варианте осуществления прежде всего получают величину временного сдвига, соответствующую точке оптимума в каждой подборке сейсмических трасс, в качестве величины Δt поправки; затем Δt используют в качестве величины поправки за нормальное приращение времени, чтобы скорректировать каждую подборку сейсмических трасс.
В еще одном варианте осуществления настоящего раскрытия способ обработки сейсмических данных содержит суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс по порядку от K-той подборки сейсмических трасс в обе стороны; и выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
В варианте осуществления вычисление точки оптимума каждой подборки сейсмических трасс по порядку от K-той подборки сейсмических трасс в обе стороны содержит скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-той подборки сейсмических трасс и (М+1)-ой подборки сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-той подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
В варианте осуществления выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума содержит выполнение остаточной коррекции за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
Далее с обращением к фиг. 3а-3с описываются положительные эффекты способа согласно варианту осуществления настоящего раскрытия.
На фиг. 3а показаны подборки сейсмических трасс после обычной сейсмической обработки (коррекции за нормальное приращение времени), на которых видно, что энергия амплитуд (черные участки) из подборки сейсмических трасс, имеющей небольшое расстояние между пунктом взрыва и геофоном (на левой стороне), значительно слабее (черные участки имеют меньшую ширину), чем энергия амплитуд из подборки сейсмических трасс, имеющей большее расстояние между пунктом взрыва и геофоном (на правой стороне), и в таком случае очевидно, что распределение энергий амплитуд не совпадает с характеристиками распределения энергии амплитуд при реальном отражении сейсмической волны. Таким образом, необходимо выполнить процесс выравнивания энергий амплитуд. На фиг. 3b показаны подборки сейсмических трасс после процесса выравнивания энергий амплитуд. Поскольку события в подборках трасс являются нерегулярными, характеристики изменения амплитуды, равно как и удаления, не могут быть правильно представлены только выполнением процесса выравнивания энергий амплитуд, и процесс остаточной коррекции за нормальное приращение времени модельной трасы также должен быть выполнен для выправления событий в подборках трасс до суммирования в горизонтальное состояние (фиг. 3с).
В вариантах осуществления настоящего раскрытия при суммировании множества подборок сейсмических трасс; вычислении коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выборе K-той подборки сейсмических трасс; вычислении точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку и выполнении коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума это решение объединяют с решением из предшествующего уровня техники, в котором реализуется процесс коррекции за нормальное приращение времени при анализе скоростей распространения сейсмических волн, так что события во всех подборках модельных трасс выправляются в горизонтальное состояние, в результате чего исключаются погрешности, обусловленные интерполяцией, в решении из предшествующего уровня техники и повышается точность коррекции до суммирования за нормальное приращение времени; поэтому качество результата суммирования сейсмических данных после оптимизации подборки трасс безусловно улучшается или повышается и вследствие этого можно более точно описывать структуру пласта, литологию и свойства флюида.
На фиг. 4 показана структурная схема устройства для обработки сейсмических данных согласно варианту осуществления настоящего раскрытия. Как показано на фиг. 4, устройство содержит запоминающее устройство 401, процессор 402, а компьютерная программа сохраняется в запоминающем устройстве 401 и выполняется в процессоре 402, при этом процессор 402 при выполнении компьютерной программы выполняет следующие операции:
суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, при этом S является целым числом;
вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-той подборки сейсмических трасс с (М+1) подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-той подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K=1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс (например, в данном случае (N+1)-ой подборки сейсмических трасс) во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс с нулевым удалением;
выполнение остаточной коррекции за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок для S подборок сейсмических трасс.
В варианте осуществления настоящего раскрытия процессор суммирует множество подборок сейсмических трасс; вычисляет коэффициент корреляции каждой подборки сейсмических трасс с модельной трассой и выбирает K-тую подборку сейсмических трасс; вычисляет точку оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; и выполняет остаточную коррекцию за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума. Это решение объединяется с решением из предшествующего уровня техники, в котором реализуется процесс коррекции за нормальное приращение времени при анализе скоростей распространения сейсмических волн, так что события во всех подборках модельных трасс выправляются в горизонтальное состояние, в результате чего исключаются погрешности, обусловленные интерполяцией, в решении из предшествующего уровня техники и повышается точность коррекции за нормальное приращение времени до суммирования.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет процесс выравнивания энергий амплитуд каждой подборки сейсмических трасс перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс, после процесса коррекции за нормальное приращение времени.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет следующие операции:
подсчет энергии Е0 амплитуд в заданном интервале временного окна:
где Q обозначает количество подсчитанных трасс, Р обозначает количество выборочных точек и aj обозначает значение амплитуды;
вычисление энергии Е1 амплитуд для каждой подборки сейсмических трасс:
вычисление коэффициента k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1;
выполнение процесса выравнивания амплитуд для каждой подборки сейсмических трасс при использовании коэффициента выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс (значений амплитуд множества выборочных точек), и обозначает значение амплитуды после процесса выравнивания амплитуд;
суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, при этом S является целым числом;
вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-той подборки сейсмических трасс с (М+1)-ой подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-той подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени (K+1)-ой подборки сейсмических трасс;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс (например, в данном случае (N+1)-ой подборки сейсмических трасс) во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс с нулевым удалением;
выполнение остаточной коррекции за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок для S подборок сейсмических трасс.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 использует нижеследующую формулу для суммирования подборок сейсмических трасс, чтобы получить модельную трассу:
где m1 и m2 обозначают порядковые номера подборок сейсмических трасс и aij обозначает значение амплитуды j-той выборочной точки из i-той трассы; m1, m2, i и j являются положительными целыми числами, 1≤m1<m2≤S, i∈[m1, m2].
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор использует нижеследующую формулу для вычисления коэффициента корреляции модельной трассы с подборкой сейсмических трасс:
где обозначает коэффициент корреляции x-вой трассы с модельной трассой, aj обозначает значение амплитуды j-той выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-той выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет следующие операции: суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, при этом S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет следующие операции: скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-той подборки сейсмических трасс с (М+1)-ой подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-той подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс; скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет остаточную коррекцию за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
Кроме того, в варианте осуществления при выполнении компьютерной программы процессор 402 выполняет следующие операции: выполнение процесса выравнивания энергий амплитуд каждой подборки сейсмических трасс перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс, после процесса коррекции за нормальное приращение времени; суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы, где S является целым числом; вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции; вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Кроме того, согласно вариантам осуществления настоящего раскрытия предусмотрен машиночитаемый носитель данных, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять следующие операции:
суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом;
вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-той подборки сейсмических трасс с максимальным коэффициентом корреляции;
вычисление точки оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку;
выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
Кроме того, в варианте осуществления компьютерная программа при выполнении побуждает процессор выполнять следующие операции:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-ой подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-той подборки сейсмических трасс с (М+1)-ой подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-той подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного кона в заданном временном интервале, вычисление коэффициента корреляции (K+1)-ой подборки сейсмических трасс с K-той подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-ой подборки сейсмических трасс; и
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-ой подборки сейсмических трасс с N-ной подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-ой подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
Кроме того, в варианте осуществления компьютерная программа при выполнении побуждает процессор выполнять остаточную коррекцию за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
Кроме того, в варианте осуществления компьютерная программа при выполнении побуждает процессор использовать нижеследующую формулу для суммирования подборок сейсмических трасс, чтобы получать модельную трассу:
где m1 и m2 обозначают порядковые номера подборок сейсмических трасс и aij обозначает значение амплитуды j-той выборочной точки из i-той трассы; m1, m2, i и j являются положительными целыми числами, 1≤m1<m2≤S, i∈[m1, m2].
Кроме того, в варианте осуществления компьютерная программа при выполнении побуждает процессор использовать нижеследующую формулу для вычисления коэффициента корреляции модельной трассы с подборкой сейсмических трасс:
где обозначает коэффициент корреляции x-вой трассы с модельной трассой, aj обозначает значение амплитуды j-той выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-той выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
Кроме того, в варианте осуществления компьютерная программа при выполнении побуждает процессор выполнять процесс выравнивания энергии амплитуд каждой подборки сейсмических трасс перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс, после процесса коррекции за нормальное приращение времени для получения модельной трассы.
Кроме того, в варианте осуществления при побуждении компьютерной программой при выполнении процессор выполняет следующие операции:
подсчет энергии Е0 амплитуд в заданном интервале временного окна:
где Q обозначает количество подсчитанных трасс, Р обозначает количество выборочных точек и aj обозначает значение амплитуды;
вычисление энергии Е1 амплитуд для каждой подборки сейсмических трасс:
вычисление коэффициента k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1;
выполнение процесса выравнивания амплитуд для каждой подборки сейсмических трасс при использовании коэффициента выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс, и обозначает значение амплитуды после процесса выравнивания амплитуд.
Согласно варианту осуществления настоящего раскрытия суммируют множество подборок сейсмических трасс; вычисляют коэффициент корреляции каждой подборки сейсмических трасс с модельной трассой и выбирают K-тую подборку сейсмических трасс; вычисляют точку оптимума каждой подборки сейсмических трасс от K-той подборки сейсмических трасс в обе стороны по порядку; и выполняют остаточную коррекцию за нормальное приращение времени подборок сейсмических трасс в соответствии с точкой оптимума. Это решение объединяют с решением из предшествующего уровня техники, в котором процесс коррекции за нормальное приращение времени реализуется при анализе скоростей распространения сейсмических волн, так что события всех подборок модельных трасс выправляются в горизонтальное состояние, в результате чего исключаются погрешности, обусловленные интерполяцией, выполняемой в решении из предшествующего уровня техники и повышается точность коррекции за нормальное приращение времени до суммирования; поэтому качество результата суммирования сейсмических данных после оптимизации подборки трасс безусловно улучшается или повышается и вследствие этого можно более точно описывать структуру пласта, литологию и свойства флюида.
Специалисту в данной области техники следует понимать, что все этапы или часть этапов согласно этим вариантам осуществления может выполняться аппаратным обеспечением или аппаратным обеспечением, выдающим релевантные команды с помощью программы, которая может сохраняться на машиночитаемом носителе данных, таком как постоянное запоминающее устройство, магнитный диск или оптический диск.
Выше описаны только предпочтительные варианты осуществления настоящего раскрытия. Любые изменения, эквивалентные замены, усовершенствования и т.д., сделанные в рамках сущности и принципов настоящего раскрытия, следует считать попадающими в объем защиты настоящего раскрытия.
Claims (72)
1. Способ обработки сейсмических данных, содержащий этапы, на которых:
суммируют подборки сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ), чтобы получить модельную трассу, при этом S является целым числом;
вычисляют коэффициент корреляции каждой подборки сейсмических трасс с модельной трассой и выбирают K-ю подборку сейсмических трасс с максимальным коэффициентом корреляции;
вычисляют точку оптимума каждой подборки сейсмических трасс от K-й подборки сейсмических трасс в обе стороны по порядку; и
выполняют остаточную коррекцию за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
2. Способ по п. 1, в котором вычисление точки оптимума каждой подборки сейсмических трасс от K-й подборки сейсмических трасс в обе стороны по порядку содержит:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-й подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-й подборки сейсмических трасс с (М+1)-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-й подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-й подборки сейсмических трасс; и
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-й подборки сейсмических трасс с N-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-й подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс с нулевым удалением.
3. Способ по п. 2, в котором выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума содержит выполнение остаточной коррекции за нормальное приращение времени S подборок сейсмических трасс при использовании вычисленных величин остаточных поправок за нормальное приращение времени для S подборок сейсмических трасс.
4. Способ по п. 1, в котором суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы содержит:
использование нижеследующей формулы для суммирования подборок сейсмических трасс для получения модельной трассы:
где m1 и m2 обозначают порядковые номера подборок сейсмических трасс и aij обозначает значение амплитуды j-й выборочной точки из i-й трассы; m1, m2, i и j являются положительными целыми числами, 1≤m1<m2≤S, i∈[m1, m2].
5. Способ по п. 1, в котором вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой содержит:
где обозначает коэффициент корреляции x-й трассы с модельной трассой, aj обозначает значение амплитуды j-й выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-й выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
6. Способ по п. 1, где перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы способ также содержит выполнение процесса выравнивания энергий амплитуд каждой подборки сейсмических трасс.
7. Способ по п. 6, в котором выполнение процесса выравнивания энергий амплитуд содержит:
подсчет энергии Е0 амплитуд в заданном интервале временного окна:
где Q обозначает количество подсчитанных трасс, Р обозначает количество выборочных точек и aj обозначает значение амплитуды;
вычисление энергии Е1 амплитуд для каждой подборки сейсмических трасс:
вычисление коэффициента k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1;
выполнение процесса выравнивания амплитуд для каждой подборки сейсмических трасс при использовании коэффициента выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс, и обозначает значение амплитуды после процесса выравнивания амплитуд.
8. Устройство для обработки сейсмических данных, содержащее запоминающее устройство, процессор и компьютерную программу, сохраняемую в запоминающем устройстве и выполняемую процессором, при этом процессор выполняет следующие операции при выполнении компьютерной программы:
суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом;
вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-й подборки сейсмических трасс с максимальным коэффициентом корреляции;
вычисление точки оптимума каждой подборки сейсмических трасс от K-й подборки сейсмических трасс в обе стороны по порядку;
выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
9. Устройство для обработки сейсмических данных по п. 8, в котором процессор также выполняет следующие операции при выполнении компьютерной программы:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-й подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-й подборки сейсмических трасс с (М+1)-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-й подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-й подборки сейсмических трасс; и
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-й подборки сейсмических трасс с N-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-й подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс с нулевым удалением.
где обозначает коэффициент корреляции x-й трассы с модельной трассой, aj обозначает значение амплитуды j-й выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-й выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
11. Устройство для обработки сейсмических данных по п. 8, в котором процессор при выполнении компьютерной программы также выполняет процесс выравнивания энергий амплитуд каждой подборки сейсмических трасс перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы.
12. Устройство для обработки сейсмических данных по п. 11, в котором процессор также выполняет следующие операции при выполнении компьютерной программы:
подсчет энергии Е0 амплитуд в заданном интервале временного окна:
где Q обозначает количество подсчитанных трасс, Р обозначает количество выборочных точек и aj обозначает значение амплитуды;
вычисление энергии Е1 амплитуд для каждой подборки сейсмических трасс:
вычисление коэффициента k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1;
выполнение процесса выравнивания амплитуд для каждой подборки сейсмических трасс при использовании коэффициента выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс, и обозначает значение амплитуды после процесса выравнивания амплитуд.
13. Машиночитаемый носитель данных, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять следующие операции:
суммирование подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени (НПВ) для получения модельной трассы, при этом S является целым числом;
вычисление коэффициента корреляции каждой подборки сейсмических трасс с модельной трассой и выбор K-й подборки сейсмических трасс с максимальным коэффициентом корреляции;
вычисление точки оптимума каждой подборки сейсмических трасс от K-й подборки сейсмических трасс в обе стороны по порядку;
выполнение остаточной коррекции за нормальное приращение времени подборок сейсмических трасс в соответствии с точками оптимума.
14. Машиночитаемый носитель данных по п. 13, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять следующие операции:
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K-1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K-1)-й подборки сейсмических трасс, где K>1 и K является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции М-й подборки сейсмических трасс с (М+1)-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для М-й подборки сейсмических трасс, где K-2≥M≥1 и М является целым числом;
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (K+1)-й подборки сейсмических трасс с K-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (K+1)-й подборки сейсмических трасс; и
скольжение временного окна в заданном временном интервале, вычисление коэффициента корреляции (N+1)-й подборки сейсмических трасс с N-й подборкой сейсмических трасс, выбор временного окна в качестве точки оптимума, когда коэффициент корреляции является максимальным, и получение величины временного сдвига, соответствующей точке оптимума, в качестве величины остаточной поправки за нормальное приращение времени для (N+1)-й подборки сейсмических трасс, где S-1≥N≥K+1 и N является целым числом, и при этом величина временного сдвига, соответствующая точке оптимума, является значением разности между временем отражения сейсмической волны из подборки сейсмических трасс во временном окне и временем отражения сейсмической волны из подборки сейсмических трасс при нулевом удалении.
где обозначает коэффициент корреляции x-й трассы с модельной трассой, aj обозначает значение амплитуды j-й выборочной точки из модельной трассы , Aa обозначает среднее значение амплитуд модельной трассы , bj обозначает значение амплитуды j-й выборочной точки из подборки сейсмических трасс, Ba обозначает среднее значение амплитуд подборки сейсмических трасс и n обозначает суммарное количество выборочных точек в подборках сейсмических трасс, где n, x и j являются положительными целыми числами и 1≤x≤S.
16. Машиночитаемый носитель данных по п. 13, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять процесс выравнивания энергий амплитуд каждой подборки сейсмических трасс перед суммированием подборок сейсмических трасс в заданном интервале из числа S подборок сейсмических трасс после процесса коррекции за нормальное приращение времени для получения модельной трассы.
17. Машиночитаемый носитель данных по п. 16, при этом на машиночитаемом носителе данных сохраняется компьютерная программа, которая при выполнении побуждает процессор выполнять следующие операции:
подсчет энергии Е0 амплитуд в заданном интервале временного окна:
где Q обозначает количество подсчитанных трасс, Р обозначает количество выборочных точек и aj обозначает значение амплитуды;
вычисление энергии Е1 амплитуд для каждой подборки сейсмических трасс:
вычисление коэффициента k выравнивания амплитуд для каждой подборки сейсмических трасс: k=E0/E1;
выполнение процесса выравнивания амплитуд для каждой подборки сейсмических трасс при использовании коэффициента выравнивания амплитуд, где обозначает последовательность, составленную из значений амплитуд каждой подборки сейсмических трасс, и обозначает значение амплитуды после процесса выравнивания амплитуд.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610204476.6 | 2016-04-01 | ||
CN201610204476.6A CN106226818B (zh) | 2016-04-01 | 2016-04-01 | 地震数据处理方法和装置 |
PCT/CN2017/078498 WO2017167191A1 (zh) | 2016-04-01 | 2017-03-29 | 地震数据处理方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2694621C1 true RU2694621C1 (ru) | 2019-07-16 |
Family
ID=57519258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2018105481A RU2694621C1 (ru) | 2016-04-01 | 2017-03-29 | Способ и устройство для обработки сейсмических данных |
Country Status (5)
Country | Link |
---|---|
US (1) | US10705238B2 (ru) |
EP (1) | EP3324216B1 (ru) |
CN (1) | CN106226818B (ru) |
RU (1) | RU2694621C1 (ru) |
WO (1) | WO2017167191A1 (ru) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106226818B (zh) | 2016-04-01 | 2018-05-04 | 中国石油天然气股份有限公司 | 地震数据处理方法和装置 |
CN107607993B (zh) * | 2017-09-07 | 2019-05-31 | 中国石油大学(北京) | 一种确定叠加速度的方法、装置及系统 |
CN107843919B (zh) * | 2017-10-27 | 2020-02-07 | 中国石油天然气集团有限公司 | 一种微地震监测数据振幅均衡方法 |
CN110609326A (zh) * | 2018-06-15 | 2019-12-24 | 中国石油化工股份有限公司 | 地震道集自动校平方法及系统 |
CN110618457B (zh) * | 2018-06-20 | 2021-02-26 | 中国石油化工股份有限公司 | 一种可控震源地震数据的联合去噪方法及其装置 |
CN110824566B (zh) * | 2018-08-10 | 2021-09-28 | 中国石油天然气股份有限公司 | 地震属性融合方法、装置及存储介质 |
CN112180446A (zh) * | 2019-07-04 | 2021-01-05 | 中国石油天然气集团有限公司 | 三维地震数据叠前道集叠加方法及装置 |
CN112230285B (zh) * | 2019-07-15 | 2023-09-26 | 中国石油天然气集团有限公司 | 地震数据激发延迟时间的校正方法及校正装置 |
CN112305593B (zh) * | 2019-07-23 | 2023-07-25 | 中国石油天然气集团有限公司 | 异常初至波识别方法及装置 |
CN112649873B (zh) * | 2019-10-10 | 2024-04-09 | 中国石油化工股份有限公司 | 地震数据分频动校正处理方法及系统 |
CN112666601B (zh) * | 2019-10-15 | 2024-03-22 | 中国石油化工股份有限公司 | 一种地震数据振幅拟合有效时窗的方法及其系统 |
CN112882102A (zh) * | 2019-11-29 | 2021-06-01 | 中国石油天然气集团有限公司 | 地震噪声衰减方法及装置 |
CN114002741B (zh) * | 2020-07-28 | 2024-01-30 | 中国石油天然气股份有限公司 | 叠前深度偏移方法及装置、计算机可读存储介质 |
CN114114418A (zh) * | 2020-08-28 | 2022-03-01 | 中国石油化工股份有限公司 | 地震信号增益方法、装置、计算机设备和存储介质 |
CN114325833A (zh) * | 2022-01-05 | 2022-04-12 | 中国海洋大学 | 一种地震记录分频带野值压制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4203161A (en) * | 1972-03-01 | 1980-05-13 | Texaco Inc. | Method of enhancing common depth point seismic data |
US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
CN102540252A (zh) * | 2011-12-15 | 2012-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于互相关的高精度中值叠加方法 |
CN104459764A (zh) * | 2014-11-20 | 2015-03-25 | 天津大学 | 一种适用于实验室研究多点地震动的连动式电磁击锤 |
WO2015106065A1 (en) * | 2014-01-10 | 2015-07-16 | Cgg Services (U.S.) Inc. | Device and method for mitigating cycle-skipping in full waveform inversion |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9726928D0 (en) * | 1997-12-19 | 1998-02-18 | Geco Prakla Uk Ltd | Method of stacking seismic signals |
CN1129008C (zh) * | 2001-06-29 | 2003-11-26 | 中国石油天然气股份有限公司 | 一种模型约束动静校正方法 |
US20030187583A1 (en) * | 2002-04-01 | 2003-10-02 | Martin Federico D. | Method and apparatus for resolving shear wave seismic data |
US7496452B2 (en) * | 2007-04-13 | 2009-02-24 | Westerngeco L.L.C. | Correction for errors caused by variation in water conditions |
US9075162B2 (en) * | 2011-11-10 | 2015-07-07 | Pgs Geophysical As | Method and system for separating seismic sources in marine simultaneous shooting acquisition |
CN102565857B (zh) * | 2011-12-16 | 2014-09-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 自动剩余动校正方法 |
CN102879821B (zh) * | 2012-09-26 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种针对地震叠前道集的同相轴精细拉平处理方法 |
US9857490B2 (en) * | 2013-12-30 | 2018-01-02 | Pgs Geophysical As | Methods and systems for optimizing generation of seismic images |
CN104181588B (zh) * | 2014-08-15 | 2016-08-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种构建高品质模型道的方法 |
CN104459794B (zh) * | 2014-12-01 | 2017-02-22 | 郑鸿明 | 共反射点道集时变时间差值的校正方法及装置 |
CN106226818B (zh) * | 2016-04-01 | 2018-05-04 | 中国石油天然气股份有限公司 | 地震数据处理方法和装置 |
-
2016
- 2016-04-01 CN CN201610204476.6A patent/CN106226818B/zh active Active
-
2017
- 2017-03-29 RU RU2018105481A patent/RU2694621C1/ru active
- 2017-03-29 WO PCT/CN2017/078498 patent/WO2017167191A1/zh unknown
- 2017-03-29 EP EP17773221.1A patent/EP3324216B1/en active Active
-
2018
- 2018-01-29 US US15/882,779 patent/US10705238B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4203161A (en) * | 1972-03-01 | 1980-05-13 | Texaco Inc. | Method of enhancing common depth point seismic data |
US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
CN102540252A (zh) * | 2011-12-15 | 2012-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于互相关的高精度中值叠加方法 |
WO2015106065A1 (en) * | 2014-01-10 | 2015-07-16 | Cgg Services (U.S.) Inc. | Device and method for mitigating cycle-skipping in full waveform inversion |
CN104459764A (zh) * | 2014-11-20 | 2015-03-25 | 天津大学 | 一种适用于实验室研究多点地震动的连动式电磁击锤 |
Also Published As
Publication number | Publication date |
---|---|
WO2017167191A1 (zh) | 2017-10-05 |
US10705238B2 (en) | 2020-07-07 |
EP3324216A1 (en) | 2018-05-23 |
EP3324216B1 (en) | 2020-08-05 |
CN106226818B (zh) | 2018-05-04 |
EP3324216A4 (en) | 2018-10-17 |
CN106226818A (zh) | 2016-12-14 |
US20180149764A1 (en) | 2018-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2694621C1 (ru) | Способ и устройство для обработки сейсмических данных | |
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
US9103935B2 (en) | Method of first arrival picking of seismic refraction data | |
CA2375495C (en) | Method and apparatus for doppler smear correction in marine seismology measurements | |
CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
MX2014003060A (es) | Sistemas y metodos de filtracion de dominio de frecuencia y discriminacion de dominio de espacio-tiempo de datos sismicos. | |
CN107272062B (zh) | 一种数据驱动的地下介质q场估计方法 | |
US20120158382A1 (en) | process for characterising the evolution of an oil or gas reservoir over time | |
RU2631407C1 (ru) | Способ и устройство для обработки сейсмических сигналов | |
US20160054465A1 (en) | Seismic data spectrum restoring and broadening | |
US7609585B2 (en) | Method for sub-salt migration velocity analysis | |
CN104181587A (zh) | 一种地震数据振幅谱的相干值获取方法及系统 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
WO2005026776A1 (en) | Wide-offset-range pre-stack depth migration method for seismic exploration | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN114966837A (zh) | 基于卷积型w2距离目标函数的全波形反演方法及装置 | |
CN109143366A (zh) | 近地表初至层析静校正方法及计算机可读存储介质 | |
CN108693560A (zh) | 一种基于互相关道的散射波成像方法及系统 | |
CN104570115A (zh) | 一种面波衰减方法及装置 | |
CN109143345B (zh) | 基于模拟退火的品质因子q非线性反演方法及系统 | |
US10613246B2 (en) | Joint estimation of electromagnetic earth responses and ambient noise | |
CN112305615B (zh) | 一种地震资料角度域共成像点道集提取方法及系统 | |
CN115170428A (zh) | 一种声波远探测成像图的降噪方法 | |
CN110824555B (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
CN112462428B (zh) | 一种多分量地震资料偏移成像方法及系统 |