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

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

Info

Publication number
EA011314B1
EA011314B1 EA200701070A EA200701070A EA011314B1 EA 011314 B1 EA011314 B1 EA 011314B1 EA 200701070 A EA200701070 A EA 200701070A EA 200701070 A EA200701070 A EA 200701070A EA 011314 B1 EA011314 B1 EA 011314B1
Authority
EA
Eurasian Patent Office
Prior art keywords
reflection
amplitude
dependence
distance
azimuthal
Prior art date
Application number
EA200701070A
Other languages
English (en)
Other versions
EA200701070A1 (ru
Inventor
Вендзи Донг
Чих-Пинг Лю
Original Assignee
Эксонмобил Апстрим Рисерч Компани
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Эксонмобил Апстрим Рисерч Компани filed Critical Эксонмобил Апстрим Рисерч Компани
Publication of EA200701070A1 publication Critical patent/EA200701070A1/ru
Publication of EA011314B1 publication Critical patent/EA011314B1/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32Transforming one recording into another or one representation into another
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/52Move-out correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]

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)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

В настоящем изобретении описаны способы количественной оценки и ослабления вызванных наклоном эффектов азимутального изменения амплитуды отражения от удаления при сейсмическом обнаружении трещин путем использования анализа азимутальной зависимости амплитуды отражения от удаления с точным учетом коррекции за геометрическое расхождение и азимутальной зависимости угла отражения. Решения предложены для трех случаев: (1) наклонные изотропные коллекторы; (2) анизотропные коллекторы с расположенными параллельно трещинами, ориентированными в произвольном направлении; и (3) анизотропные коллекторы, в которых расположенные параллельно вертикальные трещины перпендикулярны к направлению падения.

Description

В общем настоящее изобретение относится к области геофизических исследований, включая оконтуривание коллекторов геофизическим способом. В частности, изобретение представляет собой способ обработки сейсмических данных, предназначенный для получения присущих пласту свойств породы, флюида и анизотропии трещин путем анализа изменения амплитуды сейсмических волн с удалением и азимута при наличии угла падения пласта.
Уровень техники
При сейсмической разведке на углеводороды изменение амплитуды отражения с удалением (АУО) объясняется изменением амплитуды Р-волны по мере увеличения расстояния между точкой взрыва и приемником или угла отражения. Когда пласты пород являются плоскими и изотропными (изотропия означает, что упругие свойства среды остаются неизменными независимо от направлений, по которым они исследуются), амплитуда отражения при фиксированном удалении не изменяется с азимутом (углом направления линии между точкой взрыва и приемником, отсчитываемым от направления на север). Эта ситуация отражена на фиг. 1а, на которой показана траектория сейсмического луча, начинающаяся у источника (в точке сейсмического взрыва) 8, проходящая вниз с отражением от поверхности пласта в точке Ό и затем продолжающаяся в обратном направлении до приемника К, осуществляющего обнаружение. Отраженный луч имеет угол θ с нормалью к отражающей поверхности или к вертикали, поскольку отражающий горизонт предполагается плоским. На фиг. 1Ь показаны сейсмические данные, зарегистрированные с помощью приемника К для нескольких различных значений расстояния 11 между источником и приемником (называемого удалением), при этом каждая отраженная волна пришла из одной и той же отражающей точки Ό. Посредством приемника измеряется амплитуда сейсмической волны как функция времени двойного пробега сейсмической волны, распространяющейся вниз до Ό от 8 и затем обратно к К (фиг. 1а). Время пробега связано с глубиной нахождения отражающей точки через геометрию траектории луча и скорость сейсмической волны. Отражающий горизонт 12 из фиг. 1а показан на фиг. 1Ь для различных времен пробега (расстояния 8-Э-К), зависящих от удаления. Причина состоит в том, что данные еще не скорректированы с помощью стандартной процедуры обработки сейсмических данных, называемой введением поправки за нормальное приращение времени.
На фиг. 2а пояснено вышеупомянутое определение азимута, то есть угла ψ между направлением на север и линией между источником 8 и приемником К. На фиг. 2Ь показана плоская характеристика амплитуды отражения в зависимости от азимута, которая означает, что среда является азимутально изотропной. Однако амплитуда отражения становится зависимой от азимута, когда пласт под отражающей поверхностью является азимутально анизотропным вследствие, например, наличия расположенных параллельно вертикальных трещин. Тогда зависимость амплитуды отражения от азимута может быть до некоторой степени похожей на показанную на фиг. 2с. В таком случае изменение амплитуды называют азимутальным изменением амплитуды отражения с удалением. Геофизики осуществляют анализ зависимости амплитуды отражения от удаления, чтобы получать свойства породы и флюида в коллекторном пласте ниже границы отражения. (См., например, 8ткк С.С. апб С1б1оте Р.М., \Усщ1Исб 81асктд Гог госк ргореДу езйтайоп апб бсЮсбоп о! даз, Сеоркуз. Ргозр., Еиг. Аззп. 8со5ст Епд., 35, 993-1014 (1987); Ки!кегГогб 8.К. апб ^йкатз К.Н., АтрШибе-уегаш-ойке! уапакощ ίη даз запбз, Сеоркуысз 54, 680-688 (1989); Саз!адпа 1.Р., 8теап Н.\У. апб Еокег Ό.Ι., Етатетеотк йог АУО дгаб1еп! апб 1п!егсер! т!етрге1а1юп, Сеоркуысз 63, 948-956 (1998)). Как раскрыто, например, в работах Сотдап Ό., Тке еГГес! о! ахппЩка1 апщойору оп !ке уапа!юп о! теДескуку текк оГГзе!: ^откзкор оп 8е15т1с АпщоДору, 8ос. Ехр1. Сеоркуз., 41\У8А. 1645 (1990) и Кидег А., Уапа!юп о! Р-теауе теДесДуку текк оГГзе! апб ахипЩк ш апщоДорк теб1а, СеоркуДсз 63, 935-947 (1998).) Данные азимутальной зависимости амплитуды отражения от удаления могут быть инвертированы для получения свойств породы и флюида, а также трещин, когда предполагается, что в коллекторном пласте существуют трещины.
Кроме того, азимутальное изменение амплитуды отражения с удалением может наблюдаться, когда пласт является изотропным, но ограничен сверху наклонной поверхностью границы. Это явление может быть проиллюстрировано численным примером. Как показано на фиг. 3а (угловая шкала в градусах), в этом синтетическом примере сейсмические линии были взрывами по многим азимутальным направлениям в пределах простой модели геологической среды из фиг. 3Ь. На этой фигуре показаны траектории лучей от пункта 31 взрыва, обнаруживаемых многими различными приемниками в 32. Поскольку сейсмические волны расширяются или расходятся по мере того, как они распространяются, необходимы поправки на геометрическое расхождение для восстановления амплитуд отражения, чтобы представлять истинную интенсивность отражения на отражающей границе. Однако обработкой из известного уровня техники амплитуды отражения от наклонного отражающего горизонта не восстанавливаются должным образом, поскольку предполагается, что пласты пород, такие, как 33 на фиг. 3Ь, являются плоскими. Кроме того, в ней пренебрегается коррекцией амплитуд для компенсации зависимых от азимута углов
- 1 011314 отражений.
Результирующая амплитуда в дополнение к тому, что она является зависимой от удаления (или угла падения), также становится азимутально зависимой (фиг. 3с). Этим объясняется то, что амплитуда отражения при фиксированном удалении является зависимой от азимута, когда отражающий горизонт является наклонным. Изменение амплитуды такого типа может быть названо вызванным наклоном азимутальным изменением амплитуды отражения с удалением для отличия его от азимутального изменения амплитуды отражения с удалением, вызванного анизотропией из-за наличия, например, расположенных параллельно вертикальных трещин. Вызванное наклоном азимутальное изменение амплитуды отражения с удалением искажает амплитуду сейсмической волны не только тогда, когда пласты являются изотропными, но также и тогда, когда пласты являются анизотропными.
Обычно при анализе зависимости изменения амплитуды отражения от удаления не осуществляют количественную оценку и снижение вызванного наклоном азимутального изменения отражения с удалением. Это обусловлено не только небрежностью, но также и тем, что данные азимута утрачиваются при выполнении процедур обработки данных (например, при коррекции на приращение времени, обусловленное наклоном отражающей границы, и миграции). Для подавления эффекта наклона при инверсии азимутальной зависимости амплитуды отражения от удаления был предложен способ секторной миграции во временной области. В этом способе сначала разделяют множество трехмерных наборов данных на несколько подмножеств, каждое из которых содержит трассы в угловом диапазоне вокруг выбранного азимутального угла. Подмножества дополнительно обрабатывают путем осуществления сейсмической миграции до выполнения анализа азимутальной зависимости амплитуды отражения от удаления. (Сейсмическая миграция представляет собой способ обработки данных, при осуществлении которого отражения располагают в соответствии с их действительными положениями. Сейсмическую миграцию используют по причине изменяющихся скоростей сейсмических волн и наклонных горизонтов.) Однако этим способом невозможно полностью решить проблему, поскольку отражающие точки перемещаются на различную величину из сектора в сектор вследствие изменения кажущегося наклона. Результаты, получаемые способом временной секторной миграции, могут быть неоднозначными.
Другой используемый способ представляет собой инверсию азимутальной зависимости амплитуды отражения от удаления в трехмерных сейсмических данных без осуществления сейсмической миграции. В этом способе пласты породы предполагают плоскими, и поэтому применяемой коррекцией на геометрическое расхождение не компенсируется вызванное наклоном азимутальное изменение амплитуды отражения с удалением. Этот способ является предпочтительным по сравнению со способом секторной инверсии, поскольку при инверсии для повышения статистической мощности в случае оценивания атрибутов трещин в нем используют информацию об амплитуде, удалении и азимуте одновременно из всех трасс. Однако результат инверсии также может быть неоднозначным, поскольку из данных не исключается вызванное наклоном азимутальное изменение амплитуды отражения с удалением.
Поскольку газонефтяные коллекторы часто не являются плоскими, способ количественной оценки и снижения, вызванного наклоном азимутального изменения амплитуды отражения с удалением во время анализа амплитуд сейсмических волн, может быть полезным при наличии анизотропии или без нее. Для повышения надежности результатов анализа зависимости амплитуды отражения от удаления и азимутальной зависимости амплитуды отражения от удаления являются полезными аналитические соотношения, связывающие амплитуду, угол падения пласта и азимутальные углы, применение которых улучшает понимание сущности процессов, количественную оценку степени ударного сейсмического воздействия и расчетные методики, предназначенные для ослабления эффекта. Настоящее изобретение направлено на удовлетворение этих потребностей.
Сущность изобретения
Согласно одному варианту осуществления описан способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых на основании результатов сейсмического исследования наклонного коллектора, лежащего под изотропной перекрывающей породой. При выполнении способа (а) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (Ь) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных.
Согласно второму варианту осуществления описан способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного изотропного коллектора, лежащего под изотропной перекрывающей породой. Способ содержит этапы, при выполнении которых: (а) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования объема данных с четырьмя атрибутами А, В, В' и ψα, соответствующими взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка, и углового коэффициен
- 2 011314 та зависимости амплитуды отражения от удаления, соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ,, является азимутом расположенных параллельно трещин; (Ь) вычисляют углы φ падения для верхней поверхности коллектора и умножают со§р> на А для получения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зависимости амплитуды отражения от удаления для фактического наклонного коллектора; (с) вычисляют скомпенсированный за наклон атрибут В0 углового коэффициента зависимости амплитуды отражения от удаления для фактического наклонного коллектора из В0=(В+В') со5<р; и (ά) вычисляют азимут ψ0 падения фактического наклонного коллектора из ψ0α.
Согласно третьему варианту осуществления описан способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой. Этот способ содержит этапы, при выполнении которых:
(a) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема трехмерных сейсмических данных;
(b) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из объема трехмерных сейсмических данных, при этом масштабный коэффициент определяют так, чтобы сделать амплитуды из объема трехмерных сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
(c) применяют масштабный коэффициент к трехмерным сейсмическим данным последовательным выборочным способом и (ά) инвертируют уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из объема трехмерных масштабированных сейсмических данных для формирования объемов данных с четырьмя скомпенсированными за наклон атрибутами А, В, В' и ψ£ для фактического наклонного анизотропного коллектора, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ £ является азимутом расположенных параллельно трещин.
Согласно четвертому варианту осуществления описан способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой, при этом анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно к азимуту падения наклонного анизотропного коллектора. Способ содержит этапы, при выполнении которых (а) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно вертикальных трещин; (Ь) вычисляют углы φ падения для верхней поверхности наклонного анизотропного коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора путем умножения А на со»р; (с) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности фактического наклонного анизотропного коллектора и умножают сумму на со»р; (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора; (е) вычитают В0 из (В+В')со§р> для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута Βι углового коэффициента для фактического наклонного анизотропного коллектора; и (£) выполняют итерацию этапов с (а) по (е) с оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут заданный критерий сходимости.
Согласно пятому варианту осуществления описан способ оценивания магнитуды, вызванной наклоном погрешности в вычисленных значениях атрибутов отсекаемого отрезка, и углового коэффициента зависимости амплитуды отражения от удаления, извлекаемых из сейсмических данных в результате трехмерного сейсмического исследования наклонного изотропного отражающего горизонта, лежащего под изотропной покрывающей породой. Способ содержит этапы, при выполнении которых (а) вычисля
- 3 011314 ют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения с удалением, из
Магнитуда = Тап2 φ είη2 Θ1 , где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, соответственно, φ - угол падения наклонного изотропного отражающего горизонта и θ1 -угол отражения, то есть характеристика волн в сейсмических данных, отраженных от поверхности наклонного изотропного отражающего горизонта, в предположении, что наклонный изотропный отражающий горизонт является плоским; (Ь) вычисляют величину базовой амплитуды, характерную для поверхности плоского изотропного отражающего горизонта, из
Базовая амплитуда = Α00δίη2θ1; и (с) вычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
Согласно шестому варианту осуществления описан способ компенсации вызванных наклоном погрешностей. Способ включает в себя получение сейсмических данных в результате сейсмического исследования наклонного пласта; и компенсацию за расхождение энергии и угол отражения в сейсмических данных для наклонного пласта. Следует отметить, что этот способ может быть реализуемым компьютером способом, который осуществляется на компьютере, и при этом отображается или формируется отчет, основанный на результатах описанного способа.
Краткое описание чертежей
Настоящее изобретение и его преимущества станут более понятными при обращении к нижеследующему подробному описанию и сопровождающим чертежам, на которых фиг. 1а - иллюстрация взаимосвязи между траекторией сейсмического луча, точкой взрыва, прием ником, удалением и углом отражения в горизонтально-слоистом пласте;
фиг. 1Ь - иллюстрация сейсмических отражений от границы раздела слоев, изображенных в виде волн, при этом амплитуда отражения определена по вершине волны для каждого удаления;
фиг. 2а - иллюстрация определения понятия азимута линии между точкой взрыва и приемником; фиг. 2Ь - диаграмма, иллюстрирующая, что в случае границы раздела между плоскими изотропными средами амплитуда отражения при фиксированном удалении не изменяется с азимутом;
фиг. 2с - диаграмма, иллюстрирующая, что амплитуда отражения становится зависимой от азимута, когда пласт под отражающей поверхностью является азимутально анизотропным;
фиг. 3 а - вид геометрии возбуждения волн по азимутам;
фиг. 3Ь - иллюстрация числовой модели пласта и примера траекторий лучей вдоль линии взрыва; фиг. 3 с - диаграмма, иллюстрирующая, что амплитуда отражения при фиксированном удалении является зависимой от азимута в случае модели пласта из фиг. 3Ь;
фиг. 4а - иллюстрация изотропного пласта, покрывающего наклонный изотропный пласт;
фиг. 4Ь - иллюстрация изотропного пласта, покрывающего наклонный анизотропный пласт с трещинами, перпендикулярными к наклонной поверхности;
фиг. 5 - иллюстрация линейной зависимости между коэффициентом отражения и квадратом синуса угла отражения;
фиг. 6 - иллюстрация двумерной зависимости между траекторией сейсмического луча, точкой взрыва, приемником, удалением и углом отражения для случая наклонной отражающей поверхности;
фиг. 7 - трехмерная перспективная диаграмма для целей преобразования геометрии траектории луча от случая плоской отражающей границы к случаю наклонной отражающей границы;
фиг. 8 - блок-схема последовательности операций, отражающая основные этапы способа осуществ ления компенсации для наклонного изотропного коллектора согласно одному варианту осуществления настоящего изобретения (случай 1);
фиг. 9 - блок-схема последовательности операций, отражающая основные этапы способа осуществ ления компенсации для наклонного анизотропного коллектора с ориентированными вертикальными трещинами согласно одному варианту осуществления настоящего изобретения (случай 2);
фиг. 10 - блок-схема последовательности операций, отражающая основные этапы способа осущест вления компенсации для наклонного анизотропного коллектора с вертикальными трещинами, располо женными перпендикулярно направлению наклона, согласно одному варианту осуществления настоящего изобретения (случай 3);
фиг. 11а-с - три разреза по параметрам вдоль сейсмической линии, полученные в результате инверсии азимутальной зависимости амплитуды отражения от удаления; и фиг. 12 - иллюстрация процедуры сложения двух атрибутов для получения присущих пласту значений атрибута Β0 углового коэффициента зависимости амплитуды отражения от удаления вдоль сейсмической линии в кубе трехмерных сейсмических данных.
Изобретение будет описано применительно к его предпочтительным вариантам осуществления. Однако в той части, в какой нижеследующее подробное описание является характерным для конкретного
- 4 011314 варианта осуществления или конкретного использования изобретения, оно подразумевается только иллюстративным и не должно толковаться как ограничивающее объем изобретения. Напротив, оно предполагается охватывающим все варианты, модификации и эквиваленты, которые могут быть заключены в рамках сущности и объема изобретения, определенного прилагаемой формулой изобретения.
Подробное описание предпочтительных вариантов осуществления
Настоящее изобретение представляет собой способ количественной оценки и ослабления эффектов наклона при анализе зависимости амплитуды отражения от удаления и инверсии азимутальной зависимости амплитуды отражения от удаления. В способе согласно изобретению в качестве базовой точки использован ранее описанный способ инверсии азимутальной зависимости амплитуды отражения от удаления на основе общей средней точки. Уравнения для описания свойств азимутальной зависимости амплитуды отражения от удаления при наличии угла падения пласта выводятся для трех случаев: (1) случая изотропного слоя 18 (то есть покрывающей породы) поверх наклонного изотропного слоя 19 коллектора, показанного на фиг. 4а, (2) случая изотропного слоя 14 поверх наклонного и анизотропного коллектора 15 с ориентированными вертикальными трещинами 16, показанного на фиг. 4Ь, и (3) особого случая, дополнительного к случаю (2), когда азимут трещин является перпендикулярным к азимуту падения трещиноватого пласта (также показанного на фиг. 4Ь). Используемый в настоящей заявке термин «вертикальные трещины» означает трещины, перпендикулярные к наклонной отражающей поверхности над ними. Азимут трещин определяется как угол между (а) линией, образуемой пересечением плоскости трещин с горизонтальной плоскостью, и (Ь) опорным направлением, которое обычно является направлением на географический север. Азимут падения представляет собой угол между направлением падения пласта и направлением на север. Представления из этих уравнений используют для определения степени влияния эффекта наклона и для разработки простых и управляемых данными способов, предназначенных для исключения эффекта наклона из инверсии азимутальной зависимости амплитуды отражения от удаления.
Случай 1: отсутствие трещин и изотропная покрывающая порода поверх наклонного изотропного коллектора (фиг. 4а)
1. Вывод формулы
Для случая, когда пласты породы являются изотропными, 8йиеу (А ытрШтсайоп οί 111е ΖοορρπΙζсс.|иа1юп5. СеорйуШсх 50 (4), 609-614 (1985)) установил, что существует линейная зависимость между коэффициентом (К) отражения и углом (θ) отражения при угле отражения меньше, чем 30°:
^Αο+Βοδίη2#. (1)
Как показано на фиг. 5, А0 и В0 представляют собой соответственно отсекаемый отрезок и угловой коэффициент аппроксимирующей методом наименьших квадратов линии 51, проведенной на перестроенном из одной системы координат в другую графике зависимости К от §ίη2θ, полученной из сейсмических данных. Они являются параметрами зависимости амплитуды от удаления (называемыми атрибутами), используемыми для получения литологических свойств и свойств флюида нижележащего коллекторного пласта (см., например, статью Сайадпа и других, 1998). Приближенное соотношение для коэффициента К отражения и амплитуды Атр имеет вид
Άιηρ=άίν· Κ-άίν· [ Αο+,Β^ίη2#] , (2) где άίν - член расхождения, который описывает, какое количество волновой энергии сохранилось после рассеивания за счет расширения волнового фронта во время распространения волны. Величина расхождения близка и может быть аппроксимирована к обратной величине пройденного волной расстояния по траектории луча от точки 8 взрыва до приемника К через точку Ό отражения на фиг. 1а. Когда отражающая поверхность является горизонтальной (то есть не наклонной), величина άίν изменяется с расстоянием между точкой взрыва и приемником (фиг. 1а). Амплитуду отражения корректируют на геометрическое расхождение до анализа зависимости амплитуды отражения от удаления.
По установившейся практике поправку на геометрическое расхождение при анализе зависимости амплитуды отражения от удаления вычисляют в предположении, что пласты породы являются плоскими, как показанные на фиг. 1а. Однако в действительности большинство коллекторных пластов являются наклонными, и расхождение энергии волны становится зависимым от азимута вследствие того, что протяженность траектории 8-Э-К луча в случае наклонного отражающего горизонта (фиг. 6) зависит от азимута. При рассмотрении наклонного пласта как плоского пласта в две величины в уравнении (2) вносится погрешность. Во-первых, как можно видеть из сравнения протяженности траектории 8-Э-К луча на фиг. 6 с протяженностью траектории 8-Э-К луча на фиг. 1а, расхождение энергии волны полностью не компенсируется. Во-вторых, угол θ отражения в уравнении (2) вычисляют, используя геометрию из фиг. 1, тогда как вычисление должно быть основано на фиг. 6. Для осуществления анализа зависимости амплитуды отражения от удаления работающим в случае наклонного пласта (в предположении плоских слоев после коррекции за геометрическое расхождение), уравнение (1) должно быть переписано для компенсации изменений в расхождении и угле отражения. В вариантах осуществления способа настоящего изобретения во всех трех случаях, описанных в настоящей заявке осуществляется коррекция сейсмических амплитуд применительно к обоим источникам погрешности.
- 5 011314
На основе геометрии траектории луча, показанной на фиг. 7, угловые величины θ и θ1 могут быть приближенно связаны с помощью следующего выражения:
είη2 Θ-—-— [1- з!п2 φ соз2($/ - 3 3ίη2 в (3) соз φ где φ - угол падения отражающего пласта 72 (то есть угол, на который поверхность пласта отклоняется от горизонтали), ψ -азимутальный угол линии (8-К) между точкой взрыва и приемником, ψ0 - азимут падения наклонного пласта (то есть направление 73 максимального падения), θ - истинный угол отражения и θ1 -срединный угол, то есть угол отражения, который может быть получен при некорректном предположении, что наклонная плоскость 72 подземного отражающего горизонта параллельна горизонтальной плоскости 71 поверхности. На фиг. 7 точка О является средней точкой между источником 8 и приемником К, Ό - точка отражения на отражающем горизонте 72 и М, расположенная по вертикали ниже О, является точкой, где может находиться точка отражения, если отражающий горизонт 72 плоский.
Для компенсации изменения при коррекции на геометрическое расхождение вследствие различия между плоским и наклонным пластами можно получить компенсирующий расхождение член άίνοοιηρ. основанный на различии в протяженностях траекторий 8-Э-К и 8-М-К лучей (фиг. 7), чтобы учесть угол падения пласта:
Г 11 άίνεοπιρ = соз φ1--Гап2 φ зтп2^ - 3ίη2 θ. ,(4) при этом результирующая амплитуда, полностью скорректированная на геометрическое расхожде ние для случая наклонных пластов, выражается с помощью следующего уравнения:
ΑΜΡ=<3ίναοιηρ· [Доо 5ίη20] .(5)
Величина АМР эквивалентна коэффициенту К отражения в предположении наличия других факторов рассеяния энергии (например, затухания волн), не влияющих на анализ зависимости амплитуды отражения от удаления. Численные исследования этих уравнений показали, что уравнения (3) и (4) являются справедливыми при углах падения пласта до около 20° и углах отражения до 30°.
Подстановка уравнений (3) и (4) в уравнение (5) дает:
Уравнением (6) описывается взаимосвязь между амплитудой, углом падения отражающего пласта и срединным углом θ1 отражения отраженной волны. Это указывает, что, когда изотропный пласт не является плоским, амплитуда отражения на наклонной границе становится зависимой от азимута. Это уравнение может быть названо уравнением вызванного наклоном азимутального изменения амплитуды отра жения с удалением.
2. Количественная оценка вызванного наклоном азимутального изменения амплитуды отражения с удалением
Из уравнения вызванного наклоном азимутального изменения амплитуды отражения с удалением вытекает способ количественной оценки влияния наклона отражающего горизонта на амплитуду отражения. После умножения обеих частей уравнения (6) на οοκφ уравнение приобретает вид:
соз φ-АМР ® Ад + [Во + (Βυ - Ав) Ьап2 φ] З1пг θί - (Ββ - Аа) 1ап2 φ соз2(((/ - зап2 = (6-1)
Левая часть уравнения может рассматриваться как масштабированная амплитуда сейсмической отраженной волны в подборке трасс с общей средней точкой. Правая часть уравнения может быть разделена на две части. Первой частью, состоящей из первого и второго членов, характеризуется коэффициент отражения на границе пласта без угла падения (см. уравнение (1)). Второй частью, состоящей из третьего члена, определяется количественно азимутальное изменение амплитуды, обусловленное углом наклона пласта. Дальнейшее рассмотрение облегчается следующими определяющими членами:
- 6 011314
Базовая амплитуда - Αο+·Βο3ίη20ί,
Изменение амплитуды = 3ίη2(^ - ψν)·(Во - ± Ао) Сап2 φ Ξίη2 . (6.2)
Поскольку 5ίη2 (ψ-1//0) изменяется от 0 до 1, диапазон изменения амплитуды находится между 0 и (Во- Ао)сап2φ 3ίη2Θ1 .
(6.3) (6.4)
Магнитуда этого изменения выражается просто как
1 , Магнитуда = [В,, - - Ао) Сап (Ρ3ίη θί .
Этой величиной определяется максимальная степень изменения амплитуды при наличии вызванного наклоном азимутального изменения амплитуды отражения с удалением. Ее значение определяется не только средней точкой и азимутальным углом, но также и различием в присущих породе свойствах на протяжении границы отражения (то есть А0 и В0).
В сущности член изменение амплитуды может считаться аномалией амплитуды отражения, обусловленной падением пласта. Это можно заметить, когда уровень амплитуды выше уровня случайного шума в сейсмических данных. Для сравнения по интенсивности этой аномалии (то есть вызванного наклоном азимутального изменения амплитуды отражения с удалением) и сейсмического шума, можно задать следующий количественный параметр:
, Магнитуда % изменения = 10Ох----------------Базовая амплитуда
Этот количественный параметр может быть вычислен для условий конкретного коллектора и выполнено количественное сравнение его с уровнем шума в сейсмических данных, если вызванное наклоном азимутальное изменение амплитуды отражения с удалением может быть обнаружено (то есть больше, чем отношение шума к сигналу) в сейсмических данных. Если вызванное наклоном азимутальное изменение амплитуды отражения с удалением может быть обнаружено, становится необходимым ослабление влияния наклона. Этот подход иллюстрируется ниже в примере.
В табл. 1 приведен пример применения уравнения (6.4) для четырех различных условий коллектора, обычных при разведке на нефть и газ. В предположении, что отношение шума к сигналу в сейсмических данных составляет меньше 5% (это отношение может быть измерено на основании сейсмических данных), в примере показано, что вызванное наклоном азимутальное изменение амплитуды отражения с удалением является значительным в случаях 1 и 2 (то есть % изменения больше, чем уровень шума), но не в случаях 3 и 4. Поэтому снижение вызванного наклоном азимутального изменения амплитуды отражения с удалением является полезным для коллекторов, относящихся случаям 1 и 2.
Таблица 1. Максимальное процентное изменение вызванного наклоном азимутального изменения амплитуды отражения с удалением для четырех условий глинистой породы/коллектора
Случай Тип коллектора Порис- тость Угол паде- ния Сре- динный угол Ао Во % изменения
1 Песча- ник 20% 15° 26’ 0,15 -0,4 9,63
2 Песчаник 25% 15’ 26’ 0,03 -0,28 16,39
3 Песча- 34% 15’ 26’ -0,1 -0,16 1,19
НИК
4 Извест- няк 20% 15’ 26° -0,111 -0,142 3,41
3. Способ исключения влияния падения пласта при анализе зависимости амплитуды отражения от удаления
Можно заметить, что уравнение (6) вызванного наклоном азимутального изменения амплитуды отражения с удалением имеет такой же вид, как и уравнение вызванного наклоном азимутального изменения амплитуды отражения с удалением, полученного Кидег (1998):
если в уравнении (7) заменить ψ,, азимутальным углом ψ0 падения. Величина ψα в уравнении Кидег представляет собой азимут, обусловленный анизотропией, например, азимут трещины, ψ представляет собой азимутальный угол линии между точкой взрыва и приемником. Параметры А и В эквивалентны
- 7 011314 отсекаемому отрезку и угловому коэффициенту в уравнении изменения амплитуды отражения с удалением, то есть А0 и В0 в уравнении (1), а В' относится к интенсивности анизотропии. Уравнение (7) является базовым уравнением, используемым при инверсии вызванной наклоном азимутальной зависимости амплитуды отражения от удаления для получения этих параметров анизотропии на основании сейсмических данных.
Выравнивающие коэффициенты подобных членов в уравнениях (6) и (7) выражаются следующими соотношениями:
А = —ί--Ао , соз φ
I 1 ϊ
В =-----[Во + (Во - - Еап φ], соз φ 2
В =-----(Во - - Ао) Сап2 φ , соз φ 2 ί^ο= ψ&.
Дальнейшее преобразование этих уравнений приводит к следующим соотношениям между параметрами, характерными для наклонного изотропного пласта, и параметрами, получаемыми в результате обычной инверсии, вызванной наклоном азимутальной зависимости амплитуды отражения от удаления в трехмерных сейсмических данных:
Αο=Α·οο3(Ρ,(8.1)
В0=(В+В')-созр,(8.2) и ιρο = (8.3)
Из уравнений (8.1)-(8.3) видно, что в случае изотропных наклонных пластов присущие пластам параметры могут быть получены путем инверсии азимутальной зависимости амплитуды отражения от удаления в сейсмических данных. Основываясь на понимании этого, можно легко получить простую, управляемую данными процедуру для ослабления влияния падения пласта при анализе зависимости амплитуды отражения от удаления. Эта процедура, сведенная в блок-схему последовательности операций на фиг. 8, может быть сформулирована следующим образом:
1) Применение инверсии азимутальной зависимости амплитуды отражения от удаления к трехмерным сейсмическим данным для формирования кубов данных с четырьмя атрибутами (то есть параметрами) А, В, В' и ψα, этап 81. Теперь, как показано уравнением (8.3), атрибут ψ,, является азимутом падения наклонного изотропного пласта.
2) Извлечение угла падения целевого отражающего горизонта из сейсмических данных после обычной обработки (например, из данных, суммированных по общей средней точке, или данных после миграции и суммирования), этап 82.
3) Как представлено уравнением (8.1), умножение косинуса угла падения (то есть οοκφ) на атрибут А, получаемый при каждом положении общей средней точки вдоль целевого отражающего горизонта, для нахождения собственного (то есть скорректированного на наклон) параметра А0 пласта, этап 83.
4) Как представлено уравнением (8.2), суммирование атрибутов В и В' при каждом положении общей средней точки вдоль целевого отражающего горизонта и умножение суммы на косинус угла падения для этих положений с целью нахождения собственного параметра В0 пласта (этап 84). Полученные в результате параметры А0 и В0 могут быть использованы для определения свойства породы и типа флюида в наклонном коллекторе с помощью известных из уровня техники процедур, выполненных, например, согласно упомянутой выше статье Сайадиа и других, 1998.
О наличии трещин в коллекторе обычно свидетельствуют пробы из скважины (например, керн, изображения пласта и т.д.), эксплуатационные характеристики скважины или расшифровка сейсмических данных. Когда трещины не являются очевидными, предшествующий простой способ может быть применен для извлечения собственных параметров А0 и В0 пласта с целью интерпретации литологии и флюида в случае плоских или наклонных коллекторов. Этот способ предпочтительно применять к пласту с углом падения меньше, чем 20°. Угол падения может быть вычислен по трехмерным сейсмическим изображениям (например, по сейсмическим данным после миграции) и использован при обработке. Когда трещины являются очевидными, может быть применен нижеследующий способ.
Случай 2: наличие трещин и изотропная покрывающая порода поверх наклонного анизотропного коллектора (фиг. 4Ь)
1. Вывод формулы
В этом случае коэффициент отражения описывается уравнением азимутального изменения амплитуды отражения с удалением, но модифицируется членом компенсации на геометрическое расхождение: ΑΜΡ^άίνσοίηρ· [А+В 3ίη2^-Β ' соз2 ( ψ- ψμ здп2#] . (9)
- 8 011314
Отметим, что угол θ относится к истинному углу отражения на наклонной поверхности (фиг. 2а). Для разделения влияния угла падения и влияния трещин уравнение (3) может быть переписано 3ίη2£=Ο· 3ίη20ί,(10) где
С = —-— [1 - 5ίη2 φ соз2(^ -.
соз φ
Подстановка уравнения (10) в уравнение (9) дает
АМР - С-άί + В ат η2 &ί + в' соз2($/ - ψ 3ΐη2 .(11)
Можно видеть, что теперь уравнение (9) выражено в членах срединного угла θ1, а величина в квадратных скобках имеет форму уравнения инверсии азимутальной зависимости амплитуды отражения от удаления, то есть уравнения (7). Кроме того, в этом уравнении С и άίνοοιηρ могут быть вычислены для каждой отраженной волны. Эти две величины могут быть объединены и определены как масштабный коэффициент 8сйт:
зсРрг =-----?-----. (12)
Οα?ϊ усотр
Подстановка уравнения (12) в уравнение (11) и перемещение масштабного коэффициента в левую часть уравнения дает
АМР-5СГ£Г = + В 31П2 θι + В' СО52(^ - ψ.] 5ΪΠ* θί . (13)
Из этого уравнения следует, что благодаря масштабированию сейсмической амплитуды с помощью масштабного коэффициента зс£(г уравнение амплитуд имеет такую же форму, как и уравнение азимутального изменения амплитуды отражения с удалением даже в случае, когда анизотропный пласт имеет угол падения. Исходя из этого понимания, может быть разработан простой способ для исключения влияния падения пласта из сейсмических данных в случае обнаружения трещин.
2. Способ исключения влияния падения пласта при анализе азимутальной зависимости амплитуды отражения от удаления
Уравнением (13) обеспечивается простой путь ослабления влияния падения в случае наклонных анизотропных пластов. Процедура, сведенная в блок-схему последовательности операций на фиг. 9, включает в себя следующее:
1) Вычисление азимута ψ0 и угла φ падения последовательным выборочным способом (то есть для каждого пространственно-временного положения, при котором имелось отражение амплитуды или энергии) на основании просуммированных трехмерных сейсмических данных (этап 91).
2) Вычисление масштабного коэффициента зсПг в скорректированных на нормальное приращение времени подборках трасс сейсмических данных с общей средней точкой в соответствии с уравнениями (4), (10) и (12) и применение масштабного коэффициента для масштабирования амплитуды в каждой выборке отражения (этап 92). Подборка трасс с общей средней точкой представляет собой набор сейсмических трасс, которые имеют общую среднюю точку (точку О на фиг. 7).
3) Инвертирование уравнения (13) по масштабированным подборкам трасс с общей средней точкой для получения параметров А, В, В' и ψ£ при предпочтительном использовании инверсии методом наименьших квадратов (этап 93).
Специалистам в данной области техники должно быть понятно, каким образом осуществляется ин версия методом наименьших квадратов.
Параметры А и В могут быть использованы для определения свойств породы и флюида в коллекторе, тогда как В' и ψ£ указывают на интенсивность и ориентацию трещин в коллекторе.
Этот способ включает в себя большее количество вычислений в процессе обработки данных, чем способ в случае 1. Он применим, когда наличие трещин предполагается или ожидается. Совместно с углом падения, полученным по трехмерным сейсмическим изображениям, масштабные коэффициенты могут быть найдены и применены к подборкам трасс с общей средней точкой. При инверсии азимутальной зависимости амплитуды отражения от удаления по масштабированной амплитуде получают атрибуты А и В свойств породы и атрибуты В' и ψ£ свойств трещин.
Случай 3: наличие трещин, а азимут трещин перпендикулярен к азимуту падения трещиноватого пласта (фиг. 4Ь)
1. Вывод формулы
Замечено, что в трещиноватых пластах азимут расположенных параллельно (вертикальных) трещин обычно перпендикулярен к направлению падения пласта. Это особый случай случая 2. В такой ситуации ψί=ψ0±90°. При подстановке уравнений (3) и (4) в уравнение (9) и замене ψ£ на ψ0+90° получаем
- 9 011314
Здесь А0, В0 и Βι обозначают соответственно отсекаемый отрезок, угловой коэффициент на графике изменения амплитуды отражения с удалением и параметр анизотропии при азимутальном изменении амплитуды отражения с удалением. Как и уравнение (6) в случае 1, уравнение (14) имеет такую же форму, как уравнение азимутального изменения амплитуды отражения с удалением (альтернативную форму уравнения (7)):
за исключением того, что азимут ψ,, трещин заменен на азимут ψ0 падения. N представляет собой разность между измеренным коэффициентом отражения (то есть амплитудой после коррекции на геометрическое расхождение) и коэффициентом отражения, вычисленным на основании уравнения (7). Точно так же, как и в случае 1, может быть найден простой способ нахождения атрибутов свойств породы и атрибутов трещин для этого специального случая.
2. Способ исключения влияния падения пласта при анализе азимутальной зависимости амплитуды отражения от удаления
С помощью выравнивающих коэффициентов подобных членов в уравнениях (14) и (15) получаются следующие соотношения:
А =-----Ао ,
СОЗ (р в = —[Во + (Во - | Ао) Сап2 φ], сов φ2 в' = —ί— [В, - (Во А0)1:ап2 <р\,(16) соз φ2 + ε .
Теперь оценка азимутального угла ψ'0 имеет небольшую погрешность ε относительно истинного азимутального угла падения или трещин. Это вызвано пренебрежением членов высшего порядка, представленных N. Точно так же, как в случае 1, этими соотношениями предлагается простой способ получения свойств породы и трещин без вредного влияния анизотропии, вызванной падением пласта. Процедура, сведенная в блок-схему последовательности операций на фиг. 10, изложена ниже для одного варианта осуществления способа настоящего изобретения и включает в себя:
1) Применение инверсии методом наименьших квадратов к трехмерным сейсмическим данным для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ'0. Атрибут ψ'0 представляет собой оцененный азимут падения трещиноватого наклонного пласта (этап 101) .
2) Вычисление углов падения и их косинусов (то есть со5<р) на целевом отражающем горизонте и умножение атрибута А на соответствующий косинус для получения скомпенсированного на наклон параметра А0, этап 102.
3) Суммирование атрибутов В и Β' для каждого пространственного положения на целевом отражающем горизонте и умножение суммы на косинус угла падения для этого положения (этап 103). Этот результат является суммой неизвестных В0 и В1.
4) Оценивание В0 с помощью обычной инверсии зависимости амплитуды отражения от удаления или угловых сумм (этап 104).
5) Для каждого положения вычитание В0 из значения (В+В')со§р из третьего этапа для получения анизотропного члена Βι углового коэффициента, который служит мерой азимутального изменения амплитуды, обусловленного наличием трещин (этап 105).
6) Выполнение итерации этого процесса (этапов 1-5) с только что оцененными атрибутами в качестве ограничивающих параметров при инверсии до тех пор, пока не будет достигнута необходимая точность (этап 106).
Пример
Способ случая 1 был применен к массиву трехмерных сейсмических наборов данных для демонстрации простоты и управляемой данными сущности способа. Сначала инверсия азимутальной зависимости амплитуды отражения от удаления была применена для получения атрибутов А, В, В' и ψ0 азимутальной зависимости амплитуды отражения от удаления. Результаты для первых трех атрибутов отображены на фиг. 11а-С, при этом положение вдоль линий приема отложено по горизонтальной оси, а время
- 10 011314 двойного пробега по вертикальной оси. Интенсивность атрибутов представлена с помощью полутоновой шкалы, показанной на фиг. 11а, на которой отображен атрибут А. Атрибуты В и В' отображены на фиг. 11Ь и 11с, соответственно. Атрибуты В и В' были суммированы для получения присущего пласту атрибута В0 углового коэффициента, а результаты показаны на фиг. 12. В этом случае согласно геологическим исследованиям угол падения пласта составляет меньше, чем 5°. (Наклоны осей синфазности могут быть непосредственно определены по изображениям атрибутов в каждом из этих разрезов. Наблюдаемый наклон кажется больше, чем 5°. Это является ошибочным восприятием, поскольку для удобства отображения горизонтальное расстояние сжато по оси абсцисс.) При таком небольшом угле падения нет необходимости умножать атрибуты А и (В+В') на οοκφ, чтобы получать присущие пласту атрибуты А0 и В0 (см. уравнение (8)), поскольку при таком небольшом угле падения οοκφ очень близок к 1,0. В случае большого угла падения (5°<φ<20°) для этого преобразования необходимо вычислять угол падения.
Должно быть понятно, что конкретные варианты осуществления изобретения могут быть реализованы как часть программы или прикладной программы, функционирующей в процессорной системе, такой, как компьютерная система. Компьютерная система может включать в себя серверы и рабочие станции, сконфигурированные для сохранения в запоминающем устройстве данных и прикладных программ, предназначенных для осуществления приведенных выше способов и процессов. Безусловно, компьютерная система может формировать отчет, который может быть отображен на мониторе, подключенном к компьютерной системе, сохранен в запоминающем устройстве компьютерной системы или распечатан на принтере, подключенном к компьютерной системе. Кроме того, как должно быть понятно специалистам в данной области техники, пользователь может взаимодействовать с компьютерной системой для выполнения вариантов осуществления, или способ может выполняться в компьютерной системе без вмешательства пользователя. Вне зависимости от этого по меньшей мере часть конкретных вариантов осуществления изобретения, которые рассмотрены выше, может сохраняться на считываемом компьютером носителе данных для выполнения компьютерной системой.
Изложенное выше описание касается конкретных вариантов осуществления настоящего изобретения, предназначенных для иллюстрации его. Однако специалисту в данной области техники должно быть понятно, что возможны многочисленные модификации и изменения к описанным вариантам осуществления. Все такие модификации и изменения подразумеваются находящимися в рамках объема настоящего изобретения, определенного прилагаемой формулой изобретения.

Claims (36)

  1. (1) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема сейсмических данных;
    (1) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объема данных с четырьмя атрибутами А, В, В' и ψα, соответствующими взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом расположенных параллельно трещин;
    (1) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема сейсмических данных;
    (1) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объема данных с четырьмя атрибутами А, В, В' и ψα, соответствующими взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом расположенных параллельно трещин;
    1. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;
    при этом указанный наклонный резервуар является изотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
  2. (2) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом указанные масштабные коэффициенты определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
    (2) вычисляют углы <р падения для верхней поверхности коллектора и умножают со»р на А для по
    (2) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом указанные масштабные коэффициенты определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
    2. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (2) вычисляют углы φ падения для верхней поверхности коллектора и умножают οοκφ на А для получения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зависимости амплитуды отражения от удаления для фактического наклонного коллектора;
  3. (3) применяют указанные масштабные коэффициенты к сейсмическим данным последовательным выборочным способом и (4) инвертируют указанное уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из масштабированных сейсмических данных для формирования объемов данных с четырьмя скомпенсированными за наклон атрибутами А, В, В' и ψ£ для фактического наклонного коллектора, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ£ является азимутом расположенных параллельно трещин.
    (3) вычисляют скомпенсированный за наклон атрибут В0 углового коэффициента зависимости амплитуды отражения от удаления для фактического наклонного коллектора из В0=(В+В')соэд и (4) вычисляют азимут ψ0 падения фактического наклонного коллектора из ψ0α.
    3. Способ по п.1, в котором наклонный резервуар является анизотропным, при этом указанная анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно азимуту падения коллектора, а при учете азимутальной зависимости расхождения энергии и угла отражения:
    (a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к трехмерным сейсмическим данным для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно трещин;
    (b) вычисляют углы φ падения для верхней поверхности коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора путем умножения А на οοκφ;
    (c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности коллектора и умножают сумму на со»р;
    (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора;
    (е) вычитают В0 из (В+В')со»р для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для фактического наклонного резервуара и (1) выполняют итерацию этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
    (3) применяют указанные масштабные коэффициенты к сейсмическим данным последовательным выборочным способом и (4) инвертируют указанное уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из масштабированных сейсмических данных для формирования объемов данных с четырьмя скомпенсированными за наклон атрибутами А, В, В' и ψ£ для фактического наклонного коллектора, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ£ является азимутом расположенных параллельно трещин.
    (3) вычисляют скомпенсированный за наклон атрибут В0 углового коэффициента зависимости амплитуды отражения от удаления для фактического наклонного коллектора из В0=(В+В')ео8<р и (4) вычисляют азимут ψ0 падения фактического наклонного коллектора из ψ0α.
  4. 4. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного изотропного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объемов данных с четырьмя атрибутами А, В, В' и ψα, соответствующими соотношению амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом гипотетических вертикальных трещин;
    (b) вычисляют углы φ падения для верхней поверхности наклонного изотропного коллектора и умножают со»р на А для получения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зави
  5. 5. Способ по п.4, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    К=А+В8т20+В'со82^^а)8ш20 для А, В, В' и ψα, где В - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных.
  6. 6. Способ по п.4, в котором углы φ падения вычисляют на основании объема трехмерных сейсмических данных после обычной обработки.
  7. 7. Способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема трехмерных сейсмических данных;
    (b) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из объема трехмерных сейсмических данных, при этом масштабный коэффициент определяют так, чтобы сделать амплитуды из объема трехмерных сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
    (c) применяют масштабный коэффициент к сейсмическим данным последовательным выборочным способом и (ά) инвертируют уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из объема трехмерных масштабированных сейсмических данных для формирования объемов данных с четырьмя скомпенсированными за наклон атрибутами А, В, В' и ψτ для фактического наклонного анизотропного коллектора, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψτ является азимутом расположенных параллельно трещин.
  8. 8. Способ по п.7, в котором уравнение азимутального изменения амплитуды отражения с удалением имеет вид
    К=А+В511Г0+В' ^82(ψ-ψα)8ΐη2θ, где В - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных.
  9. 9. Способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой, при этом анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно азимуту падения наклонного анизотропного коллектора, содержащий этапы, при выполнении которых:
    (a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно вертикальных трещин;
    (b) вычисляют углы φ падения для верхней поверхности наклонного анизотропного коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора путем умножения А на со»р;
    (c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности фактического наклонного анизотропного коллектора и умножают сумму на со»р;
    (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора;
  10. 10. Способ по п.9, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    К=А+В81и20 +В'со§2( ψ-ψ3)δίη2θ +Ν для А, В, В' и ψα, где Я - коэффициент отражения, θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных, при этом ψ,, заменяют на ψ0', а N представляет собой разность между результатом измерения Я и значением, вычисленным из предшествующего уравнения без члена Ν.
  11. 11. Способ по п.9, в котором В0 оценивают, используя обычную инверсию зависимости амплитуды отражения от удаления.
    - 11 011314 (a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;
    при этом наклонный резервуар является анизотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
  12. 12. Способ по п.9, в котором В0 оценивают, используя просуммированные по углам трассы.
    - 12 011314 симости амплитуды отражения от удаления для фактического наклонного изотропного коллектора;
    (с) вычисляют скомпенсированный за наклон атрибут В0 углового коэффициента зависимости амплитуды отражения от удаления для фактического наклонного изотропного коллектора из В0=(В=В') СО8ф и (ά) вычисляют азимут ψ0 падения фактического наклонного изотропного коллектора из ψ0α.
  13. 13. Способ оценивания магнитуды вызванной наклоном погрешности в вычисленных значениях атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, извлекаемых из сейсмических данных в результате трехмерного сейсмического исследования наклонного изотропного отражающего горизонта, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения с удалением, из
    Магнитуда = ^Во - £ап2 φ ξϊπ2 О, , где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, р - угол падения наклонного изотропного отражающего горизонта и θ1 - угол отражения, то есть характеристика волн в сейсмических данных, отраженных от поверхности наклонного изотропного отражающего горизонта, в предположении, что наклонный изотропный отражающий горизонт является плоским;
    (b) вычисляют величину базовой амплитуды, характерную для поверхности плоского отражающего горизонта, из
    Базовая амплитуда = А^В^т^ и
    (С) вычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
    - 13 011314 (е) вычитают В0 из (В+В')со§<р для каждого положения, чтобы оценить скомпенсированный за наклон анизотропный атрибут В1 углового коэффициента для фактического наклонного анизотропного коллектора; и (£) выполняют итерацию этапов с (а) по (е) с оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
  14. - 14 011314 лучения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зависимости амплитуды отражения от удаления для фактического наклонного коллектора;
    14. Способ по п.13, в котором дополнительно осуществляют компенсацию атрибутов зависимости амплитуды отражения от удаления за вызванную наклоном погрешность, если отношение оценки вызванной наклоном относительной погрешности к репрезентативному отношению шума к сигналу в сейсмических данных больше, чем предварительно определенное значение.
  15. - 15 011314 в сейсмических данных для наклонного пласта, при этом применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами;
    вычисляют угол падения для верхней поверхности наклонного пласта и умножают косинус угла падения на первый объем данных с атрибутом для получения скомпенсированного за наклон атрибута отсекаемого отрезка зависимости амплитуды отражения от удаления для наклонного пласта;
    вычисляют скомпенсированный за наклон атрибут углового коэффициента зависимости амплитуды отражения от удаления для наклонного пласта и вычисляют азимут ψ0 падения для наклонного пласта.
    15. Способ по п.14, в котором предварительно определенное значение равно единице.
  16. - 16 011314 тетическому плоскому слою с расположенными параллельно вертикальными трещинами, при этом множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления, объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления, объем данных с третьим атрибутом В', являющимся атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно трещин;
    (b) вычисляют угол φ падения для верхней поверхности наклонного пласта и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для наклонного пласта путем умножения А на со»р;
    (c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности наклонного пласта и умножают сумму на со»р;
    (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для наклонного пласта;
    (е) вычитают В0 из (В+В')со»р для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для наклонного пласта и (ί) выполняют итерацию указанных выше этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
    16. Реализуемый компьютером способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;
    при этом указанный наклонный резервуар является изотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
  17. - 17 011314
    17. Реализуемый компьютером способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:
    (a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;
    при этом наклонный резервуар является анизотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
  18. - 18 011314 висимости амплитуды отражения от удаления для наклонного пласта;
    (е) вычитают В0 из (В+В')со8ф для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для наклонного пласта и (ί) выполняют итерацию указанных выше этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
    18. Способ по п.16, в котором наклонный коллектор является анизотропным, при этом указанная анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно к азимуту падения коллектора, а при учете азимутальной зависимости расхождения энергии и угла отражения:
    (a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к трехмерным сейсмическим данным для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где Α и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно трещин;
    (b) вычисляют углы φ падения для верхней поверхности коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора путем умножения А на соэд;
    (c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности коллектора и умножают сумму на соэд;
    (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора;
    (е) вычитают В0 из (В+В')соэд для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для фактического наклонного резервуара и (ί) выполняют итерацию этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
  19. 19. Способ компенсации вызванных наклоном погрешностей, содержащий этапы, при выполнении которых получают сейсмические данные в результате сейсмического исследования наклонного пласта и осуществляют компенсацию за расхождение энергии и азимутальную зависимость угла отражения
  20. 20. Способ по п.19, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;
    объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;
    объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψα, являющимся азимутом гипотетических вертикальных трещин.
  21. 21. Способ по п.20, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    К=А+В81П20 +В'сО82( ψ-ψ3)δίη2θ для А, В, В' и ψα, где В - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, а азимут ψ0 падения равен четвертому атрибуту ψα из объема данных.
  22. 22. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют угол φ падения и азимутальный угол ψ0 падения для каждой из множества выборок из сейсмических данных;
    вычисляют по меньшей мере один масштабный коэффициент для каждой из множества выборок в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом по меньшей мере один масштабный коэффициент определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
    применяют по меньшей мере один масштабный коэффициент к сейсмическим данным для каждой из множества выборок и инвертируют уравнение азимутального изменения амплитуды отражения с удалением из масштабированных сейсмических данных для формирования множества объемов данных с атрибутами.
  23. 23. Способ по п.22, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;
    объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;
    объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψ£, являющимся азимутом расположенных параллельно вертикальных трещин.
  24. 24. Способ по п.23, в котором уравнение азимутального изменения амплитуды отражения с удалением имеет вид
    Κ=Α+Βδίη2θ +В'со82( ψ-ψί·)/δίη2θ, где К. - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в сейсмических данных.
  25. 25. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения:
    (а) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют гипо
  26. 26. Способ по п.25, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    Κ=Α+Β5ίη2θ+Β^Ο52(ψ-ψ,,)5ίη2θ+Ν для А, В, В' и ψα, где К - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, при этом ψα заменяют на ψ0', а N представляет собой разность между результатом измерения К и значением, вычисленным из предшествующего уравнения без члена Ν.
  27. 27. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения с удалением, из где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффици ента зависимости амплитуды отражения от удаления, соответственно, φ - угол падения наклонного отражающего горизонта и θ1 - характеристика угла отражения волн в сейсмических данных, отраженных от поверхности отражающего горизонта, в предположении, что отражающий горизонт является плоским;
    вычисляют величину базовой амплитуды, характерную для плоского отражающего горизонта, из Базовая амплитуда = Α00δίη2θ1 и
    вычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
  28. 28. Способ по п.27, в котором дополнительно осуществляют компенсацию атрибутов зависимости амплитуды отражения от удаления за вызванную наклоном погрешность, если отношение оценки вызванной наклоном относительной погрешности к репрезентативному отношению шума к сигналу в сейсмических данных больше, чем предварительно определенное значение.
  29. 29. Способ компенсации вызванных наклоном погрешностей, реализуемый компьютером, содержащий этапы, при выполнении которых получают сейсмические данные в результате сейсмического исследования наклонного пласта и осуществляют компенсацию за расхождение энергии и азимутальную зависимость угла отражения в сейсмических данных для наклонного пласта, при этом применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами;
    вычисляют угол падения для верхней поверхности наклонного пласта и умножают косинус угла падения на объем данных с первым атрибутом для получения скомпенсированного за наклон атрибута отсекаемого отрезка зависимости амплитуды отражения от удаления для наклонного пласта;
    вычисляют скомпенсированный за наклон атрибут углового коэффициента зависимости амплитуды отражения от удаления для наклонного пласта и вычисляют азимут ψ0 падения для наклонного пласта.
  30. 30. Способ по п.29, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;
    объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;
    объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно вертикальных трещин.
  31. 31. Способ по п.30, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    Κ=Α+Βδίη2θ +В'со§2( ψ-ψ3)δίη2θ для А, В, В' и ψα, где Я - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, а азимут ψ0 падения равен четвертому атрибуту ψα из объема данных.
  32. 32. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют угол р падения и азимутальный угол ψ0 падения для каждой из множества выборок из сейсмических данных;
    вычисляют по меньшей мере один масштабный коэффициент для каждой из множества выборок в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом по меньшей мере один масштабный коэффициент определяют так, чтобы сделать амплитуды сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;
    применяют по меньшей мере один масштабный коэффициент к сейсмическим данным для каждой из множества выборок и инвертируют уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из масштабированных сейсмических данных для формирования множества объемов данных с атрибутами.
  33. 33. Способ по п.32, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;
    объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;
    объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии и объем данных с четвертым атрибутом ψ£, являющимся азимутом расположенных параллельно вертикальных трещин.
  34. 34. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения:
    (a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, при этом множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления, объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления, объем данных с третьим атрибутом В', являющимся атрибутом азимутального изменения амплитуды отражения с удалением, относящимся к интенсивности анизотропии, и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно трещин;
    (b) вычисляют угол р падения для верхней поверхности наклонного пласта и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для наклонного пласта путем умножения А на со»р;
    (С) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности наклонного пласта и умножают сумму на со»р;
    (ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной за
  35. 35. Способ по п.34, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравнения
    Κ=Α+Β5ίΐΓθ+ΒΑθ52(ψ-ψ,,) δίινθ+Ν для А, В, В' и ψα, где К - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, при этом ψ, заменяют на ψ0, а N представляет собой разность между результатом измерения К и значением, вычисленным из предшествующего уравнения без члена Ν.
  36. 36. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения от удаления, из Магнитуда = ± А,/| Еап2 φ з!п2 Θ, , где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффици ента зависимости амплитуды отражения от удаления соответственно, φ - угол падения наклонного отражающего горизонта и θ1 - характеристика угла отражения волн в сейсмических данных, отраженных от поверхности отражающего горизонта, в предположении, что отражающий горизонт является плоским;
    вычисляют величину базовой амплитуды, характерную для плоского отражающего горизонта, из Базовая амплитуда =Α00δίη2θ1 и
    вычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
EA200701070A 2004-11-18 2005-10-14 Способ количественной оценки и снижения вызванного наклоном азимутального изменения амплитуды отражения от удаления EA011314B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US62906704P 2004-11-18 2004-11-18
PCT/US2005/036968 WO2006055137A2 (en) 2004-11-18 2005-10-14 Method for quantification and mitigation for dip-induced azimuthal avo

Publications (2)

Publication Number Publication Date
EA200701070A1 EA200701070A1 (ru) 2007-10-26
EA011314B1 true EA011314B1 (ru) 2009-02-27

Family

ID=34956493

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200701070A EA011314B1 (ru) 2004-11-18 2005-10-14 Способ количественной оценки и снижения вызванного наклоном азимутального изменения амплитуды отражения от удаления

Country Status (8)

Country Link
US (1) US7761237B2 (ru)
EP (1) EP1812811A2 (ru)
AU (1) AU2005307047B2 (ru)
BR (1) BRPI0518362A2 (ru)
CA (1) CA2586603C (ru)
EA (1) EA011314B1 (ru)
NO (1) NO20072189L (ru)
WO (1) WO2006055137A2 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2509890C1 (ru) * 2012-11-02 2014-03-20 Общество с ограниченной ответственностью "Инситу" Способ контроля и определения координат опасного состояния массива горных пород при подземных горных работах

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7830747B2 (en) * 2006-11-01 2010-11-09 Pgs Geophysical As Method for multi-azimuth prestack time migration for general heterogeneous, anisotropic media
US8335651B2 (en) * 2008-08-01 2012-12-18 Wave Imaging Technology, Inc. Estimation of propagation angles of seismic waves in geology with application to determination of propagation velocity and angle-domain imaging
US8463551B2 (en) 2009-11-17 2013-06-11 Schlumberger Technology Corporation Consistent dip estimation for seismic imaging
US8886503B2 (en) 2011-04-19 2014-11-11 Schlumberger Technology Corporation Discontinuity detection
US9542507B2 (en) 2012-02-23 2017-01-10 Schlumberger Technology Corporation Feature detection in seismic volumes
US9234971B2 (en) * 2013-11-18 2016-01-12 Nonlinear Seismic Imaging, Inc. Direct reservoir signature using the drag wave
WO2015145257A2 (en) * 2014-03-28 2015-10-01 Cgg Services Sa Method and apparatus for characterizing azimuthal anisotropy using statistical moments
US10401514B2 (en) * 2015-06-17 2019-09-03 Conocophillips Company Seismic azimuthal gradient estimation
CN105068117B (zh) * 2015-08-25 2016-06-08 中国矿业大学(北京) 用于裂缝性介质的avo反演方法、装置和设备
CN111158053B (zh) * 2019-12-20 2023-02-03 中石化石油工程技术服务有限公司 裂缝预测方法及装置
CN113820746A (zh) * 2020-06-18 2021-12-21 中国石油化工股份有限公司 一种地震油气检测方法及系统
CN112462421B (zh) * 2020-10-30 2024-08-27 中国石油天然气集团有限公司 储层信息预测方法、装置、电子设备及存储介质
CN112682034B (zh) * 2020-12-04 2022-12-09 中国地质大学(北京) 基于致密砂岩储层的裂缝识别、倾角表征的方法及装置
CN113075727B (zh) * 2021-02-20 2023-09-26 中国石油天然气股份有限公司 基于振幅方位导数的储层裂缝分布特征确定方法及系统
CN113109872B (zh) * 2021-04-07 2022-06-21 中国煤炭地质总局物测队 一种基于覆盖次数加权的叠前振幅补偿处理方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5903515A (en) * 1997-04-04 1999-05-11 Western Atlas International, Inc. Spatial de-aliasing of the 3-D Kirchhoff DMO operator
US20020007247A1 (en) * 2000-06-30 2002-01-17 Cheng Yao C. Method for imaging discontinuities in seismic data using dip-steering
US6681184B2 (en) * 2001-05-15 2004-01-20 Input/Output, Inc. System for estimating azimuthal variations in seismic data
US20040093163A1 (en) * 2002-11-12 2004-05-13 Moshe Reshef Seismic analysis using post-imaging seismic anisotropy corrections
US20040210394A1 (en) * 2001-08-31 2004-10-21 Henning Trappe Method for the determination of local similitude from seismic 3d measured data

Family Cites Families (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2800639A (en) * 1953-12-28 1957-07-23 Texas Co Method and apparatus for analysis of seismographic records
US3242326A (en) * 1954-10-26 1966-03-22 Sun Oil Co Method and apparatus for the analysis of seismic records
US2908889A (en) * 1955-12-16 1959-10-13 Phillips Petroleum Co Computer
US2926331A (en) * 1957-02-04 1960-02-23 Phillips Petroleum Co Seismic exploration
US3014650A (en) * 1957-06-19 1961-12-26 Jersey Prod Res Co Analysis of seismic records
US3209134A (en) * 1958-12-01 1965-09-28 Exxon Production Research Co Interpretation of geophysical data
FR1265258A (fr) * 1960-05-17 1961-06-30 Geophysique Cie Gle Perfectionnements aux méthodes et dispositifs d'analyse et d'examen des enregistrements sismiques
US3209317A (en) * 1961-04-03 1965-09-28 Exxon Production Research Co Seismic exploration
US3213411A (en) * 1962-01-12 1965-10-19 Phillips Petroleum Co Method of obtaining reverse path time tie between seismic signals
US4933913A (en) * 1986-10-30 1990-06-12 Amoco Corporation Method of seismic surveying for resolving the effects of formation anisotropy in shear wave reflection seismic data
US5508973A (en) * 1994-06-06 1996-04-16 Western Atlas International, Inc. Method for determining the principal axes of azimuthal anisotropy from seismic P-wave data
US5835452A (en) * 1995-10-06 1998-11-10 Amoco Corporation Reflected shear wave seismic processes
US5719822A (en) * 1996-10-04 1998-02-17 Vector Seismic Data Processing, Inc. Seismic data radon dip moveout method
US6061301A (en) * 1997-06-30 2000-05-09 Atlantic Richfield Company Filtering of overburden azimuthal anisotropy effects from 3D seismic survey signals
US5999486A (en) * 1998-07-23 1999-12-07 Colorado School Of Mines Method for fracture detection using multicomponent seismic data
US6055482A (en) * 1998-10-09 2000-04-25 Coherence Technology Company, Inc. Method of seismic signal processing
WO2000048022A1 (en) * 1999-02-12 2000-08-17 Schlumberger Limited Uncertainty constrained subsurface modeling
US6292754B1 (en) * 1999-11-11 2001-09-18 Bp Corporation North America Inc. Vector recomposition of seismic 3-D converted-wave data
US6826486B1 (en) * 2000-02-11 2004-11-30 Schlumberger Technology Corporation Methods and apparatus for predicting pore and fracture pressures of a subsurface formation
AU2001251019A1 (en) * 2000-03-27 2001-10-08 Peter J. Ortoleva Method for simulation of enhanced fracture detection in sedimentary basins
WO2002047011A1 (en) * 2000-12-08 2002-06-13 Ortoleva Peter J Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories
US6751558B2 (en) * 2001-03-13 2004-06-15 Conoco Inc. Method and process for prediction of subsurface fluid and rock pressures in the earth
US6473696B1 (en) * 2001-03-13 2002-10-29 Conoco Inc. Method and process for prediction of subsurface fluid and rock pressures in the earth
US6678617B2 (en) * 2001-06-13 2004-01-13 Nonlinear Seismic Imaging, Inc. Mapping subsurface open fractures using elastically nonlinear interaction of two surface-generated waves
GB2381313B (en) * 2001-07-18 2005-01-05 Westerngeco Ltd A method of processing geophysical data
US6760667B1 (en) * 2001-09-14 2004-07-06 Emerald Geoscience Research Corporation Method for generating P-S and S-S seismic data and attributes from P-P seismic data
US6817218B1 (en) * 2001-09-14 2004-11-16 Emerald Geoscience Research Corp. Non-linear statistically constrained AVO inversion
US6714873B2 (en) * 2001-12-17 2004-03-30 Schlumberger Technology Corporation System and method for estimating subsurface principal stresses from seismic reflection data
WO2003060560A1 (en) * 2002-01-15 2003-07-24 Westerngeco, L.L.C. Layer stripping converted reflected waveforms for dipping fractures
US6714480B2 (en) * 2002-03-06 2004-03-30 Schlumberger Technology Corporation Determination of anisotropic moduli of earth formations
US6928367B2 (en) * 2002-09-27 2005-08-09 Veritas Dgc Inc. Reservoir fracture characterization
US6896074B2 (en) * 2002-10-09 2005-05-24 Schlumberger Technology Corporation System and method for installation and use of devices in microboreholes
US6839633B1 (en) * 2003-06-13 2005-01-04 Schlumberger Technology Corporation Methods and apparatus for imaging a subsurface fracture
US7085195B2 (en) * 2004-03-30 2006-08-01 Rdsp I, L.P. Method for detecting earth formation fractures by seismic imaging of diffractors

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5903515A (en) * 1997-04-04 1999-05-11 Western Atlas International, Inc. Spatial de-aliasing of the 3-D Kirchhoff DMO operator
US20020007247A1 (en) * 2000-06-30 2002-01-17 Cheng Yao C. Method for imaging discontinuities in seismic data using dip-steering
US6681184B2 (en) * 2001-05-15 2004-01-20 Input/Output, Inc. System for estimating azimuthal variations in seismic data
US20040210394A1 (en) * 2001-08-31 2004-10-21 Henning Trappe Method for the determination of local similitude from seismic 3d measured data
US20040093163A1 (en) * 2002-11-12 2004-05-13 Moshe Reshef Seismic analysis using post-imaging seismic anisotropy corrections

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2509890C1 (ru) * 2012-11-02 2014-03-20 Общество с ограниченной ответственностью "Инситу" Способ контроля и определения координат опасного состояния массива горных пород при подземных горных работах

Also Published As

Publication number Publication date
US7761237B2 (en) 2010-07-20
BRPI0518362A2 (pt) 2008-11-18
CA2586603C (en) 2014-04-01
WO2006055137A3 (en) 2006-08-10
EP1812811A2 (en) 2007-08-01
AU2005307047B2 (en) 2010-11-18
NO20072189L (no) 2007-08-13
CA2586603A1 (en) 2006-05-26
WO2006055137A2 (en) 2006-05-26
US20070260404A1 (en) 2007-11-08
AU2005307047A1 (en) 2006-05-26
EA200701070A1 (ru) 2007-10-26

Similar Documents

Publication Publication Date Title
EA011314B1 (ru) Способ количественной оценки и снижения вызванного наклоном азимутального изменения амплитуды отражения от удаления
Koren et al. Full-azimuth subsurface angle domain wavefield decomposition and imaging Part I: Directional and reflection image gathers
Dell et al. Common-reflection-surface-based workflow for diffraction imaging
Marfurt Robust estimates of 3D reflector dip and azimuth
RU2457513C2 (ru) Способы и системы для обработки микросейсмических данных
Dost et al. Probabilistic moment tensor inversion for hydrocarbon‐induced seismicity in the Groningen gas field, the Netherlands, part 2: Application
Shen et al. Q-model building using one-way wave-equation migration Q analysis—Part 1: Theory and synthetic test
US20160061975A1 (en) Vertical Seismic Profiling Migration Method
Bakulin et al. Building tilted transversely isotropic depth models using localized anisotropic tomography with well information
US20100118652A1 (en) Determination of depth moveout and of residual radii of curvature in the common angle domain
EA027554B1 (ru) Способ разведки углеводородов
US20140324354A1 (en) Transmission coefficient method for avo seismic analysis
Yordkayhun et al. 3D seismic traveltime tomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, Germany
Rabben et al. AVA inversion of the top Utsira Sand reflection at the Sleipner field
CN115877449A (zh) 用于在勘测区域内获得地下堆叠图像的计算机实现方法
Li et al. Analysis of RTM extended images for VTI media
Hendrickson Stacked
Zhu et al. Recent applications of turning-ray tomography
Liu et al. Tomographic velocity model building of the near surface with velocity-inversion interfaces: A test using the Yilmaz model
Khoshnavaz et al. Oriented surface passive seismic location using local slopes
Jusri et al. Advanced three‐dimensional seismic imaging of deep supercritical geothermal rocks in Southern Tuscany
Baina et al. 3D preserved-amplitude prestack depth migration and amplitude versus angle relevance
Zhou et al. Use of RTM full 3D subsurface angle gathers for subsalt velocity update and image optimization: Case study at Shenzi field
Palmer Non‐uniqueness with refraction inversion–the Mt Bulga shear zone
Xu et al. A case study of azimuthal AVO analysis with anisotropic spreading correction

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