Цифровые модели подводных трубопроводов
в северных морях
В работе представлены результаты исследований авторов по созданию пакета программ ТранSшельф, предназначенного для расчета транспортировки природного газа по морским трубопроводам с учетом специфики арктического шельфа и особенностей замерзания морской воды. Пакет позволяет рассчитать все характеристики потока и динамику оледенения (оттаивания) внешней поверхности трубопровода.
Ресурсы углеводородного сырья в Арктике, по оценке экспертов, составляют треть мировых запасов углеводородов. Разработка газоконденсатных месторождений на арктическом шельфе России является важной стратегической задачей топливно-энергетического комплекса. Ее решение невозможно без развития отечественного научно-технического потенциала.
При освоении арктических месторождений предпочтение часто отдается трубопроводной транспортировке добываемого сырья. Для анализа проектных решений по созданию новых морских трубопроводов, для обеспечения безопасности и надежности их эксплуатации необходимо проведение комплекса научно-исследовательских работ, в который входит непростая задача создания адекватной математической модели подводного трубопровода. Эта модель (цифровой двойник) в режиме реального времени должна давать информацию о процессе транспортировки сырья и о работе оборудования. Знание зависимостей всех характеристик потока от изменяющихся во времени реальных условий позволяет обеспечить безопасность и надежность транспортировки добываемого сырья.
Широкое распространение зарубежных программных комплексов OLGA, PIPESIM, Fluent, Star-CD и других создает иллюзию того, что с их помощью можно решить любую практическую задачу. В действительности это не совсем так, поскольку математические модели и вычислительные алгоритмы, используемые в этих программах, не универсальны. Адекватность математической модели в реальных сложных многофакторных задачах является предметом самостоятельного исследования. Нетривиальным моментом является также выбор эффективного численного метода решения систем уравнений модели, обеспечивающего возможность расчета в режиме реального времени.
В нормативных документах Газпрома подчеркивается, что для каждой практической задачи проектирования и расчета газопровода требуется доработка и коммерческих пакетов, и лицензионных программ с учетом специфики задачи. В последние годы использование зарубежных программ осложнилось из-за введения санкций в отношении России.
В данной статье представлена часть результатов исследований авторов по созданию пакета программ ТранSшельф, позволяющего рассчитать как характеристики смеси газов, транспортируемой по морским трубопроводам, так и динамику оледенения (оттаивания) внешней поверхности трубопроводов в северных морях.
Математическая модель течения газа по трубопроводу:
Наиболее содержательные одномерные математические модели нестационарного неизотермического течения природного газа по трубопроводу постоянного круглого сечения имеют общую структуру. Они состоят из уравнения неразрывности, уравнения движения (баланса импульса), уравнения сохранения полной энергии, уравнения состояния газовой смеси и калорического уравнения. В уравнение движения входит коэффициент гидравлического сопротивления, являющийся функцией расхода газа, внутреннего радиуса трубопровода, коэффициента относительной шероховатости внутренней стенки трубопровода и коэффициента динамической вязкости газа. Уравнения движения и энергии содержат учет рельефа трассы прокладки трубопровода. В уравнение энергии входит осредненное по углу значение радиальной составляющей вектора потока тепла на внутренней поверхности трубопровода. Величина радиальной составляющей вектора потока тепла определяется из решения сопряженной задачи теплообмена между потоком газа и окружающей средой через многослойную боковую стенку трубопровода, на внешней поверхности которого возможно образование слоя льда. Таким образом, поток тепла на внутренней поверхности трубопровода зависит от плотности, толщины, коэффициента теплопроводности и удельной теплоемкости каждого слоя обшивки трубопровода и слоя образовавшегося льда, от температуры окружающей воды и условий обтекания трубопровода, определяющих величину коэффициента теплообмена наружной стенки трубопровода с окружающей водой.
Плотность, давление, температура, скорость, массовые плотности полной и внутренней энергий газовой смеси в одномерной модели являются функциями координаты z, направленной вдоль оси трубопровода, и времени t.
Допустимость одномерного описания процессов, при котором все характеристики являются осредненными по сечению трубопровода, исследовалась. Проводились расчеты профиля скорости в турбулентном потоке сжимаемого газа для полуэмпирических моделей турбулентности Прандтля, Кармана и Новожилова–Павловского. Решение этой задачи позволило оценить погрешность замены двумерной модели течения газа одномерной моделью. Проведенное исследование привело к следующему выводу: при малых числах Маха, при малости относительных уклонов трассы, при больших числах Рейнольдса и при малых значениях коэффициента эквивалентной шероховатости допустим переход к одномерной модели течения.
Система уравнений математической модели течения газа по трубопроводу дополняется начальными и граничными условиями, соответствующими рассматриваемой задаче.
Адекватность математической модели в каждой практической задаче обеспечивается соответствующим выбором уравнения состояния транспортируемой смеси газов в заданном диапазоне изменений давления и температуры потока. Газовые компании корректируют коэффициенты уравнений состояния в соответствии с данными специально проведенных натурных экспериментов. Эти поправки нередко являются коммерческой тайной, охраняемой законодательством.
Адекватность математической модели зависит также от расчета коэффициента гидравлического сопротивления λ и суммарного коэффициента теплопередачи β. При наличии экспериментальных данных о давлении и температуре газа в ряде сечений трубопровода надежнее всего определять коэффициенты λ и β из решения задачи идентификации этих коэффициентов. Методика процедуры идентификации обобщается на произвольное число участков, удовлетворяющих условию постоянства на них λ и β.
Алгоритм решения задачи идентификации коэффициентов λ и β оформлен в виде программного комплекса ТранSшельф-rev, входящего в пакет программ ТранSшельф.
Для замыкания математической модели течения газа по трубопроводу ее необходимо дополнить моделью теплообмена потока газа с окружающей средой.
Математическая модель теплообмена при отсутствии слоя льда:
При сверхвысоких давлениях (более 20 МПа) для потока газа характерны большие числа Рейнольдса: Re ≈108. При этом интенсивность турбулентных пульсаций в потоке приводит к тому, что в радиальном направлении лимитирующей стадией теплообмена с внешней средой является теплопроводность через многослойную стенку трубопровода. В рамках одномерной модели теплообмен с окружающей средой учитывается интегрально.
В правую часть уравнения энергии входит слагаемое, имеющее размерность мощности объемного источника (стока) внутренней энергии; это слагаемое выражается через R – внутренний радиус трубопровода и через q(z, t) – осредненную по углу радиальную составляющую вектора потока тепла на внутренней поверхности трубопровода в z-ом сечении в момент t. Тепловые условия на внешней поверхности трубопровода на расстояниях порядка R допустимо считать неизменными по z и t. Дополнительный пульсационный перенос внутренней энергии газа в направлении оси z пренебрежимо мал по сравнению с конвективным переносом внутренней энергии в этом направлении, поэтому допустимо считать, что q(z, t) зависит от z и t только параметрически через зависимости от z и t температуры газа и внешних условий. Для осесимметричных задач вплоть до момента возникновения льда величина q(z, t) определяется из решения системы нестационарных одномерных уравнений теплопроводности в областях многослойной боковой стенки трубопровода и в области эффективного теплового пограничного слоя воды при соответствующих граничных условиях сопряжения и начальных условиях. В пределах теплового пограничного слоя передача тепла моделируется линейным уравнением теплопроводности. Величина δ* зависит от многих факторов, в частности от донных течений. Если трубопровод заглублен в грунт, для расчета среднего теплового потока от окружающей среды необходимо решать двухмерную задачу теплопроводности, в которой температура является функцией радиуса и угла.
Отсутствие осевой симметрии приближенно можно учесть в рамках одномерной модели за счет выбора эффективного значения параметра δ*, при этом величину δ* надежнее всего определять по экспериментальным данным с помощью процедуры идентификации λ и β и последующего расчета δ* по найденной величине β.
Необходимым условием для начала оледенения внешней поверхности трубопровода является такое падение температуры газа в потоке, при котором температура на внешней поверхности трубопровода становится ниже температуры T фазового перехода «морская вода – лед». Однако выполнение только этого условия еще недостаточно для начала оледенения. Надо, чтобы поток тепла от окружающей воды был строго меньше, чем поток тепла от внешней поверхности трубопровода внутрь его. Это условие заведомо не выполняется в теплых морях, в которых температура воды значительно выше температуры фазового перехода. Иначе обстоит дело в северных морях, например, в некоторых районах Баренцева моря придонная температура воды близка к температуре фазового перехода, равной в соленой морской воде T = 271,236 К.
Математическая модель теплообмена при оледенении трубопровода:
Уравнения теплопроводности в слоях обшивки трубопровода остаются теми же, что и в модели теплообмена при отсутствии льда. Начальным условием для температуры в каждом слое служит распределение температуры в этом слое в момент начала оледенения. Модель дополняется нестационарным одномерным уравнением теплопроводности в области слоя льда, имеющего изменяющуюся толщину y = y(t), условием неизменности температуры на границе фазового перехода и модифицированным уравнением Стефана. В отличие от пресной воды, в соленой морской воде при нарастании льда происходит отток солей из образующегося льда в слой воды, прилегающий к фронту замерзания. Кроме того, теплота плавления, плотность и теплоемкость морского льда, его коэффициент теплопроводности и температура T фазового перехода являются сложными функциями солености морского льда и его температуры. Все это приводит к искажению динамики оледенения, которая имела бы место в пресной воде. Суммарный эффект выражается в замедлении процесса оледенения. В рассматриваемой модели этот эффект учитывается введением в модифицированное уравнение Стефана дополнительного притока тепла к фронту оледенения, пропорционального скорости нарастания льда, коэффициент пропорциональности α, как и дополнительный приток тепла, в пресной воде равен нулю. Коэффициент α является эффективным параметром модели оледенения. Выбор α и выбор остальных теплофизических параметров нарастающего морского льда представляет собой самостоятельную непростую задачу, для решения которой кроме математической модели необходимо наличие экспериментальных данных. Эффективность этой методики была подтверждена при расчетах скорости оледенения поверхностей в соленой воде по экспериментальным данным, предоставленным К.Е. Сазоновым, начальником лаборатории «Морской ледотехники» ФГУП «Крыловский государственный научный центр».
Алгоритм численного решения задач оледенения поверхностей в соленой морской воде оформлен в виде программных комплексов ТранSшельф-Ste-1 и ТранSшельф-Ste-2, входящих в пакет программ ТранSшельф.
Общая модель неустановившегося неизотермического течения природного газа и оледенения морских трубопроводов в северных морях:
Эта модель состоит из уравнения неразрывности, уравнения движения, уравнения энергии, уравнения состояния Редлиха–Квонга и следующего из него уравнения связи внутренней энергии с температурой и плотностью газа, полуэмпирического уравнения Коулбрука–Уайта для расчета коэффициента гидравлического сопротивления λ и одного из блоков расчета радиальной составляющей q(z, t) вектора потока тепла на внутренней поверхности трубопровода. Модель дополняется начальными и граничными условиями.
Проведенный анализ известных численных методов решения одномерных задач газовой динамики привел к выбору для численного решения системы уравнений модели модифицированной явной двухшаговой схемы Лакса–Вендроффа, которая в этих задачах оказалась предпочтительнее других численных схем, как по скорости счета, так и по простоте реализации. Этот алгоритм численного решения системы уравнений модели транспортировки природного газа и оледенения морских трубопроводов в соленой морской воде оформлен в виде программного комплекса ТранSшельф -t.
Для выбора параметров модели, теоретический расчет которых затруднителен, предлагается использовать программный комплекс ТранSшельф-rev.
Выводы:
Результатом проведенного исследования явилось создание вычислительного алгоритма общей модели транспортировки смеси газов по морским трубопроводам, который позволяет рассчитать для неустановившихся режимов все характеристики потока и динамику нарастания (таяния) морского льда на внешней поверхности трубопровода в северных морях.
Расчет включает учет рельефа трассы, условий вдоль трассы, состава газовой смеси, конструкции трубопровода, нестационарности процессов теплообмена.
На основе этого вычислительного алгоритма компанией ООО «Питерсофтваре» был создан пакет программ ТранSшельф для персональных компьютеров, состоящий из семи программ, каждая из которых реализует численное решение определенной задачи.
Пакет ТранSшельф является отечественным программным продуктом, предназначенным для расчета транспортировки сырья по реальным трубопроводам с учетом специфики арктического шельфа и особенностей замерзания соленой морской воды. Он позволяет без обращения к коммерческим пакетам получать в режиме реального времени информацию обо всех характеристиках потока и о состоянии морского трубопровода.
источник: Статья «Цифровые модели подводных трубопроводов в северных морях» опубликована в журнале «Neftegaz.RU» (№1, Январь 2022), авторы: Курбатова Галина Ибрагимовна, Клемешев Владимир Алексеевич, Филиппов Кирилл Борисович