Математическая модель многофазной затопленной струи Текст научной статьи по специальности «Физика»

где l - удельная теплота образования гидрата.

В работе рассмотрен теоретические основы процесса всплытия пузырьков метана внутри купола и их превращение в гидратные частицы. При рассмотрении модели было принято, что интенсивность образования гидратного пузырька лимитируется интенсивностью отвода тепла.

Работа выполнена при поддержке гранта СФ БГУ, договор № В14-6.

Литература

1. Макогон Ю. Ф. Гидраты природных газов. - М.: Недра, 1974. - 285с

2. Математическое моделирование процесса отбора газа из пористой среды, частично насыщенной газовым гидратом // Кильдибаева С.Р.// Фундаментальные исследования. № 11 (часть 6) 2013, стр. 1163-1166.

3. Моделирование купола-сепаратора при разливе нефти в шельфе // Кильдибаева С.Р. // Фундаментальные исследования. 2013. № 10 (часть 5). с. 1045-1050.

4. Нигматулин Р. И. Динамика многофазных сред. Т. 1. - М.: Наука, 1987.- 464 с.

5. Gumerov N.A., Chahine G.L. Dynamics of bubbles in conditions of gas hydrate formation // Fluid Dynamics. 1992. № 5. pp. 664669.

6. Sautera E.J. et al. Methane discharge from a deep-sea submarine mud volcano into the upper water column by gas hydrate-coated methane bubbles // Earth and Planetary Science Letters. 2006. № 243(3-4). pp. 354-365.

References

1. Makogon Y.F. Natural gas hydrates. - Moscow: Nedra, 1974. - 285 р.

2. Mathematical modeling of gas extraction from a porous medium partially saturated gas hydrates // Kildibaeva S. R. // Fundamental research. Number 11 (Part 6) 2013, pp. 1163-1166.

3. Modeling dome separators for oil spills on the shelf // Kildibaeva S. R. // Fundamental research. Number 10 (Part 5) 2013, pp. 1045-1050.

4. Nigmatulin R.I. Dynamics of Multiphase Media. T. 1. - Moscow: Nauka, 1987. - 464 р.

Кильдибаева С.Р.

Аспирант, Стерлитамакский филиал Башкирского государственного университета Работа выполнена при поддержке гранта СФ БГУ, договор № В14-6.

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ МНОГОФАЗНОЙ ЗАТОПЛЕННОЙ СТРУИ

Аннотация

В статье рассматриваются математические модели многофазной затопленной струи. Струя состоит из капель нефти и пузырьков газа. В работе проведен расчеты теплофизических параметров затопленной струи, получены графики зависимости температуры и скорости от вертикальной координаты. Особенное внимание в исследовании струи уделено процессу охлаждения струи из-за «захвата» окружающей воды. С течением времени температура струи охлаждается до температуры гидратообразования. С ростом вертикальной координаты увеличение скорости струи замедляется, что вызвано увеличением массы струи.

Ключевые слова: затопленная струя, гидрат, миграция углеводородов.

Kildibaeva S.R.

Post graduate student, Sterlitamak branch of Bashkir State University MATHEMATICAL MODEL MULTIPHASE submerged jet

Abstract

In article are considered mathematical model of multiphase submerged jets. Multiphase jet consists of oil and gas bubbles. In this paper we calculated the thermal parameters of the submerged jet, obtained plots of the temperature and velocity of the vertical coordinate. Particular attention is paid to the study of the jet spray cooling process because of the "capture" of the surrounding water. With the passage of time, the temperature of the jet is cooled to a temperature of hydrate formation. With the growth of the vertical coordinate increase in jet velocity slows down, which is caused by the increase in mass of the jet.

Keywords: submerged jet, hydrate, migration of hydrocarbons.

Изучение затопленных струй актуально в связи с увеличением объемов глубоководной добычи углеводородов и проблемой охраны прибрежной и шельфовой зоны морей от загрязнения [6]. Возможные масштабы загрязнения водоемов зависят от подводных течений водоемов, от интенсивности выбросов и их теплофизических характеристик, глубины залегания трубопроводов, размеров и характера разрывов [2]. Моделирование струи используется для создания математических моделей устройств, предназначенных для сбора нефтепродуктов при их разливе [3,4].

При истечении струи в жидкость изменению подвергаются следующие параметры: скорость течения, температура, концентрация примеси, при этом возникают вихри, движущиеся вдоль и поперек потока, и между соседними струями происходит обмен конечными массами. На границе двух струй формируется область конечной толщины с непрерывным распределением скорости, температуры и концентрации примеси.

Рассмотрим следующую задачу. Пусть на дне океана из устья скважины вытекает смесь нефти и газа. Предположим, что известны объемные расходы, теплофизические параметры нефти, газа и окружающей воды. Требуется произвести расчет параметров многофазной струи.

Рассмотрим распределение температуры нефти и газа в струе. Предположим, что капли нефти, пузырьки газа и вода, вовлекаемая в струю, имеют одинаковые температуру, скорость и траекторию в каждом сечении струи. Ось z, направлена вертикально вверх. При расчетах используем методику, предложенную в [7].

nD2

Начальная площадь поперечного сечения струи: S0 = —-— , где D=2B0, B0 - радиус скважины.

4

Q

Скорость течения смеси определим как Vo = —0 , Q0 = Qg + Qg , здесь Q, Q0g - начальные объемные расходы нефти и

Sn

Удельный поток импульса на устье: Mo = Q0Vo .

Сила плавучести на единицу массы представляется в виде: g'o

g ■ Р^р,

Р

где pw - плотность окружающей воды, р - плотность истекающей смеси, g- ускорение свободного падения.

28

Удельный поток плавучести: Fo = g’o • Qt

о ■

Расстояние, на котором удельный импульс силы плавучести превышает удельный импульс, приобретаемый струей вначале:

3/1 3/

/ = м/4/ F/4.

Время для перемещения объема жидкости на расстояние z: t =

Радиус струи на высоте z представляется в виде: B

' 3,4а z

B

3а ^% Г -V*

п

V Fo j

\

+1

B

J

4

4/

• 1'Ъ .

Скорость подъема струи в сечении z: w

2 ■

Объемный расход представим в следующем виде: Q = nwB2 .

Таким образом, можно определить температуру в любом сечении струи, считая, что температура в сечении струи одинаковая и зависит только от координаты z:

т„, = Tw +(T0 -Tw)Q0,

где Tw, T0 - температуры окружающей воды и вытекающей жидкости, Q0 - начальный объемный расход.

Далее представлены зависимости параметров струи от высоты z, полученные на основе численного решения. Для параметров расчета приняты следующие значения: B0=1.7 м, Tw = 4°C, Q0 = 1м3/с, а = 0.149.

На рис. 1 представлена зависимость температуры струи от вертикальной координаты z для начального значения поступающих углеводородов T0=100 °C , T0=80 °C , T0=60 °C . Из графика видно, что на расстоянии z=10 м температура струи остывает до значения Tjgt = 8.0 °C , Tjgt = 7.2 °C , Tjgt = 6.4 °C . При давлении 150 атм. равновесная температура гидратообразования

составляет 21 °C, что говорит о том, до того как температура в струе ниже данной, происходит процесс гидратообразования.

В работе рассмотрена модель затопленной струи, проведены расчеты теплофизических и кинематических параметров, получены графики зависимости температуры струи от вертикальной координаты. Вследствие «захвата» окружающей воды струей происходит охлаждение струи до температуры окружающей среды.

Рис. 1 Зависимость температуры струи от вертикальной координаты

Литература

1. Абрамович Г.Н. Теория турбулентных струй. М.: ЭКОЛИТ, 2011.- 720 с.

2. Гильманов С.А., Шабаев Р.Б. Экспериментальное исследование струй положительной плавучести в воде // Современные проблемы физики и математики. - 2004. - Т. 2. - С. 44-48.

3. Кильдибаева С.Р. Моделирование купола-сепаратора при разливе нефти в шельфе // Фундаментальные исследования. -2013. - № 10 (часть 5). С. 1045-1050.

4. Кильдибаева С.Р., Гималтдинов И.К. Математическая модель наполнения купола нефтью. // Математическое моделирование в естественных науках: материалы XXII Всероссийской школы-конференции молодых ученых и студентов (Пермь, Пермский национальный политехнический университет, 2-5 октября 2013г.), - Пермь, 2013. - С.79-80.

5. Маликов З.М., Стасенко А.Л. Асимптотика затопленной струи и процессы переноса в ней // Труды МФТИ. Аэрогидромеханика. - 2013. - Т. 5, № 2. - С. 59-68.

6. Мельников Н.Н., Калашник А.И. Шельфовые нефтегазовые разработки западного сектора российской Арктики: геодинамические риски и безопасность // Газовая промышленность. - 2011. - № 661. - С. 46-55.

7. Lee J.H.W., Chu V.H. Turbulent jets and plumes - a Lagrangian approach. Kluwer Academic Publishers, 2003. - 390 p.

References

1. Lee J.H.W., Chu V.H. Turbulent jets and plumes - a Lagrangian approach. Kluwer Academic Publishers, 2003. - 390 p.

29