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
Links
- 230000000875 corresponding Effects 0.000 claims 4
- 238000005315 distribution function Methods 0.000 claims 3
- 238000010276 construction Methods 0.000 claims 2
- 238000001228 spectrum Methods 0.000 claims 1
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. Способ обработки сейсмограмм, включающий:
(а) сравнение допустимого порогового значения амплитудной характеристики во временном окне опорной трассы с амплитудной характеристикой анализируемой трассы в этом же временном окне; и
(б) использование не равного нулю скалярного коэффициента для обработки анализируемой трассы в данном временном окне, если ее амплитудные характеристики в этом окне выходят за пределы допустимого порогового значения амплитудной характеристики, причем применение не равного нулю скалярного коэффициента изменяет анализируемую трассу в данном временном окне таким образом, чтобы амплитуда анализируемой трассы в этом временном окне была меньше амплитуды опорной трассы;
(в) запоминание значений скалярного коэффициента вместе с соответствующей информацией о временных окнах и сейсмических трассах, где они были использованы; и
(г) восстановление первоначальных данных из обработанных сейсмограмм с использованием хранящихся в памяти значений скалярного коэффициента.
(а) сравнение допустимого порогового значения амплитудной характеристики во временном окне опорной трассы с амплитудной характеристикой анализируемой трассы в этом же временном окне; и
(б) использование не равного нулю скалярного коэффициента для обработки анализируемой трассы в данном временном окне, если ее амплитудные характеристики в этом окне выходят за пределы допустимого порогового значения амплитудной характеристики, причем применение не равного нулю скалярного коэффициента изменяет анализируемую трассу в данном временном окне таким образом, чтобы амплитуда анализируемой трассы в этом временном окне была меньше амплитуды опорной трассы;
(в) запоминание значений скалярного коэффициента вместе с соответствующей информацией о временных окнах и сейсмических трассах, где они были использованы; и
(г) восстановление первоначальных данных из обработанных сейсмограмм с использованием хранящихся в памяти значений скалярного коэффициента.
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)
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)
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 |
-
1996
- 1996-10-30 US US08/741,313 patent/US5818795A/en not_active Expired - Fee Related
-
1997
- 1997-10-22 CN CNB971984409A patent/CN1199055C/zh not_active Expired - Fee Related
- 1997-10-22 EP EP97946329A patent/EP0938687B1/en not_active Expired - Lifetime
- 1997-10-22 WO PCT/US1997/019477 patent/WO1998019182A1/en active IP Right Grant
- 1997-10-22 RU RU99106236/28A patent/RU2179732C2/ru not_active IP Right Cessation
- 1997-10-22 BR BRPI9712702-7A patent/BR9712702B1/pt not_active IP Right Cessation
- 1997-10-22 AU AU51524/98A patent/AU732854B2/en not_active Ceased
- 1997-10-29 EG EG113797A patent/EG21199A/xx active
-
1999
- 1999-03-31 SE SE9901174A patent/SE9901174D0/xx not_active Application Discontinuation
- 1999-04-29 NO NO19992087A patent/NO992087L/no not_active Application Discontinuation
-
2000
- 2000-06-05 HK HK00103408A patent/HK1024302A1/xx not_active IP Right Cessation
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) | 一种地震品质因子估算方法、装置、设备及存储介质 |