RU2308743C2 - Способ фильтрации сейсмических сигналов и устройство для его осуществления - Google Patents

Способ фильтрации сейсмических сигналов и устройство для его осуществления Download PDF

Info

Publication number
RU2308743C2
RU2308743C2 RU2004131199/28A RU2004131199A RU2308743C2 RU 2308743 C2 RU2308743 C2 RU 2308743C2 RU 2004131199/28 A RU2004131199/28 A RU 2004131199/28A RU 2004131199 A RU2004131199 A RU 2004131199A RU 2308743 C2 RU2308743 C2 RU 2308743C2
Authority
RU
Russia
Prior art keywords
signal
signals
acoustic wave
seismic
detector
Prior art date
Application number
RU2004131199/28A
Other languages
English (en)
Other versions
RU2004131199A (ru
Inventor
Гэри Э. КРУС (US)
Гэри Э. КРУС
Рассел Э. КЭППАЙУС (US)
Рассел Э. КЭППАЙУС
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
Application filed by Инпут/Аутпут, Инк. filed Critical Инпут/Аутпут, Инк.
Publication of RU2004131199A publication Critical patent/RU2004131199A/ru
Application granted granted Critical
Publication of RU2308743C2 publication Critical patent/RU2308743C2/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/181Geophones
    • G01V1/184Multi-component geophones
    • 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
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • 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
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/20Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/27Other pre-filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering
    • G01V2210/3246Coherent noise, e.g. spatially coherent or predictable

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)
  • Feedback Control In General (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

В настоящем изобретении разработан способ адаптивной фильтрации для эффективного подавления поверхностных волн, наблюдаемых при сборе данных в ходе сейсморазведки. Для осуществления способа предлагается устройство для фильтрации сейсмических сигналов для системы сбора сейсмических данных, включающее многоосный детектор и средство адаптивной фильтрации. 2 н. и 24 з.п. ф-лы, 8 ил.

Description

Область техники, к которой относится изобретение
Изобретение относится в основном к области геологической разведки и более конкретно к устройству и способу сбора и обработки сейсмических данных.
Уровень техники
Системы сейсмической разведки, использующие стандартные геофоны, акселерометры и гидрофоны, обычно содержат несколько сенсоров, генерирующих сигнал, указывающий на наличие в среде акустической энергии, обычно в форме сейсмической волны. Сейсмическая волна, как правило, генерируется энергетическим источником, таким как установленный на грузовике вибратор, взрыв или, в случае с гидрофонной системой, пневматическая пушка (пневматический источник сейсмических сигналов). Эти вызванные сейсмической волной сигналы поступают затем на схему (плату) сбора данных/преобразования сигнала. С помощью схемы сбора данных/преобразования сигнала аналоговые сигналы от одного или более удаленных сейсмических сенсоров (гидрофонов, геофонов или других сейсмоприемников) записываются и преобразовываются в ряды цифровых значений. Типичная плата сбора данных/преобразования сигнала конфигурируется таким образом, чтобы удовлетворять требуемой скорости записи данных, использовать любой цифровой фильтр или любой другой метод обработки оцифрованного сигнала, а также для проведения диагностики.
Одна или более схем сбора данных/преобразования сигнала собираются в блок сбора данных. Каждый блок сбора данных сохраняет ряды численных данных от всех сейсмических сенсоров, соединенных со схемами сбора данных/преобразования сигнала. Блок сбора данных передает информацию в системный контроллер, который может включать и сейсмозаписывающее устройство.
Сейсмические (акустические) волны содержат множество различных составляющих. Они включают Р (продольные) и S (поперечные или сдвиговые) волны, которые характеризуются различными временами вступления на детекторе (различными скоростями). Различные геологические события характеризуются сигналами разных частот. Соответственно, вертикальная и горизонтальная компоненты сигнала в конечных временных окнах не коррелируют между собой.
Другие сейсмические сигналы являются нежелательным шумом. Одна из разновидностей шума известна как "поверхностные волны" (ground roll). Поверхностные волны одновременно и с одними и теми же частотами проявляется на вертикальной и горизонтальной компонентах записи. Эти компоненты поляризованы под углом в 90° друг к другу. Таким образом, вертикальная и горизонтальная компоненты шума коррелируют друг с другом в конечном временном окне.
Поверхностные волны генерируются при включении акустического излучателя. Акустическая волна отражается от неглубоких препятствий так же, как и от границ объекта поиска на глубине. Отраженные от неглубоких препятствий волны имеют большую в сравнении с отраженными от цели на глубине амплитуду и иногда приходят к поверхностному приемнику одновременно с полезным сигналом. Многие стандартные системы имеют приемные характеристики, близкие к характеристикам шума, поэтому шум оказывает нежелательное воздействие на процесс сбора сейсмических данных. Чтобы подавить поверхностные волны, широко распространена практика монтажа группы определенно расположенных геофонов (в виде измерительной сети) для использования того факта, что частоты шума отличаются от частот полезного сигнала (по крайней мере, за счет ограничения полосы частот), кроме того, шум имеет большую величину амплитуды и дисперсии. Однако монтаж геофонов в виде измерительной сети сам по себе приводит лишь к ограниченному снижению влияния поверхностных волн. Кроме того, низкая скорость поверхностных волн предъявляет жесткие требования к пространственным характеристикам съемки при желании осуществлять подавление шума с использованием традиционных методов (цифровой) фильтрации, таких как фильтрация методом "частота - волновое число" (FK).
Поляризационные фильтры относятся к широкой категории фильтров, в которых распознавание и отделение поверхностных волн основано на том факте, что в этих волнах частицы твердого тела совершают движение по эллиптическим траекториям. Однако применение даже нескольких поляризационных фильтров дает ограниченный эффект из-за дисперсной природы поверхностных волн и из-за погрешностей в допущениях о локальной геологической структуре среды.
Адаптивные системы, как правило, характеризуются способностью к переменной во времени настройке с использованием динамической оценки шума на входе. Однако известные методы адаптивной фильтрации не включают возможности использования вертикальной и горизонтальной компонент сигнала от одного и того же многоосного сенсора. Соответственно, типичные методы фильтрации не позволяют отличить поверхностные волны от других волн и могут отфильтровать полезный сигнал вместе с поверхностными волнами. И наконец, типичные системы не дают возможности понижать влияние поверхностных волн от станции к станции (station-by-station), поскольку используют в основном пространственные данные и монтаж приемников в виде пространственной измерительной сети.
Раскрытие изобретения
Настоящее изобретение направлено на решение вышеперечисленных проблем обычных систем сбора сейсмических данных обеспечения системы, содержащей множество полевых блоков, включающих один или более сейсмических излучателей, один или более многоосных сейсмических сенсоров для регистрации акустических волн, генерируемых излучателями и отраженных подземными структурами, блок сбора данных, приемник/излучатель и процессор центрального контроллера, устанавливаемый на удалении от одного или более полевых блоков. Процессор использует алгоритм адаптивной фильтрации, позволяющий обрабатывать вертикальную и горизонтальную компоненты сигнала с сенсора с использованием начальных коэффициентов фильтра, адаптированных для оценки изменений шума. Подбор коэффициентов осуществляется опытным путем с помощью последовательного анализа выборок данных (sample-by-sample). В предпочтительном исполнении начальные коэффициенты фильтрации хранятся в запоминающем устройстве (средстве хранения) контроллера.
В соответствии с настоящим изобретением способ обработки данных включает два этапа. Первый этап заключается в создании фильтра, при котором определяются зависимые от сдвига (фаз) временные окна (окна приложений) на основе мгновенного значения разницы фаз вертикальной и горизонтальной компонент сигнала. Полоса частот может быть выбрана с использованием фильтра полосы пропускания частот. Второй этап заключается в применении адаптивного фильтрующего процесса (алгоритма) с использованием метода наименьший-средний-квадрат (LMS - least-mean-square) к данным внутри окна и к временному интервалу, заключенному в окне.
Соответственно предлагается способ фильтрации сейсмических сигналов в системе сейсморазведки, в котором принимают (регистрируют) посредством детектора акустическую (звуковую) волну, передают с детектора первый и второй сигналы, которые представляют собой пространственные компоненты принятой акустической волны, и осуществляют обработку первого сигнала в процессе адаптивной фильтрации, включающем формирование сигнала обратной связи, представляющего собой выбранную часть первого сигнала, сложенную со вторым сигналом.
Определение выбранной части первого сигнала предпочтительно осуществляют, по крайней мере частично, посредством использования предварительно определенных параметров окна для выходного сигнала процесса адаптивной фильтрации или, по крайней мере частично, посредством использования частотного фильтра для выходного сигнала процесса адаптивной фильтрации. При определении параметров окна осуществляют выбор зависимого от сдвига фаз временного окна на основе, по крайней мере частично, разницы фаз вертикальной и горизонтальной компонент акустической волны.
Первый и второй сигналы предпочтительно представляют собой сигналы с углом сдвига фаз около 90°. Первый сигнал является горизонтальной компонентой акустической волны, а второй сигнал является вертикальной компонентой акустической волны.
В предпочтительном варианте осуществления изобретения в процессе адаптивной фильтрации из сейсмического сигнала удаляют поверхностные волны.
Предварительно определяемые параметры окна включают параметры фазы и времени.
Предпочтительно, если при осуществлении способа используют адаптивный фильтр, основанный, по крайней мере частично, на начальных коэффициентах фильтрации, и при этом дополнительно осуществляют корректировку начальных коэффициентов фильтрации при использовании адаптивного фильтра в отношении успешно регистрируемых сейсмических волн.
Одна из целей настоящего изобретения заключается в том, чтобы уменьшить или исключить измерительные сети геофонов и подавить поверхностные волны. Предлагаемая в соответствии с настоящим изобретением система не зависит от азимутальных параметров, связанных с конструкцией измерительной сети, содержащей большое число приемников. Настоящее изобретение предоставляет возможность сконструировать систему, использующую как основу конструкции полезный сигнал вместо шума.
Настоящее изобретение обеспечивает эффективное отфильтровывание поверхностных волн без использования сетей сейсмоприемников и продумывания их пространственного расположения, а также понижает зависимость результата от поляризации шума. Устройство для фильтрации сейсмических сигналов в соответствии с настоящим изобретением снабжено детектором, установленным на земле с возможностью приема акустической волны и формирования выходного сигнала, содержащего первый и второй сигналы, представляющие собой пространственные компоненты принятой акустической волны, процессором, соединенным с детектором с возможностью приема первого и второго сигналов и содержащим адаптивный фильтр для обработки первого сигнала в процессе адаптивной фильтрации, включающем формирование сигнала обратной связи, представляющего собой выбранную часть первого сигнала, сложенную со вторым сигналом.
В качестве детектора желательно использовать многоосный детектор, содержащий многоосный акселерометр системы электромеханических микроакселерометров. Предпочтительно также использование частотного фильтра, предназначенного для фильтрации выходного сигнала адаптивного фильтра. В свою очередь частотный фильтр может быть использован в качестве фильтра полосы пропускания частот, например в виде запирающего фильтра.
Устройство может содержать выход сигнала, сформированного посредством процесса удаления поверхностной волны из сейсмического сигнала.
Устройство может также содержать средство хранения начальных коэффициентов фильтрации, используемых, по крайней мере частично, в адаптивном фильтре. При этом начальные коэффициенты фильтрации имеют значения, корректируемые процессором при использовании адаптивного фильтра в отношении успешно регистрируемых сейсмических волн.
Краткое описание чертежей
Новые признаки данного изобретения, так же как и изобретение само по себе, будут наилучшим образом поняты из прилагаемых чертежей и их последующего описания, в которых используется единое обозначение ссылок и в которых:
Фиг.1 представляет схематическое изображение системы сбора сейсмических данных в одном из исполнений в соответствии с настоящим изобретением.
Фиг.2 представляет многоосный сенсор для использования в системе, показанной на фиг.1.
Фиг.3 представляет исполнение сенсорного модуля с фиг.2.
Фиг.4 представляет альтернативное исполнение сенсорного модуля с фиг.3.
Фиг.5 представляет обычный адаптивный фильтр.
Фиг.6 представляет адаптивный фильтр в соответствии с настоящим изобретением.
Фиг.7 представляет блок-схему предлагаемого способа сбора сейсмических данных.
Фиг.8 представляет блок-схему альтернативного варианта осуществления способа сбора сейсмических данных.
Осуществление изобретения
Фиг.1 представляет схематическое изображение системы сбора сейсмических данных в одном из возможных вариантов осуществления настоящего изобретения. Система 100 включает блок 102 записи и первичного контроля (центральный контроллер) для передачи команд на один или более удаленных блоков системы 100, а также для получения с них регистрируемых данных. Как показано, блок 102 записи и первичного контроля представляет собой передвижной блок, известный в среде специалистов в данной области техники как "грузовик" (truck) или "собачья будка" (dog house). Блок записи и первичного контроля может быть на выбор установлен на любом автомобиле или транспортном средстве или может быть установлен полустационарно в заранее определенном месте.
Блок 102 записи и первичного контроля (также именуемый в дальнейшем центральным контроллером 102) соединяется с одним или более приемными блоками - ресиверами 104. Блоки ресиверов могут быть установлены вместе с центральным контроллером 102 или ресиверы 104 могут быть установлены в поле, например, под открытым небом. Ресиверы 104 соединяются через двусторонние линии связи с множеством блоков 10а-10b сбора данных (DAU - data acquisition unit). Блоки 106 (106a-106b) сбора данных включают схему для получения сигнала с измерительной сети 108 сейсмических сенсоров (воспринимающих элементов) 108а. Желательно, чтобы каждый сейсмоприемник имел цифровой выход, отражающий регистрируемую сейсмическую активность. В альтернативном исполнении выходной сигнал с сейсмоприемников является аналоговым, и аналого-цифровое преобразование сигнала осуществляется в блоке DAU.
Далее, система 100 включает сейсмический излучатель (источник сейсмических сигналов) 110, такой как группа из одного или более грузовиков с вибраторами, взрывчатое вещество (например, динамит), пневматические пушки или тому подобное. Излучатель 110 и центральный контроллер 102 соединены через двустороннюю линию связи 122. Линия связи 122 может быть как аналоговой, так и цифровой. Желательно, чтобы передаваемые данные включали контрольные установки, пусковые команды, информацию о местоположении (относительном и/или глобальном), информацию о работоспособности излучателей, марку излучателей и любую другую информацию, полезную для управления разведкой.
В процессе сейсмической съемки излучатель 110 передает сейсмическую энергию на поверхность 114 земли, индуцируя тем самым акустическую волну 116, распространяющуюся в основном вниз. Акустическая волна 116 частично отражается от сейсмических неоднородностей, таких как структурная граница 120, что приводит к появлению отраженной волны 118, распространяющейся в основном в обратную сторону по направлению к поверхности 114.
Отраженная волна 118 регистрируется измерительной сетью сейсмоприемников 108 через отдельные сенсоры 108а. Система может использовать любое число воспринимающих акустическую энергию сенсоров, таких как геофоны, микрофоны, детектор на базе акселерометра, гидрофон при работе на заболоченных территориях и т.д. В соответствии с выбором измерительная сеть может состоять из единственного элемента или включать много элементов. Детекторы 108а могут иметь аналоговый или цифровой выход, может применяться и комбинация детекторов в измерительной сети 108, например комбинация типов детекторов и/или типов выходных сигналов.
В предпочтительном исполнении детекторы 108а являются многоосными детекторами, использующими систему электромеханических микроакселерометров (MEMS). Детекторы MEMS детально обсуждаются в патентной заявке US 09/936634, поданной 12 сентября 2001 г., и соответствующей заявке PCT/USOO/07310, поданной 17 марта 2000 г., и включенных в данное описание в качестве ссылки.
Выходные сигналы измерительной сети 108 могут представлять собой как одиночные сигналы от индивидуальных детекторов 108а, так и сигналы от всех детекторов. В одном из исполнений сигналы предварительно обрабатываются с использованием процессора (отдельно не показан), находящегося в каждом блоке сбора данных 10а-10b. Данные затем передаются на ресивер 104. Для простоты ресивер 104 изображен в виде единственного ресивера. Специалисты в данной области техники без труда поймут, что на практике для получения данных с измерительной сети 108 в системе может быть установлено любое число ресиверов. Сигналы могут обрабатываться как в соответствии с настоящим исполнением, так и в любом другом подходящем месте (узле) системы 100, например в ресивере 104 или центральном контроллере 102, используя адаптивный фильтр, предложенный в настоящем изобретении.
В данном исполнении соединения 112 между центральным контроллером 102 и ресивером 104, соединения 124 между ресивером 104 и блоками DAU 106, между каждым блоком DAU и измерительной сетью 108, а также соединения 122 между центральным контроллером 102 и излучателями 110 представляют собой электропроводящие провода 112а, 112b, 122a, 122b, 124a, 124b. Эти провода обозначаются коллективно как соединения 112/124, а в случае необходимости отдельные соединения будут обсуждаться отдельно с использованием соответствующего ссылочного номера. Провода, включающие соединения 112/124, выбираются для адекватной передачи электропитания и данных в соответствии с конкретными требованиями. При отсутствии специальных требований соединения 112/124 могут представлять собой любые провода, пригодные для передачи электропитания и данных. Нет никакой необходимости в том, чтобы на полном пути сигнала одно соединение, например соединение 112а, было бы того же типа, как любое другое, например 112b. Отдельные соединения могут представлять собой любые комбинации проводов, радиочастот, оптического волокна или электромагнитных устройств. Для блоков с автономными источниками питания, например ресиверов с электрическими батареями или излучателей 110 с автономными генераторами, соединения 112/124 могут также включать радиосвязь. В этих альтернативных исполнениях трансиверы (комбинации приемника/передатчика) устанавливаются в блоке, с которым они связаны. Для иллюстрации этого положения соединения 112/124 между центральным контроллером 102 и ресивером 104 показаны пунктиром.
Сейсмоприемник 108а, в соответствии с представленным на фиг.2 вариантом осуществления, включает корпус 205, закрытый первой торцовой крышкой 210 и второй торцовой крышкой 215, а также сенсорный модуль 305. Корпус 205 закрывается торцовыми крышками 210 и 215, которые крепятся винтами 310, 315, 320 и 325.
Уплотнения 330а-d уплотняют соединение корпуса 205 с первой торцовой крышкой 210, а уплотнения 335a-d уплотняют соединение корпуса 205 со второй торцовой крышкой 215. Уплотнения 330а-d и 335a-d могут представлять собой эластомерные кольца, способные сжиматься до достижения необходимого уплотнения.
Желательно, чтобы корпус 205 содержал полость 340 и плоскую поверхность 345. В одном из исполнений корпус 205 представляет собой металлическую трубу, изготовленную из высокопрочного материала для того, чтобы выдерживать нежелательное, но возможное давление.
Сенсорный модуль 305 устанавливается на плоскую поверхность 345 в полости 340 корпуса 205 и соединяется с первой торцовой крышкой 210 с помощью соединения 355, применяющегося для печатных плат.
В представленном на фиг.3. предпочтительном варианте исполнения сенсорный модуль 305 включает сенсорные контейнеры 405а-с и опору 410. Каждый сенсорный контейнер 405а-с имеет свою ось чувствительности 415, 420 и 425. Каждая ось чувствительности 415, 420 и 425 примерно параллельна осям х, у и z соответственно.
В представленном на фиг.4. предпочтительном варианте исполнения сенсорный модуль 305 включает сенсорные контейнеры 405, опору 410 и монолитный корпус (кронштейн) 1902. Сенсорные контейнеры 405 установлены в монолитном корпусе 1902. Сенсорный модуль 305 включает первый сенсорный контейнер 405а, второй сенсорный контейнер 405b и третий сенсорный контейнер 405с. Первый сенсорный контейнер 405а имеет ось чувствительности 415, примерно параллельную оси х, и установлен в монолитном корпусе 1902 для удерживания оси чувствительности 415 параллельно оси х. Второй сенсорный контейнер 405b имеет ось чувствительности 420, примерно параллельную оси у, и установлен в монолитном корпусе 1902 для удерживания оси чувствительности 420 параллельно оси у. Третий сенсорный контейнер 405с имеет ось чувствительности 425, примерно параллельную оси z, и установлен в монолитном корпусе 1902 для выдерживания оси чувствительности 425 параллельно оси z.
Устройство в соответствии с настоящим изобретением содержит адаптивный фильтр для удаления нежелательных шумов, подобных "поверхностным волнам", перед (компьютерной) обработкой сигнала. Фиг.5 представляет обычный известный адаптивный фильтр. Аналогичные фильтры используются для удаления из сейсмических сигналов неоднородно распределенного шума, такого как шум работающего механизма. Пример применения такого фильтра описан в патенте США №4890264 на имя Крюс и др. Обычные фильтры применялись и для удаления поверхностных волн. Пример применения обычного фильтра можно найти в патенте США №4556962 на имя Уидроу. Конструкция Уидроу использует сигнал на входе фильтра n'(k), либо получаемый с отдельного геофона, используемого для записи поверхностных волн, либо представляющий собой сигнал напрямую с акустического излучателя. Запись n'(k) обозначает зашумленный входной сигнал, x(k) - общее обозначение сигнала на входе адаптивного фильтра, а у(k) - сигнала на выходе. Вторым входным сигналом системы является рабочий сигнал d(k), который представляет собой смесь шума n(k) с акустическим сигналом s(k), зарегистрированным обычным акустическим сенсором. Входной сигнал d(k) и выходной сигнал с адаптивного фильтра у(k) затем суммируются для определения погрешности e(k). Погрешность e(k) в режиме обратной связи передается на адаптивный фильтр для корректировки выходного сигнала y(k) и успешной обработки данных. В такой типичной системе поверхностные волны подавляются за счет суммирования шума со средним сигналом сети геофонов. Во всех описанных случаях имеются в виду одноосные геофоны.
Фиг.6 представляет адаптивный фильтр 500 в соответствии с настоящим изобретением. Адаптивный фильтр принимает сигнал с многоосного сейсмоприемника, передающего на фильтр вертикальную и горизонтальную (или радиальную) компоненты сигнала. Поверхностные волны демонстрируют движение частиц твердого тела по эллиптическим траекториям в сторону, противоположную направлению движения волны, что приводит к появлению вертикальной и горизонтальной компонент сигнала трехосного сейсмического приемника, как описано выше и показано на фиг.2-4. В дальнейшем поляризационная природа поверхностных волн позволяет отличить поверхностные волны от других сигналов и/или шумов, в которых наблюдается корреляция между горизонтальной и вертикальной компонентами.
Проблемой является низкая скорость поверхностных волн в сочетании с его высокой дисперсией. Эти факторы, действующие совместно, делают подавление поверхностных волн при обработке полученного сигнала трудным делом. Адаптивная фильтрация в соответствии с настоящим изобретением использует корреляцию между компонентами записываемого сигнала для идентификации и подавления шума, что достигается применением нестационарных фильтров, обрабатывающих сигнал с единственного приемника.
Когда поверхностные волны имеют значительную амплитуду и создают помеху для чтения сейсмических данных, они ясно различимы на коротких по времени записях как высокоамплитудный, низкочастотный, дисперсный шумовой фон. Однако существуют и другие типы волн, такие как отраженные сдвиговые S-волны и волны Лява, времена вступления которых могут почти совпадать с вступлением поверхностных волн и которые могут быть ошибочно приняты за поверхностные волны. Конструкция фильтра, предлагаемого в настоящем изобретении, выделяет поверхностные волны по времени и по частоте перед фильтрацией, так что фильтрации подвергаются только сейсмические данные, зашумленные поверхностными волнами.
Фиг.6 показывает вертикальную компоненту полученных данных как полезный сигнал плюс шум s(k)+n(k), что составляет входной сигнал системы d(k). Радиальная волна также регистрируется, и зашумленный сигнал n'(k) подается на вход адаптивного фильтра системы. В соответствии с настоящим изобретением выходной сигнал фильтра в дальнейшем подвергается очистке с помощью создания "окна", использующего временные параметры для выделения определенной части радиальных данных. Вспомогательный фильтр полосы пропускания частот может быть также использован для обуживания полосы пропускания частот сигнала. Таким образом инвертированный выходной сигнал у(k) становится более точным перед суммированием с d(k). Это снижает выходную погрешность системы e(k), которая затем возвращается на адаптивный фильтр в режиме обратной связи.
Выделение поверхностных волн осуществляется путем демонстрации на дисплее графиков разности между мгновенными значениями фаз вертикальной и горизонтальной компонент сигнала от одного и того же трехосного сейсмоприемника. Результирующая разница дает возможность легко идентифицировать поляризованный шум.
Система 100, описанная выше и показанная на фиг.1-6, может использоваться в предпочтительном варианте способа в соответствии с настоящим изобретением. Акустический излучатель применяется для генерации в земле акустического (звукового) сигнала. Акустическая энергия расходуется также и на генерацию поверхностных волн, которые достигают многоосных сенсоров примерно в то же самое время, в какое приходит и полезный сигнал, отраженный от глубинного объекта, являющегося целью поиска. Многоосный детектор преобразует полученный акустический сигнал в вертикальную и горизонтальную компоненты выходного сигнала. Вертикальная и горизонтальная компоненты с многоосного детектора 108а передаются на процессор, расположенный, например, в блоке сбора данных. Блок сбора данных содержит в запоминающем устройстве инструкции для осуществления алгоритма фильтрации, применяемого к вертикальной и горизонтальной компонентам выходного сигнала. Алгоритм (процесс) фильтрации осуществляется адаптивным фильтром, под который подбираются входные коэффициенты, адаптируемые к изменяющимся входным условиям. Начальные коэффициенты фильтрации могут сохраняться в запоминающем устройстве для использования процессором, осуществляющим адаптивную фильтрацию. Запоминающее устройство и процессор могут не быть расположены вместе. В одном из исполнений процессор и запоминающее устройство размещены в центральном контроллере. В противоположность этому коэффициенты могут храниться в запоминающем устройстве DAU или в любом другом удобном месте. Процессор также может и не располагаться в центральном контроллере. Процессор может располагаться в любом удобном месте, например в блоке DAU. Способ обработки данных описан ниже и представлен на фиг.7 и 8.
На фиг.7 представлена блок-схема процесса обработки данных, применяемого в предпочтительном варианте осуществления настоящего изобретения - способе 600 адаптивной фильтрации для подавления поверхностных волн в выбранных сигналах. В предпочтительном исполнении инструкции данного способа хранятся в памяти системы, описанной выше и представленной на фиг.1. Эти инструкции выполняются процессором системы для удаления поверхностных волн из вертикальной компоненты сигнала, включающей отраженную сейсмическую волну, т.е. полезный сигнал, и поверхностные волны.
Как описано выше, параметры окна определяются 602 с использованием сигналов, полученных с горизонтального и вертикального выходов многоосного сейсмического сенсора. Параметры окна выбираются исходя из отличительных характеристик поверхностных волн. Например, такой характеристикой является разница фаз вертикальной и горизонтальной компонент сигнала приблизительно в 90°.
Горизонтальная компонента выходного сигнала с сенсора передается на адаптивный фильтр (шаг 604), а вспомогательный фильтр полосы пропускания частот применяется на шаге 606 для выбора конкретной ширины полосы частот. Горизонтальная компонента сигнала затем сравнивается на шаге 608 с параметрами окна и те частоты, которые соответствуют критериям параметров окна, суммируются на шаге 610 с вертикальной компонентой сигнала, включающей поверхностные волны. Результирующий выходной сигнал итерационно передается обратно на адаптивный фильтр для дальнейшей фильтрации, так что вертикальная компонента сигнала в финальной форме оказывается в значительной степени свободна от поверхностных волн.
Обратимся теперь к фиг.8, где показан альтернативный вариант способа адаптивной фильтрации 700 в соответствии с настоящим изобретением. Этот способ в значительной степени аналогичен способу, показанному на фиг.7, однако вертикальная и горизонтальная компоненты выходного сигнала сенсора поменяны местами. Этот способ позволяет фильтровать горизонтальную компоненту сигнала, если это необходимо. Параметры окна определяются на шаге 702 с использованием сигналов, полученных с вертикального и горизонтального выходов многоосного сейсмического сенсора. Параметры окна выбираются исходя из отличительных характеристик поверхностных волн. Например, такой характеристикой является разница фаз вертикальной и горизонтальной компонент сигнала приблизительно в 90°.
Вертикальная компонента выходного сигнала сенсора передается на адаптивный фильтр (шаг 704), а вспомогательный фильтр полосы пропускания частот применяется на шаге 706 для выбора конкретной ширины полосы частот. Вертикальная компонента сигнала затем сравнивается на шаге 708 с параметрами окна и те частоты, которые соответствуют критериям параметров окна, суммируются на шаге 710 с горизонтальной компонентой сигнала, включающей поверхностные волны. Результирующий выходной сигнал итерационно передается обратно на адаптивный фильтр для дальнейшей фильтрации, так что горизонтальная компонента сигнала в финальной форме оказывается в значительной степени свободна от поверхностных волн.
Для удобства иллюстрирования и объяснений представленное выше описание относится конкретному исполнению настоящего изобретения. Однако любому специалисту, знакомому с данной областью техники, будет очевидно, что в рамках настоящего изобретения возможно множество различных модификаций и изменений. Следующие пункты формулы изобретения должны интерпретироваться как имеющие цель охватить все такие модификации и изменения.

Claims (26)

1. Способ фильтрации сейсмических сигналов в системе сейсморазведки, отличающийся тем, что принимают посредством детектора акустическую волну, передают с детектора первый и второй сигналы, которые представляют собой пространственные компоненты принятой акустической волны, и осуществляют обработку первого сигнала в процессе адаптивной фильтрации, включающем формирование сигнала обратной связи, представляющего собой выбранную часть первого сигнала, сложенную со вторым сигналом.
2. Способ по п.1, отличающийся тем, что определение выбранной части первого сигнала осуществляют, по крайней мере частично, посредством использования набора предварительно определенных параметров окна для выходного сигнала процесса адаптивной фильтрации.
3. Способ по п.1, отличающийся тем, что определение выбранной части первого сигнала осуществляют, по крайней мере частично, посредством использования частотного фильтра для выходного сигнала процесса адаптивной фильтрации.
4. Способ по п.3, отличающийся тем, что в качестве частотного фильтра используют фильтр полосы пропускания частот.
5. Способ по п.1, отличающийся тем, что первый и второй сигналы представляют собой сигналы с углом сдвига фаз около 90°.
6. Способ по п.1, отличающийся тем, что первый сигнал является горизонтальной компонентой акустической волны, а второй сигнал является вертикальной компонентой акустической волны.
7. Способ по п.1, отличающийся тем, что в процессе адаптивной фильтрации из сейсмического сигнала удаляют поверхностные волны.
8. Способ по п.2, отличающийся тем, что предварительно определяемые параметры окна включают параметры фазы и времени.
9. Способ по п.1, отличающийся тем, что используют адаптивный фильтр, основанный, по крайней мере частично, на начальных коэффициентах фильтрации.
10. Способ по п.9, отличающийся тем, что дополнительно осуществляют корректировку начальных коэффициентов фильтрации при использовании адаптивного фильтра в отношении успешно регистрируемых сейсмических волн.
11. Способ по п.2, отличающийся тем, что при определении параметров окна осуществляют выбор зависимого от сдвига фаз временного окна на основе, по крайней мере частично, разницы фаз вертикальной и горизонтальной компонент акустической волны.
12. Способ по п.1, отличающийся тем, что в качестве детектора используют многоосный детектор.
13. Способ по п.12, отличающийся тем, что для приема акустической волны используют многоосный детектор системы электромеханических микроакселерометров.
14. Устройство для фильтрации сейсмических сигналов, отличающееся тем, что оно снабжено детектором, установленным на земле с возможностью приема акустической волны и формирования выходного сигнала, содержащего первый и второй сигналы, представляющие собой пространственные компоненты принятой акустической волны, процессором, соединенным с детектором с возможностью приема первого и второго сигналов и содержащим адаптивный фильтр для обработки первого сигнала в процессе адаптивной фильтрации, включающем формирование сигнала обратной связи, представляющего собой выбранную часть первого сигнала, сложенную со вторым сигналом.
15. Устройство по п.14, отличающееся тем, что детектор выполнен в виде многоосного детектора.
16. Устройство по п.15, отличающееся тем, что многоосный детектор содержит многоосный акселерометр системы электромеханических микроакселерометров.
17. Устройство по п.14, отличающееся тем, что оно дополнительно содержит частотный фильтр, предназначенный для фильтрации выходного сигнала адаптивного фильтра.
18. Устройство по п.17, отличающееся тем, что частотный фильтр содержит фильтр полосы пропускания частот.
19. Устройство по п.14, отличающееся тем, что выбранная часть первого сигнала определена, по крайней мере частично, посредством набора предварительно определенных параметров окна для выходного сигнала процесса адаптивной фильтрации.
20. Устройство по п.14, отличающееся тем, что первый и второй сигналы представляют собой сигналы с углом сдвига фаз около 90°.
21. Устройство по п.14, отличающееся тем, что первый сигнал является горизонтальной компонентой акустической волны, а второй сигнал является вертикальной компонентой акустической волны.
22. Устройство по п.14, отличающееся тем, что оно содержит выход сигнала, сформированного посредством процесса удаления поверхностной волны из сейсмического сигнала.
23. Устройство по п.19, отличающееся тем, что предварительно определяемые параметры окна включают параметры фазы и времени.
24. Устройство по п.14, отличающееся тем, что оно дополнительно содержит средство хранения начальных коэффициентов фильтрации, используемых, по крайней мере частично, в адаптивном фильтре.
25. Устройство по п.24, отличающееся тем, что начальные коэффициенты фильтрации имеют значения, корректируемые процессором при использовании адаптивного фильтра в отношении успешно регистрируемых сейсмических волн.
26. Устройство по п.19, отличающееся тем, что параметры окна определены посредством выбора зависимого от сдвига фаз временного окна на основе, по крайней мере частично, разницы фаз вертикальной и горизонтальной компонент акустической волны.
RU2004131199/28A 2002-03-20 2003-03-19 Способ фильтрации сейсмических сигналов и устройство для его осуществления RU2308743C2 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US36611602P 2002-03-20 2002-03-20
US60/366,116 2002-03-20

Publications (2)

Publication Number Publication Date
RU2004131199A RU2004131199A (ru) 2005-05-20
RU2308743C2 true RU2308743C2 (ru) 2007-10-20

Family

ID=28454752

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2004131199/28A RU2308743C2 (ru) 2002-03-20 2003-03-19 Способ фильтрации сейсмических сигналов и устройство для его осуществления

Country Status (9)

Country Link
US (1) US6961283B2 (ru)
EP (1) EP1485733B1 (ru)
CN (1) CN1295520C (ru)
AT (1) ATE537469T1 (ru)
AU (1) AU2003218240A1 (ru)
CA (1) CA2479318A1 (ru)
NO (1) NO334060B1 (ru)
RU (1) RU2308743C2 (ru)
WO (1) WO2003081283A2 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656127A (zh) * 2015-02-04 2015-05-27 安徽吉思勘仪器科技有限公司 一种适用于拖曳式履带阵列地震采集的检波器触地头
RU2669268C1 (ru) * 2018-02-09 2018-10-09 Анна Борисовна Шмелева Способ фильтрации сейсмических сигналов в системе сейсморазведки жидкости

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7823689B2 (en) * 2001-07-27 2010-11-02 Baker Hughes Incorporated Closed-loop downhole resonant source
WO2007130551A2 (en) 2006-05-03 2007-11-15 Exxonmobil Upstream Research Company Data adaptive vibratory source acquisition method
CA2654887A1 (en) * 2006-06-09 2007-12-21 Ion Geophysical Corporation Operating state management for seismic data acquisition
US20080021658A1 (en) * 2006-06-10 2008-01-24 Input/Output, Inc. Apparatus and Method for Integrating Survey Parameters Into a Header
US7843765B2 (en) * 2007-08-09 2010-11-30 Westerngeco L.L.C. Attenuating a surface seismic wave
US8077541B2 (en) * 2007-10-19 2011-12-13 Westerngeco L.L.C. Testing a sensor to produce a filter for noise attenuation
WO2009120430A1 (en) * 2008-03-28 2009-10-01 Exxonmobil Upstream Research Company Method for performing constrained polarization filtering
EP2269094A4 (en) * 2008-03-28 2017-06-14 Exxonmobil Upstream Research Company Characterizing spatial variablility of surface waves in seismic processing
US7864628B2 (en) * 2008-07-09 2011-01-04 Ion Geophysical Corporation Flexural wave attenuation
US8321134B2 (en) * 2008-10-31 2012-11-27 Saudi Arabia Oil Company Seismic image filtering machine to generate a filtered seismic image, program products, and related methods
US8199611B2 (en) * 2009-02-05 2012-06-12 Westerngeco L.L.C. Deriving tilt-corrected seismic data in a multi-axis seismic sensor module
US9829595B2 (en) * 2009-02-06 2017-11-28 Westerngeco L.L.C. Particle motion sensor-based streamer positioning system
US8838392B2 (en) * 2009-10-05 2014-09-16 Westerngeco L.L.C. Noise attenuation in passive seismic data
US8712694B2 (en) * 2009-10-05 2014-04-29 Westerngeco L.L.C. Combining seismic data from sensors to attenuate noise
US9110187B2 (en) * 2009-10-05 2015-08-18 Westerngeco L.L.C. Sensor assembly having a seismic sensor and a divergence sensor
US20110085417A1 (en) * 2009-10-12 2011-04-14 Daniel Ronnow String of Sensor Assemblies Having a Seismic Sensor and Pressure Sensor
GB2474303B (en) * 2009-10-12 2012-02-08 Geco Technology Bv Noise filtering of seismic data
US8520469B2 (en) * 2009-10-12 2013-08-27 Westerngeco L.L.C. Sensor assembly having a seismic sensor, pressure sensor, and processor to apply first and second digital filters
US8654606B2 (en) * 2010-07-27 2014-02-18 Westerngeco L.L.C. Obtaining a response based on differencing of outputs of sensors
GB2482538B (en) * 2010-08-06 2013-10-16 Geco Technology Bv Seismic acquisition and filtering
US9091783B2 (en) 2010-11-04 2015-07-28 Westerngeco L.L.C. Computing a calibration term based on combining divergence data and seismic data
US20120250460A1 (en) * 2011-04-04 2012-10-04 Pascal Edme Noise attenuation using rotation data
US11774616B2 (en) 2011-08-29 2023-10-03 Seismic Innovations Method and system for microseismic event location error analysis and display
US9945970B1 (en) 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
WO2013072861A1 (en) * 2011-11-18 2013-05-23 Geco Technology B.V. Coherent noise attenuation
WO2013169937A1 (en) 2012-05-08 2013-11-14 Octave Reservoir Technologies, Inc. Microseismic event localization using both direct-path and head-wave arrivals
WO2014074133A1 (en) * 2012-11-12 2014-05-15 Conocophillips Company Fracture characterization from refraction travel time data
US9594174B2 (en) 2013-02-01 2017-03-14 Westerngeco L.L.C. Computing rotation data using a gradient of translational data
WO2015109175A1 (en) 2014-01-17 2015-07-23 Westerngeco Llc Seismic sensor coupling
CN103869361B (zh) * 2014-03-07 2016-08-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 自适应低频异常振幅压制的方法
CN204556849U (zh) * 2015-04-27 2015-08-12 魏继东 测耦检波器
CN104932017B (zh) * 2015-05-26 2018-05-04 中国石油天然气集团公司 一种异常振幅压制方法及装置
US10338246B1 (en) 2015-08-31 2019-07-02 Seismic Innovations Method and system for microseismic event wavefront estimation
US20180372900A1 (en) * 2016-03-09 2018-12-27 Cgg Services Sas System and method for acquiring seismic data with flotilla of seismic sources
US11662492B2 (en) 2018-10-25 2023-05-30 Saudi Arabian Oil Company Seismic random noise attenuation
CN113050167B (zh) * 2019-12-27 2022-11-04 中国石油天然气集团有限公司 地震数据的噪声压制方法及装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4458341A (en) * 1981-04-23 1984-07-03 Goebel Vaughn S Method and apparatus for seismic exploration
US4723231A (en) * 1985-03-04 1988-02-02 Texas Instruments Incorporated Source tunable array for ground roll reduction
US4757480A (en) * 1986-04-24 1988-07-12 Amoco Corporation Method and apparatus for filtering seismic data
US4890264A (en) * 1988-03-21 1989-12-26 Atlantic Richfield Company Seismic exploration method and apparatus for cancelling non-uniformly distributed noise
US4858199A (en) * 1988-09-06 1989-08-15 Mobile Oil Corporation Method and apparatus for cancelling nonstationary sinusoidal noise from seismic data
GB2273359B (en) * 1992-12-12 1997-01-15 Schlumberger Ltd Method for improving signal to noise ratio
FR2699347B1 (fr) * 1992-12-14 1995-02-10 Commissariat Energie Atomique Procédé et dispositif d'extraction d'un signal utile d'extension spatiale finie à chaque instant et variable avec le temps.
US5365492A (en) * 1993-08-04 1994-11-15 Western Atlas International, Inc. Method for reverberation suppression
US5448531A (en) * 1994-05-05 1995-09-05 Western Atlas International Method for attenuating coherent noise in marine seismic data
FR2737308B1 (fr) * 1995-07-26 1997-09-19 Inst Francais Du Petrole Methode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu
GB2309082B (en) * 1996-01-09 1999-12-01 Geco As Noise filtering method
GB9623945D0 (en) * 1996-11-15 1997-01-08 Geco Prakla Uk Ltd Detection of ground roll cone
GB9806847D0 (en) * 1998-03-30 1998-05-27 Geco Prakla Uk Ltd Ground roll attenuation method
GB2337591B (en) * 1998-05-20 2000-07-12 Geco As Adaptive seismic noise and interference attenuation method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656127A (zh) * 2015-02-04 2015-05-27 安徽吉思勘仪器科技有限公司 一种适用于拖曳式履带阵列地震采集的检波器触地头
RU2669268C1 (ru) * 2018-02-09 2018-10-09 Анна Борисовна Шмелева Способ фильтрации сейсмических сигналов в системе сейсморазведки жидкости

Also Published As

Publication number Publication date
AU2003218240A8 (en) 2003-10-08
WO2003081283A2 (en) 2003-10-02
ATE537469T1 (de) 2011-12-15
CA2479318A1 (en) 2003-10-02
EP1485733A2 (en) 2004-12-15
AU2003218240A1 (en) 2003-10-08
US6961283B2 (en) 2005-11-01
US20030198132A1 (en) 2003-10-23
EP1485733A4 (en) 2008-04-02
CN1295520C (zh) 2007-01-17
NO20044420L (no) 2004-10-19
EP1485733B1 (en) 2011-12-14
NO334060B1 (no) 2013-12-02
RU2004131199A (ru) 2005-05-20
CN1646938A (zh) 2005-07-27
WO2003081283A3 (en) 2004-01-15

Similar Documents

Publication Publication Date Title
RU2308743C2 (ru) Способ фильтрации сейсмических сигналов и устройство для его осуществления
EP0534648B1 (en) Method of marine seismic exploration
AU2009311497B2 (en) Variable timing ZENSEISTM
US5235554A (en) Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations
EP2507654B1 (en) Extraction of discrete records from continuous seismic recordings
US6381544B1 (en) Deterministic cancellation of air-coupled noise produced by surface seimic sources
AU711341B2 (en) A method for correcting dual sensor data for imperfect geophone coupling using production seismic data
GB2397886A (en) Reducing noise in dual sensor seismic data
AU750647B2 (en) Weighted backus filter method of combining dual sensor traces
US6021092A (en) Method for deriving surface consistent reflectivity map from dual sensor seismic data
US4930113A (en) Suppression of air-coupled noise produced by seismic vibrators
US2654874A (en) Seismic exploration system
Levchenko et al. Recording of ocean-bottom seismic noise and of a strong earthquake in the Himalayas by broadband digital OBS installed on the Mid-Atlantic Ridge
CN109983365A (zh) 用于地震传感器响应修正的系统和方法
Alcudia et al. Microphone experiments and applications in exploration seismology
CA2291333A1 (en) Noise reduction on geophone data using microphone records
Leon et al. Microphone and geophone data analysis for noise characterization and seismic signal enhancement

Legal Events

Date Code Title Description
PD4A Correction of name of patent owner
PC4A Invention patent assignment

Effective date: 20101006

MM4A The patent is invalid due to non-payment of fees

Effective date: 20170320