СПОСОБ ОПРЕДЕЛЕНИЯ ОТНОСИТЕЛЬНЫХ ФАЗОВЫХ ПРОНИЦАЕМОСТЕЙ ПЛАСТА
METHOD FOR DETERMINING RELATIVE PHASE PERMEABILITIES OF FORMATION

Изобретение относится к технологиям нефтедобычи, а именно к способам гидродинамического моделирования залежей и проектирования на их основе разработки месторождений.
Наиболее важными критериями при построении и настройке гидродинамической модели пласта являются экспериментальные кривые относительной фазовой проницаемости (ОФП), обычно получаемые путем стендовых лабораторных исследований на керне (см., например, изобретение по патенту РФ №2097740, 27.11.1997) или определяемые по аналогии с другими схожими месторождениями.
Несмотря на то что способ определения ОФП на кернах является прямым методом исследований, однако он характеризуется микроуровнем оценки параметров пласта и позволяет исследовать лишь ничтожную часть от реального объема нефтяного резервуара.
Более перспективным является способ, основанный на гидродинамических испытаниях нефтяных необводненных скважин, при которых оценивают интегральные свойства пласта в межскважинном пространстве, что соответствует оценке ОФП на макроуровне.
Такой способ, в частности, описан в работе: Закиров С.Н., Индрупский И.М., др. Новый подход к исследованиям скважин и пластов. Нефтяное хозяйство, 2002, №6, с.113-115.
В соответствии с описанным способом предложено проводить несколько циклов гидродинамических исследований скважин (ГДИС) на ранней стадии разработки путем принудительного создания в пласте разнонаправленных двухфазных фильтраций при различной текущей нефтенасыщенности.
но
н
во
во
н
но
Исследования по этому способу проводят в следующей последовательности: сначала осуществляют долговременный режим стабильного отбора нефтью, после этого скважину останавливают и проводят цикл ГДИС по технологии регистрации кривой восстановления давления (КВД) (это исследование позволяет определить проницаемость пласта k при фильтрации чистой нефти, при насыщении остаточной водой K=K), затем производят массированную закачку воды в пласт, после чего скважину останавливают и проводят цикл ГДИС по технологии регистрации кривой падения давления (КПД) (это исследование позволяет определить проницаемость пласта k при фильтрации воды, при насыщении остаточной нефтью K=K), далее производят долговременную отработку пласта, при которой постепенно, по мере очистки пласта, растет доля нефти в продукции.
в
н
В процессе отработки осуществляют несколько остановок скважины и регистрацию циклов КВД. Каждой остановке соответствуют свои значения фазовых проницаемостей по воде k и по нефти k, соответствующие конкретной нефтенасыщенности на этот момент.
во
но
но
н
во
Данная серия исследований позволяет определить фазовые проницаемости, соответствующие конкретной текущей насыщенности, то есть получают значения фазовой проницаемости для всех характерных точек ОФП: краевых, соответствующих остаточной водонасыщенности K и остаточной нефтенасыщенности K, и промежуточных, когда нефтенасыщенность находится в диапазоне K<K<1-K. При этом, если для каждого из циклов ГДИС известны текущие насыщенности пласта, то полученные результаты полностью характеризуют кривые ОФП (зависимости текущей фазовой проницаемости от насыщенности).
В соответствии с этим способом определяют зависимость ОФП от водонасыщенности пласта (по данным КВД) математическим путем.
Недостатком данного способа является отсутствие критериев для распространения полученных ОФП на соседние скважины, отдельные зоны пласта или пласт в целом.
Задачей изобретения является повышение надежности и объективности воспроизведения ОФП путем обеспечения возможности распространения полученных для конкретной скважины ОФП на другие участки пласта.
Поставленная задача решается следующим образом.
в_отн
н_отн
Выделяют участки (блоки) пласта, с различной литолого-фациальной характеристикой, в пределах каждого блока определяют обобщенные ОФП, для чего выполняют следующие операции: исходя из критерия отсутствия интервалов заколонных перетоков, негерметичностей обсадной колонны, подтягивания пластовой воды, а также прорывов нагнетаемых вод, определяют представительную опорную выборку скважин, далее в каждой из скважин опорной выборки проводят гидродинамические исследования, по результатам которых оценивают текущие фазовые проницаемости по воде k и по нефти k.
При этом в начальный период эксплуатации (когда водонасыщенность пластов и обводненность продукции скважины еще не претерпели существенного изменения) выполняют следующие операции:
в_нач
в_тек
в_нач
1) в каждой из скважин опорной выборки оценивают начальную водонасыщенность K по ГИС открытого ствола, которая в рассматриваемый период эксплуатации соответствует текущей водонасыщенности K=K;
н_отн
в_нач
2) строят кросс-плот, на который наносят точки k, K (количество точек на кросс-плоте соответствует числу скважин), точки на кросс-плоте аппроксимируются, например, известной степенной зависимостью:
в_отн
в_нач
3) строят кросс-плот, на который наносят точки k, K (количество точек на кросс-плоте соответствует числу скважин), точки на кросс-плоте аппроксимируются, например, известной степенной зависимостью:
но
во
выт
н
в
4) методом наименьших квадратов, с учетом известных по данным петрофизических исследований значений остаточной нефте- и водонасыщенности K, K, и коэффициента вытеснения β определяют коэффициенты зависимостей (1) и (2) α, α и β, по которым рассчитывают кривые фазовых проницаемостей.
На стадии эксплуатации, характеризуемой изменением средней обводненности продукции по скважинам не менее 30% по сравнению с первоначальной (то есть при снижении нефтенасыщенности пласта вследствие его выработки), выполняют следующие операции:
1) определяют представительную опорную выборку скважин, исходя из критериев:
в_тек
в_тек
а) отсутствия интервалов заколонных перетоков, негерметичностей обсадной колонны, интервалов подтягивания пластовой воды и прорывов нагнетаемых вод (при выполнении этого критерия значение проницаемости по гидродинамическим исследованиям характеризует исключительно перфорированный пласт, а текущая обводненность продукции φ меняется исключительно за счет изменения водонасыщености пласта К);
в_тек
в_тек
в_тек
б) нахождения текущей обводненности в диапазоне от 0.3<φ<0.9 (этот критерий согласно теории Баклея-Леверетта определяет диапазон обводненности продукции, при которой параметры обводненности φ и К наиболее тесно увязаны друг с другом);
в_тек
2) в каждой из скважин опорной выборки оценивают текущую обводненность продукции φ;
н_отн
в_тек
3) строят кросс-плот, на который наносят точки k, φ (количество точек на кросс-плоте соответствует числу скважин), точки на кросс-плоте аппроксимируют, например, известной степенной зависимостью:
в_отн
в
4) строят кросс-плот, на один из которых наносят точки k, φ (количество точек на кросс-плоте соответствует числу скважин), точки на кросс-плоте аппроксимируют, например, известной степенной зависимостью:
но
во
выт
н
в
5) методом наименьших квадратов, с учетом известных по данным петрофизических исследований значений остаточной нефте- и водонасыщенности К, К и коэффициента вытеснения β определяют коэффициенты зависимостей (3) и (4) α, α и β, по которым рассчитывают условные кривые фазовых проницаемостей.
в_тек
в_тек
Условные кривые фазовых проницаемостей при выполнении вышеуказанных критериев соответствуют реальным, поскольку в этом случае текущая насыщенность пласта К и обводненность продукции φ согласно теории Баклея-Леверетта изменяются пропорционально друг другу.
Таким образом, обеспечивается возможность распространения полученных для конкретной скважины ОФП на другие участки пласта.
Эффективность предлагаемого метода показали комплексные испытания, проведенные в ОАО «Газпромнефть».
Практический пример реализации способа показан на фиг.1 и 2.
Успешности решения задачи построения обобщенных ОФП способствовало то, что рассматриваемый участок залежи характеризовался малой фациальной изменчивостью и отсутствием существенных макронеоднородностей.
н_отн
в_нач
в_отн
в_нач
Кроме того, на начальном этапе эксплуатации по данным ГИС обнаружено существенное изменение начальной нефтенасыщенности по простиранию пласта. Это позволило получить точки на кросс-плотах k, K и k, K практически во всем диапазоне возможного изменения насыщенности, то есть получить представительные кривые ОФП.
Блок залежи был практически полностью охвачен как гидродинамическими, так и геофизическими исследованиями, реализованными на основе специальной обработки результатов качественных долговременных измерений датчиками на приеме насоса, что позволило определить фильтрационные свойства пласта с высокой точностью.
Значения проницаемости по ГДИС в условиях двухфазной фильтрации были получены по опорной выборке из 34 скважин рассматриваемого участка месторождения.
ΣГДИС
в_нач
ΣГДИС
Σ
в_нач
нГДИС
н
в_нач
вГДИС
в
в_нач
В результате обобщения данных о проницаемости, полученных при интерпретации ГДИС на начальном этапе эксплуатации по всем скважинам выбранного участка К, и насыщенности K, оцененной по ГИС, были определены закономерности изменения от насыщенности суммарной проницаемости пласта по ГДИС k=k(K) фазовой проницаемости по нефти k=k(K) и по воде k=k(K), которые могут быть пересчитаны в кривые относительных фазовых проницаемостей.
ΣГДИС
Σ
в_нач
Зависимость k=k(K) показана на фиг.1.
ΣГДИС
в_нач
-2.38
Факт, что реальная связь проницаемости от насыщенности получилась достаточно тесной (коэффициент корреляции более 0.85) подтверждает правомерность предположения о том, что исследуемый объект можно приближенно описать едиными кривыми ОФП. Осредненную суммарную функцию фазовых проницаемостей можно получить осреднением полученных данных аналитической зависимостью К=3400·К . Локальные отклонения от данной зависимости связаны с изменением абсолютной проницаемости пласта по простиранию.
Для иллюстрации возможности реализации предлагаемого способа для случая, когда текущая насыщенность отличается от первоначальной и не может быть оценена по результатам ГИС в открытом стволе, были выбраны скважины, где диапазон изменения обводненности находился между 0.3 и 0.9. Для этих скважин был построен кросс-плот, связывающий величины водонасыщенности по ГИС и обводненности по промысловым данным (фиг.2).
Точки на данном кросс-плоте фактически иллюстрируют поведение функции Баклея-Леверетта для рассматриваемого участка месторождения. На данный кросс-плот нанесены также зависимости Баклея-Леверетта, рассчитанные по формуле:
(тонкие линии на фиг.2). Расчеты были выполнены для определенных значений фазовых проницаемостей при различных соотношениях вязкостей воды и нефти.
Это свидетельствует, что в рассматриваемых условиях обводненность пласта с достаточной степенью точности характеризует его водонасыщенность. То есть перечисленные параметры, при приближенной оценке ОФП, могут быть взаимозаменяемы.