ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ТЕПЛОВЫХ ПРОЦЕССОВ В СКВАЖИНЕ С ПРИЗАБОЙНЫМ НАГРЕВАТЕЛЕМ
- Авторы: Пинягин Д.С1, Костарев Н.А1, Труфанова Н.М1
- Учреждения:
- Пермский национальный исследовательский политехнический университет
- Выпуск: № 36 (2020)
- Страницы: 48-62
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/elinf/article/view/2453
- DOI: https://doi.org/10.15593/2224-9397/2020.4.03
- Цитировать
Аннотация
По причине снижения природных запасов нефти в нефтедобыче остро стоит вопрос освоения месторождений с высоковязкой нефтью. Одним из способов добычи высоковязкой нефти является прогрев призабойной зоны скважины при помощи нагревателя. Цели: разработка и исследование модели процесса тепломассопереноса в нефтяной скважине с призабойным нагревателем путем сравнения различных моделей турбулентного и ламинарного течения и выявление оптимальной. Методы: в работе рассмотрена двухмерная осесимметричная математическая модель процесса тепломассопереноса в нефтяной скважине с призабойным нагревателем, который используется для скважин с высоковязкой нефтью. Использование локального нагревателя в призабойной области позволяет снизить вязкость нефти, поступающей в насос, обеспечивая тем самым бесперебойную работу скважины. Рассматривалось продольное сечение нефтяной скважины глубиной 105 м, окружённой грунтом. Нагреватель расположен над забоем. Численная реализация разработанной математической модели осуществлялась методом конечных элементов в ANSYS. Для построения геометрии и разбиения на сетку конечных элементов использовался препроцессор ICEM CFD. Поставленная задача решалась численным методом конечных элементов в среде инженерных расчетов Fluent. Результаты: разработана математическая модель нефтяной скважины с турбулентным и ламинарным течением нефти с учетом призабойного нагревателя. В результате исследования получены поля скоростей и температур во всем объеме скважины. Построены распределения скоростей и температур по сечению насосно-компрессорных труб, произведено сравнение различных моделей турбулентного течения и k -epsilon модели. Исследовано влияние величины вязкости на температуру перед входом в электрический центробежный насос. Выбрана наиболее оптимальная модель течения при данных условиях. Практическая значимость: полученные результаты позволяют выбрать оптимальную модель течения для расчета данной задачи и принять необходимые меры по снижению вязкости нефти на заранее известном участке, что существенно снизит материальные затраты при разработке месторождений.
Ключевые слова
призабойный нагреватель, нефтяная скважина, турбулентное течение, высокая вязкость, тепловой процесс.
Полный текст
Введение. В связи со снижением природных запасов нефти все большее внимание начинают привлекать способы разработки месторождений с высоковязкой нефтью. Проблема усиливается тем фактом, что большая часть месторождений высоковязкой нефти требует особых методов ее добычи. Решению проблемы добычи высоковязкой нефти посвящено большое количество работ [1-3], в которых рассматриваются различные технологии перекачки нефти. Авторами исследуется как турбулентное [4-6], так и ламинарное [7-9] течение, а также проводится сравнительный анализ результатов расчета. Стоит отметить, что в качестве одного из эффективных способов добычи высоковязкой нефти рассматривается нагрев призабойной зоны скважины при помощи электрических нагревателей, что делает задачу проработки данной технологии достаточно актуальной. Постановка цели и задач. Целью данного исследования является разработка математической модели процессов тепломассопереноса [10-15] в нефтяной скважине с призабойным нагревателем для снижения вязкости нефти перед входом в электрический центробежный насос. В рамках исследования приняты следующие допущения: задача - осесимметричная, стационарная; течение нефти - турбулентное. Теплофизические свойства твердых материалов не зависят от температуры, бесконечный массив земли заменен ограниченной областью. На границе раздела разнородных сред выполняется условие идеального контакта. Схематическое изображение исследуемой области приведено на рис. 1. Рис. 1. Граничные условия Математическая модель движения и теплопереноса в нефтяной скважине основывается на законах сохранения энергии, сохранения количества движения и сохранения массы [6, 16-19]. Система уравнений имеет вид [20], представленный ниже. Уравнения переноса (уравнения Навье-Стокса, осреднённые по Рейнольдсу): (1) (2) Уравнение переноса для кинетической энергии турбулентности k: (3) Уравнение переноса для скорости диссипации кинетической энергии турбулентности e: (4) Турбулентная вязкость mt: (5) Уравнение неразрывности: (6) Уравнение энергии для жидких элементов: (7) Уравнение энергии для твердых элементов: (8) Уравнение энергии для источника: (9) Граничные условия задачи [20] в соответствии с рис. 1 приведены в таблице. Граничные условия задачи моделирования Положение Обозначение Величина Нагреватель Г1 Источник, Вт/м3 Центр скважины Г2 Ось симметрии Вход в скважину Г3 Дебит жидкости 50 м3/сут., температура 20 °С Выход из скважины Г4 Дебит жидкости 50 м3/сут Граница грунта Г5 Геотерма Нижняя часть Г6 Температура 20 °С Верхняя часть Г7 Конвективный теплообмен Численное моделирование процессов в скважине. Рассмотрим основные результаты проведенных исследований в соответствии с разработанными моделями процесса тепломассопереноса в нефтяной скважине с призабойным нагревателем. 1. Исследование сходимости. Исследование процесса сходимости решения (рис. 2, 3) показало, что установившееся решение может быть получено при использовании 700 000 элементов и 100 итераций. 500 000 2 500 000 4 500 000 6 500 000 Количество элементов, шт 21,1 21 20,9 20,8 20,7 20,6 Рис. 2. Зависимость средней температуры на выходе из насосно-компрессорной трубы от количества элементов а Рис. 3. Зависимость температуры (а) и скорости (б) на выходе из насосно-компрессорной трубы от числа итераций б Рис. 3. Окончание 2. Сравнение моделей течения. Результаты численного исследования для случая ламинарного и турбулентного течений представлены на рис. 4 и 5. Рис. 4. Распределение температур по сечению НКТ перед ЭЦН: турбулентное (1) и ламинарное (2) течение На рис. 4 показано, что температура распределяется равномерно по всему сечению насосно-компрессорной трубы для турбулентного течения. Это связано с тем, что элементы жидкости совершают неупорядоченные, неустановившиеся движения по сложным траекториям, что приводит к интенсивному перемешиванию между слоями движущейся жидкости. В случае ламинарного течения перепад температур по сечению составляет 6 °С. Поля продольной скорости для двух видов течения приведены на рис. 5. Максимальное отличие не превосходит 5 %. Рис. 5. Распределение скоростей по сечению НКТ перед ЭЦН: турбулентное (1) и ламинарное (2) течение 3. Сравнение различных турбулентных моделей. Рассмотрены три наиболее используемые в численном анализе математические модели турбулентного течения: Sapalart-Allmaras, k-epsilon, k-omega. Распределение скоростей и температур представлено на рис. 6 и 7. Как видно из рис. 6, все модели турбулентности показали близкие друг к другу результаты по значениям скоростей и температур. Отличие как по одному, так и по другому параметру не превосходит доли процента. Кривые практически совпадают. Таким образом, для описания течения в скваженном канале можно выбрать наиболее распространенную модель k-epsilon. 4. Сравнение различных k-epsilon моделей. Поскольку существует ряд разновидностей k-epsilon моделей, отличающихся формой записи турбулентной вязкости, проведено сравнение наиболее используемых математических моделей турбулентности: Realzable, Standard, RNG. Результаты сравнения моделей приведены на рис. 7. а б Рис. 6. Сравнение турбулентных моделей при расчете значений скорости (а) и температуры (б) на выходе НКТ: 1 - Sapalart-Allmaras; 2 - k-epsilon; 3 - k-omega а б Рис. 7. Сравнение k-epsilon моделей: 1 - Realzable, 2 - Standard, 3 - RNG Анализ данных на рис. 7 позволяет сделать вывод, что все k-epsilon модели показали близкие друг к другу результаты по значениям показателей скорости и температуры. Максимальное отклонение по скорости не превысило 2 %. Поскольку из трех рассмотренных моделей наибольшей скоростью сходимости обладает модель realzable, то для дальнейших исследований была выбрана k-epsilon realzable модель. 5. Исследование влияния величины вязкости на поле температур и скоростей. Исследование проводилось как при высоких, так и при низких значениях вязкости для случаев ламинарного и турбулентного течения. Рис. 8. Распределение температур по сечению НКТ при различной вязкости: сплошная линия - ламинарное течение, штриховая линия - турбулентное; начальная вязкость: синий - 10 Па×с, зеленый - 5 Па×с, оранжевый - 1 Па×с, красный - 0,5 Па×с, фиолетовый - 0,35 Па×с, желтый - 0,2 Па×с, розовый - 0,1 Па×с, черный - 0,001 Па×с Рассмотрены 8 различных ситуаций, относящихся к значениям вязкости нефти от 0,001 до 10 Па×с. Как видно из рис. 8 и 9, наибольшие перепады температур по сечению НКТ возникают при вязкости в 0,2 и 0,1 Па×с и составляют 11 и 9 ºС. Рис. 9. Распределение скоростей по сечению НКТ при различной вязкости: сплошная линия - ламинарное течение, штриховая линия - турбулентное; начальная вязкость: синий - 10 Па×с, зеленый - 5 Па×с, оранжевый - 1 Па×с, красный - 0,5 Па×с, фиолетовый - 0,35 Па×с, желтый - 0,2 Па×с, розовый - 0,1 Па×с, черный - 0,001 Па×с Наибольший перепад скоростей возникает при вязкости 0,1 Па×с и составляет 16,5 %. В остальных случаях перепад по значениям температуры не превышает 8 ºС и 4 % по значениям скоростей. Анализ представленных на рисунках результатов также свидетельствует о том, что особенности турбулентного течения начинают проявляться со значения вязкости в 0,2 Па×с. Выводы. Анализ турбулентного и ламинарного течения показал, что разница в результатах по значениям скоростей на выходе НКТ составляет не более 5 %, а разница по значениям температурам - около 6 °C. На основании сравнения различных моделей турбулентности в качестве оптимальной модели для данной задачи выбрана k-epsilon realizable модель. Рассмотрены 8 видов вязкости нефти от 0,001 до 10 Па×с. Наибольшие перепады температур по сечению возникли при вязкости в 0,2 и 0,1 Па×с и составили 11 и 9 ºС. Наибольший перепад скоростей возник при вязкости 0,1 Па×с и составил 16,5 %. В остальных случаях перепад по температурам не превышает 8 ºС и 4 % по значениям скоростей. Показано, что особенности турбулентного течения начинают проявляться при значении вязкости в 0,2 Па×с. Разработанная модель успешно справляется с задачей по определению наиболее рационального способа прогрева призабойной зоны скважины, что позволяет снизить вязкость нефти перед входом в электрический центробежный насос как при ламинарном, так и при турбулентном течении нефти. Использование данной модели позволяет принять необходимые меры по снижению вязкости нефти на заранее известном участке, что существенно снизит материальные затраты при разработке месторождений.Об авторах
Д. С Пинягин
Пермский национальный исследовательский политехнический университет
Н. А Костарев
Пермский национальный исследовательский политехнический университет
Н. М Труфанова
Пермский национальный исследовательский политехнический университет
Список литературы
- Тарасюк В.М. Высоковязкие нефти и природные битумы: проблемы и повышение эффективности разведки и разработки // Берегиня.777.сова. - 2014. - № 2. - С. 121-125.
- Башкирцева Н.Ю. Высоковязкие нефти и природные нефти // Вестник Казан. технолог. ун-та. - 2014. - С. 296-299.
- Шандрыгин А.Н. Разработка залежей тяжелой нефти и природного битума методом парогравитационного дренажа (SAGD) // Нефтяное хозяйство. - 2006. - № 7. - С. 92-96.
- Исследование процессов турбулентного движения однородной жидкости в стволе горизонтальной скважины при различных типах перфорационных отверстий / М.А. Фатхлисламов, И.В. Владимиров, О.П. Торопчин, С.А. Кротов // Автоматизация, телемеханизация и связь в нефтяной промышленности. - 2009. - № 9. - С. 36-40.
- Исследование процессов турбулентного, вязкого движения двухфазной жидкости в полого направленном стволе скважины при различных свойствах коллектора / М.А. Фатхлисламов, Т.Г. Казакова, О.П. Торопчип, А.С. Кротов // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2010. - № 1. - С. 79-84.
- Филиппов А.И., Ахметова О.В., Родионов А.С. Температурное поле турбулентного потока в скважине // Теплофизика высоких температур. - 2013. - Т. 51, № 2. - С. 277.
- Харламов С.Н., Терещенко Р.Е. Теплообмен при ламинарном течении нефти и нефтепродуктов с аномальной реологией в трубопроводах большой протяженности // Вестник Рос. акад. естеств. наук. Западно-Сибирское отделение. - 2013. - № 15. - С. 96-100.
- Ананьев Д.В., Халитова Г.Р. Интенсификация теплопереноса при ламинарном течении высокопарафинистой нефти в круглой трубе // Труды Академэнерго. - 2014. - № 2. - С. 7-16.
- Стрункин С.И., Григорьев А.В., Хузин Л.И. Применение установок прогрева призабойной зоны пласта // Инженерная практика. - 2015. - Вып. 12.
- Ибрагимов А.М. Тепломассоперенос при граничных условиях второго и третьего рода // Промышленное и гражданское строительство. - 2006. - № 9. - С. 58-59.
- Труфанова Н.М., Костарев Н.А. Исследование процесса тепломассопереноса в нефтяной скважине // Научно-технический вестник Поволжья. - 2015. - № 4. - С. 126-129.
- Конюхов В.М., Конюхов И.В., Чекалин А.Н. Математическое моделирование и параллельные вычисления процессов тепломассопереноса при эксплуатации единого комплекса «Нефтяной пласт - система скважин - электроцентробежные насосы» // Вестник Казан. технолог. ун-та. - 2014. - Т. 17. - № 10. - С. 195-202.
- Файзуллин А.Э., Труфанова Н.М. Анализ тепловых процессов в нефтяной скважине // Научно-технический вестник Поволжья. - 2018. - № 7. - С. 150-153.
- Ковригин Л.А., Кухарчук И.Б. Конвективный теплообмен в межтрубном пространстве нефтяной скважины с греющим кабелем // Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. - 2012. - № 6. - С. 43-47.
- Фертиков М.Г., Дятлов И.Я., Труфанова Н.М. Исследование тепловых режимов греющего кабеля различных конструкций // Автоматизированные системы управления и информационные технологии: материалы всерос. науч.-техн. конф.: в 2 т. - Пермь: Изд-во Перм. нац. исследов. политехн. ун-та, 2018. - Т. 2. - С. 115-120.
- Костарев Н.А., Труфанова Н.М. Исследование процессов тепломассопереноса в нефтяной скважине с учетом процесса парафинообразования // Автоматизированные системы управления и информационные технологии: материалы всерос. науч.-техн. конф.: в 2 т. - Пермь: Изд-во Перм. нац. исследов. политехн. ун-та, 2017. - Т. 2. - С. 136-142.
- Ming-Jing Du, Yu-Lan Wang, Chao-LuTemuer. Reproducing kernel method for numerical simulation of downhole temperature distribution // Applied Mathematics and Computation. - 2017. - 297. - P. 19-30.
- Tarom N., Hossain M.M. Using ANSYS to Realize a Semi-Analytical Method for Predicting Temperature Profile in Injection/Production Well // World Academy of Science, Engineering and Technology. - 2012. - 72.
- Труфанова Н.М., Костарев Н.А. Моделирование процесса тепломассопереноса при промывке нефтяной скважины через полые штанги // XXI Зимняя школа по механике сплошных сред: тез. докл., (г. Пермь, 18-22 февр. 2019 г.) / Урал. отд-ние Рос. акад. наук, Ин-т механики сплошных сред УрО РАН - филиал Федер. гос. бюджет. учр. науки «Перм. федер. исслед. центра» УрО РАН. - Пермь: Новопринт, 2019. - С. 158.
- Пинягин Д.С., Костарев Н.А., Труфанова Н.М. Анализ процессов тепломассопереноса в нефтяной скважине при использовании призабойных нагревателей // Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. - 2019. - № 30. - C. 211-226.
Статистика
Просмотры
Аннотация - 51
PDF (Russian) - 24
Ссылки
- Ссылки не определены.