СПОСОБ ОЦЕНКИ ТРЕЩИННОЙ ПОРИСТОСТИ ПО ДАННЫМ СКВАЖИННОЙ СЕЙСМОРАЗВЕДКИ
METHOD OF ESTIMATING FRACTURE POROSITY BASED ON BOREHOLE SEISMIC DATA

Изобретение относится к нефтедобывающей промышленности, а именно - к оценке площади нефтегазоносности залежей с трещинными коллекторами методами определения пористости.
Для решения задачи определения трещинной пористости привлекаются методы корреляционной зависимости между скоростями, частотами и амплитудами продольных сейсмических волн. Результат не всегда однозначен, поскольку упругие параметры среды находятся в сложной взаимосвязи с ее геологическими параметрами, что может привести к ошибочным выводам.
В настоящее время широко применяется метод оценки структуры порово-трещинного пространства на основе измеренных кинематических параметров разнополяризованных волн различных типов (обратная геофизическая задача). Метод заключается в интерпретации соотношений между микровключениями (порами, кавернами и трещинами) с учетом их ориентации в пространстве и характером заполнения по результатам математического моделирования трещинной пористости исследуемого коллектора. По результатам этой интерпретации производится оценка относительной плотности каждого из участвующих в решении микровключений и оценка абсолютной пористоти и плотности породы, содержащей эти микровключения (Кузнецов P.M., Жуков А.П., Шнеерсон М.Б. Введение в сейсмическую анизотропию: теория и практика. - Тверь, ООО «Издательство ГЕРС», 2006, с.119-135).
Известный метод обеспечивает прогноз пористости и проницаемости, а соответственно - и продуктивности изучаемого интервала. При этом данные о системе преимущественной трещиноватости обеспечивают выбор оптимального направления при горизонтальном бурении с целью повышения продуктивности коллекторов.
Однако известный метод является трудоемким. Для обеспечения достаточно точных расчетов требуется учет типа флюидов, заполняющих трещины. Если в залежи сочетаются различные типы коллекторов (например, межзерновой и трещинный), то подсчет производится отдельно для каждого типа и затем суммируется. При этом точность решения обратной задачи по данным ВСП ограничена отсутствием информации о геометрии включений, физических свойств их материала и характера взаимодействия включений между собой.
Наиболее близким аналогом предлагаемого изобретения является способ определения трещинной пористости пород (патент RU 2012021, 1994), в котором определение трещинной пористости выполняют по данным волнового акустического и гамма-гамма каротажа по результатам измерений при двух значениях плотности заполняющей скважину жидкости.
Этот способ не может быть использован при скважинной сейсморазведке, так как расстояния между источниками и приемниками в скважинной сейсморазведке на несколько порядков больше, чем в каротаже (в каротаже - доли метра, в скважинной сейсморазведке - тысячи метров), а частоты сигналов, соответственно, ниже. В результате разных масштабов измерений в скважинной сейсморазведке свойства жидкости, заполняющей скважину, не оказывают влияний на результат измерений. Кроме того, для реализации известного способа применяется трудоемкий процесс замены скважинной жидкости. На практике известный способ может быть использован лишь в небольших объемах для оценки пористости в отдельных интервалах в отдельных скважинах.
Задачей настоящего изобретения является повышение информативности и экспрессности анализа трещинной пористости пласта непосредственно по данным многолучевого НВСП.
Поставленная задача решается тем, что согласно Способу оценки трещинной пористости пород по данным скважинной сейсморазведки, заключающемуся в проведении в изучаемом разрезе скважины волнового акустического и плотностного гамма-гамма каротажа, выявлении на основе полученных данных интервала трещиноватых пород и вычислении коэффициента трещинной пористости выявленного интервала трещиноватых пород, согласно изобретению на основе полученных данных дополнительно проводят скважинную сейсморазведку, для чего в выявленном интервале ориентированной трещиноватости один из источников излучения сейсмических волн ориентируют в направлении, близком к азимутальному направлению трещиноватости, а другой источник излучения сейсмических волн ориентируют ортогонально направлению трещиноватости, затем по анализу изменения скоростей Vp, Vs от источников излучения сейсмических волн и их отношения Vs/Vp осуществляют выделение интервалов пористых коллекторов с применением выражения
2
2
βэф=3/ρ(3Vp-4Vs),
где βэф - эффективное значение коэффициента сжимаемости; ρ - плотность,
и производят оценку трещинной пористости коллекторов порово-трещинного или порово-кавернозно-трещинного типа по коэффициенту трещинной пористости Кпт, определяемому:
а) - при хаотичной трещиноватости - по формуле:
Кпт=(βэф-βм)/βт,
где βэф - эффективное значение коэффициента сжимаемости; βм - коэффициент сжимаемости пористой/порово-кавернозной матрицы, определяемый в нетрещиноватом интервале; βт - коэффициент сжимаемости трещин, определяемый на участке с известной трещинной пористостью по формуле:
βт=(βэф-βм)Кпт;
б) - при ориентированной трещиноватости - по формуле:
a
2
Кпк=(σэф/σм-1)/2sinφ,
2
2
2
2
a
где σэф=(0.5Vp/Vs-1)/(Vp/Vs-1) - эффективный коэффициент Пуассона, определяемый по данным источника излучения сейсмических волн, ориентированного ортогонально направлению трещиноватости; σм - коэффициент Пуассона пористой/порово-кавернозной матрицы, определяемый по данным источника излучения сейсмических волн, ориентированного в направлении, наиболее близком к азимутальному направлению трещиноватости;  - среднее значение отношения раскрытости трещин к их длине; φ - угол между направлением трещин и направлением распространения волны источника излучения сейсмических волн, использованного при определении σэф.
Как известно, одним из критериев трещиноватых коллекторов служит значительное понижение скорости Vp, определяемой по данным ВСП относительно вычисленной по среднему времени с использованием коэффициента пористости Кпт, определяемого по НГК или ГГК.
Поскольку каверны в связи с их низкой сжимаемостью существенно меньше влияют на скорости Vp и Vs, трещиноватые коллекторы даже при небольшой трещиноватости породы выделяются повышенными значениями сжимаемости βэф.
Предложенный способ существенно отличается от известных тем, что обеспечивает быструю оценку трещинной пористости по аналитическим выражениям с использованием преимущественно данных ВСП и не связан с необходимостью определения физических свойств материала коллекторов и характера взаимодействия включений между собой, а также с необходимостью замены скважинной жидкости.
На практике предложенный способ оценки трещинной пористости пород по данным скважинной сейсморазведки реализуется следующим образом.
В изучаемом разрезе пласта предварительно выполняют измерения в нетрещиноватом интервале и в интервале с известной трещинной пористостью породы. Эти измерения могут быть проведены любым известным способом промысловой геофизики (например, по данным волнового акустического и гамма-гамма каротажа, как у прототипа). При этом, при ориентированной трещиноватости породы, один из источников излучения сейсмических волн ориентируют в направлении, наиболее близком к азимутальному направлению трещиноватости породы, а другой источник излучения сейсмических волн ориентируют ортогонально направлению трещиноватости породы. В процессе анализа изменения скоростей продольных Vp и поперечных Vs волн, определяемых по данным многолучевого НВСП и изменения параметра отношения этих скоростей Vs/Vp, определяют эффективное значение коэффициента сжимаемости βэф с использованием выражения
по которому выявляют интервалы пористых коллекторов.
По понижению скоростей продольных Vp и поперечных Vs волн и их отношения Vs/Vp, а также по резкому повышению сжимаемости производят оценку трещинной пористости коллекторов и рассчитывают коэффициент трещинной пористости Кпт.
При хаотичной трещиноватости коэффициент трещинной пористости Кпт коллекторов порово-трещинного или порово-кавернозно-трещинного типа определяют по формуле:
где βм - коэффициент сжимаемости пористой/порово-кавернозной матрицы;
βт - коэффициент сжимаемости трещин.
Коэффициент сжимаемости пористой/порово-кавернозной матрицы βм определяют по коэффициенту сжимаемости βэф в нетрещиноватом интервале. Наличие каверновой составляющей пористости определяют по повышению скорости Vp, вычисленной по данным ВСП относительно скорости, вычисленной по среднему времени с использованием коэффициента пористости Кпт, определяемого по НГК или ГГК. Коэффициент сжимаемости трещин βт на участке с известной трещинной пористостью Кпт определяют по формуле: βт=(βэф-βм)/Кпт.
При ориентированной трещиноватости
- в выделенных трещиноватых интервалах определяют преобладающее направление субвертикальной трещиноватости. По данным НВСП от источника излучения сейсмических волн, ориентированного в направлении, близком к азимутальному направлению трещиноватости, определяют скорости продольной волны Vp и быстрой поперечной волны Vs||. С использованием полученных значений указанных скоростей и значения плотности ρ (по данным плотностного каротажа) по выражению (1) и выражению (3)
рассчитывают эффективное значение коэффициента сжимаемости βэф и коэффициента Пуассона σэф.
Поскольку в выбранной системе наблюдений скорости Vp и Vs|| практически не зависят от трещиноватости породы, определенные коэффициент эффективной сжимаемости βэф и коэффициент Пуассона σэф соответствуют значениям пористой/пористо-кавернозной матрицы βм и σм.
Аналогичным образом определяют эффективное значение коэффициента сжимаемости βэф и коэффициента Пуассона σэф по данным НВСП от источника излучения сейсмических волн, ориентированного в направлении, наиболее отличающемся от направления трещиноватости (ортогонально направлению трещиноватости).
Коэффициент трещинной пористости Кпк породы в этом случае определяют с использованием полученных значений σм и σэф по формуле:
a
где  - среднее значение отношения раскрытости трещин к их длине; φ - угол между направлением трещин и направлением распространения волны от источника излучения сейсмических волн, использованного при определении σэф.
a
Для повышения точности результата аспектное отношение  определяют по экспериментальным данным в близких условиях на участках с известной трещинной пористостью по формуле (4) или с использованием выражения (5)
Таким образом, предложенный способ позволяет выполнить оперативную оценку трещинной пористости породы за счет использования эффективных значений упругих параметров излученных волн и избежать трудоемкую операцию математического моделирования с определением физических свойств материала коллекторов и характера взаимодействия включений между собой.