EA011315B1 - Оценка шума на одной частоте посредством замера шума на других частотах - Google Patents

Оценка шума на одной частоте посредством замера шума на других частотах Download PDF

Info

Publication number
EA011315B1
EA011315B1 EA200701744A EA200701744A EA011315B1 EA 011315 B1 EA011315 B1 EA 011315B1 EA 200701744 A EA200701744 A EA 200701744A EA 200701744 A EA200701744 A EA 200701744A EA 011315 B1 EA011315 B1 EA 011315B1
Authority
EA
Eurasian Patent Office
Prior art keywords
noise
signal
frequency
frequencies
interest
Prior art date
Application number
EA200701744A
Other languages
English (en)
Other versions
EA200701744A1 (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
Application filed by Эксонмобил Апстрим Рисерч Компани filed Critical Эксонмобил Апстрим Рисерч Компани
Publication of EA200701744A1 publication Critical patent/EA200701744A1/ru
Publication of EA011315B1 publication Critical patent/EA011315B1/ru

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/26—Measuring noise figure; Measuring signal-to-noise ratio
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/083—Controlled source electromagnetic [CSEM] surveying

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Noise Elimination (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

Описываются способ, устройство и компьютерная программа для улучшения отношения сигнал-шум сигнала S(t), причем S(t) содержит сигнал и шум. Выполняется измерение S(t) на представляющей интерес частоте. Выполняются измерения шума в S(t) на одной или нескольких шумовых частотах, где, как ожидается, часть сигнала в S(t) является малой. Оценивается шум на представляющей интерес частоте, используя измерения шума на одной или нескольких шумовых частотах. Оцененный шум вычитается из измерения S(t) на представляющей интерес частоте.

Description

Настоящая заявка притязает на приоритет по предварительной заявке на патент США № 60/653427, поданной 16 февраля 2005 г., и предварительной заявке на патент США № 60/654595, поданной 18 февраля 2005 г.
Уровень техники
Электромагнитное зондирование, включающее в себя электромагнитные исследования с управляемым источником (УИЭМ), проводится посредством передачи электромагнитного сигнала, обычно низкочастотной периодической формы волны, в геологическую среду и измерения электромагнитного отклика. Патент США № 6603313, выданный 8гпка, и публикация патента США № 2003/0050759 (публикация РСТ (Договор о патентной кооперации) № ХУО 03/025803) 8гпка и Сагаххопс описывают способы использования УИЭМ-измерений для разведки нефти и газа и оконтуривания известных объектов разведки.
Как изображено на фиг. 1, которая показывает пример оборудования, участвующего при выполнении морской УИЭМ-съемки, один конец антенны 105 прикреплен к буксируемому телу 110, которое опускается до требуемой глубины воды при помощи подводного буксирного троса 115. Буксируемое тело 110 представляет собой нечто большее, чем просто якорная точка для буксирного троса. Оно обеспечивает место, содержащее электрические компоненты, необходимые для генерирования электромагнитной исходной волны, и также может содержать системы связи, системы позиционирования, устройства измерения скорости звука, высотомеры и т.п., которые полезны при съемке. Лебедка (не показана), прикрепленная к надводному судну 120, управляет буксирным тросом 115.
Антенна 105 представляет собой преобразователь электромагнитных полей, используемых для УИЭМ-съемки. Морская УИЭМ-съемка обычно использует горизонтальные электрические диполи (ГЭД), которые могут быть выполнены следующим образом. Два изолированных проводника отходят от двух выходных клемм генератора мощности, способного подавать электрическую мощность с требуемой частотой и формой волны. Другой конец каждого изолированного проводника соединен с электродом. Альтернативно, изоляция может быть снята с конца изолированного проводника, и неизолированный проводник становится электродом. Два электрода поддерживаются на фиксированном расстоянии друг от друга. Ось диполя поддерживается в горизонтальном положении в случае ГЭД. Токовая петля между двумя электродами замыкается при морском применении водой, морским дном и, возможно, воздухом над водой.
Антенна 105 вырабатывает изменяющееся во времени электромагнитное поле 130, которое в примере, показанном на фиг. 1, проходит сквозь морское дно 125 до формации 135. Изменяющееся во времени электромагнитное поле вызывает протекание изменяющегося во времени тока 140 по формации 135. Протекание изменяющегося во времени тока 140 через формацию 135 наводит другое электромагнитное поле 145. Группа датчиков 150, обычно расположенная на морском дне, принимает, обнаруживает и анализирует электромагнитное поле 145 и сохраняет результирующие данные или передает их на поверхность для анализа. Характеристики принятого электромагнитного поля 145 зависят от характеристик переданного электромагнитного поля 130, которое известно в некоторой степени, характеристик формации 135, характеристик других особенностей и формаций геологической среды и шума. Можно определить некоторые из характеристик формации 135 посредством анализа принятого электромагнитного поля 145 в этом контексте.
Обычный датчик 150, изображенный на фиг. 2, включает в себя электронный модуль 205, соединенный с балластом 210. Электронный модуль 205 включает в себя четыре антенны 215, расположенные примерно симметрично по его периферии. Четыре антенны образуют два электрических диполя, как описано выше. Одна или более вертикальных антенн (не показаны) также могут быть включены для обнаружения вертикально ориентированного электромагнитного излучения. Антенны 215 принимают электромагнитное поле 145, и оборудование в электронном модуле 205 обнаруживает, анализирует и регистрирует данные, относящиеся к фазе и амплитуде электромагнитного поля 145. Когда будет зарегистрировано достаточное количество данных или после прошествия некоторого периода времени, электронный модуль 205 освобождается от балласта 210 и всплывает на поверхность, где его подбирают. На поверхности данные восстанавливаются из электронного модуля 205. Данные анализируются.
Съемка состоит из данных, собранных в результате одного или нескольких проходов антенны 105 по площади морского дна, содержащего один или несколько датчиков 150. Обычно каждый из проходов известен как «интервал буксирного каната».
Одним из наиболее существенных источников шума для УИЭМ-применения являются магнитотеллурические помехи 155, показанные на фиг. 1 в виде множества линий со стрелками. Кроме того, шум может происходить от океанических потоков морского дна, инициирующих механические колебания антенн («вибрации»), или из-за дефектов в электронике датчиков.
Раскрытие изобретения
В общих чертах, в одном аспекте изобретение характеризуется способом улучшения отношения сигнал-шум сигнала 8(1), причем 8(1) содержит сигнал и шум. Способ включает в себя (а) выполнение измерения 8(1) на представляющей интерес частоте, (Ь) выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах, где, как ожидается, часть сигнала 8(1) является малой, (с) оценку шума на
- 1 011315 представляющей интерес частоте, используя измерения шума на одной или нескольких шумовых частотах, и (ά) вычитание оцененного шума из измерения 8(ΐ) на представляющей интерес частоте.
Реализации изобретения могут включать в себя одно или несколько из следующего. Способ может дополнительно включать в себя повторение (а), (Ь), (с) и (ά). Оценка шума на представляющей интерес частоте Ν(Τ) может включать в себя минимизирование
Σ |5(Т) - лг(т)|г, где ^Т)=С1П1(Т)+с2П2(Т)+сзПз(Т)+...; и где сь с2, с3 ... представляют собой комплексные коэффициенты; щ, п2, п3, ... представляют собой измерения шума на шумовых частотах; и Т1 и Т2 определяют временной период, когда в 8(1) присутствует малый сигнал или отсутствует сигнал. Выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах может включать в себя выбор шумовых частот, так что они легко будут моделировать шум на представляющей интерес частоте. Выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах может включать в себя: выбор шумовых частот, которые находятся вблизи представляющей интерес частоты. Выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах может включать в себя выполнение измерений 8(1) в момент времени, когда, как ожидается, часть сигнала в 8(1) является малой. Измерения могут включать в себя данные. Данные могут включать в себя смещение, амплитуду и фазу. Оценка шума может включать в себя сортировку данных по бинам, причем каждый бин ассоциируется с соответствующим диапазоном смещений; преобразование данных в каждом бине в частотную область; выбор из каждого бина данных, ассоциированных с представляющей интерес частотой; выбор из каждого бина данных, ассоциированных с одной или несколькими шумовыми частотами; оценку, для каждого бина, шума на представляющей интерес частоте из выбранных данных о шумовых частотах; и вычитание, на побиновой основе, оцененного шума на представляющей интерес частоте из данных, ассоциированных с представляющей интерес частотой. Выполнение измерения 8(1) на представляющей интерес частоте может включать в себя выполнение измерений комплексной амплитуды 8(1) на представляющей интерес частоте. Выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах может включать в себя выполнение измерений комплексной амплитуды 8(1) на одной или нескольких шумовых частотах. Оценка шума может включать в себя оценку комплексной амплитуды шума на представляющей интерес частоте, используя комплексные амплитуды 8(1). полученные на одной или нескольких шумовых частотах. Вычитание оцененного шума может включать в себя вычитание комплексной амплитуды оцененного шума на представляющей интерес частоте из комплексной амплитуды 8(ΐ), полученной на представляющей интерес частоте.
Необходимо отметить, что зависящая от времени величина ΐ в 8(ΐ) относится к времени как мере положения источника относительно положения приемника. Она не является зависящей от времени величиной, которая преобразуется в частоту, когда данные в каждом бине преобразуются в частотную область. Преобразованная Фурье, зависящая от времени величина представляет собой, например, зависящую от времени величину, представленную на горизонтальной оси фиг. 5, тогда как зависящая от времени величина в 8(1) обозначает момент времени в центре бина после спектрального разложения (преобразования в частотную область).
Способ может включать в себя генерирование УИЭМ-сигнала, имеющего большое количество энергии на представляющей интерес частоте и малое количество энергии на множестве частот с малым сигналом. Способ дополнительно может включать в себя передачу УИЭМ-сигнала и прием сигнала 8(ΐ). Способ дополнительно может включать в себя выбор частоты для УИЭМ-сигнала. Генерирование УИЭМ-сигнала может включать в себя генерирование УИЭМ-сигнала, в котором энергия УИЭМ-сигнала концентрируется по временным частотам, пригодным для различения нефтегазоносных коллекторов. Способ дополнительно может включать в себя выбор шумовых частот, совпадающих с поднабором частот с малым сигналом. Способ дополнительно может включать в себя выбор шумовых частот так, чтобы избежать частотных составляющих передаваемого УИЭМ-сигнала. Генерирование УИЭМ-сигнала может включать в себя генерирование прямоугольной волны. Генерирование УИЭМ-сигнала может включать в себя генерирование волны с тремя пиками.
В общих чертах, в другом аспекте изобретение характеризуется компьютерной программой, хранимой на материальном носителе, для улучшения отношения сигнал-шум сигнала 8(ΐ), причем 8(ΐ) содержит сигнал и шум. Программа включает в себя исполняемые инструкции, которые вызывают компьютером (а) выполнение измерения 8(ΐ) на представляющей интерес частоте; (Ь) выполнение измерений шума в 8(ΐ) на множестве шумовых частот, где, как ожидается, часть сигнала в 8(ΐ) является малой; (с) оценку шума на представляющей интерес частоте, используя измерения шума на одной или нескольких шумовых частотах; и (ά) вычитание оцененного шума из измерения 8(ΐ) на представляющей интерес частоте.
В общих чертах в другом аспекте изобретение характеризуется УИЭМ-устройством для оценки шума на одной частоте в сигнале 8(ΐ), который содержит сигнал и шум, посредством взятия замера шума в 8(ΐ) на других частотах. Устройство включает в себя блок оценки шума для определения оценки шума в 8(ΐ) на представляющей интерес частоте, используя измеренный шум в 8(ΐ) на некоторых шумовых час
- 2 011315 тотах, и блок вычитания шума для вычитания оценки шума из 8(1) на представляющей интерес частоте.
Реализации изобретения могут включать в себя одно или несколько из следующего. УИЭМустройство дополнительно может включать в себе одну или несколько антенн; один или несколько аналого-цифровых преобразователей, соединенных с антеннами, причем каждый вырабатывает данные, представляющие 8(1); и преобразователь данных для преобразования данных из временной области в частотную область. УИЭМ-устройство дополнительно может включать в себя регистратор, соединенный с аналого-цифровым преобразователем, причем регистратор регистрирует выходной сигнал аналогоцифрового преобразователя. УИЭМ-устройство дополнительно может включать в себя блок разделения данных по бинам, соединенный с аналого-цифровым преобразователем для сортировки своего выходного сигнала по бинам, основываясь на смещении. УИЭМ-устройство дополнительно может включать в себя приводящие в определенное состояние компоненты, соединенные с антеннами. Блок оценки шума может минимизировать
где М(Т)=С1П1(Т)+с2П2(Т)+сзПз(Т)+...; Οι, с2, с3 ... представляют собой комплексные коэффициенты; щ, п2, п3 ... представляют собой измерения шума на шумовых частотах; и Т1 и Т2 определяют временной период, когда в 8(1) присутствует малый Сигнал или отсутствует Сигнал. Блок оценки шума может оценивать шум, принимаемый одной антенной. Блок оценки шума может оценивать шум, принимаемый двумя или более антеннами.
Краткое описание чертежей
Фиг. 1 иллюстрирует рабочую среду для УИЭМ-съемки.
Фиг. 2 иллюстрирует датчик УИЭМ-съемки.
Фиг. 3 представляет собой блок-схему излучателя для использования при УИЭМ-съемке.
Фиг. 4 представляет собой блок-схему устройства для приема и обработки УИЭМ-данных.
Фиг. 5 иллюстрирует идеальную прямоугольную волну.
Фиг. 6 иллюстрирует идеальную волну с тремя пиками.
Фиг. 7 изображает часть спектра идеальной прямоугольной волны.
Фиг. 8 изображает часть спектра идеальной волны с тремя пиками.
Фиг. 9 иллюстрирует область практической прямоугольной волны высокой мощности.
Фиг. 10 представляет собой блок-схему последовательности операций, описывающую методику повышения отношения сигнал-шум в УИЭМ-данных.
Фиг. 11-15 иллюстрируют пример улучшения отношения сигнал-шум, создаваемого осуществлением на практике методики, описанной в данном документе.
Осуществление изобретения
Чтобы избежать путаницы, слово сигнал, когда оно написано с прописной буквы, относится к сигнальной составляющей сигнала (без прописной буквы), который включает в себя как сигнал, так и шум. УИЭМ-система повышает отношение сигнал-шум УИЭМ-данных посредством передачи УИЭМэлектромагнитного сигнала, который имеет известные промежутки в своем спектре; приема УИЭМэлектромагнитного сигнала, который в линейной бесшумной среде, как ожидается, имел бы такой же спектральный состав, что и переданный УИЭМ-электромагнитный сигнал; использования шума, принятого в известные промежутки, для оценки шума на представляющих интерес частотах, где, как ожидается, обнаруживается сигнал, и вычитания оцененного шума из принятого УИЭМ-электромагнитного сигнала на представляющих интерес частотах.
Пример устройства для генерирования передаваемого УИЭМ-электромагнитного сигнала 130, которое обычно располагается в буксируемом теле 110, изображен на фиг. 3. Оно включает в себя генератор 305 формы волны, который генерирует форму волны, имеющую требуемые характеристики. Генератор 305 формы волны соединен с УИЭМ-излучателем 310, который передает сгенерированную форму волны антенной 105 и создает переданное электромагнитное поле 130.
При обычном УИЭМ-применении форма волны выбирается такой, чтобы концентрировать доступную мощность излучателя по нескольким выбранным временным частотам, которые выбираются так, чтобы наилучшим образом различать нефтегазоносные коллекторы в геологической среде. Ток излучателя принимает форму, которая повторяется во времени, такую как прямоугольная волна, показанная на фиг. 5, или форму волны с тремя пиками, показанную на фиг. 6. Форма волны с тремя пиками, показанная на фиг. 6, относится к типу, описанному в заявке на патент США № 60/572694, озаглавленной Ьодап11пшс 8рес1гит ТтапктШет ХУауеГогт Гог Соп1то11еб-8оигсе Е1ес1готадпе11с 8игуеутд, Ьи и 8гпка, поданной 20 мая 2004 г.
Из теории анализа Фурье известно, что несинусоидальные формы волны, такие как показанные на фиг. 5 и 6, эквивалентны сумме синусоидальных форм волны, причем каждая представляет конкретную временную частоту. После анализа Фурье амплитуда каждой синусоиды представляет относительный вклад ее частоты в несинусоидальную форму волны. Самая низкая такая частота, в основном, соответствует периоду, в течение которого форма волны повторяется. Например, если повторяется симметричная
- 3 011315
8-секундная прямоугольная волна, показанная на фиг. 5, результирующая непрерывная форма волны состоит из частот (2*Ν+1)/8 Гц, где N=0, 1, 2, ..., или 1/8, 3/8, 5/8 Гц и т.д. Эти частоты и амплитуда сигнала на каждой этой частоте изображены на фиг. 7. Следующие табл. 1 и 2 описывают непрерывную прямоугольную волну и ее несколько первых частотных составляющих, где Т представляет собой период прямоугольной волны (например, 8 секунд на фиг. 5):
Таблица 1. Переходы прямоугольной волны
Переход Время
из 1 в -1 Т/2
из -1 в 1 Т
Таблица 2. Несколько первых частотных составляющих симметричной прямоугольной волны
Частота Амплитуда Фаза
1/Т 1,2732 1 '° 1 'о' 1 ί 1 1 1
з/т 0,4244 0,0
5/Т 0,2546 0,0
7/Т 0,1819 0,0
9/Т 0,1415 0,0
11/Т 0,1157 0,0
13/Т 0,0979 0,0
15/Т 0,0849 0,0
17/Т 0, 0749 0,0
19/Т 0,0670 0,0
21/Т 0,0606 0,0
Табл. 3 и 4 ниже описывают непрерывную форму волны с тремя пиками и ее несколько первых частотных составляющих, где Т представляет собой период формы волны с тремя пиками. Эти частоты и амплитуда сигнала на каждой этой частоте изображены на фиг. 8.
Таблица 3. Переходы волны с тремя пиками
Таблица 4. Несколько первых частотных составляющих волны с тремя пиками
Частота Амплитуда Фаза
1/Т 0,6212 0,0
2/Т 0,6010 -90, 0
4/Т 0,6064 -90,0
7/Т 0,1183 180,0
10/Т 0,2966 90,0
14/Т 0,0801 90, 0
16/Т 0,1596 -90,0
20/Т 0,0756 -90,0
Реальные излучатели высокой мощности не вырабатывают идеальные прямоугольные волны, как показано на фиг. 5, или идеальные волны с тремя пиками, как показано на фиг. 6. Они вырабатывают
- 4 011315 более сложные формы волны, такие как показанные на фиг. 9, в которых положительные амплитуды состоят из положительных полуволн выпрямленного переменного тока, и отрицательные амплитуды состоят из отрицательных полуволн выпрямленного переменного тока. Частотный состав таких сложных форм волны не имеет чистого спектра, показанного на фиг. 7 и 8.
Пример устройства для обнаружения и обработки принятого УИЭМ-электромагнитного сигнала 145, изображенного на фиг. 4, включает в себя антенны 215. Антенны 215 соединены с усилителями и блоками 405 приведения в определенное состояние, которые усиливают и приводят в определенное состояние сигнал от антенн. Приведение в определенное состояние может включать в себя фильтрацию, ослабление или задержку части или всего принятого сигнала. Усилители и блоки 405 приведения в определенное состояние соединены с аналого-цифровым (А/Ц) преобразователем 410, который преобразует аналоговый сигнал в цифровое представление. А/Ц преобразователь имеет достаточно широкую полосу частот и динамический диапазон для регистрации сигнала для анализа. В частности, для устройства улучшения отношения сигнал-шум, описанного в данном документе, А/Ц преобразователь имеет достаточную полосу частот для захвата частот, где ожидается малый сигнал или отсутствие сигнала, как описано ниже. Примерный А/Ц преобразователь работает на частоте 31,125 Гц и обеспечивает 24-битовый выход. Другой примерный А/Ц преобразователь работает на частоте 50 Гц.
Как можно видеть из табл. 4 и на фиг. 8, большая часть энергии в форме волны с тремя пиками с 8секундным периодом приходится на дискретные частоты 1/8 Гц, 2/8 Гц и 4/8 Гц. И наоборот, А/Ц преобразователь производит выборку со значительно меньшим временным интервалом, таким как 0,032 с, позволяя А/Ц преобразователю надежно захватывать частоты от 0 Гц (постоянный ток) до частоты отсечки Найквиста, равной 15,625 Гц. В результате, выходные данные А/Ц преобразователя содержат много частот, отличных от тех, которые были переданы. Например, в случае волны с тремя пиками и с 8секундным периодом (Т=8), выход А/Ц преобразователя, производящий выборку с интервалами 0,032 с, содержит не только представляющие интерес частоты, показанные в табл. 4 (т.е. 1/8, 2/8, 4/8, 7/8, 10/8, 14/8 и 20/8 Гц), но также 3/8, 5/8, 6/8 и т.д. Гц, где ожидается обнаружение только энергии шума.
Возвращаясь обратно к фиг. 4, выход аналого-цифрового преобразователя 410 необязательно соединяется с регистратором 415, который регистрирует цифровые данные для последующей обработки. В некоторых конфигурациях регистратор 415 не используется. В некоторых конфигурациях данные анализируются в реальном времени.
Зарегистрированные данные затем подаются на блок 420 разделения данных по бинам. Блок 420 разделения данных по бинам извлекает временной сегмент зарегистрированных данных, соответствующий интервалу буксирного каната, и делит этот сегмент по бинам. Каждый бин ассоциируется с интервалом времени, простирающимся, в основном, от 2 до 128 с. В случае движущегося источника этот интервал времени может соответствовать диапазону смещений, причем этот диапазон, как правило, попадает между 50 и 600 м. Определяется, что смещение представляет собой расстояние со знаком, которое иногда выражается как время и иногда как физическое расстояние, от антенны 105 до датчика 150, который принял данные, разделяемые на бины. Таким образом, каждый интервал буксирного каната ассоциировал с ним набор бинов, каждый из которых содержит временной сегмент данных от этого интервала.
Данные в каждом бине затем преобразуются преобразователем 425 данных из временной области в частотную область. Преобразователь 425 данных также собирает результирующие данные о комплексной амплитуде (т.е. амплитуда и фаза) со всех бинов для одной или нескольких частот, подлежащих использованию в последующем анализе. Например, преобразователь 425 данных может собирать данные об амплитуде и фазе на частоте 1/8 Гц из каждого бина. Результатом будет совокупность данных, такая как показанная на фиг. 11. Как можно видеть, фиг. 11 имеет две диаграммы для данных, собранных на частоте 1/16 Гц. Обе диаграммы имеют смещение, измеряемое в юлианских днях, в качестве горизонтальной оси. На одной диаграмме вертикальная ось представляет собой амплитуду энергии при соответствующем смещении на горизонтальной оси при 1/16 Гц. На другой диаграмме вертикальная ось относится к фазе энергии при соответствующем смещении на горизонтальной оси при 1/16 Гц. Каждая индивидуальная точка на кривых, показанная на двух диаграммах на фиг. 11, соответствует данным в одном бине.
Данные, собранные преобразователем 425 данных, подаются на блок 430 оценки шума. Блок оценки шума оценивает шум на частотах, где ожидается мощность сигнала, используя шум, собранный на частотах, где не ожидается мощность сигнала. Например, используя спектр, изображенный на фиг. 7 в качестве примера, шум на частотах 1/8, 3/8 и 5/8 Гц может оцениваться с использованием шума, измеренного на частотах 1/16, 3/16, 1/4 и 1/2 Гц.
В одной примерной системе шум на частоте сигнала (такой как 1/8 Гц на фиг. 7) моделируется линейной комбинацией данных на шумовых частотах (таких как 1/16, 3/16, 1/4 и 1/2 Гц на фиг. 7). Т.е. моделируемый шум N на любой сигнальной частоте и время Т бина определяются:
(Г)+с2 пг (Т)+с3и3 (Т)+· · · (1) где С1, С2, ... представляют собой комплексные коэффициенты и п1, п2, ... представляют собой зарегистрированные данные на выбранных шумовых частотах. С; определяются посредством минимизирова- 5 011315 ния в смысле наименьших квадратов разности между зарегистрированным сигналом 8(1) и моделированным шумом в течение временного периода, когда источник или является неактивным, или находится на достаточном расстоянии от приемника, чтобы вносить вклад в малый зарегистрированный сигнал. Другими словами, с1 определяются минимизированием:
где {Т1, Т2} представляет собой временной период, когда присутствует малый сигнал или отсутствует сигнал, и квадрат, как понятно, относится к комплексной величине. Временное окно {Т1,Т2} обычно выбирается с использованием некоторой меры того, как хорошо шум N может моделироваться из п1. В одной примерной системе временное окно выбирается тогда, когда сигнал на представляющих интерес частотах и на шумовых частотах, в основном, находится в фазе, и сигнал является слабым.
При выборе шумовых частот п1 необходимо следить за тем, чтобы слабые, но преднамеренно передаваемые частоты по невнимательности не рассматривались шумом. В частности, форма волны с тремя пиками, показанная на фиг. 6, содержит некоторые гармоники, которые, как правило, считаются слабыми, но являются достаточно сильными, чтобы исказить оценку шума, если они будут включены в процесс оценки шума. В частности, 8-секундная форма волны с тремя пиками содержит очень сильные сигналы на частотах 1/8, 2/8 и 4/8 Гц. Она содержит значительный сигнал на частотах 7/8, 10/8, 14/8 и 16,8 Гц и других гармониках. Сильный и значительный сигнал обозначается на фиг. 8 сплошными кружками. Еще другие, более слабые гармоники (3/8, 5/8, 9/8, 12/8, 13/8, ... Гц) могут рассматриваться слишком слабыми, чтобы обеспечить полезный сигнал, но могут быть достаточно сильными, чтобы исказить оценку шума, если они будут включены в вычисления. Такой более слабый и на уровне шума сигнал обозначается на фиг. 8 незаштрихованными кружками. Сведения о форме волны излучателя и его спектре служат в качестве четкого руководства для выбора шумовых частот. Например, при заданном спектре излучателя, показанном на фиг. 8 для 8-секундной формы волны с тремя пиками, блок 430 оценки шума может выбирать частоты 41/16 и 43/16 Гц, где не ожидается сигнал, для оценки шума на частоте 21/8 Гц, где ожидается сигнал.
Эта методика может быть применена независимо к измерениям, зарегистрированным другой антенной 215 на УИЭМ-датчике 150. Альтернативно, методика может быть применена к любой комбинации данных от различных антенн. В частности, она может быть применена к линейной комбинации измерений, предназначенных для получения составляющей электромагнитного поля, коллинеарной или перпендикулярной антенне излучателя. Чтобы наилучшим образом исследовать аппаратурный шум, способ применяется независимо к различным каналам регистрации в УИЭМ-датчике 150. Канал 402 регистрации может включать в себя одну антенну 215, набор усилителей и блоков 405 приведения в определенное состояние, А/Ц преобразователь 410 и регистратор 415. Альтернативно, каждая из антенн 215 может мультиплексироваться на одни и те же усилители/блок 405 приведения в определенное состояние, А/Ц преобразователь 410 и регистратор 415. Электронный модуль 205 может включать в себя один или несколько каналов регистрации.
Эта методика может быть применена более одного раза к одним и тем же данным, как в случае, когда она сначала применяется для подавления магнитотеллурических помех и затем применяется второй раз для подавления шума от вибрации антенны.
Эта методика может быть применена как для сухопутной, так и для морской УИЭМ-съемки.
В одной примерной системе шумовые частоты выбираются так, чтобы они были достаточно близки к представляющей интерес частоте, чтобы эффективно моделировать шум на представляющей интерес частоте. Наилучший выбор шумовых частот изменяется от набора данных к набору данных в зависимости от спектрального состава и других характеристик шума.
Для специалиста в данной области техники очевидно, что модель оценки шума может быть обобщена с линейной моделью, обсужденной выше, до включения других математических операторов. В частности, свертка или фильтр могут быть применены к значениям данных в некотором диапазоне времени для оценки шума в один момент времени. Альтернативно, метод наименьших квадратов может быть обобщен на другие математические методики оптимизации, временное окно, используемое для разработки коэффициентов, может быть обобщено на включение двух или более временных окон, или коэффициенты с1 могут быть обобщены на включение линейного тренда времени, как в с1+41*Т.
Возвращаясь к фиг. 4, если шум оценивается блоком 430 оценки шума, оценки шума подаются на блок 435 вычитания шума, который вычитает оцененный шум из сигнала для получения оценки сигнала.
В примере использования, изображенном на фиг. 10, УИЭМ-система получает УИЭМ-данные (блок 1005). Как описано выше, получение УИЭМ-данных включает в себя прием принимаемого УИЭМсигнала, наводимого передаваемым УИЭМ-сигналом, где передаваемый УИЭМ-сигнал имеет известную форму волны с известными количествами энергии на представляющих интерес частотах и с малой энергией или отсутствием энергии на других частотах. В линейной бесшумной среде ожидается, что принятый УИЭМ-сигнал будет иметь такой же частотный состав, что и передаваемый УИЭМ-сигнал. Практически, принятый УИЭМ-сигнал содержит как сигнал, так и шум. Выборка принятого УИЭМ-сигнала
- 6 011315 производится с достаточной степенью детализации (т.е. частотой выборки и динамическим диапазоном) для захвата представляющих интерес частот и частот, где, как ожидается, обнаруживается малый сигнал или отсутствие сигнала.
УИЭМ-система затем выбирает по меньшей мере одну сигнальную частоту (блок 1010). Сигнальные частоты выбираются из числа представляющих интерес частот в переданной форме волны. В большинстве случаев выбирается меньшее количество, чем все представляющие интерес частоты.
УИЭМ-система затем идентифицирует одну или несколько шумовых частот (блок 1015). Шумовые частоты выбираются из числа частот, на которых в переданной форме волны было передано мало энергии или отсутствие энергии. В большинстве случаев выбирается меньшее количество, чем все частоты с малой энергией или отсутствием энергии в переданной форме волны.
УИЭМ-система затем измеряет энергию на выбранной шумовой частоте или частотах со смещением, где шум на представляющих интерес частотах легко может быть оценен из шума на шумовых частотах, например, там, где фаза шума на шумовых частотах, вероятно, близка к фазе шума на представляющих интерес частотах. УИЭМ-система использует эти измерения для оценки энергии шума на выбранных сигнальных частотах (блок 1020). Например, может использоваться описанная выше методика оценки.
УИЭМ-система затем вычитает оцененный шум из измеренного сигнала на представляющих интерес частотах (блок 1025), чтобы прийти к оценке сигнала на представляющих интерес частотах. В одной примерной системе вычитание выполняется на побиновой основе.
Пример, иллюстрирующий принцип действия методики, описанной выше, изображен на фиг. 11-15. В примере был передан практический сигнал, имеющий форму волны с тремя пиками, показанную на фиг. 6, с Т=8. Исходные УИЭМ-данные показаны на фиг. 11. Как упомянуто выше, фиг. 11 включает в себя две диаграммы: диаграмму, отражающую амплитуду УИЭМ-данных (фиг. 11 А), и диаграмму, отражающую фазу УИЭМ-данных (фиг. 11В). Горизонтальная ось обеих диаграмм представляет собой смещение, измеряемое во времени, и, конкретно, в юлианском представлении даты. Юлианское представление даты находится в диапазоне от примерно 184,4 до примерно 185,3, который представляет собой период, когда данные были собраны. При других анализах смещение может измеряться расстоянием.
Вертикальная ось амплитудной диаграммы на фиг. 11А представляет собой амплитуду собранных данных в вольтах на метр, измеренных на частоте 0,0625 Гц, которая представляет собой одну из представляющих интерес частот в переданной форме волны данного примера. Вертикальная ось использует логарифмическую шкалу, простирающуюся от примерно 10-6 до примерно 10-12 В/м. В идеальной бесшумной среде данные будут следовать плавной кривой. Как можно видеть на фиг. 11А, данные содержат изрядное количество шума, особенно ниже 184,8 дней и выше 185,1 дней.
Вертикальная ось фазовой диаграммы на фиг. 11В представляет собой косинус фазы собранных данных. В идеальной бесшумной среде данные будут равномерно распределяться от -1 до 1. Как можно видеть на фиг. 11В, данные содержат изрядное количество шума, особенно ниже 184,8 дней и выше 185,1 дней.
Фиг. 12 изображает разность между фазой сигнала, измеренной на одной из представляющих интерес частот, 0,0625 Гц, и фазой сигнала, измеренной на трех частотах, где ожидается малый сигнал или отсутствие сигнала, 0,03125 Гц (фиг. 12А), 0,09375 Гц (фиг. 12В) и 0,15625 Гц (фиг. 12С), по диапазону смещений.
Был выбран диапазон смещений, в котором разность фазы мала, что означает, что шум на представляющей интерес частоте должен легко моделироваться из шума на частотах, где ожидается малый сигнал или отсутствие сигнала. В примере, показанном на фиг. 11-15, данные в диапазоне 184,5-184,6 были выбраны для оценки шума.
Фиг. 13 изображает УИЭМ-сигнал после подавления шума. Сравнивая фиг. 13 с фиг. 11, можно видеть, что уровень шума понизился как на амплитудной диаграмме (фиг. 13А), так и на фазовой диаграмме (фиг. 13В). Снижение шума иллюстрируется на фиг. 14 и 15, которые объединяют диаграмму сигнала перед подавлением шума и диаграмму сигнала после подавления шума. На фиг. 14 данные ниже 184,65 юлианских дней, где расположена линия 1405, представляют собой исходные данные. Данные выше 184,65 юлианских дней представляют собой данные после того, как была применена методика подавления шума. Снижение шума очевидно из снижения средней амплитуды на фиг. 14А. Фиг. 14В изображает соответствующие данные о фазе.
На фиг. 15 данные свыше 185,15 юлианских дней, где располагается линия 1505, представляют собой исходные данные. Данные ниже 185,15 юлианских дней представляют собой данные после того, как была применена методика подавления шума. Снова, снижение шума очевидно из снижения средней амплитуды на фиг. 15А. Соответствующие данные о фазе показаны на фиг. 15В.
Хотя настоящее изобретение было описано со ссылкой на его примерный вариант осуществления, специалист в данной области техники узнает о различных изменениях в форме, которые могут быть сделаны без отступления от сущности и объема заявленного изобретения, определенного в прилагаемой формуле изобретения. Например, специалист в данной области техники осознает, что могут быть использованы другие методики для моделирования шума, отличные от той, которая показана в уравнениях
- 7 011315 и 2. В качестве другого примера, очевидного для специалиста в данной области техники, преобразователь 105 может быть заменен магнитным диполем, что антенна 215 может быть заменена магнитной антенной, или что как преобразователь 105, так и антенна 215 могут быть заменены магнитными устройствами. Все такие изменения, как полагается, включены в нижеследующую формулу изобретения.

Claims (38)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Способ улучшения отношения сигнал-шум сигнала 8(ΐ), причем 8(1) содержит сигнал и шум, при этом способ содержит:
    (a) выполнение измерения 8(1) на представляющей интерес частоте;
    (b) выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах, где часть сигнала в 8(ΐ), как ожидается, является малой;
    (c) оценку шума на представляющей интерес частоте, используя измерения шума на одной или нескольких шумовых частотах; и (ά) вычитание оцененного шума из измерения 8(1) на представляющей интерес частоте.
  2. 2. Способ по п.1, дополнительно содержащий повторение (а), (Ь), (с) и (ά).
  3. 3. Способ по п.1, в котором оценка шума на представляющей интерес частоте Ν(Τ) содержит минимизирование
    Σ |5(Т) - М(Т)2 , где ^Т)=С1П1(Т)+с2П2(Т)+сзИз(Т)+...; и где Οι, с2, с3 ... представляют собой комплексные коэффициенты; щ, п2, п3 ... представляют собой измерения шума на шумовых частотах; и Т и Т2 определяют временной период, когда в 8(1) присутствует малый сигнал или отсутствует сигнал.
  4. 4. Способ по п.1, в котором выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах содержит выбор шумовых частот, так что они легко моделируют шум на представляющей интерес частоте.
  5. 5. Способ по п.1, в котором выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах содержит выбор шумовых частот, которые находятся вблизи представляющей интерес частоты.
  6. 6. Способ по п.1, в котором выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах содержит выполнение измерений 8(1) в момент времени, когда часть сигнала в 8(ΐ), как ожидается, является малой.
  7. 7. Способ по п.1, в котором измерения содержат данные, данные содержат смещение, амплитуду и фазу, и оценка шума содержит сортировку данных по бинам, причем каждый бин ассоциируется с соответствующим диапазоном смещений;
    преобразование данных в каждом бине в частотную область;
    выбор из каждого бина данных, ассоциированных с представляющей интерес частотой;
    выбор из каждого бина данных, ассоциированных с одной или несколькими шумовыми частотами; оценку для каждого бина шума на представляющей интерес частоте из данных о выбранных шумовых частотах и вычитание на побиновой основе оцененного шума на представляющей интерес частоте из данных, ассоциированных с представляющей интерес частотой.
  8. 8. Способ по п.1, в котором выполнение измерения 8(1) на представляющей интерес частоте содержит выполнение измерений комплексной амплитуды 8(1) на представляющей интерес частоте;
    выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах содержит выполнение измерений комплексной амплитуды 8(1) на множестве шумовых частот;
    оценка шума содержит оценку комплексной амплитуды шума на представляющей интерес частоте, используя комплексные амплитуды 8(ΐ), полученные на одной или нескольких шумовых частотах; и вычитание оцененного шума содержит вычитание комплексной амплитуды оцененного шума на представляющей интерес частоте из комплексной амплитуды 8(ΐ), полученной на представляющей интерес частоте.
  9. 9. Способ по п.1, дополнительно содержащий генерирование сигнала для электромагнитных исследований с управляемым источником (УИЭМсигнала), имеющего большое количество энергии на представляющей интерес частоте и малое количество энергии на множестве частот с малым сигналом;
    передачу УИЭМ-сигнала и прием сигнала 8(ΐ).
  10. 10. Способ по п.9, дополнительно содержащий выбор частоты для УИЭМ-сигнала.
  11. 11. Способ по п.9, в котором генерирование УИЭМ-сигнала содержит генерирование УИЭМсигнала, в котором энергия УИЭМ-сигнала концентрируется по временным частотам, пригодным для
    - 8 011315 различения нефтегазоносных коллекторов.
  12. 12. Способ по п.9, дополнительно содержащий выбор шумовых частот, совпадающих с поднабором частот с малым сигналом.
  13. 13. Способ по п.9, дополнительно содержащий выбор шумовых частот так, чтобы избежать частотных составляющих переданного УИЭМ-сигнала.
  14. 14. Способ по п.9, в котором генерирование УИЭМ-сигнала содержит генерирование прямоугольной волны.
  15. 15. Способ по п.9, в котором генерирование УИЭМ-сигнала содержит генерирование волны с тремя пиками.
  16. 16. Компьютерная программа, хранимая на материальном носителе, для улучшения отношения сигнал-шум сигнала 8(1), причем 8(1) содержит сигнал и шум, при этом программа содержит исполняемые инструкции, которые вызывают выполнение компьютером следующих операций:
    (a) выполнение измерения 8(1) на представляющей интерес частоте;
    (b) выполнение измерений шума в 8(1) на одной или нескольких шумовых частотах, где, как ожидается, часть сигнала в 8(1) является малой;
    (c) оценку шума на представляющей интерес частоте, используя измерения шума на одной или нескольких шумовых частотах; и (ά) вычитание оцененного шума из измерения 8(1) на представляющей интерес частоте.
  17. 17. Компьютерная программа по п.16, причем программа дополнительно содержит исполняемые инструкции, которые вызывают компьютером повторение (а), (Ь), (с) и (ά).
  18. 18. Компьютерная программа по п.16, в которой при оценке шума Ν(Τ) на представляющей интерес частоте компьютер минимизирует
    Σ |5(Т) - О2 ,
    Ге{Т,,Т2) где ВДТ)=С1П1(Т)+С2П2(Т)+СзПз(Т)+···; и где С1, с2, с3 ... представляют собой комплексные коэффициенты; щ, п2, п3 ... представляют собой измерения шума на шумовых частотах; и Τι и Т2 определяют временной период, когда в 8(1) присутствует малый сигнал или отсутствует сигнал.
  19. 19. Компьютерная программа по п.16, в которой при выполнении измерений шума в 8(1) на одной или нескольких шумовых частотах компьютер выбирает шумовые частоты, так что они легко моделируют шум на представляющей интерес частоте.
  20. 20. Компьютерная программа по п.16, в которой при выполнении измерений шума в 8(1) на одной или нескольких шумовых частотах компьютер выбирает шумовые частоты, которые находятся вблизи представляющей интерес частоты.
  21. 21. Компьютерная программа по п.16, в которой при выполнении измерений шума в 8(1) на одной или нескольких шумовых частотах компьютер выполняет измерения 8(1) в момент времени, когда часть сигнала в 8(1), как ожидается, является малой.
  22. 22. Компьютерная программа по п.16, в которой измерения содержат данные, данные содержат смещение, амплитуду и фазу, и при оценке шума компьютер сортирует данные по бинам, причем каждый бин ассоциируется с соответствующим диапазоном смещений;
    преобразует данные в каждом бине в частотную область;
    выбирает из каждого бина данные, ассоциированные с представляющей интерес частотой;
    выбирает из каждого бина данные, ассоциированные с одной или несколькими шумовыми частотами;
    оценивает для каждого бина шум на представляющей интерес частоте из данных о выбранных шумовых частотах и вычитает на побиновой основе оцененный шум на представляющей интерес частоте из данных, ассоциированных с представляющей интерес частотой.
  23. 23. Компьютерная программа по п.16, в которой при выполнении измерения 8(1) на представляющей интерес частоте компьютер выполняет измерения комплексной амплитуды 8(1) на представляющей интерес частоте;
    при выполнении измерений шума в 8(1) на одной или нескольких шумовых частотах компьютер выполняет измерения комплексной амплитуды 8(1) на одной или нескольких шумовых частотах;
    при оценке шума компьютер оценивает комплексную амплитуду шума на представляющей интерес частоте, используя комплексные амплитуды 8(1), полученные на одной или нескольких шумовых частотах; и при вычитании оцененного шума компьютер вычитает комплексную амплитуду оцененного шума на представляющей интерес частоте из комплексной амплитуды 8(1), полученной на представляющей интерес частоте.
  24. 24. Компьютерная программа по п.16, причем программа дополнительно содержит исполняемые
    - 9 011315 инструкции, которые вызывают компьютером генерирование сигнала для электромагнитных исследований с управляемым источником (УИЭМсигнала), имеющего большое количество энергии на представляющей интерес частоте и малое количество энергии на множестве частот с малым сигналом;
    передачу УИЭМ-сигнала и прием сигнала 8(1).
  25. 25. Компьютерная программа по п.24, причем программа дополнительно содержит исполняемые инструкции, которые вызывают компьютером выбор частоты для УИЭМ-сигнала.
  26. 26. Компьютерная программа по п.24, в которой при генерировании УИЭМ-сигнала компьютер генерирует УИЭМ-сигнал, в котором энергия УИЭМ-сигнала концентрируется по временным частотам, пригодным для различения нефтегазоносных коллекторов.
  27. 27. Компьютерная программа по п.24, причем программа дополнительно содержит исполняемые инструкции, которые вызывают компьютером выбор шумовых частот, совпадающих с поднабором частот с малым сигналом.
  28. 28. Компьютерная программа по п.24, причем программа дополнительно содержит исполняемые инструкции, которые вызывают компьютером выбор шумовых частот так, чтобы избежать частотных составляющих переданного УИЭМ-сигнала.
  29. 29. Компьютерная программа по п.24, в которой при генерировании УИЭМ-сигнала компьютер генерирует прямоугольную волну.
  30. 30. Компьютерная программа по п.24, в которой при генерировании УИЭМ-сигнала компьютер генерирует волну с тремя пиками.
  31. 31. УИЭМ-устройство для оценки шума на одной частоте в сигнале 8(1), который содержит сигнал и шум, посредством взятия замера шума в 8(1) на других частотах, причем устройство содержит блок оценки шума для определения оценки шума в 8(1) на представляющей интерес частоте, используя измеренный шум в 8(1) на некоторых шумовых частотах; и блок вычитания шума для вычитания оценки шума из 8(1) на представляющей интерес частоте.
  32. 32. УИЭМ-устройство по п.31, дополнительно содержащее одну или несколько антенн;
    один или несколько аналого-цифровых преобразователей, соединенных с антеннами, причем каждый вырабатывает данные, представляющие 8(1); и преобразователь данных для преобразования данных из временной области в частотную область.
  33. 33. УИЭМ-устройство по п.32, дополнительно содержащее регистратор, соединенный с аналого-цифровым преобразователем, причем регистратор регистрирует выходной сигнал аналого-цифрового преобразователя.
  34. 34. УИЭМ-устройство по п.32, дополнительно содержащее блок разделения данных по бинам, соединенный с аналого-цифровым преобразователем для сортировки его выходного сигнала по бинам, основываясь на смещении.
  35. 35. УИЭМ-устройство по п.32, дополнительно содержащее компоненты приведения в определенное состояние, соединенные с антеннами.
  36. 36. УИЭМ-устройство по п.31, в котором блок оценки шума минимизирует где Ы(Т)=С1П1 (Т)+С2П2(Т)+С3П3(Т)+.
    с1, с2, Сз ... представляют собой комплексные коэффициенты; п1, п2, п3 ... представляют собой изме рения шума на шумовых частотах; и Т1 и Т2 определяют временной период, когда в 8(1) присутствует малый сигнал или отсутствует сигнал.
  37. 37. УИЭМ-устройство по п.31, в котором блок оценки шума оценивает шум, принятый одной антенной.
  38. 38. УИЭМ-устройство по п.31, в котором блок оценки шума оценивает шум, принятый двумя или более антеннами.
EA200701744A 2005-02-16 2006-01-17 Оценка шума на одной частоте посредством замера шума на других частотах EA011315B1 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US65342705P 2005-02-16 2005-02-16
US65459505P 2005-02-18 2005-02-18
PCT/US2006/001555 WO2006088591A2 (en) 2005-02-16 2006-01-17 Estimating noise at one frequency by sampling noise at other frequencies

Publications (2)

Publication Number Publication Date
EA200701744A1 EA200701744A1 (ru) 2007-12-28
EA011315B1 true EA011315B1 (ru) 2009-02-27

Family

ID=36916904

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200701744A EA011315B1 (ru) 2005-02-16 2006-01-17 Оценка шума на одной частоте посредством замера шума на других частотах

Country Status (12)

Country Link
US (2) US7593815B2 (ru)
EP (1) EP1849116A4 (ru)
CN (1) CN101107608B (ru)
AU (1) AU2006214688B2 (ru)
BR (1) BRPI0607222A2 (ru)
CA (1) CA2593030C (ru)
EA (1) EA011315B1 (ru)
MA (1) MA29318B1 (ru)
MX (1) MX2007009855A (ru)
MY (1) MY142291A (ru)
NO (1) NO20073107L (ru)
WO (1) WO2006088591A2 (ru)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2006297766B2 (en) 2005-09-19 2011-04-28 Exxonmobil Upstream Research Company Method for performing controlled source electromagnetic surveying with multiple transmitters
WO2007086993A1 (en) 2006-01-30 2007-08-02 Exxonmobil Upstream Research Company Method for spatial filtering of electromagnetic survey data
WO2008092205A1 (en) * 2007-01-31 2008-08-07 Andrew Duncan Method and apparatus for analysing geological features
WO2009062236A1 (en) * 2007-11-12 2009-05-22 Commonwealth Scientific And Industrial Research Organisation Method and apparatus for detecting marine deposits
US7865322B2 (en) * 2008-04-14 2011-01-04 Dh Technologies Development Pte. Ltd. Relative noise
CA2733989C (en) * 2008-09-24 2017-09-26 Exxonmobil Upstream Research Company Systems and methods for subsurface electromagnetic mapping
US8729901B2 (en) * 2009-07-06 2014-05-20 Merlin Technology, Inc. Measurement device and associated method for use in frequency selection for inground transmission
GB2479195A (en) * 2010-04-01 2011-10-05 Wireless Fibre Systems Ltd Buoy supported underwater radio antenna
US20120179372A1 (en) 2010-07-22 2012-07-12 Alexander Edward Kalish Collecting Control Source Electromagnetic Signals
BR112013011732A2 (pt) 2010-11-12 2016-08-16 Halliburton Energy Services Inc método e sistema para realizar medidas ambientais, e, meio de armazenagem legível por máquina
US9383469B2 (en) 2012-04-30 2016-07-05 Pgs Geophysical As Methods and systems for noise-based streamer depth profile control
US20140191760A1 (en) * 2014-01-17 2014-07-10 Bentsion Zinger Method and apparatus for suppression of the airwave in subsea exploration
RU2593384C2 (ru) * 2014-12-24 2016-08-10 Федеральное государственное бюджетное учреждение науки "Морской гидрофизический институт РАН" Способ дистанционного определения характеристик морской поверхности
CN105629317B (zh) * 2016-04-08 2019-02-05 中国矿业大学(北京) 一种基于站间传递函数的大地电磁噪声压制方法
US10482582B2 (en) 2017-02-24 2019-11-19 Sunnybrook Research Institute Systems and methods for noise reduction in imaging
CN107621656B (zh) * 2017-08-31 2019-07-16 湖南五维地质科技有限公司 瞬变电磁法的抽道时窗方法
WO2020069579A1 (en) * 2018-10-04 2020-04-09 Woodside Energy Technologies Pty Ltd A sensor device
EP3935581A4 (en) 2019-03-04 2022-11-30 Iocurrents, Inc. DATA COMPRESSION AND COMMUNICATION USING MACHINE LEARNING
CN117113909B (zh) * 2023-10-25 2024-01-09 湖南大学 一种以振动抑制为导向的变压器优化设计方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030050759A1 (en) * 2001-09-07 2003-03-13 Exxonmobil Upstream Research Company Method of imaging subsurface formations using a virtual source array
US20050251340A1 (en) * 2004-05-06 2005-11-10 Michael Tompkins Electromagnetic surveying for hydrocarbon reservoirs
US20060197534A1 (en) * 2005-03-07 2006-09-07 Exxonmobil Upstream Research Company Method for identifying resistivity anomalies in electromagnetic survey data

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4339720A (en) * 1980-04-04 1982-07-13 The Anaconda Company Telluric noise cancellation in induced polarization prospecting
US4344150A (en) 1980-07-10 1982-08-10 Newmont Mining Corporation Coherent noise cancelling filter
US4812766A (en) * 1987-11-02 1989-03-14 Atlantic Richfield Company Telluric noise cancellation in induced polarization prospecting
US4996484A (en) * 1988-12-29 1991-02-26 Atlantic Richfield Company Method and apparatus for cancelling powerline noise in geophysical electromagnetic exploration
US5448531A (en) * 1994-05-05 1995-09-05 Western Atlas International Method for attenuating coherent noise in marine seismic data
US6359923B1 (en) * 1997-12-18 2002-03-19 At&T Wireless Services, Inc. Highly bandwidth efficient communications
GB9726928D0 (en) * 1997-12-19 1998-02-18 Geco Prakla Uk Ltd Method of stacking seismic signals
MY131017A (en) 1999-09-15 2007-07-31 Exxonmobil Upstream Res Co Remote reservoir resistivity mapping
GB2396013B (en) * 2002-12-04 2006-03-08 Westerngeco Seismic Holdings Processing seismic data
CN100338480C (zh) * 2003-06-26 2007-09-19 埃克森美孚上游研究公司 用于从近海频域受控源电磁数据中除去空气波效应的方法
WO2005117326A2 (en) 2004-05-20 2005-12-08 Exxonmobil Upstream Research Company Logarithmic spectrum transmitter waveform for controlled-source electromagnetic surveying

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030050759A1 (en) * 2001-09-07 2003-03-13 Exxonmobil Upstream Research Company Method of imaging subsurface formations using a virtual source array
WO2003025803A1 (en) * 2001-09-07 2003-03-27 Exxonmobil Upstream Research Company Method of imaging subsurface formations using a virtual source array
US20050251340A1 (en) * 2004-05-06 2005-11-10 Michael Tompkins Electromagnetic surveying for hydrocarbon reservoirs
US20060197534A1 (en) * 2005-03-07 2006-09-07 Exxonmobil Upstream Research Company Method for identifying resistivity anomalies in electromagnetic survey data

Also Published As

Publication number Publication date
EA200701744A1 (ru) 2007-12-28
MA29318B1 (fr) 2008-03-03
WO2006088591A3 (en) 2007-03-29
US20080059072A1 (en) 2008-03-06
WO2006088591A2 (en) 2006-08-24
MY142291A (en) 2010-11-15
CN101107608B (zh) 2011-01-12
AU2006214688B2 (en) 2010-06-17
BRPI0607222A2 (pt) 2009-08-18
US7933741B2 (en) 2011-04-26
AU2006214688A1 (en) 2006-08-24
CA2593030A1 (en) 2006-08-24
EP1849116A4 (en) 2017-07-26
US20090276189A1 (en) 2009-11-05
CN101107608A (zh) 2008-01-16
MX2007009855A (es) 2007-09-04
NO20073107L (no) 2007-07-10
EP1849116A2 (en) 2007-10-31
US7593815B2 (en) 2009-09-22
CA2593030C (en) 2014-02-25

Similar Documents

Publication Publication Date Title
EA011315B1 (ru) Оценка шума на одной частоте посредством замера шума на других частотах
US7504829B2 (en) Methods and apparatus for subsurface geophysical exploration using joint inversion of steady-state and transient data
RU2428719C2 (ru) Способ картирования коллектора углеводородов и устройство для осуществления этого способа
US7667464B2 (en) Time segmentation of frequencies in controlled source electromagnetic (CSEM) applications
JP6806305B2 (ja) 海底資源探査システム、信号処理方法、電気探査法及び電磁探査法
US20110066379A1 (en) survey system for locating geophysical anomalies
US20100026304A1 (en) Method and Apparatus for Analysing Geological Features
Steele et al. Theory and application of calibration techniques for an NDBC directional wave measurements buoy
CA2531801A1 (en) Geophysical data acquisition system
US11567227B2 (en) Multi-mode dispersion energy imaging device and method for a four-component marine interface wave of an ocean bottom seismometer
WO2007003203A1 (en) Method of marine electric logging of oil and gas fields and arrangement of apparatuses 've-so-tem' therefor
US20090072831A1 (en) Source Monitoring for Electromagnetic Surveying
US10429538B1 (en) Underwater electromagnetic field measurement that factors in ocean dynamics
US10209386B2 (en) Processing methods for time division CSEM data
CN105380645B (zh) 一种肺磁图的检测方法与装置
RU48645U1 (ru) Исследовательский комплекс "vesotem" для морской электроразведки нефтегазовых месторождений
RU2010456C1 (ru) Способ измерения уровня звукового давления гидроакустического источника в естественном водоеме и излучатель опорных сигналов для его осуществления
RU2055402C1 (ru) Способ оценки уровня шумов удаленного судоходства в мелком море
CN113433226A (zh) 一种基于数字编码信号检测岩石物性的方法及设备
Farmer Breaking Waves, Bubbles and Langmuir Circulation: a Measurement Strategy for Studying Bubble Enhanced Air-Sea Gas Transfer
Constable et al. Marine magnetotellurics for sub salt exploration 2. A marine equipment system
Murugan et al. Hardware Implementation of Kalman Least Mean Square-Based Adaptive Algorithm for Denoising Ambient Noises in Shallow Water Region
Constable Natural Resource Exploration Using Marine Controlled-Source Electromagnetic Sounding
Norvill The use of distributed sensor arrays in electrical and electromagnetic imaging
NO323242B1 (no) Fremgangsmate og anordning for hydrokarbonidentifisering

Legal Events

Date Code Title Description
MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): AM AZ BY KZ KG MD TJ TM RU