RU2010135895A - Преобразование радона волнового фронта - Google Patents

Преобразование радона волнового фронта Download PDF

Info

Publication number
RU2010135895A
RU2010135895A RU2010135895/28A RU2010135895A RU2010135895A RU 2010135895 A RU2010135895 A RU 2010135895A RU 2010135895/28 A RU2010135895/28 A RU 2010135895/28A RU 2010135895 A RU2010135895 A RU 2010135895A RU 2010135895 A RU2010135895 A RU 2010135895A
Authority
RU
Russia
Prior art keywords
data array
multidimensional data
collection
processing
parameters
Prior art date
Application number
RU2010135895/28A
Other languages
English (en)
Other versions
RU2460095C2 (ru
Inventor
Алекс БЕРКОВИЧ (IL)
Алекс БЕРКОВИЧ
Игорь БЕЛЬФЕР (IL)
Игорь БЕЛЬФЕР
Роман КАЗИННИК (IL)
Роман КАЗИННИК
Original Assignee
Геомейдж (2003) Лтд (Il)
Геомейдж (2003) Лтд
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Геомейдж (2003) Лтд (Il), Геомейдж (2003) Лтд filed Critical Геомейдж (2003) Лтд (Il)
Publication of RU2010135895A publication Critical patent/RU2010135895A/ru
Application granted granted Critical
Publication of RU2460095C2 publication Critical patent/RU2460095C2/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/56De-ghosting; Reverberation compensation

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

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

Claims (38)

1. Способ обработки данных, выполняемый на компьютере, предусматривающий следующие шаги:
- берут первую коллекцию трасс, показывающих изменение во времени принимаемых сигналов, обусловленное отражениями сейсмических волн от подповерхностных структур;
- задают преобразование Радона относительно набора параметров фронта сейсмических волн, при этом преобразование определяет сложение амплитуд сейсмических волн вдоль траекторий, определяемых параметрами волнового фронта;
- применяют преобразование Радона к первой коллекции трасс, преобразуя ее в многомерный массив данных, определяемый как функция по меньшей мере двух из параметров волнового фронта;
- обработкой многомерного массива данных получают вторую коллекцию трасс, качество изображения в которой выше, чем в первой коллекции;
- обработкой второй коллекции трасс строят сейсмический снимок подповерхностных структур с повышенным качеством изображения.
2. Способ по п.1, отличающийся тем, что преобразование, обработку многомерного массива данных и второй коллекции производят безотносительно к расположению первых точек, в которых порождались сейсмические волны, вторых точек, в которых принимались сигналы, а также средних между первыми и вторыми точек.
3. Способ по п.1, отличающийся тем, что многомерный массив данных имеет три или более измерения, отвечающие по меньшей мере двум из параметров волнового фронта и времени пробега сейсмических волн.
4. Способ по любому из пп.1-3, отличающийся тем, что в число параметров волнового фронта входят угловой параметр и по меньшей мере один радиус кривизны фронта волны, исходящей от подповерхностного отражателя.
5. Способ по п.4, отличающийся тем, что угловой параметр есть угол выхода волны, а по меньшей мере один радиус кривизны представляет собой соответственно радиус волнового фронта с отражением от общей точки или с отражением от общего участка.
6. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных предусмотрена коррекция искажений, внесенных в массив в результате применения преобразования Радона.
7. Способ по п.6, отличающийся тем, что для коррекции искажений предусмотрено частотное преобразование многомерного массива данных с коррекцией искажений в частотных координатах.
8. Способ по п.6, отличающийся тем, что искажения корректируют в пространственно-временных координатах.
9. Способ по п.6, отличающийся тем, что коррекция искажений предполагает итеративную корректирующую обработку многомерного массива данных.
10. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных различают первые и вторые сейсмические события, которые трансформируются соответственно в первые и вторые области массива.
11. Способ по п.10, отличающийся тем, что разделение первого и второго событий предусматривает обнуление первой или второй области в многомерном массиве.
12. Способ по п.10, отличающийся тем, что первое сейсмическое событие вызвано однократным отражением, а второе - многократными отражениями.
13. Способ по п.10, отличающийся тем, что для разделения первого и второго событий задают первую и вторую области в массиве как функцию по меньшей мере двух параметров волнового фронта.
14. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных предусмотрена фильтрация по меньшей мере части данных массива.
15. Способ по любому из пп.1-3, отличающийся тем, что при обработке массива предусмотрена интерполяция данных в нем для увеличения числа трасс во второй коллекции по сравнению с первой.
16. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных предусмотрена компенсация различия скоростей сейсмических волн.
17. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных предусмотрено применение к нему миграции.
18. Способ по любому из пп.1-3, отличающийся тем, что при обработке второй коллекции предусмотрено интегрирование входящих в нее трасс для получения композитной трассы, представляющей подповерхностные структуры.
19. Способ по любому из пп.1-3, отличающийся тем, что при обработке многомерного массива данных предусмотрено применение к нему обратного преобразования Радона для получения второй коллекции трасс.
20. Устройство обработки данных, содержащее:
интерфейс, предназначенный для приема первой коллекции трасс, показывающих изменение во времени принимаемых сигналов, обусловленное отражениями сейсмических волн от подповерхностных структур;
процессор сигналов, предназначенный определять метрику преобразования Радона относительно набора параметров фронта сейсмических волн, причем преобразование определяет сложение амплитуд сейсмических волн вдоль траекторий, определяемых параметрами волнового фронта,
а также предназначенный применять преобразование Радона к первой коллекции трасс, преобразуя ее в многомерный массив данных, определяемый как функция по меньшей мере двух из параметров волнового фронта,
а также предназначенный обрабатывать многомерный массив данных для получения второй коллекции трасс, качество изображения в которой выше, чем в первой коллекции,
а также предназначенный обрабатывать вторую коллекцию трасс для построения сейсмического снимка подповерхностных структур с повышенным качеством изображения.
21. Устройство по п.20, отличающееся тем, что процессор сигналов предназначен производить преобразование, обработку многомерного массива данных и второй коллекции безотносительно к расположению первых точек, в которых порождались сейсмические волны, вторых точек, в которых принимались сигналы, а также средних между первыми и вторыми точек.
22. Устройство по п.20, отличающееся тем, что многомерный массив данных имеет три или более измерения, отвечающие по меньшей мере двум из параметров волнового фронта и времени пробега сейсмических волн.
23. Устройство по любому из пп.20-22, отличающееся тем, что в число параметров волнового фронта входят угловой параметр и по меньшей мере один радиус кривизны фронта волны, исходящей от подповерхностного отражателя.
24. Устройство по п.23, отличающееся тем, что угловой параметр есть угол выхода волны, а по меньшей мере один радиус кривизны представляет собой соответственно радиус волнового фронта с отражением от общей точки или с отражением от общего участка.
25. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен производить коррекцию искажений, внесенных в многомерный массив данных в результате применения преобразования Радона.
26. Устройство по п.25, отличающееся тем, что процессор сигналов предназначен производить частотное преобразование многомерного массива данных с коррекцией искажений в частотных координатах.
27. Устройство по п.25, отличающееся тем, что процессор сигналов предназначен производить коррекцию искажений в пространственно-временных координатах.
28. Устройство по п.25, отличающееся тем, что процессор сигналов предназначен производить коррекцию искажений путем итеративной корректирующей обработки многомерного массива данных.
29. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен различать первые и вторые сейсмические события, которые трансформируются соответственно в первые и вторые области многомерного массива данных.
30. Устройство по п.29, отличающееся тем, что процессор сигналов предназначен обнулять первую или вторую область в многомерном массиве данных.
31. Устройство по п.29, отличающееся тем, что первое сейсмическое событие вызвано однократным отражением, а второе - многократными отражениями.
32. Устройство по п.29, отличающееся тем, что процессор сигналов предназначен задавать первую и вторую области в многомерном массиве данных как функцию по меньшей мере двух параметров волнового фронта.
33. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен производить фильтрацию по меньшей мере части данных многомерного массива.
34. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен производить интерполяцию данных в многомерном массиве данных для увеличения числа трасс во второй коллекции по сравнению с первой.
35. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен путем обработки многомерного массива данных производить компенсацию различия скоростей сейсмических волн.
36. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен производить миграцию в многомерном массиве данных.
37. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен интегрировать входящие во вторую коллекцию трассы для получения композитной трассы, представляющей подповерхностные структуры.
38. Устройство по любому из пп.20-22, отличающееся тем, что процессор сигналов предназначен при обработке многомерного массива данных применять к нему обратное преобразование Радона для получения второй коллекции трасс.
RU2010135895/28A 2008-04-01 2009-02-22 Преобразование радона волнового фронта RU2460095C2 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/080,218 2008-04-01
US12/080,218 US8209126B2 (en) 2008-04-01 2008-04-01 Wavefront-defined Radon transform

Publications (2)

Publication Number Publication Date
RU2010135895A true RU2010135895A (ru) 2012-05-10
RU2460095C2 RU2460095C2 (ru) 2012-08-27

Family

ID=41118413

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2010135895/28A RU2460095C2 (ru) 2008-04-01 2009-02-22 Преобразование радона волнового фронта

Country Status (4)

Country Link
US (1) US8209126B2 (ru)
EP (1) EP2260328A4 (ru)
RU (1) RU2460095C2 (ru)
WO (1) WO2009122392A2 (ru)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2692947A1 (en) * 2007-07-06 2009-01-15 Schlumberger Canada Limited Methods and systems for processing microseismic data
US7630865B2 (en) * 2007-09-11 2009-12-08 Geomage (2003) Ltd Complex analysis of kinematics for non-hyperbolic moveout corrections
US8880348B2 (en) * 2010-09-21 2014-11-04 Baker Hughes Incorporated Radon migration of acoustic data
CN102053271B (zh) * 2010-11-17 2012-09-26 中国石油天然气股份有限公司 地震采集质量分析方法及装置
US9103935B2 (en) 2013-02-04 2015-08-11 King Fahd University Of Petroleum And Minerals Method of first arrival picking of seismic refraction data
MX2016001025A (es) * 2013-07-23 2016-10-05 Cgg Services Sa Metodo para la designacion de datos sismicos adquiridos utilizando una fuente en movimiento.
US20150198729A1 (en) * 2014-01-13 2015-07-16 Cgg Services Sa Regularization of spatially aliased seismic data
US20170248716A1 (en) * 2014-11-14 2017-08-31 Cgg Services Sas Device and method for weighted sparse inversion for seismic processing
US10032077B1 (en) 2015-10-29 2018-07-24 National Technology & Engineering Solutions Of Sandia, Llc Vehicle track identification in synthetic aperture radar images
CN106597540B (zh) * 2016-12-30 2017-10-31 中国科学院地质与地球物理研究所 高斯束偏移成像方法及装置
US11016212B2 (en) 2017-04-11 2021-05-25 Saudi Arabian Oil Company Compressing seismic wavefields in three-dimensional reverse time migration
CN108181648B (zh) * 2017-11-23 2019-07-09 中国石油天然气集团公司 连续记录混叠地震采集数据的质控方法和装置
US10684382B2 (en) * 2018-01-23 2020-06-16 Saudi Arabian Oil Company Generating target-oriented acquisition-imprint-free prestack angle gathers using common focus point operators
US11561312B2 (en) * 2019-12-16 2023-01-24 Saudi Arabian Oil Company Mapping near-surface heterogeneities in a subterranean formation
US11353609B2 (en) 2019-12-20 2022-06-07 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using angle domain gathers sampled in a spiral coordinate space
CN111220158B (zh) * 2020-01-15 2022-09-09 哈尔滨工程大学 一种基于时间方位历程图的线谱目标运动参数估计方法
US11656378B2 (en) 2020-06-08 2023-05-23 Saudi Arabian Oil Company Seismic imaging by visco-acoustic reverse time migration
CN113687416B (zh) * 2021-08-19 2022-08-30 吉林大学 一种针对倾斜金属矿体成像方法

Family Cites Families (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4760563A (en) 1986-01-09 1988-07-26 Schlumberger Technology Corporation Seismic exploration using exactly invertible discrete transformation into tau-p space
IL83306A (en) 1987-07-23 1990-12-23 Gelchinsky Boris Method of analysing the structure of a medium,particularly useful for seismic prospecting
EP0331585B1 (en) * 1988-03-04 1993-12-29 Societe De Prospection Electrique Schlumberger Method for deconvolution of unknown source signatures from unknown waveform data
US4935905A (en) 1989-08-04 1990-06-19 Terra Linda Canada, Inc. Method and apparatus for seismic exploration
IL92132A (en) 1989-10-27 1994-07-31 Gelchinsky Boris Homeomorphic character method for analyzing the structure of the medium
DE69117947D1 (de) * 1990-04-20 1996-04-18 Schlumberger Technology Bv Verfahren und vorrichtung zur rohrwellenmessung eines bohrlochs unter verwendung von diskreten frequenzen
US5383114A (en) 1992-10-05 1995-01-17 Western Atlas International, Inc. Method for displaying a volume of seismic data
US5563949A (en) 1994-12-12 1996-10-08 Amoco Corporation Method of seismic signal processing and exploration
RU2187130C2 (ru) * 1995-10-06 2002-08-10 Корэ Лэбораторис Глобал Н.В. Способ и устройство для обработки сейсмического сигнала и проведения разведки полезных ископаемых
US5831935A (en) 1996-03-05 1998-11-03 Chevron U.S.A. Inc. Method for geophysical processing and interpretation using seismic trace difference for analysis and display
US5724309A (en) 1996-03-06 1998-03-03 Chevron U.S.A. Inc. Method for geophysical processing and interpretation using instantaneous phase and its derivatives and their derivatives
CA2220274C (en) 1996-04-12 2005-06-28 Amoco Corporation Method and apparatus for seismic signal processing and exploration
US5920828A (en) 1997-06-02 1999-07-06 Baker Hughes Incorporated Quality control seismic data processing system
US5987388A (en) 1997-12-26 1999-11-16 Atlantic Richfield Company Automated extraction of fault surfaces from 3-D seismic prospecting data
US6018498A (en) 1998-09-02 2000-01-25 Phillips Petroleum Company Automated seismic fault detection and picking
US6018500A (en) * 1999-02-19 2000-01-25 Baker Hughes Incorporated Radon moveout
US6263284B1 (en) 1999-04-22 2001-07-17 Bp Corporation North America Inc. Selection of seismic modes through amplitude characteristics
US6636809B1 (en) 1999-05-03 2003-10-21 Compagnie Generale De Geophysique High resolution radon transform seismic traces processing
US6989841B2 (en) * 2001-05-29 2006-01-24 Fairfield Industries, Inc. Visualization method for the analysis of prestack and poststack seismic data
US6850845B2 (en) * 2001-07-20 2005-02-01 Tracy Joseph Stark System for multi-dimensional data analysis
US7069149B2 (en) 2001-12-14 2006-06-27 Chevron U.S.A. Inc. Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume
US6763304B2 (en) * 2002-05-21 2004-07-13 Pgs Americas, Inc. Method for processing seismic data to attenuate multiples
FR2841344B1 (fr) 2002-06-19 2005-04-29 Tsurf Procede, dispositif et produit-programme de lissage d'une propriete de subsurface
US6691039B1 (en) * 2002-08-30 2004-02-10 John M. Robinson Removal of noise from seismic data using improved radon transformations
US6904368B2 (en) 2002-11-12 2005-06-07 Landmark Graphics Corporation Seismic analysis using post-imaging seismic anisotropy corrections
NZ543753A (en) 2003-04-24 2008-11-28 Shell Int Research Thermal processes for subsurface formations
US6832161B1 (en) 2003-08-15 2004-12-14 Westerngeco, L.L.C. Method for attenuating water layer multiples
US7085195B2 (en) 2004-03-30 2006-08-01 Rdsp I, L.P. Method for detecting earth formation fractures by seismic imaging of diffractors
US7197399B2 (en) 2005-01-13 2007-03-27 Bp Corporation North America, Inc. Method of multiple attenuation
US7149630B2 (en) 2005-01-13 2006-12-12 Bp Corporation North America Inc. Method of DMO calculation for use in seismic exploration
US7561491B2 (en) * 2005-03-04 2009-07-14 Robinson John M Radon transformations for removal of noise from seismic data
US7295074B2 (en) * 2005-03-31 2007-11-13 Beceem Communications Inc. Doherty power amplifier with phase compensation
US7725266B2 (en) * 2006-05-31 2010-05-25 Bp Corporation North America Inc. System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling
US7719923B2 (en) * 2006-06-05 2010-05-18 Westerngeco L.L.C. Low frequency model estimation
WO2008111037A2 (en) 2007-03-12 2008-09-18 Geomage 2003 Ltd A method for identifying and analyzing faults/fractures using reflected and diffracted waves
US7630865B2 (en) 2007-09-11 2009-12-08 Geomage (2003) Ltd Complex analysis of kinematics for non-hyperbolic moveout corrections

Also Published As

Publication number Publication date
US8209126B2 (en) 2012-06-26
US20090248313A1 (en) 2009-10-01
RU2460095C2 (ru) 2012-08-27
WO2009122392A2 (en) 2009-10-08
EP2260328A4 (en) 2017-03-22
WO2009122392A3 (en) 2010-03-11
EP2260328A2 (en) 2010-12-15

Similar Documents

Publication Publication Date Title
RU2010135895A (ru) Преобразование радона волнового фронта
CA2597598C (en) Method and aparatus for true relative amplitude correction of seismic data for normal moveout stretch effects
RU2192658C2 (ru) Способ формирования сборок с общим удалением и общим азимутом при трехмерных сейсмических исследованиях и способ изучения вариаций характеристик отраженных сигналов
WO2008112036A1 (en) Imaging of multishot seismic data
CN109765616A (zh) 一种保幅波场延拓校正方法及系统
CN104533396A (zh) 一种远探测声波的处理方法
CN107783191B (zh) 多维空间时空时频峰值滤波消减地震勘探随机噪声的方法
CN109597021A (zh) 一种波达方向估计方法及装置
CN111175727B (zh) 一种基于条件波数谱密度的宽带信号方位估计的方法
CN106950600B (zh) 一种近地表散射面波的去除方法
EP0300477A3 (en) Method of analysing the structure of a medium, particularly useful for seismic prospecting
CN107515422B (zh) 一种高精度浅地层剖面采集系统
CN111239814B (zh) 基于同相轴分频追踪平滑的浅剖数据机械干扰压制方法
CN111257938A (zh) 基于小波互相关时移地震虚拟震源波场重构方法和系统
AU2010314924B2 (en) System and method for seismic beam formation that accounts for equipment misalignment
CN102998702A (zh) 保幅平面波叠前深度偏移方法
CN114578290A (zh) 一种基于频率修正的麦克风阵列声源定位方法
CN111766631B (zh) 一种面波压制方法及装置
WO2013159836A1 (en) A computerized method and a computer program product for determining a resulting data set representative of a geological region of interest
CN110007296B (zh) 一种基于引导信号修正的时域干扰抵消方法
CN116125535B (zh) 三维vsp成像的方法及装置
CN104793243B (zh) 基于n次根叠加的定向地震波数据处理方法
JP2544535B2 (ja) 移動音源の計測装置及び計測方法
Tang Combining source direction vectors with wavefield decomposition to calculate angle gathers from reverse time migration
Liu et al. Multiple attenuation using multichannel predictive deconvolution in radial domain