[go: up one dir, main page]

EA008163B1 - Способ подавления влияния воздушной волны от электромагнитных данных прибрежного частотно-регулируемого источника - Google Patents

Способ подавления влияния воздушной волны от электромагнитных данных прибрежного частотно-регулируемого источника Download PDF

Info

Publication number
EA008163B1
EA008163B1 EA200600106A EA200600106A EA008163B1 EA 008163 B1 EA008163 B1 EA 008163B1 EA 200600106 A EA200600106 A EA 200600106A EA 200600106 A EA200600106 A EA 200600106A EA 008163 B1 EA008163 B1 EA 008163B1
Authority
EA
Eurasian Patent Office
Prior art keywords
source
model
receiver
air
data
Prior art date
Application number
EA200600106A
Other languages
English (en)
Other versions
EA200600106A1 (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 EA200600106A1 publication Critical patent/EA200600106A1/ru
Publication of EA008163B1 publication Critical patent/EA008163B1/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric 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/083Controlled source electromagnetic [CSEM] surveying

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Electromagnetism (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

Способ подавления влияния воздушной волны от прибрежного частотно-регулируемого электромагнитного источника поисковых данных. Исследуемая область моделируется (111) с и (113) без верхнего воздушного слоя. Электромагнитное поле, обусловленное источником, вычисляется для наблюдаемых положений приемника для каждой позиции источника для обеих моделей. Различие между рассчитанными полями (116) для двух моделей - эффект воздушной волны, который вычитается (117) из данных поля для соответствующей геометрии источник-приемник.

Description

Область техники, к которой относится изобретение
Это изобретение относится к области геофизической разведки, и в частности к управляемому электромагнитному источнику разведки углеводородов. Более конкретно, изобретение относится к способу подавления влияния воздушной волны, возникающего в данных, получаемых в частотной области.
Уровень техники
Геофизические исследования с прибрежным управляемым электромагнитным источником (С8ЕМ) используют искусственные электрические и магнитные источники для возбуждения электромагнитных полей в земле и развертывания инструментов/приемников в океане, на морском дне и в буровых скважинах для измерения электрических и магнитных полей. На фиг. 1 показана схематическая диаграмма такого наблюдения с электромагнитным источником 11, соединенным кабелем с судном и приемниками 12, расположенными в океане и частично на морском дне 13. Измеряемые поля анализируются для исследования строения придонных земных структур. Эта технология применима не только для изучения океанической тектоники, но также при прибрежных углеводородных и минеральных исследованиях (Ά.Ό. СНасе с1 а1., в Е1сс1гот;щпсБс МеШобз ίπ ЛррНеб Сеорйузюз 2, 931-966 8оае!у о£ Ехр1огайоп Сеорйузю1з!з (1988); 8. Сопз!аЫе апб С.8. Сох, 1. Сеорйз. Кез. 101, 5519-5530 (1996); Ь. МасСгедог е! а1., Сеорйу. 1. 1п!., 146, 217-236 (2001); 8. ЕШпдзглб е! а1., Т1е Ьеабтд Ебде, 972-982 (2002); Т. Е16езто е! а1., Игз! Вгеак 20.3, 144-152 (2002)).
Электромагнитные сигналы, зарегистрированные приемниками, состоят из электромагнитных полей 21, 22 и 23, прошедших через землю 33, морскую воду 32 и воздух 31, соответственно, как показано на фиг. 2. Сигнал 23, прошедший частично через воздух, называют воздушной волной. В прибрежных электромагнитных геофизических исследованиях с управляемым источником обычно используются частоты ниже 1,0 кГц. Известно, что в этом квазистатическом режиме частоты проникновение электромагнитных волн в среду изменяется обратно пропорционально частоте волны и проводимости среды. Этот результат следует из теории явления скин-эффекта (ГЛ. 8!гайоп, Е1ес1готадпейс Тйеогу, р. 504, МасСгатеНй1 (1941)). Так как проводимость морской воды намного выше, чем воздуха и земли, электромагнитные сигналы в морской воде затухают намного быстрее, чем в воздухе и земле. Так, при расстоянии между источником и приемником приблизительно вдвое большем, чем глубина морской воды, зарегистрированные электромагнитные поля, главным образом, распространяются через воздух и землю. Однако только сигналы, проходящие через землю, несут информацию о строении придонных земных структур. Для глубоких придонных объектов 34 электромагнитные поля должны быть генерированы на низких частотах для гарантирования, что переданные электромагнитные сигналы 25 проникают на целевую глубину. К сожалению, для мелкой глубины воды относительно целевой глубины и на низких частотах сигнал воздушной волны может быть доминирующим в приемниках 12 при большом расстоянии до источника 11 так, что целевой сигнал становится едва различим. Очевидно, условия являются лучшими для С8ЕМ разведки, когда сигнал 25 доминирует над объединенными эффектами сигналов 21, 22 и 23.
Интерференция воздушной волны является проблемой, когда измерения осуществляются в частотной области, то есть, когда источник продолжает передавать свои сигналы, в то время как данные собираются в приемниках. Самый простой исходный сигнал - синусоидальный сигнал с выделенной частотой. Для эксплуатационной эффективности многократные частоты могут быть переданы в то же самое время в виде сложной формы волны, такой как прямоугольная волна. Дополнением к частотной области С8ЕМ является временная область С8ЕМ. Во временной области С8ЕМ источник включен и затем выключен после того, как передана желательная форма волны (например, импульс, прямоугольник или ступенчатая функция). Воздушная волна не может быть проблемой во временной области С8ЕМ, потому что воздушная волна будет зарегистрирована в более раннее время и отделится от целевых сигналов. Однако преимущества, предлагаемые в частотной области С8ЕМ в более сложном моделировании и программном обеспечении инверсии, более понятных результатах и в более высококачественных данных, делают частотную область С8ЕМ более широко используемой в геофизических исследованиях, чем временную область С8ЕМ. Поскольку квалифицированным специалистам в этой области техники будет понятно, что, несмотря на предыдущие замечания, все С8ЕМ данные фактически получены во временной области, то есть они собраны регистрирующим устройством как более или менее непрерывный поток чисел, с независимой переменной во времени, что отличием частотной области С8ЕМ являются способ, посредством которого проводится эксперимент (непрерывный источник), и способы, используемые для анализа и интерпретации данных, посредством чего данные разделяются на индивидуальные частотные компоненты, например Фурье-анализ.
Эффект воздушной волны может быть легко иллюстрирован с использованием простой одномерной (1Ό) слоистой модели. Как показано на фиг. 3, сверху вниз, модель состоит из пяти слоев: непроводящий воздух 31, морская вода 32 (проводимость=3,0 сименс/м, глубина которой будет изменяться согласно примерам ниже), илистые породы 33 (1,0 сименс/м, 1,0 км толщиной), резистивный слой продуктивного пласта 34 (0,01 сименс/м, 100,0 м толщиной) и геологический фундамент 35 (1,0 сименс/м). Если резистивный слой 34 является целевым и будет удален из этой модели, новая модель может быть определена как фоновая для первоначальной модели. Единичный горизонтальный электрический дипольный источник 11, направленный по оси X (НЕЭ.х). буксируется в направлении оси X на расстоянии 50 м от дна.
- 1 008163
Донный приемник 12 расположен прямо ниже серединой точки исходной линии буксировки (не показанный на фиг. 3).
Фиг. 4А-4С - графики амплитуды х-компоненты электрических полей (Ех) в зависимости от разнесения источника-приемника в х-направлении для этих Ю моделей и для ее фоновой модели. Глубина морской воды - 5,0 км на фиг. 4А, 1,0 км на фиг. 4В и 100 м на фиг. 4С. Фиг. 4Ό-4Ε показывают соответствующую развернутую фазу для тех же самых трех глубин морской воды. Развернутая фаза получена для изменения абсолютного скачка больше, чем в дополнении от π до 2π. Кривые из маленьких кругов представляют данные от 1Ό, модель и сплошные линии - от фоновой модели. Для глубины морской воды 5 км (фиг. 4А и 4Ό) наблюдается незначительное влияние воздушный волны на данные от обеих моделей для всех разнесений источника и приемника, представленных на фигуре. Большое различие между ГО модельными кривыми 41 и 43 и кривыми фоновой модели 42 и 44 указывает, что сигнал от резистивного слоя, утопленного на 1,0 км ниже морского дна, существенен, когда разнесение источника-приемника становится большим чем ~2 или 3 км. (Недостаточное различие между модельными и фоновыми кривыми для маленького интервала источник-приемник обусловлено соответственно низким ослаблением водного участка 22 и донного участка 21. Вклад от этих двух сигналов доминирует над принятым сигналом для приемников с небольшим смещением (разнесением источник-приемник), даже с целевым слоем в модели.) Когда глубина морской воды уменьшается до 1,0 км (фиг. 4В и 4Е), различие между этими двумя кривыми значительно снижается из-за влияния воздушной волны, то есть на участке 23 на фиг. 2, теперь за счет прохождения через воду, количество которой намного меньше, следует уменьшение ослабления нежелательной воздушной волны. Этот эффект увеличивается с увеличением расстояния. При расстояниях больше чем ~ 6 или 7 км влияния воздушной волны доминируют в принимаемом сигнале для фоновой модели. Это может быть видно на кривых фоновой модели 46 и 48 на фиг. 4В и 4Е, в частности, по разрыву в наклоне амплитудной кривой 46 в районе 6 км и постоянству фазы на фазовой кривой 48 свыше ~7 км. Однако никакие такие особенности не появляются в данных кривых 45 и 47 для Ю модели с утопленным резистивным слоем, потому что сигнал от утопленного резистивного слоя все еще сильнее, чем влияние воздушной волны для Ю модели с 1,0-километровой глубиной воды. Этого больше не происходит, когда глубина морской воды 100 м, для чего на фиг. 4С представлены данные модели с и без резистивного утопленного слоя продуктивного пласта, едва различимые по амплитуде для всех расстояний. (Различие между двумя кривыми на фиг. 4Е для расстояния больше чем ~3 км, - это прежде всего эффект бесконечного продолжения Ю модели, используемого для получения сигнала от цели; фиг. 9В показывает, по существу, устранение этого эффекта с более реалистической моделью. Материалы были бы значительно хуже для данных поля с шумом. Это подразумевает, что эффект влияния воздушной волны доминирует над полученным сигналом даже при том, что сигнал от поддонной цели сильный. Результаты этого примера ясно демонстрируют проблему влияния воздушной волны.
Влияние воздушной волны было исследовано СНауе и Сох в их теоретической цифровой модели для изучения прибрежного СЕ8М при исследовании с горизонтальным электрическим дипольным источником (Α.Ώ. Сйауе апб С.8. Сох, 1. Оеорйук. Век. 87, 5327-5338 (1982)). Сйауе и Сох поняли, что эффект глубины морской воды будет значителен при больших разнесениях источник-приемник, низких частотах или в относительно мелкой воде. Они указали на то, что эффект может быть включен в теорию, если глубина воды и положение источника точно определены, но они не раскрыли никакого метода вычисления эффекта или удаления его из С8ЕМ данных.
Е1бекшо и другие в статье ΕιγκΙ Вгеак, процитированной ранее, не только описали особенности влияния воздушной волны на амплитуду и фазу, но также и наблюдаемый диапазон, в котором воздушная волна доминирует в отклике, и потерю информации о сопротивлении морского дна, увеличивающуюся с уменьшением частоты и глубины воды. Эффект влияния воздушной волны, таким образом, может быть минимизирован выбором соответствующих частот передачи и планированием изысканий на глубокой воде, в которой цель находится на относительно небольшой глубине ниже морского дна. Однако тактика работы вне диапазона доминирования воздушной волны с тщательным выбором частоты передачи и перспективы обзора не может использоваться для изысканий, где влияние воздушной волны неизбежно для частот, которые необходимы для возбуждения цели для формирования регистрируемых сигналов.
При теоретическом изучении электромагнитного исследования морского дна, используя вертикальный магнитный диполь (УМО), Соддоп и Могпкоп установили, что в случае слабопроводящего морского дна полные горизонтальные магнитные поля являются результатом энергии, распространяющейся двумя основными путями: через основание (морское дно) и через воздух только выше морской поверхности. Они также вычислили вклад воздуха и сравнили его с фактическим эффектом на границе море/воздух для того, чтобы продемонстрировать, что полный бесфазовый отклик магнитного поля является приблизительно простой суммой вкладов воздушного и донного пути (ЕН. Соддоп апб Н.Е. Моткоп, Оеорйукюк 35, 476-489 (1970)). Однако это вычисление вклада воздушной волны выполнено, прежде всего, для проверки концепции о двух основных путях распространения энергии без увеличения целевого сигнала для удаления влияния воздушной волны на измеряемые данные, как позволяет изобретение. В заключении авторы об использовании своего исследования подводят итог в цитате: На практике часто желательно
- 2 008163 работать с Ό/В выше этого минимума так, чтобы изменениями в глубины моря можно было пренебречь. Ό - глубина воды, и В - разнесение источника-приемника. Таким образом, подобно Е1бе8то и другим, Соддоп и Моглкоп показывают, как избежать условий типа мелкой воды или больших разнесений источника-приемника, которые имеют тенденцию делать воздушные волны серьезным шумовым источником в С8ЕМ данных.
Соответственно, имеется потребность в надежном способе подавления влияния воздушной волны в частотной области С8ЕМ данных, применяющемся, где такой шум не устраним. Настоящее изобретение удовлетворяет эту потребность.
Сущность изобретения
В одном варианте осуществления настоящее изобретение относится к способу подавления влияния воздушной волны от прибрежного частотно-регулируемого электромагнитного источника поисковых данных, который включает следующие этапы: (а) формирование модели области, имеющей верхний воздушный слой, средний морской водный слой и слой земли основания, с моделью, отражающей известную глубину области и известную проводимость воздуха, морской воды и земли; (Ь) использование модели для вычисления электромагнитного поля во всех местоположениях приемника для каждого положения источника; (с) замена воздушного слоя в модели с большим количеством морской воды для формирования модели без воздуха; (б) вычисление полей для такой же геометрии источник-приемник для модели без воздуха; (е) нормализация параметров приемника и источника между двумя вычисленными откликами и поисковыми данными; (1) вычисление влияния воздушной волны путем вычитания безвоздушного поля из соответствующего поля для модели с воздухом; и (д) коррекция данных поля путем вычитания вычисленного эффекта воздушной волны для каждого местоположения приемника при каждом положении источника.
Краткое описание чертежей
Настоящее изобретение и его преимущества будут лучше понятны при обращении к следующему детальному описанию и прилагаемым чертежам, где фиг. 1 - схематичный чертеж С8ЕМ наблюдений;
фиг. 2 - схематическая диаграмма путей распространения электромагнитного сигнала для С8ЕМ наблюдений;
фиг. 3 иллюстрирует одномерную слоистую модель земли в океаническом местоположении, также показывая С8ЕМ положение приемника и источника;
фиг. 4А-С - графики амплитуды х-компоненты электрического поля Ех в зависимости от разнесения источника-приемника для различной глубины воды, рассчитанных для модели фиг. 3;
фиг. 4Ό-Ε - графики развернутой фазы в зависимости разнесения источника-приемника, рассчитанного для модели для тех же самых трех глубин воды;
фиг. 5А и 5В иллюстрируют фундаментальные модели с и без верхнего слоя воздуха;
фиг. 6А-С иллюстрируют модели, используемые в одномерном примере, на 6 А показана полная модель, на 6В показано удаление целевого слоя (продуктивного пласта) (с фоновой на основании воздуха моделью), и на 6С показано удаление целевого и воздушного слоев (фоновая модель без воздуха);
фиг. 7Α-Ό показывают смоделированные в частотной области С8ЕМ результаты для 1-Ό примера, сравнение скорректированной амплитуды (7 А) и фазы (7В) с скорректированной амплитудой (7С) и фазой (7Ό);
фиг. 8А-С иллюстрируют модели, используемые в трехмерном примере, на 8А показана полная модель, на 8В показано удаление целевого слоя продуктивного пласта, и на 8С показано отсутствие целевого слоя или воздушного слоя;
фиг. 9Α-Ό показывают моделирование результатов С8ЕМ для трехмерного примера, сравнение некорректированной амплитуды (9 А) и фазы (9В) с скорректированной амплитудой (9С) и фазой (9Ό);
фиг. 10 А и 10В показывают теоретические результаты для амплитуды Ех (10Α) и развернутой фазы (10В) для трехмерного примера; и фиг. 11 показывает блок-схему основных этапов в способе согласно настоящему изобретению.
Изобретение будет описано в связи с его предпочтительным вариантом осуществления. Однако степень этого последующего детального описания является определенной к специфическому варианту осуществления или частному использованию изобретения, это предназначено только для иллюстрации и не должно рассматриваться как ограничение возможностей изобретения. Напротив, оно предназначено охватить все альтернативы, модификации и эквиваленты, которые могут быть включены в сущность и объем изобретения, как определено в представленной формуле изобретения.
Подробное описание предпочтительного варианта осуществления изобретения
Настоящее изобретение относится к способу подавления эффекта воздушной волны в прибрежных частотных С8ЕМ данных. Это требует последовательного ввода информации, подлежащей измерению, рассчитанной или заранее известной или оцененной:
a. вычисление электромагнитных полей для данной модели и каждой предписанной геометрии источника и приемника;
b. измерение глубин области наблюдений;
- 3 008163
с. измерение профиля проводимости морской воды в области наблюдений;
б. измерение (или получение информации иным способом) проводимости морского дна в области наблюдений;
е. измерение амплитуды и фазы электрического излучаемого сигнала одним или более источниками для каждой выбранной позиции;
ί. измерение электромагнитных сигналов одно- или многокомпонентными приемниками, которые расположены в предварительно выбранных положениях;
д. измерение положения электромагнитных приемников для каждой из выбранных позиций; и
11. измерение положения источников (концы или/и поток проводной линии) для каждой из выбранных позиций.
В предпочтительном варианте осуществления настоящего изобретения положение и ориентация источника и приемника измеряются с использованием таких технических методов, как акустические методы суперкороткой базовой линии (88ВЬ), глобальная система позиционирования (СР8), магнитный компас, инерционная навигация, другими техническими методами, известными в позиционировании и навигации.
Для данных поля с эффектом воздушной волны (ЕИАА) основные этапы настоящего способа согласно изобретению выполняются в последовательности (номера ссылок относятся к блок-схемам фиг. 11): 111 - построение фоновой модели с воздушным слоем (МАА) с той же глубиной, проводимостью морской воды и проводимостью дна в области наблюдений (фиг. 5А); 112 - вычисление электромагнитного отклика (ЭАА) для фоновой модели с воздушным слоем для каждого положения приемника и источника, для заданной геометрии в области наблюдения; 113 - построение фоновой модели без воздушного слоя (ΜΝΑ), с заменой воздушного слоя в верхней части фоновой модели с воздушным слоем (МАА) слоем морской воды с проводимостью поверхностного морского слоя (фиг. 5В); 114 - вычисление электромагнитного отклика (ΌΝΑ) для фоновой модели без воздушного слоя для каждого положения приемника и источника, для заданной геометрии в области наблюдения; 115 - нормализация (в случае необходимости) параметров приемника и источника (таких, как длинна антенны и мощность источника) между двумя вычисленными откликами и полевыми данными; 116 - вычисление эффекта воздушной волны (ААЕ) посредством вычитания отклика воздушной модели без воздушного слоя из отклика фоновой модели с воздушным слоем, ΑАЕ=^АΑ-^NΑ; и 117 - коррекция посредством вычитания воздушной волны из полевых данных для каждой заданной геометрии источника и приемника, Ε^NΑ=Ε^АΑ-ΑАЕ. Результирующими данными будут скорректированные данные с удаленным влиянием воздушной волны.
Квалифицированные специалисты в данной области техники осознают, что этапы 112 и 114 приводят к решению системы уравнений Максвелла для определенного положения источника и приемника и заданных параметров модели и с непрерывным излучением волны, используемой для получения данных при режиме обработки в частотной области. Аналитическое решение уравнений возможно для одномерной модели, возбуждаемой дипольным источником (1.А. Копд, Е1ес1тотадпейс Аауе Тйеогу, 2пб Еб., 312321, 1ойп, Абеу & 8оп§, 1пс (1990)). Аналитические решения не существуют для двухмерной или трехмерной модели, кроме некоторых простых конфигураций типа сфер и цилиндров, где границы проводимости соответствуют поверхности постоянной координаты. Поэтому для многомерных моделей используются численные методы (С.А. Нойтапп, ш Е1ес1тотадпейс Ме1йоб§ ш Аррйеб Сеорйумск, 313-363, 8ос1е1у о! Ехр1ота1юп СеорйубсМх (1988)).
С8ЕМ источник может, без ограничения, иметь один из следующих типов или их комбинаций: (1) горизонтальный электрический диполь; (2) вертикальный электрический диполь; (3) горизонтальный магнитный диполь; и (4) вертикальный магнитный диполь. Квалифицированные специалисты в данной области поймут, что термин диполь не используется здесь в самом строгом смысле, который относится к бесконечно малому источнику. Источник может буксироваться на любой глубине или стационарно находиться в воде или на морском дне. Стационарные источники дают более точные результаты измерений, но в результате при низкой эффективности получения данных. Как правило, источник буксируется на 20-80 м выше морского дна, при скорости 1-4 узлов. Такие более медленные скорости предпочтительны для лучшего позиционирования источника и незначительного искажения формы волны из-за движения источника, например эффекта Допплера. Как в случае источника, приемники могут буксироваться, стационарно устанавливаться на дно или в буровых скважинах. Стационарные приемники будут иметь низкий шум и более точно позиционироваться.
Точность результатов настоящего изобретения будет зависеть от того, насколько хорошо фоновые модели (МАА и М^А) представляют фактические структуры электрической проводимости ниже морского дна в области наблюдения. Проводимости ниже морского дна могут быть получены известными спосабами, включая (а) каротаж; (й) магнитотеллурические измерения; или (с) результат инверсии собранных С8ЕМ данных.
Вышеописанный способ согласно настоящему изобретению эффективно подавляет воздушную волну и увеличивает целевой сигнал. Изобретение было проверено на синтетических данных. Далее приведены результаты для двух тестовых моделей.
Примеры
Во всех последующих примерах принято, что единичной силы передатчик генерирует излучение на
- 4 008163 частоте 0,25 Гц. Это значение выбрано только для иллюстративных целей, и тем, кто специализируется в данной области, понятно, что оно никоим образом не ограничивает настоящее изобретение.
Модель 1Ό, показанная в фиг. 6 А, использовалась для формирования набора данных, используемых как данные поля, которые имеют эффект воздушной волны (номера ссылки на фиг. 5А-5В и 6А-6С определены как для фиг. 3). Эти данные (ЕЭ\УЛ) графически показаны на фиг. 7А (амплитуда в зависимости от разнесения) и 7В (фаза в зависимости от разнесения). Эти данные являются кривыми из маленьких кругов, обозначенные как 71 и 72. Модель, по существу, та же самая, как на фиг. 3, за исключением того, что толщина слоя морской воды 250,0 м для этого примера. Как показано на фиг. 3 (но не показано на фиг. 6), направленный по оси х горизонтальный электрический дипольный источник буксируется на 50,0 м выше морского дна в направлении оси х от -15,0 до 15,0 км. Приемник расположен на морском дне непосредственно ниже середины линии буксировки источника. Соответствие фоновой модели с воздушным слоем (М\УА) и фоновой модели без воздушного слоя (ΜΝΑ) показано на фиг. 6В и 6С соответственно. Отклики этих двух моделей вычислены и представлены графически. Кривые 73 и 74 на фиг. 7 А и 7В являются откликом Э\УА. произведенным от модели М\УА согласно фиг. 6В. Отклик ΌΝΑ от модели ΜΝΑ представлен кривыми 75 и 76 на фиг. 7С и 7Ό. Эффект воздушной волны (А\УЕ) вычислен вычитанием ΌΝΑ из Э\УА. Нормализация осуществлялась установкой параметров приемника и источника таких же, как для сгенерированных данных поля и вычисленных откликов от двух фоновых моделей. Тогда влияние воздушной волны удаляется из данных поля вычитанием А\УЕ из ЕЭХУА, и скорректированные данные (ΕΌΝΑ) показаны на фиг. 7С и 7Ό как кривые 77 и 78. Может быть отмечено, что на фиг. 7А и 7В показаны нескорректированные данные поля по сравнению с фоновыми, очень похожие на результаты на мелкой воде на фиг. 4С и 4Е, где воздушная волна доминирует, и, таким образом, нет никакого разделения по амплитуде между целевыми данными и фоновыми данными. Напротив, скорректированные данные поля согласно фиг. 7 С и 7Ό очень похожи на результаты при глубокой воде, показанные в фиг. 4А и 4Ό, где есть хорошее разделение между целевыми данными и фоновыми данными, потому что влияние воздушной волны незначительно для глубокой воды. Эффект воздушной волны также незначителен на фиг. 7С и 7Ό, потому что он был подавлен способом согласно настоящему изобретению.
Фиг. 8А показывает трехмерную модель, которая использовалась для получения синтетических данных поля (ΕΌνΑ) с эффектом воздушной волны. Эти данные представлены (данные в виде маленьких кружков) на кривых 91 и 92 на фиг. 9А и 9В. Модель состоит из воздуха 31, морской воды 32, резистивной квадратной плиты конечных размеров 81 и однородной земли 33. Верхняя часть модели - проводящее воздушное полупространство, в то время как основание - однородная земля с проводимостью 1,0 сименс/м. Резистивная плита - 4,0 км х 4,0 км х 200 м с проводимостью 0,01 сименс/м и утоплена на 1,0 км ниже морского дна. Слой морской воды толщиной 250,0 м, и его проводимость увеличивается ступенчато от 3,0 сименс/м в основании до 5,0 сименс/м на поверхности моря. Горизонтальный электрический дипольный источник, направленный по оси х 11, буксируется на 50 м выше морского дна вдоль оси х от -15 до 15 км. Замкнутые изогнутые линии, выходящие из источника, представляют электромагнитное поле, сформированное источником. Пять приемников 12 развернуты на морском дне. Один приемник расположен непосредственно выше центра плиты и ниже средней точки линии буксировки источника. Другие четыре приемника разнесены с 1,0-километровом интервалом по линии буксировки (положительный х), однако, только данные от одного приемника, расположенного непосредственно выше края резистивной плиты, используются в этом примере. Соответствующая фоновая модель с воздушным слоем (ΜνΑ) и фоновая модель без воздушного слоя (ΜΝΑ) показаны на фиг. 8В и 8С соответственно. Отклики этих двух моделей (названные как Ω\νΑ и ΌΝΑ соответственно) вычислены. Ω\νΑ показано на фиг 9А и 9В как кривые 93 и 94. ΌΝΑ представлена на фиг. 9С и 9Ό как кривые 95 и 96. Эффект воздушной волны (Α\νΕ) вычислен посредством, вычитания ΌΝΑ из Ω\νΑ. Эффект воздушной волны тогда удаляется из данных поля вычитанием Α\νΕ из ΕΩ\νΑ. Данные с скорректированной амплитудой и фазой (ΕΌΝΑ) показаны как кривая 97 на фиг. 9С и кривая 98 на фиг. 9Ό. Фиг. 9А и 9В показывают, что эффект воздушной волны очень затрудняет прогнозирование наличия цели, даже при том, что цель является очень большой. Однако фиг. 9С и 9Ό показывают, что на скорректированных данных четко разделяется целевой сигнал от фона. Фиг. 10А и 10В показывают теоретические данные без эффекта воздушной волны 101 и 102. Теоретический результат рассчитан для модели, показанной на фиг. 8А с воздушным слоем, замененным морской водой тем же самым способом, которым электрические поля рассчитаны на этапах 112 и 114 на фиг. 11. Сплошные кривые 95 и 96 - данные ΌΝΑ, т.е. рассчитанные для модели согласно фиг. 8С. Сравнение кривой 101 и кривой 95 на фиг. 10А и сравнение кривой 102 и кривой 96 на фиг. 10В указывают, что коррекция является очень эффективной. Этот пример показывает эффективность настоящего изобретения в подавлении эффекта воздушной волны для трехмерной модели.
Предшествующее описание направлено на раскрытие вариантов осуществления настоящего изобретения с иллюстративной целью. Однако для специалистов в данной области техники будет очевидно, что возможны многие модификации и разновидности описанного применения. Например, изобретение обсуждалось с использованием горизонтального электрического диполя как примера электромагнитного источника; однако, изобретение одинаково применимо к любому Ο8ΕΜ источнику, например к вертикальному магнитному диполю. Также изобретение работает в пресной воде так же хорошо, как и в соленой
- 5 008163 воде, и слово море должно интерпретироваться соответственно. Все подобные модификации и варианты входят в объем настоящего изобретения, как определено в представленной формуле изобретения.

Claims (5)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Способ подавления влияния воздушной волны от частотно-регулируемого электромагнитного источника данных наблюдения, собираемых на областях, покрытых морем, содержит следующие этапы:
    (a) формирование первой модели области, имеющей сверху вниз воздушный слой, слой морской воды и слой морского основания, указанная модель отражает известную глубину и проводимость;
    (b) вычисление электромагнитных полей, вызванных источником для всех положений источника и приемника в области наблюдений, заданных для первой модели;
    (c) формирование второй модели области путем замещения воздушного слоя в первой модели морской водой и затем вычисление электромагнитных полей для второй модели при тех же самых положениях источника и приемника;
    (б) нормализация параметров приемника и источника между вычисленными полями и наблюдаемыми данными;
    (е) вычисление эффекта воздушной волны путем вычитания поля, полученного для второй модели, из поля, полученного для первой модели, в каждом положении приемника для каждого положения источника; и (ί) вычитание вычисленного эффекта воздушной волны производится из фактических данных наблюдений в каждом положении приемника для каждого положения источника.
  2. 2. Способ по п.1, в котором электромагнитные поля вычисляются путем решения уравнений Максвелла.
  3. 3. Способ по п.1, в котором параметры приемника и источника относятся к длине антенны и интенсивности источника.
  4. 4. Способ по п.1, в котором источник относится к одному из следующих типов: горизонтальный электрический диполь; горизонтальный магнитный диполь; вертикальный электрический диполь; вертикальный магнитный диполь; комбинация предыдущих типов источников.
  5. 5. Способ по п.1, в котором этапы вычислений используют измерения положения и ориентации для источника и приемников, которые получены с использованием по меньшей мере одного из следующих средств: (а) акустических способов; (Ь) системы глобального позиционирования; (с) магнитного компаса; и (б) инерционной навигации.
EA200600106A 2003-06-26 2004-02-20 Способ подавления влияния воздушной волны от электромагнитных данных прибрежного частотно-регулируемого источника EA008163B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US48268103P 2003-06-26 2003-06-26
PCT/US2004/005024 WO2005010560A1 (en) 2003-06-26 2004-02-20 Method for removing air wave effect from offshore frequency domain controlled-source electromagnetic data

Publications (2)

Publication Number Publication Date
EA200600106A1 EA200600106A1 (ru) 2006-06-30
EA008163B1 true EA008163B1 (ru) 2007-04-27

Family

ID=34102646

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200600106A EA008163B1 (ru) 2003-06-26 2004-02-20 Способ подавления влияния воздушной волны от электромагнитных данных прибрежного частотно-регулируемого источника

Country Status (14)

Country Link
US (1) US7277806B2 (ru)
EP (1) EP1636613A4 (ru)
CN (1) CN100338480C (ru)
AU (1) AU2004260003B2 (ru)
BR (1) BRPI0411967A (ru)
CA (1) CA2526266C (ru)
CO (1) CO5660315A2 (ru)
EA (1) EA008163B1 (ru)
MA (1) MA27972A1 (ru)
MX (1) MXPA05013685A (ru)
MY (1) MY135573A (ru)
NO (1) NO20055918L (ru)
TR (1) TR200505213T1 (ru)
WO (1) WO2005010560A1 (ru)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101107608B (zh) 2005-02-16 2011-01-12 埃克森美孚上游研究公司 通过采样其它频率处的噪声估算一个频率处的噪声的方法和装置
GB2423370B (en) 2005-02-22 2007-05-02 Ohm Ltd Electromagnetic surveying for resistive or conductive bodies
CA2611006C (en) * 2005-06-09 2014-09-30 Exxonmobil Upstream Research Company Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys
US7884612B2 (en) 2005-12-22 2011-02-08 Westerngeco L.L.C. Multi-component field sources for subsea exploration
GB2450641B (en) 2006-01-30 2010-06-09 Exxonmobil Upstream Res Co Method for spatial filtering of electromagnetic survey data
CA2643057C (en) 2006-02-21 2015-05-26 Exxonmobil Upstream Research Company Method for electromagnetic air-wave suppression by active cancellation and shielding
NO327007B1 (no) * 2006-05-24 2009-04-06 Norsk Hydro As Fremgangsmate for elektromagnetisk geofysisk kartlegging av undersjoiske bergartsformasjoner
GB2454835B (en) * 2006-07-13 2011-03-02 Exxonmobil Upstream Res Co Removing air wave noise from electromagnetic survey data
US7860655B2 (en) 2006-07-14 2010-12-28 Westerngeco L.L.C. Electromagnetically detecting thin resistive bodies in shallow water and terrestrial environments
US7657391B2 (en) 2006-07-14 2010-02-02 Westerngeco L.L.C. Electromagnetically detecting thin resistive bodies in shallow water and terrestrial environments
GB0616870D0 (en) * 2006-08-25 2006-10-04 Mtem Ltd Improvements In Marine EM Exploration
US7400977B2 (en) * 2006-10-12 2008-07-15 Schlumberger Technology Corporation Computing values for surveying a subterranean structure based on measurements according to different electromagnetic survey techniques
US7430474B2 (en) * 2006-10-31 2008-09-30 Schlumberger Technology Corporation Removing sea surface-related electromagnetic fields in performing an electromagnetic survey
GB0623279D0 (en) * 2006-11-22 2007-01-03 Statoil Asa Air wave modeling for MCSEM/SBL surveying
GB0702661D0 (en) * 2007-02-12 2007-03-21 Mtem Ltd Improvements in marine EM exploration
US7554881B2 (en) * 2007-07-31 2009-06-30 Westerngeco L.L.C. Determining seawater conductivity for the compensation of electromagnetics-based survey data
US7949470B2 (en) * 2007-11-21 2011-05-24 Westerngeco L.L.C. Processing measurement data in a deep water application
BRPI0909440A2 (pt) * 2008-04-17 2015-12-15 Exxonmobil Upstream Res Co métodos para planejamento de desenvolvimento de reservatório, para suporte de decisão com relação de desenvolvimento de recursos petrolíferos, para planejamento de desenvolvimento de otimização para um reservatório com base em um computador, e para produzir hidrocarbonetos de um reservatório subterrâneo.
EP2291761A4 (en) * 2008-04-18 2013-01-16 Exxonmobil Upstream Res Co MARKOV DECISION-MAKING DECISION-MAKING ASSISTANCE TOOL FOR TANK DEVELOPMENT PLANNING
BRPI0910333A2 (pt) * 2008-04-21 2015-10-06 Exxonmobil Upstream Res Co métodos para planejamento de desenvolvimento de um reservatório, para suporte à decisão que considera o desenvolvimento de recursos de petróleo, para otimização do planejamento de desenvolvimento com base em computador para um reservatório de hidrocarboneto, e para produzir hidrocarbonetos, e, produto de programa de computador
US20090290447A1 (en) * 2008-05-25 2009-11-26 Erik Vigen Measuring Electromagnetic Source Geometry
US8063642B2 (en) * 2008-06-11 2011-11-22 Mtem Ltd Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources
US7795873B2 (en) 2008-07-15 2010-09-14 Mtem Ltd Method for attenuating air wave response in marine electromagnetic surveying
US8010291B2 (en) 2008-11-08 2011-08-30 Westerngeco L.L.C. Processing measurement data that involves computing a derivative of the measurement data that represents a subterranean structure
US8258791B2 (en) 2009-01-27 2012-09-04 Mtem Ltd. Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources to impart electromagnetic signals into a subsurface formation and thereby determining a formation response to each signal
US8131522B2 (en) * 2009-06-26 2012-03-06 Pgs Geophysical As Method for estimating and removing air wave response in marine electromagnetic surveying
EP3092516B1 (en) * 2014-01-08 2021-07-14 Equinor Energy AS Inversion method
BR112017003401B1 (pt) * 2014-08-29 2022-05-03 Pgs Geophysical As Estimativa conjunta de respostas eletromagnéticas de terra e ruído ambiente
CN104391332B (zh) * 2014-11-14 2017-06-23 吉林大学 浅海双频可控源电磁勘探方法
US10365394B2 (en) * 2015-08-03 2019-07-30 Groundmetrics, Inc. Workflow to use borehole completion as antenna for geophysical applications
CN105204073B (zh) * 2015-09-18 2017-07-25 中南大学 一种张量视电导率测量方法
CN114764151B (zh) * 2021-01-13 2023-06-23 中国石油化工股份有限公司 一种大地电磁分频层析反演方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792761A (en) * 1987-04-06 1988-12-20 King Thomas C Geophysical prospecting with collimated magnetotelluric fields
US6603313B1 (en) * 1999-09-15 2003-08-05 Exxonmobil Upstream Research Company Remote reservoir resistivity mapping
US6765383B1 (en) * 2001-12-26 2004-07-20 Anthony R. Barringer Magnetotelluric geophysical survey system using an airborne survey bird

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2328026B (en) * 1997-06-16 2002-01-30 Ultra Electronics Ltd Electrical polarisation distribution of a vessel
GB2382875B (en) * 2001-12-07 2004-03-03 Univ Southampton Electromagnetic surveying for hydrocarbon reservoirs
GB2385923B (en) 2002-05-24 2004-07-28 Statoil Asa System and method for electromagnetic wavefield resolution
US6842006B2 (en) * 2002-06-27 2005-01-11 Schlumberger Technology Corporation Marine electromagnetic measurement system
GB2390904B (en) 2002-07-16 2004-12-15 Univ Southampton Electromagnetic surveying for hydrocarbon reservoirs
US7023213B2 (en) * 2002-12-10 2006-04-04 Schlumberger Technology Corporation Subsurface conductivity imaging systems and methods

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792761A (en) * 1987-04-06 1988-12-20 King Thomas C Geophysical prospecting with collimated magnetotelluric fields
US6603313B1 (en) * 1999-09-15 2003-08-05 Exxonmobil Upstream Research Company Remote reservoir resistivity mapping
US6765383B1 (en) * 2001-12-26 2004-07-20 Anthony R. Barringer Magnetotelluric geophysical survey system using an airborne survey bird

Also Published As

Publication number Publication date
BRPI0411967A (pt) 2006-08-29
CA2526266A1 (en) 2005-02-03
CN100338480C (zh) 2007-09-19
MA27972A1 (fr) 2006-07-03
EA200600106A1 (ru) 2006-06-30
MY135573A (en) 2008-05-30
US20070061078A1 (en) 2007-03-15
US7277806B2 (en) 2007-10-02
EP1636613A1 (en) 2006-03-22
NO20055918L (no) 2006-01-04
TR200505213T1 (tr) 2006-08-21
EP1636613A4 (en) 2006-12-27
CA2526266C (en) 2012-01-03
CN1809763A (zh) 2006-07-26
AU2004260003A1 (en) 2005-02-03
WO2005010560A1 (en) 2005-02-03
AU2004260003B2 (en) 2009-04-02
MXPA05013685A (es) 2006-02-24
CO5660315A2 (es) 2006-07-31

Similar Documents

Publication Publication Date Title
EA008163B1 (ru) Способ подавления влияния воздушной волны от электромагнитных данных прибрежного частотно-регулируемого источника
US8437961B2 (en) Time lapse analysis with electromagnetic data
CA2654442C (en) Method for acquiring and interpreting seismoelectric and electroseismic data
CA2383931C (en) Remote reservoir resistivity mapping
US7884612B2 (en) Multi-component field sources for subsea exploration
US8538699B2 (en) Rapid inversion of electromagnetic reconnaissance survey data
AU2006302907B2 (en) Method and apparatus for utilizing time division multiple waveform transmitting
CA2657475A1 (en) Geophysical prospecting with reduced air wave noise
USRE40321E1 (en) Remote reservoir resistivity mapping
EP1779147A1 (en) Processing electromagnetic data
WO2005096021A1 (en) Electromagnetic wavefield analysis
AU2010200948A1 (en) Method for determining resistivity anisotropy from earth electromagnetic responses
RU2349935C2 (ru) Способ и устройство для обработки электромагнитных данных
RU2794698C1 (ru) Способ обнаружения электромагнитных геофизических возмущений от движущихся источников
Birsan Remote sensing of 3-D conducting objects in a layered medium using electromagnetic surface waves
RU2277251C2 (ru) Способ для определения характера подземных резервуаров и способ поиска углеводородсодержащих подземных резервуаров

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): AZ KZ TM

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): RU