CC BY
34Добыча, транспорт и переработка нефти и газа
УДК 681.5
ИССЛЕДОВАНИЕ ГИДРОДИНАМИКИ МАГИСТРАЛЬНЫХ НЕФТЕПРОВОДОВ1
А.А. Афиногентов, Ю.А. Багдасарова, Ю.А. Тычинина
Самарский государственный технический университет 443110, г. Самара, ул. Молодогвардейская, 244
Е-mail: pondex@yandex.ru
Рассмотрена математическая модель, описывающая нестационарное движение нефти по магистральным трубопроводам. Проведен анализ пусковых режимов работы нефтепровода. Установлена область максимальных значений амплитуды ударного давления по длине трубопровода.
Ключевые слова: математическое моделирование, объект с распределенными параметрами, краевая задача, трубопроводный транспорт нефти и нефтепродуктов.
Введение. В настоящее время для решения задач технологического расчета трубопровода в системах магистрального транспорта нефти и нефтепродуктов широко используются методики, позволяющие определить требуемые параметры для заранее заданных режимов работы трубопровода [1]. Такой подход основан на приближенных методиках определения основных параметров и не позволяет оценить работу системы в нестационарных режимах, при изменении производительности, а также решать задачи оптимизации режимов работы нефтепроводов.
В статье предлагается подход к решению задач технологического расчета трубопровода на основе решения дифференциальных уравнений движения капельной сжимаемой жидкости в трубах с учетом гидравлического сопротивления, описанных в работах [2-4].
Краевая задача. Взаимосвязь основных параметров (давления P и скорости т), характеризующих движение нефтепродукта плотностью р по трубопроводу постоянного диаметра D длиной L в любой точке I по направлению движения потока и момент времени т, может быть описана системой двух уравнений, в общем случае нелинейных [2], вида
:Работа выполнена в рамках Федеральной целевой программы «Научные и научнопедагогические кадры инновационной России на 2009-2013 гг.».
Александр Александрович Афиногентов, кандидат технических наук, ассистент. Юлия Александровна Багдасарова, кандидат педагогических наук, старший преподаватель.
Юлия Александровна Тычинина, кандидат технических наук, доцент.
дР
~ Р
дт + 2•Б
дР 2 дю
— = с ■ р--------,
дт д/
здесь Р(/, т)= р- g ■
(Н (/,,)-2(1)) , где Н (/, т) - напор жидкости в трубопроводе, г{/) -превышение центра тяжести сечения трубы над произвольной горизонтальной плоскостью, g - ускорение свободного падения, А,(ю(/, т)) - коэффициент гидравлического сопротивления, с - скорость звука в нефти, текущей в стальной трубе с толщиной стенки 5, определяется по формуле Жуковского:
с = , 1 , (2)
Р + Р- °
'К 5 ■ Е
где К - модуль объемного сжатия жидкости (для нефти К~1500 МПа), Е - модуль упругости материала трубы (для стали Е~2-105 МПа).
Положим, что до момента времени т = 0 движение жидкости в трубопроводе
отсутствовало или было стационарным, тогда начальные условия можно предста-
вить в виде
Р(/,0) = Р0(/); (3)
ю(/,0) = ю0. (4)
Сформулируем краевую задачу в следующем виде: к одному концу трубопровода подключен насосный агрегат, изменяющий расход Q(т) нефтепродукта по известному закону от времени, а на другом конце предполагается известным давление (при наличии резервуаров на конечном пункте). Тогда граничные условия могут быть представлены в виде
(0, т)=Щ; (5)
я-Б 2
Р(^ т) = (РаШ +Р-g^) . (6)
Далее будем считать атмосферное давление одинаковым в начальной и конечной точках трубопровода и положим при расчетах РаШ = 0 .
Если полагать в задаче известным закон изменения напора Н(т) (давления), создаваемого насосом, то граничное условие (6) можно заменить:
Р(0, т) = р- g ■ Н н (т). (7)
Коэффициент гидравлического сопротивления А,(ю) определяется в зависимости от режима движения жидкости в трубопроводе, который характеризуется числом Рейнольдса [1]
Ке = —; (8)
при ламинарном режиме течения
, ч 64• V
Яв < 2320, Х(ю) =------------- ; (9)
ю- Б
при турбулентном режиме течения в зоне гидравлически гладких труб по формуле Стокса
2320 < Яе < — ;
^ (10)
Л(ю) = 0,3164 ■
\ 0,25 V 1
(11)
а -Б)
в зоне смешанного трения по формуле Альтшуля
10 „ 500
— < Яе <------;
А« = 0Д1.(,+ ,
в зоне квадратичного трения по формуле Шифринсона
Яе >—, Х(ю)= 0,11 -в°’25 , (12)
в
где V - кинематическая вязкость, г - относительная шероховатость труб [1].
Анализ нестационарных режимов работы трубопровода. Интегрирование системы нелинейных уравнений (1) с начальными условиями (3, 4) и граничными условиями (5, 6) либо (6, 7) проводится численными методами в программе МаШ-САБ с использованием стандартной функции Pdesolve.
Приведем решение рассматриваемой задачи для нестационарного режима работы трубопровода диаметром Б = 1 м, длиной Ь = 100 км = 105 м и профилем трассы
г{/) (рис. 1) при перекачке нефти, имеющей плотность р = 840 кг/м3 , кинематиче-
6 2 /
скую вязкость у = 7,65 -10_ м /с, скорость распространения звука в системе «нефть - труба» с^998 м/с.
Положим, что в начальный момент времени произошло мгновенное изменение напора или подачи. Причем до этого момента движение в трубопроводе отсутствовало, т. е. ю0 = 0 м/с, Р0 (/) = р- g -(Н0 - г(/)) , где Н0 - напор в сечении / = 0 трубопровода. Законы изменения скорости или напора в рассматриваемой задаче будут иметь вид:
Q(0, т) = Qcm = 7000 мУц ; (13)
Нн (0, т) = НСт = 600 м . (14)
Вид граничного условия (13, 14) определяет характер протекания переходного режима; графики изменения скорости движения и давления нефти в трубопроводе представлены на рис. 1.
Постоянной подаче насоса соответствует граничное условие (13), переходный процесс имеет колебательную составляющую в соответствии с рис. 1, а. При постоянном
напоре (14) максимальные значения достигаются в первом колебании, затем проходит
фаза выхода на стационарный режим (рис. 1, б).
Анализ изменения кривых переходных процессов в зависимости от координаты по длине трубопровода представляет отдельный интерес. В обоих случаях изменение
давления происходит с запаздыванием (Тф), пропорциональным расстоянию от начальной точки трубопровода, величина максимального отклонения давления от стационарного возрастает по мере удаления координаты от начала трубопровода, длительность интервала действия максимального давления сокращается. Изменение давления принимает характер гидравлической ударной волны, амплитуда которой может превышать рабочее давление в данной точке трубопровода в несколько раз.
Изменение давления нефти в трубопроводе
Изменение скорости потока нефти в трубопроводе
Распределение напора по длине трубопровода
а б
Рис. 1. Изменение скорости (ю) и давления (Р) нефти во времени (в пяти сечениях по длине трубопровода) и напора (Н) по длине трубопровода (в различные моменты времени) при пуске:
а - Q(0,т)=const; б - H(0,т)=const
Момент максимального отклонения (удара) при этом зависит только от расстояния участка от точки приложения возмущения - начала трубопровода.
Представим амплитуду давления разностью между значением давления, достигаемого в любой точке I трубопровода в произвольный момент времени т в переходном режиме, и значением давления в этой точке в стационарном режиме
Руст (1)-Р(/Д П ) :
А(/,т)= Р(1,х)~ Руст (О,те[0,тя ] , (15)
где хП - длительность нестационарного (переходного) режима работы трубопровода, определяемая как момент времени, в который отклонения давления в любой точке I трубопровода не будут превышать 5 % от значения давления в этой точке при новом стационарном режиме.
На рис. 2 приведено распределение амплитуды давления по длине трубопровода в различные моменты времени нестационарного режима его работы. Время переходного процесса хП составляет порядка 320 сек.
а б
Рис. 2. Распределение давления нефти по длине трубопровода в нестационарном режиме:
а - Q(0,т)=const; б-H(0,т)=const
Амплитуда изменения давления зависит от напора, развиваемого насосным агрегатом. В режиме пуска трубопровода с открытым в атмосферу (резервуар) правым концом на последних нескольких километрах трубопровода происходит затухание амплитуды ударной волны, которая достигает своих максимальных значений на участке (0,8+0,9)^Ь.
Указанные выше наблюдения имеют практическую значимость в связи с тем, что поверочные расчеты по толщине стенки трубопровода, как правило, учитывают только значение рабочего давления на рассматриваемом участке, которое значительно ниже ударного на участках (0,8^0,9)^Ь. Следовательно, разработка методик учета амплитуды и длительности ударного давления при проведении поверочных расчетов при проектировании и эксплуатации нефтепроводов, а также разработка мероприя-
тий по снижению негативных факторов пусковых режимов в системах магистрального транспорта нефти являются актуальными задачами.
Исследование распределения по длине трубопровода значения скорости движения фронта волны давления. В напорном трубопроводе при внезапном изменении скорости движения жидкости (мгновенная остановка или появление движения) возникает волна давления, сопровождающаяся резким понижением либо повышением давления. Ударная волна распространяется по жидкости с постепенным затуханием колебаний.
Анализ полученных зависимостей и графиков показывает, что скорость движения фронта волны давления изменяется в зависимости от протяженности пройденного участка трубопровода. Скорость движения фронта волны давления увеличивается по мере приближения к конечной точке участка трубопровода.
Введем два определения для скорости распространения фронта волны давления на произвольном участке трубопровода: мгновенная и усредненная.
*
Под мгновенной скоростью Уф будем понимать скорость, с которой фронт волны давления, т. е. точка Хф, проходит рассматриваемый участок I трубопровода, т. е.
d хф
УФ
(1 ) =
хф
d х
(16)
Под осредненной скоростью Уф движения фронта волны давления будем понимать отношение пути 1ф = I - ¡0 , пройденного точкой Хф , ко времени Хф , за которое
точка Хф этот путь преодолела:
ф
(1 ) = ¡-¡°-хф
(17)
Полученные зависимости мгновенной и осредненной скоростей движения фронта волны давления Уф от координаты 1 трубопровода представлены на рис. 3.
(I В ДОЛЯХ ОТ V)
а
(I в долях от V)
б
Рис. 3. Зависимость мгновенной УФ (а) и усредненной Уф (б) скорости движения фронта волны давления от координаты 1 трубопровода
Анализ представленных зависимостей позволяет сделать вывод, что волна давления распространяется по трубопроводу с неравномерной скоростью, ускоряясь на одних участках и замедляясь на других, с увеличением осредненной скорости по направлению к конечному сечению трубопровода.
Заключение. В работе рассмотрена краевая задача математического моделирования гидродинамических процессов, происходящих в магистральных нефтепроводах. Полученные в результате решения краевой задачи зависимости изменения скоростей и распределений напора по длине трубопровода в переходных режимах его работы позволили осуществить анализ характера распространения волн давления (скорости и амплитуды) в трубопроводе. Такой подход позволяет в дальнейшем формулировать и решать задачи управления нестационарными режимами работы магистральных нефтепроводов с учетом их волновой динамики.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Типовые расчеты при проектировании и эксплуатации нефтебаз и нефтепроводов / П.И. Тугунов,
B.Ф. Новосёлов, А.А.Коршак, А.М. Шамазов. - Уфа: Дизайн-ПолиграфСервис, 2002.- 658 с.
2. Чарный И.А. Неустановившееся движение реальной жидкости в трубах. - М.: Недра, 1975. - 296 с.
3. Афиногентов А.А., Плешивцева Ю.Э., Снопков А.С. Математическое моделирование управляемых гидродинамических процессов трубопроводного транспорта жидких углеводородов // Вестник Самарского государственного технического ун-та. Сер. Технические науки. - Самара: СамГТУ, 2010. - Вып. 7(28). - С. 137-144.
4. Афиногентов А.А., Плешивцева Ю.Э., Ефимов А.П. Оптимальное по быстродействию управление переходными режимами работы магистрального нефтепровода // Вестник Самарского государственного технического ун-та. Сер. Технические науки. - Самара: СамГТУ, 2011. - Вып. 3(31). -
C. 6-13.
Статья поступила в редакцию 2 марта 2012 г.
HYDRODYNAMIC ANALYSIS OF TRUNK OIL PIPELINES A.A. Afinogentov, Yu.A. Bagdasarova, Yu.A. Tychinina
Samara State Technical University
244, Molodogvardeyskaya str., Samara, 443100
The mathematical model describing non-stationary movement of the trunk oil pipelines is considered. The analysis of starting operating modes of the oil pipeline is made. The area of the maximum values of amplitude of shock pressure on length of the pipeline is established.
Keywords: mathematical simulation, object with distributed parameters, boundary problem, petroleum and petrochemical products transportation in trunk pipelines.
Alexander A. Afinogentov (Ph.D. (Techn.)), Assistant. Yulia A. Bagdasarova (Ph.D. (Techn.)), Senior Lecturer. Yulia A. Tychinina (Ph.D. (Techn.)), Associate professor.
