EA007587B1 - Устройство и способ определения удельного сопротивления с помощью направленной электромагнитной волны - Google Patents

Устройство и способ определения удельного сопротивления с помощью направленной электромагнитной волны Download PDF

Info

Publication number
EA007587B1
EA007587B1 EA200400573A EA200400573A EA007587B1 EA 007587 B1 EA007587 B1 EA 007587B1 EA 200400573 A EA200400573 A EA 200400573A EA 200400573 A EA200400573 A EA 200400573A EA 007587 B1 EA007587 B1 EA 007587B1
Authority
EA
Eurasian Patent Office
Prior art keywords
logging tool
transmitting
resistivity
directional
measurements
Prior art date
Application number
EA200400573A
Other languages
English (en)
Other versions
EA200400573A1 (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 EA200400573A1 publication Critical patent/EA200400573A1/ru
Publication of EA007587B1 publication Critical patent/EA007587B1/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/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G01V3/28Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
    • 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/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/30Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

Предложен новый способ оперативной обработки данных, используемый для выделения сигналов из азимутального изменения направленных измерений, проводимых каротажным инструментом внутри ствола скважины. Соответствующие сигналы границ, анизотропии и разрывов выделяются из отклика пласта посредством аппроксимации азимутального изменения измеренных напряжений несколькими синусоидальными функциями. В результате, также получают ориентацию наслоения. Выделенные направленные сигналы используются для получения расстояний до границ и принятия решений по геоуправлению. Можно использовать два способа, предусматривающие инверсию и построение кроссплотов, в зависимости от характера границы. Частью системы является графический интерфейс пользователя (ГИП), облегчающий гибкое определение задач инверсии, что позволяет улучшить результаты инверсии и визуализировать модель пласта, а также инверсионные замеры.

Description

Изобретение относится в основном к области каротажа скважин. Более конкретно, изобретение относится к усовершенствованным способам, при осуществлении которых приборы, оснащенные антенными системами, представленными в виде поперечного или наклоненного щелевого симметричного вибратора, используются для электромагнитных измерений подземных пластов и для надлежащего расположения скважин относительно геологических границ в коллекторе. Изобретение находит основное применение в области каротажа скважин, а конкретно используется на операции каротажа в процессе бурения.
Предшествующий уровень техники
В области разведки и добычи углеводородов хорошо известны различные способы каротажа скважин. При осуществлении этих способов обычно используют приборы или инструменты, оснащенные источниками, предназначенными для излучения энергии в подземный пласт, который пронизан стволом скважины. В этом описании термины «прибор» и «инструмент» будут употребляться как взаимозаменяемые, чтобы указать, например, электромагнитный прибор (или инструмент), инструмент (или прибор), спускаемый в скважину на тросе, либо инструмент (или прибор) для каротажа в процессе бурения. Излученная энергия взаимодействует с подземным пластом, вследствие чего формируются сигналы, которые затем обнаруживаются и измеряются посредством одного или более датчиков. Путем обработки данных, соответствующих обнаруженным сигналам, получают профиль свойств пласта.
Хорошо известны способы каротажа с использованием индукции распространения электромагнитных волн (ЭМ). Каротажные приборы располагают внутри ствола скважины для измерения удельной электропроводности (или обратной ей величины - удельного сопротивления) пластов грунта, окружающих ствол скважины. В данном описании любое указание удельной электропроводности следует считать относящимся также к обратной ей величине удельному сопротивлению, и наоборот. Типичный электромагнитный инструмент для измерения удельного сопротивления содержит передающую антенну и одну или более (в типичном случае - пару) приемных антенн, расположенных на некотором расстоянии от передающей антенны вдоль оси инструмента (см. фиг. 1) .
Индукционные инструменты измеряют удельное сопротивление (или удельную электропроводность) пласта путем изменения напряжения, индуцируемого в приемной антенне (приемных антеннах) в результате появления магнитного потока, индуцируемого токами, протекающими по излучающей антенне (или передающей антенне). Инструмент каротажа на основе распространения ЭМ волн работает аналогичным образом, но, как правило, на более высоких частотах, чем индукционные инструменты, для сравнимых расстояний между антеннами (примерно 106 Гц для инструментов на основе распространения ЭМ волн по сравнению с примерно 104 Гц для индукционных инструментов). Типичный инструмент на основе распространения ЭМ волн может работать в диапазоне частот от 1 кГц до 2 МГ ц.
Обычные передающие и приемные антенны представляют собой катушки, содержащие один или более витков провода в виде изолированного проводника, намотанного на опору. Эти антенны обычно работают как передатчики и/или приемники. Специалистам в данной области техники должно быть понятно, что одну и ту же антенну можно использовать как передающую в один момент времени и как приемную - в другой. Следует также понять, что описываемые здесь передающие и приемные конфигурации являются взаимозаменяемыми ввиду принципа взаимности, т.е. передающая антенна может быть использована как приемная, и наоборот.
Принцип работы антенн состоит в том, что катушка, проводящая ток (например, передающая катушка) генерирует магнитное поле. Электромагнитная энергия из передающей антенны передается в окружающий пласт, и эта передача индуцирует вихревые токи, протекающие в пласте вокруг передатчика (см. фиг. 2А). Эти вихревые токи, индуцируемые в пласте, которые являются функциями удельного сопротивления пласта, генерируют магнитное поле, которое, в свою очередь, генерирует электрическое напряжение в приемных антеннах. Если используется пара отстоящих друг от друга приемников, то индуцируемые напряжения в двух приемных антеннах будут иметь разные фазы и амплитуды ввиду геометрического разноса и поглощения окружающим пластом. Разность фаз (фазовый сдвиг, Ф) и отношение амплитуд (затухание, А) обоих приемников можно использовать для получения удельного сопротивления пласта. Определяемые фазовый сдвиг (Ф) и затухание (А) зависят не только от разнесения двух приемников и расстояния между передатчиком и приемниками, но также и от частоты ЭМ волн, генерируемых передатчиком.
В обычных приборах индукционного и основанного на распространении ЭМ волн каротажа передающие антенны и приемные антенны установлены своими осями вдоль продольной оси прибора. Таким образом, эти инструменты реализуются с антеннами, представленными в виде продольного щелевого симметричного вибратора (ПродЩСВ). Известным способом в области каротажа скважин является использование приборов, включающих в себя антенны, имеющие наклоненные или поперечные катушки, т.е. такие, в которых ось катушки не параллельна продольной оси инструмента. Таким образом, эти инструменты реализованы с поперечной или наклоненной щелевой симметричной вибраторной (ПЩСВ или НЩСВ) антенной. Специалистам в данной области техники известно, что существуют различные способы осуществления наклона антенны. Каротажные приборы, оснащенные ПЩСВ и НЩСВ антенна
- 1 007587 ми, описаны, например, в патентах США №№ 6163155, 5115198, 4319191, 5508616, 5757191, 5781436, 6044325 и 6147496.
На фиг. 2А представлено упрощенное отображение вихревых токов и ЭМ энергии, протекающих от каротажного прибора, расположенного в части или сегменте ствола скважины, который пронизывает подземный пласт в направлении, перпендикулярном слоям осадочных пород. Однако, это не является точным отображением всех многочисленных сегментов, образующих ствол скважины, в частности, когда ствол скважины получен направленным бурением, что описано ниже. Таким образом, сегменты ствола скважины часто пронизывают слои пластов под углом, отличающимся от 90°, как показано на фиг. 2В. Когда это происходит, говорят, что плоскость пласта имеет некоторое относительное падение. Угол θ относительного падения определяется как угол между осью ствола скважины (осью инструмента), ВА, и нормалью N к плоскости Р интересующего слоя пласта.
Известные в данной области техники способы бурения включают в себя бурение стволов скважин из выбранного географического местоположения на поверхности грунта по выбранной траектории. Эта траектория может проходить в другие выбранные географические местоположения на конкретных глубинах в пределах ствола скважины. Эти способы известны под объединяющим названием способов «направленного бурения». Одним приложением направленного бурения является бурение сильно отклоняющихся (от вертикали) или даже горизонтальных стволов скважин в пределах и вдоль относительно тонких углеводородоносных пластов (называемых «продуктивными зонами») на значительные расстояния. Эти сильно отклоняющиеся стволы скважин предназначены для значительного увеличения отбора углеводородов из продуктивной зоны по сравнению с «обычными» стволами скважин, которые пронизывают продуктивную зону «вертикально» (по существу, перпендикулярно залеганию пласта, как показано на фиг. 2А) .
При бурении сильно отклоняющихся или горизонтальных стволов скважин в пределах продуктивной зоны важно поддерживать траекторию ствола скважины таким образом, чтобы она оставалась в конкретном положении в продуктивной зоне. В данной области техники хорошо известны системы направленного бурения, в которых используются «гидравлические забойные двигатели» и «кривые переходники», а также другие средства для управления поворотом траектории ствола скважины по географическим опорным параметрам, таким, как магнитный север, направление силы земного притяжения (вертикаль) и скорость вращения Земли (относительно инерциального пространства). Однако залегание пластов может быть таким, что продуктивная зона не пролегает вдоль заранее заданной траектории в географических местоположениях, отстоящих от местоположения на поверхности ствола скважины. В типичном случае, оператор буровой скважины использует информацию (например, буровые журналы каротажа в процессе бурения (КвПБ)) , полученную во время бурения ствола скважины, для поддержания траектории ствола скважины в пределах продуктивной зоны, а также для того, чтобы убедиться в том, что бурение действительно идет в пределах продуктивной зоны.
Способы поддержания траектории, известные в данной области техники, описаны, например, в документе ТпЬе с1 а1., Ргсс15с Уе11 Р1асешеп1 иыпд Ко1агу 81еетаЬ1е 8у51еш5 апб БЛУО Меакитешеп!, 8ос1е1у о£ Ре1то1еиш Епдшеетк, Рарег 71396, 8ер1ешЬет 30, 2001. Способ, описанный в этом документе, основан на откликах датчиков удельной электропроводности, определяемой посредством КвПБ. Если, например, удельная электропроводность продуктивной зоны известна до проникновения в нее ствола скважины, и если удельные электропроводности вышележащей и нижележащей зон создают значительный контраст по отношению к продуктивной зоне, то замер удельной электропроводности пласта во время бурения можно использовать в качестве критерия «управления» стволом скважины, чтобы оставаться в пределах продуктивной зоны. Более конкретно, если измеренная удельная электропроводность значительно отклоняется от удельной электропроводности продуктивной зоны, то это указывает, что ствол скважины приближается к границе раздела с вышележащим или нижележащим пластом грунта, или даже прошел эту границу. В качестве примера отметим, что удельная электропроводность нефтенасыщенного песка может быть значительно меньшей, чем удельная электропроводность типичных вышележащего и нижележащего сланцев. Указание, что удельная электропроводность около ствола скважины увеличивается, можно интерпретировать как означающее, что ствол скважины подходит к вышележащему или нижележащему слою пласта (в этом примере - к сланцу). Способ направленного бурения с использованием замера свойств пласта в качестве направляющего параметра для коррекции траектории обычно называют «геоуправлением».
Помимо ЭМ измерений для геоуправления также используют измерения акустических свойств и радиоактивности. Снова пользуясь примером нефтепродуктивной зоны с вышележащим и нижележащим сланцами, отметим, что естественная гамма-радиоактивность в продуктивной зоне обычно значительно ниже, чем естественная активность гамма-лучей сланцевых пластов, лежащих выше и ниже продуктивной зоны. В результате, увеличение измеряемой естественной активности гамма-лучей, получаемой из датчика гамма-лучей для КвПБ, будет свидетельствовать, что ствол скважины отклоняется от центра продуктивной зоны и приближается к поверхности раздела с вышележащим или нижележащим сланцем, или даже проходит ее.
- 2 007587
Если, как и в предыдущих примерах, удельная электропроводность и естественная радиоактивность вышележащего и нижележащего сланцевых пластов совпадают, то вышеописанные способы геоуправления показывают лишь то, что ствол скважины покидает продуктивную зону, но не показывают, выходит ли ствол скважины из продуктивной зоны через верх этой зоны или через низ этой зоны. Это создает затруднения для бурового оператора, который должен корректировать траекторию ствола скважины, поддерживая ее выбранное положение в продуктивной зоне.
Приборы ЭМ индукционного каротажа весьма удобны для приложений, связанных с геоуправлением, потому что характерная для них поперечная (радиальная) глубина исследования за счет проникновения в пласты, окружающие скважину, является относительно большой, особенно - по сравнению с ядерными приборами. Возможность более глубокого радиального исследования позволяет индукционным приборам «видеть» на значительном поперечном (или радиальном) расстоянии от оси ствола скважины. В приложениях, связанных с геоуправлением, эта увеличенная глубина исследования обеспечивает обнаружение приближения к границам слоев пластов на более значительных поперечных расстояниях от ствола скважины, что дает буровому оператору дополнительное время для внесения любых необходимых поправок в траекторию. Обычные приборы на основе распространения ЭМ волн выполнены с возможностью разрешения осевых и поперечных (радиальных) изменений удельной электропроводности пластов, окружающих прибор, но характеристика этих приборов обычно не разрешает проводить азимутальные изменения удельной электропроводности пластов, окружающих прибор. Помимо этого, такие приборы не способны воспринимать анизотропию в вертикальных скважинах.
Две экономически важных области отставания делают устранение этих недостатков более срочной проблемой. Первая область отставания связана с возрастающей потребностью в точном размещении скважин, которое требует проведения направленных измерений для принятия решений по управлению для оптимального размещения ствола скважины в коллекторе. Вторая область отставания связана с продуктивными зонами малого удельного сопротивления в слоистых пластах, где точная идентификация и характеристика запасов углеводородов невозможна без знания анизотропии удельного сопротивления. Во многих патентах, полученных за последние годы, описаны способы и устройства для проведения направленных измерений и получения анизотропии удельного сопротивления. Для приложений, связанных с каротажем во время бурения, в патенте США № 5508616 на имя 8а1о и др., описан инструмент индукционного типа с двумя катушками, наклоненными в разных направлениях, не совпадающих с продольной осью инструмента. Направленность измерения иллюстрируется посредством простого утверждения, что функция чувствительности обеих наклоненных катушек нарастает в направлении к области перекрытия зон чувствительности каждой катушки. Указано, что посредством вращения инструмента можно получить изображение азимутального удельного сопротивления пласта на глубине. Однако, в этом патенте не раскрывается, как можно получить азимутальное удельное сопротивление, и не раскрыты никакие методы обнаружения и/или получения характеристик границ, необходимые для принятия количественного решения по геоуправлению.
Патент США № 6181138 на имя НадЩага и 8опд, развивает идею 8а1о и др., связанную с одиночными катушками фиксированной направленности, воплощая ее в трех собранных в блок ортогональных индукционных катушках в местах нахождения передатчика и приемника. О необходимости вращения инструмента не говорится, поскольку направление фокусирования можно настраивать, придерживаясь произвольной ориентации, посредством комбинации откликов ортогональных катушек. Неясно, имеется ли экранирующая конструкция, которая обеспечит прохождение требуемых ЭМ составляющих без сильных неконтролируемых искажений сигнала в случае приложений, связанных с «измерениями в процессе бурения».
В патенте США № 6297639 на имя С1агк и др. и переуступленном обладателю прав на настоящее изобретение, описаны способы и устройства для проведения направленных измерений с использованием различных экранирующих конструкций для обеспечения избирательного затухания энергии ЭМ волн в случаях осевых, наклоненных и поперечных антенных катушек. В этом патенте, помимо прочих особенностей, описаны измерения посредством обычного направленного индукционного каротажа и каротажа на основе распространения ЭМ волн с помощью наклоненных катушек и соответствующих экранов, наряду с нетривиальным способом проведения компенсации ствола скважины при таких измерениях. Описана комбинация одной осевой передающей катушки и одной наклоненной приемной катушки, а также применение этой комбинации для обнаружения направлений границ пластов путем наблюдения азимутального изменения индуцируемого сигнала при вращении инструмента. Азимутальное изменение связи можно использовать для управления поворотом скважин в процессе бурения. С момента выдачи упомянутого патента выдано еще несколько патентов, связанных с экранирующими конструкциями, включая патент США № 6351127 на имя Во81йа1 и др. и патент США № 6566881 на имя Ошегадю и др., причем оба эти патента переуступлены обладателю прав на настоящее изобретение.
Патент США № 6476609 на имя В1йаг, развивает ранее выданный патент США № 6163155, в котором предусматривается использование анизотропии и описываются как передатчики, так и приемники, возможно, имеющие некоторый угол наклона применительно к геоуправлению. Отклик пластов для устройства на основе индукции или распространения ЭМ волн, использующего катушки, наклоненные
- 3 007587 вверх и/или вниз, описывается посредством разности или отношения сигналов при двух разных ориентациях, но об экранировании не упоминается. Не описано и влияние анизотропии или падения пласта. Также ничего не сказано о том, как использовать эти измерения, чтобы получить точное расстояние до границы пласта. В патенте № 6476609 явно предполагается, что ориентация пластов точно известна при вычислении отклика в направлении вверх и/или вниз. Однако, не предложен никакой способ установления точного направления вверх или вниз перед расчетом сигналов, направленных вверх или вниз.
В публикации № 2003/0085707 заявки на патент США на имя МшегЬо и др. и переуступленной обладателю прав на настоящее изобретение, описаны конфигурации инструментов и способы симметризации, которые упрощают получение отклика направленных измерений до момента, когда этот отклик становится почти независимым от анизотропии или угла падения. Отклики на расстояние до границы пласта при разных падении и анизотропии, по существу, перекрываются, за исключением зоны около границы пласта. Для достижения этого упрощения можно симметризовать измерения в двух случаях - индукционного каротажа с помощью двух катушек (одного передатчика и одного приемника: «ПЕ-ПР») и каротажа на основе распространения ЭМ волн с помощью трех катушек (одной передающей и двух приемных: «ПЕ-ПР-ПР»). Симметризацию проводят между двумя наклоненными парами ПЕ-ПР, расположенными в одном и том же промежутке, но при измененном угле наклона передающих катушек и измененном угле наклона приемных катушек. Рассматриваются только те случаи, когда магнитные моменты передающих и приемных катушек оказываются в одной и той же плоскости. Недостаток такого подхода заключается в том, что не удается обеспечить сигнал, необходимый для геоуправления, на протяжении всего времени скольжения, как в случае размещения скважин, при котором приходится учитывать наличие гидравлического забойного двигателя во время угловой коррекции траектории. Если магнитный момент инструмента оказывается параллельным наслоению во время скольжения, то генерируемый сигнал, направленный вверх или вниз, будет нулевым, независимо от расстояния до границы. Таким образом, оперативный контроль расстояния до границы оказывается невозможным.
В публикации № 2003/0200029 заявки на патент США на имя Отегащс и др. и тоже переуступленной обладателю прав на настоящее изобретение, описаны направленные измерения при каротаже на основе распространения ЭМ волн для определения анизотропии в почти вертикальных скважинах с компенсацией стволов скважин. Для получения показаний свойств анизотропных пластов также используют способы инверсии. В публикации № 2003/0184302 заявки на патент США на имя Отегащс и Ектегаоу и переуступленной обладателю прав на настоящее изобретение, также описаны способы прогнозирования в связи с направленными измерениями.
В публикациях №№ 2004/0046560А1 и 2004/0046561А1 заявок на патент США на имя Й/косю/ и др., описано использование квадрупольных антенн и связи индукционных измерений, проводимых с помощью поперечных симметричных вибраторных и квадрупольных антенн, имеющих аналогичные характеристики направленности, с откликом в плоскости ΧΖ обычной поперечной вибраторной антенны. Практическая реализация таких антенн на металлической переходной муфте и их адекватное экранирование не раскрыты. Кроме того, влияние ствола скважины на такие измерения и взаимодействие и/или связь этого влияния с влиянием границ может отличаться от тех характеристик, которые получены путем измерения в плоскости ΧΖ.
Ни в одном из вышеупомянутых патентов не описано использование подробно охарактеризованных азимутальных откликов измеренного сигнала, или способы влияния на такие отклики. В этих документах также ничего не сказано о том, как применять измерение направления для определения расстояний до границ с целью использования при геоуправлении. Упоминается лишь так называемое измерение в направлении вверх или вниз, при котором получают разность измеряемого сигнала для случаев инструмента, сфокусированного непосредственно по направлению к некоторому слою пласта и от него. Точная информация о падении наслоения и азимуте обычно неизвестна до бурения и часто изменяется во время перспективных проработок ситуаций размещения скважин, при которых необходимо геоуправление. Воспользовавшись предварительно определенным направлением наслоения вверх или вниз, можно получить наименее искаженное измерение, что в худшем случае может привести лишь к неточным решениям по геоуправлению при внезапном изменении азимута наслоения. В принципе, возможно азимутальное регулирование измерений в скважине с помощью накопителей. Этот способ имеет ряд недостатков, включая затруднения при точном выравнивании верхнего и нижнего накопителей в соответствии с ориентацией слоев пласта и невозможность использовать данные, которые не находятся в верхнем и нижнем накопителях (т.е. потерю таких данных). В результате, требуется также большой объем памяти для записи азимутальных данных с достаточной точностью.
Более важно то, что достигнутый уровень техники геоуправления посредством измерений направления обеспечивает возможность работы только при управлении поворотами вверх и вниз. Однако, возникает множество случаев, в которых приходится вносить азимутальные поправки в траекторию ствола скважины во избежание выхода из продуктивной зоны.
Поэтому существует потребность в способах и методах обнаружения и анализа азимутальной зависимости направленных каротажных измерений с использованием замеров, получаемых при любых ази
- 4 007587 мутальных углах, для характеристики пласта грунта и для управления поворотом скважин во время бурения с повышенной точностью.
Также существует потребность в определении азимута наслоения по результатам направленных измерений, а также генерирования замеров, которые можно использовать для размещения скважин при управлении поворотом скважин вверх, вниз или по азимуту.
Также существует потребность в способах использования этих направленных измерений в реальном масштабе времени для получения расстояний до границ и для получения точных моделей грунта, при наличии которых можно принимать точные решения по геоуправлению при размещении скважин.
Также существует потребность в способе обнаружения присутствия анизотропии удельного сопротивления в слоях пластов рядом с почти вертикальными скважинами.
Кроме того, существует потребность в эффективной системе, которая обеспечивает такие направленные измерения, анализирует их по скважине и передает соответствующую информацию на поверхность для облегчения геоуправления вверх или вниз или по азимуту во время размещения скважин. Было бы также выгодно, если бы такая система могла выдавать информацию о расстоянии до границы во время фаз скольжения при бурении (т.е., когда бурильная колонна не вращается), а также в случае, когда система и/или инструмент вращается.
Определения
Определения некоторых терминов приводятся по всему тексту описания при их первом употреблении, а определения некоторых других терминов, употребляемых в этом описании, приведены ниже.
«Схождение» означает условие, когда значения, получаемые в процессе итеративных вычислений, приближаются к наблюдаемым значениям или конечным пределам с увеличением количества циклов итерации.
«Кроссплот» (график взаимной зависимости) означает график или построение такого графика, который показывает взаимную зависимость между двумя разными измерениями одного и того же объекта или выборки.
«Инверсия» или «инверсная» означает получение модели (известной под названием «инверсной модели») на основании измеренных данных (например, данных каротажа), которая описывает подземный пласт и согласуется с измеренными данными.
«Торец (торцевая поверхность) инструмента» относится к угловой ориентации инструмента относительно его продольной оси и представляет угол, образованный между выбранной поверхностью отсчета на корпусе инструмента (например, на утяжеленной бурильной трубе) и либо стенкой ствола скважины, крайней сверху в соответствии с направлением силы тяжести, либо географическим севером.
В том смысле, в каком он употребляется в данном описании, термин «симметрия» или «симметричная» относится к конфигурации, при которой наборы передающих и приемных устройств установлены в противоположных ориентациях вдоль продольной оси инструмента, так что эти наборы передатчиков и приемников можно коррелировать со стандартной симметричной операцией (например, переносом, отображением в зеркальной плоскости, инверсией и вращением) относительно точки на оси инструмента или плоскости симметрии, перпендикулярной оси инструмента.
Краткое изложение сущности изобретения
В одном аспекте, в настоящем изобретении предложен новый способ оперативной обработки данных для выделения сигналов на основании азимутального изменения направленного измерения, подходящих для характеристики пластов и использования в целях геоуправления. Вместо помещения данных каротажа в элементы памяти азимутальных измерений (накопители) и последующего определения значений, связанных с направлениями вверх и вниз, как это делается при традиционном отображении и предлагается другими исследователями, настоящее изобретение обладает преимуществом простоты физического представления отклика каротажа. Более конкретно, соответствующие сигналы, характеризующие границы, анизотропию и разрыв, выделяются из отклика пласта посредством аппроксимации азимутального изменения измеренных напряжений несколькими синусоидальными функциями. В результате этого также получается ориентация наслоений. Эта оперативная обработка повышает точность измерений, потому что используются данные, полученные во всех азимутальных направлениях или при любых азимутальных углах, а не осуществляется квантование в соответствии с конкретными углами. Такая объемная обработка становится возможной посредством целочисленных вычислений в цифровом процессоре сигналов (ЦПС), что также является новшеством в контексте рассматриваемой реализации.
За счет выделения релевантных связей по напряжению в соответствии с их порядком азимутальной зависимости, в настоящем изобретении обеспечивается возможность вычислений замеров в режиме распространения ЭМ волн с помощью лишь одной пары передатчика и приемника. Это также допускает суперпозицию замеров, полученных с помощью разных пар ПЕ-ПР, даже имеющих разную азимутальную ориентацию, для генерирования других замеров с уникальным свойством, таких, как симметризованные или антисимметризованные замеры.
В настоящем изобретении также предложена новая последовательность измерений для увеличения толерантности к неупорядоченному вращению, накоплению гармоник и даже скачкообразному перемещению. Используется быстрая последовательность операций, и эта последовательность рандомизируется
- 5 007587 в каждом цикле сбора данных для сокращения циклического просмотра при конкретных скоростях вращения.
Следовательно, настоящее изобретение можно описать как способ получения характеристик подземного пласта, начинающийся этапом, на котором размещают подходящий каротажный прибор внутри ствола скважины. Этот каротажный прибор оснащен, по меньшей мере, первыми передающей и приемной антеннами, отстоящими на некоторое первое расстояние друг от друга. По меньшей мере, одна из первых антенн является щелевым симметричным вибратором, наклоненным относительно продольной оси прибора. Первые антенны ориентированы относительно оси каротажного прибора таким образом, что, по меньшей мере, один наклоненный щелевой симметричный вибратор соответствует первому азимутальному углу. Каротажный прибор выполнен с возможностью азимутального вращения внутри ствола скважины, например, посредством вращения утяжеленной бурильной трубы или инструмента, расположенного на бурильной колонне и содержащего прибор. Когда каротажный прибор вращается, включается первая передающая антенна для передачи электромагнитной энергии в пласт. Кроме того, при вращении каротажного прибора с помощью первой приемной антенны осуществляется направленное измерение группы первых сигналов напряжения, связанных с передаваемой электромагнитной энергией, в зависимости от азимутальной ориентации каротажного прибора. Эти направленные измерения определяют азимутальное изменение измеренных первых сигналов напряжения. Это азимутальное изменение подгоняется к аппроксимирующим функциям. Этапы включения, измерения и аппроксимации можно повторять для выполнения последующих циклов сбора данных.
В конкретном варианте осуществления этап аппроксимации проводят во время измерения первых сигналов напряжения, и прекращают эту аппроксимацию, когда достигнуто схождение. Аппроксимирующие коэффициенты предпочтительно определяют с помощью быстрого преобразования Фурье.
В конкретном варианте осуществления, аппроксимирующие функции являются синусоидами, определяемыми путем связывания компонентов векторов ориентации щелевого симметричного вибратора первой передающей антенны и первой приемной антенны. Коэффициенты аппроксимирующих составляющих предпочтительно являются функциями параметров пласта грунта, включая, по меньшей мере, один из таких параметров, как удельное сопротивление слоев пласта, место нахождения каротажного прибора, отклонение ствола скважины, азимутальный угол в месте нахождения каротажного прибора, а также их комбинацию. Аппроксимирующие коэффициенты преимущественно включают в себя постоянные, а также члены, содержащие δίηφ. сокф, §ίη2φ и со§2ф, которые определяют итеративный алгоритм аппроксимации, используемый для определения азимутальной зависимости направленных измерений.
Как упоминалось выше, настоящее изобретение можно приспособить к суперпозиции замеров, получаемых с помощью разных пар передатчиков и приемников («ПЕ-ПР»). Поэтому в конкретном варианте осуществления каротажный прибор также оснащен вторыми передающей и приемной антеннами, отстоящими на первое расстояние друг от друга. Вторая передающая система представляет собой щелевой симметричный вибратор, наклон которого соответствует наклону первой приемной антенны, а вторая приемная антенна представляет собой щелевой симметричный вибратор, наклон которого соответствует наклону первой передающей антенны, так что по меньшей мере одна из вторых антенн представляет собой наклоненный щелевой симметричный вибратор. Вторые передающие и приемные антенны ориентированы относительно оси каротажного прибора таким образом, что, по меньшей мере, один наклоненный щелевой симметричный вибратор соответствует второму азимутальному углу. Таким образом, когда каротажный прибор вращается, включается вторая передающая антенна для передачи электромагнитной энергии в пласт, и с помощью второй приемной антенны осуществляется направленное измерение второй группы сигналов напряжения в зависимости от азимутальной ориентации каротажного прибора. Эти направленные измерения определяют азимутальное изменение измеренных вторых сигналов напряжения. Как и в случае измеренных первых сигналов напряжения, азимутальное изменение измеренных вторых сигналов напряжения подгоняется под аппроксимирующие функции.
В конкретном варианте осуществления второй азимутальный угол отличается от первого азимутального угла, по существу, на 90°. В альтернативном варианте, второй азимутальный угол может быть, по существу, равным первому азимутальному углу.
В конкретном варианте осуществления, аппроксимирующие функции являются синусоидами, определяемыми путем связывания компонентов векторов ориентации щелевого симметричного вибратора первой передающей антенны и первой приемной антенны, а также путем связывания компонентов векторов ориентации щелевого симметричного вибратора второй передающей антенны и второй приемной антенны. Коэффициенты аппроксимирующих составляющих предпочтительно являются функциями таких параметров пласта грунта, как удельное сопротивление слоев пласта, место нахождения каротажного прибора, отклонение ствола скважины и азимутальный угол в месте нахождения каротажного прибора. Аппроксимирующие коэффициенты преимущественно включают в себя постоянные, а также члены, содержащие δίηφ. сокф, §ίη2φ и со§2ф, которые определяют итеративный алгоритм аппроксимации, используемый для определения азимутальной зависимости направленных измерений. Измеренные первые и вторые сигналы напряжения предпочтительно являются комплексными сигналами напряжения. Соответ
- 6 007587 ственно, в этом конкретном варианте осуществления предлагаемый способ дополнительно включает в себя этапы, на которых вычисляют значения фазового сдвига и затухания на основании аппроксимирующих коэффициентов для измеренных первых и вторых сигналов напряжения, и объединяют вычисленные значения фазового сдвига и затухания для измеренных первых и вторых сигналов напряжения, формируя симметризованный или антисимметризованный замер. Значения фазового сдвига и затухания можно получать, беря логарифм отношения комплексных сигналов напряжения, получаемых из выражения аппроксимации, при двух азимутальных углах, предпочтительно - при азимутальных углах, которые составляют 0 и 180° от определенного азимута наслоения.
В другом аспекте, в настоящем изобретении обеспечивается получение характеристик шума в направленных измерениях. Так, в конкретных вариантах осуществления, в которых лишь одна из первых антенн является наклоненным щелевым симметричным вибратором, можно определить характеристики шума измеренных первых и вторых сигналов напряжения с использованием коэффициентов гармоник второго порядка. В конкретных вариантах, в которых каждая из первых антенн является либо наклоненным, либо поперечным щелевым симметричным вибратором, а аппроксимирующие коэффициенты включают в себя коэффициенты гармоник третьего порядка, и можно определить характеристики шума измеренных первых и вторых сигналов напряжения с использованием коэффициентов гармоник третьего порядка. В других случаях, можно определить характеристики шума измеренных первых и вторых сигналов напряжения посредством объединения первых и вторых измеренных сигналов напряжения.
Настоящее изобретение также можно приспособить к условиям, в которых каротажный прибор не вращается, например, во время остановки вращения бурильной колонны в процессе направленного бурения с использованием узла гидравлического забойного двигателя. При таких условиях, вторые антенны предпочтительно являются симметричными относительно первых антенн. Азимут интересующего слоя пласта определяют путем совместного использования связей первых и вторых антенн (описанных выше) и определения постоянных и коэффициентов гармоник первого порядка на основании первых и вторых сигналов напряжения, снятых, когда прибор не вращается. Определенные коэффициенты можно впоследствии использовать для выполнения этапа аппроксимации, когда прибор снова будет вращаться.
В другом аспекте, в настоящем изобретении предложены способ и система для использования выделенных направленных сигналов с целью получения расстояний до границ и принятия решений по геоуправлению. Выделенные направленные сигналы используются для выделения расстояний до границы двумя способами. При простом сценарии с одной границей или с известным профилем удельного сопротивления пласта используют простой график взаимной зависимости (кроссплот), тогда, как при более сложных ситуациях, а также для построения соответствующих им структурных моделей, используют способы инверсии. Двумерный кроссплот охватывает измерения, чувствительные к удельному сопротив лению, и измерение, которое чувствительно к расстоянию, если удельное сопротивление уступа в месте искривления скважины известно. Альтернативой является использование кроссплота двух направленных измерений для определения расстояния и удельного сопротивления уступа в месте искривления скважины. Можно также генерировать трехмерный кроссплот для одной ситуации с границами, в которой можно определить каждый из таких параметров, как удельное сопротивление пласта, удельное сопротивление уступа в месте искривления скважины и расстояние до границы. При способе инверсии используют множество инвертируемых измерений для различных моделей и определения наилучшей аппроксимации. Инверсия может быть проведена внутри скважины или на поверхности.
Частью системы визуализации информационной модели, а также результатов измерения и инверсии, является графический интерфейс пользователя (ГИП). ГИП облегчает интерактивное определение параметров инверсии, выбор измерений и моделей для улучшения интерпретации и генерирования удовлетворительных структурных моделей.
Алгоритм итеративной аппроксимации предпочтительно включает в себя этапы, на которых выполняют следующие действия: инициализировать Р0 и и0; для т = от 1 до N выборок „ „ Л, , · , · гт , · Рт , п >__ п /и-1 т—1 т—1____т—1 ' Έη-1 , η Т + Ги_! Л-! -ГтЧ •Л
следующее т; возврат (и) ;
где N выборок - суммарное количество выборок, полученных в одном цикле;
М - размер вектора аппроксимирующих функций (количество аппроксимирующих функций);
- 7 007587 и - вектор аппроксимирующих коэффициентов, имеющий размер М;
г - вектор значений аппроксимирующих функций в каждом положении измерения, имеющий размер М; и
Р - матрица размером М х М.
В одном варианте осуществления итеративный алгоритм аппроксимации обеспечивает определение того, меньше ли ошибка аппроксимации, чем предварительно определенное пороговое значение, и сходится ли и к значению, которое представляется аппроксимирующими коэффициентами.
В одном варианте осуществления, в итеративном алгоритме аппроксимации используется целочисленная реализация.
Целочисленная реализация предпочтительно используется тогда, когда каротажный прибор претерпевает азимутальное вращение с относительно большими скоростями и значительное количество каналов отклика требуют аппроксимацию.
В одном варианте осуществления, аппроксимирующие коэффициенты используются для определения ориентации слоя пласта. Измеренные первые и вторые сигналы напряжения предпочтительно являются комплексными сигналами напряжения. Ориентация слоя пласта относительно опорного значения азимутального угла для каждого канала направленного измерения определяется в соответствии с формулой С15(^'^)|),
С\с(вПЕПР) где С18 - действительная или мнимая часть коэффициента δίπφ, а С - действительная или мнимая часть коэффициента еозф из аппроксимирующего выражения. Общий азимутальный угол для первых и вторых сигналов напряжения можно вычислить с использованием взвешенного усреднения аппроксимирующих коэффициентов для действительных и мнимых частей измеренных сигналов напряжения. Амплитуду и фазу измеренного сигнала напряжения можно вычислить в предполагаемом направлении нормали к интересующей границе слоя. Фазовый сдвиг и затухание можно определить путем проведения измерений распространения для двух азимутальных углов, наприме, Фслоя и Фслоя+180о. Сигналы, основанные на аппроксимирующих коэффициентах для первых и вторых измерений напряжения, предпочтительно объединяют, чтобы получить сигналы, необходимые для определения расстояния до интересующих границ слоев.
Другие аспекты настоящего изобретения относятся к получению характеристик пластов посредством использования кроссплотов. Один способ включает в себя этап построения кроссплота двух направленных каротажных измерений, осуществляемых с помощью прибора, расположенного в стволе скважины, пересекающем пласт, для получения расстояния, по меньшей мере, до одной границы пласта и удельного сопротивления, по меньшей мере, одного слоя пласта. Построение кроссплота проводят с использованием модели с одной границей. Получаемое удельное сопротивление является удельным сопротивлением на границе раздела «слой - уступ в месте искривления скважины», а полученное расстояние является наименьшим расстоянием до границы раздела «слой - уступ в месте искривления скважины».
Еще один способ включает в себя этап построения кроссплота удельного сопротивления и направленного измерения, определяемых с помощью прибора, расположенного в стволе скважины, пересекающем пласт, для получения расстояния по меньшей мере до одной границы пласта и удельного сопротивления по меньшей мере одного слоя пласта. Построение кроссплота проводят с использованием модели с одной границей. Получаемое удельное сопротивление является удельным сопротивлением на границе раздела «слой - уступ в месте искривления скважины», а полученное расстояние является наименьшим расстоянием до границы раздела «слой - уступ в месте искривления скважины».
Другой способ включает в себя этап построения кроссплота удельного сопротивления и двух направленных каротажных измерений, осуществляемых с помощью прибора, расположенного в стволе скважины, пересекающем пласт, для получения расстояния по меньшей мере до одной границы пласта и удельного сопротивления по меньшей мере двух слоев пласта. Построение кроссплота проводят с использованием модели с одной границей. Получаемые удельные сопротивления являются удельными сопротивлениями слоя и границы раздела «слой - уступ в месте искривления скважины», а полученное расстояние является наименьшим расстоянием до границы раздела «слой - уступ в месте искривления скважины». Определенные расстояние до границы и удельное сопротивление слоя можно использовать для принятия решений по бурению.
Еще один аспект настоящего изобретения относится к использованию способа инверсии для интерпретации направленных измерений в приложениях, связанных с геоуправлением. В данном случае используется итеративный алгоритм аппроксимации для выбранных направленных измерений в реальном масштабе времени, находящий применение при геоуправлении. Для выбранных направленных измерений в реальном масштабе времени выбирают подходящую инверсную модель. Сразу же после проверки выбранной модели на предмет согласования с другой информацией, эту модель используют для принятия решений по бурению.
- 8 007587
Этап выбора модели предпочтительно включает в себя просмотр совокупности моделей следующих типов:
однородная изотропная (единственный параметр: удельное сопротивление);
однородная анизотропная (два параметра: Кг и Кв);
описывающая изотропный пласт с единственной границей, расположенной выше или ниже (три параметра: Кслоя, Куступа и расстояние до границы);
описывающая анизотропный пласт с единственной границей, расположенной выше или ниже (четыре параметра: Кслоя_г, Кслоя_в, Куступа и расстояние до границы);
описывающая три параметра изотропного пласта с двумя границами (пять параметров: Кслоя, Куступа_вверху, Куступа_внизу и расстояния до границ, расположенных выше и ниже инструмента);
описывающая три параметра анизотропного пласта с двумя границами (шесть параметров: Кслоя_г, Кслоя_в, Куступа_вверху, Куступа_внизу и расстояния до границ, расположенных выше и ниже инструмента).
В предпочтительном варианте этап выбора модели дополнительно включает в себя визуализацию выбранных направленных измерений.
В конкретном варианте осуществления, этап выбора модели предусматривает идентификацию известных параметров пласта, интерактивный выбор моделей, которые инвертируют выбранные направленные измерения, и выбор простейшей модели, которая аппроксимирует известную информацию.
Этап проверки модели предусматривает сравнение выбранной модели с известными географическими характеристиками и другими параметрами пласта, а также обновление выбранной модели, если выбранная модель не согласуется с известной информацией.
Этап обновления предусматривает уточнение выбранной модели на основании одной разновидности таких сведений, как тенденции, ранее накопленная информация, внешняя информация, а также их комбинация. Выбирают подходящие параметры инверсии и определяют диапазоны для выбранных параметров. Инверсную модель предпочтительно обновляют путем добавления большего количества слоев пластов. Некоторые из направленных измерений в реальном масштабе времени можно повторно взвешивать или исключать, а оставшиеся направленные измерения в реальном масштабе времени можно повторно инвертировать в обновленную модель.
Способы построения графиков взаимной зависимости (кроссплотов) в соответствии с настоящим изобретением предпочтительно предусматривают определение подходящей модели, выбор подходящих направленных измерений, ввод выбранных измерений в определенную модель для построения кроссплота, и создание визуального представления кроссплота. Кроссплот можно обновлять с помощью последующих измерений в реальном масштабе времени.
Еще один аспект настоящего изобретения относится к системе для измерения характеристик пластов грунта, окружающих ствол скважины. Это устройство включает в себя каротажный прибор, выполненный с возможностью размещения внутри ствола скважины. Этот каротажный прибор имеет продольную ось и оснащен парами первых и вторых передающих и приемных антенн. Пара первых передающей и приемной антенн включает в себя первую передающую антенну, представляющую собой щелевой симметричный вибратор, ориентированный в первом направлении относительно продольной оси каротажного прибора, и первую приемную антенну, находящуюся на первом расстоянии от первой передающей антенны и представляющую собой щелевой симметричный вибратор, ориентированный во втором направлении, причем первое и второе направления являются разными. Щелевые симметричные вибраторы первых передающей и приемной антенн определяют плоскость, в которой лежит продольная ось каротажного прибора. Пара вторых передающей и приемной антенн включает в себя вторую передающую антенну, представляющую собой щелевой симметричный вибратор, ориентированный во втором направлении относительно продольной оси каротажного прибора, и вторую приемную антенну, находящуюся на первом расстоянии от второй передающей антенны и представляющую собой щелевой симметричный вибратор, ориентированный в первом направлении. Щелевые симметричные вибраторы вторых передающей и приемной антенн определяют плоскость, в которой лежит продольная ось каротажного прибора. Устройство также включает в себя датчик угла торца инструмента для непрерывной индикации азимутальной ориентации каротажного прибора, а также контроллер для управления парами первых и вторых передающих и приемных антенн, обеспечивающий избирательную передачу электромагнитной энергии в пласт и измерение сигналов напряжения, связанных с передаваемой электромагнитной энергией в функции азимутальной ориентации каротажного прибора.
В конкретном варианте осуществления, пары вторых передающей и приемной антенн ориентированы под первым азимутальным углом (например, 90°) относительно пар первых передающей и приемной антенн вокруг продольной оси каротажного прибора.
В различных конкретных вариантах осуществления, измеряемые информационные характеристики включают в себя удельное сопротивление, а также информацию о геометрии пластов грунта, такую как падение, азимут и толщина слоя.
- 9 007587
Первое и второе направления изменяются и могут оказаться, например, по существу коллинеарными с продольной осью каротажного прибора, или образующими угол, по существу, 45° от продольной оси каротажного прибора.
В конкретном варианте осуществления, пары первых и вторых передающих и приемных антенн находятся в одних и тех же физических положениях на каротажном приборе.
В конкретном варианте осуществления каждый из передатчиков и приемников обладает функциональными возможностями приемопередатчика.
Датчик угла торца инструмента может предусматривать наличие магнетометров для индикации азимутальной ориентации каротажного прибора относительно земного магнитного севера, или гравитационные датчики для индикации азимутальной ориентации каротажного прибора относительно вектора силы земного притяжения.
Предлагаемое устройство предпочтительно включает в себя центральный процессор (ЦП) для обработки измеренных сигналов напряжения внутри ствола скважины, телеметрическое устройство для передачи измеренных сигналов и результатов обработки посредством ЦП из ствола скважины на поверхность, и расположенную на поверхности систему для дальнейшей обработки измеренных сигналов наряду с другими замерами для генерирования и отображения выбранных параметров согласованной модели грунта.
Краткое описание чертежей
Вышеизложенные признаки и преимущества настоящего изобретения поясняются ниже в подробном описании изобретения, кратко охарактеризованного выше, со ссылками на конкретные варианты его осуществления, проиллюстрированные на прилагаемых чертежах. Однако, следует отметить, что прилагаемые чертежи иллюстрируют лишь типичные варианты осуществления этого изобретения, так что не следует считать их ограничивающими объем его притязаний, поскольку изобретение может допускать другие, столь же эффективные варианты осуществления.
На фиг. 1 показаны условные изображения известных инструментов для индукционного каротажа или каротажа на основе распространения ЭМ волн.
На фиг. 2А и 2В - проекции, иллюстрирующие вихревые токи, индуцируемые каротажным инструментом в стволе скважины, пронизывающем пласт, при наличии и отсутствии относительного падения, соответственно.
На фиг. 3 - представление в виде вертикальной проекции обычной вращающейся бурильной колонны, в которой можно с выгодой применить настоящее изобретение.
На фиг. 4 - условное представление базового инструмента для направленных измерений, имеющего симметричные пары передающих и приемных антенн.
На фиг. 5А - условное представление инструмента для направленных измерений, имеющего конфигурацию «передатчик - приемник - приемник» («ПЕ-ПР-ПР»), которая не чувствительна к анизотропии при любом угле падения, в соответствии с одним аспектом настоящего изобретения.
На фиг. 5В показаны графики отклика направленного распространения для трехслойного пласта, полученные с помощью каротажного инструмента, соответствующего фиг. 5А.
На фиг. 6 - графическое представление способа накопления и способа оперативной аппроксимации в соответствии с одним аспектом настоящего изобретения.
На фиг. 7 - графики, представляющие схождение и ошибку в наборе целевых коэффициентов, получаемые способом аппроксимации в соответствии с одним аспектом настоящего изобретения, реализуемым посредством целочисленного алгоритма в цифровом процессоре сигналов в соответствии с другим аспектом настоящего изобретения.
На фиг. 8 - кроссплот, представляющий преобразованное измерение удельной электропроводности и направленное измерение, используемые для получения удельной электропроводности слоя пласта и расстояния до границ слоя при известном сопротивлении на границе раздела «слой - уступ в месте искривления скважины».
На фиг. 9 - результаты применения способа инверсии для получения удельного сопротивления и поточечного построения пограничных мест слоя пласта.
На фиг. 10 проиллюстрирован способ инверсии на основе кроссплота применительно к интерпретации направленных измерений удельного сопротивления.
На фиг. 11 представлена схема последовательности рабочих операций геоуправления в соответствии с одним аспектом настоящего изобретения.
На фиг. 12 - схема последовательности рабочих операций построения кроссплота, используемого для определения расстояний до границ слоев и применения таких расстояний при геоуправлении в реальном масштабе времени в соответствии с другими аспектами настоящего изобретения.
На фиг. 13 - блок-схема, представляющая структуру графического интерфейса пользователя (ГИП) для отображения моделей слоев пласта, построенных на основе инверсии.
На фиг. 14 - генерируемая компьютером визуализация инверсии геоуправления, связанная с ГИП, показанным на фиг. 13.
- 10 007587
Подробное описание изобретения
На фиг. 3 изображена обычная буровая установка и бурильная колонна, в которой можно с выгодой применить настоящее изобретение. Агрегат 10 наземной платформы и буровой вышки расположен над стволом 11 скважины, пронизывающим подземный пласт П. В изображенном варианте осуществления, ствол 11 скважины образован путем вращательного (роторного) бурения хорошо известным способом. Вместе с тем, специалистам в данной области техники, на основе информации из данного описания, должно быть понятно, что настоящее изобретение также находит такое же применение в приложениях, связанных с направленным бурением, как и в приложениях, связанных с вращательным бурением, и не ограничивается применением наземных буровых установок.
Внутри ствола 11 скважины подвешена бурильная колонна 12, которая включает в себя буровое долото 15, расположенное на ее нижнем конце. Вращение бурильной колонны 12 осуществляется посредством роторного стола 16, электропитание на который подают средства, не показанные на чертеже, и который сочленен с ведущей бурильной трубой 17, находящейся на верхнем конце бурильной колонны. Бурильная колонна 12 свешена с крюка 18, прикрепленного к перемещаемому блоку (не показан) через посредство ведущей бурильной трубы 17 и вращательного вертлюга 19, который обеспечивает вращение бурильной колонны относительно крюка.
Буровой раствор или глинистый раствор 26 хранится в колодце 27, вырытом на буровой площадке. Насос 29 подает буровой раствор 26 внутрь бурильной колонны 12 через отверстие в вертлюге 19, вынуждая протекание бурового раствора вниз по бурильной колонне 12, как показано направленной стрелкой 9. Буровой раствор выходит из бурильной колонны 12 через отверстия в буровом долоте 15, а затем циркулирует вверх через область между внешней стороной бурильной колонны и стенкой ствола скважины, называемую кольцевым пространством, как показано направленными стрелками 32. Таким образом, буровой раствор смазывает буровое долото 15 и переносит обломки выбуренной породы на поверхность, где возвращается в колодец 27 для рециркуляции.
Бурильная колонна 12 также включает в себя забойный узел, обозначенный как единое целое позицией 34 и находящийся около бурового долота 15 (иными словами, в пределах нескольких секций утяжеленной бурильной трубы от бурового долота). Забойный узел обладает функциональными возможностями измерения, обработки и хранения информации, а также осуществления связи с поверхностью. Забойный узел 34, который включает в себя, помимо прочих компонентов, устройство 36 для измерения и локальной связи, предназначенное для определения и передачи данных удельного сопротивления пласта П, окружающего ствол 11 скважины. Устройство 36 связи, также известное под названием «инструмент для каротажа по методу сопротивления», включает в себя пару первых передающей и приемной антенн ПЕ, ПР, а также пару вторых передающей и приемной антенн ПЕ', ПР' . Пара вторых антенн ПЕ', ПР' симметрична относительно пары первых антенн ПЕ, ПР, как подробнее описано ниже. Инструмент 36 для каротажа по методу сопротивления также включает в себя контроллер для управления сбором данных, как известно в данной области техники.
Забойный узел (ЗУ) 34 также включает в себя приборы, заключенные внутри утяжеленных бурильных труб 38, 39, для выполнения различных других измерительных функций, таких, как измерение естественной радиации, плотности (посредством гамма- или нейтронного каротажа) и порового давления пласта П. По меньшей мере, некоторые из утяжеленных бурильных труб оснащены стабилизаторами 37, как известно в данной области техники.
ЗУ 34 также включает в себя субузел 40 связи с поверхностью и локальной связи, расположенный как раз над утяжеленной бурильной трубой 39. Подузел 40 включает в себя тороидальную антенну 42, используемую для локальной связи с инструментом 36 для каротажа по методу сопротивления (хотя, для этого с выгодой могут быть использованы и другие известные средства локальной связи), и акустическую телеметрическую систему известного типа, осуществляющую связь с аналогичной системой (не показана) на поверхности Земли посредством сигналов, передаваемых в буровом растворе или глинистом растворе. Таким образом, телеметрическая система в субузле 40 включает в себя акустический передатчик, который генерирует в буровом растворе акустический сигнал (известный под названием «импульс бурового раствора»), отображающий измеряемые скважинные параметры.
Генерируемый акустический сигнал передается на поверхность измерительными преобразователями, обозначенными позицией 31. Измерительные преобразователи, например - пьезоэлектрические измерительные преобразователи, преобразуют принимаемые электрические сигналы в электронные сигналы. Выход измерительных преобразователей 31 подключен к наземной принимающей подсистеме 90, которая демодулирует переданные сигналы. Далее, выход принимающей подсистемы 90 подключен к процессору 85 компьютера и регистрирующему прибору 45. Процессор 85 можно использовать для определения профиля удельного сопротивления пласта (помимо прочих характеристик) «в реальном масштабе времени» во время каротажа или после него посредством осуществления доступа к зарегистрированным данным из регистрирующего прибора 45. Процессор компьютера подключен к монитору 92, в котором используется графический интерфейс пользователя (ГИП), посредством которого измеренные скважинные параметры и конкретные результаты, полученные на их основании (например, профили удельного сопротивления), предоставляются пользователю в графической форме.
- 11 007587
Предусмотрена также наземная передающая подсистема 95 для приема команд, вводимых пользователем (например, посредством ГИП в мониторе 92), выполненная с возможностью избирательного прерывания работы насоса 29 так, как предписывают измерительные преобразователи 99 в субузле 40. Следовательно, между субузлом 40 и наземным оборудованием осуществляется двухсторонняя связь. Подходящий субузел 40 более подробно описан в патентах США №№ 5235285 и 5517464, причем оба они переуступлены обладателю прав на настоящее изобретение. Специалистам в данной области техники должно быть понятно, что для осуществления связи с поверхностью можно применять альтернативные акустические способы, а также другие телеметрические средства (например, электромеханические и электромагнитные).
Азимутальная зависимость направленного измерения и новый эквивалент каротажа на основе распространения ЭМ волн
Для реализации измерений с направленной чувствительностью используют катушечные антенны двух типов. В антеннах одного типа присущая им направленная чувствительность достигается либо за счет сдвига антенн, например, от центра продольной оси каротажного инструмента, либо за счет их частичного складывания. Направленные измерения также можно проводить с помощью антенны, конфигурация которой такова, что ее магнитный момент не ориентирован вдоль продольной оси инструмента, на котором установлена антенна. Настоящее изобретение связано с направленно чувствительной антенной второго типа.
На фиг. 4 условно изображен базовый инструмент 36 для каротажа по методу сопротивления с целью направленного измерения посредством электромагнитной (ЭМ) волны. Инструмент 36 включает в себя передающую антенну ПЕ, которая излучает ЭМ волну некоторой частоты £, и приемную антенну ПР, находящуюся на некотором расстоянии Ь от передающей антенны. Также показана симметричная пара (ПЕ', ПР'), соответствующая положениям публикации № 2003/0085707 заявки на патент США (МтегЬо и др.), переуступленной обладателю прав на настоящее изобретение. Для ясности и простоты, нижеследующие рассуждения будут ограничены передающей антенной ПЕ и приемной антенной ПР, хотя они применимы и к симметричной паре антенн ПЕ' , ПР' . Следует отметить, что, хотя наклоненные моменты обеих симметричных пар находятся в одной и той же плоскости, как показано на фиг. 4, это не является обязательным признаком настоящего изобретения. Как очевидно из нижеследующего описания, для достижения эквивалентных результатов можно также ввести сигналы из обеих пар, моменты которых лежат в разных плоскостях, если на операции симметризации используются выделенные коэффициенты или направленный фазовый сдвиг или направленное затухание.
В процессе работы приемная антенна ПР регистрирует напряжение VПР-ПЕ, индуцируемое ЭМ волной от передающей антенны ПЕ и ее вторичными токами, генерируемыми в пласте, пронизываемом стволом скважины, содержащим каротажный прибор 36. Обе антенны - ПЕ и ПР - закреплены на инструменте 36 и поэтому вращаются вместе с инструментом. Можно предположить, что ориентации антенн соответствуют углу 0ПЕ для передающей антенны ПЕ и 0ПР для приемной антенны ПР. Тогда азимутальное изменение напряжения связи во время вращения инструмента может быть выражено членами связи в декартовых координатах щелевых симметричных вибраторов следующим образом:
Упр-пе(Ф) = [УггсозвПЕсозвПР+^- (УхХ+Ууу) 3ίηθΠΕ3ίηθΠρ] +
- Со (&ПЕ/ θχρ) +С1с (Θπε г θπρ) СОзф+СΠΕ г @Пр) 31Пф+ + Сгс(бпЕгдпр) соз2ф+С2В(@пЕгвпр) 3Ϊη2φ. (1.1.) где набор комплексных коэффициентов С0, С, С18, С, С определяется таким образом, что они представляют амплитуды разных составляющих измеренного отклика пласта. Следовательно, комплексные коэффициенты определяются таким образом:
- 12 007587
С (&Т1Ег @Пр) ~ [УхгЗϊηθπΕ^ΟεθπρΊ'УгхСО 3 Θπε3 ΐ ^θπρ ] /
С1з (ΘΠΕ, ΘΠΡ) = [νγζΞίηθπΕ^0ΞθΠρ+νζγ003θΠΕ3ίηθΠρ] ;
Сгс (Θπε, θηρ) = [ — (Цхх-Цуу) 3ίηθπΕ3Ϊηθπρ] г
.(..)
В соответствии с одним аспектом настоящего изобретения предполагается, что эти коэффициенты являются функциями удельного сопротивления пласта, отклонения ствола скважины и азимутального угла в месте нахождения инструмента.
При операции симметризации, т.е. (Θπε θπρ) происходит упрощение уравнений (1.1), принимающих следующий вид:
V (ф) = Упе-пр(Ф, Θπε/ θπρ)-Упр-пе (ф, Θπρ, ΘΠε) — = 2 [νχζζχ] 3ίη (Θπε-Θπρ) созф+2 [νγζζγ] 3ίη (ΘπΕ-Θπρ) 3ίηφ = = С ία (Θπε, θπρ) созф+С 13 (Θπεг Θπρ) 3Ϊηφ. (1·3)
Все гармоники второго порядка (С, С) исчезают после вычитания, потому что они симметричны относительно изменения углов наклона передатчиков и приемников. Таким образом, симметризация упрощает азимутальное изменение антисимметризованного сигнала.
На этой стадии точка отсчета азимутального угла является произвольной. Применительно к плоской геометрии, если мы выберем точку отсчета угла ф на направлении нормали к плоскости наслоения, то будем иметь Ууг=0 в результате симметрии, а функция V (ф) должна зависеть только от созф. В реальном приложении, ориентация наслоения неизвестна. Однако, задавшись любой плоскостью отсчета, можно вычислить ориентацию наслоения по формуле
При повороте на угол фнаслоения, вектор х будет направлен по нормали к наслоению, поэтому величина V (ф) будет точно равна [Ухг-Угх], за исключением постоянного множителя 2з1п (Θπεπ^ .
Сразу же после определения напряжения в каждой приемной катушке, индуцируемого каждой из передающих катушек, можно определить общий замер путем суммирования напряжений в случае индукционного инструмента или путем получения комплексного отношения напряжений в случае инструмента для каротажа на основе распространения ЭМ волн. Например, в случае показанного на фиг. 4 устройства для каротажа на основе распространения ЭМ волн, абсолютное значение напряжения в каждом приемнике можно получить в виде корня квадратного из суммы квадратов действительных и мнимых частей комплексного напряжения (выражаемого уравнением 1.1), а отношение абсолютных значений дает затухание, на основании которого можно получить удельное сопротивление Кзатух в функции затухания (удельное сопротивление пластов на относительно большой глубине исследования вокруг приемников). Фазу для каждого приемника получают, находя арктангенс отношения действительной и мнимой частей комплексного напряжения, а фазовый сдвиг представляет собой разность фаз между двумя приемниками. В таком случае, можно получить удельное сопротивление Кфс в функции фазового сдвига (удельное сопротивление пластов на относительно небольшой глубине исследования вокруг приемников).
В случае измерений посредством каротажа на основе распространения ЭМ волн находят разность (или отношение) логарифмических значений напряжений между двумя измерениями. В соответствии с теоретическими положениями, приведенными в патенте МшегЬо и др., находим амплитуду азимутального отклика, т.е. разность значений фазового сдвига и затухания при измерении под углами ф и (ф+180) , оцениваемую на максимуме отклика по напряжению. Это приводит к следующему приближенному уравнению, выведенному из уравнений (1.1-1.2):
УпР-ПЕ (ф) /УпР-ПЕ (180 + ф) - [Со (ΘΠΕ, θπρ) + С]_с (Θπε,Θπρ) созф + + С13ПЕ,0ПР) ΞΪηφ + СΠΕ,θπΡ) соз2ф + С (Θρε^Θπρ) 3ίη2φ] / [Со (Θπε/ΘΠρ) -С1с (Θπε,Θπρ) созф-С13 (Θπε/Θπρ) 3Ϊηφ + С (Θπε, Θπρ) соз2ф + С (θρρ, ΘΠΡ) 3Ϊη2φ] = 1 + 2 [Со (Θπε, Θπρ) + СосΠΕ, Θπρ) созф + + С13ηΕΠρ) 3ίηφ] / [Со (Θπε/Θπρ) + СΠΕ г θπρ) соз2ф + + С (Θπε,θΠρ) 3Ϊη2φ] = = 1 + 2 {[Ухгз1п0ПЕсозвПр + Угхсоз0ПЕз1п0пр] созф + + [νγζ3ίηθΠΕοο3θπρ+νζγαο3θΠΕ3ίηθπρ] з!пф} / {Уггсоз0пЕ31п0пр + + [Ухх + Ууу] ΞΪηθπεΞίηθπρ + [νγχ + νχγ] 31ηθπΕ3ίηθΠρ3ίη2φ + + [νχχγγ] 3ίηθΠΕ3ίηθΠρσο32φ] } . (1.5)
Максимум величины |ν| достигается при φ=0, если вектор х выбран вдоль направления нормали к наслоению. Если проводить оценку при угле φ=0, то уравнение (1.5) дает
УпР-ПЕ (0)/νΠρ-πε (180) — 1 + 2 [Со (Θπεг Θπρ) + С1С (Θπε, Θπρ) ] / [Со (ΘΠΕ, Θπρ) + С (Θπεг Θπρ) ] = = 1 + 2 [νχζ3ίηθΠΕοο3θπρ + Усоз0пез1п0пр] / [Уг2соз0пЕсоз0Пр + + Ухх31П0пе31П0пр] . (1.6)
Однако, это не «чистый» χζ-ζχ тип откликов, которые желательны, т.е. которые нечувствительны к анизотропии наслоений и углу падения.
Настоящее изобретение относится к направленным измерениям, которые нечувствительны к анизотропии пласта в широком диапазоне углов падения и в широком диапазоне частот. Как упоминалось выше, конкретные варианты осуществления изобретения основаны на антисимметризованных конфигурациях или системах антенн. Теперь, при проведении процедуры симметризации (0ПЕ 0ПР), описанной МтегЬо и др., имеем [Упр-пе (0 , Θπε, Θπρ) /Упр-пе (180, Θπε, Θπρ) ] ~ ~ [Упр-пе (180 , θπρ, Θπε)/Упр-пе (0, θπρ, Θπε) 1 = 1+2 [ [УХ2] 3ίη (Θπε~ θΠρ) 1 /[Уг2соз0ПЕсоз0ПР + Уххз1п0Пе31п0пр] . (1.7)
Это выражение вновь подобно отклику прибора индукционного типа, хотя в знаменателе попрежнему есть члены, которые не представляют простой тип [χζ-ζχ]. Это доказывает, что процедура симметризации для измерения посредством каротажа на основе распространения ЭМ волн может давать отклики, аналогичные откликам в случае симметризованного измерения посредством индукционного каротажа, но - не «чистого» типа. Также справедливо утверждение, что измерение посредством каротажа на основе распространения ЭМ волн можно проводить при двух произвольных ориентациях в азимутальном отклике.
Окончательный отклик, получаемый при этом анализе, содержит составляющие, полученные вследствие двух разных типов связей, которые на первый взгляд кажутся нежелательными в свете обычных методов каротажа. Однако, этот отклик позволяет осуществить усовершенствованный способ измерения посредством каротажа на основе распространения ЭМ волн, который прост и значительно ближе к индукционному каротажу. Предлагаемое направленное измерение реализуется за счет того, что разные типы связей естественно разделяются на разные азимутальные зависимости. Для проведения «более чистого» направленного измерения используются коэффициенты С0, С, С18.
Можно предположить значение фнаслоен некоторой определенной ориентации наслоения с помощью уравнения (1.4), так что все углы φ и X, у упоминаются в связи с этим направлением. В этом случае уравнение (1.1) упрощается, принимая следующий вид:
- 14 007587 УПР-ПЕ (^) ~ ΪΥζζ С08 &ПЕ С08 &ПР + 2 + УуУ ) 8^П ^ПЕ 5^П @ПР ] + + [УХ2 ЗШ ΘΠΕ СОЗ ΘΠΡ + Уа СОЗ ΘΠΕ 81П θπρ ] СОЗ φ + + [| (Αχ - Ууу ) 8ίη ΘΠΕ δίη ΘΠΡ1008 Ц> = = Со (ΘΠΕ ,0ПР) + С1с(вПЕПР)со8ф + С 2с (ΘΠΕ, ΘΠΡ) СО5 2φ. (1.8)
Тогда измерение посредством каротажа на основе распространения ЭМ волн можно просто охарактеризовать следующим образом:
[Со (Θπε г θπρ) + С ία (Θπεζ θπρ) ] / [С о (Θπε г Θπρ) ~С 1С (Θπε г Θπρ) ] (1.9) и
[С о (&пег Θπρ) +С 2с(Θπεг Θπρ) 1 / [ С о (Θπεг Θπρ) ~С 2с(Θπε/ θΠρ) ] ·
Теперь, после симметризации первого выражения, получаем [С о (ΘΠεζ θΠρ) +С ια(Θπεг Θπρ) ] / [ С о (Θπεг Θπρ) ~С ιο(Θπεг Θπρ) ] ' [Со (Θπρ г Θπε) + С1Сαρ, Θπε) ) / [ С о (θπρ, Θπε) ~С ια (θπρ, Θπε) ] = = 1+2 [ [νχζζχ] 3ΐη (Θπε~Θπρ) ] / [У22соз0пеСО30пр+^ [νχχγγ] 3ίηθΠΕ3ίηθπρ}, (1.10) которое лишь незначительно отличается от уравнения (1.7).
Таким образом, в случае двухкатушечных («ИЕ-ПР») измерений при обоих типах каротажа - индукционном и основанном на распространении ЭМ волн - нужно проводить анализ разных составляющих на уровне напряжения УПР-ПЕ (φ). Это дает точные поведения членов, содержащих δίπφ, созф, δίπ2φ, со§2ф, которые можно выделить с помощью нового алгоритма обработки, описанного в следующем разделе.
И вновь, следует отметить, что если либо 0ПЕ=О, либо 0ПР=О, то оба члена, содержащие δίπ2φ и со§2ф, исчезают. Зависимость напряжения основана лишь на членах, содержащих δίπφ и ^δφ.
Один важный аспект измерений фазового сдвига и затухания посредством каротажа на основе распространения ЭМ волн заключается в том, что они естественным образом подходят в качестве измерений «в процессе бурения», во время которых трудно получить подробные характеристики теплового дрейфа электронной аппаратуры в скважинных условиях. Направленные измерения фазового сдвига и затухания, о которых идет речь в данном описании, дают преимущество обычного инструмента для каротажа на основе распространения ЭМ волн по методу сопротивления с компенсацией в стволе скважины, заключающееся в том, что характеристики приемной антенны и передающей антенны и дрейф электронных компонентов приемника не оказывают влияние на измерение.
Следует также отметить, что в этом изобретении процесс симметризации осуществляется с выделенными коэффициентами, которые можно обрабатывать независимо от реального относительного азимута обеих пар ПЕ-ПР. Если ориентация пары ПЕ-ПР физически изменяется посредством вращения на некоторый дополнительный угол φ0 вокруг оси инструмента, то отклик описывается точно тем же выражением, за исключением следующей замены угла φ: φ φ+φο· Это имеет важные последствия для реальных конфигураций катушек. За счет использования двух пар ПЕ-ПР с разными азимутальными ориентациями для измерений, процесс симметризации в случае вращения можно осуществлять точно так же, как в случае, если бы обе пары находились в одной и той же плоскости. Вместе с тем, когда инструмент скользит, измеренный сигнал, полученный при наличии двух ориентации, можно использовать для построения требуемых направленных измерений, просто предполагая, что азимут инструмента и наслоения не изменился с последнего момента, когда инструмент еще вращался. Таким образом, можно использовать уравнение 1.3 для получения амплитуд С и ϋ, необходимых для построения измерений посредством каротажа на основе распространения ЭМ волн.
Эти анализы можно без изменения распространить на традиционный тип измерений с помощью конфигураций ПЕ-ПР-ПР, описанных МтпегЬо и др. Специалистам в данной области техники должно быть очевидно, что эта процедура дает, по существу, такой же отклик, как указанный выше, но с удвоением сигнала, поскольку промежуток в паре приемных антенн значительно меньше по сравнению с промежутком в паре ПЕ-ПР. Направленные сигналы из двух приемных антенн просто складываются.
На фиг. 5А показана конфигурация ПЕ-ПР-ПР, которая нечувствительна к анизотропии при любом угле падения, а на фиг. 5В показаны отклики, соответствующие этой конфигурации. Осуществляется подача электропитания на передающую антенну ПЕ1 и измерение фазового сдвига и затухания, обуславливаемых приемными антеннами ПР11 и ПР12. Затем осуществляется подача электропитания на пере
- 15 007587 дающую антенну ПЕ2 и измерение фазового сдвига и затухания, обуславливаемых приемными антеннами ПР21 и ПР22. Показание инструмента соответствует разностям между этими двумя наборами замеров. Поскольку отдельные замеры в однородной среде идентичны при любом угле и при любой анизотропии, показание инструмента является нулевым в однородной среде при любом падении.
Отклики измерений в трехслойном анизотропном пласте показаны на фиг. 5В. Показание инструмента является нулевым на расстоянии от границы при любом падении, и имеет место незначительная чувствительность к анизотропии вблизи границы. Различие откликов возникает вследствие того, что отклики распространения окажутся не симметричными, если поменять местами передатчик и приемник. Проведение измерений в направлениях вверх и вниз предусматривает получение только информации в этих направлениях, даже вблизи границы. Следует отметить, что отклики затухания практически перекрываются для разных углов падений, если все антенны находятся в одной и той же среде, что просто идеально для индукционного измерения типа χζ-ζχ (описанного выше). Измерения фазового сдвига также перекрываются, хотя отклики в проводящем слое (1 См/м) удваиваются.
Алгоритмы цифрового процессора сигналов (ЦПС) для выделения сигнала на основании азимутальной зависимости
Из предыдущего анализа очевидно, что наилучший способ выделения направленных измерений заключается в анализе сигналов напряжения и выделении подходящих измерений с последующим их синтезом и симметризацией для конечного использования при проведении геоуправления. Традиционный способ обработки азимутальных данных заключается в помещении их в малые элементы памяти (накопители), соответствующие одинаковому азимутальному диапазону (см. левую часть фиг. 6), как делалось со скважинными изображениями удельного сопротивления или плотности. Затем можно идентифицировать верхний и нижний накопители, определяя местонахождение максимального и минимального значений. Однако этот метод имеет следующие недостатки:
1. Накопление снижает угловое разрешение измерения;
2. Накопление приводит к потере данных, поскольку данные, не попадающие ни в один из верхнего и нижнего накопителей, не используются для вычисления расстояния до границы;
3. В случаях заедания и скольжения накопление неинформативно, вследствие чего верхний и нижний накопители могут быть пустыми или иметь малые выборки, что обуславливает ошибку в измерениях в направлениях вверх и вниз и, возможно, неправильную идентификацию пиковых значений в наихудшем случае;
4. Накопление занимает больше памяти.
Один аспект настоящего изобретения имеет принципиальное отношение к определению коэффициентов членов, содержащих зтф, созф, δίη2φ и соз2ф, которые определяют новый алгоритм, используемый вместо обычного накопления. Этот предлагаемый метод называется «оперативной» аппроксимацией азимутального отклика в соответствии с выделениями подходящих членов, содержащих зт и соз, для направленных измерений, причем эта аппроксимация проводится итеративно (см. построенные точки на фиг. 9) . Такой алгоритм аппроксимации реализуется в ЦПС посредством целочисленного алгоритма, так что он достаточно быстро выполняется для всех каналов в пределах времени выборки, составляющего 4 мс. Точное использование информации об азимутальных углах и рандомизация последовательностей сбора данных делают этот алгоритм робастным к допустимому неравномерному вращению инструмента, а также к заеданию и скольжению в тяжелых условиях бурения. Таким образом, для получения сигнала в направлении вверх или вниз, используются все данные, а не только данные в обоих накопителях, что увеличивает отношение «сигнал - шум» при измерении. Использование точных углов азимута также позволяет точнее определять ориентацию наслоения.
Подробный алгоритм можно описать следующим образом.
Реализация в режиме с плавающей запятой: начинают с начального значения матрицы Р0 и вектора и0, затем выполняют алгоритм, описанный ниже (который также представлен графически в правой части фиг. 6), с измерением у(ф;) и базисом г = (1 зтф;, созф;, зт2ф;, соз2ф;)Т, где Р - матрица размером МхМ, а и и г - векторы размером М. М - это размер базисной функции. После итерации по N и будет сходиться к значению, которое представляет коэффициенты выражения. Этот алгоритм устойчив, так что сходимость обычно достигается за 10-15 итераций.
Подробный алгоритм показан ниже:
инициализировать Р0 и и0;
для т = от 1 до N выборок
Рт Рщ-1 ~ (Рщ-1'^Тп1-1'^т-1'Рт-1) / (1 + ^т-1'Рт-1'^Тщ-1) ит-1 ~ Рт'Г т-1’ (Ут-1 ~ Е т-1’-Т т-1) следующее т;
возврат (и) ;
- 16 007587 где N выборок - суммарное количество выборок, полученных в одном цикле;
М - размер вектора аппроксимирующих функций (количество аппроксимирующих функций);
и - вектор аппроксимирующих коэффициентов, имеющий размер М;
г - вектор значений аппроксимирующих функций в каждом положении измерения, имеющий размер М; и
Р - матрица размером М х М.
Целочисленные алгоритмы в ЦПС
Во многих случаях, реализация в режиме с плавающей запятой будет слишком дорогостоящей, чтобы ее можно было осуществить в доступных в настоящее время скважинных ЦП, потому что может потребоваться аппроксимация по сотням каналов, а сбор данных для каждого азимутального угла придется проводить за довольно короткое время (исчисляемое миллисекундами), чтобы получить точный угол на повышенной скорости вращения. В этой ситуации можно применить целочисленную реализацию, несколько модифицировав ее, чтобы повысить точность (например, можно использовать 32 разряда для умножения), провести масштабирование и избежать переполнения, а также ускорить сходимость. Значения базисной функции также можно генерировать заранее и хранить в памяти, чтобы позже интерполировать их для получения значения истинного угла ф1. В целом, целочисленные алгоритмы основаны на 16-разрядном целочисленном представлении данных с 32-разрядными целочисленными промежуточными переменными;
более эффективны, но менее точны;
приспособляемы к стратегиям повышения точности и скорости сходимости: масштабированию, инициализации и повторной инициализации;
приспособляемы к стратегиям, которые в свою очередь приспособлены к среде ЦПС: делению, округлению.
Протекание процесса сходимости для целочисленной версии алгоритма показано на фиг. 7. Следует отметить, что аппроксимация является довольно точной, а ошибки в общем случае составляют менее 1%. Очень важным моментом аппроксимации является то, что почти все данные используются для получения коэффициентов (сигналов направленного измерения), что значительно увеличивает отношение «сигнал шум». Например, если используются только 32 элемента памяти (накопителя), то при реализации от верхнего до нижнего накопителей используется лишь 1/16 часть данных. В отличие от этого, при использовании оперативной аппроксимации учитываются почти все данные (помимо части первоначальной сходимости).
Поскольку при осуществлении метода аппроксимации (описанного выше) выделяются только релевантные сигналы, нужно сохранять только полезные коэффициенты. Так, в этом случае нужно сохранять только 5 коэффициентов по сравнению с 32 в случае, если бы накапливались все данные в варианте с 32 накопителями. Специалистам в данной области техники должны быть понятны преимущества предлагаемого метода, которые включают в себя точность выделенного сигнала и конкретное увеличение точности азимутального угла.
В еще одном аспекте, изобретение обеспечивает получение характеристик шума при направленных измерениях. Так, в конкретных вариантах осуществления, в которых только пара первых антенн («ПЕПР») представляет собой наклоненный щелевой симметричный вибратор, характеристики шума измеряемых сигналов напряжения можно получить с использованием коэффициентов гармоник второго порядка. В вариантах осуществления, в которых каждая из первых антенн представляет собой наклоненный либо поперечный щелевой симметричный вибратор, а аппроксимирующие коэффициенты включают в себя коэффициенты гармоник третьего порядка, шум измеряемых сигналов напряжения можно охарактеризовать коэффициентами гармоник третьего порядка. В других случаях шум измеряемых сигналов напряжения можно охарактеризовать путем объединения упомянутых сигналов.
Настоящее изобретение также можно приспособить к условиям, в которых каротажный прибор не вращается, например, таким, которые возникают при остановке вращения бурильной колонны, когда направленное бурение проводится с помощью узла гидравлического забойного двигателя. При таких условиях пара вторых антенн (ПЕ'-ПР') предпочтительно симметрична относительно первой пары антенн. Азимут интересующего слоя пласта определяют путем объединения этих связей антенн (как описано выше) и определения постоянных коэффициентов и коэффициентов гармоник первого порядка измеренных сигналов напряжения, получаемых, когда прибор не вращается. Определенные коэффициенты затем используют для проведения итеративной аппроксимации, когда каротажный прибор снова вращается.
Расстояние до границы, используемое при геоуправлении
В соответствии с еще одним аспектом настоящего изобретения, определенные коэффициенты можно использовать для получения ориентации наслоения пласта. Ориентацию наслоения для каждого канала направленных измерений можно определить с помощью уравнения 1.4. Кроме того, можно отметить, что среднее значение для многих каналов, взвешенное посредством функции относительного уровня сигнала каждого канала, который содержит замер, также можно использовать для повышения точности,
- 17 007587 поскольку эта ориентация должна быть одинаковой для всех каналов. Тогда можно вычислить амплитуду и фазу сигнала напряжения, Упр-пе (фнаслоен) в предполагаемом направлении нормали к слою.
Направленный фазовый сдвиг и затухание можно вычислить с помощью уравнения 1.10. Затем, симметризацию осуществляют с целью получения конечного сигнала, необходимого для нахождения расстояния до границы. Вместе с тем, следует понять, что порядок проведения этих этапов можно изменять, получая при этом аналогичные или идентичные отклики.
Чтобы получить расстояние до границы в реальном масштабе времени, можно использовать два способа. В случае простых моделей (где есть только одна граница), кроссплот двух направленных измерений позволяет получить расстояние до границы и удельное сопротивление одного из слоев. Характерный кроссплот показан на фиг. 8, где используется зависимость отклика затухания при 84 дюймах, 100 кГц, от удельного сопротивления при 28 дюймах и фазовом сдвиге 2 МГц (удельное сопротивление уступа в месте искривления скважины: К.г=0,8 Ом, Вв=3,6 Ом). В данном случае симметризация, по существу, позволяет исключить из рассмотрения другие параметры, такие, как анизотропия и падение. Графики для симметризованного измерения просты. Использование других комбинаций парных измерений для получения согласованной картины увеличит достоверность интерпретации.
На фиг. 9 иллюстрируется использование инверсии на основе кроссплота для интерпретации направленных измерений. Для интерпретации измерений в случае модели с одной границей и при фиксированном удельном сопротивлении слоя на уровне 100 Ом, а также при изменяемом удельном сопротивлении на границе раздела «слой - уступ в месте искривления скважины» (ниже инструмента) и изменяемом расстоянии, используются два направленных измерения при 84 дюймах и 400 кГ ц. Отклик нечувствителен к удельному сопротивлению слоя, а сигнал определяется главным образом удельным сопротивлением на границе раздела «слой - уступ в месте искривления скважины» и расстоянием. Точки, представленные на отображенном на экране изображении, показывают замеры удельного сопротивления на границе раздела «слой - уступ в месте искривления скважины» и расстояния. Значения удельного сопротивления и расстояния считываются из кроссплота и выводятся на экран.
Для случая одной границы слоев, когда удельные сопротивления слоев по обе стороны от границы известны, можно использовать три входных сигнала, по меньшей мере, один из которых представляет собой замер удельного сопротивления среды, окружающей инструмент (например, инструмент 36, показанный на фиг. 3).
Для более сложных моделей, предусматривающих большее количество границ, применяется программа инверсии, основанная на способе, предложенном в патенте США № 6594584. Траекторию ствола скважины проецируют на первоначальную модель пласта. Каротажный инструмент перемещают в пределах сегмента траектории ствола скважины и измеряют отклики инструмента во время движения вдоль упомянутого сегмента. Также определяют ожидаемые отклики инструмента в соответствии с моделью. Разности, полученные при сравнении ожидаемых и измеренных откликов вдоль сегмента, используют затем для коррекции модели, и цикл сравнения и коррекции повторяют до тех пор, пока упомянутые разности не станут меньше некоторого выбранного значения. Робастность инверсной модели увеличивается за счет использования нескольких начальных точек и физических критериев установления различий между решениями.
Затем применяют алгоритм инверсии. Этот алгоритм предусматривает использование любых замеров в качестве входных переменных, с последующим нахождением модели, наиболее согласованной с данными. Осуществляется просмотр нескольких моделей, а наилучшая модель выбирается автоматически. Реконструкция типичного профиля удельного сопротивления показана на фиг. 10А, тогда как реконструкция соответствующей границы (структуры) изображена на фиг. 10В. Эти результаты показывают, что решения в случае приближения к проводящим слоям точнее, чем в случае приближения к резистивным слоям. Это ожидаемый результат, поскольку отклики направленных измерений характеризуют главным образом проводящий слой, и поэтому они более чувствительны к изменению удельного сопротивления в проводящих слоях, расположенных выше и ниже, чем к изменению удельного сопротивления на границах раздела «резистивный слой - уступ в месте искривления скважины».
Еще один аспект настоящего изобретения предусматривает использование определенных расстояний до границ с целью принятия решений по геоуправлению. Общий алгоритм геоуправления можно описать со ссылкой на фиг. 11. Сначала на этапе 110 выбирают направленные измерения в реальном масштабе времени способом оперативной аппроксимации, описанным выше. Затем на этапе 120 визуализируют эти направленные измерения в виде ряда логарифмических откликов или с помощью кроссплотов, как показано на фиг. 8 и 10 (и описано выше).
После этого на этапе 130 принятия решения поток обработки направляется в зависимости от того, идентифицирована ли подходящая модель. Если структура информации или какие-либо иные параметры известны, то этот подход позволяет пользователю зафиксировать известные параметры. Например, может быть известно, что удельное сопротивление или пролегание верхнего слоя на уступе в месте искривления скважины стабильны, и с высокой степенью определенности может быть известно, что ниже границы нет. Такая информация позволяет, например, выбрать модель с единственной границей и фиксированным удельным сопротивлением на границе раздела «слой - уступ в месте искривления скважины».
- 18 007587
Математически это означает, что нужно инвертировать только три параметра - горизонтальное и вертикальное удельное сопротивление слоя (Кг и Кв) и расстояние до верхней границы. Это пример инверсии полной модели (см. этап 140), которая гарантирует удовлетворительную интерпретацию и позволяет избежать путаницы в случаях ограничений, когда может быть выбрана какая-либо не физическая модель, обеспечивающая лучшую математическую аппроксимацию измерения.
В альтернативном варианте, на этапе 150, проводят быстрые аппроксимированные инверсии с множеством моделей наряду с осуществлением на этапе 160 алгоритма автоматического выбора модели. При осуществлении этого алгоритма просматривают разные модели, от простых (отсутствие границ, изотропный пласт) до более сложных (два расстояния и анизотропный пласт), включая однородную изотропную (единственный параметр: удельное сопротивление);
однородную анизотропную (два параметра: Кг и КВ);
описывающую изотропный пласт с единственной границей, расположенной выше или ниже (три параметра: Кслоя, Куступа и расстояние до границы);
описывающую анизотропный пласт с единственной границей, расположенной выше или ниже (четыре параметра: Кслоя_г, Кслоя_в, Куступа и расстояние до границы);
описывающую три параметра изотропного пласта с двумя границами (пять параметров: Кслоя, Куступа_вверху, Куступа_внизу и расстояния до границ, расположенных выше и ниже инструмента);
описывающую три параметра анизотропного пласта с двумя границами (шесть параметров: Кслоя_г, Кслоя_в, Куступа_вверху, Куступа_внизу и расстояния до границ, расположенных выше и ниже инструмента).
Решение получают для каждой модели.
При осуществлении алгоритма выбора модели используются физически обоснованные ограничения для направленных измерений и накладываются условия, которые обуславливают выбор «простейшей модели, которая аппроксимирует данные» (критерий Байеса). Для исключения моделей ввиду их сложности можно использовать классические критерии - информационный критерий Акайке (ИКА) или байесовский информационный критерий.
После этого на этапе 180 принятия решения осуществляют проверку модели, чтобы определить, согласуется ли выбранная модель с уже известной информацией о геологическом образовании или с другими измерениями, такими, как измерения, проведенные посредством гамма-каротажа, или с иными измерениями, доступными в реальном масштабе времени, что отображается этапом 170. Если модель согласуется с другими данными, то она принимается как входная (этап 190) для принятия решения по бурению (этап 200), и на этапе 210 выдается соответствующее сообщение. Реализация решения по бурению приведет к дальнейшим направленным измерениям в реальном масштабе времени (этап 100), которые снова вводятся на этапы 110 и 120 аппроксимации и визуализации. Если выбранная инверсная модель не согласуется с другими измерениями (см. этап 170), то следует провести обновление модели. В этом случае на этапе 220 применяется субалгоритм интерактивного режима, согласованный с программой комплексной инверсии, описанной в патенте США № 6594584.
Автоматическая инверсия зачастую может приводить к неудовлетворительной интерпретации. Это может быть вызвано шумом измерения, вызываемым электронными средствами, а также «шумом модели», т.е. тем, что реальная модель отличается от моделей, просматриваемых при инверсии. Например, в модели не включаются параметры вторжения в ствол скважины, тонких слоев, косой или волнистая слоистости, а также конечные размеры инструмента, что может вызвать их несоответствие в процессе аппроксимации.
Гибкость и возможность для пользователя выбирать в интерактивном режиме общую модель для сегмента данных являются основополагающими факторами для успешной интерпретации измерений. Программное обеспечение, определяющее субпоследовательность действий в интерактивном режиме (этап 220), имеет особенности, которые обеспечивают: уточнение моделей (этап 221) на основании тенденций, ранее известной информации или сведений из внешнего источника информации; ограничение или фиксацию некоторых параметров (этап 222); исключение некоторых измерений, которые могут подвергаться большему влиянию окружающей среды, не учитываемому в модели; и повторную обработку данных (этап 223).
На фиг. 12 изображена схема последовательности операций построения кроссплота направленных измерений, которые обеспечивают определение расстояний до границ, описываемое в других местах данного описания. Создание или модификация кроссплота представляется субпоследовательностью 20 операций, которая инициируется определением подходящей модели (на этапе 30). Затем выбираются подходящие направленные измерения, соответствующие возможности определения и/или прогнозирования удельного сопротивления на уступе в месте искривления скважины (этап 40). Если удельное сопротивление на границе раздела «слой - уступ в месте искривления скважины» достоверно известно, что часто бывает на крупных месторождениях, где бурят много скважин, то интерпретация основывается на одном удельном сопротивлении и одном направленном измерении (этап 50), чтобы определить истинное удельное сопротивление слоя (с поправкой, учитывающей характеристики на границе раздела «слой - уступ в месте искривления скважины») и расстояние до границы. Если удельное сопротивление на границе раз
- 19 007587 дела «слой - уступ в месте искривления скважины» неизвестно, то рекомендуется использовать кроссплот двух направленных измерений (этап 60) , как показано на фиг. 9, чтобы определить удельное сопротивление на границе раздела «слой - уступ в месте искривления скважины» и расстояние до границы. Выбранные измерения вводятся в определенную модель для генерирования откликов на этапе 70, причем эти отклики на этапе 80 могут визуализироваться как показано на фиг. 8 и 9. В случаях, когда ни одно удельное сопротивление неизвестно, а известен диапазон их разности, могут строиться трехмерные кроссплоты на основе объединения одного удельного сопротивления и двух направленных измерений.
После построения или обновления кроссплота (что иллюстрируется субпоследовательностью 20 операций), его можно непрерывно обновлять путем введения дополнительных измерений в реальном масштабе времени (этап 100) и визуализации (этап 120'). Обработка кроссплота направленных измерений дает расстояние до одной или более границ пласта и удельное сопротивление одного или более слоев (этап 220), которое можно выдавать и визуализировать (этап 230) в качестве инструкции для принятия решения по бурению на этапе 200.
В соответствии с еще одним аспектом настоящего изобретения, создан графический интерфейс пользователя (ГИП), предназначенный для облегчения использования программы и отображения инверсных моделей наслоений с целью визуализации их для пользователя. Блок-схема, иллюстрирующая структуру ГИП и разные способы обработки данных, показана на фиг. 13. Соответственно, ГИП позволяет просматривать опережающую модель для заданного наслоения;
строить кроссплоты для анализа чувствительности измерений;
осуществлять инверсию данных в реальном масштабе времени, характеризующуюся:
проведением инверсии на основе кроссплота применительно к одному удельному сопротивлению и одному расстоянию (для заданного удельного сопротивления уступа в месте искривления скважины или удельного сопротивления пласта);
автоматической интерпретацией с использованием быстрой инверсии на основании суперпозиции откликов одиночной границы и автоматического выбора модели, а также физически обоснованных ограничений, в сочетании с использованием информационного критерия Акайке или смешением моделей;
интерактивной интерпретацией, основанной на полной инверсии любого параметра (до 6) из трехмерной модели (слой и два уступа в местах искривления скважины - два расстояния и четыре удельных сопротивления), просмотром полной модели слоистой среды в контуре инверсии, обеспечением уточнений модели, предусматривающим опции выбора и отмены выбора параметров инверсии и ограничения их и/или выбора и повторного взвешивания полученных замеров;
переключением между быстрым и интерактивным режимами для уточнения модели; визуализацией измерения и инверсии.
Пример экрана для интерпретации показан на фиг. 14. На экране отображаются оба измерения, которые используются при инверсии, а также траектория ствола скважины. Графически отображаются инвертированные результаты, включающие в себя расстояние до верхней и нижней границ, удельные сопротивления в слое и двух границах раздела «слой - уступ в месте искривления скважины». Результаты показаны от точки к точке и основаны на восстановленных положениях границ. Имеется возможность делать выводы о стратиграфическом падении. Следует отметить, что, хотя опережающая модель, просматриваемая в контуре инверсии, представляет одномерную слоистую среду, применяемый подход позволяет строить и более сложные информационные изображения, включающие в себя непараллельные границы;
определение сбросов, включая субсейсмические сбросы.
Оценка анизотропии
При наличии конфигурации ПЕ-ПР-ПР с наклонной передающей антенной и, по меньшей мере, одной из приемных антенн, можно проводить измерения анизотропии в соответствии с положениями публикации № 2003/0200029 заявки на патент США. Для определения анизотропии, объединяют обнаруженные ЭМ сигналы, связанные с разностью фаз или отношением амплитуд. Для такого измерения в указанной заявке предложен новый способ компенсации ствола скважины.
Направленную часть таких объединенных измерений можно анализировать так же, как описано выше. Соответственно, можно использовать коэффициенты составляющих для вывода характеристики анизотропии посредством инверсии методом, описанным в упоминавшемся выше патенте США № 6594584. Этот вывод справедлив для всех углов падения.
Специалистам в данной области техники должно быть понятно, что в сильно отклоненных стволах скважин традиционное измерение методом распространительного каротажа с использованием осевых катушек уже обеспечивает приемлемую чувствительность к анизотропии, что исключает необходимость способа компенсации ствола скважины в этих стволах скважин.
Для специалистов в данной области техники должно быть очевидно, что настоящее изобретение может быть реализовано с помощью одного или более подходящих компьютеров общего назначения, имеющих надлежащее аппаратное обеспечение и запрограммированных на проведение процессов согласно изобретению. Программирование можно реализовать посредством использования одного или бо
- 20 007587 лее устройств для запоминания программ, считываемых процессором компьютера, и кодирования одной или более программ команд, выполняемых компьютером, для проведения вышеописанных операций. Устройство для запоминания программ может принимать форму, например, одного или более флоппидисков, СО-КОМ или другого оптического диска, магнитной ленты, микросхемы постоянного запоминающего устройства (ПЗУ), а также другие формы, которые хорошо известны в данной области техники или будут разработаны впоследствии.
Программа команд может представлять собой «объектный код», т.е. может быть выражена в двоичной форме, исполняемой более или менее непосредственно компьютером, «исходный код», который требует компиляции или интерпретации перед выполнением, или может быть представлена в некоторой промежуточной форме, такой, как частично компилированный код. Точные формы устройства для запоминания программ и кодирования команд в данном случае несущественны. Таким образом, эти средства обработки можно реализовать в оборудовании, которое расположено на поверхности, в инструменте, или они могут быть реализованы и в упомянутом оборудовании, и в инструменте, как известно в данной области техники. Способы, согласно изобретению, можно использовать совместно с любого типа системой для каротажа скважин, например, с инструментами, спускаемыми в скважину на тросе, инструментами для каротажа в процессе бурения и/или скважинных исследований в процессе бурения (КвПБ/СИвПБ), или инструментов для каротажа в процессе спуска (КвПС).
Из вышеизложенного описания должно быть понятно, что в предпочтительных и альтернативных вариантах осуществления настоящего изобретения могут осуществляться различные модификации и изменения в рамках сущности изобретения.
Это описание приведено лишь в целях иллюстрации и не должно считаться носящим ограничительный характер. Объем изобретения должен определяться только текстом нижеследующей формулы изобретения. Термин «содержащий(ая, ее)» в рамках формулы изобретения следуют считать означающим «включающий(ая, ее) в себя, по меньшей мере», так что приводимый перечень элементов в некотором пункте формулы изобретения представляет собой открытую группу. Употребление терминов в единственном числе следует считать распространяющимся и на формы множественного числа таких терминов, если специально не указано иное.

Claims (31)

1. Способ получения характеристик подземного пласта, включающий в себя этапы, на которых размещают внутри ствола скважины каротажный прибор, оснащенный по меньшей мере первыми антеннами, одна из которых является передающей, а другая приемной, отстоящими на некоторое первое расстояние друг от друга, причем по меньшей мере одна из первых антенн представляет собой щелевой симметричный вибратор, наклоненный относительно продольной оси прибора, при этом первые антенны ориентированы относительно оси каротажного прибора таким образом, что по меньшей мере один наклоненный щелевой симметричный вибратор соответствует первому азимутальному углу;
осуществляют азимутальное вращение каротажного прибора внутри ствола скважины;
когда каротажный прибор вращается, включают первую передающую антенну для передачи электромагнитной энергии в пласт;
когда каротажный прибор вращается, осуществляют направленное измерение первых сигналов напряжения, связанных с передаваемой электромагнитной энергией, с помощью первой приемной антенны в зависимости от азимутальной ориентации каротажного прибора для определения азимутального изменения измеренных первых сигналов напряжения; и осуществляют аппроксимацию азимутального изменения измеренных первых сигналов напряжения, чтобы получить функции, являющиеся его приближением.
2. Способ по п.1, при котором этап аппроксимации проводят во время измерения первых сигналов напряжения.
3. Способ по п.2, дополнительно включающий в себя этап, на котором прекращают аппроксимацию, когда достигнута сходимость.
4. Способ по п.1, при котором аппроксимирующие функции являются синусоидами, определяемыми путем связывания компонентов векторов ориентации щелевого симметричного вибратора первой передающей антенны и первой приемной антенны.
5. Способ по п. 4, при котором коэффициенты аппроксимирующих составляющих являются функциями параметров пласта грунта, включая по меньшей мере один из таких параметров, как удельное сопротивление слоев пласта, место нахождения каротажного прибора, отклонение ствола скважины, азимутальный угол в месте нахождения каротажного прибора, а также их комбинацию.
6. Способ по п.5, при котором аппроксимирующие коэффициенты включают в себя постоянные, а также члены, содержащие δίηφ, сокф, §ш2ф и со§2ф, которые определяют итеративный алгоритм аппроксимации, используемый для определения азимутальной зависимости направленных измерений.
7. Способ по п.6, при котором используют итеративный алгоритм аппроксимации для выбранных направленных измерений в реальном масштабе времени, применяемых при геоуправлении.
- 21 007587
8. Способ по п.1, при котором каротажный прибор также оснащен вторыми антеннами, одна из которых является передающей, а другая приемной, отстоящими на некоторое первое расстояние друг от друга, причем вторая передающая антенна представляет собой щелевой симметричный вибратор, наклон которого соответствует наклону первой приемной антенны, а вторая приемная антенна представляет собой щелевой симметричный вибратор, наклон которого соответствует наклону первой передающей антенны, так что по меньшей мере одна из вторых антенн представляет собой наклоненный щелевой симметричный вибратор, при этом вторые передающая и приемная антенны ориентированы относительно оси каротажного прибора таким образом, что по меньшей мере один наклоненный щелевой симметричный вибратор соответствует второму азимутальному углу, и при этом способ также включает в себя этапы, на которых когда каротажный прибор вращается, включают вторую передающую антенну для передачи электромагнитной энергии в пласт;
когда каротажный прибор вращается, осуществляют направленное измерение второй группы сигналов напряжения, передающих электромагнитную энергию с помощью второй приемной антенны в зависимости от азимутальной ориентации каротажного прибора, чтобы определить азимутальное изменение измеренных вторых сигналов; и осуществляют аппроксимацию азимутального изменения измеренных вторых сигналов напряжения, чтобы получить функции, являющиеся его приближением.
9. Способ по п.8, при котором второй азимутальный угол отличается от первого азимутального угла, по существу, на 90°.
10. Способ по п.8, при котором аппроксимирующие функции являются синусоидами, определяемыми путем связывания компонентов векторов ориентации щелевого симметричного вибратора первой передающей антенны и первой приемной антенны, а также путем связывания компонентов векторов ориентации щелевого симметричного вибратора второй передающей антенны и второй приемной антенны.
11. Способ по п.10, при котором коэффициенты аппроксимирующих составляющих являются функциями таких параметров пласта грунта, как удельное сопротивление слоев пласта, место нахождения каротажного прибора, отклонение ствола скважины и азимутальный угол в месте нахождения каротажного прибора.
12. Способ по п. 11, при котором измеренные первые и вторые сигналы напряжения являются комплексными сигналами напряжения, а способ дополнительно включает в себя этапы, на которых вычисляют значения фазового сдвига и затухания на основании аппроксимирующих коэффициентов для измеренных первых и вторых сигналов напряжения и объединяют вычисленные значения фазового сдвига и затухания для измеренных первых и вторых сигналов напряжения для формирования симметризованного или антисимметризованного результата измерений.
13. Способ по п.11, при котором лишь одна из первых антенн имеет наклоненный щелевой симметричный вибратор, а способ дополнительно включает в себя этап, на котором определяют характеристики шума измеренных первых и вторых сигналов напряжения с использованием коэффициентов гармоник второго порядка.
14. Способ по п.11, при котором каждая из первых антенн представляет собой наклоненный либо поперечный щелевой симметричный вибратор, а аппроксимирующие коэффициенты включают в себя коэффициенты гармоник третьего порядка, причем способ дополнительно включает в себя этап, на котором определяют характеристики шума измеренных первых и вторых сигналов напряжения с использованием коэффициентов гармоник третьего порядка.
15. Способ по п.11, при котором определяют характеристики шума измеренных первых и вторых сигналов напряжения посредством объединения первых и вторых измеренных сигналов напряжения.
16. Способ по п.11, при котором вторые антенны являются симметричными относительно первых антенн, а способ дополнительно включает в себя этапы, на которых временно останавливают вращение каротажного прибора;
определяют азимут интересующего слоя путем объединения связей первых и вторых антенн;
определяют постоянные коэффициенты и коэффициенты гармоник первого порядка на основании первых и вторых сигналов напряжения, снятых, когда прибор не вращается; и используют постоянные коэффициенты для выполнения этапа аппроксимации, когда каротажный прибор вращается.
17. Способ по п.16, дополнительно включающий в себя этапы, на которых обновляют азимутальные изменения измеренных первых и вторых сигналов напряжения и повторно вычисляют эквивалентные напряжения, когда инструмент находится в плоскости наслоения.
18. Способ по п.1, дополнительно включающий в себя этап, на котором вычисляют амплитуду и фазу измеренного сигнала напряжения в предполагаемом направлении нормали к границе интересующего слоя.
- 22 007587
19. Способ по п.18, дополнительно включающий в себя этап, на котором определяют фазовый сдвиг и затухание, проводя измерения посредством каротажа на основе распространения ЭМ волн для двух азимутальных углов.
20. Способ определения характеристик подземного пласта, включающий в себя этапы, на которых генерируют график, используя два типа направленных измерений в качестве области значений и области определения графика и имея удельные сопротивления пласта и расстояние до границы слоя, нанесенной на график, формируют вслед за тем координатную сетку;
выбирают два направленных каротажных измерения, получаемых от прибора, размещенного в стволе скважины, пересекающей пласт;
определяют расположение точки на графике при помощи нанесения на график двух направленных каротажных измерений в качестве пары координат; и получают расстояние до границы слоя и удельное сопротивление пласта при помощи формирования графика взаимной зависимости для соответствующих значений координатной сетки в точке, расположение которой определено.
21. Способ по п.20, при котором этап формирования графика взаимной зависимости осуществляют с использованием модели с одной границей, при этом полученное удельное сопротивление является удельным сопротивлением на границе раздела «слой - уступ в месте искривления скважины», а полученное расстояние является наименьшим расстоянием до границы раздела «слой - уступ в месте искривления скважины».
22. Способ по п.20, при котором полученные удельное сопротивление и расстояния используются для принятия решений по бурению.
23. Способ определения характеристик подземного пласта, включающий в себя этапы, на которых генерируют график, используя удельное сопротивление и направленное измерение в качестве области значений и области определения графика и имея удельные сопротивления пласта и расстояние до границы слоя, нанесенной на график, формируют вслед за тем координатную сетку;
выбирают удельное сопротивление и направленное измерение, получаемые от прибора, размещенного в стволе скважины, пересекающей пласт;
определяют расположение точки на графике при помощи нанесения на график удельного сопротивления и направленного измерения в качестве пары координат; и получают расстояние до границы слоя и удельное сопротивление пласта при помощи формирования графика взаимной зависимости для соответствующих значений координатной сетки в точке, расположение которой определено.
24. Способ определения характеристик подземного пласта, включающий в себя этапы, на которых генерируют график, используя удельное сопротивление и направленное измерение в качестве области значений и области определения графика и имея удельные сопротивления пласта и расстояние до границы слоя, нанесенной на график, формируют вслед за тем координатную сетку;
выбирают удельное сопротивление и два направленных каротажных измерения, получаемых от прибора, размещенного в стволе скважины, пересекающей пласт;
определяют расположение двух точек на графике при помощи нанесения на график удельного сопротивления и каждого из двух направленных каротажных измерений в качестве пары координат; и получают расстояние до границы слоя и удельное сопротивление пласта для по меньшей мере двух слоев пласта при помощи формирования графика взаимной зависимости для соответствующих значений координатной сетки в точках, расположение которых определено.
25. Способ по п.24, дополнительно включающий в себя этап, на котором выбирают подходящую инверсную модель для выбранных направленных измерений в реальном масштабе времени;
проверяют, согласуется ли выбранная модель с другой информацией; и используют проверенную модель для принятия решений по бурению.
26. Способ по п.25, при котором этап выбора модели включает в себя создание визуализации выбранных направленных измерений.
27. Способ по п.25, при котором этап выбора модели включает в себя идентификацию параметров известной информации;
интерактивный выбор моделей, с помощью которых следует инвертировать выбранные направленные измерения; и выбор простейшей модели, которая аппроксимирует известную информацию.
28. Устройство для измерения характеристик пластов грунта, окружающих ствол скважины, содержащее каротажный прибор, выполненный с возможностью размещения внутри ствола скважины, причем этот каротажный прибор имеет продольную ось и оснащен парами первых и вторых передающих и приемных антенн;
при этом пара первых передающей и приемной антенн содержит
- 23 007587 первую передающую антенну, представляющую собой щелевой симметричный вибратор, ориентированный в первом направлении относительно продольной оси каротажного прибора; и первую приемную антенну, находящуюся на первом расстоянии от первой передающей антенны, представляющую собой щелевой симметричный вибратор, ориентированный во втором направлении, причем первое и второе направления различаются;
при этом щелевые симметричные вибраторы первых передающей и приемной антенн определяют плоскость, в которой лежит продольная ось каротажного прибора;
при этом пара вторых передающей и приемной антенн содержит вторую передающую антенну, представляющую собой щелевой симметричный вибратор, ориентированный во втором направлении относительно продольной оси каротажного прибора; и вторую приемную антенну, находящуюся на первом расстоянии от второй передающей антенны и представляющую собой щелевой симметричный вибратор, ориентированный в первом направлении, причем первое и второе направления являются разными;
причем щелевые симметричные вибраторы вторых передающей и приемной антенн определяют плоскость, в которой лежит продольная ось каротажного прибора, а устройство также содержит датчик угла торца инструмента для непрерывной индикации азимутальной ориентации каротажного прибора; и контроллер для управления парами первых и вторых передающей и приемной антенн, обеспечивающий избирательную передачу электромагнитной энергии в пласт и измерение сигналов напряжения, связь которых с передаваемой электромагнитной энергией представляет собой функцию азимутальной ориентации каротажного прибора.
29. Устройство по п.28, в котором пары вторых передающей и приемной антенн ориентированы под первым азимутальным углом относительно пар первых передающей и приемной антенн вокруг продольной оси каротажного прибора.
30. Устройство по п.28, в котором каждая из передающей и приемной антенн обладает функциональными возможностями приемопередающей антенны.
31. Устройство по п.28, в котором пары первых и вторых передающих и приемных антенн находятся в одних и тех же физических положениях на каротажном приборе.
EA200400573A 2003-05-22 2004-05-21 Устройство и способ определения удельного сопротивления с помощью направленной электромагнитной волны EA007587B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US47268603P 2003-05-22 2003-05-22
US10/709,212 US7382135B2 (en) 2003-05-22 2004-04-21 Directional electromagnetic wave resistivity apparatus and method

Publications (2)

Publication Number Publication Date
EA200400573A1 EA200400573A1 (ru) 2004-12-30
EA007587B1 true EA007587B1 (ru) 2006-12-29

Family

ID=32511770

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200400573A EA007587B1 (ru) 2003-05-22 2004-05-21 Устройство и способ определения удельного сопротивления с помощью направленной электромагнитной волны

Country Status (9)

Country Link
US (1) US7382135B2 (ru)
CN (2) CN1573013B (ru)
BR (1) BRPI0401805B1 (ru)
DE (1) DE102004024969A1 (ru)
EA (1) EA007587B1 (ru)
FR (2) FR2855272B1 (ru)
GB (1) GB2402489B (ru)
MX (1) MXPA04004692A (ru)
NO (1) NO335727B1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2655279C1 (ru) * 2017-06-19 2018-05-24 Публичное акционерное общество "Татнефть" имени В.Д. Шашина Способ определения геомеханических параметров горных пород

Families Citing this family (148)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7659722B2 (en) 1999-01-28 2010-02-09 Halliburton Energy Services, Inc. Method for azimuthal resistivity measurement and bed boundary detection
US6163155A (en) * 1999-01-28 2000-12-19 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations
GB2407170B (en) * 2002-12-31 2005-06-22 Schlumberger Holdings System and method for locating a fracture in an earth formation
US7202670B2 (en) * 2003-08-08 2007-04-10 Schlumberger Technology Corporation Method for characterizing a subsurface formation with a logging instrument disposed in a borehole penetrating the formation
GB2417783B (en) * 2003-08-08 2006-09-20 Schlumberger Holdings Methods and apparatus for characterising a subsurface formation
US7663363B2 (en) * 2004-02-09 2010-02-16 Baker Hughes Incorporated Method and apparatus for high signal-to-noise ratio NMR well logging
US7423426B2 (en) * 2004-02-09 2008-09-09 Baker Hughes Incorporated Selective excitation in earth's magnetic field nuclear magnetic resonance well logging tool
US7239145B2 (en) * 2004-03-29 2007-07-03 Schlumberger Technology Center Subsurface electromagnetic measurements using cross-magnetic dipoles
US7027926B2 (en) * 2004-04-19 2006-04-11 Pathfinder Energy Services, Inc. Enhanced measurement of azimuthal dependence of subterranean parameters
US7848887B2 (en) * 2004-04-21 2010-12-07 Schlumberger Technology Corporation Making directional measurements using a rotating and non-rotating drilling apparatus
US7786733B2 (en) * 2004-07-14 2010-08-31 Schlumberger Technology Corporation Apparatus and system for well placement and reservoir characterization
US8736270B2 (en) 2004-07-14 2014-05-27 Schlumberger Technology Corporation Look ahead logging system
US7755361B2 (en) * 2004-07-14 2010-07-13 Schlumberger Technology Corporation Apparatus and system for well placement and reservoir characterization
US20060017443A1 (en) * 2004-07-23 2006-01-26 Baker Hughes Incorporated Deep reading propagation resistivity tool for determination of distance to a bed boundary with a transition zone
US7630872B2 (en) * 2004-09-16 2009-12-08 Schlumberger Technology Corporation Methods for visualizing distances between wellbore and formation boundaries
CA2584068C (en) * 2004-10-22 2011-05-10 Baker Hughes Incorporated Magnetic measurements while rotating
US7471088B2 (en) * 2004-12-13 2008-12-30 Baker Hughes Incorporated Elimination of the anisotropy effect in LWD azimuthal resistivity tool data
US7536261B2 (en) * 2005-04-22 2009-05-19 Schlumberger Technology Corporation Anti-symmetrized electromagnetic measurements
CN100429529C (zh) * 2006-01-10 2008-10-29 中国石油天然气集团公司 利用全波列、偶极横波测井资料确定气层的方法
US8931335B2 (en) * 2006-04-07 2015-01-13 Baker Hughes Incorporation Processing of multi-component induction data in the presence of borehole abnormalities
AU2007248114B2 (en) * 2006-05-04 2010-12-16 Shell Internationale Research Maatschappij B.V. Method of analyzing a subterranean formation using time dependent transient response signals
CN101467070B (zh) * 2006-06-19 2014-04-09 哈里伯顿能源服务公司 井下管内的天线挖切部
JP5060555B2 (ja) * 2006-07-11 2012-10-31 ハリバートン エナジー サービシーズ,インコーポレーテッド モジュラー式ジオステアリング用ツールアセンブリ
MX2008014830A (es) * 2006-07-12 2009-03-05 Halliburton Energy Serv Inc Metodo y aparato para construir una antena inclinada.
US8126650B2 (en) * 2006-07-25 2012-02-28 Exxonmobil Upstream Research Co. Method for determining physical properties of structures
US8593147B2 (en) 2006-08-08 2013-11-26 Halliburton Energy Services, Inc. Resistivity logging with reduced dip artifacts
US7979211B2 (en) * 2006-08-24 2011-07-12 Exxonmobil Upstream Research Co. Electromagnetic data processing system
WO2008033184A2 (en) * 2006-09-13 2008-03-20 Exxonmobil Upstream Research Company Rapid inversion of electromagnetic reconnaissance survey data
US8466683B2 (en) * 2006-12-14 2013-06-18 Schlumberger Technology Corporation Determining properties of earth formations using the electromagnetic coupling tensor
EP2066866B1 (en) * 2006-12-15 2018-09-12 Halliburton Energy Services, Inc. Antenna coupling component measurement tool having rotating antenna configuration
WO2008115229A1 (en) * 2007-03-16 2008-09-25 Halliburton Energy Services, Inc. Robust inversion systems and methods for azimuthally sensitive resistivity logging tools
US8049508B2 (en) * 2007-03-16 2011-11-01 Baker Hughes Incorporated Method and apparatus for determining formation boundary near the bit for conductive mud
US7751280B2 (en) * 2007-03-27 2010-07-06 Schlumberger Technology Corporation Determining wellbore position within subsurface earth structures and updating models of such structures using azimuthal formation measurements
US8497685B2 (en) 2007-05-22 2013-07-30 Schlumberger Technology Corporation Angular position sensor for a downhole tool
US7657377B2 (en) 2007-05-31 2010-02-02 Cbg Corporation Azimuthal measurement-while-drilling (MWD) tool
US7558675B2 (en) * 2007-07-25 2009-07-07 Smith International, Inc. Probablistic imaging with azimuthally sensitive MWD/LWD sensors
US7912648B2 (en) * 2007-10-02 2011-03-22 Baker Hughes Incorporated Method and apparatus for imaging bed boundaries using azimuthal propagation resistivity measurements
US8932879B2 (en) * 2007-11-06 2015-01-13 Ambergen, Inc. Methods and compounds for phototransfer
US8906700B2 (en) 2007-11-06 2014-12-09 Ambergen, Inc. Methods and compositions for phototransfer
EP2223161B1 (en) * 2007-11-30 2019-01-16 Exxonmobil Upstream Research Company Method for reservoir fracture and cross beds detection using tri- axial/multi-component resistivity anisotropy measurements
CA2703588C (en) 2007-12-12 2015-12-01 Exxonmobil Upstream Research Company Method and apparatus for evaluating submarine formations
US20090157361A1 (en) * 2007-12-12 2009-06-18 Toghi Farid Method of well placement modeling and geosteering
JP5662804B2 (ja) * 2007-12-18 2015-02-04 シュルンベルジェ ホールディングス リミテッドSchlnmberger Holdings Limited 表面電磁探査を改善するシステム及び方法
GB2468734B (en) 2008-01-18 2012-08-08 Halliburton Energy Serv Inc Em-guided drilling relative to an existing borehole
CN101525999B (zh) * 2008-03-06 2013-04-24 中国石油化工股份有限公司 电磁随钻测量系统的适应性分析方法
US8278931B2 (en) * 2008-07-14 2012-10-02 Baker Hughes Incorporated Anisotropy orientation image from resistivity measurements for geosteering and formation evaluation
US7991555B2 (en) * 2008-07-30 2011-08-02 Schlumberger Technology Corporation Electromagnetic directional measurements for non-parallel bed formations
US8278928B2 (en) * 2008-08-25 2012-10-02 Baker Hughes Incorporated Apparatus and method for detection of position of a component in an earth formation
US8427162B2 (en) * 2008-08-25 2013-04-23 Baker Hughes Incorporated Apparatus and method for detection of position of a component in an earth formation
WO2010039357A2 (en) * 2008-10-01 2010-04-08 Schlumberger Canada Limited Logging tool with antennas having equal tilt angles
WO2010074678A2 (en) 2008-12-16 2010-07-01 Halliburton Energy Services, Inc. Azimuthal at-bit resistivity and geosteering methods and systems
US8089268B2 (en) * 2009-03-24 2012-01-03 Smith International, Inc. Apparatus and method for removing anisotropy effect from directional resistivity measurements
US8195400B2 (en) * 2009-05-08 2012-06-05 Smith International, Inc. Directional resistivity imaging using harmonic representations
US8159227B2 (en) * 2009-05-11 2012-04-17 Smith International Inc. Methods for making directional resistivity measurements
US7990153B2 (en) * 2009-05-11 2011-08-02 Smith International, Inc. Compensated directional resistivity measurements
US8497673B2 (en) * 2009-09-28 2013-07-30 Schlumberger Technology Corporation Directional resistivity antenna shield
US8466682B2 (en) * 2009-09-29 2013-06-18 Schlumberger Technology Corporation Apparatus and method for downhole electromagnetic measurement while drilling
US8433518B2 (en) * 2009-10-05 2013-04-30 Schlumberger Technology Corporation Multilevel workflow method to extract resistivity anisotropy data from 3D induction measurements
US8271199B2 (en) * 2009-12-31 2012-09-18 Smith International, Inc. Binning method for borehole imaging
BRPI1013305B1 (pt) 2010-01-22 2019-09-10 Halliburton Energy Services Inc sistema para medir uma resistividade de uma formação, método para determinar uma resistividade,e, broca de perfuração instrumentada
WO2011129828A1 (en) * 2010-04-15 2011-10-20 Halliburton Energy Services, Inc. Processing and geosteering with a rotating tool
CN103250072B (zh) * 2010-04-27 2017-02-22 哈里伯顿能源服务公司 通过使用钻头地震对断裂的干涉测量钻头成像、时间反演成像来进行断裂表征、和经由时间反演的声学和电震对断裂产生的监测
BR112012029525A2 (pt) * 2010-05-21 2016-12-06 Landmark Graphics Corp dispositivos e método para a correlação e geodirecionamento de poços horizontais
US8600115B2 (en) 2010-06-10 2013-12-03 Schlumberger Technology Corporation Borehole image reconstruction using inversion and tool spatial sensitivity functions
US9372276B2 (en) 2010-06-10 2016-06-21 Schlumberger Technology Corporation Combinations of axial and saddle coils to create the equivalent of tilted coils for directional resistivity measurements
US8749243B2 (en) * 2010-06-22 2014-06-10 Halliburton Energy Services, Inc. Real time determination of casing location and distance with tilted antenna measurement
US9933541B2 (en) * 2010-06-22 2018-04-03 Schlumberger Technology Corporation Determining resistivity anisotropy and formation structure for vertical wellbore sections
US9115569B2 (en) 2010-06-22 2015-08-25 Halliburton Energy Services, Inc. Real-time casing detection using tilted and crossed antenna measurement
US8844648B2 (en) 2010-06-22 2014-09-30 Halliburton Energy Services, Inc. System and method for EM ranging in oil-based mud
US8756018B2 (en) 2010-06-23 2014-06-17 Schlumberger Technology Corporation Method for time lapsed reservoir monitoring using azimuthally sensitive resistivity measurements while drilling
US8558548B2 (en) 2010-07-28 2013-10-15 Schlumberger Technology Corporation Determining anisotropic resistivity
US9360582B2 (en) 2010-07-02 2016-06-07 Halliburton Energy Services, Inc. Correcting for magnetic interference in azimuthal tool measurements
WO2012015542A1 (en) 2010-07-27 2012-02-02 Exxonmobil Upstream Research Company Inverting geophysical data for geological parameters or lithology
CA2806874C (en) 2010-08-16 2016-12-06 Exxonmobil Upstream Research Company Reducing the dimensionality of the joint inversion problem
US9273517B2 (en) 2010-08-19 2016-03-01 Schlumberger Technology Corporation Downhole closed-loop geosteering methodology
BR112013007048A2 (pt) * 2010-09-17 2016-06-14 Baker Hugues Inc aparelho e método de perfuração de poços determinando furos de sondagem existentes usando dispositivos de indução
CN101982734A (zh) * 2010-10-29 2011-03-02 南昌航空大学 一种用于地下磁导航的计算方法
US8536871B2 (en) 2010-11-02 2013-09-17 Schlumberger Technology Corporation Method of correcting resistivity measurements for toll bending effects
US9658360B2 (en) 2010-12-03 2017-05-23 Schlumberger Technology Corporation High resolution LWD imaging
US9377555B2 (en) * 2010-12-23 2016-06-28 China Petroleum & Chemical Corporation Apparatus and method for well logging and data processing device
US9103936B2 (en) 2011-01-20 2015-08-11 Halliburton Energy Services, Inc. Method and apparatus for determining and displaying geosteering information
US8626446B2 (en) 2011-04-01 2014-01-07 Schlumberger Technology Corporation Method of directional resistivity logging
WO2012144981A1 (en) * 2011-04-18 2012-10-26 Halliburton Energy Services, Inc. Method for real-time downhole processing and detection of bed boundary for geosteering application
US9453929B2 (en) 2011-06-02 2016-09-27 Exxonmobil Upstream Research Company Joint inversion with unknown lithology
EP2721478A4 (en) 2011-06-17 2015-12-02 Exxonmobil Upstream Res Co FREEZING OF DOMAINS IN A CONNECTION VERSION
US9494711B2 (en) 2011-07-21 2016-11-15 Garrett M Leahy Adaptive weighting of geophysical data types in joint inversion
CN103717833B (zh) * 2011-07-29 2016-10-19 兰德马克绘图国际公司 将测量测井与预测测井相关的方法和系统
US10539009B2 (en) * 2011-08-10 2020-01-21 Scientific Drilling International, Inc. Short range data transmission in a borehole
CA2844111C (en) 2011-08-18 2016-11-08 Halliburton Energy Services, Inc. Improved casing detection tools and methods
US9239403B2 (en) * 2011-08-29 2016-01-19 Hallibburton Energy Services, Inc. Apparatus and methods of controlling recordation of resistivity-related readings in determining formation resistivity
US10823871B2 (en) * 2011-11-15 2020-11-03 Halliburton Energy Services, Inc. Enhanced resistivity measurement with at-bit resistivity sensor
US10371852B2 (en) * 2011-12-21 2019-08-06 Schlumberger Technology Corporation Formation properties from conductivity tensor
WO2013106318A1 (en) 2012-01-09 2013-07-18 Zimmer, Inc. Composite device that combines porous metal and bone stimuli
US9540922B2 (en) 2012-03-29 2017-01-10 Schlumberger Technology Corporation Electromagnetic method for obtaining dip azimuth angle
WO2014003701A1 (en) 2012-06-25 2014-01-03 Halliburton Energy Services, Inc. Resistivity logging systems and methods employing ratio signal set for inversion
WO2014003702A1 (en) 2012-06-25 2014-01-03 Halliburton Energy Services, Inc. Tilted antenna logging systems and methods yielding robust measurement signals
RU2598003C1 (ru) * 2012-08-10 2016-09-20 Хэллибертон Энерджи Сервисиз, Инк. Способы и системы для прямого моделирования скважинного изображения свойств пласта
US9182510B2 (en) 2012-08-20 2015-11-10 Landmark Graphics Corporation Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models
NO2895950T3 (ru) * 2012-11-13 2018-07-07
CN103015995B (zh) * 2012-12-14 2015-06-24 中国电子科技集团公司第二十二研究所 通过带倾斜线圈的电阻率设备的地层倾角测量方法及装置
BR112015012050A2 (pt) 2012-12-31 2019-12-17 Halliburton Energy Services Inc imageamento de formação com antenas multipolo
RU2628000C2 (ru) * 2012-12-31 2017-08-14 Халлибертон Энерджи Сервисез, Инк. Глубокая азимутальная система с применением многополюсных датчиков
US10591638B2 (en) 2013-03-06 2020-03-17 Exxonmobil Upstream Research Company Inversion of geophysical data on computer system having parallel processors
US10324219B2 (en) * 2013-03-15 2019-06-18 Halliburton Energy Services, Inc. Identifying unconventional formations
US11480705B2 (en) * 2013-04-01 2022-10-25 Oliden Technology, Llc Antenna, tool, and methods for directional electromagnetic well logging
US10345475B2 (en) 2013-04-02 2019-07-09 Schlumberger Technology Corporation Extended 1D inversion of electromagnetic measurements for formation evaluation
US9846255B2 (en) 2013-04-22 2017-12-19 Exxonmobil Upstream Research Company Reverse semi-airborne electromagnetic prospecting
CA2915348C (en) 2013-06-12 2023-05-02 Well Resolutions Technology Apparatus and methods for making azimuthal resistivity measurements
US9207053B2 (en) * 2013-06-21 2015-12-08 Rosemount Aerospace Inc. Harmonic shuttered seeker
WO2015050985A1 (en) 2013-10-01 2015-04-09 Schlumberger Canada Limited Methods, computer-readable media, and systems for applying 1-dimensional (1d) processing in a non-1d formation
US10190408B2 (en) * 2013-11-22 2019-01-29 Aps Technology, Inc. System, apparatus, and method for drilling
CN103670387A (zh) * 2013-12-12 2014-03-26 贝兹维仪器(苏州)有限公司 一种地层定向电阻率测量方法及装置
CN103711474B (zh) * 2013-12-19 2016-08-17 天津大学 一种正交偶极子声电组合测井仪器
US9765613B2 (en) 2014-03-03 2017-09-19 Aps Technology, Inc. Drilling system and electromagnetic telemetry tool with an electrical connector assembly and associated methods
WO2015137921A1 (en) 2014-03-11 2015-09-17 Halliburton Energy Services, Inc. Multi-component induction logging systems and methods using blended-model inversion
US10215878B2 (en) 2014-03-29 2019-02-26 Schlumberger Technology Corporation Gain compensated directional propagation measurements
US9581721B2 (en) 2014-03-29 2017-02-28 Schlumberger Technology Corporation Method for making downhole electromagnetic logging while drilling measurements
EP3126626A4 (en) * 2014-03-30 2017-11-15 Services Pétroliers Schlumberger Gain compensated measurements using tilted antennas
MX369871B (es) 2014-04-03 2019-11-25 Halliburton Energy Services Inc Metodos y sistemas de registro de induccion con componentes multiples utilizando inversion de frecuencia seleccionada.
US9790784B2 (en) 2014-05-20 2017-10-17 Aps Technology, Inc. Telemetry system, current sensor, and related methods for a drilling system
CN105464650A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 一种随钻测井解释方法
US10401529B2 (en) 2014-10-17 2019-09-03 Halliburton Energy Services, Inc. Fast-changing dip formation resistivity estimation
US10267945B2 (en) * 2014-10-20 2019-04-23 Schlumberger Technology Corporation Use of transverse antenna measurements for casing and pipe detection
CA2960154C (en) * 2014-11-13 2019-06-11 Halliburton Energy Services, Inc. Resistivity logging tools with tilted ferrite elements for azimuthal sensitivity
US9784880B2 (en) * 2014-11-20 2017-10-10 Schlumberger Technology Corporation Compensated deep propagation measurements with differential rotation
DE112015005897T5 (de) * 2015-01-07 2017-10-12 Halliburton Energy Services, Inc. Funktionale Erdmodellparametrierung zur Widerstandsinvertierung
US9976413B2 (en) 2015-02-20 2018-05-22 Aps Technology, Inc. Pressure locking device for downhole tools
CN104732562A (zh) * 2015-03-24 2015-06-24 中国石油集团川庆钻探工程有限公司 一种电成像测井图像相似度比对方法
US11307323B2 (en) * 2015-09-30 2022-04-19 Schlumberger Technology Corporation Methods and systems to analyze bed boundary detection
WO2017132098A1 (en) * 2016-01-25 2017-08-03 Schlumberger Technology Corporation Drilling with information characterizing lateral heterogeneities based on deep directional resistivity measurements
CA3005858C (en) * 2016-02-16 2021-07-27 Halliburton Energy Services, Inc. Generating an earth model from spatial correlations of equivalent earth models
CN106324689B (zh) * 2016-06-24 2018-05-11 杭州迅美科技有限公司 一种水平井地层环境下电阻率各向异性识别方法
CN107762494A (zh) * 2016-08-19 2018-03-06 中国石油集团长城钻探工程有限公司 随钻偶极声波测井仪器以及声波测井方法
WO2018144029A1 (en) * 2017-02-06 2018-08-09 Halliburton Energy Services, Inc. Multi-layer distance to bed boundary (dtbb) inversion with multiple initial guesses
US9995840B1 (en) * 2017-04-17 2018-06-12 Nabors Drilling Technologies Usa, Inc. Azimuthal minor averaging in a wellbore
US11391859B2 (en) 2018-06-29 2022-07-19 Halliburton Energy Services, Inc. Determining formation properties in a geological formation using an inversion process on a modified response matrix associated with a downhole tool
US11015432B2 (en) * 2018-07-10 2021-05-25 Halliburton Energy Services, Inc. Relative azimuth correction for resistivity inversion
CN109668938B (zh) * 2018-12-13 2024-01-09 长江水利委员会长江科学院 堤防渗漏通道三维磁测电阻率探测装置及方法
WO2020149824A1 (en) * 2019-01-14 2020-07-23 Halliburton Energy Services, Inc. Azimuthal borehole rendering of radioelement spectral gamma data
WO2020159484A1 (en) * 2019-01-29 2020-08-06 Halliburton Energy Services, Inc. Iterative borehole correction
CN113188437B (zh) * 2019-07-17 2022-12-23 伟志股份公司 一种高精度城市建筑基坑常规变形监测系统及其监测方法
EP4111240A4 (en) * 2020-02-28 2024-03-27 Novamera Inc GROUND penetrating radar device and method
US20210285297A1 (en) * 2020-03-13 2021-09-16 Baker Hughes Oilfield Operations, Llc Automated geosteering based on a distance to oil-water contact
CN111827986B (zh) * 2020-09-03 2023-10-10 北京贝威通能源科技集团有限公司 一种近钻头无线短传系统和方法
CN112963142B (zh) * 2021-03-26 2024-02-02 北京吉星恒大能源科技有限公司 一种双制式无线传输高低电阻率测量近钻头系统
CN113359199B (zh) * 2021-07-06 2022-07-08 北京航空航天大学 一种基于聚焦磁场的井周电阻率测量方法
CN113885085B (zh) * 2021-08-05 2023-06-09 中煤科工集团西安研究院有限公司 一种井下直流轴向偶极动源超前探测方法
CN114526063A (zh) * 2022-02-14 2022-05-24 北京工业大学 获取探边电磁波测井仪器结构参数的方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU998995A1 (ru) * 1981-09-11 1983-02-23 Всесоюзный научно-исследовательский институт нефтепромысловой геофизики Устройство дл электромагнитного каротажа
US4766384A (en) * 1986-06-20 1988-08-23 Schlumberger Technology Corp. Well logging apparatus for determining dip, azimuth, and invaded zone conductivity
US5406206A (en) * 1991-05-28 1995-04-11 Schlumberger Technology Corporation Method of evaluating a geological formation using a logging tool including slot antenna having two nonparallel elements
RU2069878C1 (ru) * 1992-02-10 1996-11-27 Кузьмичев Олег Борисович Способ электромагнитного каротажа скважин
RU2107313C1 (ru) * 1996-07-12 1998-03-20 Дворецкий Петр Иванович Способ геофизических исследований скважин сложной конфигурации, основанный на применении направленных широкополосных электромагнитных импульсов, возбуждаемых щелевой цилиндрической антенной решеткой
US6297639B1 (en) * 1999-12-01 2001-10-02 Schlumberger Technology Corporation Method and apparatus for directional well logging with a shield having sloped slots

Family Cites Families (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4319191A (en) 1980-01-10 1982-03-09 Texaco Inc. Dielectric well logging with radially oriented coils
US4689569A (en) * 1984-12-17 1987-08-25 Southwest Research Institute Directional antenna system for use in a borehole incorporating antenna dipole elements
US5115198A (en) 1989-09-14 1992-05-19 Halliburton Logging Services, Inc. Pulsed electromagnetic dipmeter method and apparatus employing coils with finite spacing
US4980643A (en) * 1989-09-28 1990-12-25 Halliburton Logging Services, Inc. Induction logging and apparatus utilizing skew signal measurements in dipping beds
US5235285A (en) 1991-10-31 1993-08-10 Schlumberger Technology Corporation Well logging apparatus having toroidal induction antenna for measuring, while drilling, resistivity of earth formations
JP2534193B2 (ja) 1993-05-31 1996-09-11 石油資源開発株式会社 指向性インダクション検層法および装置
EP0665958B1 (en) * 1993-07-21 1999-01-13 Western Atlas International, Inc. Method of determining formation resistivity utilizing combined measurements of inductive and galvanic logging instruments
US5517464A (en) 1994-05-04 1996-05-14 Schlumberger Technology Corporation Integrated modulator and turbine-generator for a measurement while drilling tool
US5757191A (en) 1994-12-09 1998-05-26 Halliburton Energy Services, Inc. Virtual induction sonde for steering transmitted and received signals
CA2256771C (en) 1996-07-01 2005-04-26 David R. Beard Electrical logging of a laminated earth formation
US5781436A (en) 1996-07-26 1998-07-14 Western Atlas International, Inc. Method and apparatus for transverse electromagnetic induction well logging
US5992228A (en) * 1996-10-23 1999-11-30 Dunham; Lanny L. Method for determining resistivity derived porosity and porosity derived resistivity
US6044325A (en) 1998-03-17 2000-03-28 Western Atlas International, Inc. Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument
US6068957A (en) 1998-06-29 2000-05-30 Eastman Kodak Company Lubricating layer in photographic elements
US6476609B1 (en) 1999-01-28 2002-11-05 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone
US6163155A (en) 1999-01-28 2000-12-19 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations
US6181138B1 (en) 1999-02-22 2001-01-30 Halliburton Energy Services, Inc. Directional resistivity measurements for azimuthal proximity detection of bed boundaries
US6442488B2 (en) * 1999-03-08 2002-08-27 Baker Hughes Incorporated Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well
US6304086B1 (en) 1999-09-07 2001-10-16 Schlumberger Technology Corporation Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers
US6594584B1 (en) 1999-10-21 2003-07-15 Schlumberger Technology Corporation Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument
US6351127B1 (en) 1999-12-01 2002-02-26 Schlumberger Technology Corporation Shielding method and apparatus for selective attenuation of an electromagnetic energy field component
US6566881B2 (en) 1999-12-01 2003-05-20 Schlumberger Technology Corporation Shielding method and apparatus using transverse slots
US6727705B2 (en) 2000-03-27 2004-04-27 Schlumberger Technology Corporation Subsurface monitoring and borehole placement using a modified tubular equipped with tilted or transverse magnetic dipoles
US6836218B2 (en) 2000-05-22 2004-12-28 Schlumberger Technology Corporation Modified tubular equipped with a tilted or transverse magnetic dipole for downhole logging
US6393363B1 (en) * 2000-06-28 2002-05-21 Schlumberger Technology Corporation Method and apparatus for the measurement of the electrical resistivity of geologic formations employing modeling data
US6509738B1 (en) * 2000-07-14 2003-01-21 Schlumberger Technology Corporation Electromagnetic induction well logging instrument having azimuthally sensitive response
US6573722B2 (en) 2000-12-15 2003-06-03 Schlumberger Technology Corporation Method and apparatus for cancellation of borehole effects due to a tilted or transverse magnetic dipole
US6541979B2 (en) 2000-12-19 2003-04-01 Schlumberger Technology Corporation Multi-coil electromagnetic focusing methods and apparatus to reduce borehole eccentricity effects
US6636045B2 (en) * 2001-04-03 2003-10-21 Baker Hughes Incorporated Method of determining formation anisotropy in deviated wells using separation of induction mode
US6584408B2 (en) 2001-06-26 2003-06-24 Schlumberger Technology Corporation Subsurface formation parameters from tri-axial measurements
US6969994B2 (en) 2001-09-26 2005-11-29 Schlumberger Technology Corporation Directional electromagnetic measurements insensitive to dip and anisotropy
US6556015B1 (en) 2001-10-11 2003-04-29 Schlumberger Technology Corporation Method and system for determining formation anisotropic resistivity with reduced borehole effects from tilted or transverse magnetic dipoles
US6819110B2 (en) 2002-03-26 2004-11-16 Schlumberger Technology Corporation Electromagnetic resistivity logging instrument with transverse magnetic dipole component antennas providing axially extended response
US6690170B2 (en) 2002-03-29 2004-02-10 Schlumberger Technology Corporation Antenna structures for electromagnetic well logging tools
US6667620B2 (en) 2002-03-29 2003-12-23 Schlumberger Technology Corporation Current-directing shield apparatus for use with transverse magnetic dipole antennas
US6930652B2 (en) 2002-03-29 2005-08-16 Schlumberger Technology Corporation Simplified antenna structures for logging tools
US6998844B2 (en) 2002-04-19 2006-02-14 Schlumberger Technology Corporation Propagation based electromagnetic measurement of anisotropy using transverse or tilted magnetic dipoles
US6794875B2 (en) 2002-05-20 2004-09-21 Halliburton Energy Services, Inc. Induction well logging apparatus and method
US6937022B2 (en) 2002-09-06 2005-08-30 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
US6903553B2 (en) 2002-09-06 2005-06-07 Baker Hughes Incorporated Method and apparatus for a quadrupole transmitter for directionally sensitive induction tool
US6819111B2 (en) * 2002-11-22 2004-11-16 Baker Hughes Incorporated Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electro-magnetic antenna combination and orientation with additional rotation and position measurements
US6937021B2 (en) * 2002-12-09 2005-08-30 Schlumberger Technology Corporation Method and apparatus for determining the presence and orientation of a fraction in an earth formation
US6924646B2 (en) * 2002-12-31 2005-08-02 Schlumberger Technology Corporation System and method for locating a fracture in an earth formation
US7202670B2 (en) 2003-08-08 2007-04-10 Schlumberger Technology Corporation Method for characterizing a subsurface formation with a logging instrument disposed in a borehole penetrating the formation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU998995A1 (ru) * 1981-09-11 1983-02-23 Всесоюзный научно-исследовательский институт нефтепромысловой геофизики Устройство дл электромагнитного каротажа
US4766384A (en) * 1986-06-20 1988-08-23 Schlumberger Technology Corp. Well logging apparatus for determining dip, azimuth, and invaded zone conductivity
US5406206A (en) * 1991-05-28 1995-04-11 Schlumberger Technology Corporation Method of evaluating a geological formation using a logging tool including slot antenna having two nonparallel elements
RU2069878C1 (ru) * 1992-02-10 1996-11-27 Кузьмичев Олег Борисович Способ электромагнитного каротажа скважин
RU2107313C1 (ru) * 1996-07-12 1998-03-20 Дворецкий Петр Иванович Способ геофизических исследований скважин сложной конфигурации, основанный на применении направленных широкополосных электромагнитных импульсов, возбуждаемых щелевой цилиндрической антенной решеткой
US6297639B1 (en) * 1999-12-01 2001-10-02 Schlumberger Technology Corporation Method and apparatus for directional well logging with a shield having sloped slots

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2655279C1 (ru) * 2017-06-19 2018-05-24 Публичное акционерное общество "Татнефть" имени В.Д. Шашина Способ определения геомеханических параметров горных пород

Also Published As

Publication number Publication date
MXPA04004692A (es) 2004-11-25
EA200400573A1 (ru) 2004-12-30
NO335727B1 (no) 2015-02-02
GB0410076D0 (en) 2004-06-09
US20050140373A1 (en) 2005-06-30
FR2855272A1 (fr) 2004-11-26
CN1573013A (zh) 2005-02-02
CN101082276B (zh) 2012-03-28
FR2859539B1 (fr) 2013-01-18
BRPI0401805A (pt) 2005-01-18
FR2855272B1 (fr) 2013-02-08
FR2859539A1 (fr) 2005-03-11
BRPI0401805B1 (pt) 2016-07-26
CN101082276A (zh) 2007-12-05
DE102004024969A1 (de) 2006-01-05
GB2402489B (en) 2005-11-30
US7382135B2 (en) 2008-06-03
NO20042106L (no) 2004-11-23
CN1573013B (zh) 2010-09-01
GB2402489A (en) 2004-12-08

Similar Documents

Publication Publication Date Title
EA007587B1 (ru) Устройство и способ определения удельного сопротивления с помощью направленной электромагнитной волны
US11119239B2 (en) Measuring petrophysical properties of an earth formation by regularized direct inversion of electromagnetic signals
US10451765B2 (en) Post-well reservoir characterization using image-constrained inversion
US9547102B2 (en) Resistivity logging systems and methods employing ratio signal set for inversion
US8364404B2 (en) System and method for displaying data associated with subsurface reservoirs
US9678240B2 (en) Detecting boundary locations of multiple subsurface layers
US10358911B2 (en) Tilted antenna logging systems and methods yielding robust measurement signals
US9110188B2 (en) Detecting bed boundary locations based on gradients determined from measurements from multiple tool depths in a wellbore
US20050083063A1 (en) Electromagnetic method for determining dip angles independent of mud type and borehole environment
EP3410160A1 (en) Method for real-time downhole processing and detection of bed boundary for geosteering application
US10330818B2 (en) Multi-component induction logging systems and methods using real-time OBM borehole correction
CN105074505A (zh) 真岩层电阻率的确定
BR112018017328B1 (pt) Método e sistema para caracterização de uma formação subterrânea
BR112018008755B1 (pt) Método e sistema para caracterização de uma formação subterrânea
US20210373190A1 (en) Evaluation and visualization of azimuthal resistivity data
GB2417328A (en) Methods of characterising earth formations
BR112019014132B1 (pt) Método de avaliação de uma formação de terra atravessada por um poço inacabado
BR112018072967B1 (pt) Método para caracterizar uma formação geológica atravessada por um poço de exploração
GB2417783A (en) Method for characterising a subsurface formation

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 BY KG MD TJ

QB4A Registration of a licence in a contracting state
QB4A Registration of a licence in a contracting state
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