Способы интерпретации диффузионных-т2 карт, полученных с использованием ямр данных

Способы интерпретации диффузионных-т2 карт, полученных с использованием ямр данных

Иллюстрации

Показать все

Реферат

Область изобретения

Как вода, так и углеводороды в земных пластах формируют детектируемые ЯМР (ядерно-магнитный резонанс) сигналы. Желательно, чтобы сигналы от воды и от углеводородов были разделяемыми так, чтобы можно было идентифицировать зоны, содержащие углеводороды. Однако не всегда легко отличить, какие сигналы получены от воды, а какие от углеводородов. Для идентифицирования сигналов от воды и от углеводородов были предложены различные способы.

Способы дифференциального спектра (DSM, СДС) и смещенного спектра (SSM, ССС), предложенные Akkurt и соавт. в публикации “NMR Logging of Natural Gas Reservoirs” Paper N. Transactions of the Society of Professional Well Log Analysts (SPWLA) Annual Logging Symposium, 1995, сравнивают распределения T₂, выведенные из двух измерений по Carr-Purcell-Meiboom-Gill (CPMG), выполненных с различными временами поляризации (DSM, СДС) или разнесениями эхо-сигналов (SSM, ССС). Модификация указанных способов, известная как анализ временных выборок (TDA, АВВ), позже была введена Prammer и соавт. в публикации “Lithology-Independent Gas Detection by Gradient-NMR Logging,” SPE paper 30562, 1995. В способе АВВ "различные" данные вычисляются непосредственно во временном представлении путем вычитания одного набора измеренных амплитуд из другого. Затем предполагается, что разностный набор данных охватывает только легкую нефть и/или газ. В способе АВВ относительные вклады от легкой нефти или газа выводятся путем выполнения линейного анализа способом наименьших квадратов разностных данных, используя предполагаемые ЯМР отклики от указанных флюидов. Как в способе СДС, так и в способе АВВ предполагается, что сигнал от воды имеет время T₁ релаксации, существенно более короткое, чем таковое от углеводородов. Однако данное предположение оказывается не всегда справедливым. В особенности указанное предположение не оправдывается в пластах, в которых имеются большие поры или в которых углеводороды обладают средней или высокой вязкостью. ССС способ или его преемник, способ повышенной диффузии (EDM, СПД), предложенный Akkurt и соавт. в публикации “Enhanced Diffusion: Expanding the Range of NMR Direct Hydrocarbon Typing Applications”, Paper GG. Transactions of the Society of Professional Well Log Analysts (SPWLA) Annual Logging Symposium, 1998, разделяет вклады от газа, нефти и воды на основе изменений в распределениях T₂, которые происходят в результате изменений в разнесении эхо-сигналов измерений CPMG. Способы применимы в ограниченном диапазоне сред, и точность результата значительно страдает от неполного разделения сигналов воды и углеводородов в выборке T₂. Более того, указанные способы рассчитаны таким образом, чтобы функционировать с CPMG последовательностями. Однако, с диффузионными способами, CPMG последовательности импульсов обеспечивают плохие отношения сигнал/шум вследствие уменьшенного числа эхо-сигналов, которые могут быть измерены. Стратегия комбинирования и выбора указанных различных ЯМР способов недавно описана Coates и соавт. в публикации патента США № 6366087.

Последовательность импульсов с корректированной диффузией (DE, КД), предложенная Hiirlimann, обеспечивает другой подход. См.: Hiirlimann и соавт., “Diffusion-Editing: New NMR Measurement of Saturation and Pore Geometry,” paper presented at the 2002 Annual Meeting of the Society of Professional Well Log Analysts, Osio, Japan, June 2-5; см. также патент США № 6570382, зарегистрированный 28 ноября 2000 Hiirlimann. Данный патент передан тому же приемнику патентовладельца, что и настоящее изобретение, и тем самым включен здесь ссылкой.

КД последовательности импульсов подобны CPMG последовательностям, за исключением того, что воспринимаются два исходных эхо-сигнала с более длинными разнесениями эхо-сигналов, а третий и последующие эхо-сигналы воспринимаются с более короткими разнесениями эхо-сигналов. В КД последовательностях импульсов, информация диффузии кодируется во время восприятия первых двух эхо-сигналов, тогда как третий и последующие эхо-сигналы обеспечивают объемную и поверхностную информацию времен релаксации с относительно малым затуханием сигнала, обусловленным диффузией. Использование известной CPMG последовательности для кодирования информации диффузии требует длительного разнесения между эхо-сигналами, что приводит к недостаточной объемной и поверхностной информации времени релаксации, вследствие того, что диффузионный спад ослабляет сигнал после незначительного количества отражений. Впоследствии серия данных, воспринятых с помощью КД последовательностей, обеспечивает лучшую информацию диффузии и лучшее отношение сигнал/шум в данных спинового эхо по сравнению с аналогичной серией, воспринятой с помощью CPMG последовательностей. Следовательно, КД последовательности могут обеспечить более точные и надежные вычисления T₂ распределений рапы и нефти, чем CPMG последовательности.

Для улучшения детектирования углеводородов в дополнение к КД последовательностям были разработаны специализированные способы интерпретации для ЯМР данных. Указанные способы обычно применяют прямое моделирование к сериям ЯМР данных, воспринятым с различными параметрами. Серия ЯМР данных обычно воспринимается с различным временем разнесения эхо-сигналов (TE, ВРЭС) или временем поляризации (WT, ВП), а иногда воспринимается с различными градиентами (G) магнитного поля. КД последовательности являются примером такого сбора данных. Два иллюстративных способа включают: способ MACNMR, предложенный Slijkerman и соавт. в публикации SPE paper 56768, “Processing of Multi-Acquisition NMR Data”, 1999, и способ описания характера флюида с помощью магнитного резонанса (MRF, ОХФМР), раскрытый в патенте США № 6229308 на имя Freedman и переданный приемнику патентовладельца настоящего изобретения ("патент Freedman"). Патент Freedman тем самым включен здесь ссылкой.

ОХФМР способ (способ описания характера флюидов магнитным резонансом) способен получать раздельные распределения T₂ нефти и воды. Данный способ использует модель составных вязкостей (CVM, МСВ), которая соотносит время релаксации и скорости диффузии с составными вязкостями, геометрическое среднее которых идентично макроскопической вязкости флюида. При использовании ОХФМР способа оценки для объемов воды и углеводородов получают посредством применения прямой модели для того, чтобы смоделировать ЯМР отклики серии ЯМР измерений, воспринятых с различными параметрами. В частности, ОХФМР метод основан на установленных физических законах, которые, для учета ЯМР откликов флюидов нисходящей скважины, калибруются эмпирически. Используя реалистические модели флюидов, ОХФМР стремится минимизировать число регулируемых параметров, чтобы они были сравнимыми с информационным содержимым данных ЯМР каротажа. Поскольку параметры модели по проекту соотносятся с объемами и свойствами индивидуальных флюидов, определение значений параметров (то есть, аппроксимация данных) ведет непосредственно к оценкам искомых нефтефизических величин.

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

Другой подход, разработанный Schlumberger, основанный на принципе максимальной энтропии (MEP, ПМЭ), состоит из общего, не зависящего от модели способа для анализа комплексных данных флюидов, воспринятых с помощью инструментов ЯМР каротажа, и представляет результаты в визуально привлекательном и легком для понимания формате, который называется здесь как диффузионно-релаксационные карты, или D-T2 карты. Указанные карты уже использовали для понимания случаев, в которых модельный анализ дает неудовлетворительные результаты вследствие отклонений ЯМР свойств от "идеального" поведения, предполагаемого в моделях. Указанные ситуации могут возникнуть из-за аномальных взаимодействий флюид/порода, таких как ограниченная диффузия, смешанная смачиваемость и внутренние градиенты. Отклонения от предполагаемых свойств также наблюдались для некоторых сырых нефтепродуктов, приводя к неточным прогнозам в анализе модели. Посредством использования D-T2 карт, ПМЭ подход обеспечивает простое графическое представление данных, которые могут быть использованы для идентификации откликов флюидов во всех средах. Диффузионно-релаксационные карты описываются в патентах США №№ 6570382 и 6462542.

Хотя указанные способы, известные из уровня техники, полезны в прогнозировании наличия углеводородов в пластах, желательно иметь более простые способы, которые на основе ЯМР данных могут прогнозировать присутствие углеводородов в пластах и которые, вообще говоря, применимы к ЯМР данным, воспринятым с различными последовательностями импульсов. Более того, хотя двух- и трехмерная визуализация была разработана для получения, в первую очередь, качественной информации, желательно иметь количественные методы интерпретации, которые могут обеспечивать точные результаты описания характера флюидов.

Сущность изобретения

Согласно одному аспекту раскрытого объекта изобретения, описан способ для интерпретации многомерных данных ядерно-магнитного резонанса, полученных на образце земного пласта. В частности, собирается набор ЯМР данных для образца флюида, расположенного либо в стволе скважины, либо в лабораторных условиях. Из набора ЯМР данных вычисляется многомерное распределение с использованием математической инверсии, которая является независящей от свойств образца флюида, известных до проведения исследований. Многомерное распределение графически отображается на многомерной карте. Каждое проявление флюида или артефакт, видимый на графике, идентифицируется как представление вероятного существования детектируемого флюида. Одно или несколько решений количественных оценок пласта для одного или нескольких проявлений флюидов вычисляются на основе многомерного распределения, ассоциированного с соответственным проявлением флюида.

Согласно другому аспекту, решения количественной оценки пласта определяются из многомерного распределения ЯМР данных путем первоначального определения набора параметров модели, которые представляют аспекты многомерного распределения. Затем, для вычисления свойств флюидов, применяется зависящая от модели инверсия.

Согласно другому аспекту, решения количественной оценки пласта определяются из многомерного распределения ЯМР данных способом указания и щелчка. Используя мышь компьютера или прикладную программу автоматического выбора края, из ЯМР данных выбирается один или несколько артефактов флюидов. Для того чтобы определить свойства флюида, ассоциированного с выбранной областью, амплитуда интегрируется в ее пределах.

Согласно другому аспекту, решения количественной оценки пласта определяются из многомерного распределения ЯМР данных посредством определения среднего значения диффузии по области диффузионно-T2-релаксационного распределения. Для определения свойств флюида, ассоциированного с выбранной областью, используется средняя диффузия.

Краткое описание чертежей

В дальнейшем изобретение поясняется описанием конкретных вариантов реализации изобретения со ссылками на сопровождающие чертежи, на которых:

фиг.1 изображает общую схему иллюстративной системы сбора данных ядерно-магнитного резонанса нисходящей скважины,

фиг.2 изображает более подробную схему системы согласно фиг.1,

фиг.3 изображает многомерную карту или график для отображения ЯМР данных,

фиг.4 изображает карту ЯМР данных, иллюстрирующую неидеальные эффекты на ЯМР данных,

фиг.5 изображает другую карту ЯМР данных, показывающую неидеальные эффекты на ЯМР данных,

фиг.6 изображает другую карту ЯМР данных, показывающую неидеальные эффекты на ЯМР данных,

фиг.7 изображает другую карту ЯМР данных, иллюстрирующую неидеальные эффекты на ЯМР данных,

фиг.8 изображает двухмерный график согласно способу инверсии согласно известному уровню техники,

фиг.9 изображает набор карт ЯМР данных, иллюстрирующих коррекцию неидеальных эффектов согласно аспекту описанного объекта изобретения,

фиг.10 иллюстрирует коррекцию способа инверсии, известного из уровня техники, согласно набору карт ЯМР данных по фиг.9,

фиг.11 изображает двухмерный график согласно способу инверсии из известного уровня техники,

фиг.12 изображает набор карт ЯМР данных, показывающих коррекцию неидеальных эффектов согласно аспекту описанного объекта изобретения,

фиг.13 иллюстрирует коррекцию способа инверсии из уровня техники согласно набору карт ЯМР данных фиг.12,

фиг.14 изображает набор карт и графиков ЯМР данных согласно одному из способов оценки пласта описанного объекта изобретения,

фиг.15 изображает блок-схему способа, иллюстрируемого на фиг.14,

фиг.16 изображает многомерную карту или график другого способа оценки пласта раскрытого объекта изобретения,

фиг.17 изображает набор карт и графиков ЯМР данных согласно способу, иллюстрируемому на фиг.16,

фиг.18 изображает блок-схему способа, иллюстрируемого на фиг.16 и 17,

фиг.19 изображает другой набор карт и графиков ЯМР данных способа оценки пласта описанного объекта изобретения,

фиг.20 изображает набор корректированных графиков ЯМР данных согласно способу, иллюстрируемому на фиг.19,

фиг.21 изображает снимок экрана, иллюстрирующий примерные решения оценки пласта, обеспеченные способом описанного объекта изобретения.

Подробное описание предпочтительных вариантов реализации изобретения

Описанный объект изобретения описывает количественные способы для интерпретации двухмерных карт ядерно-магнитного резонанса (ЯМР), полученных из обычных оценочных ЯМР измерений пласта. Хотя могут использоваться другие значения, предпочтительный вариант реализации изобретения, обсуждаемый здесь в первую очередь, основан на картах диффузии в функции от T₂(D-T2). Согласно настоящему изобретению, D-T2 карты могут использоваться для того, чтобы помочь в выборе параметров для применения к существующим кодам инверсии, основанным на модели. Далее, полные нефтефизические решения (пористость, проницаемость, объемы флюидов, насыщение, вязкость нефти и т.д.) могут быть получены непосредственно из D-T2 карт. Чтобы воспользоваться преимуществом визуального представления карт, предложенные способы являются интерактивными и, согласно варианту реализации изобретения, состоят из последовательных процедур указания и щелчка.

Сбор данных в ЯМР измерениях согласно вариантам реализации изобретения может быть выполнен с помощью различных способов ЯМР измерений, известных из уровня техники. Например, измерения могут быть выполнены в лаборатории с использованием образца, изъятого из земного пласта. Альтернативно, ЯМР измерения могут быть выполнены в процессе операции каротажа, с использованием инструмента проводной линии связи, инструмента каротажа во время бурения или измерения во время бурения, или же тестера пласта. Фиг.1 иллюстрирует схему системы ЯМР каротажа. На фиг.1 показан инструмент 30 ЯМР каротажа, предназначенный для исследования земных пластов 31, через которые проходит ствол 32 скважины. Инструмент 30 ЯМР каротажа подвешивается в стволе 32 скважины на армированном кабеле 33, длина которого по существу определяет относительную осевую глубину устройства 30. Длина кабеля управляется соответствующим средством, находящимся на поверхности, таким, как барабан и лебедочный механизм 8. Поверхностное оборудование 7 может быть обычного типа и может содержать подсистему процессора, которая сообщается с оборудованием нисходящей скважины, содержащим устройство 30 ЯМР каротажа.

Устройство 30 ЯМР каротажа может быть любым подходящим каротажным прибором ядерно-магнитного резонанса; оно может быть устройством для использования в прикладных задачах каротажа проводной линии связи, как показано на фиг.3, или же оно может быть из тех, которые могут использоваться в прикладных задачах каротажа в процессе бурения (LWD, КПБ) или измерения в процессе бурения (MWD, ИПБ). Кроме того, устройство 30 ЯМР каротажа может быть частью тестера пласта, известного в технике, например, такого, который продается под торговым названием MDT™ технологической корпорации Schlumberger (Хьюстон, TX). Устройство 30 ЯМР каротажа обычно содержит средство для создания статического магнитного поля в пластах и средство радиочастотной (РЧ) антенны для создания импульсов магнитного поля в пластах и для приема спиновых эхо от пластов. Средство для создания статического магнитного поля может содержать постоянный магнит или матрицу магнитов, а средство радиочастотной антенны для создания импульсов магнитного поля в пластах и для приема спиновых эхо от пластов может содержать одну или несколько РЧ антенн.

Фиг.2 иллюстрирует схему некоторых компонентов одного типа устройства 30 ЯМР каротажа. Фиг.2 изображает первый центральный магнит или матрицу 36 магнитов и РЧ антенну 37, которая может быть соответственно ориентированной катушкой или катушками. Фиг.2 также иллюстрирует общее представление близко расположенных тонких цилиндрических оболочек 38-1, 38-2… 38-N, которые могут быть выбраны по частоте в операции многочастотного каротажа. Одно такое устройство описано в патенте США № 4710713. На фиг.2 показан другой магнит или матрица 39 магнитов. По мере того, как устройство 30 каротажа поднимается в стволе скважины в направлении стрелки Z, для предварительной поляризации земного пласта перед областью исследования может использоваться матрица 39 магнитов. Примеры таких устройств раскрыты в патентах США № 5055788 и 3597681.

На фиг.3 показана иллюстративная D-T2 карта с ЯМР данными спинового эхо, представленными как амплитуды в функции от диффузии (D) и релаксации (T2). Карта, показанная в левой панели, представляет собой трехмерный вид в перспективе. Правая панель обеспечивает более практичное представление D-T2 карты в виде двухмерной карты. Однако, следует отметить, что раскрытые способы могут применяться к набору данных, имеющему несколько измерений: 2-D, 3-D, 4-D и т.д. Кроме того, следует отметить, что, хотя D-T2 карты обсуждаются здесь для иллюстративных целей, раскрытые способы могут быть одинаково эффективными в получении решений количественной оценки пласта на основе многих других комбинаций свойств ЯМР данных (D, T1, T2, T1/T2 и т.д.).

В контексте двухмерной D-T2 карты амплитуда диффузии представлена согласно схеме кодирования цветом. Различия диффузионных свойств между газом, водой и нефтепродуктами с различными вязкостями охватываются D-T2 картой и изображаются как отдельные и различимые пики. В частности, группирование цвета в области А, также упоминаемое здесь как артефакт или проявление флюида, представляет вероятное детектирование первого флюида. Аналогично, группирование более светлого цвета или проявлений флюидов в областях B, C и D также представляет вероятное детектирование трех дополнительных флюидов. Для того чтобы помочь интерпретации, на карты наложены теоретические отклики воды, нефти и газа. Таким образом, для группировки или проявления A, вероятно, что флюид представляет собой газ, поскольку пик находится вблизи теоретического значения диффузии газа. Для группировок B и C, вероятно, что флюиды имеют варьирующиеся вязкости или фазы нефти, лежащие вдоль теоретической линии диффузии нефти. Наконец, вероятно, что группировка или проявление D представляет собой воду, подверженную ограниченной диффузии (обсуждается ниже).

Согласно одному из вариантов реализации изобретения, сводный график фиг.3 не зависит от модели. Это означает, что для получения результатов ЯМР измерений из данных спинового эхо, в вычисления не вводятся заданные значения или пределы диффузии. Хотя существуют и другие примеры, один из примеров вычисления, независящего от модели, представляет собой ранее упоминаемый подход ПМЭ. На фиг.3, чтобы содействовать интерпретации ЯМР данных, перекрытие теоретических откликов воды, нефти и газа является полезным.

Как упоминалось, предыдущие попытки определить решения количественной оценки пластов были основаны на модели заданных значений параметров диффузии и времен T2 релаксации флюидов. Кроме того, требовалось выбирать модель флюидов, основанную на лучшей приблизительной оценке, относительно которой будут осуществлять детектирование флюидов. Разумеется, любые неточности в начальных оценках модели флюида и параметров флюидов создают неточности в конечных решениях. Фиг.4 графически иллюстрирует некоторые исключения для общих идеальных моделей: 1) внутренний градиент тянет данные в направлении стрелки вверх (так, что вода могла бы быть ошибочно принята за газ), 2) ограниченная диффузия тянет данные в другом направлении, показанном стрелкой вниз (так, что вода могла бы ошибочно быть принята за нефть), 3) смешанная смачиваемость тянет данные нефти влево (так, что объем/насыщение нефти вычисляются слишком низкими), 4) высокий газовый фактор (GOR, ГФ) заставляет сдвигаться нефть вдоль северо-восточного направления (так, что нефтепродукты могли бы ошибочно быть приняты за воду или газ).

Например, фиг.5 иллюстрирует действие эффекта внутреннего градиента на ЯМР отклик флюида. Внутренний градиент прибавляется к градиенту инструмента и, следовательно, диффузия оказывается выше ожидаемой. Образцом является богатый глиной песчаник в известной водонесущей зоне скважины, пробуренной с водным буровым раствором. Также, расчетные значения диффузии оказываются выше ожидаемого результата, показанного перекрытием линии диффузии воды.

Другой пример показан на фиг.6, иллюстрирующей действие ограниченной диффузии в известном пласте с доминирующим карбонатом. Свободная вода подвергается неограниченной диффузии в больших порах и согласуется с теоретическим откликом. Связанная вода задерживается в более мелких порах и поэтому испытывает ограниченную диффузию. Результатом является один пик, представляющий свободную воду, на ожидаемом перекрытии диффузии воды, и второй пик, представляющий связанную воду, вытянутый поперек диапазона значений диффузии, даже пересекая значения, показательные о наличие углеводородов.

Еще один пример, показанный на фиг.7, иллюстрирует газ в месторождениях песчаника с высокой проницаемостью на Ближнем Востоке. Большие размеры пор заставляют воду релаксировать с параметрами, близкими к ее объемному значению, накладываясь на сигнал газа и приводя к размазанному пику на длительном T2.

Согласно варианту реализации настоящего изобретения, такие D-T2 карты генерируются посредством кодов инверсии, независящей от модели, как, например, ПМЭ. Указанные инверсии, независящие от модели, не требуют априорного ввода значений диффузии или T2 флюидов. Согласно одному из вариантов реализации, инверсии, независящие от модели, не требуют ни априорного знания свойств флюидов, ни того, какие флюиды присутствуют. Для отображения результирующих многомерных ЯМР данных по двум, трем или более осям в легко считываемой форме из этих инверсий генерируются D-T2 карты или графики.

Согласно одному из вариантов реализации, для улучшения результатов инверсии, независящей от модели, как, например, ОХФМР анализ, используются D-T2 карты. В частности ПМЭ и другие D-T2 карты, полученные независимо от модели, обеспечивают несмещенное представление ЯМР данных. Также, перекрытие теоретических откликов трех наиболее часто встречающихся флюидов, т.е. воды, нефти и газа, показывает, нужна ли для откликов настройка для использования в инверсии, основанной на модели.

Например, ОХФМР модель (для нефти, газа и воды) формулирует, что коэффициенты диффузии воды и газа не зависят от T2, а зависят от температуры T и давления P (для газа):

Для нефти коэффициент диффузии линейно пропорционален T2,

Из уравнений (1)-(3) следует, что две горизонтальные линии (то есть при постоянных значениях D), представляющие теоретические отклики воды и газа, и диагональная линия, представляющая теоретический отклик нефти, могут быть наложены на D-T2 карту. Отклонения от идеальных откликов флюидов будут видны на картах как сигналы, расположенные вдали от линий перекрытия. Будучи известными, указанные отклонения применяются к инверсии, основанной на модели. Результат, полученный из инверсии, обеспечивает решение с улучшенной точностью на основе исследований независящих от модели D-T2 карт.

Фиг.8-10 иллюстрирует коррекцию инверсии, зависящей от модели, в водном буровом растворе, формате цезия, а фиг.11-13 иллюстрируют коррекцию в образце сырой первичной нефти. Фиг.8 изображает иллюстративный анализ, зависящий от модели, как, например, ОХФМР подход, для раствора формиата цезия. Фиг.8 изображает анализ, зависящий от модели, имеющий пик 200 нефти и пик 201 воды. Анализ, основанный на модели, использующий коэффициент диффузии воды, взятый по умолчанию, показывает главным образом нефть для бурового раствора формиата цезия. На фиг.9, слева, изображена D-T2 карта для таких же данных, имеющих проявление или артефакт 206, лежащий вдоль линии 204 диффузии нефти по умолчанию. В частности, изображенная слева D-T2 карта показывает, что детектированный формиат Ce, например 206, имеет намного более низкий коэффициент диффузии, чем чистая вода (по умолчанию), обозначенная горизонтальной линией 202. То есть можно определить, что диффузия 202 идеальной воды является неправильной для образца формиата Ce. Таким образом, для вычисления точных значений насыщения флюидов эффективный коэффициент диффузии "воды" (точнее, формиата Ce) должен быть уменьшен до 40% от его значения по умолчанию, представленного горизонтальной линией 208. Затем указанное скорректированное значение диффузии программируется в инверсию, основанную на модели. Фиг.10 изображает повторно обработанные ОХФМР результаты, полученные с использованием коэффициента диффузии воды, уменьшенного до 40% от его первоначального значения. Теперь ОХФМР результаты правильно показывают преобладающую воду, показанную на пике 210, и очень низкий сигнал нефти, показанный на пике 212. Также, могут быть более точно вычислены определения оценок ассоциированных пластов, такие как объем и насыщение воды.

Фиг.11 представляет анализ, зависящий от модели, известного образца сырой нефти. Результаты показывают наличие нефти, показанной пиком 220, но также неправильно показывают значительное насыщение воды, показанное пиком 221. Фиг.12 изображает D-T2 карту для тех же данных. Слева, изображена D-T2 карта с линиями перекрытия флюидов по умолчанию, с диффузией воды, представленной горизонтальной линией 222, и с диффузией нефти, представленной линией 224. Сигнал или проявление 226 появляется над D-T2 линией 224 корреляции нефти, показывающей, что данная нефть имеет необычно высокое отношение диффузии к T2. Таким образом, можно определить, что существующая D-T2 линия 224 нефти является некорректной и требует регулировки. Такая же карта построена справа с новой D-T2 линией 228 корреляции нефти, которая делит пополам главный пик сигнала. Затем проводится ОХФМР анализ, использующий увеличенное значение λ (см. Ур.(3)). Результаты повторно обработанных данных показаны на фиг.13. Как ожидалось, теперь ОХФМР анализ правильно прогнозирует преобладающую нефть на пике 232, при этом показывая очень низкое присутствие воды на кривой 230.

В дополнение к приведенным выше вспомогательным методам интерпретации, зависящей от модели, непосредственно из двухмерных карт могут быть получены полные количественные нефтефизические решения. В частности, согласно одному из применений изобретения, из D-T2 карт могут быть получены количественные измерения пористости, проницаемости, объемов флюидов, насыщения, вязкости нефти и других величин. По сравнению с предыдущими способами, для того, чтобы получить больше, чем ранее используемая качественная информация, необходима дополнительная интерпретация. Согласно раскрытому объекту изобретения, для получения количественных решений используются два подхода, подход указания и щелчка, а именно указания объекта на экране и выбор отрывистым нажатием кнопки мыши и подход логарифмического среднего значения диффузии.

Согласно одному из варианту реализации, обеспечивается подход указания и щелчка, который позволяет пользователю взаимодействовать с D-T2 картой посредством фокусировки на конкретных артефактах, графически изображенных на карте. Амплитуда сигнала A из серии ЯМР последовательностей импульсов может быть выражена как

где WT, TE, t представляют собой время ожидания, разнесение эхо-сигналов и время ЯМР последовательностей импульсов; (i, j k) представляет собой индексы T2, D и T1/T2 распределений, f(i,j,k) - амплитуда трехмерной составляющей в T2, D, T1/T2 пространстве, H(WT, TE, t,i,j, k) - ядро этой составляющей, а δ - шумовая составляющая.

Из вышеупомянутого уравнения (4) можно видеть, что D-T2 карта является представлением амплитуд сигналов, проинтегрированных по k размерности (T1/T2) результата ПМЭ инверсии. Следовательно, в благоприятном случае, когда распределения флюидов в D-T2 пространстве хорошо разделены, как видно на фиг.14, их соответствующие объемы могут быть непосредственно получены путем интегрирования амплитуд в D и T2 окнах, заданных каждым типом флюидов. Тогда соответствующие насыщения могут быть получены путем деления соответствующих объемов флюидов на полный объем флюидов. Для обеспечения правильных решений, отдельные объемы должны быть скорректированы на водородный индекс согласно известным способам.

Фиг.14 иллюстрирует, согласно одному из вариантов реализации, пример подхода указания и щелчка, имеющий четыре артефакта флюидов 240, 242, 244 и 246, которые хорошо разделены в D-T2 пространстве (левая верхняя часть). Заметим, что четыре флюида не разрешены ни в T2 пространстве (левая нижняя часть), ни в D пространстве (правая верхняя часть). В частности, графики T2 и D пространства показывают три вместо четырех флюидов, имеющих различимые свойства. В нижнем левом T2 графике пик 248 соответствует артефакту 246, а пик 250 соответствует артефакту или проявлению флюида 244. Однако, пик 252 не может разрешить артефакты 240 и 242. В верхнем правом распределении диффузии пик 254 соответствует артефакту 240, а пик 256 соответствует проявлению флюида 242. Однако, пик 258 не может разрешить артефакты 244 и 246.

Согласно предшествующим способам, количественные решения флюидов могли бы быть получены, используя T2 графики. Однако, как показано, в некоторых случаях T2 карты не могут полностью разрешить многочисленные флюиды, имеющие различные диффузионные свойства. Известные из уровня техники способы сделали шаг к дальнейшему использованию D-T2 карты, с тем, чтобы оценить точность (качественно) решений, полученных из T2 графиков. Согласно варианту реализации настоящего изобретения, визуальный способ указания и щелчка использует D-T2 карту для определения количественных соответствующих объемов каждого из четырех флюидов, путем интегрирования амплитуд сигналов вдоль T2 и D измерений в окнах, определенных прямоугольниками 260, 262, 264 и 266. Заметим, можно аналогично применять другие формы для очерчивания исследуемой области карты, например многоугольники или круги. Таким образом, раскрытые способы содействуют в продвижении состояния уровня техники в части разрешения многочисленных флюидов, имеющих подобные T2 распределения, а также посредством количественного определения некоторых решений оценки пласта, как только идентифицированы артефакты флюидов.

Согласно одному применению способа указания и щелчка, интерпретатор может легко выбрать область интегрирования D-T2 карты, используя, например, мышь компьютера или цифровое перо. Выбор области также можно осуществлять автоматически посредством реализации алгоритма программного обеспечения, основанного на заданном пороге амплитуды. Также заметим, что интерпретация типа флюида управляется путем перекрытия теоретических откликов газа, нефти и воды на D-T2 карте, как видно в левой верхней части. Указанный этап также может осуществляться автоматически посредством прикладной программы программного обеспечения, основанной, например, на близости точки максимальной амплитуды к теоретическим откликам газа, нефти и воды.

Фиг.5 иллюстрирует блок-схему подхода указания и щелчка. На этапе 500 начинается процедура указания и щелчка с D-T2 картой. Как упоминалось, указанная D-T2 карта предпочтительно генерируется с использованием подхода инверсии, независящей от модели. Далее, на этапе 502, пользователя просят ввести тип модели флюида, которую нужно использовать в вычислении. Обычно, указанная процедура осуществляется процессом осмотра карты. Сначала пользователь принимает решение о числе флюидов. Например, если присутствуют две группы амплитуд (яркие пятна), то выбирается подход двух флюидов. Во-вторых, каждый тип флюида интерпретируется с помощью перекрытия теоретических откликов воды, нефти и газа. На этапе 504 пользователь выбирает область флюида, используя мышь компьютера. Согласно одному из вариантов реализации, отображается рамка, которая управляется так, чтобы охватить по существу весь артефакт на D-T2 карте. Согласно варианту реализации, окончательная невыбранная область, имеющая положительную ненулевую амплитуду, накапливается и отображается для обеспечения пользователю индикации того, были ли пропущены какие-либо существенные части. Согласно другому варианту реализации, прикладная программа компьютера или программного обеспечения выбирает область, по меньшей мере, частично окружающую каждое проявление флюида на основе пороговой амплитуды.

Как только артефакт был выбран, на этапе 506 вычисляются объем и насыщение путем интегрирования по выбранной области D-T2 карты. Поскольку интегрирование по полной области карты дает полную пористость, значения насыщения отдельных флюидов могут быть вычислены путем деления объемов флюидов на полную пористость. Здесь уместно отметить, что хотя раскрытый способ указания и щелчка обсуждается в контексте D-T2 карт для иллюстративных целей, практически любая многомерная карта может быть использована для того, чтобы определить искомую величину. Далее на этапе 510, как только были определены величины объема флюида и насыщения или другие базовые значения, вспомогательные вычисления могут быть затребованы пользователем или автоматически алгоритмом. Например, на этапе 512 может быть определена вязкость путем вычисления среднего T2 в окне нефти, и посредством использования опубликованных диаграмм вязкости-релаксации можно оценить вязкость. Другой пример заключается в том, чтобы использовать объем связанного флюида, полученный из карты, для вычисления проницаемости Timur-Coates согласно уравнению k_Timur=a*phit^b*((phit - bfv)/bfv)^c, где а, b, c - константы, phit - полная пористость, а bfv - объем связанного флюида.

Второй подход для непосредственного определения количественных результатов из D-T2 карты применяет определение среднего логарифмического значения диффузии (D_LM). При стандартном ОХФМР анализе необработанные данные аппроксимируются непосредственно, используя ограничения уравнений 1, 2 и 3. Подобные ограничения накладываются для других инверсий, независящих от модели, уровня техники. Альтернативный подход, раскрытый здесь, заключается в том, чтобы использовать сами карты в качестве входных данных для получения решения, которое пытается получить ОХФМР. Поскольку информация, содержащаяся в картах, по существу идентична таковой у исходных данных, должны быть сопоставимы два способа решения. Однако, на практике данных часто не хватает в информации диффузии, и, следовательно, D-T2 амплитуды различных флюидов размазываются по большим площадям карты (разрешение). Такая ситуация противоположна идеальной ситуации для применения способа указания и щелчка, в котором каждый артефакт флюида по существу является отделенным от других проявлений флюида. Тогда проблема состоит в повторном присвоении разброса