EA011314B1 - Способ количественной оценки и снижения вызванного наклоном азимутального изменения амплитуды отражения от удаления - Google Patents
Способ количественной оценки и снижения вызванного наклоном азимутального изменения амплитуды отражения от удаления Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 96
- 230000000116 mitigating effect Effects 0.000 title abstract 2
- 238000011002 quantification Methods 0.000 title 1
- 230000015572 biosynthetic process Effects 0.000 claims description 64
- 230000008859 change Effects 0.000 claims description 62
- 239000011435 rock Substances 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 8
- 238000012217 deletion Methods 0.000 claims description 2
- 230000037430 deletion Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 3
- 238000011156 evaluation Methods 0.000 claims 1
- 230000005855 radiation Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 9
- 238000012937 correction Methods 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 9
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000007598 dipping method Methods 0.000 abstract 1
- 238000005755 formation reaction Methods 0.000 description 43
- 238000004880 explosion Methods 0.000 description 12
- 230000001419 dependent effect Effects 0.000 description 9
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 239000012530 fluid Substances 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 102100029095 Exportin-1 Human genes 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 235000013601 eggs Nutrition 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 108700002148 exportin 1 Proteins 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000005070 ripening Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
- G01V2210/632—Amplitude 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 -угол отражения, то есть характеристика волн в сейсмических данных, отраженных от поверхности наклонного изотропного отражающего горизонта, в предположении, что наклонный изотропный отражающий горизонт является плоским; (Ь) вычисляют величину базовой амплитуды, характерную для поверхности плоского изотропного отражающего горизонта, из
Базовая амплитуда = Α0+Β0δίη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 в1Г (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) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема сейсмических данных;(1) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объема данных с четырьмя атрибутами А, В, В' и ψα, соответствующими взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом расположенных параллельно трещин;(1) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема сейсмических данных;(1) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объема данных с четырьмя атрибутами А, В, В' и ψα, соответствующими взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом расположенных параллельно трещин;1. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;при этом указанный наклонный резервуар является изотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
- (2) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом указанные масштабные коэффициенты определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;(2) вычисляют углы <р падения для верхней поверхности коллектора и умножают со»р на А для по(2) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом указанные масштабные коэффициенты определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;2. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(2) вычисляют углы φ падения для верхней поверхности коллектора и умножают οοκφ на А для получения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зависимости амплитуды отражения от удаления для фактического наклонного коллектора;
- (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. Способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного изотропного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объемов данных с четырьмя атрибутами А, В, В' и ψα, соответствующими соотношению амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψα является азимутом гипотетических вертикальных трещин;(b) вычисляют углы φ падения для верхней поверхности наклонного изотропного коллектора и умножают со»р на А для получения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зави
- 5. Способ по п.4, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияК=А+В8т20+В'со82^^а)8ш20 для А, В, В' и ψα, где В - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных.
- 6. Способ по п.4, в котором углы φ падения вычисляют на основании объема трехмерных сейсмических данных после обычной обработки.
- 7. Способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) вычисляют угол φ падения и азимутальный угол ψ0 падения последовательным выборочным способом на основании объема трехмерных сейсмических данных;(b) вычисляют масштабный коэффициент для каждой выборки в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из объема трехмерных сейсмических данных, при этом масштабный коэффициент определяют так, чтобы сделать амплитуды из объема трехмерных сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;(c) применяют масштабный коэффициент к сейсмическим данным последовательным выборочным способом и (ά) инвертируют уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из объема трехмерных масштабированных сейсмических данных для формирования объемов данных с четырьмя скомпенсированными за наклон атрибутами А, В, В' и ψτ для фактического наклонного анизотропного коллектора, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψτ является азимутом расположенных параллельно трещин.
- 8. Способ по п.7, в котором уравнение азимутального изменения амплитуды отражения с удалением имеет видК=А+В511Г0+В' ^82(ψ-ψα)8ΐη2θ, где В - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных.
- 9. Способ компенсации вызванных наклоном погрешностей в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного анизотропного коллектора, лежащего под изотропной покрывающей породой, при этом анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно азимуту падения наклонного анизотропного коллектора, содержащий этапы, при выполнении которых:(a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к объему трехмерных сейсмических данных для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где А и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно вертикальных трещин;(b) вычисляют углы φ падения для верхней поверхности наклонного анизотропного коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора путем умножения А на со»р;(c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности фактического наклонного анизотропного коллектора и умножают сумму на со»р;(ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного анизотропного коллектора;
- 10. Способ по п.9, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияК=А+В81и20 +В'со§2( ψ-ψ3)δίη2θ +Ν для А, В, В' и ψα, где Я - коэффициент отражения, θ - угол отражения для каждой отраженной волны в объеме трехмерных сейсмических данных, при этом ψ,, заменяют на ψ0', а N представляет собой разность между результатом измерения Я и значением, вычисленным из предшествующего уравнения без члена Ν.
- 11. Способ по п.9, в котором В0 оценивают, используя обычную инверсию зависимости амплитуды отражения от удаления.- 11 011314 (a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;при этом наклонный резервуар является анизотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
- 12. Способ по п.9, в котором В0 оценивают, используя просуммированные по углам трассы.- 12 011314 симости амплитуды отражения от удаления для фактического наклонного изотропного коллектора;(с) вычисляют скомпенсированный за наклон атрибут В0 углового коэффициента зависимости амплитуды отражения от удаления для фактического наклонного изотропного коллектора из В0=(В=В') СО8ф и (ά) вычисляют азимут ψ0 падения фактического наклонного изотропного коллектора из ψ0=ψα.
- 13. Способ оценивания магнитуды вызванной наклоном погрешности в вычисленных значениях атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления, извлекаемых из сейсмических данных в результате трехмерного сейсмического исследования наклонного изотропного отражающего горизонта, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения с удалением, изМагнитуда = ^Во - £ап2 φ ξϊπ2 О, , где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, р - угол падения наклонного изотропного отражающего горизонта и θ1 - угол отражения, то есть характеристика волн в сейсмических данных, отраженных от поверхности наклонного изотропного отражающего горизонта, в предположении, что наклонный изотропный отражающий горизонт является плоским;(b) вычисляют величину базовой амплитуды, характерную для поверхности плоского отражающего горизонта, изБазовая амплитуда = А^В^т^ и(С) вычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.- 13 011314 (е) вычитают В0 из (В+В')со§<р для каждого положения, чтобы оценить скомпенсированный за наклон анизотропный атрибут В1 углового коэффициента для фактического наклонного анизотропного коллектора; и (£) выполняют итерацию этапов с (а) по (е) с оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
- - 14 011314 лучения скомпенсированного за наклон атрибута А0 отсекаемого отрезка зависимости амплитуды отражения от удаления для фактического наклонного коллектора;14. Способ по п.13, в котором дополнительно осуществляют компенсацию атрибутов зависимости амплитуды отражения от удаления за вызванную наклоном погрешность, если отношение оценки вызванной наклоном относительной погрешности к репрезентативному отношению шума к сигналу в сейсмических данных больше, чем предварительно определенное значение.
- - 15 011314 в сейсмических данных для наклонного пласта, при этом применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами;вычисляют угол падения для верхней поверхности наклонного пласта и умножают косинус угла падения на первый объем данных с атрибутом для получения скомпенсированного за наклон атрибута отсекаемого отрезка зависимости амплитуды отражения от удаления для наклонного пласта;вычисляют скомпенсированный за наклон атрибут углового коэффициента зависимости амплитуды отражения от удаления для наклонного пласта и вычисляют азимут ψ0 падения для наклонного пласта.15. Способ по п.14, в котором предварительно определенное значение равно единице.
- - 16 011314 тетическому плоскому слою с расположенными параллельно вертикальными трещинами, при этом множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления, объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления, объем данных с третьим атрибутом В', являющимся атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно трещин;(b) вычисляют угол φ падения для верхней поверхности наклонного пласта и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для наклонного пласта путем умножения А на со»р;(c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности наклонного пласта и умножают сумму на со»р;(ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для наклонного пласта;(е) вычитают В0 из (В+В')со»р для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для наклонного пласта и (ί) выполняют итерацию указанных выше этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.16. Реализуемый компьютером способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;при этом указанный наклонный резервуар является изотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
- - 17 01131417. Реализуемый компьютером способ компенсации вызванных наклоном погрешностей в атрибутах зависимости амплитуды отражения от удаления или в атрибутах азимутальной зависимости амплитуды отражения от удаления, извлекаемых из объема трехмерных сейсмических данных, получаемых в результате сейсмического исследования наклонного коллектора, лежащего под изотропной покрывающей породой, содержащий этапы, при выполнении которых:(a) учитывают азимутальную зависимость расхождения энергии в сейсмической волне, отраженной от наклонного коллектора; и (b) учитывают азимутальную зависимость угла отражения волны, отраженной от наклонного коллектора, в объеме трехмерных сейсмических данных;при этом наклонный резервуар является анизотропным, а при учете азимутальной зависимости расхождения энергии и угла отражения:
- - 18 011314 висимости амплитуды отражения от удаления для наклонного пласта;(е) вычитают В0 из (В+В')со8ф для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для наклонного пласта и (ί) выполняют итерацию указанных выше этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии азимутальной зависимости амплитуды отражения от удаления до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.18. Способ по п.16, в котором наклонный коллектор является анизотропным, при этом указанная анизотропия в значительной степени обусловлена трещинами, ориентированными, по существу, перпендикулярно к азимуту падения коллектора, а при учете азимутальной зависимости расхождения энергии и угла отражения:(a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к трехмерным сейсмическим данным для формирования объемов данных с четырьмя атрибутами А, В, В' и ψ0', соответствующими гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, где Α и В в азимутальной зависимости амплитуды отражения от удаления являются эквивалентами атрибутов отсекаемого отрезка и углового коэффициента зависимости амплитуды отражения от удаления соответственно, В' является атрибутом азимутальной зависимости амплитуды отражения от удаления, относящимся к интенсивности анизотропии, и ψ0' является оцененным азимутом расположенных параллельно трещин;(b) вычисляют углы φ падения для верхней поверхности коллектора и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора путем умножения А на соэд;(c) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности коллектора и умножают сумму на соэд;(ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной зависимости амплитуды отражения от удаления для фактического наклонного коллектора;(е) вычитают В0 из (В+В')соэд для каждого положения с целью оценивания скомпенсированного за наклон анизотропного атрибута В1 углового коэффициента для фактического наклонного резервуара и (ί) выполняют итерацию этапов с (а) по (е) с только что оцененными атрибутами в качестве ограничивающих параметров для инверсии до тех пор, пока не будет достигнут предварительно определенный критерий сходимости.
- 19. Способ компенсации вызванных наклоном погрешностей, содержащий этапы, при выполнении которых получают сейсмические данные в результате сейсмического исследования наклонного пласта и осуществляют компенсацию за расхождение энергии и азимутальную зависимость угла отражения
- 20. Способ по п.19, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψα, являющимся азимутом гипотетических вертикальных трещин.
- 21. Способ по п.20, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияК=А+В81П20 +В'сО82( ψ-ψ3)δίη2θ для А, В, В' и ψα, где В - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, а азимут ψ0 падения равен четвертому атрибуту ψα из объема данных.
- 22. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют угол φ падения и азимутальный угол ψ0 падения для каждой из множества выборок из сейсмических данных;вычисляют по меньшей мере один масштабный коэффициент для каждой из множества выборок в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом по меньшей мере один масштабный коэффициент определяют так, чтобы сделать амплитуды из сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;применяют по меньшей мере один масштабный коэффициент к сейсмическим данным для каждой из множества выборок и инвертируют уравнение азимутального изменения амплитуды отражения с удалением из масштабированных сейсмических данных для формирования множества объемов данных с атрибутами.
- 23. Способ по п.22, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψ£, являющимся азимутом расположенных параллельно вертикальных трещин.
- 24. Способ по п.23, в котором уравнение азимутального изменения амплитуды отражения с удалением имеет видΚ=Α+Βδίη2θ +В'со82( ψ-ψί·)/δίη2θ, где К. - коэффициент отражения, а θ - угол отражения для каждой отраженной волны в сейсмических данных.
- 25. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения:(а) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют гипо
- 26. Способ по п.25, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияΚ=Α+Β5ίη2θ+Β^Ο52(ψ-ψ,,)5ίη2θ+Ν для А, В, В' и ψα, где К - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, при этом ψα заменяют на ψ0', а N представляет собой разность между результатом измерения К и значением, вычисленным из предшествующего уравнения без члена Ν.
- 27. Способ по п.19, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения с удалением, из где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффици ента зависимости амплитуды отражения от удаления, соответственно, φ - угол падения наклонного отражающего горизонта и θ1 - характеристика угла отражения волн в сейсмических данных, отраженных от поверхности отражающего горизонта, в предположении, что отражающий горизонт является плоским;вычисляют величину базовой амплитуды, характерную для плоского отражающего горизонта, из Базовая амплитуда = Α0+Β0δίη2θ1 ивычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
- 28. Способ по п.27, в котором дополнительно осуществляют компенсацию атрибутов зависимости амплитуды отражения от удаления за вызванную наклоном погрешность, если отношение оценки вызванной наклоном относительной погрешности к репрезентативному отношению шума к сигналу в сейсмических данных больше, чем предварительно определенное значение.
- 29. Способ компенсации вызванных наклоном погрешностей, реализуемый компьютером, содержащий этапы, при выполнении которых получают сейсмические данные в результате сейсмического исследования наклонного пласта и осуществляют компенсацию за расхождение энергии и азимутальную зависимость угла отражения в сейсмических данных для наклонного пласта, при этом применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют взаимосвязи амплитуда-удаление-азимут гипотетического плоского слоя с расположенными параллельно вертикальными трещинами;вычисляют угол падения для верхней поверхности наклонного пласта и умножают косинус угла падения на объем данных с первым атрибутом для получения скомпенсированного за наклон атрибута отсекаемого отрезка зависимости амплитуды отражения от удаления для наклонного пласта;вычисляют скомпенсированный за наклон атрибут углового коэффициента зависимости амплитуды отражения от удаления для наклонного пласта и вычисляют азимут ψ0 падения для наклонного пласта.
- 30. Способ по п.29, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии; и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно вертикальных трещин.
- 31. Способ по п.30, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияΚ=Α+Βδίη2θ +В'со§2( ψ-ψ3)δίη2θ для А, В, В' и ψα, где Я - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, а азимут ψ0 падения равен четвертому атрибуту ψα из объема данных.
- 32. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют угол р падения и азимутальный угол ψ0 падения для каждой из множества выборок из сейсмических данных;вычисляют по меньшей мере один масштабный коэффициент для каждой из множества выборок в скорректированных за нормальное приращение времени подборках трасс с общей средней точкой из сейсмических данных, при этом по меньшей мере один масштабный коэффициент определяют так, чтобы сделать амплитуды сейсмических данных удовлетворяющими уравнению азимутального изменения амплитуды отражения с удалением для гипотетического плоского анизотропного слоя, содержащего расположенные параллельно трещины;применяют по меньшей мере один масштабный коэффициент к сейсмическим данным для каждой из множества выборок и инвертируют уравнение азимутального изменения амплитуды отражения с удалением по подборкам трасс с общей средней точкой из масштабированных сейсмических данных для формирования множества объемов данных с атрибутами.
- 33. Способ по п.32, в котором множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления;объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления;объем данных с третьим атрибутом В', являющимся в азимутальной зависимости амплитуды отражения от удаления атрибутом, относящимся к интенсивности анизотропии и объем данных с четвертым атрибутом ψ£, являющимся азимутом расположенных параллельно вертикальных трещин.
- 34. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения:(a) применяют инверсию азимутальной зависимости амплитуды отражения от удаления к сейсмическим данным для формирования множества объемов данных с атрибутами, которые соответствуют гипотетическому плоскому слою с расположенными параллельно вертикальными трещинами, при этом множество объемов данных с атрибутами содержит объем данных с первым атрибутом А, который в азимутальной зависимости амплитуды отражения от удаления является эквивалентом атрибутов отсекаемого отрезка зависимости амплитуды отражения от удаления, объем данных со вторым атрибутом В, являющимся в азимутальной зависимости амплитуды отражения от удаления эквивалентом атрибутов углового коэффициента зависимости амплитуды отражения от удаления, объем данных с третьим атрибутом В', являющимся атрибутом азимутального изменения амплитуды отражения с удалением, относящимся к интенсивности анизотропии, и объем данных с четвертым атрибутом ψα, являющимся азимутом расположенных параллельно трещин;(b) вычисляют угол р падения для верхней поверхности наклонного пласта и получают скомпенсированный за наклон атрибут А0 отсекаемого отрезка азимутальной зависимости амплитуды отражения от удаления для наклонного пласта путем умножения А на со»р;(С) суммируют атрибуты В и В' для каждого пространственного положения на верхней поверхности наклонного пласта и умножают сумму на со»р;(ά) оценивают скомпенсированный за наклон атрибут В0 углового коэффициента азимутальной за
- 35. Способ по п.34, в котором инверсию азимутальной зависимости амплитуды отражения от удаления выражают с помощью уравненияΚ=Α+Β5ίΐΓθ+ΒΑθ52(ψ-ψ,,) δίινθ+Ν для А, В, В' и ψα, где К - коэффициент отражения, θ - угол отражения для каждой отраженной волны в сейсмических данных, при этом ψ, заменяют на ψ0, а N представляет собой разность между результатом измерения К и значением, вычисленным из предшествующего уравнения без члена Ν.
- 36. Способ по п.29, в котором при осуществлении компенсации за расхождение энергии и азимутальную зависимость угла отражения вычисляют величину магнитуды, представляющую максимальную степень изменения амплитуды при вызванном наклоном азимутальном изменении амплитуды отражения от удаления, из Магнитуда = ± А,/| Еап2 φ з!п2 Θ, , где А0 и В0 являются представляющими значения атрибутов отсекаемого отрезка и углового коэффици ента зависимости амплитуды отражения от удаления соответственно, φ - угол падения наклонного отражающего горизонта и θ1 - характеристика угла отражения волн в сейсмических данных, отраженных от поверхности отражающего горизонта, в предположении, что отражающий горизонт является плоским;вычисляют величину базовой амплитуды, характерную для плоского отражающего горизонта, из Базовая амплитуда =Α0+Β0δίη2θ1 ивычисляют оценку вызванной наклоном относительной погрешности путем деления величины магнитуды на величину базовой амплитуды.
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2509890C1 (ru) * | 2012-11-02 | 2014-03-20 | Общество с ограниченной ответственностью "Инситу" | Способ контроля и определения координат опасного состояния массива горных пород при подземных горных работах |
Families Citing this family (15)
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)
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)
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 |
-
2005
- 2005-10-14 WO PCT/US2005/036968 patent/WO2006055137A2/en active Application Filing
- 2005-10-14 EP EP05812478A patent/EP1812811A2/en not_active Withdrawn
- 2005-10-14 CA CA2586603A patent/CA2586603C/en not_active Expired - Fee Related
- 2005-10-14 EA EA200701070A patent/EA011314B1/ru not_active IP Right Cessation
- 2005-10-14 US US11/666,579 patent/US7761237B2/en active Active
- 2005-10-14 BR BRPI0518362-6A patent/BRPI0518362A2/pt not_active IP Right Cessation
- 2005-10-14 AU AU2005307047A patent/AU2005307047B2/en not_active Ceased
-
2007
- 2007-04-27 NO NO20072189A patent/NO20072189L/no not_active Application Discontinuation
Patent Citations (5)
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)
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 |