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

ВЛИЯНИЕ СТЕПЕНИ ЗАПОЛНЕНИЯ ГОРИЗОНТАЛЬНЫХ ЦИЛИНДРИЧЕСКИХ РЕЗЕРВУАРОВ НА ВРЕМЯ ИХ ОПОРОЖНЕНИЯ ПРИ АВАРИИ

Д.В. Седов,

доцент кафедры автотехнической экспертизы и автоподготовки ФГКОУ ВПО ВСИ МВД России, кандидат технических наук

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

In the petrochemical industry for the storage and transportation ofpetroleum products are widely used horizontal cylindrical tanks. The scale of accidents resulting from the violation of their integrity, depend primarily on the rate of discharge of oil. Speed depends, in part, on the degree of filling of the reservoir, which determines the initial height of the liquid column. However, the regulatory impact assessment methodologies do not take into account the failures of the receptacle, that is assumed that the initial height of the liquid column is equal to the diameter of the tank. This leads to incorrect results, unnecessarily overstate or understate the effects of the disaster. In this paper, we propose a simple method of accounting for the value of degree completion .

Sedov D. Influence degree completion horizontal cylindrical tanks while they emptied by accident

При хранении и транспортировке нефтепродуктов существует опасность возникновения аварии из-за нарушения герметичности обечайки резервуаров (образование трещин, расхождение швов, износ уплотнений, расслоение биметалла, коррозия, пробои при падении резервуара с опор или обрушении каких-либо конструкций и т.д.). Тяжесть ее последствий зависит от скорости истечения нефтепродуктов на прилегающую территорию. На скорость истечения влияет не только диаметр аварийного отверстия, но и степень заполнения резервуара s, определяемая технологическим регламентом (она определяет начальную высоту столба жидкости Н0, от которой зависит скорость истечения). Если для вертикальных резервуаров постоянного сечения Н0 определяется просто (объем жидкости в резервуаре ^ = sV, где V - объем резервуара, следовательно, Н0 = sHpe3, где Нрез - высота резервуара), то для горизонтальных цилиндрических резервуаров, у которых сечение по высоте не постоянно, формулу Н0 = sHpe3 не применяют. Поэтому зачастую при оценке последствий аварии степень заполнения s не учитывают совсем, то есть принимают, что Н0 равно диа-

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

Время опорожнения горизонтального цилиндрического резервуара определяется по известной формуле (см., например, [1]):

3

X =____________,

оп 3Mo^dLJ2g ’

где L - длина резервуара;

D - диаметр резервуара;

Но - коэффициент расхода через аварийное отверстие (0,62 для отверстия в тонкой стенке);

dome - диаметр аварийного отверстия в днище резервуара (возникновение отверстия в днище - наиболее опасный вариант);

g - ускорение свободного падения.

Однако выражение (1) применимо в случае, когда резервуар заполнен на 100 %. При частичном заполнении резервуара необходимо изменить граничные условия дифференциального уравнения, из которого получена формула (1) (см. [1]):

ird 2

^0-^V2gHdx = -2Цн (D - H)dH,

где Н - высота столба жидкости в резервуаре.

В левой части уравнения (2) -объем нефтепродукта, вытекшего из резервуара за время dr, в правой -объем резервуара, освободившийся от нефтепродукта за это же время.

Уравнение (1) получено из (2) при следующих граничных условиях: в начальный момент времени (при т1 = 0) высота столба жидкости равна диаметру резервуара (Hj = D), в момент опорожнения (при т2 = топ) высота столба Н2 = 0.

В случае же частичного заполнения при т] = 0 начальная высота столба Н] будет равна не D, а некоторой величине Н0 < D. После интегрирования (2) в новых границах получаем выражение для времени опорожнения частично заполненного горизонтального цилиндрического резервуара:

16L

т =

оп

D2-

3^отв

Величина Н0, как правило, неизвестна и зависит только от степени заполнения s резервуара, которая задается технологическим регламентом. Зависимость Н0 = fs) можно получить из уравнения (2), записанного в виде:

где dV - объем нефтепродукта, вытекшего из резервуара за время dr.

В начальный момент времени Н] = Н0, а объем вытекшего нефтепродукта V] = 0. В конечный момент времени Н2 = 0, а V2 = ^ (Vm - объем жидкости в резервуаре). Интегрируя (4) в данных границах, получаем:

о Но

J dV = -2L |д/H (D - H)dH;

Vx: о

Но _________

Vx = 2L JVDН - H2 dH.

о

Интеграл

DН - H2 dH

явля-

ется интегралом вида Jv a + bx - ox2 dx

(см., например, [2]), решение которого выглядит следующим образом:

dV = -2U\H

JV a + bx - cx2dx =

2cx - b I : 2 b2 + 4ac

-----Va + bx - cx +-._

4c 8^7

2cx - b

arcsin^ + const.

Следовательно, искомый интеграл составит:

H о ____ 2 H - D i-------

JV DН - H2 dH = о--VDH о - H о

D

2

+-----arcsin

8

yjb + 4ac

. 2H о - D tD 2

D

+

16

2

а выражение (5) примет вид:

F = L

2Hо D DH - Hn 2 +

V

D 2

— arcsin-----—

4 D

2H - D tD

+

2

8

2

TD

Так как Vж = sV = s—^~L , то вместо (6) получаем:

= 2(2Я D)jвио _ но2 + I

tD t

+— arcsin-t

2 H - D

D

+ о,5.

s

Выражение (7) и есть зависимость Н0 = f(s) (в явном виде выразить Н0 не удается).

На примере горизонтального цилиндрического резервуара размерами L = 5 м и D = 2,5 м построен график данной функции (рис. 1). Из рисунка видно, что график функции (7) можно приближенно заменить прямой линией:

H 0 = sD

(прослеживается аналогия с вертикальными резервуарами, для которых Н0 = sHpe3). Как будет показано ниже, данная замена дает все же более точные результаты, чем расчеты, не учитывающие s

По графику функции (7) (см. рис. 1) определяем, что если s = 0,9, то Н0 = 2,15 м. Теперь по формуле (3), зная диаметр аварийного отверстия dome, можем найти время опорожнения резервуара.

Рис. 1. Зависимость начальной высоты столба жидкости Н0 от степени заполнения е горизонтального цилиндрического резервуара: 1 -функция (7);

2 - зависимость Н0 = eD

Рассмотрим наиболее часто образующиеся при авариях отверстия (dome = 0,03...0,1 м [3]). Так, при dome = 0,04 м нефтепродукт вытечет из резервуара через 2 ч. Расчеты же без учета s, то есть по фор-

муле (1), дают завышенный результат: 2,12 ч. Результаты с учетом (8) более точны: 2,06 ч.

На рис. 2 показан график зависимости времени опорожнения резервуара (3) от степени его заполнения

при dome = 0,04 м. Как видно из рисунка, если степень заполнения уменьшается, расхождения в результатах заметно увеличиваются. Например, при s = 0,8 (Н0 = 1,86 м) время выхода

нефтепродукта по формуле (3) с учетом зависимости (7) составляет 1,85 ч. Если вместо (7) использовать (8), то получается 1,93 ч.

Рис. 2. Зависимость времени опорожнения резервуара (3) от степени его заполнения (при dome = 0,04 м): 1 - с учетом зависимости (7);

2 - с учетом приближенной замены

На рис. 3 показан график зависимости времени опорожнения (3) от степени его заполнения s при различных значениях dome. Из рисунка видно, что если не учитывать s, то расхождения в результатах увеличиваются с уменьшением dome. Например, при s = 0,9 время опорожнения резервуара с учетом (7) для отверстий диаметром 0,1 м, 0,05 м и 0,03 м соответственно составляет 0,32 ч, 1,27 ч и 3,53 ч. Без учета s, то есть по формуле (1) время опорожнения равно 0,34 ч, 1,36 ч и 3,77 ч.

Однако относительная ошибка вычислений

8 = ^оп ~Топ (s)l .100

Г *оп

(топ - время опорожнения либо без учета s, либо с использованием приближения (8); топ(ё) - время опорожнения с учетом s по зависимости (7)) от размера отверстия не зависит и определяется только точностью учета степени заполнения (рис. 4).

Рис. 3. Зависимость времени опорожнения резервуара от степени его заполнения при различных диаметрах аварийного отверстия: 1 - 0,03 м; 2 - 0,04 м; 3 - 0,05 м;

4 - 0,07 м; 5 - 0,1 м; пунктиром показаны результаты соответствующих расчетов без

учета s, по формуле

Рис. 4. Относительные ошибки вычислений времени опорожнения горизонтального цилиндрического резервуара: 1 - по формуле (3) с учетом приближенной замены (8);

2 - без учета s, по формуле

Как видно из рис. 4, относительные ошибки расчетов с использованием замены (8) не превышают 4,5 % при s = 0.. .0,5. При меньших s (менее 0,4) погрешности превышают 5 % и быстро увеличиваются. Относительная погрешность вычислений без учета s уже при s = 0,9 достигает 7 %, при s = 0,8 она составляет 14,7 %, а при s= 0,5 составляет 54,7 %.

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

Увеличение абсолютной ошибки при уменьшении размера отверстия и степени заполнения резервуара наглядно прослеживается по графикам на рис 5.

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

Так, при образовании в рассматриваемом резервуаре пробоя с диаметром отверстия dome = 0,1 м время его опорожнения без учета степени заполнения пр^ = 0,9...0,8 оказывается завышенным на 0,02-0,04 ч.

При малых значениях dome и снижении s абсолютная погрешность, напротив, резко увеличивается. Так, если dome = 0,03 м, то при s = 0,9 время опорожнения резервуара становится больше на 0,23 ч, а при s = 0,8 время опорожнения будет завышено на 0,48 ч.

Рис. 5. Абсолютные ошибки вычислений времени опорожнения горизонтального цилиндрического резервуара при различных диаметрах аварийного отверстия: 1

ПРОБЛЕМЫ ПОЖАРНОЙ БЕЗОПАСНОСТИ

2(65) 2013

и 4 - 0,03 м; 2 и 5 - 0,05 м; 3 и 6 - 0,1 м; сплошная линия - по формуле

Таким образом, при исследовании влияния степени заполнения горизонтальных цилиндрических резервуаров на время их опорожнения получена точная аналитическая зависимость, которая позволяет избежать завышения или занижения теоретических последствий аварий. Обосновано применение приближенной зависимости Н0 = sD, которая, тем не менее, дает более точные результаты, чем расчеты без учета степени заполнения s. На примере горизонтального цилиндрического резервуара объемом 25 м3 показано, что неучет данной величины приводит к тому, что при наиболее распространенных размерах аварийных отверстий (0,03-0,1 м) и степенях заполнения резервуаров (0,9-0,8) время опорожнения завышается на 0,23-0,48 ч, относительная ошибка расчетов составляет 5,1-14,7 %. При увеличении объема резервуаров погрешность будет уменьшаться.

ПРИМЕЧАНИЯ

1. Мальцев Е. Д., Бубырь Н. Ф., Во-ротынцев Ю. П. и др. Гидравлика и пожарное водоснабжение. - М.: ВИПТШ МВД СССР, 1976. - 448 с.

2. Выгодский М. Я. Справочник по высшей математике. - 7-е изд. - М.: Наука,

1964. - 870 с.

3. Руководство по определению зон воздействия опасных факторов аварий со сжиженными газами, горючими жидкостями и аварийно-химически опасными веществами на объектах железнодорожного транспорта: Введ. в действ. 24.11.97 Указанием МПС России № Г-1362у / МПС России;

МЧС России. - М.: Гос. ун-т техникоэкономических изысканий и проектирования ж.-д. транспорта МПС России, 1997. - 198 с.