INFLUENCE OF LAYERED NONUNIFORMITY OF THE RESERVOIR ON THE PARAMETERS OF EXTRACTION AT NON-ISOTHERMAL WATER-OIL DISPLACEMENT
Rubrics: TECHNICAL SCIENCES
Authors:
1.
Dagestan State University
(Doctor of Technical Sciences, Professor; Professor of the Department of Information Technology and Security of Computer Systems)
2.
Dagestan State University
(Candidate of Physics and Mathematics, Assistant Professor; Head of the Department of Information Technologies and Security of Computer Systems)
Russian Federation
Russian Federation
3.
Astrakhan State Technical University
Type:
Article
Pages:
from 14
to 21
Status:
Published
Received:
10.04.2016
Accepted:
15.04.2016
Published:
15.04.2016
Subject area:
UDK 622.276.031:53
GRNTI 20.01 Общие вопросы информатики
GRNTI 26.03 Общественно-политическая мысль
GRNTI 43.01 Общие вопросы естественных и точных наук
GRNTI 44.01 Общие вопросы энергетики
GRNTI 45.01 Общие вопросы электротехники
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 62.01 Общие вопросы биотехнологии
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 70.01 Общие вопросы водного хозяйства
GRNTI 73.34 Водный транспорт
GRNTI 20.01 Общие вопросы информатики
GRNTI 26.03 Общественно-политическая мысль
GRNTI 43.01 Общие вопросы естественных и точных наук
GRNTI 44.01 Общие вопросы энергетики
GRNTI 45.01 Общие вопросы электротехники
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 62.01 Общие вопросы биотехнологии
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 70.01 Общие вопросы водного хозяйства
GRNTI 73.34 Водный транспорт
Language:
Russian
Keywords:
reservoir, filtering, paraffin oil, pressure, saturation, indicators, development
Abstract (English):
Improvement and creation of essentially new highly effective technologies of oil refining require the deep analysis of the complex mechanisms of intra sheeted processes. Influence on the results of production of the layers at non-isothermal filtration of layered heterogeneity of the collectors on permeability is shown. Many oil fields typically contain paraffinic and asphaltene components. At the fields of high-wax oils with a saturation temperature of the paraffin close to the reservoir temperature, a small reduction can lead to crystallization of paraffin, cooling and plugging of steam channels. At the fields of highly viscous and waxy crude oil the proportion of low permeability formations, exposed primarily to temperature effects, i. e., cooling by injecting cold water in the range from 0.01 to 0.08 µm2, sometimes up to 40 %.
Improvement and creation of essentially new highly effective technologies of oil refining require the deep analysis of the complex mechanisms of intra sheeted processes. Influence on the results of production of the layers at non-isothermal filtration of layered heterogeneity of the collectors on permeability is shown. Many oil fields typically contain paraffinic and asphaltene components. At the fields of high-wax oils with a saturation temperature of the paraffin close to the reservoir temperature, a small reduction can lead to crystallization of paraffin, cooling and plugging of steam channels. At the fields of highly viscous and waxy crude oil the proportion of low permeability formations, exposed primarily to temperature effects, i. e., cooling by injecting cold water in the range from 0.01 to 0.08 µm2, sometimes up to 40 %.
Keywords:
reservoir, filtering, paraffin oil, pressure, saturation, indicators, development
reservoir, filtering, paraffin oil, pressure, saturation, indicators, development
Введение Нефти многих месторождений обычно содержат парафинистые и асфальтосмолистые компоненты. На месторождениях высокопарафинистых нефтей с температурой насыщения парафином, близкой к пластовой, небольшое ее снижение может привести к кристаллизации парафина, его охлаждению и закупорке паровых каналов. Целью исследования в этом случае становится существенно важное реальное строение пласта. Постановка задачи. Пласт предполагается слоистым, число пропластков и их расположение произвольно, пропластки имеют различную проницаемость. Они считаются гидродинамически изолированными (т. е. нет перетоков), но термически контактирующими. В каждом из пропластков вытеснение происходит автономно. Распределение давлений, насыщенностей и температур в пропластках различно. Система расположения нагнетательных и эксплуатационных скважин предполагается двоякопериодической, имеющей прямоугольный элемент симметрии. Каждая из скважин вскрывает все пропластки многослойного пласта, а давление на забое для всех пропластков поддерживается одинаковым. В нагнетательных и добывающих скважинах поддерживаются фиксированные значения давления, которые могут быть изменены в выбранные моменты времени на новые значения. Добывающая скважина может быть закрыта по определённому пропластку, если этот пропласток в данной скважине имеет высокую обводнённость. Температура нагнетаемой воды отлична от температуры пластовой и может менять свои значения в определённые моменты времени или после прокачивания определённого объёма воды. Теплоперенос в пропластках осуществляется за счёт конвекции тепла с фильтрующейся жидкостью и теплообмена с окружающими пропластками через глинистые перемычки. Крайние пропластки обмениваются теплом с кровлей и подошвой. При постановке задачи приняты некоторые ограничительные предположения, широко применяемые в теоретических исследованиях [1, 2]. Метод исследования В ходе исследования принимается схема Ловерье - схема учёта обмена теплом с кровлей и подошвой (глинистые перемычки, пропластки и окружающая порода считаются теплопроводными лишь в направлении, перпендикулярном напластованию). Низкопроницаемые пропластки при их охлаждении могут исключаться из разработки в двух случаях: если перепад давлений мал и недостаточен для преодоления градиента сдвига и если вокруг нагнетательной скважины пласт остудился до температуры застывания нефти до прохождения фронта вытеснения. Для высокопарафинистых нефтей закон фильтрации отличается от закона Дарси, и в данном случае скорость фильтрации фаз выражается следующим образом: где kв и kн - относительная проницаемость воды и нефти соответственно; x - структурный множитель, который при пластовой температуре близок к единице, а со снижением температуры становится существенно меньше единицы. Относительная фазовая проницаемость нефти и воды была получена с учётом возможного выпадения парафина при охлаждении пласта [3, 4]: ; ; ; . При T ³ T0 принимают R = 0; при T £ Tкр (Tкр - температура застывания нефти в пласте) нефть застывает и kн = 0. При пониженных значениях температуры, вплоть до температуры застывания, R вычисляется по формуле . Для определения вязкости воды и нефти были приняты известные аналитические выражения: , где с - концентрация соли в воде, % Рассмотрим двумерную задачу неизотермического вытеснения нефти водой в тонком, однородном по мощности пласте. В условиях несжимаемости фаз и отсутствия их теплового расширения для законов сохранения масс воды и нефти при их фильтрации сквозь пористый коллектор имеем следующие выражения [2]: ; (1) Принято пользоваться суммой уравнений (1) для определения поля давлений: . По известному полю давлений находим поле скоростей, а затем и приращения насыщенностей, пользуясь упрощённым уравнением переноса массы воды: . Уравнение переноса тепла в каждом из пропластков записываем в следующем дивергентном виде: , где qk - тепловой поток, уходящий из пропластка с номером k; Hk - мощность пропластка; ck - теплоёмкость пропластка с присоединённой половиной толщины hk слоя глинистой перемычки. Удельная теплоёмкость на единицу площади k-го пропластка вычисляется с учётом присоединения глинистых перемычек: . Тепловое сопротивление между пропластками больше, чем тепловое сопротивление глинистой перемычки, т. к. под температурой пропластка подразумевается средняя по её мощности температура, а при линейном распределении значений температуры, например, станет необходимым добавлять ещё и половины тепловых сопротивлений соседних пропластков: . К этим уравнениям присоединяются также уравнения теплопроводности для окружающих горных пород: , (2) где координата z перпендикулярна напластованию, а под c подразумевается объёмная теплоёмкость кровли или подошвы. Уравнение (2) означает, что для теплообмена с кровлей и подошвой принята гипотеза Ловерье, т. е. не учитывается теплопроводность по напластованию пластов. Уравнение переноса тепла с учётом конвективной и кондуктивной составляющих потока тепла для элемента пористой среды имеет вид , (3) где с - объёмная теплоёмкость пористой среды; св и сн - объёмная теплоёмкость воды и нефти соответственно; lэ - эффективная кондуктивная теплопроводность пористой среды с учётом ее неоднородности и движения жидкости. Для больших периодов времени, какими являются периоды разработки пластов, кондуктивная теплопроводность не играет существенной роли, поэтому правую часть выражения (3) часто заменяют нулём. Кроме того, фронт тепла намного отстаёт от фронта вытеснения, изменение температуры происходит в зоне, где скорость фильтрации нефти настолько мала, что ею можно пренебречь. Можно пользоваться уравнением теплопереноса вида [2, 3]: . Объёмная теплоёмкость пласта определяется как сумма теплоёмкостей фаз и скелета породы: . Семиточечная система площадного заводнения получается в том случае, если каждый отдельный ряд нагнетательных скважин чередуется с двумя рядами добывающих скважин. Элементом симметрии данной системы заводнения является правильный шестиугольник, в вершинах которого расположены добывающие скважины, а в центре - нагнетательная (что и отразилось на названии системы). В трёхрядной системе расстановки скважин одной нагнетательной скважине отнесены две эксплуатационные скважины с дебитами q2 и q3, если дебит нагнетательной скважины q1 = q2 + q3. Для двумерных задач вытеснения выражение для давления удовлетворяет эллиптическому уравнению [2]: (4) где k - проницаемость; H - мощность пласта; fв и fн - относительная проницаемость воды и нефти, зависящая от водонасыщенности; mв и mн - вязкость воды и нефти. При разностной аппроксимации уравнения (4) необходимо учитывать наличие особенностей в скважинах, которые заменяются источниками и стоками. Для узла (i, j) разностной сетки, когда соседние с ним узлы не являются скважинами, условие баланса запишется в следующем виде [2, 4]: . (5) При , используя формулы (4), имеем для определения давления в обычном узле: (6) Если в узле (i, j + 1) оказывается добывающая скважина, то вместо (5) имеем: , где qсек выражается через разность давлений и гидропроводность. В кольцевой зоне возле нагнетательной скважины суммарную скорость фильтрующихся фаз представим в виде . (7) Эта формула справедлива для любого пропластка. Из выражения (7) находим градиент давления. Дебит скважины на единицу мощности q равен: , где рнаг соответствует давлению в узле (1,1), т. е. за начало координат принята точка с индексами (1,1). Для равномерной сетки концентрических колец фильтрационное сопротивление вычисляется по формуле . Для добывающей скважины, окрестность которой разбита на четыре сектора, дебит определяется суммированием расходов жидкости по каждому сектору: где Rс - радиус скважины; - расстояние до диагональных соседних узловых точек. Вблизи добывающих скважин температура и насыщенность меняются достаточно плавно, и их значения в секторах можно принять за постоянные, а потерянный перепад можно совсем не учитывать: . При реализации расчётов в программе для определения поля давлений в узлах сетки (6) применялся метод простых итераций с поточечной релаксацией с коэффициентом 1,75. Вместе с давлениями итерировались и значения дебитов нагнетательных скважин. Для определения насыщенности ячеек, удалённых от скважин, используется уравнение баланса воды, которое записывается в дивергентной форме: , где Fв - доля воды в потоке. Термогидродинамические расчёты, основанные на моделях 2- и 3-слойных пластов, могут искажать картину фильтрации и приводить к существенным ошибкам при расчёте показателей разработки парафинистых нефтей [5]. Результаты исследования Результаты расчётов были апробированы для элемента семиточечной системы заводнения многослойного пласта (до 10 пропластков) при фиксированных значениях большинства определяющих параметров. Теплофизические параметры брались близкими к типичным параметрам месторождений высокопарафинистых нефтей бывшего Советского Союза (Ромашкинское, Узеньское, Астраханское и т. д.): температура застывания нефти - 40 °С; пластовая температура - 65 °С; температура нагнетаемой воды - 15, 65 и 90 °С. Рассмотрим пласт эффективной толщины 19 м, состоящий из 10 пропластков мощностью, м: 2,4; 2,5; 3,4; 0,6; 1,5; 1,5; 2,0; 1,5; 1,5; 2,1 и проницаемостью, мкм2: 0,36; 0,01, 0,45; 0,04; 0,25; 0,2; 0,05; 0,1; 0,08; 0,03. Расстояние между рядами равно 250 м, а между скважинами в ряду - 600 м. Элемент симметрии 500 ´ 300 м (рис. 1, 2). При холодном вытеснении 2-й пропласток застыл в течение первого года разработки, 4-й пропласток застыл к 18-му году. Нефтеотдача пласта при обводненности добывающей скважины 60,1 % составляет 24,1 %. При создании горячей оторочки воды с температурой 90 °С на начальном этапе разработки в течение 2-х лет, нефтеотдача к 35-му году составляет 28,8 % при обводнённости 77,7 %. При нагнетании горячей воды с температурой равной пластовой, т. е. 65 °С, нефтеотдача к 35-му году составила 31,9 % при обводнённости 78,6 %. Рис. 1. Зависимость нефтеотдачи от времени разработки для 10-слойного пласта при различных значениях температуры закачиваемой воды (TVC) Рис. 2. Зависимость общей обводнённости пласта от времени разработки для 10-слойного пласта при различных значениях температуры закачиваемой воды Таким образом, согласно расчетам и данным графиков, нагнетание горячей воды в пласт ведет к отсутствию застывших и закупорившихся за время разработки пропластков, увеличению нефтеотдачи и обводнённости разрабатываемого пласта, что обеспечивает улучшение общих показателей разработки месторождения. Выводы Разработка нефтяных месторождений, содержащих вязкие и парафинистые нефти, с применением тепловых методов воздействия на пласты связана с существенными затратами. Хотя результаты разработки при поддержании пластовой температуры выше, чем при создании горячей оторочки, добыча нефти термическими методами должна строиться с учётом организации энергозатрат для конкретного месторождения. Из приведённых расчётов и графиков зависимости нефтеотдачи и обводнённости для 10-слойного пласта при температуре 15, 65, 90 °С в течение первых 2-х лет и 15 °С далее видно, что нагнетание горячей воды в пласт ведёт как к отсутствию застывших и закупорившихся за время разработки пропластков, так и к увеличению нефтеотдачи и обводнённости разрабатываемого пласта и, как следствие, к улучшению общих показателей разработки месторождения.
1. Aziz H. Matematicheskoe modelirovanie plastovyh sistem / H. Aziz, E. Settari. M.: Nedra, 1982. 430 s.
2. Alishaev M. G. Konechno-raznostnyy raschet neizotermicheskogo vytesneniya nefti vodoy v mnogosloynom plaste / M. G. Alishaev, S. A. Ahmedov, G. O. Hazamov // Mezhvuz. nauch.-tehn. seminar «Funkcional'nyy analiz: teoriya funkciy i ih prilozheniya». Mahachkala, 1985. S. 42-49.
3. Alishaev M. G. Raschet temperaturnogo polya plasta pri inzhekcii zhidkosti dlya ploskogo fil'tracionnogo techeniya / M. G. Alishaev // Izv. AN SSSR. Mehanika zhidkosti i gaza. 1980. № 1. S. 67-75.
4. Alishaev M. G. Neizotermicheskoe vytesnenie parafinistoy nefti pri trehryadnoy sisteme razmescheniya skvazhin s uchetom mezhsloynogo obmena / M. G. Alishaev, S. A. Ahmedov // Neftyanoe hozyaystvo. 1998. № 11. S. 31-32.
5. Ahmedov S. A. Komp'yuternoe modelirovanie processa neizotermicheskogo vytesneniya vyazkoplastichnoy nefti v mnogosloynom plaste / S. A. Ahmedov, R. M. Alisultanov // Neftyanoe hozyaystvo. 1999. № 12. S. 38-41.
