Середюк М. Д. Вплив зміни швидкості на змішування різносортних рідин у процесі їх послідовного перекачування трубопроводом // Міжнародний науковий журнал "Інтернаука". — 2018. — №20.
Технічні науки
УДК 622.692.4
Середюк Марія Дмитрівна
доктор технічних наук, професор,
професор кафедри газонафтопроводів та газонафтосховищ
Івано-Франківський національний технічний університет нафти і газу
Середюк Мария Дмитриевна
доктор технических наук, профессор,
профессор кафедры газонефтепроводов и газонефтехранилищ
Ивано-Франковский национальный технический университет нефти и газа
Serediuk Mariya
Doctor of Technical Sciences, Professor of the
Department Oil and Gas Pipelines and Storage Facilities
Ivano-Frankivsk National Technical University of Oil and Gas
ВПЛИВ ЗМІНИ ШВИДКОСТІ НА ЗМІШУВАННЯ РІЗНОСОРТНИХ РІДИН У ПРОЦЕСІ ЇХ ПОСЛІДОВНОГО ПЕРЕКАЧУВАННЯ ТРУБОПРОВОДОМ
ВЛИЯНИЕ ИЗМЕНЕНИЯ СКОРОСТИ НА СМЕСЕОБРАЗОВАНИЕ РАЗНОСОРТНЫХ ЖИДКОСТЕЙ В ПРОЦЕССЕ ИХ ПОСЛЕДОВАТЕЛЬНОЙ ПЕРЕКАЧКИ ПО ТРУБОПРОВОДУ
THE FLOW VELOCITY CHANGER INFLUENCE ON THE DIFFERENT KINDS OF LIQUIDS MIXING DURING ITS THE PIPELINES IN SERIES TRANSPORTATION
Анотація. Розроблено метод і програмне забезпечення, які дають змогу врахувати вплив зміни швидкості на інтенсивність сумішоутворення у процесі послідовного перекачування різносортних нафт або нафтопродуктів у трубопроводі. Запропоновано розрахункові формули для розподілу концентрацій рідин довжиною зони суміші, для об’єму суміші та її розкладання в кінцевому пункті трубопроводу.
Ключові слова: нафтопровід, послідовне перекачування, зміна швидкості, рівняння турбулентної дифузії, ефективний коефіцієнт змішування, об’єм суміші
Аннотация. Разработан метод и программное обеспечение, позволяющие учитывать влияние изменения скорости на интенсивность смесеобразования в процессе последовательной перекачки разносортных нефтей или нефтепродуктов по трубопроводу. Предложены расчетные формулы для распределения концентраций жидкостей по длине зоны смешения, для объема смеси и ее раскладки на конечном пункте трубопровода.
Ключевые слова: нефтепровод, последовательная перекачка, изменение скорости, уравнение турбулентной диффузии, эффективный коэффициент смешения, объем смеси.
Summary. The methodology and computer software are provided, which can to indicate the influence of flow velocity changes on the intensiveness of mixture creation during different sorts of oil and oil products transportation by pipelines in series. The applied formulas are proposed to estimate the concentration profiles of mixture components along the mixture zone length, volume of mixture and to characterize the process of it's separation at the end-point of the pipeline.
Key words: pipeline, transportation in series, velocity change, equation of turbulent diffusion, the effective coefficient of mixing, volume of mixture.
Вступ. Магістральні нафтопроводи України сьогодні працюють зі значним недовантаженням, що суттєво знижує енергоефективність їх експлуатації. Одним із шляхів збільшення завантаження нафтопроводів і зменшення енерговитратності трубопровідного транспорту є впровадження прогресивної технології послідовного перекачування різносортних нафт. Зазначена технологія широко застосовується у світовій практиці як при трубопровідному транспорті різносортних нафтопродуктів, так і при транспортуванні нафт різних сортів.
В Україні послідовне перекачування світлих нафтопродуктів в кінці двадцятого століття здійснювалось системою магістральних нафтопродуктопроводів. За останні роки на магістральних нафтопроводах України проведено низку успішних послідовних перекачувань партій нафт з альтернативних джерел постачання.
При значних перевагах технології послідовного перекачування вона має недолік, а саме обов’язкове змішування різносортних рідин у зоні їх контакту у трубопроводі. Об’єм суміші різносортних рідин залежить від багатьох чинників, у тому числі від режимних параметрів експлуатації трубопроводу. Основним завданням при реалізації послідовного перекачування різносортних нафт або нафтопродуктів є досягнення мінімального об’єму суміші. Це визначає теоретичне та практичне значення дослідження закономірностей сумішоутворення різносортних рідин у трубопроводі за їх послідовного перекачування [1; 2].
Аналіз досліджень і публікацій. Дослідженням впливу різних чинників на процес сумішоутворення різносортних рідин присвячено низку робіт. Так у роботах [3; 4] встановлено вплив зміни температури послідовно транспортованих нафт на процес утворення суміші. Роботи [5; 6] присвячено дослідженням ефективності застосування рідинних розділювачів з метою покращення якості суміші і полегшення її розкладання у кінці трубопроводу. У роботах [7; 8] розроблено методи визначення пропускної здатності нафтопроводів за неусталених режимів послідовного перекачування різносортних нафт.
Інтенсивність сумішоутворення послідовно транспортованих рідин суттєво залежить від режиму їх руху в трубопроводі, що визначається значенням числа Рейнольдса. За сталого діаметра трубопроводу та в’язкості транспортованих рідин режим руху залежить від витрати, що пропорційна швидкості руху рідин. У процесі пересування зони контакту різносортних рідин по довжині трубопроводу швидкість руху є величиною змінною. Це впливає на розподіл концентрацій різносортних рідин по довжині зони суміші, а отже, і на об’єм суміші. Не дивлячись на значну кількість робіт щодо реалізації технології послідовного перекачування особливості сумішоутворення різносортних рідин зі змінною швидкістю руху в трубопроводі вимагають додаткового розгляду.
Метою досліджень є встановлення впливу зміни швидкості руху послідовно транспортованих рідин на закономірності їх сумішоутворення у трубопроводі.
Для досягнення мети необхідно розробити метод і програмне забезпечення для розрахунку з урахуванням динаміки зміни швидкості руху в трубопроводі таких параметрів:
Виклад основного матеріалу. Приймаємо, що на початку процесу послідовного перекачування вся порожнина трубопроводу заповнена однією із транспортованих рідин. Позначимо першу у контакті рідину – рідиною А. У певним момент часу розпочинається закачування у трубопровід другої рідини – рідини B. У зоні їх контакту утворюється суміш двох рідин, об’єм якої зростає при переміщенні від початку до кінця трубопроводу.
Процес змішування різносортних рідин при послідовному перекачуванні трубопроводом зазвичай описують рівнянням одномірної турбулентної дифузії, яке відносно рухомої системи координат, що переміщається зі сталою швидкістю потоку, має вигляд
(1)
де - миттєва концентрація другої у контакті рідини;
- час змішування, що відлічується від моменту зміни сортів рідин на початку трубопроводу;
- відстань від початку рухомої системи координат до перерізу, в якому необхідно знайти концентрацію рідин;
- ефективний коефіцієнт змішування рідин.
Як засвідчують результати теоретичних досліджень [7-9], процес зміни послідовно транспортованих рідин у трубопроводі при їх послідовному перекачуванні супроводжується помітною зміною швидкості. Нами розроблені метод та комп’ютерна програма, які враховують зазначений фактор і дають змогу визначити пропускну здатність магістрального нафтопроводу та питомі витрати електроенергії при реалізації технології послідовного перекачування нафт різних сортів. Методика включає такі елементи:
Блок гідравлічного розрахунку лінійної частини нафтопроводу з урахуванням переміщення зони контакту двох сортів нафти базується на застосуванні рівняння балансу тисків та модифікованої формули Колбрука для обчислення коефіцієнта гідравлічного опору [10].
Блок розрахунку режиму роботи нафтоперекачувальної станції (НПС) включає запропоновані нами математичні моделі для напірної характеристики та кривої коефіцієнта корисної дії (ККД) нафтових насосів за довільного значення обертової частоти ротора насоса [11; 12].
Для модельного нафтопроводу довжиною 100 км з внутрішнім діаметром 0,702 м досліджено гідродинамічні параметри послідовного перекачування двох сортів нафти з різними значеннями розрахункової густини і в’язкості.
За результатами багатоваріантних розрахунків для різних положень зони контакту двох сортів нафти на трасі нафтопроводу одержано графічні залежності пропускної здатності та питомих витрат електроенергії на реалізацію послідовного перекачування двох сортів нафти. Розрахунки засвідчили що переміщення зони контакту двох сортів нафти порожниною нафтопроводу спричинює помітну зміну пропускної здатності, яка може досягати 15-20 %. Залежність пропускної здатності нафтопроводу від положення зони контакту різносортних нафт зазвичай близька до лінійної.
Таким чином, методом математичного моделювання доведено, що за умов послідовного перекачування різносортних нафт магістральними нафтопроводами України залежність між витратою нафти та лінійною координатою розміщення зони контакту нафт адекватно описується так:
(2)
де - коефіцієнти математичної моделі, значення яких залежать від геометричних характеристик трубопроводу, особливостей траси, характеристик насосного обладнання, властивостей послідовно транспортованих рідин тощо; відстань від початку трубопроводу до середини зони суміші різносортних рідин.
Коефіцієнт у формулі (2) відповідає пропускній здатності нафтопроводу за умови перекачування рідини А.
Помітна зміна робочої витрати, а отже і швидкості руху в трубопроводі у процесі переміщення зони контакту різносортних рідин спричинює зміну інтенсивності утворення суміші.
Для випадку змінної у часі швидкості руху послідовно транспортованих трубопроводом рідин рівняння одномірної турбулентної дифузії набуває вигляду [1-2; 9]
(3)
де - відстань від початку трубопроводу до перерізу, в якому необхідно знайти концентрацію послідовно транспортованих рідин.
Із рівняння нерозривності потоку рідини можна записати:
(4)
де - внутрішній діаметр трубопроводу.
Формула (4) з урахуванням виразу (2) набуває вигляду
(5)
Розв’язуємо диференціальне рівняння (5), у результаті одержуємо аналітичну залежність між часом переміщення зони контакту двох різносортних рідин та відстанню від початку трубопроводу до середини зони суміші
(6)
Виразимо відстань від початку трубопроводу до середини зони суміші послідовно транспортованих нафт як функцію часу їх змішування у трубопроводі
(7)
Використовуючи рівняння (4) з урахуванням виразу (7), одержуємо закономірність зміни швидкості руху послідовно транспортованих рідин як функцію часу переміщення середини зони суміші двох різносортних рідин
(8)
Вводимо позначення
(9)
(10)
З урахуванням (9) і (10) рівняння турбулентної дифузії набуває вигляду
(11)
Шляхом уведення нових змінних
(12)
(13)
зводимо диференціальне рівняння (11) до такого вигляду
(14)
Розв’язок рівняння (14) має вигляд [1; 2]
(15)
де - інтеграл ймовірностей;
(16)
- аргумент інтегралу ймовірностей,
(17)
Для ефективного коефіцієнта змішування застосуємо запропоновану нами у роботі [3] узагальнену математичну модель, в яку вводимо швидкість руху замість робочої витрати
(18)
де - кінематична в’язкість суміші послідовно транспортованих рідин у співвідношенні 1:1 за об’ємом;
- коефіцієнти математичної моделі.
Перекачування різносортних нафт відбувається зазвичай у зоні гідравлічно гладких труб турбулентного режиму. У цьому випадку маємо такі значення коефіцієнтів:
Малов’язкі нафтопродукти транспортуються трубопроводами зазвичай у зоні змішаного тертя турбулентного режиму. У цьому випадку будемо мати такі значення коефіцієнтів математичної моделі (18)
де - абсолютна еквівалентна шорсткість внутрішньої поверхні труби.
Із формули ( 10) випливає, що аналітичний вираз для змінної у рівнянні турбулентної дифузії співпадає з виразом для відстані від початку трубопроводу до середини зони суміші різносортних рідин. У результаті підстановки формули (7) у вираз (8) після математичних перетворень одержуємо швидкість руху рідин як функцію змінної
(19)
З урахуванням (18) і (19) вираз для визначення змінної в рівнянні турбулентної дифузії (14) набуває вигляд
(20)
Після інтегрування виразу (20), одержуємо
(21)
Для випадку послідовного перекачування різносортних нафт в зоні гідравлічно гладких труб формула (21) матиме вигляд
(22)
Закономірності розподілу концентрацій послідовно транспортованих рідин в трубопроводі дають змогу одержати формулу для об’єму суміші.
Позначимо відстань від рухомої системи координат до довільного перерізу зони суміші різносортних рідин. Тоді довжина суміші рідин з урахуванням знаку при за певних граничних концентрацій дорівнює
(23)
Враховуючи вираз (19), одержуємо формулу для визначення об’єму суміші різносортних рідин за змінної швидкості їх руху у трубопроводі
(24)
де - аргументи інтеграла ймовірностей, що відповідають граничним концентраціям на початку і в кінці суміші рідин.
Для вирішення питання розкладання суміші послідовного транспортованих рідин у кінці трубопроводу необхідно знати, який об’єм домішок (рідин іншого сорту) міститься у суміші і надійде у резервуари з товарними продуктами.
Об’єм рідини B, що надійде у вигляді домішки у резервуари з рідиною А, у випадку розділення суміші на дві частини в перерізі трубопроводу з концентрацією рідини В знайдемо шляхом інтегрування такого виразу:
(25)
З урахуванням виразів (15) і (23), після інтегрування одержуємо
(26)
Запропонований метод урахування впливу зміни швидкості руху послідовно транспортованих рідин на інтенсивність їх сумішоутворення реалізований у комп’ютерній програмі. Програма дає змогу будувати криву розподілу концентрацій з урахуванням всіх впливових чинників для довільного положення зони контакту різносортних рідин в трубопроводі. Розраховуються об’єм суміші різносортних рідин, а також об’єми домішок, які надходять у резервуари з товарними продуктами на кінцевому пункті трубопроводу.
З метою апробації методики виконані розрахунки параметрів послідовного перекачування двох сортів нафти з такими розрахунковими властивостями:
Як модельний нафтопровід вибрана ділянка діючого магістрального нафтопроводу довжиною 100 км з внутрішнім діаметром 0,702 м. Нафтоперекачувальна станція оснащена магістральними насосами марки НМ 3600-230.
За результатами багатоваріантних розрахунків одержано, що пропускна здатність модельного нафтопроводу ( м3/с) так залежить від лінійної координати зони контакту нафт (м):
(27)
За комп’ютерною програмою проведено розрахунки концентрації другої у циклі нафти для чотирьох положень зони контакту на трасі нафтопроводу та різних значень відстаней від початку рухомої системи координат.
За результатами обчислень побудовано графіки розподілу концентрацій нафти В по довжині зони суміші (рисунок 1).
Визначено об’єм суміші різносортних нафт для граничних концентрацій = 1–99 % для різних положень зони контакту на трасі нафтопроводу (рисунок 2).
Обчислено об’єм нафти В, що надійде у резервуари з нафтою А як домішка для різних варіантів розподілу суміші на дві частини у кінці нафтопроводу (рис. 3).
Для порівняння виконано розрахунки параметрів сумішоутворення за зазначених вище умов без урахування зміни швидкості руху послідовно транспортованих рідин. У результаті одержано, що урахування зміни швидкості дає змогу на 5-10 % уточнити як об’єм суміші, так і об’єми рідин, які як домішка надходять у резервуари з товарними продуктами на кінцевому пункті трубопроводу.
Рис. 1. Зміна концентрації нафти сорту В по довжині зони змішування для різних положень початку рухомої системи координат
Рис. 2. Залежність об’єму суміші послідовно транспортованих нафт від розташування середини зони змішування у нафтопроводі
Рис. 3. Залежність об’єму нафти сорту В, який як домішка надходить у резервуари з нафтою А, від концентрації розділення суміші на дві частини в кінцевому пункті нафтопроводу
Висновки
Література