Full Text

Введение В настоящее время нефтедобывающая промышленность сталкивается с проблемой активного истощения большинства месторождений традиционного углеводородного сырья. Поддержание и увеличение достигнутых темпов добычи может быть осуществлено путем интенсивного вовлечения в разработку запасов тяжелой нефти и природных битумов, потенциальные ресурсы которых превышают традиционные залежи [1-4]. Основной особенностью данного вида углеводородов является их малая подвижность, обусловленная высокой вязкостью. По этой причине высокие показатели нефтеотдачи могут быть достигнуты только путем проведения специальных геолого-технических мероприятий [5, 6]. Одной из наиболее эффективных технологий разработки залежей сверхвязких нефтей является парогравитационное дренирование (ПГД) [7, 8]. Данный метод предполагает бурение двух параллельных горизонтальных скважин, расположенных друг над другом. Верхняя скважина используется для нагнетания перегретого пара в продуктивный пласт с целью создания и развития паровой камеры с температурой выше температуры парообразования. По мере повышения температуры в пласте подвижность нефти увеличивается, и она под действием гравитации стекает к нижней добывающей скважине вместе с пластовой водой и конденсированным паром. Рост паровой камеры обусловлен повышенным поровым давлением и теплообменом между продуктивным пластом и закачиваемым горячим паром, включающим в себя процессы теплопроводности, конвективного теплопереноса и выделения скрытой теплоты фазового перехода. Происходящее при этом изменение напряженно-деформированного состояния продуктивного пласта может приводить к разрушению кровли [1, 9-11], срезу обсадных колон нагнетательной и добывающей скважин [12] и изменению его фильтрационно-емкостных характеристик. В работе [1] отмечается, что при высоком давлении нагнетания может возникнуть вертикальный разлом, распространение которого приводит к нарушению целостности кровли и прорыву пара в верхние слои. Другим геомеханическим эффектом является сжатие пористого скелета кровли в результате теплового воздействия, что грозит падением горизонтального напряжения и как следствие возникновением вертикального гидроразрыва. В работах [9-11] проводится численное моделирование процесса ПГД с целью исследования прочности кровли продуктивного пласта на основе критериев Мора-Кулона и Друкера-Прагера с учетом анизотропии механических свойств. В работе [13] на основе термопороупругости показано, что воздействие высоких температур приводит к увеличению сжимающих напряжений, что может вызывать разрушение обсадных колон скважин. Исследованию влияния на устойчивость обсадных колон сжимающих и сдвиговых напряжений, возникающих при повышении давления и температуры, посвящены работы [14, 15]. В работе [16] отмечается, что одной из основных причин возникновения повреждений обсадных колон в ходе ПГД является ползучесть околоскважинного пространства пласта. Обширные экспериментальные исследования показывают, что существенное влияние на прочностные и деформационные характеристики горной породы оказывает наличие в ней пор и характер их распределения [17, 18]. В ходе технологического процесса ПГД изменение пористости происходит вследствие различных механизмов, связанных с увеличением температуры и порового давления. В работе [19] на основе данных мониторинга процесса ПГД, проведенного с использованием радиолокационной интерферометрии, было установлено, что при высоких давлениях нагнетания пара возникает сдвиговая дилатансия, которая существенно увеличивает как пористость, так и проницаемость продуктивного пласта. В работах [20-22] этот эффект был исследован на основе численного моделирования. В свою очередь, в работе [12] отмечается, что при высоком поровом давлении нагрев может привести к уменьшению пористости и проницаемости вследствие расширения зерен породного скелета. В работе [23] с использованием численного моделирования показано, что отказ от учета геомеханических эффектов приводит к некорректной оценке нефтеотдачи. Для описания нефтедобычи методом ПГД используются аналитические и численные модели. В работе [24] на основании закона Дарси и законов сохранения массы и энергии было получено одно из первых аналитических соотношений для расчета дебита добывающей скважины. Дальнейшее развитие этого подхода связано с учетом дополнительных факторов, таких как форма паровой камеры и отток тепла в кровлю продуктивного пласта [25, 26], нестационарность теплообмена на границе паровой камеры [27], предельного градиента давления, достижение которого необходимо для течения сверхвязкой нефти [28]. Несмотря на значительный прогресс в развитии аналитических моделей, они не позволяют в полной мере описать конвективный теплоперенос и фазовый переход, возникающий при конденсации пара и испарении воды, а также не учитывают геомеханические эффекты. По этой причине для описания сложного взаимодействия физических процессов при ПГД широкое распространение получило численное моделирование. Большинство моделей нефтедобычи с использованием паротеплового воздействия включают в себя закон сохранения массы для трехфазного потока, закон сохранения энергии и уравнение равновесия [29-31]. Влияние порового давления и температуры на механические напряжения описывается с использованием теорий Био и термоупругости. В то же время существует большое разнообразие в выборе определяющих соотношений для описания изменения пористости [32]. В наиболее простой модели пористость линейно зависит от порового давления с коэффициентом пропорциональности, равным коэффициенту поровой сжимаемости [21, 33]. В [34, 35] с использованием теории гидростатической пороупругости в предположении изменения объемов породного скелета и порового пространства была получена зависимость пористости от средних напряжений в приращениях с коэффициентом, линейно зависящим от текущего значения пористости. При малой сжимаемости породного скелета пористость может быть связана со средними эффективными напряжениями через композицию дробно-линейной функции и экспоненты [36]. В работе [37] на основе лабораторных экспериментов по сжатию песчаника была предложена другая экспоненциальная зависимость пористости от средних эффективных напряжений. В [38, 39] на основе теории Био и теоремы Бетти приращение пористости было выражено через линейную комбинацию приращений объемных деформаций и порового давления. В модификации этого соотношения для случая малой сжимаемости породного скелета зависимость между пористостью и объемными деформациями устанавливается также через экспоненту [40, 41]. В то же время в работах [31, 42, 43] предложена дробно-линейная зависимость пористости от объемных деформаций. Проведенный анализ литературы показывает, что выделяют три основных параметра, влияющих на эволюцию пористости: изменение порового давления, объемных деформаций и эффективных напряжений. Данная работа посвящена анализу ключевых моделей пористости применительно к ПГД, поскольку сложность и физическая связанность протекающих при этом процессов может приводить к принципиально различным результатам. Для описания ПГД была разработана модель, которая включает в себя уравнения баланса массы для трехфазной жидкости, закон сохранения энергии с учетом фазового перехода, возникающего при испарении пара и конденсации воды, а также закон Дарси для описания фильтрационного течения [44]. 1. Связанная термогидромеханическая модель фильтрации многофазной жидкости в пористой среде Описание всего многообразия явлений, возникающих при парогравитационном дренировании, является достаточно сложной задачей. Тепло- и массоперенос в деформируемой пористой среде характеризуется взаимодействием температурных, механических и гидродинамических процессов и трех различных фаз (твердой, жидкой и газообразной). В связи с этим в работе принимаются следующие гипотезы: · флюид, заполняющий поровое пространство, состоит из трех различных несмешиваемых компонентов (вода, пар, нефть); · фазовый переход обусловлен только температурными изменениями; · вследствие высокой пористости пород-коллекторов влияние капиллярного давления не учитывается; · течение каждой из составляющих флюида подчиняется закону Дарси; · продуктивный пласт изотропен, а его деформации малы. Уравнения, используемые для описания связанных термогидромеханических процессов, условно можно разделить на две группы: законы сохранения массы, энергии, момента импульса и определяющие соотношения (уравнения состояния), к которым относятся законы Дарси, Фурье, Гука и т.д. Таким образом, математическое описание процесса парогравитационного дренирования с учетом принятых выше гипотез включает в себя: · уравнение баланса массы воды (1) · уравнение баланса массы пара ; (2) · уравнение баланса массы нефти ; (3) · замыкающее соотношение ; (4) · уравнение баланса энергии с учетом конвективного теплопереноса и фазового перехода, вызванного испарением и конденсацией воды (5) · закон Дарси для описания фильтрации каждого из компонентов: ; (6) ; (7) ; (8) · уравнение равновесия ; (9) · определяющее соотношение для описания механического поведения пористой среды с учетом температурного расширения частиц грунта и влияния порового давления ; (10) · геометрическое соотношение для описания линейного тензора полной деформации . (11) В уравнениях (1)-(11): - пористость; , , - водонасыщенность, паронасыщенность, нефтенасыщенность (отношение части порового объема, занимаемой фазой, к общему объему пор); - время; - оператор Гамильтона; , , - скорости фильтрации воды, пара и нефти; , - массовые источники, вызванные фазовым переходом; - поровое давление; - ускорение свободного падения; , , - относительные проницаемости для воды, пара и нефти; , , - динамические вязкости воды, пара и нефти; - абсолютная температура; - тепловой источник, вызванный фазовым переходом; , , , - плотности воды, пара, нефти и частиц грунта; - эффективная плотность, ; , , , - удельные теплоемкости воды, пара, нефти и частиц грунта; , , , - коэффициенты теплопроводности воды, пара, нефти и частиц грунта; - эффективный коэффициент теплопроводности, ; - абсолютная проницаемость продуктивного пласта; - тензор напряжений; - тензор упругих постоянных, который в случае применения изотропной линейной теории упругости сводится к двум упругим постоянным (модулю всестороннего сжатия и модулю сдвига); - тензор малых деформаций; - коэффициент теплового расширения грунта; - начальное значение температуры; - коэффициент Био; - единичный тензор; - вектор перемещений. Для определения относительных проницаемостей использовались следующие соотношения [45]: (12) (13) (14) где , , , , , - эмпирические параметры; , , - остаточные значения водо-, паро- и нефтенасыщенности. Источники, появляющиеся в уравнениях (1), (2), (5) вследствие фазовых переходов, вызывающих испарение и конденсацию воды, имеют вид [46] , (15) (16) где - скрытая теплота парообразования; - эмпирический параметр, характеризующий скорость фазового перехода; - температура фазового перехода. Для описания изменения пористости горных пород разрабатываются различные модели, которые, как правило, сводятся к установлению взаимосвязи между пористостью и следующими параметрами: объемные деформации, поровое давление или эффективное напряжение. В данной работе анализируются некоторые ключевые модели, приведенные ниже. Одним из общепринятых соотношений, связывающим пористость с поровым давлением флюида, является [33] , (17) где - начальная пористость; - сжимаемость горной породы; - начальное давление. Для учета изменения пористости, вызванного влиянием объемной деформации и порового давления, Y. Wang и B. Lu [39] использовали следующее уравнение: . (18) Появление объемной деформации вызывает изменение пористости, которое определяется согласно соотношению [40] , (19) где - объемные деформации. Для одновременного учета порового давления и механических напряжений применяют уравнение [35] , (20) где - сжимаемость пор; ; . На основе лабораторных экспериментов в работе [37] было предложено соотношение для связи пористости с эффективными напряжениями , (21) где - остаточная пористость; - эмпирический параметр. В предположении постоянства объема частиц грунта может быть также получена следующая формула [36]: . (22) Система уравнений (1)-(11) дополнялась начальными и граничными условиями: , (23) , (24) , (25) , (26) , (27) , (28) , (29) , (30) , (31) , (32) , (33) , (34) , (35) , (36) , (37) , (38) где - начальное значение давления; - начальное значение температуры; - начальная водонасыщенность; - граница нагнетающей скважины; - давление нагнетаемого пара; - граница добывающей скважины; - давление в добывающей скважине; - вектор нормали; - вектор теплового потока; , - правая и левая границы расчетной области; , - верхняя и нижняя границы расчетной области. Таким образом, замкнутая система уравнений для описания парогравитационного дренирования включает уравнения (1)-(11), выражения (12)-(16), одно из соотношений (17)-(22), а также начальные и граничные условия (23)-(38). 2. Численная реализация модели Для решения системы уравнений был разработан численный алгоритм [44], основные этапы которого приведены ниже: · преобразование уравнений (1)-(4), (6)-(8) с использованием общей скорости потока : , (39) , (40) , (41) где ; ; ; · запись уравнений (39)-(41) в слабой форме при помощи умножения на тестовую функцию и применения формулы интегрирования по частям; · добавление к каждому из уравнений слагаемого с искусственной диффузией для стабилизации осцилляций численного решения: (42) (43) где , , , - тестовые функции; , , - параметры искусственной диффузии; · преобразование уравнения (5) к виду (44) · решение уравнений методом конечных элементов в пакете Comsol Multiphysics® с использованием следующих модулей: Weak Form PDE Interface для (42), (43), Heat Transfer для (44), Structural Mechanics (9)-(11). Для дискретизации по времени использовалась неявная схема Эйлера первого порядка. 3. Результаты численного моделирования эволюции пористости в процессе парогравитационного дренирования Трехмерное моделирование паротеплового воздействия на нефтяной коллектор проводилось для прямоугольной области, имеющей следующие размеры: толщина 5 м, высота 24 м, ширина 30 м. Добывающая скважина располагалась на расстоянии 10 м от подошвы пласта, нагнетающая - на 5 м выше. Скважины имели диаметр 0,178 м. Начальные и граничные условия, используемые при моделировании, даны в табл. 1. Физико-механические свойства продуктивного пласта, нефти, воды и пара представлены в табл. 2, 3. Зависимость динамической вязкости нефти от температуры изображена на рис. 1 [47]. Расчетная область разбивалась на тетраэдральные конечные элементы, общее число которых составляло 11 142. Размер элементов варьировался от минимального (вблизи скважины), равного 0,09 м, до максимального (на периферии), равного 1 м. Таблица 1 Начальные и граничные условия Table 1 Initial and boundary conditions Свойство Символ Значение Размерность Начальное значение нефтенасыщенности So0 0,88 Начальное значение водонасыщенности Sw0 0,12 Начальная температура продуктивного пласта T0 281,15 К Температура закачиваемого пара Tb 496 К Начальное значение давления P0 2·106 Па Давление закачиваемого пара Pb 3,575·106 Па Давление в добывающей скважине Pw 1,7·106 Па Таблица 2 Механические свойства Table 2 Mechanical properties Свойство Символ Значение Размерность Модуль Юнга E 3·109 Па Коэффициент Пуассона ν 0,25 Коэффициент теплового расширения грунта αт 5·10-6 1/К Коэффициент Био αB 0,44 Таблица 3 Теплофизические свойства Table 3 Thermophysical properties Свойство Символ Значение Размерность Абсолютная проницаемость продуктивного пласта K 2,68·10-12 м2 Вязкость воды µw 1004·10-6 Па∙с Вязкость пара µs 16·10-6 Па∙с Остаточная нефтенасыщенность Sro 0,68 Остаточная водонасыщенность Srw 0,0 Остаточная паронасыщенность Srs 0,0 Плотность нефти ρo 933 кг/м3 Плотность воды ρw 1012 кг/м3 Плотность скелета ρr 2100 кг/м3 Плотность пара ρs 0,3 кг/ м3 Теплопроводность нефти λo 0,14 Вт/(м∙К) Теплопроводность воды λw 0,58 Вт/(м∙К) Теплопроводность пара λs 23·10-2 Вт/(м∙К) Теплопроводность скелета λr 2,325 Вт/(м∙К) Теплоемкость нефти co 2090 Дж/(кг∙К) Теплоемкость пара cs 2000 Дж/(кг∙К) Теплоемкость воды cw 4200 Дж/(кг∙К) Теплоемкость скелета cr 1050 Дж/(кг∙К) Пористость n 0,27 Эмпирический параметр a1 75·10-2 Эмпирический параметр a2 40·10-3 Эмпирический параметр a3 12·10-3 Эмпирический параметр m1 3,45 Эмпирический параметр m2 3 Эмпирический параметр m3 2,1 Рис. 1. Зависимость динамической вязкости нефти от температуры [47] Fig. 1. Dynamic viscosity of oil versus temperature [47] На рис. 2-17 представлены результаты расчетов основных параметров, характеризующих процесс ПГД, согласно разработанной модели (1)-(16) и соотношений (17)-(22). На рис. 3, 4, 8, 9, 13, 14 показаны характерные распределения паро- и нефтенасыщенностей. Максимальная скорость распространения пара сосредоточена внутри паровой камеры. Значение нефтенасыщенности в паровой камере совпадает с остаточным. Нефть за пределами паровой камеры является неподвижной вследствие высокой вязкости. Результаты расчетов пористости согласно формулам (17)-(22) показали, что существует три качественно различных картины ее распределения по расчетной области. На рис. 2-6 показаны результаты, полученные согласно модели (18). В этом случае распределение пористости (рис. 5) определяется слагаемым, зависящим от порового давления (рис. 6), несмотря на присутствие в формуле объемной деформации. Аналогичный результат можно получить с использованием модели (17). Второй характерный вид представлен на рис. 7-11. В этом случае пористость (рис. 10) рассчитывалась согласно (19), и ее вид соответствует распределению объемной деформации (рис. 11). К третьему типу, показанному на рис. 12-16, относятся соотношения (20)-(22). На рис. 15 представлено распределение пористости, полученное согласно модели (21). Особенностью таких моделей является зависимость от эффективных механических напряжений, вид которых представлен на рис. 16. На рис. 15, 16 можно наблюдать образование «кольца», соответствующего фронту фазового перехода. В его пределах возникает повышенное уплотнение грунта. На начальном этапе ПГД в области, соответствующей паровой камере, происходит увеличение пористости вследствие влияния порового давления (рис. 15, а). В дальнейшем уплотнение грунта приводит к уменьшению порового давления и росту эффективных сжимающих напряжений (рис. 16, б, в), в результате вся расчетная область находится в состоянии сжатия (рис. 15, б, в). Модели пористости по-разному влияют на основные параметры, характеризующие ПГД, что может приводить к разным оценкам дебита нефти. При использовании формулы (19) паровая камера имеет минимальный размер. Соответствующие ей величины нефтенасыщенности и температуры приведены на рис. 7 и 9. Использование модели (21) приводит к повышению скорости ПГД (рис. 12-14). В этом случае паровая камера к 100-му дню достигает кровли пласта, что соответствует переходу от начального этапа роста паровой камеры ко второму этапу, связанному с расширением паровой камеры в горизонтальном направлении. Наиболее быстрый процесс ПГД наблюдается при использовании формулы (18) (рис. 2-4). Результаты показывают, что уже на 60-й день паровая камера распространяется вдоль кровли пласта более чем наполовину. а б в Рис. 2. Распределение температуры по расчетной области, полученное при использовании модели пористости (18): а - 20-й день; б - 60-й день; в - 100-й день Fig. 2. Evolution of temperature field obtained with the use of porosity model (18) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 3. Распределение паронасыщенности совместно с наложенным полем скорости по расчетной области, полученное при использовании модели пористости (18): а - 20-й день; б - 60-й день; в - 100-й день Fig. 3. Evolution of steam saturation together with the velocity field obtained with the use of porosity model (18) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 4. Распределение нефтенасыщенности по расчетной области, полученное при использовании модели пористости (18): а - 20-й день; б - 60-й день; в - 100-й день Fig. 4. Evolution of oil saturation obtained with the use of porosity model (18) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 5. Распределение изменения пористости по расчетной области, полученное при использовании модели пористости (18): а - 20-й день; б - 60-й день; в - 100-й день Fig. 5. Evolution of porosity change obtained with the use of porosity model (18) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 6. Распределение порового давления по расчетной области, полученное при использовании модели пористости (18): а - 20-й день; б - 60-й день; в - 100-й день Fig. 6. Evolution of pore pressure obtained with the use of porosity model (18) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 7. Распределение температуры по расчетной области, полученное при использовании модели пористости (19): а - 20-й день; б - 60-й день; в - 100-й день Fig. 7. Evolution of temperature field obtained with the use of porosity model (19) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 8. Распределение паронасыщенности совместно с наложенным полем скорости по расчетной области, полученное при использовании модели пористости (19): а - 20-й день; б - 60-й день; в - 100-й день Fig. 8. Evolution of steam saturation together with the velocity field obtained with the use of porosity model (19) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 9. Распределение нефтенасыщенности по расчетной области, полученное при использовании модели пористости (19): а - 20-й день; б - 60-й день; в - 100-й день Fig. 9. Evolution of oil saturation obtained with the use of porosity model (19) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 10. Распределение изменения пористости по расчетной области, полученное при использовании модели пористости (19): а - 20-й день; б - 60-й день; в - 100-й день Fig. 10. Evolution of porosity change obtained with the use of porosity model (19) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 11. Распределение величины объемной деформации по расчетной области, полученное при использовании модели пористости (19): а - 20-й день; б - 60-й день; в - 100-й день Fig. 11. Evolution of volumetric strain obtained with the use of porosity model (19) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 12. Распределение температуры по расчетной области, полученное при использовании модели пористости (21): а - 20-й день; б - 60-й день; в - 100-й день Fig. 12. Evolution of temperature field obtained with the use of porosity model (21) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 13. Распределение по расчетной области паронасыщенности совместно с наложенным полем скорости, полученное при использовании модели пористости (21): а - 20-й день; б - 60-й день; в - 100-й день Fig. 13. Evolution of steam saturation together with the velocity field obtained with the use of porosity model (21) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 14. Распределение нефтенасыщенности по расчетной области, полученное при использовании модели пористости (21): а - 20-й день; б - 60-й день; в - 100-й день Fig. 14. Evolution of oil saturation obtained with the use of porosity model (21) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 15. Распределение изменения пористости по расчетной области, полученное при использовании модели пористости (21): а - 20-й день; б - 60-й день; в - 100-й день Fig. 15. Evolution of porosity change obtained with the use of porosity model (21) for (a) - 20th day; (b) - 60th day; (c) - 100th day а б в Рис. 16. Распределение эффективного напряжения по расчетной области, полученное при использовании модели пористости (21): а - 20-й день; б - 60-й день; в - 100-й день Fig. 16. Evolution of effective stress obtained with the use of porosity model (21) for (a) - 20th day; (b) - 60th day; (c) - 100th day На рис. 17 приведено распределение параметра (разность между значением пористости на 100-й день и начальным значением) вдоль горизонтальной прямой, расположенной на расстоянии одного метра от кровли пласта (рис. 17, а), и центральной вертикальной прямой, соединяющей кровлю пласта с нагнетательной скважиной (рис. 17, б). Область вблизи кровли пласта представляет наибольший интерес, поскольку считается, что начальный период формирования паровой камеры заканчивается по достижении ее границы. Полученные данные показывают качественно и количественно различные результаты. Модели (20), (22) предсказывают малое изменение пористости. Наибольшее отличие по сравнению с начальным значением дает модель (21), которая показывает, что в окрестности паровой камеры существует область, в которой изменение пористости принимает постоянное значение, максимальное по абсолютной величине. Возрастающие ветви на рис. 17, а соответствуют тем частям расчетной области, куда паровая камера еще не успела распространиться. Возрастающий участок на рис. 17, б вызван наличием вокруг паровой камеры области повышенного порового давления радиусом 3 м, которое стремится расширить поры. Похожий, но менее ярко выраженный характер имеют результаты, полученные при помощи формулы (22). Соотношения (17), (18) дают противоположную формулам (21), (22) картину. Их применение позволяет описывать увеличение пористости вследствие влияния порового давления. При этом на рис. 17, а наибольшее значение пористости соответствует отрезку, пересекающему паровую камеру, а на рис. 17, б максимум совпадает с границей скважины, где задано давление нагнетаемого пара. Таким образом, соотношения (17)-(22) дают принципиально разные, как качественно, так и количественно, результаты. х, м а у, м б Рис. 17. Изменение пористости на 100-й день после начала ПГД, рассчитанное согласно различным моделям: 1 - формула (17); 2 - формула (18); 3 - формула (19); 4 - формула (20); 5 - формула (21); 6 - формула (22)); вдоль- горизонтальной линии, расположенной на расстоянии 1 м от кровли пласта (а); вдоль центральной вертикальной линии, соединяющей кровлю пласта с нагнетательной скважиной (б) Fig. 17. Change of porosity for 100th day obtained by various models: 1 - formula (17); 2 - formula (18); 3 - formula (19); 4 - formula (20); 5 - formula (21); 6 - formula (22)) along the horizontal line located 1 m from the upper boundary of the reservoir (а); along central vertical line which connects the upper boundary of the reservoir with the injection well (b) Проведенное исследование моделей пористости позволяет заключить, что для корректного описания дебита нефти в произвольном пласте-коллекторе необходимы дополнительные лабораторные испытания, позволяющие определить преобладающий механизм изменения пористости (температурное расширение, поровое давление или сжатие пласта под действием нагнетания пара). Выбор формулы для расчета пористости будет зависеть от установленного механизма. Заключение В работе предложена связанная термогидромеханическая модель фильтрации многофазной жидкости в пористой среде для описания парогравитационного дренирования. Модель описывает процессы конвективного теплообмена, фазового перехода, теплового расширения грунта и влияние порового давления на напряженно-деформированное состояние продуктивного пласта. Для ее реализации был разработан алгоритм, основанный на решении системы уравнений относительно неизвестных величин давления, паро-, водо- и нефтенасыщенностей. На основе разработанной модели проведено исследование особенностей различных соотношений для описания эволюции пористости. В результате было получено три характерных, качественно различных картины ее распределения по расчетной области. Модели первого типа основаны на установлении закономерностей между пористостью и поровым давлением. В результате их применения образуются две области, одна из которых расположена вблизи нагнетательной скважины (в ней пористость увеличивается), а другая - вблизи добывающей (пористость уменьшается). При связи пористости с объемной деформацией можно наблюдать увеличение пористости вследствие теплового расширения грунта. Для моделей третьего типа характерна зависимость от эффективных механических напряжений. В этом случае при паротепловом воздействии на продуктивный пласт образуется «кольцо» повышенного уплотнения грунта, которое соответствует фронту фазового перехода. Количественные расхождения между моделями могут достигать нескольких порядков. Результаты расчетов также показали, что применение различных моделей пористости приводит к изменению размеров паровой камеры, что влияет на оценку дебита нефти. Таким образом, проведенное исследование моделей пористости позволяет заключить, что для корректного описания дебита нефти в произвольном пласте-коллекторе необходимы дополнительные лабораторные испытания, позволяющие определить преобладающий механизм изменения пористости (температурное расширение, поровое давление или сжатие пласта под действием нагнетаемого пара).

About the authors

A A Kostina

Institute of Continuous Media Mechanics UB RAS

M S Zhelnin

Institute of Continuous Media Mechanics UB RAS

O A Plekhov

Institute of Continuous Media Mechanics UB RAS

References

Statistics

Views

Abstract - 460

PDF (Russian) - 311