Новый метод расчета многопродуктовых магистральных трубопроводов Текст научной статьи по специальности «Химические технологии»

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

УДК 622.691.4

А.С. Дидковская1, e-mail: didkovskaya.a@gubkin.ru

1 Федеральное государственное автономное образовательное учреждение высшего образования «Российский государственный университет нефти и газа (Национальный исследовательский университет) имени И.М. Губкина» (Москва, Россия).

Новый метод расчета многопродуктовых магистральных трубопроводов

Особенность эксплуатации многопродуктовых трубопроводов или трубопроводов последовательной перекачки состоит в том, что по одной трубе перекачивают несколько продуктов, различающихся по свойствам (составу). Такой метод широко используется для транспортирования по трубам светлых нефтепродуктов, а также некоторых видов газопродуктов, которые перекачиваются в жидком состоянии (например, широкой фракции легких углеводородов). При этом перекачка осуществляется «прямым контактированием»: каждый следующий продукт вытесняет предыдущий и, в свою очередь, вытесняется последующим. В зонах контактирования продуктов происходит их смесеобразование, требующее постоянного контроля: необходимо знать не только протяженность образовавшейся смеси, но и то, как меняется ее состав по длине трубопровода.

В статье представлены новые алгоритмы, позволяющие рассчитать концентрацию отельных компонентов нефте-и газопродуктов в любом сечении многопродуктового трубопровода. Научная новизна методики состоит в том, что она позволяет рассчитывать смесеобразование многокомпонентых продуктов, используя произвольное начальное распределение концентрации каждого продукта в смеси, а также учитывать множественные путевые подкачки продукта произвольного состава. Приведены окончательные уравнения для расчета многопродуктового трубопровода, подробно изложено, как эти уравнения получены. Кроме того, демонстрируется применение предлагаемой методики расчета на конкретном примере трубопровода, транспортирующего широкую фракцию легких углеводородов, в котором закачиваемый продукт имеет переменный состав, в том числе за счет путевой подкачки. По результатам расчета построены графики, отражающие компонентный состав широкой фракции легких углеводородов в пункте закачки в трубопровод, пункте подкачки и в конечном пункте трубопровода, что позволяет наглядно сравнить изменение состава продукта в процессе перекачки.

Ключевые слова: нефтепродукт, нефтепродуктопровод, последовательная перекачка, газопродукт, широкая фракция легких углеводородов, путевые подкачки, смесеобразование, концентрация, компонентный состав, математическое моделирование.

A.S. Didkovskaya1, e-mail: didkovskaya.a@gubkin.ru

1 National University of Oil and Gas "Gubkin University" (Moscow, Russia).

Current Analysis of Multi-Product Trunk Pipelines

The feature of multi-product or products pipeline operation is that a single pipe conveys several products different in composition. This technique is widely used for piping of light oil products, as well as some gas products piped in liquid state (e.g. wide light hydrocarbon fraction). Pumping is carried out by "direct contact": each following product displaces a preceding one, and in turn, is displaced by the subsequent. Besides, mixing of products takes place in their contact zones and this process requires monitoring: the spread of the mixture formed is not the only thing one must know, its behavior with changes of pipeline lenght is also important.

The article presented new algorithms which enable to compute the concentration of some components in oil-and-gas products in any multi-product pipeline passage. The scientific novelty of the procedure is that it allows to compute mixing of multi-component products using random initial distribution of each product concentration in the mixture, as well as to take into account multiple way injections of a product of arbitrary composition. The article gives final equations to design a multi-product pipeline and detailed description of the way these equations were derived. Besides it demonstrates the application of the proposed design procedure by the example of the pipeline for wide light hydrocarbon fraction in which composition of an injected product is variable due to way injections as well. Analytic data were used to construct graphs reflecting blend composition of wide light hydrocarbon fraction at the stations of pumping into the pipeline, and of way injections, as well as at the pipeline terminal enabling to compare product behavior in pumping visually.

Keywords: oil product, oil-products pipeline, products pipeline operation, gas product, wide light hydrocarbon fraction, way injections, mixing, concentration, blend composition, mathematical modeling.

68

№ 9 сентябрь 2018 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS TRANSPORTATION AND STORAGE

Технология последовательной перекачки светлых нефтепродуктов, когда партии разных продуктов перекачиваются по трубе одна за другой, и методика расчета смесеобразования в зоне контакта партий достаточно широко освещались в научной литературе [1-5]. Примером использования многопродуктового трубопровода является транспортировка широкой фракции легких углеводородов (ШФЛУ) - продукта первичной переработки попутного нефтяного газа и нефтегазового конденсата, состоящего из различных высших углеводородов. По аналогии с термином «нефтепродукт» для таких продуктов, как ШФЛУ, применим термин «газопродукт», хотя он и не настолько популярен. ШФЛУ транспортируют по магистральным трубопроводам в жидком агрегатном состоянии. На сегодняшний день в России лидером по протяженности и объемам перекачки ШФЛУ является ПАО «СИБУР Холдинг». Несколько меньшие объемы транспортируют ПАО «Газпром», ПАО «НК «Роснефть», ПАО «Сургутнефтегаз» и ПАО «Татнефть». Типичная ШФЛУ имеет следующий состав: метан-этановая (СН4 - С2Н6) фракция - 2-5 %, пропан-бутановая (С3Н8 - СД0) фракция - 20-40 %, гек-сановая (С6Н14) фракция - 5-30 %, на пентановую (С5Н12) фракцию приходится остаток. При этом конкретный состав ШФЛУ может сильно варьироваться в зависимости от состава исходного сырья и метода переработки. Если бы транспортируемая ШФЛУ имела постоянный, т. е. не изменяющийся во времени компонентный состав, то, естественно, продукт, поступающий в пункт приема в конце трубопровода, имел бы тот же самый состав. Однако массовые концентрации компонентов ШФЛУ в пункте закачки в трубопровод постоянно варьируются, и взаимное вытеснение порций ШФЛУ различного состава происходит в условиях непрерывного путевого перемешивания. Кроме того, на трубопроводах ШФЛУ осуществляются путевые подкачки продукта.

Компонентный состав определяет давление упругости насыщенных паров ШФЛУ, которое имеет ключевое значение для выбора безопасных режимов ее перекачки по трубопроводу. Также в зависимости от состава ШФЛУ регулируется технологический процесс ее переработки. Поэтому возникает вопрос, как по данным о компонентном составе продукта в пункте закачки выяснить, как изменится его состав в пункте приема, находящемся зачастую на значительном расстоянии.

ОСНОВНЫЕ УРАВНЕНИЯ МОДЕЛИ

Согласно теории продольного перемешивания осредненная по сечению концентрация с(х, ^ любой «пассивной» примеси в потоке жидкости в трубе удовлетворяет одномерному уравнению типа уравнения теплопроводности:

öc + i) дс-кд2с

(1)

сг, при 0 < t < t2; с2, npnt^i«^;

cn,npvitn_1<t<tn.

(2иx-v2t)

c(x,t) = е sr^0(x, t).

(2)

Тогда слагаемые уравнения (1) можно записать через функцию 0(х, t):

и

4 К

D (2DX - uzi)

дС = dt дс =

дх

Э2с = Эх2 4/С2'

(2ux-uzt)

4К

0(х, t) + e

(2«-ift) 00 "эр

4-К

г к

и'

4 К

(2мЭ0

0(х, t) + в « -дх,

(2их - и2£)

В 4К ■

в(Х, t) +

V 30 d20

К дХ дХ*

Подставляя полученные соотношения в (1) и приводя подобные слагаемые, получим для функции 0(х, t) уравнение:

dt Эх2.

(3)

Таким образом, наша задача свелась к необходимости найти решение уравнения (3), удовлетворяющее при х = 0 краевому условию:

где и - средняя по сечению скорость перекачки, м/с; х - координата вдоль оси трубопровода; t - время, с; К - эффективный коэффициент продольного перемешивания, м2/с. Решим краевую задачу. Пусть при х = 0 (в начале трубопровода) известно, как изменяется концентрация с() любого продукта в смеси с течением времени I > 0. Граничная функция с0(£) может быть задана в виде многоступенчатой функции:

u2t

0o(i) = e«-co(t).

(4)

Искомое решение может быть представлено в виде [4]:

К

2Vti 0 [K-(t - т)]3/2

* e"«(t-x)-0 (т)-£/т.

(5)

Функция (5) называется интегралом Дюгамеля.

Вернемся к неизвестной функции с(х, £). Заменим в уравнении (5) функцию 0(х, t) согласно зависимости (2), а 00(т) подставим из (4). Получаем:

Требуется найти концентрацию с(х, t) этого продукта в любом сечении х по длине трубопровода в любой момент времени, т. е. найти решение уравнения (1) в области 0 « х « L, t > 0. В частности, требуется найти изменение концентрации продукта с(Ь) с течением времени в конечном сечении х = L трубопровода. Введем новую функцию 0(х, V) согласно равенству:

2VÄ 0 [K-(t - т)]3/2 *е «(t-г) -сДт^-от.

(6)

Введем безразмерный параметр Пекле Pe = dL/K и перейдем к безразмерным координате х, = x/L и времени t, = t/T процесса, где T = L/d - время движения жидкости по трубопроводу, с. Тогда уравнение (6) запишется в безразмерном виде:

Для цитирования (for citation):

Дидковская А.С. Новый метод расчета многопродуктовых магистральных трубопроводов // Территория «НЕФТЕГАЗ». 2018. № 9. С. 68-73. Didkovskaya A.S. Current Analysis of Multi-Product Trunk Pipelines. Territorija «NEFTEGAS» = Oil and Gas Territory, 2018, No. 9, P. 68-73. (In Russ.)

TERRITORIJA NEFTEGAS - OIL AND GAS TERRITORY No. 9 September 2018

69

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

с(х n = VPe Р*_, i

*e * (f.-ч) -cjiJ-OT.

(7)

По формуле (7) можно найти распределение концентрации любого из компонентов смеси по длине трубопровода (0 « x, « 1) в любой момент времени

(Ь > 0).

Если необходимо найти зависимость cL(t)) концентрации продукта в смеси, приходящей на конечный пункт трубопровода, т. е. в сечении x = L, от времени, то в формулу (7) нужно подставить x) = 1:

"2VitJo (í,-T)3/2

Pe (t.-t-l)2

(t.-т) -c0(x)-dx.

(8)

Формула (8) позволяет получить закон €(:) изменения концентрации нефтепродукта в смеси от времени в конце трубопровода в зависимости от закона с0(^ изменения концентрации того же нефтепродукта в смеси в начале трубопровода.

Численное решение конкретных задач с использованием формул (7)-(8) может быть найдено путем составления несложных компьютерных программ. Проводилось сравнение результатов расчета смесеобразования по предлагаемому алгоритму и «классическим» формулам. Был осуществлен расчет длины смеси, приходящей на конечный пункт нефтепродуктопровода последовательной перекачки при условии, что в начальном пункте распределение концентрации имеет вид:

/хчГО, 0 <t<T; ti, t>T.

То есть, как это принято в «классической» теории последовательной перекачки, считалось, что переключение резервуаров с одного продукта на другой происходит мгновенно в момент времени t = Т. При этом достигнута очень хорошая сходимость результатов по двум методикам.

В то же время возможности новой методики гораздо шире, чем «классической». Она позволяет рассчитывать смесеобразование при последовательной перекачке при абсолютно любом

начальном распределении с0(£), в том числе пригодна для расчета компонентного состава ШФЛУ, поступающей на приемные пункты трубопроводной сети.

Как указывалось, состав ШФЛУ, закачиваемой в трубопровод, претерпевает изменения. Можно считать, что ее транспортируют по трубопроводу небольшими партиями, отличающимися друг от друга фракционным составом. Таким образом, получаем аналогию с последовательной перекачкой светлых нефтепродуктов. И, как при последовательной перекачке, в результате процессов конвекции и турбулентной диффузии в области контакта партий ШФЛУ возникает смесь. Протяженность этой смеси непрерывно возрастает от начала к концу трубопровода. Естественно, в области смеси фракционный состав ШФЛУ плавно изменяется от состава одной партии к составу другой. Поэтому состав газопродуктов, поступающих на приемные пункты, определяется, как и в случае последовательной перекачки нефтепродуктов, исходным составом партий, протяженностью трубопроводов и режимами перекачки. На выходе из нефтеперекачивающей станции обычно устанавливается поточный хроматограф, регистрирующий компонентный состав ШФЛУ в начальном сечении x = 0 продуктопровода. При этом для определения компонентного состава ШФЛУ принято использовать термин «процентное содержание» (процентное содержание ps(x, ^ 5-го компонента ШФЛУ может быть записано через концентрацию с (х, ^ этого компонента как р5(х, ^ = с(х, ^.100 %).

МЕТОД РЕШЕНИЯ

Итак, известна функция р0 (£) изменения процентного содержания 5-го углеводородного компонента ШФЛУ в начальном сечении в зависимости от времени (5 = 1, 2, 3, ..., т, где т - число компонентов ШФЛУ). Также известны плотность р, кг/м3, и кинематическая вязкость V продукта, м2/с, диаметр d, м, и протяженность L, м, магистрального трубопровода,расход й0 перекачки, м3/ч. Задача состоит в том, чтобы спрогнозировать процентное содержание

ps(x, t) каждого компонента ШФЛУ в трубопроводе в зависимости от координаты и момента времени. Возможно, сечение x=L1 является промежуточным сечением трубопровода, в котором проводят подкачку с расходом Q1 ШФЛУ другого состава. Процентное содержание ks (t) каждого компонента в подкачиваемом ШФЛУ известно. Как в этом случае определить состав ШФЛУ, поступающей на конечный пункт x = L трубопровода?

По формуле(8) можно определить процентное содержание p1 (t) компонентов ШФЛУ, поступающих в пункт смешения x = L'

VPe^ 1

1)3**

Peo(t.-i-l)2

' (t.-z) -p (z)-ch:. (9)

Здесь безразмерное время t, = t/T1, где T = Ly^ - время прохождения ШФЛУ участка трубопровода от начального пункта до пункта подкачки, с; d0 = 4Q0/(Ttd2) - скорость перекачки ШФЛУ, м/с; а Ре0 - число Пекле для этого участка. В рассмотрение включаются только моменты времени t > Т1, то есть t, > 1, поскольку ранее ШФЛУ, закачиваемая в трубопровод в момент t = 0, не достигнет пункта подкачки. В пункте подкачки происходит смешение продуктов. Процентное содержание q, s-го компонента ШФЛУ в продукте после пункта подкачки (в начале следующего расчетного участка) может

быть определено из условия:

= (10)

ШФЛУ измененного в результате подкачки состава поступает в трубопровод начиная с момента времени t = Т. (t, = 1). Теперь ее скорость движения по трубопроводу равна u1 = 4(Q0 + + Q1)/^7td2), и время прохождения участка от пункта подкачки до конечного пункта Т2 = (L - L1)^1. Процентное содержание pL (t) компонентов ШФЛУ, поступающих в конечный пункт трубопровода, можно также определить по формуле (8). При этом нужно подставить вместо начального условия процентное содержание q1 (t) после пункта подкачки:

70

№ 9 сентябрь 2018 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS TRANSPORTATION AND STORAGE

уРе~ 1

Ре^.-т-!)'

*е< (т)-£/т, (11)

где Ре1 - число Пекле для участка трубопровода от пункта подкачки до конечного пункта. Теперь нас будут интересовать моменты времени t > (Т1 + Т2), т. е.

тг

Таким образом, с помощью уравнений (9)-(11) можно определить компонентный состав ШФЛУ, поступающей на конечный пункт трубопровода после путевой подкачки. Подкачек продукта может быть несколько. В этом случае расчет в соответствии с представленным алгоритмом повторяется нужное количество раз.

ПРИМЕР РАСЧЕТА

Рассмотрим реализацию алгоритма расчета изменения состава продукта при транспортировке по многопродуктовым трубопроводам на конкретном примере. Пусть ШФЛУ, кинематическая вязкость которой V = 6.10-7 м2/с, состоит из восьми углеводородных компонентов: метана, этана, пропана, изобутана, п-бутана, изопентана, п-пентана и гексана. ШФЛУ перекачивают по участку трубопровода длиной L1 = 400 км с внутренним диаметром d = 700 мм и абсолютной шероховатостью труб Д = 0,15 мм. Расход перекачки 00 = 500 м3/ч. Требуется рассчитать компонентный состав ШФЛУ, поступающей в конечный пункт участка трубопровода.

Изменение процентного содержания р0 (■) каждого компонента ШФЛУ (5 = 1, 2, ..., 8) в процессе закачки в трубопровод известно и для наглядности представлено в виде графика (рис. 1). График построен «с накоплением» и отражает равенство суммы процентных содержаний всех компонентов в любой момент времени 100 %. Алгоритм предварительного расчета следующий. Рассчитываем:

• скорость перекачки: о)0 = 400/^кё2) = = 4.500/(я.0,72.3 6 00) = 0,3 6 0 9 м/с;

• число Рейнольдса: Re0 = ^o0d/v = = 0,3609.0,7/(6.10-7) = 421 050;

100

Время t, ч Time t, h

метан

methane

этан

ethane

пропан

propane

i-бутан

i-butane

л-бутан

л-butane

i-пентан

i-pentane

л-пентан

n-pentane

гексан

hexane

Рис. 1. Компонентный состав p0(t) широкой фракции легких углеводородов в процессе закачки в трубопровод

Fig. 1. Blend composition p0(t) of wide light hydrocarbon fraction when pumped into the pipeline

• коэффициент гидравлического сопротивления (например, по формуле Альтшуля): \ = 0,11.(б8/Рв0 + Д/с/)°,25 = = 0,11.(68/421050 + 0,15/700)025 = 0,01532;

• коэффициент продольного перемешивания: К0 = 3,211^,.о)0 ё = = 3,2Ш0,01532.0,3609.0,7 = 0,1004 м2/с;

• число Пекле: Ре0 = о)0 L1/К0 = = 0,3609.400.103/0,1004 = 1,438.106;

• время движения ШФЛУ по участку трубопровода: Т1 = L1/г)0 = 400.103/ /(0,3609.3600) = 307,87 ч.

Теперь по формуле (9), заменяя интегрирование суммированием, можно определить процентное содержание рь (■) каждого из компонентов ШФЛУ, поступающих в сечение х = L1 трубопровода:

VPe «-1

Pen fa-Tj-l):

х е 4

■/>„>) rv

где N - число точек разбиения временнб-го интервала (например, N = 1 000 000), ^ = >Дтг Поскольку TN = ■) и TN = шаг Дт^ определяется для каждого момента времени как Д^ = t)/N. Расчет производится для ■ » Т1, то есть » 1. Расчетное процентное содержание р1 (■) компонентов ШФЛУ в продукте, приходящем в сечение х = L1 трубопровода, в зависимости от времени демонстрирует график рис. 2.

Как видно из сравнения рис. 1 и рис. 2, в процессе движения по трубопроводу на границах скачкообразного изменения содержания компонентов ШФЛУ образуется область плавного перехода от одного значения содержания к другому. «Ступеньки» сглаживаются аналогично тому, как это происходит на границе контакта партий при последовательной перекачке нефтепродуктов.

Пусть теперь в сечении L1 = 400 км рассмотренного ранее трубопровода осуществляется подкачка продукта другого состава. Процентное содержание к(£) компонентов подкачиваемой ШФЛУ в зависимости от времени отражено на рис. 3. Расход подкачки 01 = 300 м3/ч. Далее полученную смесь транспортируют на расстояние еще 200 км, т. е. общая длина трубопровода L составляет 600 км. Необходимо рассчитать состав рД^ ШФЛУ, поступающей в конец трубопровода.

Как следует из рис. 1, в пункте подкачки к ШФЛУ добавляются две партии, различающиеся по составу. Первая из них поступает в трубопровод до момента времени ■ = 316 ч, а вторая закачивается за ней вслед. В пункте подкачки происходит смешивание продуктов. Процентное содержание qL компонентов ШФЛУ после пункта подкачки в каждый момент времени определяется по формуле (10).

TERRITORIJA NEFTEGAS - OIL AND GAS TERRITORY No. 9 September 2018

71

ТРАНСПОРТ И ХРАНЕНИЕ НЕФТИ И ГАЗА

(t),%

Время t, ч Time t, h

метан

methane

этан

ethane

пропан

propane

7-бутан

7-butane

n-бутан

n-butane

7-пентан

7-pentane

n-пентан

n-pentane

гексан

hexane

Рис. 2. Компонентный состав p (t) широкой фракции легких углеводородов в конце L1 = 400 км участка трубопровода

Fig. 2. Blend composition p (t) of wide light hydrocarbon fraction at the end of L1 = 400 km of pipeline section

100

Время i Time t, h

метан

methane

этан

ethane

пропан

propane

/-бутан

7-butane

п-бутан

л-butane

7-пентан

7-pentane

л-пентан

n-pentane

гексан

hexane

Рис. 3. Компонентный состав kL (t) широкой фракции легких углеводородов подкачки Fig. 3. Blend composition kL (t) of wide light hydrocarbon fraction injected on the way

m /+\ OI

100

Время i Time t, h

метан

methane

этан

ethane

пропан

propane

7-бутан

7-butane

n-бутан

п-butane

7-пентан

7-pentane

п-пентан

n-pentane

гексан

hexane

Рис. 4. Компонентный состав pL(t) широкой фракции легких углеводородов в конечном пункте L = 600 км трубопровода

Fig. 4. Blend composition pL(t) of wide light hydrocarbon fraction at the terminal of L = 600 km of pipeline

Для участка трубопровода от пункта подкачки до конца участка можно рассчитать:

• скорость перекачки: ^ = 4(00 + 01)/ /(я/2) = 4.(500 + 300)/(л.0,72.3600) = = 0,5774 м/с;

• число Рейнольдса: Re1 = x>1d/v = = 0,5774.0,7/(6.10-7) = 673633;

• коэффициент гидравлического сопротивления по формуле Альтшуля:

Х1 = 0,11.(68^ + Л/с/ )°,25 = 0,11х х (68/673633 + 0,15/700)025 = 0,01466;

• коэффициент продольного перемешивания: К1 = 3,211^.1^/ = = 3,2Ш0,01466.0,5774.0,7 = 0,1571 м2/с;

• число Пекле: Ре1 = иД. - ¿1)/К1 = = 0,5 7 74.20 0.103/0,1571 = 0,735.106;

• время движения ШФЛУ по участку трубопровода: Т2 = (I - ¿1)^1 = 20 0.103/ /(0,5774.3600) = 96,22 ч.

Получаем, что ШФЛУ достигнет конечного пункта в момент t = Т1 + Т2 = 307,87 + + 96,22 = 404,09 ч (0, = 0/Т=404,09/307,87 = = 1,313).

Процентное содержание р1 (0) компонентов ШФЛУ, поступающих в конечный пункт трубопровода, определяется по формуле (11), заменяя для компьютерного расчета интегрирование суммированием:

хе"4 • ы-ч) -^(т) -Дт„,

где N - число точек разбиения временного интервала (например, N = 1 000 000), т. = 1 + >Лтг Поскольку ^ = 0, и ^ = 1 + + шаг Лтр^ определяется для каждого момента времени как Л^ = (0, - 1)/М Расчет производится для 0, » 1,313. На рис. 4 отражен результат расчета компонентного состава р. (0) ШФЛУ, приходящей в конечный пункт х = I трубопровода.

ВЫВОДЫ

Таким образом, новый универсальный алгоритм расчета смесеобразования позволяет определять концентрацию (или процентное содержание) различных компонентов нефте- и газопродуктов в процессе их последовательной перекачки по трубам. Концентрация каждого компонента в конце трубопровода

72

№ 9 сентябрь 2018 ТЕРРИТОРИЯ НЕФТЕГАЗ

OIL AND GAS TRANSPORTATION AND STORAGE

определяется как функция от времени через известную функцию от времени концентрации этого же компонента в начале трубопровода. Сопоставле-

ние графиков рис. 1-4, построенных по результатам расчета, позволяет наглядно продемонстрировать изменение компонентного состава ШФЛУ при его

движении по трубопроводу за счет конвективной и турбулентной диффузии, а также в результате подкачки ШФЛУ другого состава.

Литература:

1. Годунов Н.Н. Параметры последовательной перекачки нефтепродуктов с использованием малых противотурбулентных добавок для уменьшения объема образующейся смеси // Территория «НЕФТЕГАЗ». 2018. № 5. С. 68-72.

2. Думболов Д.У., Дроздов Д.А. Основные подходы к определению объема смеси и новый метод ее идентификации при последовательной перекачке нефтепродуктов // Территория «НЕФТЕГАЗ». 2012. № 8. С. 92-98.

3. Ишмухаметов И.Т., Исаев С. Л., Лурье M.B., Макаров С.П. Трубопроводный транспорт нефтепродуктов. М.: Нефть и газ, 1999. 300 с.

4. Лурье М.В., Марон В.И., Мацкин Л.А. и др. Оптимизация последовательной перекачки. М.: Недра, 1979. 256 с.

5. Яблонский В.С., Юфин В.А., Бударов И.П. Последовательная перекачка нефтепродуктов и нефтей по магистральным трубопроводам. М.: Гостоп-техиздат, 1959. 148 с.

6. Лурье М.В., Дидковская А.С. Изменение состава ШФЛУ при транспортировке по магистральному трубопроводу // Газовая промышленность. 2012. № 12 (683). С. 48-50.

7. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1976. 735 с.

References:

1. Golunov N.N. Parameters of the Oil Products Batching by Using Drag Reducing Agents to Decrease an Amount of Interface. Territorija "NEFTEGAS" = Oil and Gas Territory, 2018, No. 5, P. 68-72. (In Russian)

2. Dumbolov D.U., Drozdov D.A. Main Approaches to Mixture Volume Determinition and the New Method of Its Identification in Products Pipeline Operation. Territorija "NEFTEGAS" = Oil and Gas Territory, 2012, No. 8, P. 92-98. (In Russian)

3. Ishmukhametov I.T., Isayev S. L., Lurie M.V., Makarov S.P. Pipeline Bulk Oil Transportation. Moscow, Neft' i gaz [Oil and Gas], 1999, 300 p. (In Russian)

4. Lurie M.V., Maron V.I., Matskin L.A., et al. Optimization of Products Pipeline Operation. Moscow, Nedra, 1979, 256 p. (In Russian)

5. Yablonskiy V.S., Yufin V.A., Budarov I.P. Oil-Refined Products Trunk Pipeline Operation. Moscow, Gostopttekhizdat, 1959, 148 p. (In Russian)

6. Lurie M.V., Didkovskaya A.S. Change of Wide Light Hydrocarbon Fraction Composition in Trunk Pipeline Transportation. Gazovaya promyshlennost' = Gas Industry, 2012, No. 12 (683), P. 48-50. (In Russian)

7. Tikhonov A.N., Samarskiy A.A. Equations of Mathematical Physics. Moscow, Nauka [Science], 1976, 735 p. (In Russian)

19-22 НОЯБРЯ 2018

ГУБКИНСКИИ УНИВЕРСИТЕТ МОСКВА

OILANDGASHORIZONS.COM