RU99106236A - Способ снижения уровня помех на сейсмических трассах - Google Patents

Способ снижения уровня помех на сейсмических трассах

Info

Publication number
RU99106236A
RU99106236A RU99106236/28A RU99106236A RU99106236A RU 99106236 A RU99106236 A RU 99106236A RU 99106236/28 A RU99106236/28 A RU 99106236/28A RU 99106236 A RU99106236 A RU 99106236A RU 99106236 A RU99106236 A RU 99106236A
Authority
RU
Russia
Prior art keywords
time window
amplitude characteristic
amplitude
trace
threshold value
Prior art date
Application number
RU99106236/28A
Other languages
English (en)
Other versions
RU2179732C2 (ru
Inventor
Саймон Райс БАРНС
Кейт Силва ХОУКИНС
Грегори Петер Гвин ФУКС
Original Assignee
Пгс Тенсор (Юс) Инк.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US08/741,313 external-priority patent/US5818795A/en
Application filed by Пгс Тенсор (Юс) Инк. filed Critical Пгс Тенсор (Юс) Инк.
Publication of RU99106236A publication Critical patent/RU99106236A/ru
Application granted granted Critical
Publication of RU2179732C2 publication Critical patent/RU2179732C2/ru

Links

Claims (40)

1. Способ уменьшения помех на сейсмических записях включающий:
(а) сравнение допустимого порогового значения амплитудной характеристики во временном окне опорной трассы с амплитудной характеристикой анализируемой трассы сейсмограммы в этом же временном окне; и
(б) использование не равного нулю скалярного коэффициента для обработки анализируемой трассы в данном временном окне, если ее амплитудные характеристики в этом окне выходят за пределы допустимого порогового значения амплитудной характеристики, при этом применение не равного нулю скалярного коэффициента приводит амплитуду анализируемой трассы в данном временном окне к значению, меньшему амплитуды опорной трассы и выше нуля, а не равный нулю скалярный коэффициент изменяется в пределах временного окна.
2. Способ по п. 1, отличающийся тем, что не равный нулю скалярный коэффициент остается постоянным при обработке записей внутри временного окна.
3. Способ по п. 1, отличающийся тем, что применение скалярного коэффициента производят с учетом функции распределения, центрированной относительно максимума, равного среднеквадратичному значению в этом временном окне.
4. Способ по п. 1, отличающийся тем, что сейсмограмма содержит трассы с равными расстояниями источник-приемник.
5. Способ по п. 1, отличающийся тем, что сейсмограмма содержит трассы с общим каналом.
6. Способ по п. 1, отличающийся тем, что сейсмограмма содержит трассы с общей средней точкой.
7. Способ по п. 1, отличающийся тем, что сейсмограмма содержит трассы с общим пунктом приема.
8. Способ по п. 1, отличающийся тем, что допустимое пороговое значение амплитудной характеристики является функцией средних пиковых значений амплитуды во временном окне опорной трассы.
9. Способ по п. 1, отличающийся тем, что пороговое значение амплитудной характеристики является функцией среднеквадратичного значения амплитуд во временном окне опорной трассы.
10. Способ по п. 1. отличающийся тем, что пороговое значение амплитудной характеристики является функцией спрямленных средних значений амплитуд во временном окне опорной трассы.
11. Способ по п. 1, отличающийся тем, что пороговое значение амплитудной характеристики составляет от 120% до 500% амплитуды во временном окне опорной трассы.
12. Способ по п. 1, отличающийся тем, что пороговое значение амплитудной характеристики является функцией амплитуд во временном окне опорной трассы в частотной полосе между 10 Гц и 200 Гц.
13. Способ по п. 12, отличающийся тем, что пороговое значение амплитудной характеристики является функцией амплитуд во временном окне опорной трассы в частотной полосе от 40 Гц до 200 Гц.
14. Способ по п. 1, отличающийся тем, что временное окно опорной трассы меньше или равно длительности последовательности из импульсов помехи.
15. Способ по п. 1, отличающийся тем, что сравнение включает определение разности амплитудных характеристик во временном окне опорной и анализируемой трассы, при этом амплитудная характеристика во временном окне опорной трассы приблизительно равна среднему значению амплитудных характеристик всех трасс сейсмограммы в этом временном окне.
16. Способ по п. 1, отличающийся тем, что сравнение включает определение разности амплитудных характеристик во временном окне опорной и анализируемой трасс, при этом амплитудная характеристика во временном окне опорной трассы приблизительно равна альфа-тримм медиане амплитудных характеристик всех трасс сейсмограммы в данном временном окне.
17. Способ по п. 1, отличающийся тем, что сравнение включает определение разности амплитудных характеристик во временном окне опорной и анализируемой трасс. при этом амплитудная характеристика во временном окне опорной трассы практически равна альфа-тримм среднему значению амплитудных характеристик всех трасс сейсмограммы в данном временном окне.
18. Способ по п. 1, отличающийся тем, что сравнение включает определение разности амплитудных характеристик во временном окне опорной и анализируемой трасс, при этом амплитудная характеристика во временном окне опорной трассы приблизительно равна медианному значению амплитудных характеристик всех трасс сейсмограммы в данном временном окне.
19. Способ по любому из пп. 15 - 18, отличающийся тем, что сравнение проводится в частотной полосе между 50 Гц и 150 Гц.
20. Способ по любому из пп. 15 - 18, отличающийся тем, что сравнение проводится в частотной полосе между 5 Гц и 40 Гц.
21. Способ по п. 1, отличающийся тем, что не равный нулю скалярный коэффициент является по существу мультипликативной функцией отношения амплитудных характеристик опорной и анализируемой трасс в данном временном окне.
22. Способ по п. 1, отличающийся тем, что не равный нулю скалярный коэффициент является по существу функцией отношения амплитудных характеристик опорной и анализируемой трасс в данном временном окне, причем это отношение возводят в степень.
23. Способ по п. 22, отличающийся тем, что используют часть сейсмических трасс на сейсмограмме с одинаковыми расстояниями источник-приемник, расположенных в интервале профиля, равном примерно 1,5 ширины зоны Френеля для преобладающей частоты в спектре сейсмического сигнала.
24. Способ по п. 22, отличающийся тем, что пороговая амплитудная характеристика является функцией среднеквадратичного значения амплитуд во временном окне опорной трассы.
25. Способ по п. 24, отличающийся тем, что пороговое значение амплитудной характеристики составляет от 120% до 500% среднеквадратичного значения амплитуды во временном окне опорной трассы.
26. Способ по п. 22. отличающийся тем, что пороговая амплитудная характеристика является функцией среднеквадратичного значения амплитуд в первом временном окне опорной трассы в частотной полосе между 10 Гц и 200 Гц.
27. Способ по п. 26, отличающийся тем. что пороговая амплитудная характеристика является функцией среднеквадратичного значения амплитуды в первом временном окне опорной трассы в частотной полосе между 40 Гц и 200 Гц.
28. Способ по п. 1 или 22, отличающийся тем, что первое временное окно опорной трассы выбрано меньшим примерно 20% периода преобладающей частоты помехи.
29. Способ по п. 1 или 22, отличающийся тем, что он дополнительно включает выполнение операций (а)-(б) для всего множества временных окон на анализируемой трассе.
30. Способ по п. 29, отличающийся тем, что временные окна перекрываются.
31. Способ по п. 30, отличающийся тем, что перекрытие множества окон выбирают в пределах от 20% до 80% длительности окна.
32. Способ по п. 1, отличающийся тем, что не равный нулю скалярный коэффициент является, в сущности, степенью отношения корректирующего скалярного коэффициента, который уравнивает амплитудные характеристики анализируемой и опорной трасс в данном временном окне. к пороговому значению амплитудной характеристики.
33. Способ по п. 1 или 29, отличающийся тем, что он включает дополнительно: получение новой сейсмограммы с измененными параметрами, основанной на применении операций (а)-(б) к первоначальной сейсмограмме, и выполнение операций (а)-(б) для новой сейсмограммы не менее одного раза.
34. Способ по п. 1, отличающийся тем, что он включает дополнительно: выполнение стандартной обработки сейсмограмм до суммирования трасс; суммирование трасс; и выполнение стандартной последующей обработки просуммированных трасс.
35. Способ по п. 1, отличающийся тем, что обработку производят для всего множества временных окон на анализируемой трассе и эта обработка включает: вычисление поправочного коэффициента для каждого временного окна, который позволяет изменить анализируемую трассу в данном временном окне таким образом, чтобы амплитудная характеристика анализируемой трассы в этом временном окне согласовывалась с амплитудной характеристикой соответствующего окна опорной трассы; использование поправочного скалярного коэффициента, равного 1, там, где поправочный коэффициент превышает допустимое пороговое значение амплитудной характеристики плюс постоянная величина, при этом определяются окна, свободные от помех; возведение поправочного скалярного коэффициента в степень, там, где он ниже допустимого порогового значения амплитудной характеристики, при этом определяют не равный единице поправочный скалярный коэффициент, и выделяют окна, содержащие помехи; проведение линейной интерполяции между поправочным скалярным коэффициентом, используемым для окна, смежного с окном, содержащим помехи, и не равным единице поправочным скалярным коэффициентом; проведение линейной интерполяции между соседними окнами, содержащими помехи, и имеющими не равные единице поправочные скалярные коэффициенты; проведение между окнами, содержащими помехи, и смежными с ними окнами без помех, параболической интерполяции между 1 и допустимым пороговым значением амплитудной характеристики для всех точек линейной интерполяции между допустимым пороговым значением амплитудной характеристики плюс постоянная величина и допустимым пороговым значением амплитудной характеристики; применение между окнами, содержащими помехи, и смежными окнами без помех линейной интерполяции между допустимым пороговым значением амплитудной характеристики и не равным единице поправочным скалярным коэффициентом, для всех точек линейной интерполяции между допустимым пороговым значением амплитудной характеристики и не равным единице поправочным скалярным коэффициентом.
36. Способ по п. 1, отличающийся тем, что обработку производят для всего множества временных окон на анализируемой трассе и эта обработка включает: вычисление поправочного коэффициента для каждого временного окна, который позволяет изменить анализируемую трассу в данном временном окне таким образом, чтобы амплитудная характеристика анализируемой трассы в этом временном окне согласовывалось с амплитудной характеристикой соответствующего окна опорной трассы; использование поправочного скалярного коэффициента, равного 1, там, где поправочный коэффициент превышает установленное пороговое значение амплитудной характеристики плюс постоянная величина, при этом определяются окна, свободные от помех; деление поправочного коэффициента для окна с помехами на допустимое пороговое значение амплитудной характеристики, при этом определяют нормализованный поправочный коэффициент для окон с помехами; применение нормализованного коэффициента при обработке сигналов в окнах с помехами; применение-интерполяции между поправочным скалярным коэффициентом, используемым в окнах, смежных с окнами с помехами, и не равным единице, с поправочным скалярным коэффициентом.
37. Способ по п. 36, отличающийся тем, что перед применением нормализованного поправочного коэффициента его возводят в степень.
38. Способ по п. 1, отличающийся тем, что он дополнительно включает определение допустимого порогового значения амплитудной характеристики в интервале между статистическими параметрами амплитудной характеристики сигнала и статистическими параметрами, определяющими помеху.
39. Способ по п. 38, отличающийся тем, что определение допустимого порогового значения амплитудной характеристики включает: вычисление поправочного коэффициента для каждого временного окна, который позволяет изменить анализируемую трассу в данном временном окне таким образом, чтобы амплитудная характеристика анализируемом трассы в этом временном окне согласовывалось с амплитудной характеристикой соответствующего окна опорной трассы; построение функции распределения поправочных коэффициентов, связанных с полезным сигналом; построение функции распределения поправочных коэффициентов, связанных с помехами: и выбор порога для допустимого порогового значения амплитудной характеристики в интервале между распределениями полезного сигнала и помех.
40. Способ обработки сейсмограмм, включающий:
(а) сравнение допустимого порогового значения амплитудной характеристики во временном окне опорной трассы с амплитудной характеристикой анализируемой трассы в этом же временном окне; и
(б) использование не равного нулю скалярного коэффициента для обработки анализируемой трассы в данном временном окне, если ее амплитудные характеристики в этом окне выходят за пределы допустимого порогового значения амплитудной характеристики, причем применение не равного нулю скалярного коэффициента изменяет анализируемую трассу в данном временном окне таким образом, чтобы амплитуда анализируемой трассы в этом временном окне была меньше амплитуды опорной трассы;
(в) запоминание значений скалярного коэффициента вместе с соответствующей информацией о временных окнах и сейсмических трассах, где они были использованы; и
(г) восстановление первоначальных данных из обработанных сейсмограмм с использованием хранящихся в памяти значений скалярного коэффициента.
RU99106236/28A 1996-10-30 1997-10-22 Способ снижения уровня помех на сейсмических трассах RU2179732C2 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/741,313 US5818795A (en) 1996-10-30 1996-10-30 Method of reduction of noise from seismic data traces
US08/741,313 1996-10-30

Publications (2)

Publication Number Publication Date
RU99106236A true RU99106236A (ru) 2001-01-20
RU2179732C2 RU2179732C2 (ru) 2002-02-20

Family

ID=24980213

Family Applications (1)

Application Number Title Priority Date Filing Date
RU99106236/28A RU2179732C2 (ru) 1996-10-30 1997-10-22 Способ снижения уровня помех на сейсмических трассах

Country Status (11)

Country Link
US (1) US5818795A (ru)
EP (1) EP0938687B1 (ru)
CN (1) CN1199055C (ru)
AU (1) AU732854B2 (ru)
BR (1) BR9712702B1 (ru)
EG (1) EG21199A (ru)
HK (1) HK1024302A1 (ru)
NO (1) NO992087L (ru)
RU (1) RU2179732C2 (ru)
SE (1) SE9901174D0 (ru)
WO (1) WO1998019182A1 (ru)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6016462A (en) * 1997-08-29 2000-01-18 Exxon Production Research Company Analysis of statistical attributes for parameter estimation
US6148264A (en) * 1998-07-06 2000-11-14 Exxonmobil Upstream Research Company Method for removing seismic noise caused by external activity
GB0015974D0 (en) * 2000-06-30 2000-08-23 Geco As Quality control of data
US6751559B2 (en) 2002-09-10 2004-06-15 Pgs Exploration (Uk) Limited Method for suppressing noise from seismic signals by source position determination
US7778110B2 (en) * 2003-03-26 2010-08-17 Westerngeco L.L.C. Processing seismic data representative of the acceleration wavefield
US6882938B2 (en) * 2003-07-30 2005-04-19 Pgs Americas, Inc. Method for separating seismic signals from two or more distinct sources
US7869303B2 (en) * 2007-08-14 2011-01-11 Pgs Geophysical As Method for noise suppression in seismic signals using spatial transforms
US7675812B2 (en) * 2008-06-30 2010-03-09 Pgs Geophysical As Method for attenuation of multiple reflections in seismic data
US8811115B2 (en) * 2008-08-14 2014-08-19 Pgs Geophysical As Attenuating seismic interference noise using a dual sensor recording system
US8274858B2 (en) 2009-11-12 2012-09-25 Pgs Geophysical As Method for full-bandwidth deghosting of marine seismic streamer data
US20110182142A1 (en) * 2010-01-27 2011-07-28 Qinglin Liu Technique and Apparatus for Seismic Data Quality Control
US8612158B2 (en) 2010-05-06 2013-12-17 Westerngeco L.L.C. Seismic interference noise elimination
EP2703847A1 (en) * 2011-02-25 2014-03-05 Cuadra Canales, Francisco Method of dynamic seismic prospecting
CN103336940B (zh) * 2011-10-26 2016-12-21 中国石油化工股份有限公司 一种数字信号同相轴识别方法及设备
US9541659B2 (en) 2011-11-18 2017-01-10 Westerngeco L.L.C. Noise removal from 3D seismic representation
CN102879825B (zh) * 2012-09-18 2015-01-07 吉林大学 可控震源地震数据的强脉冲噪声检测及压制方法
US9360575B2 (en) 2013-01-11 2016-06-07 Fairfield Industries Incorporated Simultaneous shooting nodal acquisition seismic survey methods
US9310504B2 (en) * 2013-02-01 2016-04-12 Pgs Geophysical As Systems and methods for detecting swell noise in a seismic gather
US9329293B2 (en) 2013-03-12 2016-05-03 Pgs Geophysical, As Systems and methods for removing acquisition related effects from seismic data
US9274241B2 (en) 2013-03-14 2016-03-01 Pgs Geophysical As Method and system for suppressing swell-induced electromagnetic noise
CN104597501B (zh) * 2013-11-01 2017-07-07 中国石油天然气集团公司 炮集记录处理方法及炮集记录处理装置
US11391857B2 (en) 2013-12-30 2022-07-19 Pgs Geophysical As Methods and systems for attenuating residual acoustic energy in seismic data
US10429531B2 (en) * 2015-08-24 2019-10-01 Shreya Biswas Ley Advanced noise reduction in acoustic well logging
CN107589449B (zh) * 2017-08-29 2020-04-28 电子科技大学 基于曲线Gabor滤波的三维数据断层增强方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2733412A (en) * 1956-01-31 Seismic amplifier interference eliminator
US3344395A (en) * 1965-02-01 1967-09-26 Pan American Petroleum Corp Displaying subtractively combined seismic data
US3611279A (en) * 1969-10-13 1971-10-05 Phillips Petroleum Co Noise elimination from seismic signals
US3704444A (en) * 1970-06-01 1972-11-28 Western Geophysical Co Seismic data processing method and system
US4203161A (en) * 1972-03-01 1980-05-13 Texaco Inc. Method of enhancing common depth point seismic data
US4204279A (en) * 1972-03-01 1980-05-20 Texaco Inc. Method for enhancing seismic data
US4210968A (en) * 1975-12-16 1980-07-01 Lindseth Roy O Seismic exploration technique
NO170441C (no) * 1985-05-06 1992-10-14 Western Atlas Int Inc Fremgangsmaate for undertrykkelse av koherent stoey i marine seismikk-data
US4937794A (en) * 1985-05-06 1990-06-26 Western Atlas International, Inc. Seismic noise suppression method
US4707812A (en) * 1985-12-09 1987-11-17 Atlantic Richfield Company Method of suppressing vibration seismic signal correlation noise
US4882713A (en) * 1988-09-06 1989-11-21 Exxon Production Research Company Method for noise suppression in the stacking of seismic traces
US4910716A (en) * 1989-01-31 1990-03-20 Amoco Corporation Suppression of coherent noise in seismic data
NO166902C (no) * 1989-03-17 1991-09-11 Geco As Fremgangsmaate ved sanntids kontrollsystem for seismisk interferens.
US5182729A (en) * 1991-09-03 1993-01-26 Exxon Production Research Company Reduction of sideswipe noise from seismic data by null steering
US5237538A (en) * 1992-02-20 1993-08-17 Mobil Oil Corporation Method for removing coherent noise from an array of seismic traces
FR2698697B1 (fr) * 1992-12-02 1995-02-17 Geophysique Cie Gle Procédé de réduction des bruits industriels en sismique marine, et enregistrement obtenu par ce procédé.
GB2273359B (en) * 1992-12-12 1997-01-15 Schlumberger Ltd Method for improving signal to noise ratio
US5293352A (en) * 1993-01-07 1994-03-08 Western Atlas International, Inc. Method for removing noise due to near surface scatterers
US5365492A (en) * 1993-08-04 1994-11-15 Western Atlas International, Inc. Method for reverberation suppression

Similar Documents

Publication Publication Date Title
RU99106236A (ru) Способ снижения уровня помех на сейсмических трассах
Ferguson et al. A simple algorithm for band-limited impedance inversion
AU2006214552B2 (en) Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects
RU2179732C2 (ru) Способ снижения уровня помех на сейсмических трассах
US5173879A (en) Surface-consistent minimum-phase deconvolution
US8112147B2 (en) Method and apparatus for generating determination indexes for identifying ECG interfering signals
US5392213A (en) Filter for removal of coherent noise from seismic data
US20160091623A1 (en) Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data
AU739793B2 (en) Method of stacking seismic traces
US4964098A (en) Method for seismic trace interpolation
US4561074A (en) Computationally efficient weighting and vertical stacking methods and apparatus for improving signal-to-noise ratio of seismic data
Buttkus Homomorphic filtering—theory and practice
Evans et al. Design of a sanguine noise processor based upon world-wide extremely low frequency (ELF) recordings
US4995007A (en) Method for processing seismic data
CN104502965A (zh) 一种振幅补偿因子的反演方法
US4688198A (en) Entropy guided deconvolution of seismic signals
MY107350A (en) A method in a real time control system for seismic interference.
US6266620B1 (en) Detection of ground roll cone
EP0354813A3 (en) Seismic data processing
Hampson The discrete Radon transform: a new tool for image enhancement and noise suppression
CN110865410A (zh) 一种基于nar-tfpf压制地震勘探随机噪声的方法
Anderson et al. Automatic editing of noisy seismic data 1
USH1693H (en) Method for removing non-geologic amplitude variations from seismic data
Wenz Some periodic variations in low‐frequency acoustic ambient noise levels in the ocean
CN112782763B (zh) 一种地震品质因子估算方法、装置、设备及存储介质