EA001137B1

EA001137B1 - Способ определения жидкой фракции в земной формации

Info

Publication number: EA001137B1
Application number: EA199800826A
Authority: EA
Inventors: Роланд Йоханнес Мария Бонни; Пауль Хофстра; Вильхельмус Йоханнес Лойестейн; Роберт Карл Йосеф Сандор
Original assignee: Шелл Интернэшнл Рисерч Маатсхаппий Б.В.
Priority date: 1996-03-14
Filing date: 1997-03-13
Publication date: 2000-10-30
Also published as: CN1185507C; BR9708038A; EP0886792B1; EA199800826A1; WO1997034166A1; AU2027697A; CN1213433A; EP0886792A1; MY122012A

Description

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

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

Наблюдаемая кривая затухания намагничивания высоты эхо является существенно многоэкспоненциальной функцией и может быть представлена:

Μ(ί) = Μ(0){ (1) о

где Ρ(Τ₂)άΤ₂ представляет фракцию жидкости с временем поперечной релаксации от Т₂ до Τ2+άΤ.

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

ЕР-А-489578 раскрывает способ проведения измерения эхо-импульса ЯМР в буровой скважине. В этом способе статическое магнитное поле и градиент магнитного поля прилагаются к формации, окружающей буровую скважину, сопровождаемые серией электромагнитных импульсов ЯМР. Затухание эхо ЯМР используется, чтобы определить коэффициент диффузии Ό и поперечное время релаксации Т₂. Эта публикация раскрывает, что коэффициент диффузии Ό может применяться, чтобы определить уровни насыщения воды и углеводорода. Однако не раскрыто, как такие уровни насыщения воды и углеводорода могут быть определены. Кроме того, только упоминается единственный коэффициент диффузии Ό, который относится или только к воде, или только к углеводородной жидкости.

Целью изобретения является получение улучшенного метода определения жидкой фракции в земной формации, содержащей, по крайней мере, две жидкости.

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

a) индуцирование магнитного поля в области упомянутой земной формации;

b) проведение измерения эхо-импульса ЯМР в вышеупомянутой земной формации;

c) выбор соотношения между эхооткликом ЯМР от вышеуказанных жидкостей, фракций жидкостей и, по крайней мере, одной переменной, которая воздействует на эхоотклик ЯМР способом, зависящим от фракций вышеупомянутых жидкостей;

4) изменение вышеуказанной, по крайней мере, одной переменной в ходе измерения ЯМР, воздействующего таким образом на измеренный эхо-отклик ЯМР способом, зависящим от фракций вышеуказанных жидкостей; и

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

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

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

Вариации в величине произведения могут быть получены за счет вариаций в импульсном интервале ЯМР и/или градиенте магнитного поля.

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

Предпочтительно, чтобы, по крайней мере, два импульса временных интервала ЯМР различной длины применялись в измерении ЯМР, эти импульсные временные интервалы различной длины соответственно применяются в одиночной ЯМР каротажной последовательности, в одиночном каротажном проходе циклично или в отдельных каротажных проходах.

Подходящая ЯМР каротажная последовательность - это Сагг-Ригсе11-Ме1Ьоош-ОШ (СРМО) последовательность в соответствии с:

ТК - 90° ± х - (±ср4 - 180°у - Ц,, - есйо,) (2) где ТК - время ожидания между последовательностями;

1_ср| - Сагг-Ригсе11 интервал;

х, у - это фазы 90° и 180° импульсов соответственно; и

I - индекс эхо-сигнала.

Для одиночной жидкости кривая затухания эхо ЯМР при градиенте магнитного поля и постоянной 1_ср может быть описана Μ(ΐ)=Μ(0)ΣΑίβχρ(-ΐ/Τ_2>ί)βχρ(-1γ² ϋ О² 1ср²/3)(1-ехр(-ТВ./Тм)) (3) где Άί - жидкая фракция с временем поперечной релаксации Т₂д;

М(0) - амплитуда сигнала в момент времени ΐ = 0;

Т₂д - время поперечной релаксации жидкой фракции Ль

Тц - время продольной релаксации жидкой фракции Ль γ - гиромагнитное отношение подчиненного ядра жидкости;

О - градиент магнитного поля;

Ό - молекулярный коэффициент самодиффузии жидкости в пористой породе.

Кривая эхо-сигнала ЯМР из формации, которая содержит множество жидкостей, есть суперпозиция эхо-сигналов, произведенных индивидуальными жидкостями в соответствии с

Μ(ΐ)=ΣΜ/ι) (4) где М/ΐ) является затуханием эхо ЯМР )-той фракции жидкости, как описано в выражении (3).

Величина коэффициента диффузии Ό связана с температурой и вязкостью жидкости и может быть аппроксимирована эмпирическим соотношением

1)2.5Т 3001-| (10'⁹ м²/с) (5) где Т - температура (К);

η - вязкость (сП).

Молекулярная диффузия в пористой среде по существу ограничена, что подразумевает, что Ό не константа, а функция эффективного времени диффузии, заданного 21_ср, и геометрией пористой системы.

Будет ясно, что различие между эхосигналами от различных жидкостей может быть сделано на основании влияния различных коэффициентов диффузии на отклике ЯМР (изменяя импульсный временной интервал) и/или на основании влияния их различных времен продольной релаксации на отклике ЯМР (изменяя ТК).

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

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

Шаг е) способа согласно изобретению предпочтительно содержит применяемый метод инверсии к измеренному эхо-отклику ЯМР и вышеуказанному выбранному соотношению.

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

Пример.

Определение водного насыщения в формации горных пород, содержащей нефть со средним удельным весом и воду, применяя измерения ЯМР градиента магнитного поля на образце формации горных пород, иллюстрируется в дальнейшем со ссылкой на прилагаемые чертежи, на которых фиг. 1 -4 схематично показывают кривые затухания ЯМР для земных формаций, имеющих нефтяные насыщения соответственно 0,1 5, 0,3, 0,45 и 0,6; и фиг. 5 - отношение эхо-сигналов ЯМР для двух различных импульсных интервалов времени как функцию нефтяного насыщения при различных значениях времени.

Два времени поперечной релаксации и две соответствующие фракции объема были отобраны для модели воды. Компонент с коротким временем релаксации представляет связанную воду и остается постоянным, принимая во внимание, что компонент с длинным временем релаксации представляет подвижную воду и изменяется с нефтяным насыщением. Нефть была смоделирована одним временем поперечной релаксации и одной соответствующей фракцией объема. Значения параметров были выбраны, как следует ниже:

Т2,«,1 = 10 мс; = 0,25;

Т2,«,2 = 100 мс;

Л„₂ = 0,60; 0,45; 0,30; 0,15 соответственно; = 3,0 Е-9 м²/с;

Т₂,₀ = 50 мс;

Ао = 0,15; 0,30; 0,45; 0,60 соответственно;

= 0,2 Тл/м;

γ = 2π 42,565 МГц/Тл.

Каждая из фиг. 1-4 показывает нормализованные кривые затухания эхо ЯМР для двух импульсных интервалов (2Ц,), где в каждой из этих фигур верхняя кривая (обозначенная 1а, 1Ь, 1с, 16) представляет ЯМР отклик как функцию времени 1 для импульсного интервала времени 21_ср = 2 мс, и более низкая кривая (обозначенная цифрами 2а, 2Ь, 2с, 26) представляет отклик на функцию времени 1 для импульсного интервала времени 21_ср = 6 мс. Кроме того, нефтяное насыщение на фиг. 1 равно А₀ = 0,15, на фиг.2 - А₀= 0,3, на фиг. 3 - А₀ = 0,45 и на фиг. 4 - А₀ = 0,6.

Как ясно из этих фигур разделение между верхней кривой 1а, 1Ь, 1с, 16) и более низкой кривой (2а, 2Ь, 2с, 26) уменьшается с увеличением нефтяного насыщения. Таким образом, применяя, по крайней мере, два различных импульсных интервала времени, было определено разделение между верхней и более низкой кривой, а из разделения было определено нефтяное насыщение.

В то время как предложено, что водное насыщение (или нефтяное насыщение) было определено подбором полных кривых к соответствующей форме уравнения (4), используя соответствующую программу числовой минимизации, чувствительность способа, согласно изобретению, может быть оценена из фиг. 5, которая показывает кривые 3, 4, 5, 6, представляющие отношение В = М(1, 21_ср=6)/М(1, 21_ср=2) для выбранных значений 1, как функцию нефтяного насыщения А0. Выбранное значение 1 для кривой 3 - 1=90 мс, для кривой 4 - 1=48 мс, для кривой 5 - 1=18 мс и для кривой 6 - 1=12 мс.

Повторяя способ для диапазона практических значений для водных и нефтяных параметров, было определено, что способ, согласно изобретению, лишь слабо зависит от фактических величин нефтяных параметров. Если никакая информация относительно этих параметров не доступна, погрешности в оценке водного насыщения могут доходить до 0,1. Если вязкость нефти может быть оценена с точностью двух десятичных цифр, результирующая погреш ность в водном насыщении незначительна по сравнению с полной точностью измерения.

Способ может быть выполнен, например, на месте через буровую скважину, образованную в земной формации или в лаборатории, используя образец керна, взятый из формации.

Кроме того, способ привлекателен для определения остаточного нефтяного насыщения в случае, если течение воды или движение газов применялось к земной формации для замещения нефти в формации.

Claims

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

a) индуцирование магнитного поля в области земной формации, которое имеет градиент магнитного поля;

b) проведение измерения эхо-импульса ЯМР в упомянутой области земной формации;

c) выбор соотношения между эхооткликом ЯМР от вышеупомянутых жидкостей, фракций жидкостей, различных коэффициентов ЯМР упомянутых жидкостей и, по крайней мере, одной переменной, которая воздействует на ЯМР эхо-отклик в зависимости от фракций упомянутых жидкостей;

6) изменение упомянутой, по крайней мере, одной переменной в ходе измерения ЯМР, и таким образом воздействия на измеренный эхоотклик в зависимости от фракций упомянутых жидкостей; и

е) определение фракции выбранной жидкости путем аппроксимации выбранного соотношения на ЯМР эхо-отклик.

2. Способ по п.1, в котором одна из упомянутой, по крайней мере, одной переменной является произведением импульсного интервала ЯМР и градиента магнитного поля.

3. Способ по п.2, в котором в измерении ЯМР используют, по крайней мере, два импульсных интервала различной длины.

4. Способ по п.1, в котором, по крайней мере, одна переменная содержит время ожидания между ЯМР последовательностями импульсов.

5. Способ по любому из пп.2, 3, в котором упомянутое произведение импульсного интервала ЯМР и градиента магнитного поля, формирует первую переменную, и в котором упомянутая, по крайней мере, одна переменная, включает вторую переменную, являющуюся временем ожидания между ЯМР последовательностями импульсов.

6. Способ по любому из пп.1-5, в котором этап е) включает применение инверсионного метода к измеренному эхо-отклику ЯМР и вышеуказанному соотношению.

7. Способ по п.6, в котором измеренные данные эхо-отклика ЯМР, полученные изменением первой и второй переменных, инвертированы одновременно.

8. Способ по любому из пп.1-7, в котором вышеуказанное измерение ЯМР содержит каротажную последовательность в соответствии с

ТК - 90° ± х - (Ц^ - 180°у - - еко]), где ТК - время ожидания между последовательностями;

- интервал каротажной последовательности;

х, у - фазы 90° и 180° импульсов соответственно; и

I - индекс эхо-сигнала.