RU2674419C1 - Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса - Google Patents
Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса Download PDFInfo
- Publication number
- RU2674419C1 RU2674419C1 RU2018103623A RU2018103623A RU2674419C1 RU 2674419 C1 RU2674419 C1 RU 2674419C1 RU 2018103623 A RU2018103623 A RU 2018103623A RU 2018103623 A RU2018103623 A RU 2018103623A RU 2674419 C1 RU2674419 C1 RU 2674419C1
- Authority
- RU
- Russia
- Prior art keywords
- seismograms
- seismic
- common
- total
- images
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 229910052704 radon Inorganic materials 0.000 claims abstract description 13
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 13
- 230000005284 excitation Effects 0.000 claims abstract description 5
- 238000013508 migration Methods 0.000 claims abstract description 3
- 230000005012 migration Effects 0.000 claims abstract description 3
- 238000003384 imaging method Methods 0.000 claims 1
- 229930195733 hydrocarbon Natural products 0.000 abstract description 5
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 5
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 4
- 230000004807 localization Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 230000009466 transformation Effects 0.000 description 5
- 210000001233 cdp Anatomy 0.000 description 2
- 238000004637 computerized dynamic posturography Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000000969 carrier Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V1/364—Seismic filtering
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)
Abstract
Изобретение относится к области сейсморазведки, а именно к методам построения разрезов геологической среды по сейсмическим данным (сейсмических разрезов), позволяющий, используя различие свойств отраженных и рассеянных событий на сейсмограммах общего выноса, более устойчиво (надежно) и с меньшими затратами локализовать области рассеяния (дифракции). Способ состоит в том, что производят пересортировку сейсмограмм общей точки возбуждения в сейсмограммы общего выноса, затем осуществляют локализацию зеркальных отражений с помощью преобразования Радона на каждой из полученных сейсмограмм общего выноса и маскируют локализованные зеркальные отражения на полученных Радон образах, после чего выполняют обратное преобразование Радона над маскированными образами и получают искомые сейсмограммы общего выноса, которые подвергают обратной пересортировке в сейсмограммы общей точки возбуждения и осуществляют построение изображения с помощью одного из известных алгоритмов миграции. Технический результат – повышение информативности получаемых данных за счет обеспечения возможности обнаружения неструктурных элементов, которые зачастую могут быть ловушками углеводородов. 4 ил.
Description
Область техники
Изобретение относится к области сейсморазведки, а именно к методам построения разрезов геологической среды по сейсмическим данным (сейсмических разрезов).
Уровень техники
Существующие методы построения сейсмических разрезов традиционно ориентированы на изображение структурных объектов т.е. гладких, протяженных границ, часто являющихся границами коллекторов углеводородов (1). Системы наблюдений довольно разнообразны и зависят от целей, условий наблюдения и других факторов, но всегда данные наблюдения являются переопределенными. В случае профильных наблюдений наиболее распространенной является система наблюдений многократных перекрытий (система наблюдений ОГТ). Она осуществляется путем регистрации сейсмических колебаний при заданных координатах источника и расстановки приемников на профиле и последующих сдвигов источника и всей расстановки приемников на определенный шаг и снова регистрации. Таким образом, получают функцию трех переменных, координат источника xs приемника xr и времени t:U(xs,xr,t). При фиксированной координате источника получаем функцию Us(xr,t), называемую сейсмограммой общей точки возбуждения (ОТВ). Часто вместо переменных, xs, xr удобнее использовать переменные xm=(xr+xs)/2 и xh=(xr-xs)/2, называемые координатой средней точки и координатой выноса, соответственно. При фиксированной переменной xh получаем функцию Um(xm,t), называемую сейсмограммой общей средней точки (ОСТ), при фиксации же переменной xm получаем функцию Uh(xh,t), называемую сейсмограммой общего выноса (ОВ).
Существующие методы обработки сейсмических данных, ориентированные на поиск гладких синклинальных структур в геологической среде, в зависимости от геологической ситуации (наличие априорной информации о скоростях, поверхностные условия, рельеф и т.д.) работают с различными из перечисленных комбинаций данных.
Однако, установлено, что ловушками углеводородов могут служить и различные неструктурные элементы разреза, такие как разломы, рифы, диапиры, области раскарстования и др. Носителями информации для таких объектов являются рассеянные (дифрагированные) волны, которые при традиционной обработке частично или полностью маскируются зеркальными, гораздо более интенсивными отражениями. В последние полтора десятилетия стали обращать внимание на методы, использующие этот тип информации о геологических объектах. Вследствие малой интенсивности рассеянных волн эти методы являются более тонкими и более трудоемкими, и могут использоваться только после проведения комплекса стандартной обработки данных.
Заявляемое изобретение рассчитано на устранение указанного недостатка, что может помочь геологу-интерпретатору при определении коллекторов углеводородов и проектировании бурения.
Наиболее близким к заявленному техническому решению является метод, опубликованный в работе (2), где для фокусировки отраженных сигналов на сейсмограммах ОТВ используется обобщенное гиперболическое преобразование Радона, в связи с близостью их к гиперболам или к гиперболоидам в трехмерном случае. Однако годографы рассеянных волн в этом случае также представляются гиперболами и разделение отраженных и рассеянных сигналов не всегда производится надежно.
Раскрытие изобретения
Задачей, на решение которой направлено заявляемое изобретение, является создание способа, позволяющего построение такого изображения геологической среды, на котором подавляются зеркальные отражатели и, наоборот, подчеркиваются рассеивающие элементы геологического разреза.
Технический результат состоит в возможности обнаружения неструктурных элементов, которые зачастую могут быть ловушками углеводородов.
В заявляемом методе используются сейсмограммы общего выноса (ОВ), на которых сигналы, отраженные от границ с малыми кривизнами, представляются почти плоскими событиями, а сигналы, порождаемые мелкомасштабными, рассеивающими объектами, остаются близкими к кривым второго порядка. Такое различие в кривизне сигналов является наиболее благоприятной ситуацией для их разделения.
Преобразование Радона
осуществляющее интегрирование по всем прямым в плоскости определения двумерной функции, является идеальным инструментом для разделения упомянутых событий, где x, t - переменные, определяющие исходную сейсмограмму, k - тангенс угла наклона прямой, τ - сдвиг прямой.
После применения указанного преобразования к полученной сейсмограмме общего выноса, события с осями синфазности близкими к прямолинейным окажутся локализованными в ограниченной области Радон образа, в то время как другие события окажутся распределенными по всей области. Это позволит путем маскирования областей наибольшей концентрации исключить основную энергию отраженных событий в пространстве Радон образа оставив почти неизмененными другие события.
Существует обратное преобразование Радона
u(x,t)=R-1 (k, τ),
которое приближенно можно вычислить, например, используя преобразование Фурье или линейно алгебраическим способом.
Применение обратного преобразования Радона к полученному отфильтрованному образу должно существенно усилить рассеянную компоненту волнового поля по сравнению с зеркальной компонентой.
Осуществление способа
Для демонстрации работоспособности метода ниже приводятся результаты теста, проведенного по упрощенной схеме расчета на примере данных, полученных по стандартной системе ОГТ для 2-d акустической модели среды. Упрощение заключается в том, что вместо обратного преобразования Радона было использовано более простое с точки зрения программирования и вычислительных затрат сопряженное преобразование, также называемое оператором обратного проектирования:
где тильда обозначает некоторое приближение к обратному оператору. Это преобразование хорошо известно в геофизике, является кинематически эквивалентным обратному преобразованию и часто дает вполне приемлемое приближение.
Тестовая модель, приведенная на рис. 1 содержит три слоя в последнем (тонком) из которых расположена область со случайно распределенными кавернами, которые хорошо видны на рис. 2, изображающем увеличенный фрагмент модели.
Данные для системы многократных перекрытий были рассчитаны конечно-разностным методом.
Общепринятый метод миграции Кирхгофа построения изображения среды, использующий сейсмограммы ОТВ данных многократного перекрытия рассчитанным для приведенной модели (рис. 3) дает едва уловимое изображение области рассеяния, тогда как изображение рассеянных элементов (рис. 4), построенное с помощью заявляемого метода, на тех же данных, более четко локализует эту область.
Таким образом, предлагаемый способ в связи с более сильным различием свойств отраженных и рассеянных событий, именно, на сейсмограммах ОВ должен позволить более устойчиво (надежно) и с меньшими затратами локализовать области рассеяния (дифракции).
Использованные источники информации
1. Р. Шерифф, Л. Гелдарт. Сейсморазведка, Москва, "Мир", 1987, т. 1, 2.
2. Khaidukov V., Landa Е., Moser T.J. Diffraction imaging by focusing-defocusing: An outlook on seismic super resolution //, Geophysics, 2004, v. 69, p. 1478-1490.
Claims (1)
- Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса, состоящий в том, что производят пересортировку сейсмограмм общей точки возбуждения в сейсмограммы общего выноса, затем осуществляют локализацию зеркальных отражений с помощью преобразования Радона на каждой из полученных сейсмограмм общего выноса и маскируют локализованные зеркальные отражения на полученных Радон образах, после чего выполняют обратное преобразование Радона над маскированными образами и получают искомые сейсмограммы общего выноса, которые подвергают обратной пересортировке в сейсмограммы общей точки возбуждения и осуществляют построение изображения с помощью одного из известных алгоритмов миграции.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2018103623A RU2674419C1 (ru) | 2018-01-30 | 2018-01-30 | Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2018103623A RU2674419C1 (ru) | 2018-01-30 | 2018-01-30 | Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2674419C1 true RU2674419C1 (ru) | 2018-12-07 |
Family
ID=64603857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2018103623A RU2674419C1 (ru) | 2018-01-30 | 2018-01-30 | Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2674419C1 (ru) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909655A (en) * | 1995-06-10 | 1999-06-01 | Elf Aquitaine Production | Method for processing reflection seismic traces recorded for variable offsets |
US5970023A (en) * | 1999-02-19 | 1999-10-19 | Baker Hughes Incorporated | Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms |
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
-
2018
- 2018-01-30 RU RU2018103623A patent/RU2674419C1/ru active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909655A (en) * | 1995-06-10 | 1999-06-01 | Elf Aquitaine Production | Method for processing reflection seismic traces recorded for variable offsets |
US5970023A (en) * | 1999-02-19 | 1999-10-19 | Baker Hughes Incorporated | Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms |
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
Non-Patent Citations (5)
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Iterative deblending of simultaneous-source seismic data using seislet-domain shaping regularization | |
Kühl et al. | Least-squares wave-equation migration for AVP/AVA inversion | |
US9726771B1 (en) | Time migration diffraction imaging systems and methods | |
Wang et al. | Simultaneous reverse time migration of primaries and free-surface related multiples without multiple prediction | |
US9625593B2 (en) | Seismic data processing | |
US11249207B2 (en) | Method of redatuming geophysical data | |
Silvestrov et al. | Poststack diffraction imaging using reverse‐time migration | |
Liu et al. | Reverse time migration of internal multiples for subsalt imaging | |
EA017379B1 (ru) | Способ идентификации и удаления многократных отражений для построения изображения с помощью пучков | |
Staring et al. | Three‐dimensional Marchenko internal multiple attenuation on narrow azimuth streamer data of the Santos basin, Brazil | |
Wang et al. | Waveform inversion based on wavefield decomposition | |
Lin et al. | Accurate diffraction imaging for detecting small-scale geologic discontinuities | |
Protasov et al. | 3D diffraction imaging of fault and fracture zones via image spectral decomposition of partial images | |
GB2499096A (en) | Simultaneous joint estimation of P-P and P-S residual statics | |
Dafni et al. | Scattering and dip angle decomposition based on subsurface offset extended wave-equation migration | |
Métivier et al. | A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform | |
CA3117329A1 (en) | Seismic random noise attenuation | |
RU2674419C1 (ru) | Способ обнаружения неструктурных элементов геологического разреза по сейсмограммам общего выноса | |
van der Neut et al. | Point‐spread functions for interferometric imaging | |
AU2019243987B2 (en) | Wavefield propagator for tilted orthorhombic media | |
Mildner et al. | Target-oriented velocity analysis using Marchenko-redatumed data | |
Wang et al. | Iterative deblending using unsupervised learning with double-deep neural networks | |
Li et al. | Angle domain common image gathers from reverse time migration by combining the Poynting vector with directional decomposition | |
Popov et al. | Reverse time migration with Gaussian beams and velocity analysis applications | |
Zheng et al. | Eliminating artifacts in migration of surface-related multiples: An application to marine data |