Задачник по трубопроводному транспорту нефти нефтепродуктов и газа
..pdf211
р* - ра в сечении утечки рассчитывается так, как если бы ее
не было. Имеем: |
|
|
|
Н =150+ 4,5-106 |
= 677,26 м, Н |
=100+ |
0,3-106 = 135,15 м, |
870*9,81 |
' |
ж |
870*9,81 |
i = i k ^ k = 6 7 7 '2 6 -1 3 5 ,1 5 s |
, ш -з |
L120000
Н. = Н„ - i - L . =677,26 - 4,5176-10‘3 -80000 = 315,85 м,
AH = H . - z , = 315,85 -50 = 265,85 м,
q=H • s • V2g-AH = 0,62 ■0,000Ц/2 • 9,81 • 265,85 = 4,478 • Ю'3 м3/с.
Объем V вытекшей за 6 часов нефти составляет:
4,478 • 10‘3 • 6 • 3600 = 96,7 м3.
85. В данной задаче в отличие от предыдущей разност
АН напоров, входящую в формулу q = |Х • s • >/2g*AH , нельзя вычислять предполагая, что появившаяся утечка не изменя
ет режима перекачки, |
ведь площадь |
отверстия в стенках |
|
трубопровода ( s = 25 |
Л |
велика |
и истечение |
см ) достаточно |
|||
жидкости через отверстие делает скорости v, |
течения неф |
ти до сечения утечки и v2 - после нее не равными друг дру гу. Для решения задачи необходимо составить и решить полную систему уравнений, получающуюся на основе урав нения Бернулли, примененного к двум сегментам участка - до и после сечения утечки:
f |
4 5 • 106 |
^ |
|
80000 |
V, |
|
150+ —— — — |
-(5 0 + А Н ) = Х, |
0,7 |
2*9,81 |
|
|
870*9,81 |
|
|
||
|
( |
|
0,3-Ю6 Л_ ^ |
40000 |
2 |
(5 0 + А Н )- 100+ |
v2 |
||||
|
V |
|
870*9,81 |
0,7 |
2*9,81 |
|
|
|
|
|
|
(v ,- v2) •3,14 0-’ - |
= 0,62-0,0025• V2-9.81-AH. |
||||
|
4 |
|
|
|
|
14*
212
Таким образом, имеется система трех нелинейных алгеб раических уравнений для определения трех неизвестных: v ,,v 2 и АН.
Полученные уравнения решаем методом последова тельных приближений. Сложив друг с другом первые два уравнения системы, получим:
5285-Х,v,2 +2912,5-X2v22 =542,1,
а из третьего уравнения -
V , - V 2 = 0 , 0 1 7 8 5 - V A H .
Первое приближение. Сначала разность напоров АН в сечении утечки находим так, как если бы возникшая утечка не изменила режим течения в трубопроводе (см. решение предыдущей задачи): АН = 265,85 м. Тогда
v i —V 2 =0,291 => v2 = v, -0,291 м/с. 5285-X,v,2 + 2912,5-X2(v, -0 ,2 9 l)2 =542,1.
1) Положим X,(l) = X2(l) = 0,02. Тогда из квадратного
уравнения 5285-0,02-v ,2 + 2912,5-0,02-(v, -0 .291)2 =542,1
или 3v,2 -0,582- v, -9 ,2 2 = 0 находим: v, = 1,853 м/с. |
И |
далее: v2 = 1,562 м/с. Проверяем, правильно ли выбраны |
А., |
и Х2: |
|
v, = l,853;/fe, = 1,853-0,7/(15• 10"6) = 86473; X, =0,0185; |
|
v, = 1,562 ;Ле2 =1,562-0,7/(15 -10“б) = 72893 ; Х2 =0,0193;
Полученные результаты говорят о том, что необходима кор рекция.
2) Положим теперь А./2^ =0,0185, |
А.2^ =0,0193 Тогда |
|
из квадратного уравнения |
|
|
5285 • 0,0185 |
• v ,2 + 2912,5 • 0,0193 • (v, - 0,291)2 = 542,1 |
|
или 3,043-V,2 - |
0 ,6 0 7 -V ,-9,973 = 0 |
находим: v, = 1,913 |
213
м/с. И далее: v2 = 1,622 м/с. Проверяем, правильно ли вы
браны Х 1 и Х2:
v, |
= |
1,913;Ле, = 1,913-0 ,7 /(l5 -10_6) = 89273; X, =0,0183; |
v, |
= |
1,622;Ле2 =1,622-0,7/(15-10‘6) = 75693; Х2 =0,0191. |
Видим, что достигнута достаточная точность. Итак, v, =1,913 м/с, v, = 1,622 м/с.
Подставив полученный результат в исходную систему
уравнений, получим: |
|
|
||
^150 + |
6 |
Л |
80000 |
1,9132 |
4’5 1 ° |
- (50 + ЛН) = 0,0183 |
|||
\ |
870-9,81 |
у |
0,7 |
2-9,81 |
откуда найдем: АН = 237,2 м. Этот результат меньше при нятого АН = 265,85 м. Найдем второе приближение.
Второе приближение. Полагаем АН = 237,2 м. Тогда
vi —V 2 = 0,01785-^237,2 => |
v2 = v, -0,275 |
м/с. |
|
5285-X,v,2 +2912,5-X2(v, -0 |
.275)2 = 542,1. |
|
|
Положим X,(0 =0,0183, X2*2) = 0,0191. Тогда из квад- |
|||
ратного уравнения |
|
|
|
5285 0,0183-v,2 + 2912,5-0,0191 |
(v, -0,275)2 = 542,1 |
||
или 3,043-v,2 -0,574 -V ,-9,982 |
= 0 |
находим: |
v, = 1,908 |
м/с. И далее: v2 = 1,633 м/с. Проверяем, правильно ли вы
браны и Х2:
v, = 1,908;Ле,=1,908-0,7/(15-10'6) = 89040; X, =0,0183;
v, = 1,633;Ле2 = 1,633-0,7/(15 -10_6) = 76207 ;Х2 =0,0190.
Полученные результаты говорят о том, что имеется доста
точная точность. |
|
Итак, v, - v 2 =0,275 м/с. Это означает, |
что расход q |
утечки равен: S0 -(v, - V2) = 3 ,1 4 - 0 ,7 2/ 4 - 0 , 2 7 5 |
= 0 , 1 0 5 8 M 3/ C . |
214
2
За 6 часов из трубопровода вытечет: 0,1058-6-3600=2285 м , то есть более чем в 20 раз больше, чем в предыдущей зада че.
86. Чертеж к задаче представлен на рис. 2.5.
Рис. 2.5. К решению задачи № 86
Между сечениями х = 23 и х = 24 км находится сечение х , чья высота равна 170 м плюс вакуумметрическая высота hB, соответствующая разряжению, возникающему в трубо
проводе:
h _ (0,1-0,03)-106
в ~ |
= 9,77 м. |
730-9,81 |
|
Для определения х, составляем пропорцию: |
|
х, - 2 3 |
17 0 + 9 ,7 7 -2 0 0 |
2 4 -2 3 |
х. = 23,405 км. |
730-9,81 |
Последнее означает, что между этими километрами опо рожнилось 405 м трубы.
Между сечениями 28 и 30 км также имеется сечение х2, чья геодезическая высота равна 180 м плюс hB. Для опреде
ления х2 составляем пропорцию: |
|
||
х2- 2 8 = 180 + 9 ,7 7 -7 5 |
= 29,996 км. |
||
3 0 - 2 8 _ |
1 9 0 -7 5 |
||
Х |
215
Это означает, здесь опорожнилось всего 4 м трубы.
Кроме того, опорожняются 2 км трубопровода от 25 до 27 км, т.е. 2000 м.
Таким образом, всего опорожнилось 2409 м. Объем Vn образовавшейся полости
V„ = (3,14• 0,3612/4 ) • 2409 = 246,45 м3.
87. Если бы отвод к нефтебазе был заполнен полно стью, то избыточное давление перед нефтебазой было бы не меньше давления, определяемого столбом жидкости между наивысшей точкой (х = 1,5 км) отвода и его концом, то есть
р в > p g (z15 - z 6) = 735-9,81 (180-50) = 0,937 МПа.
Однако по условию задачи избыточное давление в конце отвода составило 0,45 МПа, что свидетельствует о наличии в отводе пустот, рис. 2.6.
Образовавшиеся в отводе пустоты заполнены парами бензина с давлением 70 кПа, то есть в них существует ваку ум. Согласно формуле гидростатического распределения давления имеем (см. рис. 2.6):
рк = 70000+ 735 *9,81* hK= 98100+ 450000 .
Отсюда находим: h K=66,31 м. Следовательно, уровень z2 бензина в отводе имеет отметку 50 + 66,31 = 116,31 м.
0
216
Найдем соответствующую z2 координату х2 зеркала
бензина в отводе. Для этого составим пропорцию:
х2-3 ,0 |
11631 -150 |
3 3 -3 ,0 |
х , = 3,084 км. |
7 0 -1 5 0 |
Следовательно, здесь опорожнилось примерно 84 м трубы. Кроме того, опорожнилась часть трубы на участке меж
ду сечениями 1,5 и 2,0 км. По закону сообщающихся сосу
дов z„ = 150 |
м. Для определения соответствующей этой |
высоте координаты х, составляем пропорцию: |
|
хп- 1 3 |
150—180 |
2 ,0 -1 3 |
х, = 1,636 км. |
7 0 -1 8 0 |
Следовательно, здесь опорожнилось примерно 136 м трубы. Итого, за счет несанкционированного отбора топлива из отвода в нем опорожнилось 84 + 136 = 220 м. Это составля
ет примерно 220■ (3,14 ■ОД442/4 ) =3,58 м3.
8&. Сначала находим гидравлический уклон i на рас сматриваемом участке. Имеем:
Н„ = 100+ 4Л Ю ‘/(840 9,81)£ 646,1м,
Н . = 6 0 + 0 3 -1 0 S/(840 -9,81)= 96,4 M,
i = - Н, )/L = (646,1 -96,4)/125000 £ 4,4 10 3
шта 4,4 W/KM.
Затем рассчитываем напор Н# в месте аварии:
Н. = 646Д—4„4-56 = 399,7 м,
ишротидитаиорНд, = z , +pame/p g :
Нф= 180+ 10I30q f(840 -931)= 192,3 м.
Таким образом, на отверстии имеется перепад АН напоров, ршшыш 399,7 - 1923 = 207,4 м.
По формуле (53) рассчитываем расход q утечки: q =0,62-4-10"® -/2-9,81-207,4 = U82-10"* м’/с.
За 20 сутаж выплетет 1582-10* 3600-24-20=27337 м*.
217
89. Чертеж к этой задаче представлен на рис. 2.7.
Рис. 2.7. К решению задачи № 89
Судя по профилю рассматриваемого трубопровода, топливо могло вытекать только из двух участков: 1) 5,5 < х < 6,0 км и 2) 9,0 < х < 10,0 км.
Высотная отметка z, пробитого сечения находится ли
нейной интерполяцией точек профиля:
z„ —180 |
5,5 -6 ,0 |
160 -180 |
откуда z, = 176,7 м. |
3 ,0 -6 ,0 ’ |
Найдем сначала объем дизельного топлива, которое вы текло из трубопровода за первые 20 мин, в течение которых ПНУ еще работала. Расход q, через отверстие можно рас считать по формуле (5.3):
Р » -Р а
Pg
в которой р, - давление в месте аварии. Записав уравнение Бернулли для участка трубопровода между ПНУ и местом аварии, получим:
|
|
218 |
|
|
|
|
Г150-ь |
1’6 ' 1° 6 1 |
176,7 + Р* Ра™ |
= Х- |
5,5 IQ3 |
v2 |
|
|_ |
840-9,8 lj L |
Рё |
J |
|
0,146 |
2-9,81 |
где v = 4Q/rcd2 = 4 -80/(3600-3,14 |
0,1462) = 1,328 |
м/с. По |
||||
скольку |
число РейнольдсаRe = 1,328 0,14б/6-10"6 = 32315, |
то коэффициент X гидравлического сопротивления нахо дится по формуле Блазиуса:
Х = в’3164 5 0,0236. </32315
Подставив значения X и v в уравнение Бернулли, сначала найдем: (р* - р атм )/pg = 87,55 м, а затем и расход q , :
3,14 |
0,008 2 \ |
q, =0,62- |
• V2- 9,81-87,55 = 0,00129 м3/с. |
За 20 мин (или 1200 с) вытекло 1200 0,00129 = 1,548 м3топ лива.
Чтобы определить, как вытекало дизельное топливо по сле остановки ПНУ, используем ту же формулу (53), только в качестве движущего напора возьмем разность z ( t ) - z , ,
где z(t) — уровень жидкости на участке CD , а затем |
на |
||
участке АВ. Имеем: |
|
|
|
q2 =M -s-v/2g[z(t)-z,] = - S ~ |
= - S c t g a ~ , |
|
|
|
at |
at |
|
где S = 3,14 0,1462/4 = 0,0167 M 2 - площадь сечения трубо |
|||
провода; |
ctg a = 1000/(190-180) = 100-котангенс угла |
a |
|
наклона трубопровода на участке CD. В итоге получаем |
|||
дифференциальное уравнение |
|
|
|
dt |
S c tg a |
|
|
решение которого определяется интегралом
|
|
219 |
J |
dz |
_ |xs-y/2g { |
i l o j z - z . |
S c tg a |
Вычисляя интеграл в левой части этого равенства, по лучаем:
V z ( t ) - z . - V l 9 0 - z . = - W * * -t, 2 S c tg a
то есть зависимость z = z (t). Подставляя в полученную за
висимость t = 3600-6 = 2160 с, |
z, = 176,7 |
м, |1 = 0,62, |
||
S= 0,0167 м2, |
s = 5,024 10~5 м2, находим, |
что |
z(2160) = 184,29 |
|
м. Поскольку |
184,29 > z(C) = 180 |
м, то |
можно заключить, |
что зеркало опускающейся жидкости в трубопроводе не вы ходит за пределы участка CD.
Координата х зеркала находится из уравнения:
1 0 0 0 -х
----------------- = ctg a = 100. 190-184,29
Отсюда находим: х = 9429 м, то есть опорожнилось 571 м трубы на участке CD.
Объем опорожненной части трубопровода составляет: 0,0167 *571 = 9,5 м3. Таким образом, из трубы вытекло: 1,54 + 9,54 = 11,08 м3 топлива.
90. Профиль сегмента нефтепродуктопровода приведен на рис. 2.8. Сразу же после аварии в наивысшей точке про филя (х = 20км) происходит разрыв сплошности столба жидкости и в вершине образуется вакуум (упругостью на сыщенных паров керосина пренебрегаем). Воздух не может проникнуть в трубу через образовавшееся отверстие, по скольку на участке между 12 и 13 км имеется спуск.
Уравнение Бернулли, записанное для двух сечений неф тепродуктопровода: х, = 12 км и х = x (t), с учетом условий
при х, =12 км: z, =56 м, р, =0,10 МПа,
|
|
220 |
|
при х = x (t): z = z(t), p = 0,0 МПа |
|
||
дает следующее уравнение: |
|
|
|
flOOOOO |
ГЛ /Ч1 |
. х -1 2 0 0 0 |
v -Ivi |
^780-9,81 |
) |
0,205 |
2-9,81 |
в котором |
v - скорость |
столба нефтепродукта (v < 0), |
|
H = -v- |
|
|
|
Для того, чтобы определить высотную отметку z(t), со ответствующую свободной поверхности керосина в сечении x(t) учитывается, что наклон трубопровода на участке меж ду 14-м и 20-м километрами равен 0,002, вследствие чего уравнение его оси имеет вид:
z = 0,002 ( х - 14000)+ 60.
Исключив z из уравнения Бернулли для разности напоров, получим исходное расчетное уравнение:
0,008х -147,87 , 2 |
(*) |
----------------------= Л,- v |
х-12000
а.Сначала рассчитываем скорость жидкости в началь ный момент времени, когда ее свободная поверхность нахо
дится в точке х = 20000 м. Из (*) следует, что
к • v2 =0,0015, откуда методом последовательных прибли