Аспекты применения интегрального лагранжевого метода для описания особенностей глубоководного истечения нефтепродуктов при повреждениях морских трубопроводов
Исследование особенностей распространения нефтепродуктов при возможных утечках в воды Мирового океана. Приводится математическая модель, описывающая распространение многофазной затопленной струи, в результате повреждения нефтедобывающей конструкции.
Рубрика | Математика |
Вид | статья |
Язык | русский |
Дата добавления | 14.08.2020 |
Размер файла | 660,0 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Размещено на http://www.allbest.ru/
АСПЕКТЫ ПРИМЕНЕНИЯ ИНТЕГРАЛЬНОГО ЛАГРАНЖЕВОГО МЕТОДА ДЛЯ ОПИСАНИЯ ОСОБЕННОСТЕЙ ГЛУБОКОВОДНОГО ИСТЕЧЕНИЯ НЕФТЕПРОДУКТОВ ПРИ ПОВРЕЖДЕНИЯХ МОРСКИХ ТРУБОПРОВОДОВ
Кильдибаева С.Р.
Стерлитамакский филиал ФГБОУ ВО Башкирский государственный университет, Стерлитамак, Россия
Аннотация. В работе приводится математическая модель, описывающая распространение многофазной затопленной струи, возникновение которой обусловлено повреждением нефтедобывающей конструкции или морского трубопровода, через который транспортировались углеводороды. Исследованы особенности течения многофазных и однофазных (газовых) струй, рассмотрено влияние гидратообразования на параметры струи. Для расчета теплофизических параметров использован интегральный Лагранжевый метод. В результате расчетов получены зависимости температуры и скорости струи от вертикальной координаты. Анализ данных зависимостей позволит спрогнозировать поведение струи, а также исследовать возможность использования устройств для ликвидации утечек.
Ключевые слова: нефтяной разлив, шельф, миграция пузырьков метана, гидратообразование, устранение утечек.
математическая модель утечка нефтепродукты
Abstract. The author uses a mathematical model describing the propagation of a multiphase flooded jet caused by the damage to the oil-producing structure or the marine pipeline used for the transportation of hydrocarbons. The flow features of multiphase and single-phase (gas) jets are studied. The author also studied the effect of hydrate formation on the jet parameters. The integral Lagrange method is used to calculate the thermophysical parameters. As a result of these calculations, the author obtained the dependences of temperature and velocity on the vertical coordinate. Data analysis of dependencies allows predicting jet behavior and study the possibilities of application of devices for leak detection.
Keywords: oil spill, shelf, methane bubble migration, hydrate formation, leak elimination.
Введение. Исследованию особенностей ликвидации глубоководных техногенных разливов уделяется большое внимание в связи с увеличением объемов добычи нефти на шельфе Мирового океана. Данное направление нефтедобычи активно развивается в связи с практически нетронутыми запасами углеводородов, но на практике в случае аварии возрастает опасность разлива нефтепродуктов в водоем. Как показали случаи таких разливов в Мексиканском заливе [1], технологии, доступные на сегодняшний день, обладают рядом негативных аспектов: продолжительное время ликвидации разлива, отсутствие методов сбора углеводородов и учета влияния гидратообразования.
В связи с этим увеличивается интерес исследователей к моделированию устройств для сбора углеводородов при их разливе. Один из способов ликвидации утечки - установка устройства в форме купола над местом утечки. Такое устройство может производить сбор и дальнейшую транспортировку углеводородов. Купол крепится ко дну специальными тросами, к нему подключаются специальные трубки, которые позволяют откачивать нефтепродукты. Более подробная схема купола приводится в работах [2], [3]. Для моделирования работы устройства необходимо рассмотреть распространение струйных течений, характер течения которых зависит от состава струи, начальной скорости углеводородов и окружающей среды, а также температуры углеводородов и термобарических параметров окружающей воды [4].
Целью работы является математическое моделирование многофазных затопленных струй. Анализ динамики изменения теплофизических параметров струи позволит спрогнозировать поведением струи и расширить теоретические знания о распространении затопленных струй.
Актуальность работы связана с необходимостью исследования особенностей распространения нефтепродуктов при возможных утечках в воды Мирового океана. Прогноз миграции нефтепродуктов позволит расширить теоретические знания о процессе распространения углеводородов и уменьшить время устранения возможной утечки.
Постановка задачи и основные допущения
Рассмотрим на дне водоема поврежденную скважину или трубопровод, из которого вытекает смесь нефти и газа. Известны следующие параметры истечения: радиус, состав, температура, плотность, объемный расход. Также известны характеристики окружающей среды. Течение затопленной струи (схема приведена на рис.1) может быть описано с использованием интегрального Лагранжевого метода. Суть метода заключается в следующем: струя представляется в виде элементарных объемов. Каждый такой объем характеризуется значениями температуры, скорости, плотности, координатами в пространстве, а также объемными содержаниями каждого типа углеводорода в струе. Так как параметры элементарного объема соответствуют параметрам струи, то получив зависимости температуры, скорости, плотности, координат и объемных содержаний в контрольном объеме, будем знать все аналогичные параметры затопленной струи.
Подробнее описание метода приведено в работе [4]. Траектория струи может искривляться из-за подводных течений. При моделировании распространения примем, что нефть и газ распространяются соответственно в виде капель и пузырьков. В зависимости от характеристик окружающей среды, где истекает струя, существует два варианта течения струи: с фазовыми превращениями и без. В первом случае температура и давление на глубине распространения струи таковы, что на пузырьках газа начинается процесс образования гидратной корки. Во втором случае значения температуры и давления не соответствуют условиям гидратообразования, и гидратная корка не образуется. В данной работе рассмотрим первый случай, соответствующий глубоководным разливам.
Процесс гидратообразования оказывает существенное влияние на параметры струи: например, так как гидратообразование является экзотермическим процессом, температура струи для двух описанных выше случаев может различаться. Интенсивность выделения тепла напрямую соответствует интенсивности гидратообразования и содержания газа в струе.
Для процесса гидратообразования будет характерна следующая динамика изменения пузырька: сначала он будет полностью газовым, затем на его поверхности начнет появляться гидрат (композитный пузырёк), затем пузырек полностью станет гидратным, что станет завершающим этапом.
Рис. 1 - Схема течения многофазной затопленной струи
Запишем систему уравнений для многофазной струи:
где - масса и плотность компонентов в струе; - объемный расход «захваченной» в струю воды; - объемный расход пузырьков, «покидающих» контрольный объем; - радиус и плотность пузырька; - теплоёмкость и температура КО, - теплоёмкость и температура воды, - теплоёмкость и массовое содержание компоненты в КО. Здесь и далее нижние индексы o - нефть, g - газ, h - гидрат, w - вода.
Система уравнений (1) решается численно для параметров, соответствующих разливу нефти в Мексиканском заливе.
Результаты расчетов
Начальные параметры системы:
На рисунке 2 представлена зависимость температуры от вертикальной координаты. Здесь и далее сплошная линия ? случай, когда струя многофазная (состоит из газа и нефти), а штриховая линия ? случай однофазной газовой струи ( ). Начальная температура струи , окружающей среды . При течении струи происходит турбулентное вовлечение окружающей более холодной воды в струю, поэтому с увеличением вертикальной координаты температура струи уменьшается. Для случая газовой струи уменьшение температуры происходит более быстро, это связано с меньшим объемным содержанием углеводородов в струе.
Рис. 2 - Зависимость температуры струи по вертикальной координате
На рисунке 3 представлен график изменения скорости струи с увеличением вертикальной координаты. Начальное значение соответствует скорости, с которой углеводороды поступают из поврежденного трубопровода. По мере распространения струи происходит «захват» окружающей воды, струя становится более тяжелой. Уменьшение скорости струи связано с утяжелением струи и уменьшением начального импульса струи.
Рис. 3 - Зависимость скорости струи по вертикальной координате
Заключение. В работе приводится уточненная математическая модель течения затопленной струи с учётом образования гидрата. Исследованы особенности течения многофазных и однофазных (газовых) струй, для расчета теплофизических параметров использован интегральный Лагранжевый метод. В результате расчетов получены зависимости температуры и скорости струи от вертикальной координаты.
Полученные в работе зависимости теплофизических параметров позволят спрогнозировать процесс миграции углеводородов, проанализировать какой глубины они достигнут и за какое время, каким будет состав струи и температура. Такие данные о струе являются особенно важными для процесса ликвидации утечек и могут быть использованы для моделирования работы устройств для ликвидации утечек.
Список литературы
1. Akhmetov R.R. Oil spills from offshore drilling and development: causes and effects on plants and animals / Akhmetov R.R., Krainov S.A. // European Science. - 2017. - № 8 (30). - P. 16-21.
2. Gimaltdinov I.K. On the theory of accumulation of hydrocarbons in a dome used to eliminate a technogenic spill at the bottom of the ocean / Gimaltdinov I.K., Kildibaeva S.R. // Journal of Engineering Physics and Thermophysics. - 2018. - Vol. 91. - № 1. - P. 246-251.
3. Gimaltdinov I.K. About the theory of initial stage of oil accumulation in a dome-separator / Gimaltdinov I.K., Kildibaeva S.R. // Thermophysics and Aeromechanics. - 2015. - Vol. 22. № 3. - P. 387-392.
4. Gimaltdinov I.K. Model of a submerged jet accounting for two limiting schemes of hydrate formation / Gimaltdinov I.K., Kildibaeva S.R. // Thermophysics and Aeromechanics. -2018. -Vol. 25. № 1. -P. 75-83.
Размещено на Allbest.ru
...Подобные документы
Исследование понятия "форма" в биологии и векторной геометрии. Математическая модель формообразования и пути познания энергетических процессов в геометрии. Деление отрезка в золотом сечении. Уравнение экспансии как векторная основа формообразования.
реферат [400,8 K], добавлен 20.08.2009Математическая постановка задачи для прямоугольной пластины. Исследование спектра частот при сложных граничных условиях с помощью асимптотического метода. Определение корреляционной функции прогиба пластины. Случайная нагрузка и методы ее описания.
курсовая работа [354,2 K], добавлен 13.11.2016Определение исследования операция как применения научного метода комплексными научными коллективами для решения задач, связанных с управлением организованными (человеко-машинными) системами с целью получения решений. Анализ отличительных особенностей ИСО.
реферат [20,6 K], добавлен 27.06.2011Математическая модель линейной непрерывной многосвязной системы. Уравнение движения и общее решение неоднородной системы линейных дифференциальных уравнений. Сигнальный граф системы и структурная схема. Динамики САУ и определение ее характеристик.
реферат [55,7 K], добавлен 26.01.2009Динамическая модель как теоретическая конструкция, описывающая изменение состояний объекта. Характеристика основных подходов к построению: оптимизационный, описательный. Рассмотрение способов построения математических моделей дискретных объектов.
контрольная работа [769,7 K], добавлен 31.01.2013Идеи интегрального исчисления в работах древних математиков. Особенности метода исчерпывания. История нахождения формулы объема тора Кеплера. Теоретическое обоснование принципа интегрального исчисления (принцип Кавальери). Понятие определенного интеграла.
презентация [1,8 M], добавлен 05.07.2016Рассмотрение особенностей метода построения полного проверяющего теста для недетерминированных автоматов относительно неразделимости для модели "черного ящика" и разработка предложений по его модификации. Исследование условий усечения дерева преемников.
курсовая работа [1,3 M], добавлен 20.08.2010Синтез оптимального управления при осуществлении разворота. Разработка математической модели беспилотных летательных аппаратов. Кинематические уравнения движения центра масс. Разработка алгоритма оптимального управления, результаты моделирования.
курсовая работа [775,3 K], добавлен 16.07.2015Нечёткие системы логического вывода. Исследование основных понятий теории нечетких множеств. Операции над нечёткими множествами. Нечёткие соответствия и отношения. Описания особенностей логических операций: конъюнкции, дизъюнкции, отрицания и импликации.
презентация [191,0 K], добавлен 29.10.2013Словесная, математическая постановка исходной задачи. Исследование математической задачи на корректность. Применение метода экспертных оценок и парных сравнений основных объективных, субъективных факторов, послуживших причиной к поступлению учиться в МАИ.
курсовая работа [145,1 K], добавлен 19.12.2009Ознакомление с историей появления метода золотого сечения. Рассмотрение основных понятий и алгоритма выполнения расчетов. Изучение метода чисел Фибоначчи и его особенностей. Описание примеров реализации метода золотого сечения в программировании.
курсовая работа [416,0 K], добавлен 09.08.2015Разработка логико-формальной модели описания методики изготовления винных изделий. Разделение ингредиентов и продукции на множества. Исследование на рефлексивность, транзитивность, симметричность. Построение графа, матрицы смежности и инцидентности.
контрольная работа [165,2 K], добавлен 07.06.2010Основные этапы математического моделирования - приближенного описания класса явлений или объектов реального мира на языке математики. Методы кодирования информации. Построение устройства, которое позволяет переводить код азбуки Морзе в машинный код.
курсовая работа [507,2 K], добавлен 28.06.2011Рассмотрение особенностей сравнения рядов. Характеристика признаков сходимости Даламбера. Критерий Коши как ряд утверждений в математическом анализе. Анализ геометрической интерпретации интегрального признака. Способы определения сумы числового ряда.
контрольная работа [214,6 K], добавлен 01.03.2013Изучение методов решения уравнений математической физики, которые используются для расчётов распространения тепла, концентрации, волн. Решение уравнения теплопроводности интегро-интерполяционным методом (методом баланса), который применим во всех случаях.
курсовая работа [269,2 K], добавлен 15.11.2010Предмет, методы и понятия математической статистики, ее взаимосвязь с теорией вероятности. Основные понятия выборочного метода. Характеристика эмпирической функции распределения. Понятие гистограммы, принцип ее построения. Выборочное распределение.
учебное пособие [279,6 K], добавлен 24.04.2009Характеристика понятий "порядок", "хаос" и особенностей их применения в точных науках: математике, физике. Исследование взаимосвязи упорядоченных и хаотических явлений и методы формулировки (содержательно и математически строго) правил относительно них.
реферат [595,3 K], добавлен 29.11.2010Литералы рассуждения и вопрос об их отрицаниях. Математическая модель отрицания для рассуждения, содержащего связную совокупность суждений. Отрицания в математической логике и дополнения в алгебре множеств. Интерпретации формул математической логики.
контрольная работа [40,8 K], добавлен 03.09.2010Решение систем уравнений по правилу Крамера, матричным способом, с использованием метода Гаусса. Графическое решение задачи линейного программирования. Составление математической модели закрытой транспортной задачи, решение задачи средствами Excel.
контрольная работа [551,9 K], добавлен 27.08.2009Основные методы формализованного описания и анализа случайных явлений, обработки и анализа результатов физических и численных экспериментов теории вероятности. Основные понятия и аксиомы теории вероятности. Базовые понятия математической статистики.
курс лекций [1,1 M], добавлен 08.04.2011