Тонкослоистые пласты
.pdfvk.com/club15268505072 Практuческuе аспекты| vk.com/id446425943геофизических исследований скважин
в случаях, когда расхождения объяснить не удается, можно считать разум ным подход, при котором наиболее пессимистический набор функций при нимается в качестве низшей границы, наиболее оптимистический набор -
в качестве верхней, а посредине между ними строится базовая линия.
После определения функции, описывающей базовую линию, для пере счета значений насыщенности следует в обратном направлении с помо
щью J-функции повторно интерпретировать каротажные данные. Учи тывая, что эта функция использует только кривые для высоты над уров нем свободной воды, пористость и проницаемость в качестве входных данных (полученных из диаграммы плотностного каротажа), результаты
могут оказаться неожиданными. Имея опыт применения данной мето
дики для более чем 20 залежей, я всегда удивлялся тому, как логично
выглядят результаты и насколько хорошо они зачастую соответствуют
значениям Sw' полученным на основе данных каротажа сопротивления,
за исключением тонких прослоев, где значение Sw по данным каротажа
сопротивления обычно бывает слишком завышено. Диаграмма плотнос тного каротажа в вертикальном направлении обычно имеет разрешение
1 фут, намного превосходящее разрешение индукционного или бокового
каротажа. Результирующее увеличение средней нефтеили газонасыщен
ности составляет, как правило, порядка 20%.
Я полагаю, что хорошо откалиброванная J-функция Леверетта, при
водящая в соответствие друг с другом надежные каротажные данные по
насыщенности и данные кернового анализа, представляет собой наилуч
ший возможный способ расчета насыщенности коллектора и передачи
этой информации в статические и динамические модели. Следует также
осознавать, что в случае, когда данные кернового анализа надежны, но
местоположение уровня свободной воды неизвестно, эта функция может
быть использована для прогнозирования глубины уровня свободной во
ды, исходя из каротажных данных по насыщенности. Вообще говоря, сам
я не достиг особых успехов в применении этого метода, но он, возможно,
заслуживает дальнейшего исследования.
Упражнение 4.2. Получение J~Функции на основе данных каротажа
1.Для получения J-функции на основе данных каротажа используйте
результаты выполнения упражнений 2.1 и 3.1.
2.Проведите повторный расчет суммарных и средних значений с помо
щью J-функции, полученной на основе данных каротажа.
3.Которую из двух функций: полученную на основе кернового анализа
или на основе данных каротажа, - вы рекомендуете использовать для
будущего определения STOIIP? На чем основано это предпочтение?
vk.com/club152685050 | vk.com/id446425943
ГЛАВА 5
СОВРЕМЕННЫЕ МЕТОД".
ИНТЕРПРЕТАЦИИ ДАННЫХ rЕОФИЭИЧЕСКИХ
ИССЯЕДОВАНИЙ СКВA)I(ИН
5.1. АНАЛИЗ ХАРАКТЕРИСТИК ГЛИНИСТОГО ПЕСЧАНИКА
Наличие глин в пласте усложняет аналитическую работу петрофизика,
поскольку, как правило, они имеют высокую электропроводность и за
нижают электрическое сопротивление пород, насыщенных углеводоро
дами. Кристаллы глины связывают воду, которая адсорбируется на их
поверхности, а также катионы (например, натрия), окруженные гидрат
ной водой. Это приводит К росту электропроводности глинистого песча
ника, содержащего углеводороды, по сравнению с песчаником, в котором
глины отсутствуют.
Использование уравнения Арчи для глинистого песчаника приводит
к слишком высоким значениям воДонасыщенности Sw' вследствие чего зоны, потенциально богатые углеводородами, могут быть не замечены.
В прошлом предлагалось достаточно много уравнений для учета завы шенной электропроводности, возникающей вследствие наличия диспер
сных глин, занижающих сопротивление и приводящих к завышению ве
личины Sw при использовании уравнения Арчи. Я планирую привести
все эти уравнения, однако детально предлагаю исследовать только один
метод, а именно модификацию метода Ваксмана-Смита (Waxman-Smit).
Этот метод я успешно применял на нескольких месторождениях, и его
преимущество состоит в том, что он не требует дополнительных иссле дований керна для калибровки модели (хотя эти данные и могут быть
включены в нее).
73
vk.com/club15268505074 Пракmuческuе |аспектыvk.com/id446425943геофизических исследований скважин
(5.1.1)
где В - константа, относящаяся к температуре, а Qv - катионообменная
емкость удельного объема порового пространства. Здесь т" и n" имеют
тот же самый смысл, что и в методе Арчи, но получают их другими мето
дами исследования керна. Температура определяет добавочную проводи- мость глины (в омах) на единицу Qv(B см3 ).
Уравнение, которое может быть использовано для связи В с темпера турой пласта и удельным сопротивлением воды R." было опубликовано
Томасом (Thomas) и выглядит следующим образом:
в = |
(-1,28 + О,255хТ - о,оОО4059хТZ)Ю + (О,О4хТ - O,27)XR.,I,23), |
|
(5.1.2) |
где |
Т - температура в градусах Цельсия. |
Qv (в мг-экв/ед. порового объема в см3 ) относится К катионообменной
емкости глины (cation exchange capacity, СЕС) (в мг-экв/lоо г), измеряе
мой в лабораторных условиях. Qv может быть получена из СЕС по фор
муле: |
|
Qv = СЕС х плотность/(IОО х ф). |
(5.1.3) |
СЕС может быть измерена методом химического титрования раздроб
ленных образцов керна. Величина СЕС зависит от типа глины. Ее типич ные значения для различных глин приведены в таблице 5.1.1.
Необходимо отметить, что данные исследований раздробленных об разцов керна с большой степенью вероятности могут быть непредстави тельными, поскольку в процесс дробления проявляется гораздо больше
типов катионного обмена, чем в пласте. Более того, мне никогда не удава
лось связать Qv с каким-либо параметром, полученным на основе данных каротажа (например, спористостью, Vsh )'
Дополнительная неопределенность связана с коэффициентом В, по
лучаемым, как правило, с помощью «стандартной» корреляции, которая
Таблица 5.1.1
Типичные свойства rлины
|
СЕС, |
|
Плотность эерен, |
ВОДОРОДНЫЙ |
|
Тип rлины |
мr-экв/l |
00 r |
r/См |
З |
индекс |
|
|||||
|
|
|
|
|
|
Каолинит |
3 |
-15 |
2,64 |
0,37 |
Иллит |
10-40 |
2,77 |
0,09 |
|
|
|
|
|
|
Монтмориллонит |
80 |
-150 |
2,62 |
0,12 |
|
|
|
|
|
Хлорит |
1 |
-30 |
3,0 |
0,32 |
|
|
|
|
|
vk.com/club152685050 | vk.com/id446425943 |
75 |
Современные методы интерпретации данных геофизuческих исследований скважин |
может быть, а может и не быть применима. Гораздо лучше подход, в рам
Kax которого из каротажных данных для водонасыщенного песчаника
с известными свойствами выводится комбинированный коэффициент BQv. Ниже я даю описание этого полезного метода. Прежде всего, делает ся предположение о том, что BQv удовлетворяет уравнению вида:
BQv = (ф,- ф)/(С х ф). |
(5.1.4) |
Фактически это уравнение предполагает существование «чистой по
ристости» (фJ (clean porosity - пористость скелета породы), а диспер
сная глина распределяется в поровом пространстве этого скелета, что
приводит к уменьшению измеряемой на образцах пористости по отно
шению к скелетной. Предполагается также, что добавочная электропро водность BQv должна быть связана с этим пропорциональным уменьше
нием пористости через константу С.
Если уравнение 5.1.3 подставить в уравнение 5.1.1 и присвоить Sw зна чение 1 (для воДонасыщенных песчаников), это уравнение приобретет
следующий вид:
Cwa = ф-m"/R,= 1/~ + (фjС) х (l!ф - 1/фJ |
(5.1.5) |
Следовательно, если для водоносных песчаников построить график Cwa
Взависимости от l/ф, точки должны ложиться на прямую линию:
•Значение Cwa В начале массива данных равно 1/~.
•Значение (l/ф) в начале массива данных равно l/фс
•Точки должны ложиться на прямую с наклоном (фjС).
Отметим, что в отсутствие эффектов глинистости точки должны об
разовывать ровную горизонтальную линию, с помощью которой для Cwa
можно получить значение 1/~, если коэффициент С устремить в беско
нечность. Если данные не образуют единого тренда, то этот метод можно
считать непригодным. Пример графика с такими данными показан на
рис. 5.1.1.
Этот график был построен в предположении, что т* = 2,0. Из него
можно определить, что: ~ = 0,05 Омм; ф, = 0,13 и С = 0,01. Используя
предполагаемое значение n *, можно рассчитать насыщенность методом
Ваксмана-Смита в углеводородной части пласта.
Если доступны данные SCAL (special care analysis - специальный ана лиз керна), можно вывести т* и n* из экспериментальных данных, как
это изложено ниже. При традиционном измерении показателя цемента
ции (т) параметр пористости (Р - formation resistivity factor) наносится на график в зависимости от пористости (ф) в двойном логарифмическом масштабе. В рамках модели Арчи справедливо следующее соотношение:
Современные методы интерпретации данных геофизических исследований скважин |
77 |
vk.com/club152685050 | vk.com/id446425943 |
Следовательно, если построить график 10g(I) в зависимости от 10g(Sw)'
то его наклон определит значение n.
В случае модели Ваксмана-Смита необходимо получить величину 1*,
задаваемую уравнением:
(5.1.9)
Построение зависимости друг от друга функций 10g(I*) и 10g(Sw) позво
ляет получить уточненный показатель насыщенности n*. Как и в случае
с m*, показатель n* всегда оказывается ниже, чем n, примерно на 0,2. Значение насыщенности углеводорода не очень чувствительно к зна
чению т * при условии, если для калибровки значения ~ используется
воДонасыщенный песчаник. Следует отметить, что полученное значение ~ может не соответствовать той минерализации, которую получают при
анализе воды, поступающей из пласта при опробовании и эксплуатации
скважин. На это значение влияет не участвующая в притоке находящая
ся в глине связанная вода, которая может оказаться более пресной, чем свободная. Поэтому величина ~, как правило, оказывается выше ожида
емой.
Эффект от использования модели Ваксмана-Смита обычно оказыва ется существенным только для относительно высоких значений ~. Это происходит потому, что множитель RJ3Qv становится мал по сравнению
с единицей, если ~ мало (высокоминерализованная вода). В этой ситуа
ции расчетное значение Sw будет лишь слегка отличаться от соответству
ющего значения, рассчитанного на основе модели Арчи.
Отметим, что при использовании этого уравнения из-за того, что Sw
появляется в обеих его частях, возникают определенные вычислительные проблемы. Эти проблемы легко преодолеваются следующим образом.
Для начала сделайте предположение о том, что Sw в правой части уравне
ния 5.1.1 равно 1. Рассчитайте Sw и подставьте в правую часть уравнения его новое значение. Продолжайте эту процедуру до тех пор, пока в ходе
последовательных итераций Sw' находящаяся в левой части уравнения,
не перестанет изменяться более чем на 0,001. Обычно для этого бывает
достаточно около пяти итераций.
Другой способ применения метода Ваксмана-Смита - так называе
мый метод нормированной Qv' предложенный Иштваном Юхасом (Istvan Juhasz). Читателям рекомендуется прочесть соответствующую статью,
опубликованную Обществом петрофизиков и интерпретаторов каро тажных диаграмм, в которой подробно рассказывается об этом методе
(см. ссылки). Ее сокращенный текст приводится ниже. Юхас показывает,
что уравнение Ваксмана-Смита может быть преобразовано к виду:
|
|
Практuческuе аспекты геофизических исследований скважин |
|
vk.com/club15268505078 |
| vk.com/id446425943 |
||
|
|
|
(5.1.10) |
где |
|
|
|
Qvn = Vsh Х Фsh/Ф; |
(5.1.11) |
||
11R.,sh |
= ф-mIRSh ; |
(5.1.12) |
|
Rsh |
- |
удельное сопротивление глины; |
|
Фsh |
- |
пористость глины. |
|
Параметр т" может быть определен с помощью графика зависимости log(R,) от lоg(ф) в водоносных глинистых пластах (не чистых пластах),
поскольку наклон его линии эквивалентен т". Параметры R.,sh и R., могут
быть определены из графика зависимости Cwa от Qvn' поскольку точка пе
ресечения оси Cwa линией графика для Qvn = О дает величину 1IR." а для Qvn = 1,0 - величину R.,Sh' Другие часто используемые уравнения:
Модель «двойной ВОДЫ»
1IR,= S;/(F х R.,) + BQv х SjF, |
(5.1.13) |
где
(5.1.14)
Симанду (Siтandoux)
|
|
|
|
|
|
|
|
|
(5.1.15) |
где |
Rsh - |
удельное сопротивление глины. |
|
|
|
||||
Индонезийское уравнение |
|
|
|
|
|
||||
S |
= R(-I1n)x |
, |
ГR+ |
|
(mI2) / ~(Ax R) |
}( 2111) |
. |
|
|
V(IUI",I2) / |
<Р |
|
(5.1.16) |
||||||
w |
t |
{ |
5}, |
" l\.sh |
,v |
|
|
||
Поскольку модель Ваксмана-Смита удовлетворяет всем требованиям, которые я предъявляю к уравнению для глинистого песчаника (то есть
она вводит проводимость глины, соотнесенную с количеством глины на
каротажных диаграммах), мне очень редко в своей работе приходилось прибегать к другим уравнениям. Как уже говорилось в главе 3, я предпо чел бы всегда получать насыщенность для определения начальных геоло
гических запасов товарной нефти (stock tank oil initially in place, STOIPP)
и начальных геологических запасов газа с по-
vk.com/club152685050 | vk.com/id446425943 |
79 |
Современные методы интерпретации данных геофизических исследований скважин |
|
мощью зависимости «насыщенность - высота», откалиброванной качес
твенными измерениями параметров керна.
Упражнение 5.1. Анализ параметров rлинистоro песчаника
1. Используя данные воДонасыщенной части пласта в скважине test1,
постройте зависимость BQv от пористости. Получите с помощью это го графика значение R,.
2. Используя приведенные в приложении 2 данные кернового анализа,
рассчитайте величины т* и n *.
3. Рассчитайте насыщенности с помощью уравнения Ваксмана-Смита
и составьте новую таблицу суммарных и средних значений парамет
ров пластов.
4. Как согласуются насыщенности, полученные методом Ваксмана-Сми
та, методом Арчи и методом анализа керна с привлечением J-функ
ции?
5.2. КАРБОНАТЫ
Большинство методов интерпретации, описанных в предыдущих гла
вах, хорошо работают и для карбонатных пластов. Однако плотность зе
рен, используемая для получения пористости из диаграммы плотностно
го каротажа, должна иметь другое значение. Для известняков ее наиболее
распространенным значением является 2,71 г/см3 , хотя по возможности лучше выполнять ее непосредственные измерения на керне. Доломит
и ангидрит характеризуются более высокой плотностью скелета породы
(см. приложение 6). На многих месторождениях с высококачественными карбонатными коллекторами оценка коллекторов осуществляется про
ще, чем в случае песчаника, благодаря отсутствию проблемы глинистос
ти. Пласты же при этом обычно бывают большей толщины. Вообще го
воря, в этих условиях очень хорошо работает модель Арчи.
Однако за время моей деятельности мне приходилось работать и на та ких низкопроницаемых газо- и нефтеносных карбонатных коллекторах,
где нормальные методы не давали хороших результатов. Это я и хочу обсудить ниже. Типичные характеристики низкопроницаемых « 1 мД) карбонатных коллекторов:
•Пористость поддается оценке, но лишь в виде изолированных зон или
каверн. Это хорошо бывает видно на примере мелоподобных коллек
торов.
vk.com/club152685050 | vk.com/id446425943 |
|
80 |
Пракmuческuе аспекты геофизических исследований скважин |
•Высоты капиллярного подъема и переходные зоны могут быть исклю
чительно протяженными (до 100 м).
•Проницаемость скелета породы может быть чрезвычайно низкой, но скважина может фонтанировать благодаря присутствию в пласте ес
тественных (или созданных механическим путем) трещин. Из-за то го, что эти трещины могут быть причиной значительных поглощений промывочной жидкости при бурении, скважина может нуждаться
в выносе наружу большого количества бурового раствора, прежде чем начнут фонтанировать углеводороды пласта.
Стратегия разработки такого месторождения зачастую предусматри
BaeT бурение горизонтальных скважин перпендикулярно ориентации ес
тественных трещин, позволяющее оптимизировать добычу в скважинах.
STOIIP или GIIP будут связаны главным образом с пористостью скелета породы, но поступать в скважину нефть и/или газ будут через эти тре щины. Следовательно, при отсутствии обширной системы трещин, обес печивающих дренирование матрицы, скважины будут либо обводняться, либо очень быстро глохнуть.
Очевидно, что умение обнаруживать трещины и определять их ори
ентацию имеет чрезвычайно большое значение. Во время бурения бы
BaeT полезно контролировать уровень потерь и фиксировать глубины,
на которых эти потери начинают происходить. Это может дать ценную
информацию о наличии трещин. Нет нужды говорить о том, что меры,
предпринимаемые при бурении для устранения потерь бурового рас
твора, должны гарантировать, что трещины не будут необратимо заглу шены. К традиционным каротажным методам, способным обнаруживать
трещины, относятся:
•Акустический каротаж (наличие пропуска циклов).
•Микрокаротаж (неустойчивое поведение при пересечении трещины
прижимным башмаком каротажного зонда).
•Гамма-каротаж (пики, возникающие в тех местах, где трещины оказа
лись заполненными радиоактивными минералами).
•Кавернометрия (сечение ствола буровой скважины стремится принять
форму эллипса, главная ось которого направлена перпендикулярно
к направлению действия тектонического напряжения).
Чтобы определить тип и ориентацию трещин, необходимо использо вать средства визуализации (микросканеры). Приборы каротажа сопро
тивления предпочтительнее, поскольку позволяют легче отличать запол
ненные флюидом открытые трещины от сцементированных. Однако для
этого подходят и ультразвуковые приборы, причем они представляют
собой единственно возможный вариант при использовании бурового
раствора на нефтяной основе.
vk.com/club152685050 | vk.com/id446425943 |
81 |
Современные методы интерпретации данных геофизических исследовании скважин |
|
Заметим, что естественные трещины стремятся ориентироваться в на
правлении максимума горизонтально направленного напряжения, су
ществующего в данных отложениях. В частности, в областях, соседству
ющих со значительными тектоническими нарушениями разница между
напряжениями, действующими в различных направлениях, может быть очень большой и способствующей растрескиванию пород. Бесценными для описания любых трещин могут оказаться керны, хотя необходимо
позаботиться о том, чтобы исключить из рассмотрения техногенные трещины, возникшие в результате бурения. Необходимо также, чтобы
совместно с керноотборником работало устройство для определения
ориентации керна. После завершения анализа системы трещин полезно
получить кривую плотности расположения трещин, которая может рас
сматриваться в качестве одной из каротажных диаграмм. Такая кривая может быть увязана с горизонтальной проницаемостью kh, полученной
в процессе испытания скважины и используемой для прогноза продук
тивности будущих скважин при их вскрытии.
В отличие от коллекторов, находящихся в терригенных породах, кар
бонатные коллекторы хорошо поддаются кислотной обработке с помо щЬЮ НС! либо совместно с гидроразрывом, либо без него. Однако такой вид обработки может влиять лишь непосредственно на призабойную
зону, и не способен создать высокую проницаемость в более широкой
области резервуара. Измерение давления и отбор проб в низкопрони цаемых карбонатных коллекторах с помощью обычного пластоиспы тателя на кабеле практически всегда является безуспешным. Чтобы
получить хоть какую-то надежду на успех, необходимо использовать
устройство пакерного типа, и даже тогда вероятность успеха, как пра
вило, невысока.
В течение срока разработки месторождения, приуроченного к мелопо добным коллекторам, могут возникать осложнения, связанные с их уп
лoTHeHиeM. Несмотря на то что оно обеспечивает дополнительный меха
низм поддержания пластового давления, на стадии проектирования их
разработки, особенно особенно в шельфовой зоне, потребуются всесто ронние исследования этой проблемы.
5.3. ПОЛИМИНЕРАЛЬНЫЕ (СТАТИСТИЧЕСКИЕ) МОДЕЛИ
Как указывалось в главе 2, наиболее предпочитаемым мной методом
расчета пористости является диаграмма плотностного каротажа, исполь
зующая соответствующие плотности матрицы и флюида. Этот метод мо
жет оказаться крайне ошибочным, если в пласте в различных количест
вах присутствуют также и более тяжелые минералы. Мне приходилось
работать на некоторых месторождениях, где породы коллектора содер-
