CC BY
46
2013 Механика № 1
УДК 532.546
А.С. Чиглинцева, А.А. Русинов
Бирский филиал Башкирского государственного университета,
Бирск, Россия
К ТЕОРИИ ПРОЦЕССА ГИДРАТООБРАЗОВАНИЯ ГАЗОВЫХ ПУЗЫРЕЙ В УСЛОВИЯХ МИРОВОГО ОКЕАНА
В работе предложена и построена теоретическая модель миграции метановых пузырей в условиях образования гидрата в вертикальном канале, находящегося над газовым источником. Изучено влияние различных параметров на процесс образования гидратных пузырей. Получены минимальные значения массовых расходов газа и воды, необходимых для процесса образования гидрата.
Установлено, что при миграции газовых пузырей в вертикальном канале возможны два режима протекания процесса гидратообразования в зависимости от значения массового расхода воды. Если начальное значение массового расхода воды больше критического, то при данном массовом расходе газа, достаточного для гидратообразования, газовые пузырьки полностью превращаются в гидратные, при этом температура воды в канале не достигает равновесной температуры образования гидрата. Если начальное значение массового расхода воды меньше критического, то газовые пузырьки покрываются гидратной оболочкой, а температура воды в канале достигает равновесной и процесс гидратообразования завершается.
Получено, что с увеличением начального массового расхода воды высота всплытия гидратных пузырьков растет до некоторой максимальной величины, которая соответствует его критическому значению, а при дальнейшем увеличении - падает.
В работе проведен анализ влияния на процесс гидратообразования различных глубин, на которых расположены газовые источники: 800, 1200 и 1500 м. Получено, что если происходит процесс частичного гидратообразования, то чем глубже расположен источник, тем на большую высоту всплывают пузырьки. Если происходит процесс полного образования гидрата, то картина совершенно иная: чем выше расположен источник газа - тем больше высота гидратообразования.
Ключевые слова: газовые пузырьки, гидратообразование, миграция газовых пузырьков, вертикальный канал.
A.S. Chiglintseva, A.A. Rusinov
Birsk branch of the Bashkir State University, Birsk, Russian Federation
TO THE THEORY OF PROCESS HYDRATING OF GAS BUBBLES IN THE CONDITIONS OF THE WORLD OCEAN
The theoretical model of migration metane bubbles in the conditions of formation hydrate is offered and constructed in work. Influence of various parameters on process of formation hydrate bubbles is studied. The minimum values of mass expenses of gas and the water necessary for process hydrating are received.
Found that the migration of gas bubbles in a vertical channel, there are two modes of the process hydrate formation, depending on the mass expenses of water. If the initial value of the mass expenses greater than the critical, then for a given mass expenses of gas sufficient for hydrate formation, gas bubbles become fully hydrated, thus water temperature in the channel doesn't reach the equilibrium temperature of hydrating. If the initial value of the mass expenses is less critical, then gas bubbles are covered hydration shell, and the water temperature in the channel reaches equilibrium, and the process hydrating is completed.
Also found that the increase of the initial mass expenses, the height of ascent of hydrate bubbles grow to a certain critical expenses, and a further increase - falls.
The analysis of different depths at which gas sources are located is carried out: 800, 1200 and 1500 m. It is found that if there is a process of partial hydrating, the deeper the source is located, the greater the height of the rising bubbles. If there is a process full of hydrating, the picture is quite different: the higher up the gas source - the greater the height of hydrate formation.
Keywords: gas bubbles, hydrating, migration of gas bubbles, vertical well.
Введение
Согласно современным акустическим исследованиям в Мировом океане ежегодно происходят непрерывные выбросы метана, которые составляют порядка нескольких сотен тонн газа [1]. Во всем мире насчитывается около 10000 подводных грязевых вулканов, которые хранят до 1010-1012 м метана в газогидратах [1]. Исследования, проведенные в Охотском море на Сахалинском склоне, показали, что выходы метановых пузырей со дна моря сопровождаются образованием гидрата на их поверхности [2, 3]. Это связано с тем, что в процессе всплытия пузыря к поверхности океана в районах аварий, в зависимости от глубины, создаются термобарические условия для образования и разложения гидратной оболочки на поверхности пузыря [4], фазовая диаграмма которых представлена на рис. 1 [5]. Точкой отмечены условия, которые реализуются при эксплуатации вертикального канала на газовом источнике, находящемся на глубине 1500 м.
301
294'
287 ■
280.
г, к г< 33 + во, <1 за \S$ £^
0,0^ гидрат
i > р, МПа
0 5 10 15 20 25
Рис. 1. Фазовая диаграмма
Поскольку в 1 м3 гидрата содержится около 160 м3 газа и 0,8 м3 воды [6], это обстоятельство позволяет рассматривать газогидрат как один из перспективных источников углеводородного сырья. Поэтому одним из возможных способов получения газа из гидратов является сбор гидратных пузырей с помощью купола-ловушки и последующей транспортировки таких систем в специальные реакторы, созданные для их разложения.
1. Основные уравнения
В работе рассматривается процесс миграции систем газовых пузырьков в вертикальном канале в условиях образования гидрата (рис. 2). Согласно предлагаемой схеме над источником газа в условиях океана устанавливается цилиндрический канал, в который снизу поступает вода. Далее, в канале мигрирующие метановые гидратные пузыри попадают в ловушку, которая расположена над потоком.
Рис. 2. Технологическая схема: 1 - вода; 2 - газовый пузырек;
3 - гидратный пузырек; 4 - источник газа; 5 - осадочные породы;
6 - вертикальный канал
Ось ^ направим по оси цилиндрического канала вертикально вверх. Полагаем, что все основные параметры течения трехфазной системы, состоящей из частиц гидрата, воды и газа, однородны по сечению цилиндра. Пузырьки газа поднимаются вверх вдоль оси 2, причем на поверхности этих пузырьков образуется гидрат.
Пусть п - число пузырьков в единице объема. Дроблением
и слипанием гидратных пузырьков, поднимающихся со скоростью ,
будем пренебрегать. Тогда уравнение сохранения числа пузырьков запишется в виде [7]
Snw ) = 0, 5 = пЯ2,
^ ^ в)
где 5 - площадь сечения канала. Здесь и далее нижние индексы Ь,/,т относятся к параметрам гидрата, воды и газа.
Запишем уравнения сохранения масс соответственно для воды, газа и гидрата, при этом будем полагать, что скорости газовой и гид-ратной фазы
жм=_,, ^=Л, (1)
м1 = , = SPg^, мн = ^ ^ = Ч), (2)
+а^ +аь = 1,5 = пЯ2 ),
где Mi,р0,аг-,^(г=Ь,/,g) - массовые расходы, истинные плотности, объемные содержания и скорости фаз; Jl, Jg, Jh - интенсивности расхода воды, газа и образования гидрата.
Г азогидрат является клатратным соединением с массовым содержанием газа О. Поэтому интенсивности расхода воды, газа и образования гидрата связаны следующим образом:
Jl =(1-О)Jh, Jg = ОJh. (3)
Приведенную систему уравнений необходимо дополнить следующими кинематическими соотношениями:
а1 + а^ +аь =1, аg = 3 па3 ^, ^ = - п(4 - а3) ^, (4)
где аг - радиус газовых пузырьков; ак - радиус гидратных пузырьков.
Уравнение импульсов для трехфазного потока в рамках вышеприведенных допущений может быть записано в виде [7]
^ (М№ +(К + Мь ) ^ ) = “5 ^ («/ Р? +ag Pg +аьР0)“ 2^
где тс - сила трения между потоком и стенкой канала, отнесенная
к единице площади его поверхности.
Для замыкания системы уравнений запишем уравнение миграции для газового пузыря относительно жидкости, пренебрегая инерционными эффектами жидкостного потока, которое будет иметь вид [7]
где ть,трг - масса соответственно гидратного пузырька и присоединенная; хъ - сила трения между гидратным пузырьком и потоком жидкости, отнесенная к единице площади его поверхности; п^ - скорость
миграции газового пузыря относительно жидкости.
При записи уравнения сохранения энергии воды в канале будем пренебрегать потерями тепла через стенки канала. Тогда можем записать следующее уравнение баланса тепла, которое будет иметь вид
Здесь Т1 - температура воды в канале; с1 - теплоемкость воды (при постоянном давлении); д - интенсивность теплообмена между жидкостным потоком в канале и поверхностью гидратной оболочки пузыря, отнесенная к единице его площади.
Жидкость будем считать несжимаемой, а газ калорически совершенным.
М1сТ) = я - ^с/Ть, Я = ^ 4™ъ2д.
(5)
р° = сош^ р = .
2. Силовое и тепловое взаимодействия потока со стенками канала и гидратным пузырем
Силу трения между потоком и каналом примем в виде [8]
Здесь выражение для коэффициента гидравлического сопротивления Е,с соответствует формуле Блазиуса для гладких труб; Яе - критерий
Рейнольдса и ц/ - эффективный коэффициент динамической вязкости
для жидкостного потока.
Сила трения между гидратным пузырьком и потоком жидкости может быть записана в виде [8]
Для интенсивности теплового потока между жидкостью в канале и поверхностью пузыря будем использовать формулу [7]
где Рг и № - приведенные числа Прандтля и Нуссельта; Х1 - эффективный коэффициент теплопроводности для жидкости.
3. Описание интенсивности образования гидратных пузырьков
Интенсивность образования гидрата примем в виде [9]
где - интенсивность образования гидрата, отнесенная на единицу
площади поверхности гидратного пузырька.
Примем, что интенсивность образования гидратного пузыря лимитируется интенсивностью отвода тепла от его поверхности жидкостью. При этом температура поверхности гидратного пузыря ТЬ равна
Ш=2 + 0,46Яе°/55 Рг0 3, Рг,
^ X,
равновесной температуре образования гидрата Г, которая соответст-
Здесь /А - удельная теплота образования гидрата.
Из уравнений (1) с учетом (2), (3) и (4) получим выражения, которые выражают законы изменения соответственно радиусам гидрат-ного и газового пузырьков:
4. Критические значения масовых расходов газа и воды, необходимые для полного образования гидратного пузырька
На основе второго уравнения из (5) с учетом (1) и (3) получим интеграл энергии
Здесь с - постоянная, определяемая из условий на входе и на выходе из канала.
Пусть в канал поступает вода с некоторым начальным массовым расходом М1 = М10 при температуре Т = Т10 и газ с массовым расходом
М^ = М^0. Следовательно, на входе ( 2 = 0) и на выходе ( 2 = ) канала
можем записать следующие граничные условия:
Тогда на основе уравнения (7) с учетом (8) получим выражения для определения критических значений массовых расходов воды и газа, необходимые для полного образования гидрата
вует значению давления р для данной глубины (Тк= Г(р)). Для зависимости Т (р) будем использовать выражение
Т (р)=Т„ 0)+№ .
Р(к 0)
(6)
В соответствии с вышесказанным можно записать
^ = л_ {± _ О) ^=_ % а2
МсТ + Мь/ь +(1-О) МсТь = с.
(7)
Мь/Ъ = Мс Т -Т0), М& = ОМ,,
М* =-М^-. (9)
/ с АТО
Здесь М*, М* - критические значения массовых расходов газа и воды.
Таким образом, при миграции газовых пузырей в вертикальном канале возможны два режима протекания процесса гидратообразования. Данный процесс может завершаться при полном переходе газа в гидратное состояние или при достижении температуры воды в канале равновесной температуры образования гидрата для данного значения давления на соответствующей глубине. При этом структура пузырьков может быть различной. Если начальное значение массового расхода воды больше критического М*, то ПРИ данном массовом расходе газа
достаточного для гидратообразования, газовые пузырьки полностью превращаются в гидратные, при этом температура воды в канале не достигает равновесной температуры образования гидрата. Если же начальное значение массового расхода воды меньше критического, то газовые пузырьки покрываются гидратной оболочкой, а температура воды в канале достигает равновесной и процесс гидратообразования завершается.
5. Результаты расчетов
Для параметров, определяющих теплофизические свойства гидрата, воды и газа, были приняты следующие значения: р, = 910 кг/м ;
/, = 5-105 Дж/кг; р/ =1000 кг/м3; Х1 = 0,58 Дж/(м-с-К); с/ = 4200 Дж/(кг-К);
Х§ = 0,03 Дж/(м-с-К); ^ = 0,001 кг/(м-с); Я§ = 520 Дж/(кг-К); О = 0,12.
Для базовых параметров, определяющих геометрию и режим работы канала, принимались следующие численные величины: Я=1 м; М^0 =1 кг/с [10]; То = 277 К; р=15-106 Па; 0 = 0,1 м/с. Соответст-
вующая равновесная температура Т для данного значения давления
согласно формуле (6) равна 290,6 К.
В работе проведен анализ влияния начального значения массового расхода воды на процесс гидратообразования. Поскольку минимальное значение массового расхода воды (М/0 = 76 кг/с) для данного
начального значения массового расхода газа принимается согласно формуле (9), то задавались значения М10 = 50 кг/с, что соответствует
частичному образованию гидрата на пузырьке и М10 =100 кг/с, которое достаточно для полного образования гидрата.
На рис. 3 представлены распределения массовых расходов газа и гидрата, скоростей газовых пузырьков и воды, температуры воды вдоль канала, а также радиусов газового и гидратного пузырька, находящихся на глубине 1500 м с начальными радиусами 1 мм при полном образовании гидрата. Видно, что процесс гидратообразования заканчивается, т.е. газ в пузырьке полностью переходит в гидратное состояние на высоте порядка 0,5 м. При этом температура воды в канале не достигает равновесной температуры образования гидрата. Из рис. 3, б видно, что газовые пузырьки при заданной начальной скорости ускоряются и принимают максимальную скорость на высоте порядка 2,5 см, после чего происходит уменьшение скорости всплытия гидратных пузырьков.
в г
Рис. 3. Распределения массовых расходов газа и гидрата (а); скоростей пузырьков и воды (б); температуры воды (в) вдоль канала; радиусов газового и гидратного
пузырька (г)
На рис.4 представлены распределения массовых расходов газа и гидрата, скоростей газовых пузырьков и воды, температуры воды вдоль канала, а также радиусов газового и гидратного пузырька, находящихся на глубине 1500 м с начальными радиусами 1 мм при частичном образовании гидрата. Видно, что процесс гидратообразования заканчивается на высоте порядка 1 м. вследствие достижения температуры воды в канале равновесной температуры образования гидрата для данного значения давления. При этом газовые пузырьки не полностью переходят в гидратное состояние, а покрываются гидратной оболочкой. Из рис. 4, б видно, что скорость миграции газовых пузырьков устанавливается.
Рис. 4. Распределения массовых расходов газа и гидрата (а); скоростей пузырьков и воды (б); температуры воды вдоль канала (в); радиусов газового и гидратного
пузырька (г)
На рис. 5 представлена зависимость влияния начального массового расхода воды на процесс гидратообразования при различных значениях глубины функционирования газового источника: 1500, 1000 и 800 м, что соответствует значениям глубинного давления: р = 15, 10, 8 МПа.
Видно, что при увеличении начального массового расхода воды высота, на которой завершается процесс гидратообразования, растет до некоторого максимального значения и затем уменьшается. Это связано
с тем, что, пока 0 < М*, процесс гпдратообразованпя завершается
при достижении температуры жидкости Т1 в канале равновесной температуры Т для данного значения давления. Чем больше М(0, тем меньше Т и, следовательно, процесс гидратообразования. При М 0 > М процесс гидратообразования завершается при полном переходе газовых пузырей в гидратное состояние: чем больше М10, тем меньше Т и АТ увеличивается, следовательно, увеличивается вектор потока тепла, а это означает, что гидрат образуется интенсивнее и высота уменьшается.
Чем меньше глубина всплытия метановых пузырей, тем, соответственно, ниже давление и равновесная температура для данной глубины. Поэтому разность температур АТ=Т -Т1 становится меньше, что, в свою очередь, приводит к увеличению критического значения массового расхода воды. Например, для давления р = 15 МПа М* = 76 кг/с,
а для р = 8 МПа согласно (9) М* =140 кг/с. А это означает смещение пика вправо в зависимости от увеличения глубины всплытия.
Из рис. 5 видно, что если М(0 < М* и происходит процесс частичного гидратообразования, то чем больше донное давление, тем на большую высоту поднимается гидрат.
1,8 1,6 1,4 1,2 1,0 0,8
0,6 0,4 0,2
40 60 80 100 120 140 160 180 200
М,, кг/с
*0
Рис. 5. Зависимость процесса завершения гидратообразования от начального массового расхода воды на различных глубинах:
1 - 1500 м; 2 - 1000 м; 3 - 800 м
Если Ml 0 > м и происходит процесс полного гидратообразования,
то чем выше донное давление, тем больше AT и больше интенсивность гидратообразования, поэтому газовые пузырьки переходят в гид-ратное состояние на более меньших высотах zh.
Заключение
В работе получены минимальные значения массовых расходов газа и воды, необходимые для процесса образования гидрата. Изучено влияние различных начальных параметров на процессы миграции и образования гидратных пузырей в вертикальном канале. Установлено, что при миграции газовых пузырей в вертикальном канале возможны два режима протекания процесса гидратообразования. Данный процесс может завершаться либо при полном переходе газа в гидратное состояние, либо при достижении температуры воды в канале равновесной температуры образования гидрата для данного значения давления на соответствующей глубине.
Полученные результаты исследований расширяют теоретические представления о процессе образования гидратных пузырьков на морских глубинах, которые в дальнейшем могут быть использованы при проведении комплекса инженерно-технологических мероприятий по сбору таких систем в условиях Мирового океана и последующего извлечения из них газа.
Библиографический список
1. Methane discharge from a deep-sea submarine mud volcano into the upper water column by gas hydrate-coated methane bubbles / E.J. Sau-tera [et al.] // Earth and Planetary Science Letters. - 2006. - No. 243(3-4). -P. 354-365.
2. Gumerov N.A., Chahine G.L. Dynamics of bubbles in conditions of gas hydrate formation // Fluid Dynamics. - 1992. - No. 5. - P. 664-669.
3. Study on the kinetics of hydrate formation in a bubble column / Y.-T. Luoa, J.-H. Zhua, S.-S. Fanb, G.-J. Chena // Chemical Engineering Science. - 2007. - No. 62. - P. 1000-1009.
4. Шагапов В.Ш., Чиглинцева A.C., Сыртланов B.P. О возможности вымывания газа из газогидратного массива посредством циркуляции теплой воды // Прикладная механика и техническая физика. -2009. - Т. 50, № 4. - С. 100-111.
5. Maksimov A.O., Sosedko E.V. Dynamics of sea bubbles covered by a hydrate skin // XVI Session of the Russian Acoustical Society, Moscow, November 14-18, 2005. - Moscow, 2005. - P. 459-462.
6. Дмитриевский A.H. Газогидраты морей и океанов - источник углеводородов будущего. - М.: ИРЦ Газпром, 2009. - 416 с.
7. Нигматулин Р.И. Динамика многофазных сред. - Т. 1. - М.: Наука, 1987.- 464 с.
8. Химическая гидродинамика: отрав. пособие / А.М. Кутепов,
А.Д. Полянин, З.Д. Запрянов, А.В. Вязьмин [и др.]. - М.: Квантум, 1996. - 336 с.
9. Покусаев Б.Г. Процессы переноса в многофазной среде // Теоретические основы хим. технологии. - 2007. - Т. 41, № 1. - С. 35-43.
10. Исаев В.П. О газовом палеовулканизме на Байкале // Геология нефти и газа. - 2001. - № 5. - C. 45-50.
References
1. 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, no. 243(3-4), pp. 354-365.
2. Gumerov N.A., Chahine G.L. Dynamics of bubbles in conditions of gas hydrate formation. Fluid Dynamics, 1992, no. 5, pp. 664-669.
3. Luoa Y.-T., Zhua J.-H., Fanb S.-S., Chena G.-J. Study on the kinetics of hydrate formation in a bubble column. Chemical Engineering Science, 2007, no.62, pp. 1000-1009.
4. Shagapov V.Sh., Chiglinceva A.S., Syrtlanov V.R. O vozmozhnosti vymyvaniya gaza iz gazogidratnogo massiva posredstvom cirkulyacii teploj vody [On the possibility of washing of gas from gas hydrate solid by circulating warm water]. Prikladnaya mexanika i texnicheskaya fizika, 2009, vol. 50, no. 4, pp. 100-111.
5. Maksimov A.O., Sosedko E.V. Dynamics of sea bubbles covered by a hydrate skin. XVI Session of the Russian Acoustical Society, Moscow, 2005, pp. 459-462.
6. Dmitrievskij A.N. Gazogidraty morej i okeanov - istochnik uglevo-dorodov budushhego [Hydrates the seas and oceans is the future source of hydrocarbons]. Moscow: IRC Gazprom, 2009. 416 p.
7. Nigmatulin R.I. Dinamika mnogofaznyx sred [Dynamics of multiphase flows]. Moscow: Nauka, 1987, vol. 1. 464 p.
В. Kutepov A.M., Polyanin A.D., Zapryanov Z.D., Vyaz’min A.V., Kazenin D.A. Ximicheskaya gidrodinamika: Spravochnoe posobie [Chemical Hydrodynamics: A Reference Guide]. Moscow: Kvantum, 1996. 336 p.
9. Pokusaev B.G. Processy perenosa v mnogofaznoj srede [Transport processes in multiphase environment]. Teoreticheskie osnovy khimicheskoj. tekhnologii, 2GG7, vol. 41, no. 1, pp. 35-43.
1G. Isaev V.P. O gazovom paleovulkanizme na Bajkale [About the gas paleovulkanizme Baikal]. Geologiya nefti i gaza, 2GG1, no. 5, pp. 45-5G.
Об авторах
Чиглинцева Ангелина Сергеевна (Бирск, Россия) - кандидат физико- математических наук, доцент кафедры математического анализа и прикладной математики Бирского филиала Башкирского государственного университета (452451, Башкортостан, г. Бирск, ул. Интернациональная, 1G, e-mail: [email protected]).
Русинов Алексей Александрович (Бирск, Россия) - аспирант кафедры математического анализа и прикладной математики Бирского филиала Башкирского государственного университета (452451, Башкортостан, г. Бирск, ул. Интернациональная, 1G, e-mail: [email protected]).
About the authors
Chiglintseva Angelina Sergeevna (Birsk, Russian Federation) -Ph.D. of Physical and Mathematical Sciences, Ass. Professor, Department of Mathematical Analysis and Applied Mathematics, Birsk Branch of Bashkir State University (12, В March St., 452451, Birsk, Bashkortostan, e-mail: [email protected]).
Rusinov Alex Alexandrovich (Birsk, Russian Federation) -Assistant, Department of Mathematical Analysis and Applied Mathematics, Birsk Branch of Bashkir State University (111, Gagarin St., 452452, Birsk, Bashkortostan, e-mail: [email protected]).
Получено 7.G2.2G13
