Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

ДИПЛОМНЫЙ ПРОЕКТ

Анализ эксплуатации технологического участка Омск - Анжеро-Судженск с учетом движения партий нефти с различными физико-химическими свойствами

Введение

нефтепровод насос гидравлический

Трубопроводный транспорт является наиболее экономичным и экологически чистым видом транспорта. Этот вид транспорта обладает рядом преимуществ, такими как:

кратчайшее расстояние между пунктами;

герметичность системы

является непрерывным видом транспорта

хорошо поддается автоматизации.

Недостатками трубопроводного транспорта являются, жесткость системы и большие метало и капиталовложения.

В настоящее время создана довольно обширная и разветвленная сеть трубопроводов обеспечивающих доставку углеводородного сырья с промыслов к потребителям.

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

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

Поэтому возникает необходимость контролировать и параметры перекачиваемого продукта и своевременно корректировать режимы перекачки с учетом движения партий нефти.

На сегодняшний день, при перекачки нефти из одного пункта в другой, к трассе, подключаются /отключаются участки других месторождений. Это вызывает изменение параметров продукта по длине трубопровода. Таким образом, необходимость планировать и корректировать режимы возникает практически ежесуточно.

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

Таким же образом влияют работы по обслуживанию запорной арматуры.

Основной задачей дипломного проекта является разработка методов по расчету и корректировки режимов перекачки нефти при движении партий, с различными свойствами такими как:

- вязкости нефти;

- плотность нефти;

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

1. Краткая характеристика участка «Анжеро-Судженск-Омск»

Нефтепровод «Омск-Иркутск» условным диаметром Ду = 700 мм, совместно с нефтепроводами «Анжеро-Судженск-Красноярск-Иркутск» и «Красноярск-Иркутск» условным диаметром Ду = 1000 мм. Магистральный нефтепровод разделяется на четыре самостоятельных участка:

«Анжерская-Омск» (или «Омск-Анжерская»). Головная НПС - «Омская» ЛПДС, промежуточные: НПС «Сокур», НПС «Чулым», НПС «Барабинск» и НПС «Татарск»;

«Анжерская-Рыбное». Головная НПС - ЛПДС «Анжеро-Судженск», промежуточные НПС: «Мариинск», «Каштан», «Ачинск», «Кемчуг» и «Вознесенка»;

«Рыбинское-Нижнеудинск». Головная НПС - НПС «Рыбинское», промежуточные НПС: «Пойма», «Тайшет» и «Замзор»;

«Нижнеудинск-Ангарск». Головная НПС - ЛПДС «Нижнеудинск», промежуточные НПС: «Тулун» и «Кимильтей». Конечным пунктом является ангарский участок налива нефти.

Технологический участок «Анжеро-Судженск-Омск», имеет протяженностью 925,5 км и перепад нивелирных отметок 107 м. Основная нитка выполнена из труб диаметром 720 мм. Технологический участок имеет лупинги общей протяжённостью 603,2 км и диаметром 720 мм. Проектная пропускная способность участка составляла 12 млн. т/год, а фактическая - 10 млн. т/год.

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

Перекачка в настоящее время ведется из Омска в Анжеро-Судженск, крупным потребителем нефти в Омске является Омский нефтеперерабатывающий завод. Проведение очистки и диагностики осуществляется из Омска в Анжеро-Судженск вызвано это тем, что камеры приема-пуска СОД развернуты с учетом перекачки нефти на ВСТО.

Водные переходы нефтепровода: р. Омь - однониточный переход; р. Каргат - однониточный переход; р. Чулым - однониточный переход; р. Оеш - однониточный переход; р. Грязнуха - однониточный переход; р. Старица - двухниточный переход.

Расположение станций, относительно «Омской» ЛПДС (от нулевого километра) и высотные отметки приведены в таблице 1.8

Таблица 1.1 - Расстановка НПС по трассе нефтепровода

Название	Расположение НПС, км	Высотная отметка zi, м	Перепад высот Дzi, м
1. ЛПДС «Анжеро-Судженск»	925,5	228,4	-12,39
2. НПС «Сокур»	703,9	216,01	-79,78
3. НПС «Чулым»	525,65	136,23	-17,63
4. ЛПДС «Барабинск»	347,5	118,6	-9,7
5. НПС «Татарская»	180	108,9	12,6
6. ЛПДС «Омск»	0	121,5	-

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

1.1 Краткая характеристика НПС, участка и насосно-силового оборудования

ЛПДС «Омская» входит в состав Омского районного нефтепроводного управления (ОРНУ) акционерного общества «Транссибнефть».

Головная ЛПДС «Омск» имеет два насосных цеха. На основной насосной НПС-1 установлены 3 подпорных насоса 20НДСН и четыре магистральных агрегата НМ 2500-230 (диаметры рабочих колес 1 и 2 - 345 мм, 3 и 4 - 440 мм). На НПС -2 установлены два подпорных насоса НМП 3600 и четыре магистральных агрегата НМ 7000-210.

НПС «Татарская» имеет четыре магистральных агрегата НМ 2500-230 (диаметры рабочих колес 1 и 2 - 392 мм, 2 и 3 - 440 мм) с электродвигателями типа 4АЗМВ1-1600/6000.

НПС «Барабинск» имеет три магистральных агрегата НМ 2500-230 (диаметры рабочих колес 1 - 360 мм, 2 и 3 - 390 мм) с электродвигателями типа 4АЗМВ 2000/6000.

НПС «Чулым» имеет четыре магистральных агрегата НМ 2500-230 (диаметры рабочих колес 440 мм) с электродвигателями типа 4АЗМВ 2000/6000.

НПС «Сокур» имеет три магистральных агрегата НМ 2500-230 (диаметры рабочих колес 440 мм) с электродвигателями типа 4АЗМВ 2000/6000 и одни магистральный агрегат НМ 2500-230 (диаметр рабочего колеса 424 мм).

ЛПДС «Анжеро-Судженск» имеет два насосных цеха. На основной насосной НПС-1 установлены четыре насоса НМ 2500-230 (диаметры рабочих колес 440 мм) с электродвигателями типа 2АЗМВ1-2000/6000, мощность которых 2000 кВт и число оборотов вала 2979 об/мин. Для насосов НПС-1 подпор создается подпорной насосной, в которой установлены три подпорных насоса НМП 3600-78 с электродвигателями типа ВАО 2-560L номинальной мощностью 5000 кВт.

На основной насосной НПС-2 установлены четыре насоса НМ 10000-210 с электродвигателями типа СТД 6300, мощность которых 63000 кВт, а число оборотов вала - 3000 об/мин; подпор создается одним из четырех подпорных насосов НПВ 5000-120 с электродвигателями типа ВАОВ-800L-4.

1.2 Основные насосно-силовые агрегаты и их характеристики

Характеристики насосных агрегатов представлены в таблицах 1.1-1.6

Таблица 1.1 - Характеристика подпорного насоса 20НДСН

Q, м3/ч	0	584	1168	1460	1753	2045	2337	2629	2921
Н, м	51,2	51	47,4	45,5	44,6	43,6	42,7	39,8	36

Таблица 1.2 - Характеристика насоса НМ - 2500-230 (345)

Q, м3/ч	0	375	750	1125	1312.5	1500	1687,5	2250	2700	3700
Н, м	166	165,3	159,9	150,4	146	138	130	104	80	0

Таблица 1.3 - Характеристика насоса НМ - 2500-230 (440)

Q, м3/ч	0	500	1000	1500	2000	2500	3000	3500	4000	4500
Н, м	280	276	270	255	240	227	210	178	110	0

Таблица 1.4 - Характеристика насоса НМ - 2500-230 (424)

Q, м3/ч	0	500	1000	1500	2000	2500	3000	3500	4000	4500
Н, м	259	258	249	239	220	194	164	125	70	0

Таблица 1.5 - Характеристика насоса НМ - 2500-230 (390)

Q, м3/ч	0	443	886	1330	1551	1773	1994	2600	3500	4500
Н, м	210	207	202	196	192	186	179	159	97	0

Таблица 1.6 - Характеристика насоса НМ - 2500-230 (360)

Q, м3/ч	0	409	818	1227	1432	1636	1841	2045	2250	2455
Н, м	187	185	181	175	172	167	161	154	147	137

Определение аппроксимационных коэффициентов для насосов

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

Метод наименьших квадратов служит для оценки неизвестных величин по результатам измерений (эксперимента), содержащим случайные погрешности. Среди многих приложений метода наиболее важным является нахождение наилучшего уравнения (функциональной зависимости) определенного вида для представления опытных данных.

Процесс выражения опытных данных функциональной зависимостью с помощью метода наименьших квадратов состоит из двух этапов: на первом выбирают вид искомой формулы, на втором для данной формулы подбираются параметры [1].

Построим табличную зависимость для подпорного насоса 20НДСН, (см. рисунок 1.1)

Рисунок 1.1 - Паспортная характеристика насоса 20НДСН

В соответствии с идеей метода наименьших квадратов необходимо минимизировать сумму:

 (1.1)

где - значение функции, определяемое опытным путем в точке ;

- теоретическая зависимость

В случае зависимости функция является функцией двух переменных и сумма принимает вид:

(1.2)

Проиллюстрируем метод для нахождения коэффициентов в аналитической зависимости вида (1.2), для N экспериментальных значений тогда невязка будет равна

(1.3)

Составим систему уравнений

(1.4)

после преобразований получим:

(1.5)

(1.6)

где , число опытных величин.

Прежде чем воспользоваться формулами (1.5) и (1.6) найдем значения сумм

, ,

,

Определим коэффициенты по формулам (1.5) и (1.6)

a₀== 50,76

a₂==1,76*10^-6

Построим графики функций и паспортной характеристики (см. рисунок. 1.2.)

Аналогичным образом найдены аппроксимационные коэффициенты для остальных насосов и приведены в таблице 1.7.

Таблица 1.7 - Данные по насосам

Марка насоса	Диаметр рабочего колеса, мм	a	b
НМ - 2500-230	345	166,12	0,0000121
НМ - 2500-230	440	290,85	0,0000121
НМ - 2500-230	424	266,29	0,0000124
НМ - 2500-230	390	215,26	0,0000102
НМ - 2500-230	360	187,26	0,0000803

Рисунок 1.2 - Совмещение паспортной и аппроксимированной зависимости

Рисунок 1.3 - Совмещение паспортной и аппроксимированной зависимости насоса НМ - 2500-230 (345)



Рисунок 1.4 - Совмещение паспортной и аппроксимированной зависимости насоса НМ - 2500-230 (440)

Рисунок 1.5 - Совмещение паспортной и аппроксимированной зависимости насоса НМ - 2500-230 (424)



Рисунок 1.6 - Совмещение паспортной и аппроксимированной зависимости насоса НМ - 2500-230 (390)

Рисунок 1.7 - Совмещение паспортной и аппроксимированной зависимости насоса НМ - 2500-230 (360)

Как видно из рисунка 1.2 вычисленные аппроксимационные коэффициенты подобранны правильно - подобранная функция практические полностью совпадает с таблично, следовательно, погрешности будут минимальными и не будут превышать 0,1%.

Из рисунка 1.3 следует, что подобранная функция подходит табличной зависимости, следовательно, погрешности будут минимальными и не будут превышать 0,1%.

Из рисунка 1.4 следует, что подобранная функция подходит табличной зависимости, но при определении напора на участках подач от 0 до 2100 м³/ч и 2500 до 4000 м³/ч будет погрешность не превышающая 2%.

Из рисунка 1.5 следует, что подобранная функция практически полностью совпадает с табличной зависимостью, погрешность определения напора составит не более 1%.

Из рисунка 1.6 следует, что подобранная функция практически полностью совпадает с табличной зависимостью, погрешность определения напора составит не более 1%

Из рисунка 1.7 следует, что подобранная функция практически полностью совпадает с табличной зависимостью, погрешность определения напора составит не более 1%

1.4 Физико-химические свойства нефтей

Параметры перекачиваемой нефти взяты в лаборатории химического анализа Омского РНУ ОАО «Транссибнефть» усредненные (всесезонные)

– плотность нефти = 0,853 тн/м³ = 853 кг/м³;

– вязкость нефти = 8,5 сСт = 8,510^-6 м²/с;

– давление насыщенных паров P_s = 480 мм рт. ст. = 63984 Па.

Параметры нефти при стандартных условиях приводятся в таблице 1.9

Таблица 1.9 - Свойства перекачиваемой нефти

Показатель	Величина
Плотность при 20°С, кг/м3 Вязкость при 20°С, сСт Вязкость при 50°С, сСт Содержание хлористых солей, % Содержание воды, % Суммарное содержание балласта, % Содержание серы, % Давление насыщенных паров при температуре 37,8°С, кПа	826,5 9,8 2,6 0,0013 0,148 0,1428 0,62 62,1

За расчетную температуру перекачки принимается минимальная температура окружающей магистральный нефтепровод среды с поправкой на начальную температуру нефти, тепловыделение в трубопроводе, обусловленное трением потока, и теплопередачу тепла в окружающую среду. Для заглубенных трубопроводов это минимальная температура грунта на глубине заложения трубопровода.

Таблица 1.10 - Партии нефти, перекачиваемые по трубопроводу

Номер партии	Вязкость при 20°С, сСт	Вязкость при 50°С, сСт	Плотность при 20°С, кг/м3
1	8	2,5	842
2	9	4	793
3	8,4	3,3	850
4	9,8	2,9	822

Плотность и вязкость нефти зависят от температуры: при повышении температуры они уменьшаются, а при понижении - увеличиваются. Плотность при температуре перекачки t в ⁰С определяется по формуле [2]

(1.7)

где с₂₀ - плотность нефти при t = 20⁰C,

о - коэффициент объемного расширения.

Вязкость нефтей различных месторождений значительно отличается и в справочной литературе приводятся, как правило, значения кинематической вязкости нефти определенного состава при температуре 20 и 50 ⁰С н₂₀ и н₅₀.

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

Зависимость плотности нефти от давления представляется формулой

(1.8)

в которой - плотность нефти при стандартных условиях (атмосферном давлении P₀ и температуре +20 С); - коэффициент сжимаемости нефти, среднее значение которого составляет 0,00078 МПа^-1

В тех случаях когда одновременно отклоняются от номинальных значений и температура и давление, плотность нефти при давлении и температуре можно рассчитывать по формуле

(1.9)

где K - модуль упругости нефти и равен 1,3*10⁹ Па

Давление насыщенных паров определяется экспериментально для каждой нефти и также зависит от температуры перекачки, однако для НПС МН можно принять Р_S = 500 мм рт. ст. = 66650 Па.

Расчетную вязкость можно вычислить, например, по формуле Рейнольдса - Филонова:

(1.10)

где н₀ - вязкость при температуре t₀ = 20 или 50 ⁰С, а k - эмпирический коэффициент, равный

(1.11)

Расчетную кинематическую вязкость нефти определяем по формуле Вальтера

(1.12)

где _Т - кинематическая вязкости нефти, мм²/с;

А_н и В_н - постоянные коэффициенты, определяемые по двум значениям вязкости ₁ и ₂ при двух температурах Т₁ и Т₂:

; (1.13)

где , - расчетные значения вязкости при температурах и соответственно

; (1.14)

Произведем пересчет партии нефтей на рабочую температуру. Для этого воспользуемся зависимостями (1.4)

Рабочая температура t=5 ^oC или 278 К

Определим плотность нефти при расчетной температуре для первой партии по формуле (1.7)

кг/м³

Расчетную кинематическую вязкость нефти определяем по формуле Вальтера (1.12). Для этого определим постоянные величины при помощи зависимостей (1.13) и (1.14)

где сСт, сСт, Т₂=323К, Т₁=293К

;

сСт

Определим плотность нефти при расчетной температуре для второй партии по формуле (1.7)

кг/м³

Расчетную кинематическую вязкость нефти определяем по формуле Вальтера (1.12). Для этого определим постоянные величины при помощи зависимостей (1.13) и (1.14)



где сСт, сСт, Т₂=323К, Т₁=293К

;

сСт

Определим плотность нефти при расчетной температуре для третьей партии по формуле (1.7)

кг/м³

Расчетную кинематическую вязкость нефти определяем по формуле Вальтера (1.12). Для этого определим постоянные величины при помощи зависимостей (1.13) и (1.14)

где сСт, сСт, Т₂=323К, Т₁=293К

;

сСт

Определим плотность нефти при расчетной температуре для четвертой партии по формуле (1.7)

кг/м³

Расчетную кинематическую вязкость нефти определяем по формуле Вальтера (1.12). Для этого определим постоянные величины при помощи зависимостей (1.13) и (1.14)

где сСт, сСт, Т₂=323К, Т₁=293К

;

сСт

1.5 Определение влияния вязкости на изменение характеристик насоса

Известно, что свойства перекачиваемой жидкости могут влиять на напорную характеристику насосов, поэтому необходимо определить будут ли насосы изменять своих характеристики в исследуемом диапазоне вязкостей от 15 до 25 сСт.

Определяем границы автомодельного режима на примере магистрального насоса НМ2500-230 с рабочим колесом D₂=440 мм.

Условную границу перехода режима течения жидкости из автомодельной области в область зависящих от вязкости значений параметров насоса определяют переходное и граничное число Рейнольдса вычисляемые по формулам:



где - коэффициент быстроходности насоса

Зная число можно найти предельное значение вязкости, начиная с которой необходимо вести пересчет характеристик насоса:

сСт

где -число оборотов ротора насоса, -диаметр рабочего колеса насоса.

Следовательно, при вязкость продукта не будет влиять на напорную характеристику, так как диапазон вязкостей намного меньше следовательно перекачиваемый продукт не будет влиять на характеристику магистрального насоса НМ2500-230 (440).

Аналогично проведем расчет для всех насосов (см. таблицу 1.11)

Таблица 1.11 - Таблица переходных вязкостей

Тип насоса	Диапазон сСт	Переходная вязкость, сСт
НМ - 2500-230 (440)	15-25	128,11
НМ - 2500-230 (424)		118,96
НМ - 2500-230 (390)	15-25	100,3
НМ - 2500-230 (360)		85,2
НМ - 2500-230 (345)		78,76

Как видно из таблицы 2.5 - в исследуемом диапазоне вязкость продукта не влияет на напорную характеристику насосов.

1.6 Факторы, влияющие на распределение напоров

При расчете технологический режимов по всесезонным усредненным параметрам расчетная модель может не совпадать с фактическим распределением напоров по длине трассы.

Неравномерный характер обусловлен следующими группами факторов:

1) Топографические и климатические условия: неверная интерпретация показаний датчиков давления, обусловленная ошибочным определением нивелирных отметок и места расположения по трассе; сезонные изменения температуры и т.д.

2) Технические средства измерений: отклонение погрешности датчиков давления и т.п.

3) Конструктивные особенности МН: искажение диаметра в месте, где датчик давления показывает аномальное давление и т.п.

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

2. Расчетная часть

2.1 Сущность последовательной перекачки нефтей методом прямого контактирования

В настоящее время нефти различных промыслов и месторождений, а число таких нефтей, отличающихся по физико-химическим свойствам (плотности, вязкости, содержанию серы, минеральных солей, парафина и т.п.), может достигать несколько десятков, перекачивают методом, получившим название «последовательная перекачка прямым контактированием» [2].

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

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

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

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

На рисунке 2.1 приведена принципиальная схема трубопроводной системы для последовательной перекачки нефтей, где видны последовательно движущиеся в магистральном нефтепроводе (МНП) партии нефтей (№1, 2, 3, 4, 5 и т.д.), вытесняющие предыдущие партии и в свою очередь вытесняемые последующими, а также нефтепромыслы (П1) и (П2) с пунктами подготовки нефти к транспорту (ППНТ) резервуарные парки (РП) головной перекачивающей станции (ГПС), резервуарные парки промежуточных перекачивающих станций (ППС) и резервуарный парк конечного пункта (КП), на котором происходит прием нефтей в резервуары нефтеперерерабатывающего завода (НПЗ) или перевалочной нефтебазы для перевалки на другие виды транспорта.

В систему для последовательной перекачки нефтей входят те же основные объекты, что и в систему транспортировки однородной нефти:

- соединительные трубопроводы (СТ), соединяющие пункты подготовки нефти к транспорту с головной перекачивающей станцией;

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

- линейная часть магистрального нефтепровода с линейными сооружениями, обеспечивающими ее функционирование;

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

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

Рисунок 2.1 - Система перекачки

2.2 Гидравлические особенности расчета нефтепровода при последовательной перекачки нефтей

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

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

Если на участке нефтепровода (рисунок 2.2) происходит вытеснение одной нефти другой, а плотности и вязкости этих жидкостей отличаются друг от друга, то уравнение Бернулли в обычном виде неприменимо. Предполагая процесс замещения нефтей практически стационарным (или, как говорят, квазистационарным) и пренебрегая протяженностью зоны смеси (х_с), запишем уравнения Бернулли для участков (х₁ х_с) и (х_с, х₂) трубопровода, занятых соответственно нефтью №1 и №2 [2]:

Рисунок 2.2 - Вытеснение одной партии нефти другой нефтью

 (2.1)

Здесь p_lr p₂ - плотности нефтей; h_1-c; h_c-2 - потери напора на участках, занятых нефтью №1 и №2, соответственно; z_c - высотная отметка профиля трубопровода в месте контакта партий.

Исключив из этих уравнений давление р_с в месте контакта партий, получим обобщенное уравнение Бернулли для рассматриваемого случая:

(2.2)

Если же исключить из уравнений (2.1) высотную отметку z_c, получим еще одно уравнение:

(2.3)

В пренебрежении потерями напора на местных сопротивлениях величины h_1-c и h_c-2 обозначают здесь потери напора на участках (x₁, х_с) и (х_с, х₂) соответственно:

(2.4)

где L_i - протяженности этих участков.

Уравнение (2.2) служит для определения скорости и перекачки, а уравнение (2.3) - для определения давления р_с в месте контакта нефтей.

Уравнения (2.2) и (2.3) полезно записать в терминах давлений при этом после преобразований получим:

(2.5)

Из второго уравнения следует, что разность напоров (Н₁ - H₂) между началом и концом трубопровода не равна сумме гидравлических потерь на участках, занятых первой и второй нефтью; она содержит еще одно слагаемое, обусловленное разностью плотностей перекачиваемых нефтей (см. рисунок 2.2). Это означает, что напор Н(х), представляемый линией АСДВ гидравлического уклона, не является непрерывным; в месте контакта нефтей существует разрыв СД. Значение этого разрыва

(2.6)

называется скачком напора.

Скачок напора =0, если, т.е. плотности нефтей равны; скачок < 0 (напор скачком уменьшается), если р₂ > р_1г т.е. менее плотная нефть вытесняет более плотную нефть; скачок > 0 (напор скачком увеличивается), если р₂ < р_1г т.е. более плотная нефть вытесняет менее плотную.

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

2.3 Методика гидравлического расчета при последовательной перекачки нефтей

Исходя, из проведенных исследований влияния входных параметров на режим перекачки стало - известно, что все входные параметры необходимо учитывать при расчете.

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

Для того чтобы все необходимые параметры были известны для гидравлического расчета при последовательной перекачки, необходимо провести гидравлический расчет технологического участка по усредненным всесезонным параметрам.

В первую очередь необходимо провести перерасчет свойств нефти при стандартных условиях на температуру перекачки нефти по формулам (1.7), (1.12)

Затем зная схему включении насосов определяем расход методом итерации:

Принимая в первом приближении к-т гидравлического сопротивления запишем уравнение для определения расхода

, (2.7)

где ;- сумма всех коэффициентов магистральных насосов, работающих на станциях,

и - аппроксимационные коэффициенты подпорного насоса,

Нк - конечный напор, равный уровню взлива в резервуаре на последней станции технологического участка,

L - длина технологического участка,

D - диаметр эквивалентный.

После определения расхода в первом приближении ведется пересчет

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

(2.8)

(2.9)

после определения производим проверку сравнивая с

(2.10)

Если расхождение составляет менее 1%, то расчет прекращается и по последнему коэффициенту ведется расчет расхода, если погрешность превышает 1%, то расчет повторяется вновь, пока расхождение не будет менее 1%.

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

- для определения подпора на головной станции

; (2.11)

- для определения напора создаваемого магистральными насосами первой станции

; (2.12)

- для определения напора на первой станции

; (2.13)

Для определения подпора перед следующей станцией воспользуемся зависимостью:

, (2.14)

где - потери на трение по длине участка

(2.15)

Зависимость (2.15) не совсем удобна для расчетов т.к. здесь расход выражен через скорость, поэтому преобразуем это уравнение в другой вид с помощью формулы:

, где S - площадь трубопровода равная , тогда

уравнение (2.15) примет вид:

, тогда (2.16)

с учетом уравнения (2.14) получим:

(2.17)

Аналогично для всех остальных участков находятся напоры после станций и подпоры перед ними.

По приведенным выше напорам и подпорам можно сделать вывод о работоспособности режима, а именно про анализировать нет ли превышений по уставкам давлений.

После расчета по усредненным всесезонным параметрам ведется расчет с учетом движения партий нефти с различными физико-химическими свойствами.

Перед тем как вести расчет стало возможным определить зависимость плотности от давления, т.е. для более точного расчета необходимо знать не только как плотность меняется от температуры, но и как плотность меняется от давления. По обобщенной формуле (1.8) найдем эту зависимость для каждой партии нефти.

После пересчета партий нефти ведется основной гидравлический расчет по уравнению Бернулли

, (2.18)

где - потери на по длине трубопровода определяются по формуле Лейбензона:

, (2.19)

где - коэффициенты при определенном режиме течения.

В таблице 2.1 приведены значения коэффициентов и m для различных режимов течения жидкости. Согласно [3]

Таблица 2.1 - Зависимость коэффициентов и m от режимов течения

Режим течения	Re	m
Ламинарный	Re<2320	1	4,153
Турбулентный	Гидравлически гладкие трубы	2320<Re<	0,25	0,0246
	Смешанное трение	<Re<	0,1	0,0166
	Квадратичное трение	Re>	0	9,08910-30,25

Расчет ведется из конца трубопровода к началу с интервалом в 1 км, предполагая, что все насосы расположены в начале технологического участка. Это обусловлено тем, что расчет ведется с учетом различных нефтей с различными свойствами.

Преобразуя уравнение (2.18) получим:

(2.20)

За расчетный расход принимается расход, найденный при расчете по усредненным параметрам в качестве первого приближения. Причем расчет ведется из резервуара в конце участка в резервуар в начале участка.

Зная распределение партий нефти по длине трубопровода по уравнению (2.20) находим давления в промежуточных точках по всей длине трассы. За расчетный шаг принимается 1 км.

Для более удобного анализа необходимо перейти от единиц давления в единицы напора. Для этого воспользуемся зависимостью:

 (2.21)

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

Соответственно получиться график распределения напоров с учетом движения партий нефти с различными физико химическими свойствами. Анализируя который стало возможным наблюдать реальные распределения напоров.

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

2.4 Гидравлический расчет по усредненным всесезонным параметрам

Для разработки и демонстрации методики выбран участок «Омск-Анжеро-Судженск», исходными данными для расчета являются:

- длина участка - 925,5 км;

- высотная отметка начального участка (Омска) - 120 м;

- высотная отметка конечного участка (Анжеро-Судженск) - 227 м;

- перепад высот.- 107 м;

- вязкость усредненная - 9,5 сСт;

- плотность усредненная - 826 кг/м3;

- насосное оборудование:

2 подпорных насоса - один рабочий, второй резервный,

4 магистральных насоса на каждой станции - один резервный, остальные работают в зависимости от схемы перекачки.

Насосы во всех схемах перекачки соединяются последовательно.

Воспользуемся уравнением (2.7)

= 1565,17; = 51,3 + 30; Нк = 30 м; = 107 м; = 3.526*10^-4

=0,0000133; D=0,682 м; = 0,02 (в первом приближении), тогда

Производим проверку:

- Определяем число Re



- Определяем коэффициент гидравлического сопротивления, предполагая, что режим течения находится в зоне Блазиуса (гладкой трубы)

- После определения производим проверку, сравнивая с , используя зависимость (2.10)

Из заданного диапазона в 1% коэффициент не выходит, следовательно его можно принять для определения расхода. Таким образом, расход найден.

Необходимо провести расчет на подтверждение предполагаемого режима течения. Режимы течения и зависимости их описывающие представлены в таблице 2.2

Таблица 2.2 - Режимы течений жидкости

Режим течения	Re
Ламинарный	Re<2320
Турбулентный	Гидравлически гладкие трубы	2320<Re<
	Смешанное трение	<Re<
	Квадратичное трение	Re>

Для этого необходимо определить переходное число Re

Переходное число Re определиться из следующей зависимости

В свою очередь к-т равен

,

где - относительная шероховатость трубы равная 2 мм.

Откуда переходное число Re определиться

Как видно, что число Re намного меньше найденного переходного числа, следовательно режим течения определен верно.

Далее определяем подпор на первой станции используя зависимость (2.11)

Затем определяем напор создаваемый магистральными насосами первой станции по формуле (2.12)

Определяем напор на первой станции ГНПС «Омск» при помощи формулы (2.13)

Далее определяем напоры, создаваемые магистральными насосами на последующих станциях:

Напор на станции НПС «Татарск»

Напор на станции ГНПС «Барабинск» работают два насоса последовательно.

Напор на станции НПС «Чулым» работает один насос.

Напор на станции НПС «Сокур» работает один насос.

Определяем напоры на входе станций, которые являются подпорами, необходимыми для без кавитационной работы агрегатов, установленных на станциях. Подпоры определяются по формуле (2.17)

Определяем напор на входе в НПС «Татарск»

Определяем напор на входе в ГНПС «Барабинск»

Определяем напор на входе в НПС «Чулым»

Определяем напор на входе в НПС «Сокур»

Определяем напор на входе в ГНПС «Анжеро-Судженск»

Построим распределение напоров в виде графика (см. рисунок 2.4)



2.5 Определение влияния плотности от давления

После гидравлического расчета по усредненным параметрам стало возможным определить влияние давления на реологические свойства нефтей. В частности на плотность, т.к. при давлении менее 10 МПа вязкость практически не изменяется. Про анализируем влияние давления на плотность. За давление будем принимать то давление, которое определилось при расчетах по усредненным параметрам. Для более точного расчета будем использовать базовые точки по всей длине технологического участка с определенным интервалом. В качестве базовых точек возьмем точки по длине трассы с соответствующими давлениями, представленные в таблице 2.5

Таблица 2.5 - Базовые точки

P, МПа	L, км	P, МПа	L, км	P, МПа	L, км	P, МПа	L, км	P, МПа	L, км
3,13	0	2,94	180,3	4,46	347,5	4,29	525,65	4,12	703,9
2,67	50	2,75	200	3,98	400	4,07	550	3,69	750
2,20	100	2,29	250	3,51	450	3,61	600	3,23	800
1,74	150	1,83	300	3,05	500	3,14	650	2,77	850
1,46	180,3	1,39	347,5	2,81	525,65	2,65	703,9	2,07	925,5

Воспользуемся зависимостью (1.8) и находим плотность на первом участке в 0 км. [4]

Плотность на отметке 50 км.

Плотность на отметке 100 км.

Плотность на отметке 150 км.

Плотность на отметке 180,3 км.

Далее пересчитанные значения плотностей на других участках приведены в таблице 2.6

Таблица 2.6 - Распределение плотностей

L, км	p, кг/м3	L, км	p, кг/м3	L, км	p, кг/м3	L, км	p, кг/м3
180,3	838,73	347,5	839,70	525,65	839,59	703,9	839,48
200	838,61	400	839,39	550	839,45	750	839,21
250	838,32	450	839,09	600	839,158	800	838,92
300	838,02	500	838,80	650	838,86	850	838,62
347,5	837,75	525,65	838,6	703,9	838,54	925,5	838,18

После расчета плотностей в базовых точках необходимо пересчитать напор по формуле , т.к. напор будет отличаться в связи изменением плотности по длине трассы. Напоры представлены в таблице 2.7

Таблица 2.7 - Распределение напоров

L, км	H, м	L, км	H, м	L, км	H, м	L, км	H, м	L, км	H, м
0	380,8	180,3	357,4	347,5	542,1	525,65	521,9	703,9	500,8
50	324,7	200	335,3	400	483,3	550	494,6	750	449,2
100	268,6	250	279,2	450	427,3	600	438,7	800	393,2
150	212,5	300	223,1	500	371,3	650	382,6	850	337,2
180.3	178,4	347,5	169,7	525,65	342,6	703,9	322,2	925,5	252,5

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



Из рисунка 2.5 видно, что линии распределения напоров практически совпадают. Это позволяет сделать вывод, что зависимость плотности нефти от давления незначительно и в дальнейших расчетах не стоит это учитывать, т.к. расчет очень неудобный.

2.6 Гидравлический расчет с учетом движения партий нефти с различными физико-химическими свойствами

Исходными данными для расчета являются:

- длина участка - 925,5 км;

- высотная отметка начального участка (Омска) - 120 м;

- высотная отметка конечного участка (Анжеро-Судженск) - 227 м;

- перепад высот.- 107 м;

- насосное оборудование:

2 подпорных насоса - один рабочий, второй резервный,

4 магистральных насоса на каждой станции - один резервный, остальные работают в зависимости от схемы перекачки.

Партии нефти представлены в таблице 2.5

Таблица 2.5 - Значение свойств нефтей

Номер партии	Вязкость при 20°С, сСт	Вязкость при 50°С, сСт	Плотность при 20°С, кг/м3
1	8	2,5	842
2	9	4	793
3	8,4	3,3	850
4	9,8	2,9	822

Так как расчетная усредненная температура перекачки равна 5^о С,

то дальше в таблице 2.6 приведены характеристики нефти при данной температуре.

Таблица 2.6 - Свойства нефтей

Номер партии	Вязкость при 5°С, сСт	Плотность при 5°С, кг/м3
1	19,3	853
2	15,54	804
3	16,19	860
4	24,87	832

Распределение партий нефти по длине участка трубопровода приведена на рисунке 2.6



Рисунок 2.6 - Распределение партий нефти

Расчет ведется по уравнению Бернулли из конца технологического участка, где известен напор равный уровню взлива резервуара. За расчетный шаг принимаем 1 км.

Рисунок 2.7

Запишем уравнение баланса напоров для всего участка: [5]

, (2.22)

где - гидроуклон отдельно взятой партии нефти со своими свойствами.

Т.к. численно решить такое уравнение не представляется возможным, то разделим это уравнение на две части и будем решать их отдельно, т.е. отдельно решаем левую и правую части так, чтобы они оказались равными или их разность была равна нулю.

Для этого запишем по отдельности это уравнение, заменяя напоры на давления.

Для решения правой части уравнения, вместо левой запишем давления начальные и в качестве первого приближения, подставляем расход, найденный по усредненным параметрам.

Запишем уравнение:

(2.23)

Используя выведенную зависимость (2.23) будем рассчитывать линию распределения давлений, после чего при помощи уравнения (2.21) преобразовывать в распределение напоров.

Определим давление на границе партий нефти с конца участка 925,5 км до границы 831 км:

Определим давление на границе партий нефти с 830 км до границы 800 км:

Определим давление на границе партий нефти с 799 км до границы 730 км:

Определим давление на границе партий нефти с 729 км до границы 608 км:

Определим давление на границе партий нефти с 607 км до границы 497 км:

Определим давление на границе партий нефти с 496 км до границы 342 км:

Определим давление на границе партий нефти с 341 км до границы 194 км:

Определим давление на границе партий нефти с 193 км до границы 0 км:



Через уравнение переведем распределение из единиц давления в единицы напора.

Определим напор на 831 км:

Определим напор на 830 км:

Определим напор на 800 км:

Определим напор на 799 км:

Определим напор на 730 км:

Определим напор на 729 км:

Определим напор на 608 км:

Определим напор на 607 км:

Определим напор на 497 км:

Определим напор на 496 км:

Определим напор на 342 км:

Определим напор на 341 км:

Определим напор на 194 км:

Определим напор на 193 км:

Определим напор на 0 км:

Определим напор на резервуаре:

м

Далее построим графики распределения давлений и распределения напоров и проведем анализ линий гидроуклонов.



Рисунок 2.8 - Распределение давлений

Покажем на рисунке 2.9 распределение напоров по всему технологическому участку.

Рисунок 2.9 - Распределение напоров

Из системы координат «Расстояние - давление» в зоне контакта партий нефти скачки не наблюдаются, а заметно лишь искривление. Это связанно с изменение вязкости по длине трассы.

На рисунке 2.9 наоборот наблюдаются скачки напоров в зоне контакта различных партий и их отчетливо видно. Это связанно с различной плотностью у разных партий нефти.

На графике распределения напоров отчетливо видно, что гидроуклоны, посчитанные по усредненным параметрам и с учетом различных партий не попадают в одну точку в начале технологического участка. Это связанно с неточным определением пропускной способности т.к. при расчетах гидроуклона с различными партиями в качестве первого приближения использовался расход, посчитанный по усредненным параметрам. Расчет показал, что он не подходит. Для того что бы определить расход, нужно свести линии гидроуклона в одну точку в начале трассы. Т.к. в ручную это очень долгий и тяжелый подсчет, то расчет велся в пакете программ Excel, при помощи функции «подбор параметра». Таким образом после расчета в Excel, был найден расход, который удовлетворял всем условиям.

Q = 0,302 м³/с

Построим график распределения напоров (см. Рисунок 2.10)

2.7 Гидравлический расчет с учетом движения партий нефти с различными физико-химическими свойствами с подкачкой на НПС «Барабинск»

Исходными данными для расчета являются:

- длина участка - 925,5 км;

- высотная отметка начального участка (Омска) - 120 м;

- высотная отметка конечного участка (Анжеро-Судженск) - 227 м;

- перепад высот.- 107 м;

- насосное оборудование:

2 подпорных насоса - один рабочий, второй резервный,

4 магистральных насоса на каждой станции - один резервный, остальные работают в зависимости от схемы перекачки.

Партии нефти представлены в таблице 2.9

Таблица 2.9 - Значение свойств нефтей

Номер партии	Вязкость при 20°С, сСт	Вязкость при 50°С, сСт	Плотность при 20°С, кг/м3
1	8	2,5	842
2	9	4	793
3	5,3	4	832
4	9,8	2,9	822

Так как расчетная усредненная температура перекачки равна 5^о С, то дальше в таблице 2.10 приведены характеристики нефти при данной температуре.

Таблица 2.10 - Значение свойств нефтей

Номер партии	Вязкость при 5°С, сСт	Плотность при 5°С, кг/м3
1	19,3	853
2	15,54	804
3	10,1	845
4	24,87	832

Распределение партий нефти по длине участка трубопровода приведена на рисунке 2.11

Рисунок 2.11 - Распределение партий нефти

Для гидравлического расчета технологического участка МН с учетом движения партий нефти необходимо привести методику гидравлического расчета с подкачкой по осредненным параметрам, затем считать с различными партиями нефтей.

Рисунок 2.12 - Технологический режим с подкачкой

Запишем систему уравнений для всего технологического участка МН [5]

 (2.24)

Из системы (2.24) не определены 2 величины: Q и .

Значение q задается таким образом, чтобы обеспечить посуточный график движения нефти, поэтому это известная величина. Складываем два уравнения, тогда получим:

, (2.25)

где , а , тогда численно числено решить не представляется возможным. С помощью пакета программ Excel и функции «подбор параметра» найти Q возможно, если от правой части уравнения отнять левую и приравнять нулю.

(2.26)

После определения Q находим из любого уравнения системы (2.24) величину

(2.27)

Затем используя зависимости (2.12), (2.13), (2.17) определяем подпор на первой станции, напор, создаваемый основными насосными агрегатами, и подпоры перед следующими станциями соответственно.

Запишем уравнение (2.26) численно, если q = 3 тыс. тонн/сут. Или

q = 0,042 м³/с

Решая это уравнение через «подбор параметра» получили . При данном значении левая и правая часть равны

Затем при помощи уравнения (2.27) определим подпор перед станцией, на которой производится подкачка:

Далее используя зависимости (2.12), (2.13), (2.17), находим подпор на первой станции, напор, создаваемый основными насосными агрегатами, и подпоры перед следующими станциями соответственно.

Определяем подпор на первой станции используя зависимость (2.11), если Q - q = 0,346 - 0,042 = 0,304 м³/с

Затем определяем напор создаваемый магистральными насосами первой станции по формуле (2.12)

Определяем напор на первой станции НПС «Омск» при помощи формулы (2.13)

Далее определяем напоры, создаваемые магистральными насосами на последующих станциях

Напор на станции НПС «Татарск»

Напор на станции ГНПС «Барабинск» работают один насос.

Напор на станции НПС «Чулым» работает один насос.

Напор на станции НПС «Сокур» работает один насос.

Определяем напоры на входе станций, которые являются подпорами, необходимыми для без кавитационной работы агрегатов, установленных на станциях.

Определяем напор на входе в НПС «Татарск»

Определяем напор на входе в ГНПС «Барабинск»

Определяем напор на входе в НПС «Чулым»

Определяем напор на входе в НПС «Сокур»