CC BY
12
УДК 550.832.6
05.00.00 Технические науки
МОДЕЛИРОВАНИЕ ФИЛЬТРАЦОННО-ВОЛНОВОГО НАГРЕВА НЕФТЯНОГО ПЛАСТА
Коныс Ерканат Магзомович Студент 3-го курса, [email protected]
Кульчубаев Нурислам Булатович Студент 3-го курса
Шантиев Едиге Исмагулович Студент 3-го курса,
Хисматуллин Азат Салаватович к.ф-м.н., доцент
Филиал ««Уфимский государственный нефтяной технический университет» в г. Салавате, 453250, Россия, Республика Башкортостан, г. Салават, ул. Губкина, д. 22Б.
В данной статье разработана математическая модель термодинамических процессов при колебательном движении жидкостей в пористой среде с учетом теплового взаимодействия нефтяного пласта с окружающими породами и фазовых переходов вследствие растворения парафинов; разработка программного средства на языке высокого уровня и проведение вычислительного эксперимента, обосновывающие возможность увеличения нефтеотдачи пластов при воздействии на них периодическими полями давления. Предложено использование баротермического эффекта для увеличения температуры удаленных от скважины зон пласта. Математическая модель базируется на уравнении теплового баланса, учитывающего вклад фазовых переходов за счет парафинизации в случае, когда движение смеси происходит в однородной пористой изотропной среде с постоянной пористостью и постоянной плотностью. В работе развита теория баротермического эффекта при колебательном движении жидкостей в пористой среде с учетом теплового взаимодействия нефтяного пласта с окружающими породами и фазовых переходов вследствие растворения парафинов. Сопоставление теоретических и экспериментальных данных подтвердило возможность использования первого приближения для расчетов реальных термодинамических процессов в фильтрационно-волновых полях. Разработанная модель процесса позволяет рассчитывать температурные поля в практических условиях
Ключевые слова: ПАРАФИНИЗАЦИЯ, ТЕМПЕРАТУРНЫЙ ЭФФЕКТ, ПАРАФИНИСТАЯ НЕФТЬ, АСИМПТОМИЧЕСКИЙ МЕТОД
UDC 550.832.6 Engineering
FILTRATION-WAVE HEATING OF THE OIL POOL
Konys Yerkanat Magzomovich Student of the 3rd course, erkanat [email protected]
Kulchubayev Nurislam Bulatovich Student of the 3rd course
Shantiev Yedige Ismagulovich Student of the 3rd course
Hismatullin Azat Salavatovich Cand.Fiz-Mat.Sci., associate professor Branch of "Ufa state petroleum technological University" in the city of Salavat, 453250, Russia, Bashkortostan Republic, Salavat, St. Gubkin, D. 22B.
We have developed a mathematical model of thermo-dynamic processes at oscillatory movement of liquids in porous environment with consideration of heat interaction of an oil layer with surrounding breeds and phase transitions caused by dissolution of paraffin's dissolution. We have received analytical dependences of temperatures on time, distances to well and depth, that proves an opportunity to increase of oilrecavery when influenced by periodic fields of pressure. We have suggested a bartering effect to increase the temperature in zones far from the well formation. The mathematical model is based on the equation of heat balance, taking into account the contribution of phase transitions due to paraffinization in the case when the motion of the mixture takes place in a homogeneous isotropic porous medium with constant porosity and constant density. In the work we have developed a theory of bartering effect upon oscillatory movement of fluids in porous media considering thermal interaction between the oil reservoir and surrounding rocks and phase transitions due to the dissolution of paraffins. Comparison of theoretical and experimental data has confirmed the possibility of using the first approximation for the calculation of real thermodynamic processes in the filtration-wave fields. The developed process model allows the calculation of the temperature field in the practical conditions
Keywords: WAXING, THERMAL EFFECT, PARAFFIN OIL, ASYMPTOTIC METHOD
йо!: 10.21515/1990-4665-128-095
В работах [1-9] предложено использование баротермического эффекта для увеличения температуры удаленных от скважины зон пласта с использованием программного средства на языке высокого уровня. Дело в том, что поля давления в нефтеносных пластах при наличии фильтрацион-но-волновых процессов нестационарны, поэтому фильтрация нефти в пластах приводит к появлению баротермического эффекта - изменению температуры при течении флюида в пористой среде в нестационарном поле давления. Величина баротермического эффекта зависит от свойств пористой среды, фильтрующейся жидкости, времени, геометрии течения и т.д. Одной из наиболее важных особенностей баротермического эффекта является возрастание температуры со временем в фильтрационно-волновых полях.
Математическая модель базируется на уравнении теплового баланса, учитывающего вклад фазовых переходов за счет парафинизации в случае, когда движение смеси происходит в однородной пористой изотропной среде с постоянной пористостью т=сопБ1 и постоянной плотностью Р=СОПБ1:
1 + —/1 (Т)
с„
дт , ,
ъ + ^ о-,')
1 + /И'
+ -
Срр
дТ , , дР сж дР
— + и^ (Г, t) • еет ---ЛеТт
дг дг сп дt
1 д
г дг
где
Г ЗтЛ
г
дТ
дг
2Т (1)
+ а2
у иг у
д_Т
д2:
и/(г,0 -^РМР = «л (2)
сп т I дг сп т г
/1 (т )=А др=с
т - т
1 1 н
т - т
V т 1н у
ехр
- 2-
т - т
1 1 н
т - т
±т -'и у
(3)
2
тн, тт - температуры начала и максимума кристаллизации парафина соответственно,
С - постоянная для данной нефти величина,
Ьх - удельная энтальпия фазового превращения парафина. Скорость фильтрации подчиняется закону Дарси:
к
Я = —VР. (4)
т
В основу положено также уравнение теплопроводности флюида
дт 1 э
— = аг--
Эt г Эг
эт^ э2т
^ J + . (5)
V
Рассмотрим волновые поля давления и скорости, возникающие в пористой среде в процессе фильтрации жидкости в определенных условиях. Пласт считается однородным и изотропным.
Предположим, что твердая фаза парафина и скелет несжимаемы и поле давления в пористой среде определяются, в основном, сжимаемостью несущей жидкой фазы; температура т и связанные с ней параметры , , ^ существенно не меняются за период фильтрационно-волновых колебаний, то есть являются адиабатическими параметрами.
Тогда в адиабатическом приближении уравнение для фильтрации сжимаемой жидкости запишем в виде:
ЭР
тр-кДР = 0, (6)
дt
р0 к )
где к = 0 , н! /(?п) - коэффициент пьезопроводности.
т{\- ^^ )
Уравнение (6) называют также уравнением упругого режима. Коэффициент пьезопроводности к характеризует скорость перераспределения пластового давления при неустановившейся фильтрации упругой жидкости в упругой пористой среде. Наиболее часто встречающиеся в нефте-
промысловой практике значения коэффициента пьезопроводности заключены в пределах от 0.1 до 5 м/с.
Рассмотрим задачу об изменении температуры жидкости при колебательном движении в пористой среде с учетом фазовых переходов, обусловленных растворением парафина. Пусть горизонтальный пласт представлен в виде трех полубесконечных областей, с плоскими границами раздела, перпендикулярными оси г. Первая I и вторая II области непроницаемы; средняя область толщины 2к, расположенная горизонтально, является пористой и насыщена парафинистой нефтью. Будем рассматривать случай радиального колебательного движения парафинистой нефти в средней области -к<г<к. При описании температурной задачи примем, что температура нефти, парафина и скелета пористой среды в каждой точке совпадают, отлична от нуля только радиальная координата скорости конвективного переноса тепла, т.е. ифО, иф=0, иг=0.
Решение задачи в предельном случае нулевого приближения с учетом фазовых переходов соответствует отысканию осредненной по толщине пласта температуре. Задача для нулевого приближения с учетом парафини-зации примет вид:
ОН - ал А [г ЯП = 0 г > 1, г > 0, , > 0, (7
д( а21 Ог у Ог ) Ог2
О?
(г,
'2 (Г + /(г,1
ОГ
Г > 0, г > 0, | г |< 1;
С ( ОГ 1г2 ОГ2 \
2 [ Ог г=1 1г1 Ог г=_1)
(8)
ОГ2 аг 2 1 О
О? аг1 г Ог
ОГ2
Ог
а2 2 О Г
а21 Ог 2
= 0, г <_1, г > 0, ? > 0
(9)
Граничные и начальные условия
Tz=1 = T21 z=-l = T , (10)
T
= T
11/=0 21/=0
=T n= 0.
I/=0
(11)
limT
1,2
= 0.
r + z
(12)
За одно колебательное движение температура фильтрующейся жидкости изменяется незначительно, поэтому зависимостью скорости конвективного переноса тепла от температуры пренебрегаем. Зависимость для определения температуры фильтрующейся жидкости с учетом парафини-зации примет вид
,эг
,дГ
Ф1 (Г)— + Uef (r, t)Ф2 (Г)--D1 (r)Sin2 wt -D2 (r)cos wt = 0:
dt
dr
(13)
где Dx(r) = eef
cm к Pf
0
сп m ln2 f rA
D2 (r ) = hef P0
ln(r / R)
40/ R)'
R
Решим уравнение (13) используя программный пакет Delphi. Для этого воспользуемся методом сеток.
Основные результаты представлены в виде следующих трехмерных графиков:
Рис. 1 График зависимости температуры от радиального расстояния г,м и времени ¿,с: ю:=0.32Гц; г0:=0.1м; £=9-10-11м2/Дж; |1=11.7-10-3Па- с;
6
Я:=100м; Р0=Ь106 Па; 8ef:=0.04•10"5К/Па; nef:=0.0137^10"5 К/Па; сж= 1600000Дж/(м3^К); сп = 2120000 Дж/(м3-К).
2
Как видно из графика температура возрастает с течением времени по закону синуса, а по радиусу убывает.
Разработанная математическая модель процесса подтверждает возможность практического использования фильтрационно-волнового нагрева пластов для увеличения нефтеотдачи и позволяет рассчитывать температурные поля в практических условиях.
Выводы. В работе развита теория баротермического эффекта при колебательном движении жидкостей в пористой среде с учетом теплового взаимодействия нефтяного пласта с окружающими породами и фазовых переходов вследствие растворения парафинов. Сопоставление теоретических и экспериментальных данных подтвердило возможность использования первого приближения для расчетов реальных термодинамических процессов в фильтрационно-волновых полях.
СПИСОК ЛИТЕРАТУРЫ
1. Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. М.: Недра, 1984. 211 с.
2. Ефимова Г.Ф. Математическое моделирование температурных процессов в фильтрационно-волновых полях с учетом фазовых переходов // Диссертация на соискание ученой степени кандидата физико-математических наук. Стерлитамак, 2004
3. Филиппов А.И., Хисматуллин А.С., Михайлов П.Н. Фильтрационно-волновой нагрев нефтяного пласта // Инженерная физика М.: Научтехиздат, 2006., №5, С. 13 - 22.
4. Филиппов, А.И. Установка для исследования коэффициента температуропроводности в жидкости / А.И. Филиппов, М.Р. Минлибаев, А.С. Хисматуллин // Новые промышленные технологии. 2010. № 2. С. 62 _ 63.
5. Хисматуллин А.С., Филиппов А.И., Минлибаев М.Р., Серебренников Н.П. Определение коэффициента трансцилляторного переноса при барботаже в жидкости // Вестник Воронежского государственного технического университета. - 2010. - № 2. -С. 52 - 53.
6. Хисматуллин А. С., Сагитов К.И., Абсатаров И.Х. Увеличение нефтеотдачи пласта при помощи фильтрационно-волнового нагрева //Успехи современной науки и образования. 2016. Т. 4. № 7. С. 103-107.
7. Хисматуллин А.С., Минлибаев М.Р., Сагитов К.И., Абсатаров И.Х. Исследование зависимости работы, необходимой для извлечения жидкости из пласта, от его пористости // Успехи современной науки. 2016. Т. 4. № 7. С. 92-97.
8. Хисматуллин А.С., Кузнецов А.С., Мусин А.Р., Вильданов Р.Р., Юсупов Р.Р., Когот А.Ю. Исследования технического состояния трансформаторных подстанции 10, 35 и 110 кВ // Успехи современной науки. 2017. Т. 4. № 1. С. 55-59.
9. Хисматуллин А.С. Совместная фильтрация двух несмешивающихся жидкостей в пластах // XII Всероссийская научная конференция студентов-физиков и молодых ученыхматериалы конференции и тезисы докладов. 2006. С. 690-691.
10. Чекалюк Э.Б. Термодинамика нефтяного пласта. М.:Недра. 1965. 238 с.
11. Яворский Б. М., Детлаф А. А. Справочник по физике. - М.: Наука, 1971. - 940 с.
References
1. Barenblatt G.I., Entov V.M., Ryzhik V.M. Dvizhenie zhidkostej i gazov v prirodnyh plastah. M.: Nedra, 1984. 211 s.
2. Efimova G.F. Matematicheskoe modelirovanie temperaturnyh processov v fil'tracionno-volnovyh poljah s uchetom fazovyh perehodov // Dissertacija na soiska-nie uchenoj stepeni kandidata fiziko-matematicheskih nauk. Sterlitamak, 2004
3. Filippov A.I., Hismatullin A.C., Mihajlov P.N. Fil'tracionno-volnovoj nagrev neftyanogo plasta // Inzhenernaya fizika M.: Nauchtekhizdat, 2006g., №5, C. 13 - 22.
4. Filippov, A.I. Ustanovka dlya issledovaniya koehfficienta temperaturoprovodnosti v zhidkosti / A.I. Filippov, M.R. Minlibaev, A.S. Hismatullin // Novye promyshlennye tekhnologii. 2010. № 2. s. 62.
5. Hismatullin A.S., Filippov A.I., Minlibaev M.R., Serebrennikov N.P. Opredelenie koehfficienta transcillyatornogo perenosa pri barbotazhe v zhidkosti // Vestnik Voronezhskogo gosudarstvennogo tekhnicheskogo universiteta. - 2010. - № 2. - S. 52-53.
6. Hismatullin A.S., Sagitov K.I., Absatarov I.H. Uvelichenie nefteotdachi plasta pri pomoshchi fil'tracionno-volnovogo nagreva //Uspekhi sovremennoj nauki i obrazovaniya. 2016. T. 4. № 7. S. 103-107.
7. Hismatullin A.S., Minlibaev M.R., Sagitov K.I., Absatarov I.H. Issledovanie zavisimosti raboty, neobhodimoj dlya izvlecheniya zhidkosti iz plasta, ot ego poristosti // Uspekhi sovremennoj nauki. 2016. T. 4. № 7. S. 92-97.
8. Hismatullin A.S., Kuznecov A.S., Musin A.R., Vil'danov R.R., YUsupov R.R., Kogot A.YU. Issledovaniya tekhnicheskogo sostoyaniya transformatornyh podstancii 10, 35 i 110 kV // Uspekhi sovremennoj nauki. 2017. T. 4. № 1. S. 55-59.
9. Hismatullin A.S. Sovmestnaja fil'tracija dvuh nesmeshivajushhihsja zhidkostej v plastah / XII Vserossijskaja nauchnaja konferencija studentov-fizikov i molodyh uchenyhma-terialy konferencii i tezisy dokladov. 2006. S. 690-691.
10. Chekalyuk E.B. Termodinamika neftyanogo plasta. M.:Nedra. 1965, 238 s.
11. Yаvorskij B. M., Detlaf A. A. Spravochnik po fizike. - M.: Nauka, 1971. - 940 s.
