Руководство по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ»
(утв. Приказом Ростехнадзора от 20.04.2015 № 158)
Внимание
Документ утратил силу в связи с изданием Приказа Ростехнадзора от 02.11.2022 № 385, утвердившего новое Руководство.
Утверждено Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РУКОВОДСТВО ПО БЕЗОПАСНОСТИ «МЕТОДИКА МОДЕЛИРОВАНИЯ РАСПРОСТРАНЕНИЯ АВАРИЙНЫХ ВЫБРОСОВ ОПАСНЫХ ВЕЩЕСТВ»
I. ОБЩИЕ ПОЛОЖЕНИЯ
1. Руководство по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ» (далее — Руководство) разработано в целях содействия соблюдению требований Федеральных норм и правил в области промышленной безопасности «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств», утвержденных Приказом Ростехнадзора от 11 марта 2013 г. № 96 (зарегистрирован Минюстом России 16 апреля 2013 г., регистрационный № 28138), требований Федеральных норм и правил в области промышленной безопасности «Общие требования к обоснованию безопасности опасного производственного объекта», утвержденных Приказом Ростехнадзора от 15 июля 2013 г. № 306 (зарегистрирован Минюстом России 20 августа 2013 г., регистрационный № 29581), требований Федеральных норм и правил в области промышленной безопасности «Правила безопасности химически опасных производственных объектов», утвержденных Приказом Ростехнадзора от 21 ноября 2013 г. № 559 (зарегистрирован Минюстом России 31 декабря 2013 г., регистрационный № 30995).
2. Настоящее Руководство содержит рекомендации к расчетам зон аварийного распространения опасных веществ в атмосфере при оценке риска аварий для обеспечения требований промышленной безопасности при проектировании, строительстве, капитальном ремонте, техническом перевооружении, реконструкции, эксплуатации, консервации и ликвидации опасных производственных объектов, на которых возможны случаи выброса ОВ в атмосферу, и не является нормативным правовым актом.
3. Организации, выполняющие расчеты зон распространения ОВ в атмосфере при авариях, могут использовать иные обоснованные способы и методы, чем те, которые указаны в настоящем Руководстве.
4. В настоящем Руководстве используются обозначения и сокращения, а также термины и определения, приведенные в Приложениях № 1 и № 2 к настоящему Руководству.
5. Руководство распространяется на опасные производственные объекты, на которых возможны случаи выброса газообразных или жидких ОВ в однофазном или двухфазном состоянии с плотностью большей, чем плотность воздуха при соответствующих условиях.
II. РЕКОМЕНДУЕМЫЕ МЕТОДИЧЕСКИЕ ОСНОВЫ, ДОПУЩЕНИЯ И УСЛОВИЯ ВЫБРОСА ОПАСНЫХ ВЕЩЕСТВ
6. Расчеты распространения ОВ в атмосфере, приведенные в настоящем Руководстве, основаны на модели рассеяния «тяжелого» газа. Основными особенностями образования «тяжелого» газа являются: соотношение плотности газообразного ОВ и плотности воздуха более 1, низкая температура, наличие аэрозолей.
7. Модель «тяжелого» газа учитывает следующие процессы:
а) движение облака с учетом изменения скорости ветра по высоте;
б) гравитационное растекание облака;
в) рассеяние облака в вертикальном направлении за счет атмосферной турбулентности (подмешивание воздуха в облако);
г) рассеяние облака в горизонтальном направлении за счет подмешивания воздуха в облако, происходящего как за счет атмосферной турбулентности, так и за счет гравитационного растекания;
д) нагрев или охлаждение облака за счет подмешивания воздуха;
е) фазовые переходы ОВ в облаке («газ-жидкость» и «жидкость-газ»);
ж) теплообмен облака с подстилающей поверхностью.
8. В настоящем Руководстве приняты следующие допущения:
а) газообразное ОВ считается идеальным газом, свойства которого не зависят от температуры;
б) жидкое ОВ считается несжимаемой жидкостью, свойства которой не зависят от температуры;
в) гравитационное растекание облака ОВ учитывается с помощью эмпирической зависимости;
г) ОВ находится в аварийном оборудовании при давлении, равном или превышающем атмосферное;
д) истечение и испарение ОВ происходят с постоянной скоростью, соответствующей максимальной скорости истечения (испарения);
е) разлив жидкой фазы происходит на поверхности без впитывания;
ж) для случаев отсутствия обвалования толщина слоя разлившегося жидкого ОВ принимается равной 0,05 м <1>;
<1> При наличии достаточных обоснований допускается задание слоя разлития с глубиной, отличной от 0,05 м, в частности в соответствии с нормами пожарной безопасности «Определение категорий помещений, зданий и наружных установок по взрывопожарной и пожарной опасности», утвержденными Приказом Министерства Российской Федерации по делам гражданской обороны, чрезвычайным ситуациям и ликвидации последствий стихийных бедствий от 18.06.2003 г. № 314 (НПБ 105-03).
з) осаждение (конденсация) на подстилающую поверхность выброса газообразного ОВ и его химические превращения при рассеянии не учитываются;
и) рассматриваются только случаи аварийного истекания ОВ из оборудования, когда отверстия разгерметизации или места расположения аварийных трубопроводов всем сечением находятся выше или ниже уровня налива жидкого ОВ в оборудовании, соответственно, рассматривается истечение только газовой либо только жидкой фазы.
9. В настоящем Руководстве используются стандартные характеристики атмосферы и профили ветра, а также известные скорости подмешивания воздуха в выброс. Для описания устойчивости атмосферы используется 6 классов устойчивости — A, B, C, D, E и F (по Паскуилу). Первые три класса соответствуют неустойчивой стратификации атмосферы, последние два — устойчивой. Класс D соответствует нейтральной стратификации атмосферы. Предполагается, что в течение времени распространения облака характеристики атмосферы не меняются.
10. В настоящем Руководстве рассматриваются следующие условия изменения состояния ОВ:
а) проливы жидкой фазы имеют, как правило, форму квадрата, вдоль одной из сторон которого направлен ветер; в случае если фактические поперечные размеры пролива существенно отличаются от продольных, допускается принимать поперечный размер пролива отличным от принятого в Руководстве;
б) выброс происходит на уровне земли или площадки (этажерки), где расположено технологическое оборудование, рассеяние выброса проходит от уровня поверхности земли;
в) в начальный момент времени (на месте выброса) первичное облако ОВ имеет форму цилиндра (рисунок 3-1 Приложения № 3 к настоящему Руководству), а вертикальное сечение вторичного облака ОВ вдоль направления ветра во всех сценариях представляет собой прямоугольник (рисунок 3-2 Приложения № 3 к настоящему Руководству);
г) капли в облаках и пролив жидкости на подстилающую поверхность не «захолаживаются», то есть их температура не снижается ниже температуры кипения ОВ;
д) распространение выброса происходит над твердой ровной поверхностью, с которой нет обмена массой, а есть только обмен теплом;
е) учет подмешивания воздуха в облака на начальном этапе выброса осуществляется в соответствии с Приложением № 13 к настоящему Руководству. При отсутствии данных допускается принимать, что в начальный момент времени в облаках (первичном и вторичных) опасное вещество воздухом не разбавлено;
ж) в облаке существует фазовое равновесие газ-жидкость, это равновесие устанавливается мгновенно;
з) фазовые переходы ОВ приводят только к изменению высоты облака;
и) при определении размеров зон, где возможно горение (или детонация) топливно-воздушной смеси, предполагается, что горение (или детонация) может быть инициировано в областях со средней концентрацией от 0,5 нижнего концентрационного предела распространения пламени до верхнего концентрационного предела распространения пламени.
11. Руководство не рекомендуется применять либо рекомендуется применять с ограничениями в следующих случаях:
а) расчет рассеяния ОВ в штилевых условиях (отсутствие ветра);
б) расчет распространения выброса внутри помещений;
в) распространение выброса за пределами более 10 км от места выброса;
г) распространение выброса с массой ОВ в первичном облаке более 500 т;
д) распространение облаков от пролива ОВ со стороной квадрата пролива более 500 м;
е) наличие на пути движения облака препятствий, размеры которых превышают размеры облака;
з) выпадение конденсированной фазы на подстилающую поверхность и ее повторное поступление в атмосферу.
III. НЕКОТОРЫЕ ВОЗМОЖНЫЕ СЦЕНАРИИ ВЫБРОСА ОПАСНЫХ ВЕЩЕСТВ
12. Руководство рекомендует проводить расчеты для следующих сценариев выброса ОВ в зависимости от характера разрушения оборудования и агрегатного состояния ОВ.
Сценарий 1. Полное разрушение оборудования, содержащего газообразное ОВ.
Сценарий 2. Нарушение герметичности (частичное разрушение) оборудования, содержащего газообразное ОВ.
Сценарий 3. Полное разрушение оборудования, содержащего жидкое ОВ.
Сценарий 4. Нарушение герметичности (частичное разрушение) оборудования, содержащего жидкое ОВ.
По сценариям 1 и 3 ОВ поступает в окружающую среду мгновенно; по сценариям 2 и 4 ОВ поступает в окружающую среду через отверстие площадью S в течение продолжительного времени.
Сценарии 1 и 3 применимы только к емкостному оборудованию, сценарии 2 и 4 - как к емкостному оборудованию, так и к трубопроводам.
Возможные конфигурации оборудования и схемы его разрушения приведены в Приложении № 4 к настоящему Руководству, возможные стадии развития аварийных ситуаций приведены в Приложении № 5 к настоящему Руководству.
Приведенный перечень сценариев выброса не охватывает всего разнообразия возможных ситуаций, поэтому при выборе сценария для случаев, не перечисленных выше, рекомендуется руководствоваться соображениями физического подобия процессов.
При прогнозировании наибольших масштабов химического заражения и размеров зон, ограниченных концентрационными пределами распространения пламени ОВ, в качестве исходных данных рекомендуется принимать:
а) сценарий с полным разрушением емкости (технологической, складской, транспортной и др.), содержащей ОВ в максимальном количестве;
б) сценарий «гильотинного» разрыва трубопровода с максимальным расходом при максимальной длительности выброса;
в) метеорологические условия — класс устойчивости атмосферы — F, скорость ветра на высоте 10 м — 1 м/с.
IV. РЕКОМЕНДУЕМЫЕ ИСХОДНЫЕ ДАННЫЕ И ПОРЯДОК РАСЧЕТА ХАРАКТЕРИСТИК ВЫБРОСА ОПАСНЫХ ВЕЩЕСТВ
13. Исходными данными для расчета характеристик выбросов ОВ являются свойства ОВ, технологические данные с учетом параметров оборудования, в котором обращается ОВ, сценарии выброса ОВ в атмосферу, характеристики территории вблизи аварийного объекта, метеоусловия на момент аварии. Полный перечень данных, необходимых для проведения расчета по оценке последствий аварийных выбросов опасных веществ, а также полей концентрации и токсодоз приведен в Приложении № 6 к настоящему Руководству.
14. Основными величинами, расчет которых проводится в соответствии с настоящим Руководством, являются:
а) пространственно-временное распределение концентраций ОВ с учетом максимально достигаемой концентрации в данной точке на поверхности земли;
б) пространственные размеры зон достижения пороговой и смертельной токсодоз при ингаляционном воздействии на организм человека, размеры зон токсического поражения заданной вероятности по пробит-функциям, приведенным в таблицах № 7-1, 7-2 Приложения № 7 к настоящему Руководству;
в) размеры зон, ограниченных концентрационными пределами распространения пламени;
г) количество ОВ в облаке, способное участвовать во взрывных превращениях (горении и детонации);
д) количество ОВ, поступающее в окружающую среду в газовой и жидкой фазах;
е) количество ОВ, распространяющееся в атмосфере и выпадающее на подстилающую поверхность (при наличии жидкой фазы);
ж) площадь пролива и скорость испарения ОВ из пролива (при наличии жидкой фазы);
з) при продолжительном выбросе определяется скорость и длительность поступления ОВ в окружающую среду, масса капельных включений в облаках (первичном и вторичных), эффективные температура и плотность в облаках, геометрические характеристики облаков (эффективные высота и радиус (полуширина)), скорости распространения облаков, времена подхода и поражающего действия облаков.
15. В общем случае для выбранного i-го сценария рекомендуется рассчитывать следующие характеристики выброса:
Для сценария 4 при истечении из трубопровода, если количество ОВ в отсекаемом участке аварийного трубопровода составляет более 20% общей массы выброса, дополнительно рекомендуется рассчитывать:
q_{i}^{\text{отс}.\text{выб}},q_{\text{ж}i}^{0\text{отс}.\text{выб}},\rho_{i}^{\text{отс}.\text{выб}},t_{i}^{\text{отс}.\text{выб}},B_{i}^{\text{отс}.\text{выб}},H_{i}^{\text{отс}.\text{выб}} .
16. Для сценария 1 характеристики выброса рекомендуется рассчитывать по следующим Формулам:
если известна масса ОВ в оборудовании Q, или
Q_{1} = \frac{\mu}{R} \frac{V_{1}P_{1}}{T_{1}} , (2)
если неизвестна масса ОВ в оборудовании Q, но известны объем оборудования V1 давление в оборудовании P1 и температура в оборудовании T1:
q_{1}^{\text{ж}} = q_{1}^{\text{г}} = q_{1}^{\text{ги}} = q_{1}^{\text{и}} = q_{1}^{\text{е}} = q_{1}^{\text{отс}.\text{выб}} = q_{\text{ж}1}^{0\text{ж}} = q_{\text{ж}1}^{0\text{г}} = q_{\text{ж}1}^{0\text{ги}} = q_{\text{ж}1}^{0\text{и}} = q_{\text{ж}1}^{0\text{е}} = q_{\text{ж}1}^{0\text{отс}.\text{выб}} = 0 , <1> (4)
<1> Данные и нижеследующие аналогичные обнуления не имеют физического смысла и необходимы для формализации подстановок в формулах; равенство на какой-либо стадии нулю начальных размеров облака, плотности и расхода (или количества опасного вещества) в нем, а также длительности поступления опасного вещества в атмосферу означает отсутствие данной стадии развития аварийной ситуации.
t_{1}^{\text{ж}} = t_{1}^{\text{г}} = t_{1}^{\text{ги}} = t_{1}^{\text{и}} = t_{1}^{\text{е}} = t_{1}^{\text{отс}.\text{выб}} = 0 , (5)
\rho_{1}^{\text{выб}} = \rho_{1}(\frac{P_{0}}{P_{1}})^{\frac{1}{\gamma}} , <1> (6)
<1> При расчете по формуле (6), а также по формулам (18), (38) - (39) и (79) - (83) при наличии достоверных данных о начальной плотности выброса допускается задание иной величины с приведением соответствующих обоснований.
где \rho_{1} = Q_{1}/V_{1} — плотность газообразного ОВ в оборудовании:
\rho_{1}^{\text{ж}} = \rho_{1}^{\text{г}} = \rho_{1}^{\text{ги}} = \rho_{1}^{\text{и}} = \rho_{1}^{\text{е}} = \rho_{1}^{\text{отс}.\text{выб}} = 0 . (7)
При отсутствии данных о начальных размерах первичного облака рекомендуется принимать его радиус равным его высоте:
, (8)B_{1}^{\text{ж}} = B_{1}^{\text{г}} = B_{1}^{\text{ги}} = B_{1}^{\text{и}} = B_{1}^{\text{е}} = B_{1}^{\text{отс}.\text{выб}} = 0 , (9)
H_{1}^{\text{ж}} = H_{1}^{\text{г}} = H_{1}^{\text{ги}} = H_{1}^{\text{и}} = H_{1}^{\text{е}} = H_{1}^{\text{отс}.\text{выб}} = 0 . (10)
17. Для сценария 2 характеристики выброса рекомендуется рассчитывать по следующим формулам:
Q_{2} = Q_{\text{ж}2}^{0} = 0 (11)
докритическое истечение (при \frac{P_{0}}{P_{2}} \geq (\frac{2}{\gamma + 1})^{\frac{\gamma}{\gamma - 1}} ):
q_{2}^{\text{ги}} = 0,8S \sqrt{2 \frac{\gamma}{\gamma - 1} P_{2}\rho_{2}((\frac{P_{0}}{P_{2}})^{\frac{2}{\gamma}} - (\frac{P_{0}}{P_{2}})^{\frac{\gamma + 1}{\gamma}})} (12)
сверхкритическое истечение (при \frac{P_{0}}{P_{2}} <(\frac{2}{\gamma + 1})^{\frac{\gamma}{\gamma - 1}} ):
q_{2}^{\text{ги}} = 0,8S \sqrt{P_{2}\rho_{2}\gamma(\frac{2}{\gamma + 1})^{\frac{\gamma + 1}{\gamma - 1}}} . (13)
Если истечение происходит из трубопровода, на входе которого стоит компрессор, и S превосходит величину 0,2Sтр, то q_{2}^{\text{ги}} полагается равным расходу компрессора q^{\text{комп}} :
q_{2}^{\text{ж}} = q_{2}^{\text{г}} = q_{2}^{\text{и}} = q_{2}^{\text{е}} = q_{2}^{\text{отс}.\text{выб}} = q_{\text{ж}2}^{0\text{ж}} = q_{\text{ж}2}^{0\text{г}} = q_{\text{ж}2}^{0\text{ги}} = q_{\text{ж}2}^{0\text{и}} = q_{\text{ж}2}^{0\text{е}} = q_{\text{ж}2}^{0\text{отс}.\text{выб}} = 0 , (14)
t_{2}^{\text{ги}} = min((Q + Q_{\text{тр}})/q_{2}^{\text{ги}},t_{\text{отс}} + Q_{\text{тр}}/q_{2}^{\text{ги}},t_{\text{ликвид}}) , (15)
если известна масса ОВ в оборудовании Q, и
t_{2}^{\text{ги}} = min(\frac{\mu}{R} \frac{V_{2}P_{2}}{T_{2}q_{2}^{\text{ги}}} + Q_{\text{тр}}/q_{2}^{\text{ги}},t_{\text{отс}} + Q_{\text{тр}}/q_{2}^{\text{ги}},t_{\text{ликвид}}) , (16)
если неизвестна масса ОВ в оборудовании Q, но известны объем оборудования V2, давление в оборудовании P2 и температура в оборудовании T2:
t_{2}^{\text{ж}} = t_{2}^{\text{и}} = t_{2}^{\text{г}} = t_{2}^{\text{е}} = t_{2}^{\text{отс}.\text{выб}} = 0 , (17)
\rho_{2}^{\text{ги}} = \rho_{2}(\frac{P_{0}}{P_{2}})^{\frac{1}{\gamma}} , (18)
где \rho_{2} = \frac{\mu}{R} \frac{P_{2}}{T_{2}} — плотность газообразного ОВ в оборудовании:
\rho_{2}^{\text{выб}} = \rho_{2}^{\text{ж}} = \rho_{2}^{\text{и}} = \rho_{2}^{\text{г}} = \rho_{2}^{\text{е}} = \rho_{2}^{\text{отс}.\text{выб}} = 0 . (19)
При отсутствии данных о начальных размерах вторичного облака рекомендуется для прямоугольного сечения этого облака принимать его полуширину равной его высоте <1>:
<1> В формуле (20), а также далее в формулах (30), (34), (43), (53), (57), (60), (64), (68), (71), (73), (74), (77), (80), (81), (85), (86), (88) - (93) и в соотношениях после формулы (56) используется величина u_{0}\prime_{\text{эфф}i} , которая вычисляется согласно формуле (120).
H_{2}^{\text{ги}} = B_{2}^{\text{ги}} = \sqrt{\frac{q_{2}^{\text{ги}}}{2\rho_{2}^{\text{ги}}u_{0\text{эфф}2}^{\text{ги}}}} , (20)
R_{2} = B_{2}^{\text{ж}} = B_{2}^{\text{и}} = B_{2}^{\text{е}} = B_{2}^{\text{отс}.\text{выб}} = 0 , (21)
H_{2} = H_{2}^{\text{ж}} = H_{2}^{\text{и}} = H_{2}^{\text{е}} = H_{2}^{\text{отс}.\text{выб}} = 0 . (22)
18. Для сценария 3 характеристики выброса рекомендуется рассчитывать по следующим формулам:
Q_{3} = Q_{3}^{\text{г}} + Q_{3}^{\text{ж}} + Q_{3}^{\text{и}} + Q^{\text{г}} , (23)
Q^{\text{г}} = \alpha \frac{\mu}{R} \frac{V_{3}P_{3}}{T_{3}} , (24)
где \alpha — объемная доля оборудования, заполненная газовой фазой (формула (24) применяется, если заранее неизвестна величина Qг):
Q_{3}^{\text{г}} = Q^{\text{ж}}(1 - exp(- \frac{C_{p}(T_{3} - T_{\text{кип}} + |T_{3} - T_{\text{кип}}|)}{2\Delta H_{\text{кип}}})) , (25)
{}_{Q_{3} = min}^{\text{ж}}{_{Q_{3},Q - Q_{3}}^{\text{г ж г}}} , (26)
Q_{\text{ж}3}^{0} = Q_{3}^{\text{ж}} , (27)
Q_{3}^{\text{и}} = min\{ \frac{T_{\text{п}} - T_{\text{кип}} + |T_{\text{п}} - T_{\text{кип}}|}{\Delta H_{\text{кип}}} \sqrt{\frac{\lambda_{\text{п}}c_{\text{п}}\rho_{\text{п}}}{\pi}} \frac{F_{\text{кон}}^{2}}{F} \sqrt[\text{т}]{t_{\text{кип}}} ,Q^{\text{ж}} - Q_{3}^{\text{г}} - Q_{3}^{\text{ж}}\} , (28)
F = min\{ \frac{{}^{\cdot}Q^{\text{ж}} - Q_{3}^{\text{г}} - Q_{3}^{\text{ж}}}{0,05\rho_{\text{ж}}.} ,F_{\text{обв}}\} , <1> (29)
<1> При наличии достоверных исходных данных о геометрии и характеристиках поверхности, на которую происходит пролив, допускается задание иной площади пролива с приведением соответствующих обоснований.
где Fконт площадь контакта с твердой поверхностью, эта площадь включает как боковую поверхность обвалования, так и подстилающую поверхность; при проливе на неограниченную поверхность или внутри обвалования (без достижения проливом его границ) Fконт = F;
(30)p_{\text{н}} = 760exp(\Delta H_{\text{кип}}\mu(\frac{1}{T_{\text{кип}}} - \frac{1}{T_{\text{прив}}})R) (31)
{}_{T_{\text{прив}} = max}{_{\text{возд прол}}^{T,T}} (32)
(33)q_{3}^{\text{и}} = F \sqrt{\mu} 10^{- 6}(5,38 + 4,1u_{0\text{эфф}3}^{\text{и}})p_{\text{н}} , <1> (34)
<1> При расчете интенсивности испарения из пролива по формулам (34), (57), (60), (64) и (66), а также при использовании интенсивности испарения в качестве вспомогательной величины в формулах (30), (53), (65), (67), (68), (71), равно как и в четырех ненумерованных формулах после формулы (56), допускается при наличии достоверных данных задание иных интенсивностей испарения с приведением соответствующих обоснований.
q_{3}^{\text{ж}} = q_{3}^{\text{г}} = q_{3}^{\text{ги}} = q_{3}^{\text{е}} = q_{3}^{\text{отс}.\text{выб}} = q_{\text{ж}3}^{0\text{ж}} = q_{\text{ж}3}^{0\text{г}} = q_{\text{ж}3}^{0\text{ги}} = q_{\text{ж}3}^{0\text{и}} = q_{\text{ж}3}^{0\text{е}} = q_{\text{ж}3}^{0\text{отс}.\text{выб}} = 0 , (35)
t_{3}^{\text{и}} = (Q - Q_{3})q_{3}^{\text{и}} , (36)
t_{3}^{\text{ж}} = t_{3}^{\text{г}} = t_{3}^{\text{ги}} = t_{3}^{\text{е}} = t_{3}^{\text{отс}.\text{выб}} = 0 , (37)
(38)
(39)\rho_{3}^{\text{ж}} = \rho_{3}^{\text{г}} = \rho_{3}^{\text{е}} = \rho_{3}^{\text{ги}} = \rho_{3}^{\text{отс}.\text{выб}} = 0 . (40)
При отсутствии данных о начальных размерах первичного и вторичного облаков рекомендуется для первичного облака принимать начальный радиус равным его высоте, а полуширину вторичного облака — полуширине пролива:
, (41)B_{3}^{\text{и}} = 0,5 \sqrt{F} , (42)
H_{3}^{\text{и}} = q_{3}^{\text{и}}/(2u_{0\text{эфф}3}^{\text{и}}B_{3}^{\text{и}}\rho_{3}^{\text{и}}) , (43)
B_{3}^{\text{ж}} = B_{3}^{\text{отс}.\text{выб}} = B_{3}^{\text{г}} = B_{3}^{\text{ги}} = B_{3}^{\text{е}} = 0 , (44)
H_{3}^{\text{ж}} = H_{3}^{\text{отс}.\text{выб}} = H_{3}^{\text{г}} = H_{3}^{\text{ги}} = H_{3}^{\text{е}} = 0 . (45)
19. Для сценария 4 характеристики выброса рекомендуется рассчитывать по следующим формулам.
Если истечение происходит из трубопровода, на входе которого стоит емкость, и величина S превосходит 0,2Sтр либо истечение происходит непосредственно из емкости, то расход определяется по формуле:
q_{\text{выб}} = 0,6Ssign(Q_{H}^{\text{ж}}) \sqrt{2Hg\rho_{\text{ж}}^{2} + 2\rho_{\text{ж}}(P_{4} - p_{\text{н}}(T_{4})) + \frac{1}{K} \frac{\Delta H_{\text{кип}}^{2}\rho^{2}(T_{4},p_{\text{н}}(T_{4}))}{C_{p}T_{\text{кип}}}} , (46)
где p_{\text{н}}(T_{4}) = P_{0}exp(\Delta H_{\text{кип}}\mu(\frac{1}{T_{\text{кип}}} - \frac{1}{T_{4}})R) , (46-1)
\rho(T_{4},p_{\text{н}}(T_{4})) = \frac{\mu}{R} \frac{p_{\text{н}}(T_{4})}{T_{4}} — плотность газообразного ОВ при температуре T4 и давлении pн(T4); (46-2)
K — функция, зависящая от L — длины участка трубопровода от входа до места разгерметизации (в случае истечения непосредственно из емкости L = 0).
(47)Если истечение происходит из трубопровода, на входе которого стоит насос, а величина S превосходит 0,25Sтр, то qвыб полагается равным расходу насоса qнасос.
В остальных случаях расход определяется по формуле:
q_{\text{выб}} = sign(Q_{H}^{\text{ж}})0,6S\rho_{\text{ж}} \sqrt{2Hg + 2 \frac{P_{4} - P_{0}}{\rho_{\text{ж}}}} . (48)
В случае если истечение происходит через трубопровод и количество ОВ в отсекаемом участке аварийного трубопровода составляет более 20% общей массы выброса, то рекомендуется рассмотреть истечение после блокировки аварийного участка. Скорость выброса на этой стадии будет составлять:
. (49)Если истечение из трубопровода после отсечения аварийного участка не рассматривается, то q_{\text{выб}}^{\text{отс}} = 0 .
(50)Q_{\text{ж}4}^{0} = q_{4}\prime\prime t\prime + q_{4}\prime\prime_{\text{отс}}t_{\text{о}}\prime_{\text{ст}} , (51)
(52)
(53)
(54)t\prime = min\{t_{\text{кип}},t_{\text{отс}},t_{\text{ликвид}},t_{\text{и}}\prime_{\text{сп}}, \frac{Q_{H}^{\text{ж}} + Q_{H^{\text{тр}}}^{\text{ж}}}{q_{\text{выб}}} \} , (55)
t_{\text{о}}\prime_{\text{тс}} = min\{t_{\text{кип}} - t\prime,t_{\text{ликвид}} - t\prime,t_{\text{и}}\prime_{\text{сп}} - t\prime,t_{\text{и}}\prime\prime_{\text{сп}} - t\prime, \frac{Q_{H^{\text{отс}.\text{тр}}}^{\text{ж}}}{q_{\text{выб}}^{\text{отс}}} \} , (56)
где Fконт — площадь контакта с твердой поверхностью, эта площадь включает как боковую поверхность обвалования, так и подстилающую поверхность; при проливе на неограниченную поверхность или внутри обвалования (без достижения проливом его границ) Fконт = F';
t_{\text{и}}\prime_{\text{сп}} = \frac{(q_{\text{выб}} - q_{4}\prime - q_{4}\prime\prime) \frac{Q_{H}^{\text{ж}} + Q_{H^{\text{тр}}}^{\text{ж}}}{q_{\text{выб}}}}{F \sqrt{\mu} 10^{- 6}(5,38 + 4,1u_{0\text{эфф}4}^{0})p_{\text{н}}} при проливе в обвалование
t_{\text{и}}\prime_{\text{сп}} = \frac{0,05\rho_{\text{ж}}}{\sqrt{\mu} 10^{- 6}(5,38 + 4,1u_{0\text{эфф}4}^{0})p_{\text{н}}} при проливе на неограниченную поверхность;
при проливе в обвалованиеt_{\text{и}}\prime\prime_{\text{сп}} = \frac{0,05\rho_{\text{ж}}}{\sqrt{\mu} 10^{- 6}(5,38 + 4,1u_{0\text{эфф}4}^{0})p_{\text{н}}} при проливе на неограниченную поверхность;
где p_{\text{н}} = 760exp(\Delta H_{\text{кип}}\mu(\frac{1}{T_{\text{кип}}} - \frac{1}{T_{\text{прив}}})R) — давление насыщенных паров ОВ при приведенной температуре Tприв, которая принимается равной максимуму из температуры воздуха Tвозд и температуры жидкой фазы в проливе в начальный момент времени Tпрол: Tприв = max(Tвозд, Tпрол). Величина Tпрол принимается аналогично формуле (33): Tпрол = Tкип, при T4 > Tкип и Tпрол = T4 в противных случаях;
{}_{4}^{q^{\text{ж}} = min}{_{\text{выб }4}^{q,q}\prime_{4}^{+ q}\prime\prime^{+ F} \sqrt{\mu} {}_{0\text{эфф}4\text{н}}^{10^{- 6}(5,38 + 4,1u^{\text{ж}})p}} , (57)
q_{4}\prime = q_{\text{выб}}(1 - exp(- \frac{C_{p}(T_{4} - T_{\text{кип}} + |T_{4} - T_{\text{кип}}|)}{2\Delta H_{\text{кип}}})) , (58)
{}_{q_{4}}\prime\prime_{= q_{\text{ж}4} = min}^{0\text{ж}}{_{q_{4}}\prime_{,q_{\text{выб}} - q_{4}}\prime} , (59)
{}_{4}^{q^{\text{отс}.\text{выб}} = min}{_{\text{выб }4}^{q^{\text{отс}},q}\prime_{\text{отс }4}^{+ q}\prime\prime_{\text{отс}}^{+ F} \sqrt{\mu} {}_{0\text{эфф}4\text{н}}^{10^{- 6}(5,38 + 4,1u^{\text{отс}.\text{выб}})p}} , (60)
q_{4}\prime_{\text{отс}} = q_{\text{выб}}^{\text{отс}}(1 - exp(- \frac{C_{p}(T_{4} - T_{\text{кип}} + |T_{4} - T_{\text{кип}}|)}{2\Delta H_{\text{кип}}})) , (61)
{}_{q_{4}}\prime\prime_{{}_{\text{отс}} = q_{\text{ж}4} = min}^{0\text{отс}.\text{выб}}{_{q_{4}\prime_{\text{отс}},q_{\text{выб}} - q_{4}\prime_{\text{отс}}}^{\text{отс}}} , (62)
q_{\text{ж}4}^{0\text{г}} = q_{\text{ж}4}^{0\text{ги}} = q_{\text{ж}4}^{0\text{и}} = q_{\text{ж}4}^{0\text{е}} = 0.0 , (63)
(64)
(65)q_{4}^{\text{и}} = F \sqrt{\mu} 10^{- 6}(5,38 + 4,1u_{0\text{эфф}4}^{\text{и}})p_{\text{н}} + 5,38S_{max} \sqrt{\mu} 10^{- 6}p_{\text{н}} , (66)
q_{4}^{\text{е}} = 5,38S_{max} \sqrt{\mu} 10^{- 6}p_{\text{н}} , (67)
где Smax — эффективная площадь эмиссии из разгерметизированного оборудования (определяется с учетом максимальной площади свободной поверхности в оборудовании Sобор) находится из соотношения:
{}_{max}^{S = min}{_{\text{обор }40\text{эфф}4}^{S,\rho^{\text{е}}u^{\text{е}}S}(\sqrt{\mu} {}_{\cdot 5,38 \cdot 10p_{\text{н}}}^{- 6})} , (68)
t_{4}^{\text{ж}} = min\{ \frac{Q_{H}^{\text{ж}} + Q_{H^{\text{тр}}}^{\text{ж}}}{q_{\text{выб}}} - t\prime,t_{\text{отс}} - t\prime\} , (69)
t_{4}^{\text{отс}.\text{выб}} = min\{t_{\text{ликвид}} - t_{\text{о}}\prime_{\text{тс}} - t_{4}^{\text{ж}} - t\prime, \frac{Q_{H^{\text{отс}.\text{тр}}}^{\text{ж}}}{q_{\text{выб}}^{\text{отс}}} - t_{\text{о}}\prime_{\text{тс}}\} , (70)
(71)Если пролив происходит в обвалование, то площадь пролива F совпадет с площадью обвалования, в противном случае площадь пролива определяется по следующей формуле:
, (72)
(73)
(74)
(75)t_{max}^{\text{ги}} = t_{max}^{\text{г}} - t_{4}^{\text{г}} , (76)
(77)
, (78)
(79)
(80)
(81)\rho_{4}^{\text{г}} = \rho_{4}^{\text{ги}} = \rho_{4}(\frac{P_{0}}{P_{4}})^{\frac{1}{\gamma}} , (82)
где \rho_{4} = \frac{\mu}{R} \frac{P_{4}}{T_{4}} .
(83)При отсутствии данных о начальных размерах первичного облака рекомендуется принимать его радиус равным его высоте:
. (84)При отсутствии данных о начальных размерах вторичных облаков, формирующихся при наличии пролива, рекомендуется для облаков принимать полуширину вторичного облака равной полуширине пролива либо в отсутствие пролива:
(85)
(86)B_{4}^{\text{г}} = B_{4}^{\text{и}} = 0,5 \sqrt{F} , (87)
H_{4}^{\text{ж}} = q_{4}^{\text{ж}}/(2u_{0\text{эфф}4}^{\text{ж}}B_{4}^{\text{ж}}\rho_{4}^{\text{ж}}) , (88)
H_{4}^{\text{отс}.\text{выб}} = q_{4}^{\text{отс}.\text{выб}}/(2u_{0\text{эфф}4}^{\text{отс}.\text{выб}}B_{4}^{\text{отс}.\text{выб}}\rho_{4}^{\text{отс}.\text{выб}}) , (89)
H_{4}^{\text{г}} = q_{4}^{\text{г}}/(2u_{0\text{эфф}4}^{\text{г}}B_{4}^{\text{г}}\rho_{4}^{\text{г}}) , (90)
H_{4}^{\text{и}} = q_{4}^{\text{и}}/(2u_{0\text{эфф}4}^{\text{и}}B_{4}^{\text{и}}\rho_{4}^{\text{и}}) . (91)
При отсутствии данных о начальных размерах вторичных облаков, формирующихся в отсутствие пролива, рекомендуется для прямоугольного сечения вторичного облака принимать его полуширину равной его высоте.
H_{4}^{\text{ги}} = B_{4}^{\text{ги}} = \sqrt{\frac{q_{4}^{\text{ги}}}{2\rho_{4}^{\text{ги}}u_{0\text{эфф}4}^{\text{ги}}}} , (92)
H_{4}^{\text{е}} = B_{4}^{\text{е}} = \sqrt{\frac{q_{4}^{\text{е}}}{2\rho_{4}^{\text{е}}u_{0\text{эфф}4}^{\text{е}}}} . (93)
V. РЕКОМЕНДАЦИИ ПО РАСЧЕТУ ПОЛЕЙ КОНЦЕНТРАЦИИ И ТОКСОДОЗ
20. Для условий, в которых происходит выброс ОВ, определяются характерный размер шероховатости поверхности zпов, класс устойчивости атмосферы, характеристика профиля ветра \alpha_{\text{в}} , масштаб Монина — Обухова LМО, динамическая скорость u*, скорость подмешивания воздуха u_{\text{подм}}^{\text{верх}} и коэффициент дисперсии в поперечном направлении \sigma_{y} .
21. Если характерный размер шероховатости не может быть задан, исходя из реальных метеорологических условий с приведением соответствующих обоснований, то он задается согласно данным, приведенным в таблице № 7-3 Приложения № 7 к настоящему Руководству.
Если в результате расчета окажется, что вертикальный размер облака с требуемыми характеристиками меньше, чем характерный размер шероховатости, то это означает, что результаты расчетов по данному Руководству можно рассматривать лишь в качестве оценочных.
22. Если класс устойчивости не может быть задан, исходя из реальных метеорологических условий с приведением соответствующих обоснований, то класс устойчивости задается с использованием данных, приведенных в таблице № 7-4 Приложения № 7 к настоящему Руководству, в зависимости от скорости ветра и интенсивности теплового потока у поверхности (инсоляция и облачность).
Для расчета наихудшего варианта принимается класс устойчивости — F и скорость ветра — 1 м/с.
23. Для характерного размера шероховатости zпов, класса устойчивости и скорости ветра u10 на высоте z10 определяется коэффициент \alpha_{\text{в}} в соответствии с данными, приведенными в таблице № 7-5 <1> Приложения № 7 к настоящему Руководству.
<1> При наличии достоверных исходных данных о характере изменения скорости ветра с высотой допускается задание иной величины \alpha_{\text{в}} с приведением соответствующих обоснований.
При проведении расчетов первоначально рекомендуется выбирать \alpha_{\text{в}} для высоты до 20 м (первое значение см. таблицу № 7-5 Приложения № 7 к настоящему Руководству).
Если в результате расчета окажется, что вертикальный размер облака с требуемыми характеристиками больше 20 м, то рекомендуется провести новый расчет с измененным \alpha_{\text{в}} для высоты до или свыше 50 м в зависимости от рассчитанного вертикального размера облака.
24. Для характерного размера шероховатости zпов и выбранного класса устойчивости определяется масштаб Монина — Обухова LМО. L_{\text{МО}} = \infty для нейтральной устойчивости атмосферы (класс D). Для остальных условий устойчивости атмосферы масштаб Монина — Обухова определяется по формуле (94) с использованием данных, приведенных в таблице № 7-6 Приложения № 7 к настоящему Руководству:
L_{\text{МО}} = k_{L} \cdot z_{\text{пов}}^{p} . (94)
25. Динамическая скорость определяется по формуле:
u_{*} = \frac{k \cdot u_{10}}{(ln((z_{10} + z_{\text{пов}})/z_{\text{пов}}) - \varphi)} , (95)
где \varphi задается исходя из класса устойчивости атмосферы согласно (96):
(96)26. Скорость подмешивания воздуха в облако через верхнюю границу определяется в зависимости от характеристик облака (эффективной плотности \rho_{\text{эфф}} , эффективной высоты Hэфф, эффективной температуры Tэфф и эффективной теплоемкости облака Cэфф) <1> по следующей формуле:
<1> В следующих далее формулах используются переменные \rho_{\text{эфф}} , Hэфф, Tэфф, Eпов и Cэфф, для первичного облака эти переменные равны \rho_{\text{эфф}i} , Hэффi, Tэффi, Eповi и Cэффi, а для вторичных облаков — \rho_{\text{эфф}i}^{l} , H_{\text{эфф}i}^{l} , T_{\text{эфф}i}^{l} , E_{\text{эфф}i}^{l} и C_{\text{эфф}i}^{l} .
u_{\text{подм}}^{\text{верх}} = \frac{k \cdot u_{t}}{\text{Ф }(Ri_{*})} , (97)
, (98)Ri_{*} = g \cdot \frac{(\rho_{\text{эфф}} - \rho_{\text{возд}})}{\rho_{\text{возд}}} \cdot \frac{H_{\text{эфф}}}{u_{t}^{2}} , (99)
где u_{t} = \sqrt{u_{*}^{2} + (b \cdot w_{*})^{2}} ,b = 0,2,w_{*} = [\frac{g \cdot |E_{\text{пов}}| \cdot H_{\text{эфф}}}{\rho_{\text{эфф}} \cdot T_{\text{эфф}} \cdot C_{\text{эфф}}}]1/3 , Eпов — удельный тепловой поток (от поверхности земли в облако), описываемый ниже ((114) - (117) или (132) - (135)).
27. Коэффициент дисперсии в поперечном направлении \sigma_{y} определяется по следующей формуле с помощью данных, приведенных в таблице № 7-7 Приложения № 7 к настоящему Руководству:
\sigma_{y}(x) = \delta \cdot x \cdot (1 + \gamma_{y} \cdot x)^{- 1/2} , (100)
где \gamma_{y} = 0,0001\text{м}^{- 1} , \delta = \delta_{600} \cdot (t_{av}/600) .
28. Для каждой из стадий выброса по сценарию i-му определяются поля концентрации ОВ и максимальная концентрация ОВ на оси x.
Для первичного облака концентрация ОВ и размеры облаков при их рассеянии вычисляются по формулам (101) - (117).
29. Состояние первичного облака в каждый момент времени характеризуется следующими параметрами:
cцi, Szi, Syi, ri, Eэффi, xцi, Qжi, \rho_{\text{эфф}i} (или Vэффi) и Tэффi.
Они дополняются еще четырьмя, которые рассчитываются на основе введенных выше переменных <1>: Rэффi, Hэффi, uэффi, Qсумi.
<1> В случае если величина Hэффi меньше 0,5 м, то в формуле (103) для расчета uэффi используется величина Szi, полученная по формуле (102) в предположении, что Hэффi = 0,5 м.
R_{\text{эфф}i} = \sqrt{r_{i}^{2} + S_{yi}^{2}} (101)
H_{\text{эфф}i} = \frac{1}{\beta} \cdot \text{Г}(\frac{1}{\beta}) \cdot S_{zi} (102)
u_{\text{эфф}i} = \frac{\text{Г}[(1 + \alpha_{\text{в}})/\beta]}{\text{Г}[1/\beta]} \cdot u_{10} \cdot (\frac{S_{zi}}{z_{10}})\alpha_{\text{в}} (103)
Q_{\text{сумм}i} = \pi \cdot R_{\text{эфф}i}^{2} \cdot H_{\text{эфф}i}/V_{\text{эфф}i} (104)
Схема первичного облака изображена на рисунке 3-1 Приложения № 3 к настоящему Руководству.
30. Распределение концентрации ОВ в облаке рекомендуется описывать зависимостями:
c_{i}^{0}(x,y,z,t) = c_{\text{ц}i} \cdot exp[- [\frac{z}{S_{zi}}]^{\beta}],\text{при }(x - x_{\text{ц}i})^{2} + y^{2}<r_{i}^{2} (105)
(106)31. Для определения пространственного распределения концентрации, профиль которой задан в пункте 30 настоящего Руководства, с помощью перечисленных параметров в пункте 29 настоящего Руководства используются следующие уравнения.
Сохранение массы выброшенного вещества Qi:
Q_{i} = c_{\text{ц}i} \cdot \pi \cdot R_{\text{эфф}i}^{2} \cdot H_{\text{эфф}i} . (107)
\frac{d}{dt} [Q_{\text{сум}i}] = \pi \cdot R_{\text{эфф}i}^{2} \cdot \rho_{\text{возд}} \cdot u_{\text{подм}}^{\text{верх}} + 2 \cdot \pi \cdot R_{\text{эфф}i} \cdot H_{\text{эфф}i} \cdot \rho_{\text{возд}} \cdot \gamma_{\text{подм}} \cdot \frac{d}{dt} [R_{\text{эфф}i}] . (108)
Гравитационное растекание облака:
\frac{d}{dt} [R_{\text{эфф}i}] = u_{\text{эфф}i} \frac{d}{dx} [R_{\text{эфф}i}] = C_{e} \sqrt{g \cdot H_{\text{эфф}i} \cdot [1 - \frac{\rho_{\text{возд}}}{\rho_{\text{эфф}i}}]} . (109)
Боковое рассеяние выброса за счет атмосферной диффузии:
\frac{d}{dt} [S_{yi}] = \frac{u_{\text{эфф}i}}{S_{yi}} \cdot 2 \cdot (2/\pi)^{1/2} \cdot (r_{i} + 1/2 \cdot \pi^{1/2} \cdot S_{yi}) \cdot \frac{d}{dx} [\sigma_{y}],\text{при }r_{i}>0 (110)
S_{yi}(x) = 2^{1/2} \cdot \sigma_{y}(x + x_{t}),\text{при }r_{i} = 0 . (111)
Сохранение энергии в облаке Eэффi:
(112)Способ расчета удельного теплового потока от подстилающей поверхности в облако Eповi приведен ниже в пункте 32 настоящего Руководства.
\frac{d}{dt} [x_{\text{ц}i}] = u_{\text{эфф}i} . (113)
Величины \rho_{\text{эфф}i}(V_{\text{эфф}i}) и Tэффi, а также величина Qжi вычисляются согласно подходу, изложенному в Приложении № 8 к настоящему Руководству.
32. Определение удельного теплового потока (от поверхности земли в облако):
(114)E_{\text{пов}.\text{вын}i} = 1.22 \cdot (\frac{u_{*}^{2}}{u_{10}})^{2} \cdot \rho_{\text{эфф}i} \cdot C_{\text{эфф}i} \cdot (T_{\text{пов}.\text{рас}} - T_{\text{эфф}i}) (115)
C_{\text{эфф}i} = \frac{((Q_{i} - Q_{\text{ж}i}) \cdot C_{p\text{газ}} + Q_{\text{ж}i} \cdot C_{p} + (Q_{\text{сум}i} - Q_{i}) \cdot C_{p\text{возд}})}{Q_{\text{сум}i}} (116)
E_{\text{пов}.\text{ест}i} = 3.5 \cdot 10^{- 3} \cdot (\frac{(T_{\text{пов}.\text{рас}} - T_{\text{эфф}i})^{2}}{0.5 \cdot (T_{\text{пов}.\text{рас}} + T_{\text{эфф}i})})^{23} \cdot \frac{P_{0}}{R} \cdot (g)^{1/3} . (117)
Для длительных выбросов концентрация ОВ и размеры облаков при их рассеянии вычисляются по формулам (118) - (135).
33. Дрейф вторичного облака рассматривается для следующих шести стадий распространения выброса ОВ (см. Приложение № 5 к настоящему Руководству):
истечение жидкой фазы (в стабильном или перегретом состоянии) до отсечения аварийного участка;
истечение жидкой фазы (в стабильном или перегретом состоянии) из аварийного участка после его отсечения (для сценария 4);
истечение газа при наличии пролива жидкой фазы и испарение с пролива;
истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы;
испарение с пролива при отсутствии истечения жидкости или газа из разрушенного оборудования;
испарение из емкости при отсутствии пролива.
Для каждой из этих стадий рассчитывается свое вторичное облако. При расчете каждой из этих стадий задаются свой расход ОВ в шлейфе q_{i}^{l} и расход жидкой фазы в начальном сечении q_{\text{ж}i}^{0l} :
q_{i}^{\text{ж}},q_{i}^{\text{отс}.\text{выб}},q_{i}^{\text{г}},q_{i}^{\text{ги}},q_{i}^{\text{и}},q_{i}^{\text{е}},q_{\text{ж}i}^{0\text{ж}},q_{\text{ж}i}^{0\text{отс}.\text{выб}},q_{\text{ж}i}^{0\text{г}},q_{\text{ж}i}^{0\text{ги}},q_{\text{ж}i}^{0\text{и}},q_{\text{ж}i}^{0\text{е}} ,
свой начальный размер облака (полуширина B_{i}^{l} и высота H_{i}^{l} ):
B_{i}^{\text{ж}},B_{i}^{\text{отс}.\text{выб}},B_{i}^{\text{г}},B_{i}^{\text{ги}},B_{i}^{\text{и}},B_{i}^{\text{е}}\text{и }H_{i}^{\text{ж}},H_{i}^{\text{отс}.\text{выб}},H_{i}^{\text{г}},H_{i}^{\text{ги}},H_{i}^{\text{и}},H_{i}^{\text{е}} .
Ранее эти величины для разных сценариев были рассчитаны по формулам (4), (8) - (10), (12), (14), (20) - (22), (34) - (35), (41) - (45), (57), (60), (64) - (67).
Состояние вторичного облака в каждом поперечном сечении характеризуется следующими параметрами:
c_{\text{ц}i}^{l};S_{zi}^{l};S_{yi}^{l};b_{i}^{l};E_{\text{эфф}i}^{l},x_{\text{п}i}^{l},x_{\text{з}i}^{l},q_{\text{ж}i}^{l} .
Кроме того, облако характеризуется такими параметрами, как \rho_{\text{эфф}i}^{l}(V_{\text{эфф}i}^{l}) и T_{\text{эфф}i}^{l} . Они дополняются еще четырьмя — B_{\text{эфф}i}^{l} ; H_{\text{эфф}i}^{l} ; u_{\text{эфф}i}^{l} ; q_{\text{сум}i}^{l} , которые рассчитываются на основе введенных выше переменных по формулам (118) - (121) <1>:
<1> В случае если величина H_{\text{эфф}i}^{l} меньше 0,5 м, то в формуле (120) для расчета u_{\text{эфф}i}^{l} используется величина S_{zi}^{l} , полученная по формуле (119) в предположении, что H_{\text{эфф}i}^{l} = 0,5\text{м} .
B_{\text{эфф}i}^{l} = b_{i}^{l} + 0,5 \cdot \sqrt{\pi} \cdot S_{yi}^{l} , (118)
H_{\text{эфф}i}^{l} = \frac{1}{\beta} \cdot \text{Г}(\frac{1}{\beta}) \cdot S_{zi}^{l} , (119)
u_{\text{эфф}i}^{l} = \frac{\text{Г}[(1 + \alpha_{\text{в}})/\beta]}{\text{Г}[1/\beta]} \cdot u_{10} \cdot (\frac{S_{zi}^{l}}{z_{10}})\alpha_{\text{в}} , (120)
q_{\text{сум}i}^{l} = 2 \cdot B_{\text{эфф}i}^{l} \cdot H_{\text{эфф}i}^{l} \cdot u_{\text{эфф}i}^{l}/V_{\text{эфф}i}^{l} . (121)
Схема распространения вторичного облака изображена на рисунке 3-2 Приложения № 3 к настоящему Руководству.
34. Распределение концентрации ОВ во вторичном облаке описывается зависимостями:
c_{i}^{l}(x,y,z) = c_{\text{ц}i}^{l} \cdot exp[- [\frac{z}{S_{zi}^{l}}]^{\beta}],\text{при }|y|<b_{i}^{l}\text{и }x_{\text{з}i}^{l}<x<x_{\text{п}i}^{l} (122)
c_{i}^{l}(x,y,z) = c_{\text{ц}i}^{l} \cdot exp[- [\frac{z}{S_{zi}^{l}}]^{\beta}] \cdot exp[- [\frac{|y| - b_{i}^{l}}{S_{yi}^{l}}]^{2}],\text{при }|y| \geq b_{i}^{l}\text{и }x_{\text{з}i}^{l}<x<x_{\text{п}i}^{l} . (123)
При x_{\text{п}i}^{l}<x или x<x_{\text{з}i}^{l}c_{i}^{l}(x,y,z) = 0 .
35. Для определения пространственного распределения концентрации, профиль которой задан в пункте 34 настоящего Руководства, с помощью перечисленных параметров в пункте 33 настоящего Руководства используются следующие уравнения:
сохранение массы выброшенного вещества q_{i}^{l} :
q_{i}^{l} = 2 \cdot c_{\text{ц}i}^{l} \cdot B_{\text{эфф}i}^{l} \cdot H_{\text{эфф}i}^{l} \cdot u_{\text{эфф}i}^{l} ; (124)
изменение расхода в шлейфе q_{\text{сум}i}^{l} :
\frac{d}{dx} [q_{\text{сум}i}^{l}] = 2 \cdot B_{\text{эфф}i}^{l} \cdot \rho_{\text{возд}} \cdot u_{\text{подм}}^{\text{верх}} + 2 \cdot H_{\text{эфф}i}^{l} \cdot \rho_{\text{возд}} \cdot \gamma_{\text{подм}} \cdot \frac{d}{dt} [B_{\text{эфф}i}^{l}] ; (125)
гравитационное растекание облака:
\frac{d}{dx} [B_{\text{эфф}i}^{l}] = \frac{C_{e}}{u_{\text{эфф}i}^{l}} \sqrt{g \cdot H_{\text{эфф}i}^{l} \cdot [1 - \frac{\rho_{\text{возд}}}{\rho_{\text{эфф}i}^{l}}]} ; (126)
боковое рассеяние выброса за счет атмосферной диффузии:
\frac{d}{dx} [S_{yi}^{l}] = \frac{1}{S_{yi}^{l}} \cdot 2 \cdot (2/\pi)^{1/2} \cdot (b_{i}^{l} + 1/2 \cdot \pi^{1/2} \cdot S_{yi}^{l}) \cdot \frac{d}{dx} [\sigma_{y}],\text{при }b_{i}^{l}>0 (127)
S_{yi}^{l}(x) = 2^{1/2} \cdot \sigma_{y}(x + x_{t}),\text{при }b_{i}^{l} = 0 ; (128)
сохранение энергии в облаке E_{\text{эфф}i}^{l} :
(129)Способ расчета удельного теплового потока E_{\text{пов}i}^{l} приведен в пункте 36 настоящего Руководства.
Положение переднего края облака x_{\text{п}i}^{l} определяют по формулам:
(130)После окончания l-й стадии выброса определяется положение заднего края x_{\text{з}i}^{l} :
(131)Величины \rho_{\text{эфф}i}^{l}(V_{\text{эфф}i}^{l}) и T_{\text{эфф}i}^{l} , а также величина q_{\text{ж}i}^{l} вычисляются согласно рекомендациям, изложенным в Приложении № 8 к настоящему Руководству.
36. Определение удельного теплового потока проводится по следующим формулам:
. (132)E_{\text{пов}.\text{вын}i}^{l} = 1.22 \cdot (\frac{u_{*}^{2}}{u_{10}})^{2} \cdot \rho_{\text{эфф}i}^{l} \cdot C_{\text{эфф}i}^{l} \cdot (T_{\text{пов}.\text{рас}} - T_{\text{эфф}i}^{l}) . (133)
C_{\text{эфф}i}^{l} = \frac{((q_{i}^{l} - q_{\text{ж}i}^{l}) \cdot C_{p\text{газ}} + q_{\text{ж}i}^{l} \cdot C_{p} + (q_{\text{сум}i}^{l} - q_{i}^{l}) \cdot C_{p\text{возд}})}{q_{\text{сум}i}^{l}} . (134)
E_{\text{пов}.\text{ест}i}^{l} = 3.5 \cdot 10^{- 3} \cdot (\frac{(T_{\text{пов}.\text{рас}} - T_{\text{эфф}i}^{l})^{2}}{0.5 \cdot (T_{\text{пов}.\text{рас}} + T_{\text{эфф}i}^{l})})^{2/3} \cdot \frac{P_{0}}{R} \cdot (g)^{1/3} . (135)
37. Концентрацию ОВ в точке в момент времени t при i-ом сценарии рекомендуется определять по формуле:
(136)38. Определяются максимально возможная концентрация ОВ на расстоянии x от места аварии при i-ом сценарии и время tmax i(x, y, z), когда достигается эта концентрация. По tmax i(x, y, z) определяется стадия (или сочетания стадий распространения выброса ОВ), на которой (которых) достигаются максимальные концентрации ОВ:
c_{i}^{max}(x,y,z) = ma_{t}x(c_{i}(x,y,z,t)) . (137)
39. Путем интегрирования по времени концентрации ci(x,y,z,t) определяется поле токсодозы Di(x,y,z), а также распределение максимальной токсодозы, достигаемой на заданном расстоянии от места выброса Dimax(x). Максимальная токсодоза для заданного расстояния за все время экспозиции достигается на оси y = 0, z = 0. При необходимости определяются составляющие токсодозы, соответствующие облакам (первичному и вторичным), образовавшимся на различных стадиях аварии, D_{i}^{0}(x,y,z) и D_{i}^{l}(x,y,z) , и максимальные токсодозы, достигаемые на заданном расстоянии от облаков, образовавшихся на разных стадиях аварии, D_{imax}^{0}(x) и D_{imax}^{l}(x) .
40. Сравнением с пороговыми и смертельными токсодозами (см. таблицу № 7-1 Приложения № 7 к настоящему Руководству) определяются расстояния, соответствующие смертельному поражению и пороговому воздействию. Для оценки вероятности смертельного поражения человека рекомендуется использовать пробит-функцию Pr, по которой с использованием таблиц № 7-1, 7-2 Приложения № 7 к настоящему Руководству определяется вероятность смертельного поражения человека на открытом пространстве.
Величина Pr определяется по следующей формуле:
Pr = a_{Pr} + b_{Pr}ln(\int_{0}^{t_{\text{эксп}}} (c_{i}(x,y,z,t))^{n_{Pr}}dt) , (138)
где коэффициенты a, b, n берутся из данных таблицы № 7-1.
VI. РЕКОМЕНДАЦИИ ПО РАСЧЕТУ ГРАНИЦ ЗОН, ОГРАНИЧЕННЫХ КОНЦЕНТРАЦИОННЫМИ ПРЕДЕЛАМИ РАСПРОСТРАНЕНИЯ ПЛАМЕНИ
41. Для взрывопожароопасных выбросов в момент времени t0 определяются поверхности, ограничивающие в пространстве области ВКПР (\sum \text{ВКПР }) и 0,5 НКПР (\sum 0.5\text{НКПР }) .
Граница области ВКПР \sum \text{ВКПР} определяется уравнением:
(139)Граница зоны 0,5 НКПР \sum 0.5\text{НКПР} определяется уравнением:
(140)Область в пространстве, где возможны воспламенение и горение (детонация) пожаровзрывоопасного вещества, определяется как огибающая поверхностей \sum 0.5\text{НКПР} за все моменты времени t0 существования в пространстве концентраций выше 0,5 НКПР.
42. Для взрывопожароопасных выбросов определяются размеры зон, на которые может дрейфовать выброс, сохраняя способность к воспламенению. Полагается, что этот размер соответствует достижению средних концентраций 0,5 НКПР.
Граница зоны достижения 0,5 НКПР на уровне z = z0 в момент времени t0 определяется соотношением:
(141)Граница зоны достижения 0,5 НКПР в вертикальной плоскости, перпендикулярной ветру (x = x0), в момент времени t0 определяется соотношением:
(142)Граница зоны достижения 0,5 НКПР в вертикальной плоскости, параллельной ветру (y = y0), в момент времени t0 определяется соотношением:
(143)Линия, ограничивающая в соответствующей плоскости (z = z0, y = y0 или x = x0) область, где возможны воспламенение и горение пожаровзрывоопасного вещества, определяется как огибающая профилей линий (141) - (143) за все моменты времени t0 существования в пространстве концентраций выше 0,5 НКПР.
43. Для взрывопожароопасных выбросов в момент времени t определяется масса топлива, находящаяся во взрывоопасных пределах и способная участвовать в процессах горения или детонации. Эта масса определяется путем интегрирования концентрации по пространству, ограниченному поверхностями, образованными границей, на которой концентрация облака достигает НКПР и ВКПР:
Q_{\text{вз}} = \int_{V_{\text{взр}}} \int_{V_{\text{взр}}} \int_{\text{взр}} c(x,y,z,t_{0})dxdydz , (144)
где {}_{\text{взр}}^{V =}{_{(x,y,z)|_{\text{НКПР }0\text{ВКПР}}}^{c \leq c(x,y,z,t) \leq c}} .
Если в результате расчета по формуле (144) в первичном облаке во взрывоопасных пределах окажется масса больше 10% всей массы топлива, находящейся в первичном облаке, то масса топлива во взрывоопасных пределах первичного облака принимается равной 10% всей массы топлива, находящейся в первичном облаке.
Приложение № 1 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
СПИСОК ОБОЗНАЧЕНИЙ И СОКРАЩЕНИЙ
В настоящем Руководстве применены следующие обозначения и сокращения:
НКПР — нижний концентрационный предел распространения пламени;
ВКПР — верхний концентрационный предел распространения пламени;
B_{i}^{l} — полуширина начального сечения вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l B_{i}^{l} = B_{i}^{\text{ж}} , B_{i}^{\text{отс}.\text{выб}} , B_{i}^{\text{г}} , B_{i}^{\text{ги}} , B_{i}^{\text{и}} , B_{i}^{\text{е}} ), м;
B_{\text{эфф}i}^{l} — эффективная полуширина вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l B_{\text{эфф}i}^{l} = B_{\text{эфф}i}^{\text{ж}} , B_{\text{эфф}i}^{\text{отс}.\text{выб}} , B_{\text{эфф}i}^{\text{г}} , B_{\text{эфф}i}^{\text{ги}} , B_{\text{эфф}i}^{\text{и}} , B_{\text{эфф}i}^{\text{е}} ), м;
Ce — коэффициент, использующийся при расчете скорости гравитационного растекания, равен 1,15;
Cp — теплоемкость жидкого опасного вещества, Дж/кг/К;
Cp возд — теплоемкость воздуха при постоянном давлении, Дж/кг/К;
Cp газ — теплоемкость газообразного опасного вещества, Дж/кг/К;
Cv — теплоемкость газообразного опасного вещества, Дж/кг/К;
Cv возд — теплоемкость воздуха при постоянном объеме, Дж/кг/К;
Cэффi — эффективная теплоемкость вещества в первичном облаке в i-ом сценарии, Дж/кг/К;
C_{\text{эфф}i}^{l} — эффективная теплоемкость вещества во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l C_{\text{эфф}i}^{l} = C_{pi}^{\text{ж}} , C_{pi}^{\text{отс}.\text{выб}} , C_{pi}^{\text{г}} , C_{pi}^{\text{ги}} , C_{pi}^{\text{и}} , C_{pi}^{\text{е}} ), Дж/кг/К;
Di — ингаляционная токсодоза за рассматриваемое время экспозиции в некоторой точке в i-м сценарии, кг·с/м3;
D_{i}^{0} — составляющая ингаляционной токсодозы за рассматриваемое время экспозиции в некоторой точке при рассеянии первичного облака в i-м сценарии, кг·с/м3;
D_{i}^{l} — составляющая ингаляционной токсодозы за рассматриваемое время экспозиции в некоторой точке при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l D_{i}^{l} = D_{i}^{\text{ж}} , D_{i}^{\text{отс}.\text{выб}} , D_{i}^{\text{г}} , D_{i}^{\text{ги}} , D_{i}^{\text{и}} , D_{i}^{\text{е}} ), кг·с/м3;
Di max — ингаляционная токсодоза в некоторой точке на оси y = 0, z = 0 за рассматриваемое время экспозиции в i-м сценарии, кг·с/м3;
D_{imax}^{0} — составляющая ингаляционной токсодозы в некоторой точке на оси y = 0, z = 0 за рассматриваемое время экспозиции при рассеянии первичного облака в i-м сценарии, кг·с/м3;
D_{imax}^{l} — составляющая ингаляционной токсодозы за рассматриваемое время экспозиции в некоторой точке на оси y = 0, z = 0 при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l D_{imax}^{l} = D_{imax}^{\text{ж}} , D_{imax}^{\text{отс}.\text{выб}} , D_{imax}^{\text{г}} , D_{imax}^{\text{ги}} , D_{imax}^{\text{и}} , D_{imax}^{\text{е}} ), кг·с/м3;
Dтр — диаметр трубопровода, м;
Eпов i — удельный тепловой поток от подстилающей поверхности в первичное облако в i-ом сценарии, Дж/м2;
E_{\text{пов }i}^{l} — удельный тепловой поток от подстилающей поверхности во вторичное облако, образующееся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l E_{\text{пов }i}^{l} = E_{\text{пов }i}^{\text{ж}} , E_{\text{пов }i}^{\text{отс}.\text{выб}} , E_{\text{пов }i}^{\text{г}} , E_{\text{пов }i}^{\text{ги}} , E_{\text{пов }i}^{\text{и}} , E_{\text{пов }i}^{\text{е}} ), Дж/м2/с;
Eпов вын i _ удельный тепловой поток за счет вынужденной конвекции от подстилающей поверхности в первичное облако в i-ом сценарии, Дж/м2;
E_{\text{пов вын }i}^{l} — удельный тепловой поток за счет вынужденной конвекции от подстилающей поверхности во вторичное облако, образующееся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l E_{\text{пов вын }i}^{l} = E_{\text{пов вын }i}^{\text{ж}} , E_{\text{пов вын }i}^{\text{отс}.\text{выб}} , E_{\text{пов вын }i}^{\text{г}} , E_{\text{пов вын }i}^{\text{ги}} , E_{\text{пов вын }i}^{\text{и}} , E_{\text{пов вын }i}^{\text{е}} ), Дж/м2/с;
Eпов ест i — удельный тепловой поток за счет естественной конвекции от подстилающей поверхности в первичное облако в i-ом сценарии, Дж/м2;
E_{\text{пов ест }i}^{l} — удельный тепловой поток за счет естественной конвекции от подстилающей поверхности во вторичное облако, образующееся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l E_{\text{пов ест }i}^{l} = E_{\text{пов ест }i}^{\text{ж}} , E_{\text{пов ест }i}^{\text{отс}.\text{выб}} , E_{\text{пов ест }i}^{\text{г}} , E_{\text{пов ест }i}^{\text{ги}} , E_{\text{пов ест }i}^{\text{и}} , E_{\text{пов ест }i}^{\text{е}} ), Дж/м2/с;
Eэффi — эффективная внутренняя энергия в первичном облаке в i-ом сценарии, Дж;
E_{\text{эфф}i}^{l} — эффективный поток внутренней энергии во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l E_{\text{эфф}i}^{l} = E_{\text{эфф}i}^{\text{ж}} , E_{\text{эфф}i}^{\text{отс}.\text{выб}} , E_{\text{эфф}i}^{\text{г}} , E_{\text{эфф}i}^{\text{ги}} , E_{\text{эфф}i}^{\text{и}} , E_{\text{эфф}i}^{\text{е}} ), Дж/с;
F — площадь поверхности пролива жидкого опасного вещества, м2;
F' — площадь поверхности пролива жидкого опасного вещества при образовании первичного облака в сценарии 4, м2;
Fконт — площадь контакта жидкого опасного вещества с подстилающей поверхностью при проливе, м2;
Fобв — площадь обвалования, при неограниченном проливе принимается большой (бесконечной) величиной, м2;
\Delta H_{\text{кип}} — теплота испарения (кипения) жидкого опасного вещества, Дж/кг;
H — высота столба жидкости опасного вещества в оборудовании над уровнем отверстия, через которое происходит истечение, м; при истечении из трубопровода, на входе которого стоит насос, принимается равной 0 м;
Hi — высота первичного облака опасного вещества в начальный момент времени (на месте выброса) в i-м сценарии, м;
Hэффi — эффективная высота первичного облака в i-ом сценарии, м;
H_{\text{эфф}i}^{l} — эффективная высота вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l H_{\text{эфф}i}^{l} = H_{\text{эфф}i}^{\text{ж}} , H_{\text{эфф}i}^{\text{отс}.\text{выб}} , H_{\text{эфф}i}^{\text{г}} , H_{\text{эфф}i}^{\text{ги}} , H_{\text{эфф}i}^{\text{и}} , H_{\text{эфф}i}^{\text{е}} ), м;
H_{i}^{l} — высота начального сечения вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l H_{i}^{l} = H_{i}^{\text{ж}} , H_{i}^{\text{отс}.\text{выб}} , H_{i}^{\text{г}} , H_{i}^{\text{ги}} , H_{i}^{\text{и}} , H_{i}^{\text{е}} ), м;
H_{\text{тр}}^{\text{отс}} — высота столба жидкого опасного вещества в трубопроводе над уровнем отверстия, через которое происходит истечение, м; берется на момент отсечения аварийного участка трубопровода; если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса, принимается равной 0 м;
K — безразмерная вспомогательная величина, зависящая от длины участка трубопровода L от входа до места разгерметизации;
L — длина участка трубопровода от его начала до места разгерметизации, м;
LМО — масштаб Монина — Обухова, м;
LCt50 — смертельная токсодоза, мг·мин/л;
Pi — давление в оборудовании в i-м сценарии, Па; при истечении из трубопровода, на входе которого стоит насос (компрессор), принимается равной давлению на выходе насоса (компрессора);
P0 — давление в окружающей среде, кПа; при нормальных условиях принимается равным 101,325 кПа;
PCt50 — пороговая токсодоза, мг·мин/л;
Q — общая масса опасного вещества в оборудовании, включает массу жидкости и массу газа; при выбросах из трубопровода с насосом (компрессором) на входе задается равной бесконечной величине; при выбросе из трубопровода, на котором возможно отсечение аварийного участка, не включает массу в отсекаемом участке аварийного трубопровода, кг;
Qг — масса газообразного опасного вещества в оборудовании, кг;
Qж — масса жидкого опасного вещества в оборудовании (при истечении из трубопровода с насосом на входе равно Q_{H}^{\text{ж}} ), кг;
Q_{3}^{\text{г}} — масса опасного вещества, переходящая в газовую фазу в первичное облако при мгновенном вскипании перегретого опасного вещества в сценарии 3, кг;
Q_{3}^{\text{ж}} — масса опасного вещества, переходящая в аэрозоль в первичное облако в сценарии 3, кг;
Q_{3}^{\text{и}} — масса опасного вещества, переходящая в газовую фазу в первичное облако при кипении пролива в сценарии 3, кг;
Qi — масса опасного вещества (включая жидкую и газообразную фазы), образующая первичное облако в i-м сценарии, кг;
Q_{i}^{\text{разб}} — масса воздуха, подмешиваемого в первичное облако в i-м сценарии на начальной стадии, кг;
Q_{H}^{\text{ж}} — масса жидкого опасного вещества в оборудовании выше уровня отверстия, через которое происходит истечение, кг; задается на момент времени, соответствующий началу аварии; при истечении из трубопровода с насосом на входе принимается равной бесконечной величине; если аварийное отверстие выше уровня жидкости, то величина принимается равной 0 кг; если истечение происходит из трубопровода, присоединенного к емкости, то Q_{H}^{\text{ж}} включает и массу жидкости в трубопроводе от емкости до начала отсекаемого аварийного участка;
Q_{H^{\text{отс тр}}}^{\text{ж}} — масса жидкого опасного вещества в отсекаемом участке аварийного трубопровода выше уровня отверстия, через которое происходит истечение, кг; задается на момент отсечения аварийного участка трубопровода; при истечении из трубопровода в сценарии 4 принимается равной 0 кг, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса; если аварийное отверстие находится выше уровня жидкости, то величина принимается равной 0 кг;
Q_{H^{\text{тр}}}^{\text{ж}} — масса жидкого опасного вещества в отсекаемом участке аварийного трубопровода выше уровня отверстия, через которое происходит истечение, кг; величина задается на момент начала аварии; при истечении из трубопровода в сценарии 4 принимается равной 0 кг, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
Qвз — масса опасного вещества, находящаяся во взрывоопасных пределах (распределенная в пространстве между поверхностями с концентрациями, равными НКПР и ВКПР) и способная участвовать в процессах горения и взрыва, кг;
Qжi — масса опасного вещества, находящаяся в первичном облаке в жидкой фазе (каплях) в i-ом сценарии, кг;
Q_{\text{ж}i}^{0} — масса опасного вещества, находящаяся в первичном облаке в жидкой фазе (каплях) в i-ом сценарии в начальный момент (на месте выброса), кг;
Q_{\text{отс}}^{\text{г}} — масса газообразного опасного вещества в аварийном участке трубопровода на момент отсечения, кг; при истечении из трубопровода в сценарии 4 принимается равной 0 кг, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
Qсумi — эффективная масса первичного облака в i-ом сценарии, включает газообразное опасное вещество, опасное вещество в жидкой фазе и воздух, кг;
Qтр — общая масса опасного вещества в отсекаемом участке аварийного трубопровода, кг; включает массу жидкости в сценарии 4 или массу газа в сценарии 2; при выбросах опасного вещества из емкости задается равной нулю;
Q_{\text{тр}}^{\text{ж}} — масса жидкого опасного вещества, способного вместиться в отсекаемом участке трубопровода; величина задается на момент начала аварии; определяется перемножением объема, занимаемого жидкостью в отсекаемом участке трубопровода, на плотность жидкого опасного вещества, кг;
Q_{\text{тр}^{\text{отс}}}^{\text{ж}} — масса жидкого опасного вещества в отсекаемом участке аварийного трубопровода на момент отсечения, кг; при истечении из трубопровода в сценарии 4 принимается равной 0 кг, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
R — универсальная газовая постоянная, равная 8,3144 Дж/К/моль;
Ri — радиус первичного облака опасного вещества в начальный момент времени (на месте выброса) в i-м сценарии, м;
Rэффi — эффективный радиус первичного облака в i-ом сценарии, м;
S — площадь аварийного отверстия, м2;
Smax — эффективная площадь эмиссии из разгерметизированного оборудования, м2;
Syi — горизонтальная дисперсия при рассеянии первичного облака в i-м сценарии, м;
S_{yi}^{l} — горизонтальная дисперсия при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l S_{yi}^{l} = S_{yi}^{\text{ж}} , S_{yi}^{\text{отс}.\text{выб}} , S_{yi}^{\text{г}} , S_{yi}^{\text{ги}} , S_{yi}^{\text{и}} , S_{yi}^{\text{е}} ), м;
Szi — вертикальная дисперсия при рассеянии первичного облака в i-м сценарии, м;
S_{zi}^{l} — вертикальная дисперсия при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l S_{zi}^{l} = S_{zi}^{\text{ж}} , S_{zi}^{\text{отс}.\text{выб}} , S_{zi}^{\text{г}} , S_{zi}^{\text{ги}} , S_{zi}^{\text{и}} , S_{zi}^{\text{е}} ), м;
Sобор — максимальная площадь поверхности жидкости внутри оборудования, м2;
Sтр — площадь поперечного сечения трубопровода, м2;
T_{i}^{0} — температура в облаке в начальный момент времени (на месте выброса), К;
Ti — температура, при которой находится опасное вещество внутри оборудования в i-м сценарии, К;
Tэффi — эффективная температура среды в первичном облаке в i-ом сценарии, К;
T_{\text{эфф}i}^{l} — эффективная температура среды во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l T_{\text{эфф}i}^{l} = T_{\text{эфф}i}^{\text{ж}} , T_{\text{эфф}i}^{\text{отс}.\text{выб}} , T_{\text{эфф}i}^{\text{г}} , T_{\text{эфф}i}^{\text{ги}} , T_{\text{эфф}i}^{\text{и}} , T_{\text{эфф}i}^{\text{е}} ), К;
Tвозд — температура воздуха, К;
Tкип — температура кипения жидкого опасного вещества при давлении P0, К;
Tп — температура подстилающей поверхности, на которую происходит пролив жидкой фазы опасного вещества, К;
Tпов рас — температура подстилающей поверхности, над которой происходит рассеяние выброса, К;
Tприв — вспомогательная (приведенная) температура, используемая при расчете скорости испарения, К;
Tпрол — температура жидкой фазы в проливе в начальный момент времени, К;
Vi — объем оборудования в i-м сценарии, м3, при выбросе с трубопровода, на входе которого стоит компрессор (насос), принимается равным бесконечной величине;
Vэффi — эффективный удельный объем среды в первичном облаке в i-ом сценарии, м3/кг;
V_{\text{эфф}i}^{l} — эффективный удельный объем среды во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l V_{\text{эфф}i}^{l} = V_{\text{эфф}i}^{\text{ж}} , V_{\text{эфф}i}^{\text{отс}.\text{выб}} , V_{\text{эфф}i}^{\text{г}} , V_{\text{эфф}i}^{\text{ги}} , V_{\text{эфф}i}^{\text{и}} , V_{\text{эфф}i}^{\text{е}} ), м3/кг;
a — безразмерная вспомогательная величина при расчете динамической скорости;
aPr — коэффициент при расчете пробит-функции;
b — вспомогательный коэффициент при расчете числа Ричардсона;
b_{i}^{l} — полуширина ядра вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l b_{i}^{l} = b_{i}^{\text{ж}} , b_{i}^{\text{отс}.\text{выб}} , b_{i}^{\text{г}} , b_{i}^{\text{ги}} , b_{i}^{\text{и}} , b_{i}^{\text{е}} ), м;
bPr — коэффициент при расчете пробит-функции;
ci — концентрация опасного вещества в воздухе в некоторой точке в некоторый момент времени в i-м сценарии, кг/м3;
c_{i}^{0} — концентрация опасного вещества в некоторой точке в некоторый момент времени при рассеянии первичного облака в i-м сценарии, кг/м3;
c_{i}^{l} — концентрация опасного вещества в некоторой точке в некоторый момент времени при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l c_{i}^{l} = c_{i}^{\text{ж}} , c_{i}^{\text{отс}.\text{выб}} , c_{i}^{\text{г}} , c_{i}^{\text{ги}} , c_{i}^{\text{и}} , c_{i}^{\text{е}} ), кг/м3;
cвкпр — концентрация топлива (пожаровзрывоопасного вещества) в смеси с воздухом, соответствующая ВКПР, кг/м3;
cп — теплоемкость подстилающей поверхности, на которую проливается опасное вещество, Дж/кг/К;
cнкпр — концентрация топлива (пожаровзрывоопасного вещества) в смеси с воздухом, соответствующая НКПР, кг/м3;
cцi — концентрация опасного вещества в центре облака в некоторый момент времени при рассеянии первичного облака в i-м сценарии, кг/м3;
c_{\text{ц}i}^{l} — концентрация опасного вещества в центре (на оси) облака в некоторый момент времени при рассеянии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l c_{\text{ц}i}^{l} = c_{\text{ц}i}^{\text{ж}} , c_{\text{ц}i}^{\text{отс}.\text{выб}} , c_{\text{ц}i}^{\text{г}} , c_{\text{ц}i}^{\text{ги}} , c_{\text{ц}i}^{\text{и}} , c_{\text{ц}i}^{\text{е}} ), кг/м3;
eвозд — удельная внутренняя энергия подмешиваемого воздуха, Дж/кг;
g — ускорение свободного падения, равно 9,81 м/с2;
kL — вспомогательный коэффициент для вычисления масштаба Монина — Обухова;
P — вспомогательный коэффициент для вычисления масштаба Монина — Обухова;
pн — давление насыщенного пара опасного вещества при температуре воздуха, мм рт. ст.;
pн(T4) — давление насыщенного пара опасного вещества при температуре T4, Па;
qкомп — производительность компрессора, кг/с.
Учет производительности компрессора необходим при расчете выброса с трубопровода, на входе которого стоит компрессор, в случае, если площадь аварийного отверстия трубопровода превосходит 20% площади поперечного сечения трубопровода;
qнасос — производительность насоса, кг/с.
Учет производительности насоса необходим при расчете выброса с трубопровода, на входе которого стоит насос, в случае, если площадь аварийного отверстия трубопровода превышает 20% площади поперечного сечения трубопровода;
q_{i}^{l\text{разб}} — расход воздуха, подмешенного на месте выброса во вторичное облако, образующееся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l q_{i}^{l\text{разб}} = q_{i}^{\text{ж разб}} , q_{i}^{\text{отс}.\text{выб разб}} , q_{i}^{\text{г разб}} , q_{i}^{\text{ги разб}} , q_{i}^{\text{и разб}} , q_{i}^{\text{е разб}} ), кг/с;
q_{4}\prime — скорость поступления в атмосферу газообразного опасного вещества, образующегося при мгновенном вскипании жидкой фазы в случае истечения жидкого опасного вещества из разрушенного оборудования в сценарии 4, кг/с;
q_{4}\prime\prime — скорость поступления в атмосферу капель (аэрозоля) опасного вещества, образующихся при мгновенном вскипании жидкой фазы в случае истечения жидкого опасного вещества из разрушенного оборудования в сценарии 4, кг/с;
q_{4}\prime_{\text{отс}} — скорость поступления в атмосферу газообразного опасного вещества, образующегося при мгновенном вскипании жидкой фазы в случае истечения жидкого опасного вещества из разрушенного трубопровода после отсечения аварийного участка в сценарии 4, кг/с; принимается равной 0 кг/с, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
q_{4}\prime\prime_{\text{отс}} — скорость поступления в атмосферу капель (аэрозоля) опасного вещества, образующихся при мгновенном вскипании жидкой фазы в случае истечения жидкого опасного вещества из разрушенного трубопровода после отсечения аварийного участка в сценарии 4, кг/с; принимается равной 0 кг/с, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
q_{i}^{l} — расход опасного вещества (газообразного и находящегося в жидкой фазе) во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l q_{i}^{l} = q_{i}^{\text{ж}} , q_{i}^{\text{отс}.\text{выб}} , q_{i}^{\text{г}} , q_{i}^{\text{ги}} , q_{i}^{\text{и}} , q_{i}^{\text{е}} ), кг/с; q_{i}^{\text{отс}.\text{выб}} принимается равным 0 кг/с, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общего количества выброса;
qвыб — скорость выброса жидкого опасного вещества при истечении жидкого опасного вещества из разрушенного оборудования в сценарии 4, кг/с;
q_{\text{выб}}^{\text{отс}} — скорость выброса опасного вещества при истечении жидкого опасного вещества из разрушенного оборудования после отсечения аварийного участка трубопровода в сценарии 4, кг/с;
q_{\text{ж}i}^{l} — суммарный расход опасного вещества, находящегося в жидкой фазе (каплях), во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l q_{\text{ж}i}^{l} = q_{\text{ж}i}^{\text{ж}} , q_{\text{ж}i}^{\text{отс}.\text{выб}} , q_{\text{ж}i}^{\text{г}} , q_{\text{ж}i}^{\text{ги}} , q_{\text{ж}i}^{\text{и}} , q_{\text{ж}i}^{\text{е}} ), кг/с;
q_{\text{ж}i}^{0l} — суммарный расход опасного вещества, находящегося в жидкой фазе (каплях), на месте эмиссии вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l q_{\text{ж}i}^{0l} = q_{\text{ж}i}^{0\text{ж}} , q_{\text{ж}i}^{0\text{отс}.\text{выб}} , q_{\text{ж}i}^{0\text{г}} , q_{\text{ж}i}^{0\text{ги}} , q_{\text{ж}i}^{0\text{и}} , q_{\text{ж}i}^{0\text{е}} ), кг/с;
q_{\text{сум}i}^{l} — эффективный суммарный расход воздуха и опасного вещества (включая жидкую и газообразную фазы) во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l q_{\text{сум}i}^{l} = q_{\text{сум}i}^{\text{ж}} , q_{\text{сум}i}^{\text{отс}.\text{выб}} , q_{\text{сум}i}^{\text{г}} , q_{\text{сум}i}^{\text{ги}} , q_{\text{сум}i}^{\text{и}} , q_{\text{сум}i}^{\text{е}} ), кг/с;
k — константа Кармана, равна 0,41;
ri — радиус центральной части (ядра) первичного облака в i-ом сценарии, м;
t' — составляющая времени формирования первичного облака от начала выброса до времени отсечения аварийного участка трубопровода в сценарии 4, с;
t0 — момент времени, для которого определяются размеры зон достижения 0,5 НКПР и ВКПР, с;
tav — время осреднения при расчете дисперсии вдоль оси y(\sigma_{y}) , с; принимается в каждой точке равным времени прихода в эту точку облака, но не менее 600 с;
t_{i}^{l} — длительность поступления опасного вещества в атмосферу за счет истечения (испарения) (время формирования соответствующего вторичного облака) на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l t_{i}^{l} = t_{i}^{\text{ж}} , t_{i}^{\text{отс}.\text{выб}} , t_{i}^{\text{г}} , t_{i}^{\text{ги}} , t_{i}^{\text{и}} , t_{i}^{\text{е}} ), с; t_{i}^{\text{отс}.\text{выб}} принимается равным 0 с, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общего количества выброса;
tmax i — время достижения максимальной концентрации опасного вещества в некоторой точке в i-ом сценарии, с;
t_{max}^{\text{г}} — максимально возможная длительность истечения газообразного опасного вещества из разрушенного оборудования в сценарии 4, с;
t_{max}^{\text{ги}} — максимально возможная длительность истечения газообразного опасного вещества из разрушенного оборудования в сценарии 4 после испарения пролива, с;
tисп — длительность испарения пролива после окончания истечения жидкого опасного вещества для случая, если продолжается истечение газовой фазы, с;
t_{\text{и}}\prime_{\text{сп}} — вспомогательная величина при расчете t_{\text{о}}\prime_{\text{тс}} , с;
t_{\text{и}}\prime\prime_{\text{сп}} — вспомогательная величина при расчете t_{\text{о}}\prime_{\text{тс}} , с;
tкип — время, в течение которого опасное вещество поступает в первичное облако из-за интенсивного кипения жидкого ОВ в проливе за счет теплопритока от подстилающей поверхности, с;
tликвид — время ликвидации аварийного отверстия (разгерметизации) и пролива, с; если ликвидации не происходит, принимается равным бесконечности;
tотс — время отсечения аварийного участка трубопровода (время остановки компрессора (насоса), с;
t_{\text{о}}\prime_{\text{тс}} — составляющая времени формирования первичного облака, с;
u* — динамическая скорость, м/с;
u_{0\text{эфф}i}^{0} — эффективная скорость движения первичного облака в i-ом сценарии в начальный момент времени (на месте выброса), м/с;
u_{0\text{эфф}i}^{l} — эффективная скорость движения вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l u_{0\text{эфф}i}^{l} = u_{0\text{эфф}i}^{\text{ж}} , u_{0\text{эфф}i}^{\text{отс}.\text{выб}} , u_{0\text{эфф}i}^{\text{г}} , u_{0\text{эфф}i}^{\text{ги}} , u_{0\text{эфф}i}^{\text{и}} , u_{0\text{эфф}i}^{\text{е}} ) на месте выброса, м/с;
u10 — скорость ветра на высоте z10, м/с;
ut — вспомогательная величина при расчете числа Ричардсона, м/с;
u_{\text{подм}}^{\text{верх}} — скорость подмешивания воздуха в облако за счет диффузии в вертикальном направлении, м/с;
uэффi — эффективная скорость движения первичного облака в i-ом сценарии, м/с;
u_{\text{эфф}i}^{l} — эффективная скорость движения вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l u_{\text{эфф}i}^{l} = u_{\text{эфф}i}^{\text{ж}} , u_{\text{эфф}i}^{\text{отс}.\text{выб}} , u_{\text{эфф}i}^{\text{г}} , u_{\text{эфф}i}^{\text{ги}} , u_{\text{эфф}i}^{\text{и}} , u_{\text{эфф}i}^{\text{е}} ), м/с;
w* — вспомогательная величина при расчете числа Ричардсона, м/с;
x — пространственная переменная (координата вдоль ветра), м;
x0 — расстояние от места выброса до плоскости, перпендикулярной направлению ветра, для которой определяются размеры (по y и по z) зон достижения 0,5 НКПР и ВКПР, м;
x_{\text{з}i}^{l} — координата задней кромки вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l x_{\text{з}i}^{l} = x_{\text{з}i}^{\text{ж}} , x_{\text{з}i}^{\text{отс}.\text{выб}} , x_{\text{з}i}^{\text{г}} , x_{\text{з}i}^{\text{ги}} , x_{\text{з}i}^{\text{и}} , x_{\text{з}i}^{\text{е}} ), м;
x_{\text{п}i}^{l} — координата передней кромки вторичного облака, образующегося на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l x_{\text{п}i}^{l} = x_{\text{п}i}^{\text{ж}} , x_{\text{п}i}^{\text{отс}.\text{выб}} , x_{\text{п}i}^{\text{г}} , x_{\text{п}i}^{\text{ги}} , x_{\text{п}i}^{\text{и}} , x_{\text{п}i}^{\text{е}} ), м;
xцi — координата центра первичного облака в i-ом сценарии, м;
y — пространственная переменная (координата, перпендикулярная направлению ветра), м;
y0 — расстояние от оси выброса до перпендикулярной поверхности земли плоскости, для которой определяются размеры (по x и по z) зон достижения 0,5 НКПР и ВКПР; для плоскости, проходящей через ось выброса, y0 = 0, м;
z — пространственная переменная (координата высоты), м;
z0 — высота над поверхностью земли плоскости, для которой определяются размеры (по x и по y) зон достижения 0,5 НКПР и ВКПР, для поверхности земли z0 = 0, м;
z10 — стандартная высота, на которой задается скорость ветра, принимается равной 10 м;
zпов — характерный размер шероховатости поверхности, м;
\sum {}_{\text{вкпр}}(x,y,z) — соотношение, описывающее поверхность в пространстве с концентрацией ВКПР;
\sum {}_{0,5\text{нкпр}}(x,y,z) — соотношение, описывающее поверхность в пространстве с концентрацией 0,5 НКПР;
Ф — безразмерная вспомогательная величина при расчете скорости подмешивания воздуха в облако;
\alpha — объемная доля газовой фазы в оборудовании;
\alpha_{\text{в}} — показатель степенной зависимости скорости ветра от высоты;
\alpha_{i\text{возд}} — количество массовых долей воздуха (на единичную массовую долю ОВ), подмешивающегося в первичное облако на начальном этапе в i-ом сценарии;
\alpha_{i\text{возд}}^{l} — количество массовых долей воздуха (на единичную массовую долю ОВ), подмешивающегося во вторичное облако на начальном этапе на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии;
\beta — вспомогательный коэффициент, \beta = 1 + \alpha_{\text{в}} ;
\gamma — показатель адиабаты опасного вещества в газообразной фазе;
\gamma_{y} — вспомогательная величина при расчете \sigma_{y} ;
\gamma_{\text{подм}} — коэффициент пропорциональности при расчете воздуха в облаке при подмешивании через боковую поверхность; равен 0,63;
\delta — вспомогательная величина при расчете \sigma_{y} ;
\delta_{600} — вспомогательная величина при расчете \sigma_{y} ;
\varphi — вспомогательная величина при расчете динамической скорости;
\lambda_{\text{п}} — коэффициент теплопроводности подстилающей поверхности, на которую проливается опасное вещество, Вт/К/м;
\mu — молярная масса опасного вещества, кг/моль;
\mu_{\text{возд}} — молярная масса воздуха, кг/моль;
\mu_{\text{эфф}i} — эффективная молярная масса газообразной смеси опасного вещества с воздухом в первичном облаке в i-ом сценарии, кг/моль;
\mu_{\text{эфф}i}^{l} — эффективная молярная масса газообразной смеси опасного вещества с воздухом во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l \mu_{\text{эфф}i}^{l} = \mu_{\text{эфф}i}^{\text{ж}} , \mu_{\text{эфф}i}^{\text{отс}.\text{выб}} , \mu_{\text{эфф}i}^{\text{г}} , \mu_{\text{эфф}i}^{\text{ги}} , \mu_{\text{эфф}i}^{\text{и}} , \mu_{\text{эфф}i}^{\text{е}} ), кг/моль;
\pi — число, равное 3,14159..., отношение длины окружности к диаметру;
\rho(p_{H}(T_{4}),T_{4}) — плотность газообразного опасного вещества при температуре T4 и давлении pн(T4), кг/м3;
\rho_{i} — плотность газовой фазы опасного вещества в оборудовании в i-м сценарии, кг/м3;
\rho_{\text{возд}} — плотность воздуха в окружающей среде, кг/м3;
\rho_{\text{ж}} — плотность жидкого опасного вещества, кг/м3;
\rho_{\text{кип}} — плотность газообразного опасного вещества при температуре кипения и давлении P0, кг/м3;
\rho_{\text{п}} — плотность материала подстилающей поверхности, на которую проливается опасное вещество, кг/м3;
\rho_{\text{эфф}i} — эффективная плотность среды в первичном облаке в i-ом сценарии, кг/м3;
\rho_{\text{эфф}i}^{l} — эффективная плотность среды во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l \rho_{\text{эфф}i}^{l} = \rho_{\text{эфф}i}^{\text{ж}} , \rho_{\text{эфф}i}^{\text{отс}.\text{выб}} , \rho_{\text{эфф}i}^{\text{г}} , \rho_{\text{эфф}i}^{\text{ги}} , \rho_{\text{эфф}i}^{\text{и}} , \rho_{\text{эфф}i}^{\text{е}} ), кг/м3;
\rho_{i}^{\text{выб}} — плотность опасного вещества в первичном облаке в начальный момент (на месте выброса) времени в i-м сценарии, кг/м3;
\rho_{i}^{l} — плотность опасного вещества в начальный момент времени во вторичном облаке, образующемся на l-ой стадии поступления опасного вещества в атмосферу в i-ом сценарии (при различных l \rho_{i}^{l} = \rho_{i}^{\text{ж}} , \rho_{i}^{\text{отс}.\text{выб}} , \rho_{i}^{\text{г}} , \rho_{i}^{\text{ги}} , \rho_{i}^{\text{и}} , \rho_{i}^{\text{е}} ), кг/м3; \rho_{i}^{\text{отс}.\text{выб}} принимается равной 0 кг/м3, если количество опасного вещества в отсекаемом участке аварийного трубопровода составляет менее 20% общей массы выброса;
\sigma_{y} — дисперсия вдоль оси y (в поперечном направлении), м.
г — стадия аварии, на которой происходит истечение газообразного опасного вещества из оборудования при наличии пролива;
ги — стадия аварии, на которой происходит истечение газообразного опасного вещества из оборудования в отсутствие пролива;
е — стадия аварии, на которой происходит испарение опасного вещества, оставшегося в оборудовании (пролив отсутствует);
ж — стадия аварии, на которой происходит истечение жидкого опасного вещества из оборудования;
и — стадия аварии, на которой происходит испарение опасного вещества из пролива;
отс. выб — стадия аварии, на которой происходит истечение жидкого опасного вещества из разрушенного трубопровода после отсечения аварийного участка.
| | — модуль величины, равен самой величине, если величина больше нуля, и величине со знаком минус, если величина меньше нуля; например, 5 = |-5|;
sign — знак величины, равен 1, если величина больше нуля, минус 1, если величина меньше нуля, и 0, если величина 0; например, 1 = sign(5);
min — минимальное из всех значений, например, 3 = min{3, 5};
minsum(\alpha_{1},\alpha_{2},\alpha_{3},...,\alpha_{n}) — разность величины \alpha_{1} и суммы \alpha_{2} , \alpha_{3} , ..., \alpha_{n} , если эта разность величина положительная, и 0 в противном случае;
minsum(\alpha_{1},\alpha_{2},\alpha_{3},...,\alpha_{n}) = \frac{\alpha_{1} - \sum_{k = 2}^{n} \alpha_{i} + \alpha_{1} - \sum_{k = 2}^{n} \alpha_{i}}{2} ||;
minsum(2,1,3,4) = \frac{2 - (1 + 3 + 4) + |2 - (1 + 3 + 4)|}{2} = \frac{- 6 + 6}{2} = 0
minsum(9,1,3,4) = \frac{9 - (1 + 3 + 4) + |9 - (1 + 3 + 4)|}{2} = \frac{1 + 1}{2} = 1 .
\text{Г}(a) = \int_{0}^{+ \infty} x^{a - 1}e^{x}dx — гамма-функция (Г-функции, интеграл Эйлера второго рода) при a > 0.
Основные соотношения для Г-функции:
exp(x) — экспонента действительного числа x (число e = 2,71... в степени x, ex), например exp(-0,3) = 0,74081822068..., exp(1,3) = 3,6692966676...;
ln(x) — натуральный логарифм действительного числа x; например ln(0,740818220) = 0,3000..., exp(3,669297) = 1,3...;
arctg(x) — арктангенс действительного числа x (в радианах), например arctg(1) = 0,785398...
Приложение № 2 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
В настоящем документе используются следующие термины с соответствующими определениями:
авария — разрушение сооружений и (или) технических устройств, применяемых на опасном производственном объекте, неконтролируемый взрыв и (или) выброс опасных веществ;
время экспозиции — время, за которое набирается ингаляционная токсодоза (верхний предел интегрирования концентрации опасного вещества по времени в формуле расчета токсодозы);
вторичное облако (или шлейф) — облако опасного вещества, образующееся в результате длительного выброса газа или перегретой вскипающей жидкости, а также в результате испарения опасного вещества с подстилающей поверхности или из разгерметизированного оборудования и распространяющееся по ветру от места выброса; во вторичном облаке может существовать ядро — область пространства, в которой концентрация на заданной высоте постоянна (не изменяется при перемещении в горизонтальном направлении, перпендикулярном ветру, хотя может изменяться при перемещении по вертикали);
ингаляционная токсодоза — интеграл по времени концентрации опасного вещества в воздухе; при условно постоянной во времени концентрации опасного вещества в заданной точке — произведение концентрации опасного вещества в воздухе на время экспозиции;
класс устойчивости атмосферы — совокупность метеорологических факторов (скорость ветра, облачность, вертикальный тепловой поток), создающих определенные условия рассеяния в атмосфере. В зависимости от состояния атмосферы выделяют шесть классов: A, B, C, D, E, F;
концентрационные пределы распространения пламени — нижний (верхний) концентрационный предел распространения пламени (НКПР (ВКПР) — минимальное (максимальное) содержание горючего вещества в однородной смеси с окислительной средой, при котором возможно распространение пламени по смеси на любое расстояние от источника зажигания);
облако опасного вещества — область пространства, ограниченная поверхностями заданной концентрации опасного вещества;
обоснование безопасности — документ, содержащий сведения о результатах оценки риска аварии на опасном производственном объекте и связанной с ней угрозы, условия безопасной эксплуатации опасного производственного объекта, требования к эксплуатации, капитальному ремонту, консервации и ликвидации опасного производственного объекта;
опасный производственный объект — предприятие или его цеха, участки, площадки, а также иные производственные объекты, указанные в приложении 1 к Федеральному закону от 21 июля 1997 г. № 116-ФЗ «О промышленной безопасности опасных производственных объектов»;
опасные вещества — воспламеняющиеся, окисляющие, горючие, взрывчатые, токсичные, высокотоксичные вещества и вещества, представляющие опасность для окружающей природной среды, указанные в приложении 1 к Федеральному закону от 21 июля 1997 г. № 116-ФЗ «О промышленной безопасности опасных производственных объектов»;
отсекаемый участок аварийного трубопровода — участок трубопровода, на котором происходит выброс опасного вещества, между двумя ближайшими к месту выброса задвижками (насосами или компрессорами), которые при срабатывании (отключении) блокируют аварийный участок трубопровода;
отсечение аварийного участка трубопровода (отсечение аварийного потока) — остановка насоса (компрессора), стоящего на входе трубопровода, и (или) срабатывание запорной арматуры (задвижек), установленной на трубопроводе; прекращение подачи опасного вещества к месту выброса опасного вещества из трубопровода. При отсутствии трубопровода, но наличии истечения из емкости термин «отсечение аварийного участка» относится к прекращению истечения из отверстия, то есть к закупориванию отверстия;
первичное облако — облако опасного вещества, образующееся в результате очень быстрого (за 1 — 2 минуты) перехода в атмосферу части опасного вещества и распространяющееся по ветру от места выброса, в первичном облаке может существовать ядро — область пространства, в которой концентрация на заданной высоте постоянна;
пороговая токсодоза — наименьшая ингаляционная токсодоза опасного вещества, вызывающая у человека, не оснащенного средствами защиты органов дыхания, начальные признаки поражения организма с определенной вероятностью (табулированное значение для каждого опасного вещества);
разгерметизация оборудования — образование в оборудовании отверстий с размером, существенно меньшим, чем размеры оборудования, через которые опасное вещество в жидком или газообразном состоянии в течение некоторого времени поступает в окружающую среду;
разрушение оборудования — существенное нарушение целостности оборудования с образованием отверстий с размером, сопоставимым с размерами оборудования, при этом содержащееся в оборудовании опасное вещество в жидком или газообразном состоянии мгновенно выбрасывается в окружающую среду;
риск аварии — мера опасности, характеризующая возможность возникновения аварии на опасном производственном объекте и тяжесть ее последствий;
смертельная (летальная) токсодоза — наименьшая ингаляционная токсодоза опасного вещества, вызывающая у человека, не оснащенного средствами защиты органов дыхания, смерть с 50% вероятностью (табулированное значение для каждого опасного вещества);
«тяжелый» газ — смесь воздуха, газовой фазы и капель опасного вещества, плотность которой выше плотности окружающего воздуха.
Приложение № 3 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
СХЕМЫ РАСПРОСТРАНЕНИЯ ПЕРВИЧНОГО И ВТОРИЧНОГО ОБЛАКОВ (С УКАЗАНИЕМ ХАРАКТЕРИСТИК ОБЛАКОВ)
контуры изоконцентрации
контуры достижения заданной концентрации
Рисунок 3-1. Схема распространения первичного облака
контуры изоконцентрацииРисунок 3-2. Схема распространения вторичного облака
Приложение № 4 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
ВОЗМОЖНЫЕ КОНФИГУРАЦИИ ОБОРУДОВАНИЯ И СХЕМЫ ЕГО РАЗРУШЕНИЯ
- Для емкости с газом или жидкой фазой возможно мгновенное разрушение (разрыв) емкости и выброс всего содержимого в окружающую среду (сценарий 1 или сценарий 3 в соответствии с рисунками 4-1, 4-2).
Рисунок 4-1. Мгновенное разрушение емкости с газом (сценарий 1)
Рисунок 4-2. Мгновенное разрушение емкости с жидкой фазой (сценарий 3)
- Для емкости с газом или жидкой фазой возможны образование аварийного отверстия (разгерметизации) в стенке емкости и последующее истечение газа и (или) жидкости в окружающую среду (сценарий 2 или сценарий 4). При разгерметизации емкости с жидкой фазой возможно образование отверстия как выше, так и ниже уровня жидкости (в соответствии с рисунками 4-3 - 4-5).
Рисунок 4-3. Разгерметизация емкости с газом (сценарий 2)
Рисунок 4-4. Разгерметизация емкости с жидкой фазой выше уровня жидкости (сценарий 4)
Рисунок 4-5. Разгерметизация емкости с жидкой фазой ниже уровня жидкости (сценарий 4)
- Для емкости с газом или жидкой фазой с присоединенным трубопроводом возможно образование аварийного отверстия (разгерметизации) в стенке трубопровода либо полный разрыв трубопровода на некотором расстоянии от емкости (сценарий 2 или сценарий 4). При этом трубопровод может быть оснащен запорной арматурой, которая при срабатывании изолирует разгерметизированный (разрушенный) участок трубопровода от емкости. В этом случае в окружающую среду поступают газ и (или) жидкость и (или) двухфазный поток. На рисунках 4-6 - 4-8 отсекаемый участок аварийного трубопровода обозначен «1» и расположен справа от задвижки, которая, в свою очередь, расположена на трубопроводе около емкости.
Рисунок 4-6. Разгерметизация трубопровода с газом, присоединенного к емкости (1 — отсекаемый участок аварийного трубопровода) (сценарий 2)
Рисунок 4-7. Разгерметизация трубопровода с жидкой фазой, присоединенного к емкости (1 — отсекаемый участок аварийного трубопровода) (сценарий 4)
Рисунок 4-8. Разгерметизация трубопровода с газовой фазой, присоединенного к емкости (1 — отсекаемый участок аварийного трубопровода) (сценарий 4)
- Для трубопровода с газом или жидкой фазой с нагнетающим компрессором или насосом возможно образование аварийного отверстия (разгерметизации) в стенке трубопровода либо полный разрыв трубопровода на некотором расстоянии от компрессора (насоса) (сценарий 2 или сценарий 4). Возможны также разрушения (разгерметизации) самих компрессоров, насосов, в том числе по причине выдавливания сальников у электронасосов, неисправности запорной арматуры (вентилей) и т.д. При этом трубопровод может быть оснащен запорной арматурой, которая при срабатывании изолирует разгерметизированный (разрушенный) участок трубопровода от компрессора (насоса). В этом случае в окружающую среду поступают газ и (или) жидкость и (или) двухфазный поток. На рисунках 4-9, 4-10 отсекаемый участок аварийного трубопровода обозначен «1» и расположен справа от задвижки, которая, в свою очередь, расположена на трубопроводе около компрессора (насоса).
Рисунок 4-9. Разгерметизация трубопровода с газом, присоединенного к компрессору (1 — отсекаемый участок аварийного трубопровода) (сценарий 2)
Рисунок 4-9. Разгерметизация трубопровода с жидкой фазой, присоединенного к насосу (1 — отсекаемый участок аварийного трубопровода) (сценарий 4)
Приложение № 5 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
ВОЗМОЖНЫЕ СТАДИИ РАЗВИТИЯ АВАРИЙНОЙ СИТУАЦИИ
В качестве возможных стадий развития аварийной ситуации рекомендуется рассматривать:
разрушение оборудования и образование первичного облака;
истечение жидкой фазы до отсечения аварийного участка;
истечение жидкой фазы из аварийного участка после его отсечения;
истечение газа при наличии пролива жидкой фазы и испарение с пролива;
истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы;
испарение с пролива при отсутствии истечения жидкости или газа из разрушенного оборудования;
испарение из емкости при отсутствии пролива;
завершение аварии (ликвидация аварийного отверстия (разгерметизации) и пролива).
Каждая из вышеперечисленных стадий вносит свой вклад в суммарную массу выброса ОВ.
На каждой стадии аварии формируются свои облака ОВ в атмосфере (первичное и вторичные).
В зависимости от сценария, конфигурации оборудования, характера разрушения, свойств ОВ и действий по ликвидации аварии отдельные стадии из приведенных выше могут либо присутствовать, либо отсутствовать в той или иной аварийной ситуации.
Предполагается, что на каждой стадии процесс протекает стационарно.
Для более точных расчетов допускается разделение перечисленных стадий на отдельные подстадии меньшей продолжительности. Для каждой подстадии в этом случае определяются входные данные, указанные в Приложении № 6 к настоящему Руководству, с учетом изменений в системе, происшедших на предыдущих подстадиях.
Разрушение оборудования с выбросом всего объема ОВ, образование первичного облака, рассеяние первичного облака и воздействие на окружающую среду (в соответствии с рисунком 4-1).
Разрушение оборудования и истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы; рассеяние облака и воздействие на окружающую среду (в соответствии с рисунками 4-3, 4-6, 4-9).
При истечении газа из разрушенного трубопровода возможно отсечение аварийного участка трубопровода (либо в результате использования запорной арматуры, либо в результате остановки компрессоров, подающих ОВ в трубопровод, либо в результате и того, и другого) и истечение ОВ из него.
Возможно прекращение выброса путем ликвидации аварийного отверстия (разгерметизации).
Разрушение оборудования с жидким опасным веществом, выброс ОВ в окружающую среду, при наличии перегрева у жидкой фазы возможно ее вскипание с образованием в атмосфере газокапельного облака (в соответствии с рисунком 4-2). Часть жидкой фазы может пролиться на подстилающую поверхность — либо в обвалование, либо на неограниченную площадь. Если температура кипения жидкости при этом меньше температуры поверхности, то произойдет вскипание жидкости при ее соприкосновении с подстилающей поверхностью. Из газовой фазы, содержавшейся в оборудовании, из образовавшейся при вскипании за счет перегрева жидкой фазы газокапельной фазы и из газа, образующегося при кипении пролива, образуется первичное облако, которое рассеивается в атмосфере, воздействуя на окружающую среду.
Из пролива происходит испарение ОВ, в результате чего образуется вторичное облако, которое также рассеивается в атмосфере, воздействуя на окружающую среду.
Возможно прекращение поступления ОВ в окружающую среду путем ликвидации пролива.
Разрушение оболочки емкости выше уровня жидкости и длительное истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы (если жидкость находится в перегретом состоянии, то происходит вскипание жидкости, в результате которого в дополнение к газовой фазе, содержащейся в емкости на момент начала аварии, добавится газовая фаза, образовавшаяся при кипении), рассеяние газового облака ОВ (вторичного) и воздействие его на окружающую среду (в соответствии с рисунками 4-4, 4-5). После спада давления в емкости до атмосферного поступление ОВ в окружающую среду будет обусловлено лишь испарением ОВ с поверхности жидкости в емкости. При этом в атмосфере образуется вторичное облако, состоящее из газообразного ОВ, поступающего из разгерметизированной емкости за счет испарения ОВ из нее. Вторичное облако будет формироваться на месте аварии до тех пор, пока не испарится все ОВ из емкости. Возможно прекращение выброса путем ликвидации аварийного отверстия (разгерметизации).
Разрушение оболочки емкости ниже уровня жидкости и истечение жидкой фазы из разрушенного оборудования, образование пролива на месте выброса.
Если жидкость в емкости находилась в перегретом состоянии, то происходит вскипание жидкости сразу после ее выброса из емкости и образование в атмосфере газокапельной взвеси. Затем, если температура кипения меньше температуры поверхности, происходит кипение жидкой фазы (той ее части, что не участвовала в формировании газокапельной взвеси в атмосфере) при проливе ее на подстилающую поверхность. При этом из газа, образовавшегося при кипении пролива, а также из газокапельной взвеси выброса, поступившего из емкости в атмосферу за время кипения пролива, формируется первичное облако, которое рассеивается в атмосфере и воздействует на окружающую среду.
При истечении перегретой жидкости выброс жидкой фазы вскипает до выпадения на землю, образуя в атмосфере газокапельную взвесь ОВ. При этом в атмосфере образуется вторичное облако, состоящее из газообразного ОВ, испарившегося с пролива, и (при выбросе перегретой жидкости) из газокапельной взвеси, образующейся при вскипании ОВ сразу после выброса. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать возможность выброса жидкой фазы, то есть до момента выброса из оборудования всей жидкой фазы, находившейся выше уровня разгерметизации. Причем если в емкости находилась перегретая жидкость и если давление могло опуститься ниже давления насыщенного пара, то по мере вытекания жидкости возможно вскипание перегретого ОВ в самой емкости.
После окончания выброса жидкой фазы из емкости начинается истечение газовой фазы за счет наличия избыточного давления в емкости. При этом на стадии истечения газовой фазы продолжается испарение пролива. В результате в атмосфере образуется вторичное облако ОВ, состоящее из газообразного ОВ, испарившегося с пролива, и из газообразного ОВ, поступающего из разрушенной емкости за счет наличия избыточного давления в емкости. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать возможность выброса газовой фазы из емкости за счет избыточного давления в емкости и испарения из пролива. Причем если в емкости находилась перегретая жидкость, то в формирующееся на месте аварии вторичное облако будет поступать не только газовая фаза, находившаяся в емкости на начало аварии, но и газовая фаза, образующаяся при вскипании в емкости с перегретой жидкостью.
Если давление в емкости становится равным атмосферному, то газовая фаза из разгерметизированной емкости перестает поступать в окружающую среду за счет избытка давления. Если пролив ОВ при этом еще не испарился, то в атмосфере образуется вторичное облако ОВ, состоящее из газообразного ОВ, испарившегося с пролива. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать испарение ОВ из пролива.
Если после испарения пролива давление в емкости еще не стало равным атмосферному, то газовая фаза из разгерметизированной емкости продолжает поступать в окружающую среду за счет избыточного давления. В атмосфере при этом образуется вторичное облако ОВ, состоящее из газообразного ОВ. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать избыточное давление в емкости.
Если пролив ОВ на месте аварии испарился и давление в емкости стало равным атмосферному, то поступление ОВ в окружающую среду будет обусловлено лишь испарением ОВ с поверхности жидкости в емкости. При этом в атмосфере образуется вторичное облако, состоящее из газообразного ОВ, поступающего из разгерметизированной емкости за счет испарения. Вторичное облако будет формироваться на месте аварии до тех пор, пока не испарится все ОВ из емкости.
Возможно прекращение поступления ОВ в окружающую среду путем ликвидации аварийного отверстия и пролива.
Сценарий 4. Аварии на трубопроводах, на входе которых стоит емкость
Разрушение трубопровода, выходящего из емкости выше уровня жидкости, и истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы (в соответствии с рисунками 4-7, 4-8). Если жидкость находится в перегретом состоянии, то при падении давления может происходить вскипание жидкости, в результате которого в дополнение к газовой фазе, содержащейся в емкости на момент начала аварии, добавится газовая фаза, образовавшаяся при кипении. При этом в атмосфере происходит рассеяние вторичного газового облака ОВ и воздействие его на окружающую среду. При наличии запорной арматуры возможно отсечение аварийного участка трубопровода и истечение ОВ только из этого аварийного участка. Если давление в емкости (или в отсеченном участке трубопровода) уменьшилось и стало равным атмосферному, то поступление ОВ в окружающую среду будет обусловлено лишь его испарением с поверхности жидкости в емкости и (или) в отсеченном участке трубопровода. При этом в атмосфере образуется вторичное облако, состоящее из газа, поступающего из разгерметизированной емкости и (или) отсеченного участка трубопровода за счет испарения. Вторичное облако будет формироваться на месте аварии до тех пор, пока не испарится все ОВ из емкости.
Возможно прекращение выброса путем ликвидации аварийного отверстия (разгерметизации).
При разрушении трубопровода, выходящего из емкости ниже уровня жидкости, стадии аварии будут аналогичны описанным в разделе «Аварии на емкости», с двумя отличиями.
Во-первых, при существенном разрушении трубопровода в случае выброса перегретой жидкости вскипание жидкости может происходить не только в воздухе, но и в самом трубопроводе.
Во-вторых, при наличии на трубопроводе запорной арматуры и срабатывании ее в ходе аварии появляется этап — истечение из отсеченного участка трубопровода, который полностью совпадает по физической картине с описанным выше выбросом из емкости - сценарий 4. Аварии на емкости.
Сценарий 4. Аварии на трубопроводах, на входе которых стоит насос
При разрушении трубопровода происходит истечение жидкой фазы и образование пролива на месте выброса (в соответствии с рисунком 4-10). Если жидкость в трубопроводе находилась в перегретом состоянии, то происходит вскипание жидкости либо сразу после ее выброса из трубопровода (случай, когда при сильном разрушении трубопровода место выброса находится вблизи насоса, и случаи небольших отверстий разгерметизации), либо внутри трубопровода (случай, когда при сильном разрушении трубопровода место выброса находится вдали от насоса).
При выбросе перегретой жидкости в атмосфере образуется газокапельная взвесь, а затем происходит кипение жидкой фазы (той ее части, которая не участвовала в формировании газокапельной взвеси в атмосфере) при проливе на подстилающую поверхность, если температура кипения меньше температуры подстилающей поверхности. При этом из газа, образовавшегося при кипении пролива, а также из газокапельной взвеси выброса, поступившего из трубопровода в атмосферу за время кипения пролива, формируется первичное облако, которое рассеивается в атмосфере и воздействует на окружающую среду.
Продолжающийся выброс жидкой фазы идет в пролив, с которого происходит испарение ОВ. При выбросе перегретой жидкости часть жидкой фазы вскипает до выпадения на землю, образуя в атмосфере газокапельную взвесь ОВ. При этом в атмосфере образуется вторичное облако ОВ, состоящее из газообразного ОВ, испарившегося с пролива, и (при выбросе перегретой жидкости) из газокапельной взвеси, образующейся при вскипании ОВ, либо сразу после выброса, либо в трубопроводе. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать возможность выброса жидкой фазы (двухфазной смеси), то есть до момента отсечения аварийного участка трубопровода и выброса из отсеченного участка трубопровода всей жидкой фазы, находящейся в нем выше уровня аварийного отверстия (разгерметизации).
После окончания выброса жидкой фазы из отсеченного участка трубопровода начинается истечение газовой фазы за счет наличия избыточного давления в нем, что возможно только в случае, если в трубопроводе находится перегретая жидкость. При этом на стадии истечения газовой фазы продолжается испарение пролива. В результате в атмосфере образуется вторичное облако, состоящее из газообразного ОВ, испарившегося с пролива, и из газообразного ОВ, поступающего из разрушенного трубопровода за счет наличия избыточного давления в нем. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать возможность выброса газовой фазы из трубопровода за счет избыточного давления в нем и испарения из пролива.
Если давление в отсеченном участке аварийного трубопровода становится равным атмосферному, то газовая фаза из разгерметизированного трубопровода перестает поступать в окружающую среду за счет избытка давления. Если пролив ОВ при этом еще не испарился, то в атмосфере образуется вторичное облако ОВ, состоящее из газообразного ОВ, испарившегося с пролива. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать испарение ОВ из пролива.
Если испарился пролив, а давление в отсеченном участке аварийного трубопровода осталось выше атмосферного, то газовая фаза из разгерметизированного трубопровода продолжает поступать в окружающую среду за счет избыточного давления. При этом в атмосфере образуется вторичное облако, состоящее из газообразного ОВ, поступающего из отсеченного участка аварийного трубопровода за счет избыточного давления. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока будет существовать избыточное давление в отсеченном участке аварийного трубопровода.
Если пролив ОВ на месте аварии испарился и давление в отсеченном участке аварийного трубопровода стало равным атмосферному, то поступление ОВ в окружающую среду будет обусловлено лишь испарением ОВ с поверхности жидкости в трубопроводе. В атмосфере при этом образуется вторичное облако, состоящее из газообразного ОВ, поступающего из разгерметизированного трубопровода за счет испарения ОВ. Такое вторичное облако будет формироваться на месте аварии до тех пор, пока не испарится все ОВ из трубопровода.
Возможно прекращение поступления ОВ в окружающую среду путем ликвидации аварийного отверстия и пролива в любой момент времени.
Приложение № 6 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
ПЕРЕЧЕНЬ ДАННЫХ, НЕОБХОДИМЫХ ДЛЯ ПРОВЕДЕНИЯ РАСЧЕТА ПО ОЦЕНКЕ ПОСЛЕДСТВИЙ АВАРИЙНЫХ ВЫБРОСОВ ОПАСНЫХ ВЕЩЕСТВ
Для проведения расчетов по любому из четырех сценариев, приведенных в разделе III настоящего Руководства, рекомендуется иметь следующие данные:
тип местности, где происходит рассеяние выброса (необходим для определения zпов — характерного размера шероховатости подстилающей поверхности (таблица № 7-3);
показатель степенной зависимости скорости ветра от высоты \alpha_{\text{в}} ;
скорость ветра на высоте 10 м;
температура воздуха, плотность воздуха и удельная внутренняя энергия воздуха, атмосферное давление (при нормальных условиях принимается равным 101,325 кПа);
класс устойчивости атмосферы или уровень инсоляции и облачности (для определения класса устойчивости атмосферы);
температура поверхности, над которой происходит рассеяние выброса;
значения токсодоз и концентраций ОВ в воздухе, соответствующие тому или иному воздействию ОВ на человека, способ расчета пробит-функции токсического поражения человека;
индивидуальные характеристики ОВ в газовой и жидкой фазах (молярная масса, показатель адиабаты, температура кипения, теплота испарения (кипения), теплоемкость газовой и жидкой фаз, плотность жидкой фазы, значения НКПР и ВКПР для пожаровзрывоопасных веществ);
при наличии пролива в обвалование — площадь обвалования, ее высота и площадь возможного контакта жидкости с твердой поверхностью;
при наличии пролива вскипающей жидкости — характеристики подстилающей поверхности (плотность, теплоемкость, теплопроводность) и температура подстилающей поверхности (на месте пролива);
время экспозиции для человека, попадающего в область поражения.
Дополнительно для каждого сценария задаются следующие данные:
Дополнительные данные для каждого сценария
| Сценарий 1 | Q, V1, T1, либо Q, P1, T1, либо Q, P1, V1, либо V1, P1, T1 |
|---|---|
| Сценарий 2 | Q, V2, T2, либо Q, P2, T2, либо Q, P2, V2, либо V2, P2, T2, Qтр, S, Sтр, qкомп, qликвид |
| Сценарий 3 | P3, T3, Qг либо V3 и \alpha |
| Сценарий 4 | P4, T4, \rho_{4} , H, L, Dтр, qнасос, S, tликвид, tотс, Sобор, Qж, Q_{H}^{\text{ж}} , Q_{\text{тр}}^{\text{ж}} , H_{\text{тр}}^{\text{отс}} , Q_{H^{\text{тр}}}^{\text{ж}} , Q_{\text{тр}^{\text{отс}}}^{\text{ж}} , Q_{\text{тр}}^{\text{ж отс}} , Q_{H}^{\text{ж отс}.\text{тр}} |
Приложение № 7 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
СПРАВОЧНЫЕ ДАННЫЕ ПО СВОЙСТВАМ ОПАСНЫХ ВЕЩЕСТВ, СВЯЗЯМ ВЕРОЯТНОСТИ ПОРАЖЕНИЯ С ПРОБИТ-ФУНКЦИЕЙ, ХАРАКТЕРИСТИКАМ ПОДСТИЛАЮЩИХ ПОВЕРХНОСТЕЙ И СОСТОЯНИЯ АТМОСФЕРЫ
Связь вероятности поражения с пробит-функцией
Характерный размер шероховатости поверхности zпов в зависимости от типа местности, где происходит рассеяние
| Характерный размер шероховатости поверхности zпов, м | Местность | ||
| Природный ландшафт | Антропогенный ландшафт | ||
| 1,00·10-5 | Лед, равнина, покрытая грязью | ||
| 9,00·10-5 | Равнины, покрытые снегом, или укатанный грунт | ||
| 1,00·10-4 | Поверхность открытого моря при штиле | ||
| 1,00·10-4 — 1,00·10-3 | Обширные водные поверхности | ||
| 5,00·10-4 | Ровная поверхность пустыни | ||
| 9,00·10-4 | Поверхность моря в прибрежной зоне при ветре с моря | ||
| 2,00·10-3 | Снежная целина (с/х угодья) | ||
| 7,50·10-3 | Равнинная местность: скошенная трава (~= 30 см) | Равнинная местность, сельскохозяйственные угодья | |
| 1,00·10-2 | Равнинная местность: трава, редкие деревья (зима, без листьев) | ||
| 2,20·10-2 | Равнинная местность: некошеная трава | ||
| 2,50·10-2 | Равнинная местность: одиночные деревья | Аэропорт — летное поле | |
| 5,00·10-2 | Равнинная местность: высокая трава (60 см) | Сельскохозяйственные угодья: неубранные посевы зерновых | |
| 5,50·10-2 | Равнинная местность: редкие деревья (лето) | ||
| 8,50·10-2 | Холмистая местность | ||
| 0,15 — 0,30 | Деревья | Большое количество заборов, изгородей, редкие здания | |
| 0,40 | Лесистая местность | Окраины города <*> | |
| 0,55 | Центры малых городов <*> | ||
| 0,65 | Центры больших поселений и городов <*> | ||
| 0,90 | Леса | ||
| 1,30 | |||
| 1,35 | Центры больших городов, районы с высотными зданиями <*> | ||
| 2,00 | Крайне холмистая и гористая местность <*> | ||
| 3,00 | |||
<*> Неоднородный рельеф местности.
Класс устойчивости атмосферы (модифицированные классы устойчивости по Паскуилу)
| Скорость ветра, м/с | День <1> | Сумерки <2> | Ночь | |||||
| Интенсивность солнечного излучения (Вт/м2) | Уровень облачности (восьмые доли) | |||||||
| сильная (> 600) | умеренная (300 — 600) | слабая (< 300) | сплошная облачность | 0 — 3 | 4 — 7 | 8 | ||
| <= 2 | A | A — B | B | C | D | F | F | D |
| 2 — 3 | A — B | B | C | C | D | F | E | D |
| 3 — 5 | B | B — C | C | C | D | E | D | D |
| 5 — 6 | C | C — D | D | D | D | D | D | D |
| >= 6 | C | D | D | D | D | D | D | D |
<1> Период времени за исключением 1 часа после восхода солнца и 1 часа до заката.
<2> Сумерки включают в себя 1 час после восхода солнца и 1 час до заката.
Величина коэффициента \alpha_{\text{в}} в зависимости от класса устойчивости атмосферы для различных характерных размеров шероховатости
Примечание. Приведены значения коэффициента для трех высот распространения выброса до 20 м, до 50 м, свыше 50 м.
Величины коэффициентов kL и p для вычисления масштаба Монина — Обухова
| Класс устойчивости | A | B | C | E | F |
|---|---|---|---|---|---|
| kL | -11,4 | -26,0 | -123,0 | 123,0 | 26,0 |
| p | 0,10 | 0,17 | 0,30 | 0,30 | 0,17 |
Значения коэффициента \delta_{600} при различной устойчивости атмосферы
| Класс устойчивости | A | B | C | D | E | F |
|---|---|---|---|---|---|---|
| \delta_{600} | 0,22 | 0,16 | 0,11 | 0,08 | 0,06 | 0,04 |
Значения характеристик подстилающих поверхностей
Приложение № 8 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО ОПРЕДЕЛЕНИЮ ПАРАМЕТРОВ СМЕСИ ОПАСНОГО ВЕЩЕСТВА С ВОЗДУХОМ В ОБЛАКЕ
При определении состояния смеси ОВ с воздухом предполагается, что в воздухе нет паров и капель воды. ОВ в облаке может находиться как в газовой, так и в жидкой фазах. Химические реакции между опасным веществом и компонентами воздуха, образование растворов газовой фазы одного вещества в жидкой фазе (в каплях) другого вещества не учитываются.
Предполагается, что объем, занимаемый жидкой фазой, незначителен, при расчетах он полагается равным нулю.
Определение состояния смеси ОВ с воздухом в облаке проводится для каждой рассматриваемой точки пространства.
Исходными данными для проведения расчета являются:
общая масса (расход) (в газообразной и жидкой фазах) ОВ в облаке (Qi — для первичного облака и q_{i}^{l} — для вторичных облаков);
масса (расход) жидкой фазы в облаке в начальный момент времени ( Q_{\text{ж}i}^{0} — для первичного облака и q_{\text{ж}i}^{0l} — для вторичных облаков);
общая масса (расход) облака (Qсумi — для первичного облака и q_{\text{сум}i}^{l} — для вторичных облаков);
полная внутренняя энергия (удельный поток внутренней энергии) смеси вещества в облаке (Eэффi — для первичного облака и e_{\text{эфф}i}^{l} — для вторичных облаков);
геометрические характеристики облаков (радиус и высота для первичного облака, полуширина и высота для вторичных).
Кроме того, необходимы данные об окружающей среде и физические характеристики воздуха и ОВ.
Перечисленные исходные данные получаются в результате решения уравнений (101) - (104), (107) - (113) для первичного облака или (118) - (121), (124) - (131) для вторичных облаков. В результате определения состояния смеси ОВ и воздуха в облаке получаются:
плотность смеси ( \rho_{\text{эфф}i} — для первичного облака или \rho_{\text{эфф}i}^{l} — для вторичных облаков);
температура смеси (Tэффi — для первичного облака или T_{\text{эфф}i}^{l} — для вторичных облаков);
масса (расход) ОВ в жидкой фазе (Qжi — для первичного облака или q_{\text{ж}i}^{l} — для вторичных облаков);
изменение геометрического размера облака — эффективная высота (Hэффi — для первичного облака или H_{\text{эфф}i}^{l} — для вторичных облаков).
Далее приводятся формулы для расчета состояния смеси в первичном облаке, для вторичных облаков все приведенные ниже формулы справедливы при замене соответствующих переменных, описывающих состояние в первичном облаке, на переменные, описывающие состояния во вторичных облаках, а именно: Qжi на q_{\text{ж}i}^{l} , Qi на q_{i}^{l} , Q_{\text{ж}i}^{0} на q_{\text{ж}i}^{0l} , Qсумi на q_{\text{сум}i}^{l} , Tэффi на T_{\text{эфф}i}^{l} , \mu_{\text{эфф}i} на \mu_{\text{эфф}i}^{l} , Vэффi на V_{\text{эфф}i}^{l} , \rho_{\text{эфф}i} на \rho_{\text{эфф}i}^{l} .
Состояние в облаке определяется из условия совпадения энергии системы и ее плотности, распределенной в объеме, с энергией облака и массой облака при давлении P0 = 101,325 кПа.
В случае наличия в начальный момент времени перегретой жидкой фазы (Q_{\text{ж}i}^{0}>0) и Tвозд > Tкип сначала определяется количество жидкой фазы ОВ в облаке:
. (145)Если величина Qжi > 0, то это означает, что в облаке находится капельная взвесь массой Qжi, тогда температура в облаке Tэффi полагается равной Tкип.
Плотность (удельный объем) в облаке определяется по соотношениям:
\mu_{\text{эфф}i} = \frac{(Q_{\text{сум}i} - Q_{\text{ж}i})\mu_{\text{возд}}\mu}{(Q_{i} - Q_{\text{ж}i})\mu_{\text{возд}} + (Q_{\text{сум}i} - Q_{i})\mu} , (146)
\rho_{\text{эфф}i} = \frac{P_{0}\mu_{\text{эфф}i}}{T_{\text{эфф}i}R} + c_{\text{ц}i} - \frac{P_{0}\mu}{T_{\text{эфф}i}R} , (147)
V_{\text{эфф}i} = \frac{1}{\rho_{\text{эфф}i}} . (148)
Если величина Qжi <= 0, то это означает, что в облаке находится только газовая смесь и массу капельной взвеси Qжi рекомендуется полагать равной 0.
Температура в облаке в этом случае, а также в случае изначального отсутствия жидкой фазы в облаке определяется из соотношения:
T_{\text{эфф}i} = \frac{(Q_{i} - Q_{\text{ж}}^{0})C_{v}T_{i}^{0} + Q_{\text{ж}}^{0}C_{p}T_{i}^{0} - \Delta H_{\text{кип}}Q_{\text{ж}}^{0} + (Q_{\text{сум}i} - Q_{i})C_{v\text{возд}}T_{\text{возд}}}{(Q_{\text{сум}i} - Q_{i})C_{v\text{возд}} + Q_{i}C_{v}} . (149)
Плотность (удельный объем) в облаке определяется по соотношениям:
\mu_{\text{эфф}i} = \frac{Q_{\text{сум}i}\mu\mu_{\text{возд}}}{Q_{i}\mu_{\text{возд}} + (Q_{\text{сум}i} - Q_{i})\mu} , (150)
\rho_{\text{эфф}i} = \frac{P_{0}\mu_{\text{эфф}i}}{T_{\text{эфф}i}R} , (151)
V_{\text{эфф}i} = \frac{1}{\rho_{\text{эфф}i}} . (152)
В случае наличия в начальный момент времени перегретой жидкой фазы (Q_{\text{ж}i}^{0}>0) и Tвозд < Tкип количество жидкой фазы в облаке предполагается постоянным и равным Q_{\text{ж}i}^{0} .
Температура в облаке в этом случае определяется из соотношения:
T_{\text{эфф}i} = \frac{(Q_{i} - Q_{\text{ж}i}^{0})C_{vi}T_{i}^{0} + Q_{\text{ж}i}^{0}C_{p}T_{\text{кип}} + (Q_{\text{сум}i} - Q_{i})C_{v\text{возд}}T_{\text{возд}}}{Q_{\text{ж}i}^{0}C_{p} + (Q_{i} - Q_{\text{ж}i}^{0})C_{vi} + (Q_{\text{сум}i} - Q_{i})C_{v\text{возд}}} . (153)
После определения состояния в облаке производится определение геометрических размеров, изменение которых произошло за счет испарения жидкой фазы, а также за счет нагрева (охлаждения) газовой фазы при смешении с воздухом.
Примеры расчетов последствий аварий с выбросом ОВ при рассматриваемых выше условиях приведены в Приложении № 9 к настоящему Руководству.
Приложение № 9 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
ПРИМЕРЫ РАСЧЕТОВ ПОСЛЕДСТВИЙ ВЫБРОСОВ ОПАСНЫХ ВЕЩЕСТВ
Расчеты проводятся в соответствии со схемой, приведенной в Приложении № 10 к настоящему Руководству.
Пример 1. В опорожненном шаровом резервуаре вместимостью 2000 м3 находится хлорметан в газовой фазе при температуре окружающей среды (18 °C) и давлении 1 атм. Резервуар расположен на товарно-сырьевой базе в равнинной местности, характерный размер шероховатости 0,018 м. Происходит разрушение емкости, и выбрасывается все опасное вещество. Метеоусловия: скорость ветра на высоте 10 м — 3,2 м/с, безоблачная ночь.
Определить максимальные приземные концентрации хлорметана в атмосфере, размеры взрывопожароопасного облака смеси хлорметана с воздухом при его рассеянии до концентрации, соответствующей нижнему концентрационному пределу воспламенения и его половине, а также зависимость площади, покрываемой облаком, от времени.
1. При решении используем следующие характеристики хлорметана:
температура кипения жидкой фазы при давлении окружающей среды (101325 Па) — -24 °C;
теплота испарения жидкой фазы — 3,87·106 Дж/кг;
теплоемкость жидкой фазы — Cp = 1585 Дж/кг/К;
плотность жидкой фазы — 930 кг/м3;
2. По условиям хранения и характеру разрушения определяем тип сценария аварийной ситуации - 1.
3. По формулам (1) и (2) определяем количество хлорметана в первичном облаке Q1 = 4272 кг.
4. По формуле (6) с использованием формулы (2) определяем плотность ОВ в выбросе — 2,13 кг/м3.
5. Учитывая, что хлорметан находится в сферической емкости, начальный диаметр первичного облака полагается равным его высоте. Они соответственно равны 13,6 м.
6. Характерный размер шероховатости принимается из исходных данных zпов = 0,018 м.
7. С использованием данных таблицы № 7-4 определяем класс устойчивости атмосферы — E.
8. С использованием данных таблиц № 7-5 - 7-7 определяем коэффициенты, используемые при расчете дисперсии, скорости подмешивания воздуха в облако.
9. Для определения концентрации хлорметана в атмосфере (по формулам (105) - (106)) численно решаем систему обыкновенных дифференциальных уравнений (формулы (107) - (115)) с использованием уравнений (101) - (104).
Результаты расчета приведены на рисунке 9-1.
Рисунок 9-1. Максимальная приземная концентрация на оси облака (по направлению ветра) при выбросе хлорметана (пример 1)
10. С помощью соотношений (139) - (143) определяем размеры взрывопожароопасного облака смеси хлорметана с воздухом при его рассеянии (рисунок 9-2).
Значения концентраций, кг/м3 (плоскость XY, Z = 0.0 м)
Рисунок 9-2. Сечение зоны, ограниченной НКПР и 0,5 НКПР плоскостью, z = 0 (в направлении по ветру) при выбросе хлорметана (пример 1)
11. Так как высота облака в пределах зоны, где рассматривается рассеяние (1000 м), не превосходит 20 м, то согласно пункту 23 настоящего Руководства повторный расчет не требуется.
12. Для определения площади облака воспользуемся соотношением S_{\text{эфф}} = \pi R_{\text{эфф}}^{2} . Зависимость площади, покрываемой облаком, от времени приведена на рисунке 9-3.
Рисунок 9-3. Изменение площади облака в зависимости от времени при выбросе хлорметана (пример 1)
Пример 2. Технологический трубопровод Ду200 с хлорцианом в газовой фазе расположен на эстакаде высотой 50 см. Трубопровод расположен на территории промышленной площадки в равнинной местности, характерный размер шероховатости 0,018 м. На входе трубопровода стоит компрессор, расход которого составляет 10,3 кг/с. Давление в трубопроводе 1,3 атм. Трубопровод не теплоизолирован. Ночью в нижней образующей трубопровода в результате разгерметизации возникает отверстие диаметром 10 см. Время ликвидации утечки составляет 400 сек. Метеоусловия: скорость ветра на высоте 10 м — 2,1 м/с, температура воздуха — 30 °C, облачность — 50%. Профиль ветра аппроксимируется степенной зависимостью с показателем — 0,22.
Определить концентрации хлорциана в атмосфере, а также распределение токсодозы в зоне химического заражения, где возможно смертельное поражение незащищенных людей на открытой местности.
1. По условиям хранения и характеру разрушения определяем тип сценария аварийной ситуации - 2.
2. С учетом комментария к формуле (14) определяем расход выброса qвыб = 10,3 кг/с.
3. Определяем плотность ОВ в выбросе по формуле (18) \rho_{2}^{\text{ги}} = 2,65\text{кг}/\text{м}^{3} .
4. По формуле (16) определяем длительность выброса t_{2}^{\text{ги}} = 400\text{с} .
5. Поскольку выброс осуществляется с высоты 0,5 м и плотность выброса больше плотности воздуха, высота вторичного облака в начальном сечении была принята H_{2}^{\text{ги}} = 0,5\text{м} . Тогда по формуле (20) рассчитывается ширина сечения вторичного облака (шлейфа) в начальный момент времени (на месте выброса) — B_{2}^{\text{ги}} = 4,6\text{м} .
6. С использованием данных таблицы № 7-4 определяем класс устойчивости атмосферы — E.
7. С использованием данных таблиц № 7-5 - 7-7 определяем коэффициенты, используемые при расчете дисперсии, скорости подмешивания воздуха в облако.
8. Для определения концентрации хлорциана в атмосфере (122) - (123) численно решаем систему обыкновенных дифференциальных уравнений (124) - (132) с использованием уравнений (118) - (121) для стадии истечения газа в отсутствии пролива.
9. По формулам (136) - (137) рассчитываем поле концентрации при выбросе газообразного хлорциана из трубопровода. Соответствующие значения концентрации на оси выброса приведены на рисунке 9-4.
Рисунок 9-4. Максимальная концентрация на оси облака (в направлении по ветру) при выбросе хлорциана (пример 2)
10. В соответствии с пунктом 40 настоящего Руководства и по значению смертельной токсодозы для хлорциана 11 мг·мин/л определяем размеры зоны смертельного поражения (рисунки 9-5, 9-6). Протяженность зоны смертельного поражения составляет 905 м. Максимальная высота — 4,45 м на расстоянии 382 м. Максимальная ширина 407 м на расстоянии 566 м.
Рисунок 9-5. Размеры зоны смертельного поражения в плоскости y = 0 (вертикальное сечение через ось выброса) при выбросе хлорциана (пример 2)
Рисунок 9-6. Размеры зоны смертельного поражения на поверхности земли при выбросе хлорциана (пример 2)
11. В соответствии с пунктом 40 настоящего Руководства и по значению пороговой токсодозы для хлорциана 0,75 мг·мин/л определяем размеры зоны порогового воздействия (рисунки 9-5, 9-6). Протяженность зоны порогового воздействия — 4077 м. Максимальная высота — 22,3 м на расстоянии 1677 м. Максимальная ширина — 874 м на расстоянии 1844 м.
Пример 3. На объекте хранения ОВ, расположенном в пределах небольшого города, имеется емкость объемом 100 м3. Она на 50% заполнена сжиженным аммиаком под давлением. Аммиак в емкости имеет температуру окружающей среды — 30 °C. Емкость стоит на бетонной поверхности.
Оценить, полагая наихудшие условия рассеяния, последствия аварии в случае мгновенного разрушения емкости. При оценке токсического поражения использовать время экспозиции — 30 мин.
1. По условиям хранения и характеру разрушения определяем тип сценария аварийной ситуации - 3. Согласно пункту 22 настоящего Руководства в качестве наихудших условий рассеяния принимается скорость ветра на высоте 10 м — 1,0 м/с, устойчивость атмосферы — класс F. Выбираем тип подстилающей поверхности, над которой происходит рассеяние, — центр малых городов, с характерным размером шероховатости 0,55 м.
2. Время экспозиции принимаем равным 30 мин.
3. По формуле (24) определяем массу газовой фазы в емкости — 393,16 кг.
4. По объему емкости (100 м3), степени заполнения (50%) и плотности жидкой фазы аммиака 681 кг/м3 определяем массу жидкой фазы в емкости — 34050 кг.
5. По формуле (25) определяем массу аммиака, поступающую в первичное облако в виде газа, образовавшегося при вскипании перегретой жидкой фазы — 6577,65 кг.
6. По формуле (26) определяем массу аммиака, поступающую в первичное облако в виде аэрозоля, образовавшегося при вскипании перегретой жидкой фазы, — 6577,65 кг. Согласно (27) это же количество жидкого аммиака задаем в первичном облаке в начале рассеяния.
7. По формуле (29) определяем площадь пролива — 613,65 м2. Эту же величину согласно правилу, приведенному после формулы (29), принимаем в качестве площади контакта жидкой фазы с подстилающей поверхностью.
8. По формуле (31) давление насыщенных паров при температуре окружающего воздуха — 8718 мм. рт. ст.
9. По формуле (30) определяем время интенсивного кипения, в течение которого ОВ поступает в первичное облако из пролива, — 13 с.
10. По формуле (28) определяем массу газа, поступившего в первичное облако от пролива при его кипении, — 183,73 кг.
11. По формуле (23) определяем массу аммиака в первичном облаке — 13732,2 кг.
12. Согласно формуле (35) задаем равными 0 кг/с общие расходы аммиака и начальный расход жидкой фазы аммиака в следующих вторичных облаках:
в облаке, образующемся при истечении жидкой фазы;
в облаке, образующемся при истечении жидкой фазы после отсечения аварийного участка;
в облаке, образующемся при истечении газовой фазы при наличии пролива;
в облаке, образующемся при истечении газовой фазы в отсутствие пролива;
в облаке, образующемся при испарении из оборудования.
Согласно формуле (37) также полагаем равным 0 секунд длительности выброса аммиака в соответствующих вторичных облаках.
Согласно формулам (40), (44) и (45) также полагаем равным нулю начальные плотности аммиака в этих же вторичных облаках, их полуширину и высоту.
13. По формуле (39) определяем плотность аммиака в газовой фазе (при давлении 101325 Па (1 атм) и температуре кипения -33,6 °C) — 0,865 кг/м3. Эта же величина полагается равной начальной плотности во вторичном облаке, образующемся при испарении с пролива. Начальную температуру в первичном и вторичном облаках определяем как температуру кипения 240 К.
14. По формуле (38) определяем начальную плотность первичного облака — 1,76 кг/м3.
15. По формуле (41) определяем начальный радиус и высоту первичного облака, они равны 13,55 м.
16. По формуле (42) определяем полуширину вторичного облака, образующегося на стадии испарения аммиака из пролива в отсутствие газообразного выброса, — 12,39 м.
17. По формуле (43) определяем высоту, она равна 1,19 м.
18. По формуле (34) определяем скорость испарения с поверхности пролива и расход аммиака во вторичном облаке, образующемся на стадии испарения из пролива, она равна 4,52 кг/с. При расчете испарения использовалась начальная эффективная скорость вторичного облака, образующегося на стадии испарения из пролива, которая равна 0,18 м/с. Согласно формуле (35) расход аммиака в жидкой фазе во вторичном облаке, образующемся на стадии испарения ОВ из пролива, полагается равным нулю.
19. По формуле (36) определяем время испарения пролива — 4580 с.
20. Характерный размер шероховатости согласно пункту 24 настоящего Руководства определяем по таблице № 7-3 - 0,55 м — центры малых городов.
21. Согласно пункту 23 настоящего Руководства определяем показатель степени в степенной зависимости скорости ветра от высоты — 0.655. С помощью этого показателя будет аппроксимироваться скорость в атмосферном слое 20 м над поверхностью земли.
22. По формуле (94) с использованием таблицы № 7-6 определяем масштаб Монина — Обухова — 23,49 м.
23. По формуле (95) определяем динамическую скорость — 0,07 м/с.
24. С использованием таблицы № 7-7 определяем коэффициент для расчета рассеяния.
25. Для определения концентраций аммиака в атмосфере при распространении первичного облака по формулам (105) - (106) численно решаем систему обыкновенных дифференциальных уравнений (107) - (115) с использованием уравнений (101) - (104).
26. На стадии испарения из пролива плотность аммиака в образующемся вторичном облаке 0,87 кг/м3, что меньше плотности окружающего воздуха (1,16 кг/м3), поэтому для расчета концентраций во вторичном облаке используем соотношения из методики «ТОКСИ-2».
27. Согласно пункту 40 настоящего Руководства сравнением с пороговыми и смертельными токсодозами для аммиака определяем расстояние, в пределах которого будет наблюдаться соответствующее поражение.
На поверхности земли протяженность зоны смертельного поражения составляет 212 м в направлении ветра и 160 м в направлении против ветра. Ширина зоны смертельного поражения составляет 393 м на 40 м от места выброса.
На поверхности земли протяженность зоны порогового воздействия составляет 863 м в направлении ветра и 385 м в направлении против ветра. Ширина зоны смертельного поражения составляет 1186 м на 360 м от места выброса.
Результаты расчета зон смертельного поражения и порогового воздействия приведены на рисунках 9-7, 9-8. На рисунке 9-7 показаны границы этих зон на поверхности земли, а на рисунке 9-8 - в вертикальной плоскости, проходящей через ось выброса. Положению места выброса на этих рисунках соответствуют пространственные координаты (0,0,0).
Рисунок 9-7. Зоны смертельного поражения и порогового воздействия на поверхности земли при выбросе аммиака (пример 3)
Рисунок 9-8. Зоны смертельного поражения и порогового воздействия в вертикальной плоскости, проходящей через ось выброса, при выбросе аммиака (пример 3)
Пример 4. Происходит частичное разрушение трубопровода, идущего от емкости с жидким пропаном. Трубопровод содержит жидкую фазу. В результате происходит выброс пропана.
Емкость содержит 4 тонны жидкого пропана. Трубопровод диаметром 5 см идет от нижней части емкости. Высота уровня жидкой фазы над трубопроводом 0,5 м. Давление в емкости 9,2 атм, температура — 18,5 °C.
Аварийное отверстие (разгерметизация) в трубопроводе имеет эффективный диаметр 2 см (площадь 3,14 см2). Место разрушения трубопровода находится на расстоянии 41 метра от емкости. Предполагается, что через 60 секунд после разрушения происходит отсечение аварийного участка.
Выброс происходит на песчаную почву, рассеяние происходит над ровным полем с травой высотой до 15 см. Температура воздуха и почвы 20,5 °C. Скорость ветра на высоте 10 м — 6 м/с. Класс устойчивости атмосферы — D.
Оценить размеры пожаровзрывоопасной зоны, возникающей в результате такой аварии.
При решении используем характеристики пропана, приведенные в таблице № 7-1.
Выбираем тип подстилающей поверхности, над которой происходит рассеяние «трава до 15 см», с характерным размером шероховатости 0,01 м.
При расчете кипения пролива при выпадении жидкости на подстилающую поверхность согласно исходным данным в качестве поверхности пролива берется песок со свойствами, приведенными в таблице № 7-8.
1. По условиям хранения и характеру разрушения определяем тип сценария аварийной ситуации - 4.
2. По формуле (39) определяем плотность газовой фазы пропана при 1 атм и температуре кипения — 2,29 кг/м3.
3. В качестве вспомогательной величины определяем плотность газовой фазы пропана в оборудовании — 14,16 кг/м3.
4. По вспомогательной формуле (46) определяем давление насыщенных паров пропана при существующей температуре в оборудовании — 774463 Па.
5. По вспомогательной формуле (46) определяем плотность насыщенных паров пропана при существующей температуре в оборудовании — 14,06 кг/м3.
6. В качестве вспомогательной величины определяем плотность газовой фазы при температуре кипения и давлении насыщенного пара, соответствующем температуре жидкости в оборудовании, — 17,76 кг/м3.
7. По формуле (48) определяем общий (жидкости и газа) расход в выбросе из разрушенного оборудования (до отсечения аварийного участка) — 4,97 кг/с.
8. По формуле (49) определяем общий (жидкости и газа) расход в выбросе из разрушенного оборудования (после отсечения аварийного участка) — 4,95 кг/с.
9. По формуле (31) определяем давление насыщенных паров пропана при температуре окружающего воздуха 6125 мм. рт. ст.
10. По формуле (58) определяем расход газа, образовавшегося при вскипании перегретой жидкости, в выбросе из разрушенного оборудования (до отсечения аварийного участка) — 1,52 кг/с.
11. По формуле (58) определяем расход жидкости (в виде аэрозоля, образовавшегося при вскипании перегретой жидкости) в выбросе из разрушенного оборудования (до отсечения аварийного участка) — 1,52 кг/с. Эту же величину полагаем равной расходу пропана в жидкой фазе (в виде аэрозоля) в начальном сечении вторичного облака, образующегося на стадии истечения жидкой фазы (до отсечения аварийного участка).
12. С помощью формулы (53) выполняем предварительный расчет времени выравнивания скоростей кипения и испарения за счет тепломассообмена с воздухом — 1,8 с.
13. По формуле (61) определяем расход пропана в газовой фазе, образовавшейся при вскипании перегретой жидкости, в выбросе из разрушенного оборудования (после отсечения аварийного участка) — 1,51 кг/с.
14. По формуле (62) определяем расход жидкого пропана (в виде аэрозоля, образовавшегося при вскипании перегретой жидкости) в выбросе из разрушенного оборудования (после отсечения аварийного участка) — 1,51 кг/с. Эту же величину полагаем равной расходу пропана в жидкой фазе (в виде аэрозоля) в начальном сечении вторичного облака, образующегося на стадии истечения жидкой фазы (после отсечения аварийного участка).
15. Определяем время испарения из пролива всего выброшенного пропана, в том числе после отсечения аварийного участка, — 651 с.
16. По формуле (55) проводим предварительный расчет времени формирования первичного облака на стадии истечения жидкости из разрушенного оборудования до отсечения аварийного участка — 1,8 с.
17. По формуле (56) делаем предварительный расчет времени формирования первичного облака на стадии истечения жидкости из разрушенного оборудования после отсечения аварийного участка — 0 с.
18. По формуле (54) определяем площадь пролива пропана на стадии формирования первичного облака — 0,14 м2. После формулы (56) приравниваем эту величину к площади контакта пролива с подстилающей поверхностью на стадии формирования первичного облака.
19. С помощью формулы (53) окончательно определяем время формирования первичного облака — 0,12 с. Согласно расчету по формуле (55) первичное облако формируется на стадии истечения жидкости из разрушенного оборудования до отсечения аварийного участка. Время формирования первичного облака на стадии истечения жидкости из разрушенного оборудования после отсечения аварийного участка согласно формуле (56) равно 0 с.
20. По формуле (69) определяем длительность стадии истечения жидкости из разрушенного оборудования до отсечения аварийного участка — 59,88 с.
21. По формуле (69) определяем длительность стадии истечения жидкости из разрушенного оборудования после отсечения аварийного участка — 8,27 с.
22. По формуле (50) определяем массу пропана в первичном облаке — 0,38 кг. Согласно формуле (51) в начальный момент времени в первичном облаке находится 0.19 кг жидкого пропана. По формуле (79) рассчитываем начальную плотность первичного облака 4,56 кг/м3. По формуле (84) определяем радиус и высоту первичного облака — 0,3 м. Начальная температура в первичном облаке составляет 231 К.
23. По формуле (72) определяем площадь пролива пропана — 5,17 м2.
24. По формуле (85), первый вариант, определяем полуширину вторичного облака, образующегося на стадии истечении пропана в жидкой фазе (до отсечения аварийного участка), — 1.13698 м. По формуле (88) определяем высоту этого облака — 0,13 м.
25. По формуле (57) расход пропана во вторичном облаке, образующемся на стадии истечении жидкой фазы (до отсечения аварийного участка), составляет 3,14 кг/с. Начальный расход пропана в газовой фазе в этом вторичном облаке составляет 1,62 кг/с. Начальная скорость перемещения этого облака составляет 2,37 м/с. По формуле (80) рассчитываем начальную плотность среды в этом облаке — 4,44 кг/м3. Начальная температура во вторичном облаке составляет 231 К.
26. По формуле (86), первый вариант, определяем полуширину вторичного облака, образующегося на стадии истечении пропана в жидкой фазе (после отсечения аварийного участка), — 1,14 м. По формуле (89) определяем высоту этого облака — 0,13 м.
27. По формуле (57) расход пропана во вторичном облаке, образующемся на стадии истечении жидкой фазы (после отсечения аварийного участка), составляет 3,13 кг/с. Начальный расход пропана в газовой фазе в этом вторичном облаке составляет 1,62 кг/с. Начальная скорость перемещения этого облака составляет 2,37 м/с. По формуле (81) рассчитываем начальную плотность среды в этом облаке — 4,44 кг/м3. Начальная температура во вторичном облаке составляет 231 К.
28. По формулам (73), (75) и (78) получаем, что длительности таких этапов, как истечения из разрушенного оборудования газа при наличии пролива, истечение газа в отсутствие пролива и испарение из емкости, равны нулю, то есть эти стадии в рассматриваемой аварийной ситуации отсутствуют.
29. По формуле (87) определяем полуширину вторичного облака, образующегося на стадии испарения пропана из пролива в отсутствие выброса из оборудования, — 1,14 м. По формуле (91) определяем высоту этого вторичного облака — 0,01 м. Средняя начальная скорость перемещения этого облака составляет 1,46 м/с. По формуле (66) определяем расход пропана в этом облаке — 0,08 кг/с. Согласно формуле (63) жидкая фаза в этом вторичном облаке отсутствует. Согласно формуле (83) начальная плотность этого облака составляет 2,29 кг/м3. Начальная температура во вторичном облаке составляет 231 К.
30. По формуле (71) определяем время испарения пролива пропана после окончания истечения жидкой фазы — 1462 с.
31. Характерный размер шероховатости согласно пункту 21 настоящего Руководства определяем по таблице № 7-1 Приложения № 7 настоящего Руководства — 0,01 м.
32. Согласно пункту 23 настоящего Руководства определяем показатель степени в степенной зависимости скорости ветра от высоты — 0,19. С помощью этого показателя будет аппроксимироваться скорость в атмосферном слое 20 м над поверхностью земли.
33. По формуле (94) с использованием таблицы № 7-6 Приложения № 7 настоящего Руководства определяем масштаб Монина — Обухова — бесконечно большая величина.
34. По формуле (95) определяем динамическую скорость — 0,35 м/с.
35. С использованием таблицы № 7-7 Приложения № 7 настоящего Руководства определяем необходимый коэффициент для расчета рассеяния.
36. Для определения концентраций пропана в атмосфере при распространении первичного облака (105) - (106) численно решаем систему обыкновенных дифференциальных уравнений (107) - (115) с использованием уравнений (101) - (104).
37. Согласно разделу VI настоящего Руководства определяем размеры зон со взрывоопасной концентрацией.
Результаты расчета зон с концентрацией ВКПР, НКПР и 0,5 НКПР представлены на рисунках 9-9, 9-10. На рисунке 9-9 показаны границы этих зон на поверхности земли, а на рисунке 9-10 - в вертикальной плоскости, проходящей через ось выброса. Положению места выброса на этих рисунках соответствуют пространственные координаты (0,0,0).
Значения концентраций, кг/м3 (плоскость XY, Z = 0.0 м)
Рисунок 9-9. Зоны достижения концентрацией ВКПР, НКПР и 0,5 НКПР на поверхности земли при выбросе пропана (пример 4)
Значения концентраций, кг/м3 (плоскость XZ, Y = 0.0 м)
Рисунок 9-10. Зоны достижения концентрацией ВКПР, НКПР и 0,5 НКПР в вертикальной плоскости, проходящей через ось выброса, при выбросе пропана (пример 4)
Приложение № 10 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
БЛОК-СХЕМА РАСЧЕТА ПОСЛЕДСТВИЙ АВАРИЙНОГО ВЫБРОСА ОПАСНЫХ ВЕЩЕСТВ
Рисунок 10-1. Блок-схема расчета последствий аварийного выброса опасных веществ
Приложение № 11 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО УЧЕТУ СЛОЖНЫХ МНОГОКОМПОНЕНТНЫХ СМЕСЕЙ
Руководство применимо для случаев выброса сложных многокомпонентных смесей различных ОВ, в том числе смесей углеводородов, при этом рекомендуется рассматривать один из двух вариантов:
Используется «модельное» вещество, свойства которого корректно описывают соответствующие физические процессы, соответствуют имеющимся экспериментальным данным по испарению и кипению многокомпонентных смесей. При этом рекомендуется учитывать, что в случае многокомпонентных смесей состав жидкостной и газовой фаз может существенно различаться.
В качестве коэффициента \gamma (Cp/Cv) выбирается показатель адиабаты, соответствующий составу газовой фазы либо находящейся изначально в оборудовании (сценарии 1, 2 раздела III настоящего Руководства), либо образующейся после испарения (кипения) перегретой жидкой фазы (сценарии 3, 4 раздела III настоящего Руководства).
В качестве плотности жидкости (\rho_{\text{ж}}) выбирается плотность жидкой фазы, остающейся после вскипания жидкости, находящейся в оборудовании.
В качестве температуры кипения (Tкип) выбирается температура, соответствующая равновесному состоянию «пар — жидкость» при атмосферном давлении в случае изэнтропического расширения ОВ от условий его содержания в оборудовании.
В качестве теплоемкости жидкости (Cv) выбирается теплоемкость жидкой фазы, остающейся после вскипания жидкости, находящейся в оборудовании.
В качестве теплоты испарения (\Delta H_{\text{исп}}) выбирается величина, обеспечивающая испарение наблюдаемого на практике количества перегретой жидкости для выбранных значений температуры кипения (Tкип) и теплоемкости жидкости (Cp).
В качестве молярной массы (\mu) выбирается величина, соответствующая составу газовой фазы либо находящейся изначально в оборудовании (сценарии 1, 2 раздела III настоящего Руководства), либо образующейся после испарения (кипения) перегретой жидкой фазы (сценарии 3, 4 раздела III настоящего Руководства).
В качестве теплоемкости газа (Cv) выбирается величина, соответствующая составу газовой фазы либо находящейся изначально в оборудовании (сценарии 1, 2 раздела III настоящего Руководства), либо образующейся после испарения (кипения) перегретой жидкой фазы (сценарии 3, 4 раздела III настоящего Руководства).
При определении свойств газовой фазы рекомендуется учитывать изменение \gamma , \mu и Cv при изменении состава газовой фазы, которое происходит во времени при испарении, с одной стороны, образовавшегося пролива, а с другой стороны, аэрозольных капель, находящихся в дрейфующих облаках.
При определении свойств жидкости рекомендуется учитывать изменение \rho_{\text{ж}} , Tкип, Cp и \Delta H_{\text{исп}} при изменениях состава жидких фаз, которые происходят во времени при испарении, с одной стороны, в образовавшемся проливе, а с другой стороны, в аэрозольных каплях, находящихся в дрейфующих облаках.
Используются результаты решения различных термодинамических задач для расчета:
фазового равновесия сложных многокомпонентных систем при заданной температуре для определения состава газовой и жидкостной фракций, давления насыщения;
адиабатического расширения газов;
адиабатического расширения вскипающей жидкости;
интенсивности кипения пролива;
параметров ОВ в облаке (Приложение № 12 к настоящему Руководству).
При этом в обоих вариантах для сложных многокомпонентных смесей:
скорость испарения из проливов и из капель (находящихся в дрейфующих облаках), состав и свойства образующейся газовой фазы определяются на основе эмпирических данных либо расчетным путем с соответствующим обоснованием метода расчета (термодинамика, молекулярная динамика);
концентрационные пределы распространения пламени (НКПР и ВКПР) и показатели токсичности (пороговая токсодоза PCt50, смертельная токсодоза LCt50, коэффициенты пробит-функции) определяются на основе эмпирических данных;
стехиометрическая концентрация определяется на основе соответствующего состава газовой фазы.
Приложение № 12 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО ПРИМЕНЕНИЮ ТЕРМОДИНАМИЧЕСКИХ РАСЧЕТОВ
При использовании Руководства рекомендуется применять также термодинамические расчеты, которые позволят для системы, находящейся изначально при определенных давлении и температуре и состоящей в общем случае из множества различных веществ и фаз, а в простейшем случае из одного вещества, определить ее параметры:
состав по веществам (компонентам);
полный набор термодинамических параметров (давление, температура, плотность, внутренняя энергия, энтропия, энтальпия, энергия Гиббса и т.д.);
набор первых производных термодинамических параметров (например, скорость звука — это квадрат первой производной давления по плотности).
Термодинамические расчеты могут быть проведены корректно в случае, когда:
известен начальный состав системы (либо молекулярный — набор и количество различных индивидуальных веществ, либо атомарный — брутто формула смеси различных веществ);
известны два любых параметра смеси (например, давление и температура);
задана степень равновесия (например, полное фазовое равновесие при отсутствии химических реакций).
При расчете всегда предполагается отсутствие химических реакций.
При проведении термодинамических расчетов (Приложение № 11 к настоящему Руководству) решаются следующие термодинамические задачи:
для расчета параметров фазового равновесия сложных многокомпонентных систем решается T-P задача определения состояния при данном давлении и температуре, исходный состав полагается равным составу жидкой фазы в неразрушенном оборудовании;
для расчетов адиабатического расширения газов решается S-P задача определения состояния при заданном значении давления и энтропии, при этом энтропия полагается равной энтропии в исходном состоянии (в неразрушенном оборудовании), а давление — равным атмосферному, исходный состав полагается равным составу газовой фазы в неразрушенном оборудовании;
для расчета адиабатического расширения вскипающей жидкости решается S-P задача определения состояния при заданном значении давления и энтропии, при этом энтропия полагается равной энтропии в исходном состоянии, а давление — равным атмосферному, исходный состав полагается равный составу в неразрушенном оборудовании;
для расчета скоростей истечения со звуковой скоростью решается S-P задача определения состояния, включая скорость звука, при заданном значении давления и энтропии, при этом энтропия полагается равной энтропии в исходном состоянии, а давление — равным атмосферному, исходный состав полагается равным составу в неразрушенном оборудовании;
для расчета интенсивности кипения пролива решается H-P задача определения состояния системы при заданных значениях энтальпии (энергии) и давлении, при этом энергия полагается равной сумме исходной внутренней энергии и энергии, подводимой от подстилающей поверхности, давление полагается равным атмосферному, исходный состав полагается равным составу жидкой фазы, проливающемуся на поверхность в начальный момент времени;
для определения состояния в облаке решается H-P задача определения состояния системы при заданных значениях энтальпии (энергии) и давлении, при этом энергия (энтальпия) полагается равной сумме исходной внутренней энергии ОВ, распространяющегося в атмосфере, внутренней энергии (энтальпии) воздуха, смешивающегося с ОВ в ходе распространения с выбросом, а также энергии, поступающей за счет теплообмена с подстилающей поверхностью.
Приложение № 13 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО УЧЕТУ ПОДМЕШИВАНИЯ ВОЗДУХА НА НАЧАЛЬНОМ ЭТАПЕ ВЫБРОСА
При использовании Руководства допускается рассматривать также выбросы, в которых имеет место разбавление воздухом уже на начальном этапе распространения. Это разбавление происходит за счет турбулизации потоков вблизи границ раздела сред (ОВ — воздух). Такая турбулизация обусловлена существованием градиентов скорости и давления в высокоскоростных потоках, формирующихся на начальном этапе вблизи места выброса: струи газа и двухфазных сред, адиабатически расширяющихся объемов.
Учет разбавления производится через задание величины Q_{i}^{\text{разб}} — массы воздуха в первичных облаках в начале рассеяния, и q_{i}^{l\text{разб}} — расходы воздуха во вторичных облаках в начале рассеяния. Эти величины определяются путем умножения соответствующих масс (расходов) ОВ в первичных (Qi) или вторичных (q_{i}^{l}) облаках на коэффициенты \alpha_{i\text{возд}} для первичных облаков или \alpha_{i\text{возд}}^{l} для вторичных: Q_{i}^{\text{разб}} = \alpha_{i\text{возд}} \cdot Q_{i} ; q_{i}^{l\text{разб}} = q_{i}^{l} \cdot \alpha_{i\text{возд}}^{l} .
За счет динамических и (или) гравитационных эффектов при истечении ОВ в атмосферу может произойти смещение в пространстве положения первичного или вторичного облака по отношению к месту разгерметизации оборудования, то есть вместо начального положения облака с координатами (0,0,H/2) принимаются координаты ( \Delta X_{\text{обл}} , \Delta Y_{\text{обл}} , \Delta Z_{\text{обл}} ).
Величины коэффициентов \alpha_{i\text{возд}} и \alpha_{i\text{возд}}^{l} , а также параметры смещения облака \Delta X_{\text{обл}} , \Delta Y_{\text{обл}} и \Delta Z_{\text{обл}} обосновываются и определяются по расчетам начальных этапов распространения выбросов. Эти расчеты и обоснования могут быть сделаны:
на основе расчета по интегральным моделям теории турбулентных струй;
на основе расчета методами прямого численного моделирования (методы вычислительной газодинамики).
Зная величины Q_{i}^{\text{разб}} и q_{i}^{l\text{разб}} с использованием Приложения № 8 к настоящему Руководству, определяются параметры облака, разбавленного воздухом, непосредственно в начале рассеяния, а именно:
плотность смеси после подмешивания воздуха ( \rho_{\text{эфф}i} — для первичного облака или \rho_{\text{эфф}i}^{l} — для вторичных облаков);
температура смеси после подмешивания воздуха (Tэффi — для первичного облака или T_{\text{эфф}i}^{l} — для вторичных облаков);
масса ОВ в жидкой фазе после подмешивания воздуха (Qжi — для первичного облака или q_{\text{ж}i}^{l} — для вторичных облаков).
При расчете вместо используемых в формулах (145) - (153) величин Qсумi для первичных облаков и q_{\text{сум}i}^{l} — для вторичных облаков рекомендуется использовать рассчитанные вышеуказанным образом величины Q_{i}^{\text{разб}} + Q_{i} для первичных и q_{i}^{l\text{разб}} + q_{i}^{l} для вторичных облаков.
Эти же величины массы первичных (Q_{i}^{\text{разб}} + Q_{i}) и расходов вторичных (q_{i}^{l\text{разб}} + q_{i}^{l}) облаков с учетом подмешенного на начальном этапе воздуха рекомендуется использовать и для расчета начальных размеров облаков в формулах (8), (20), (39), (41), (82), (83), (84), (86), (87), (88), (89), (90) и (91), заменяя ими соответствующие величины массы без учета подмешивания воздуха.
Также для расчета начальных размеров облаков в (8), (20), (39), (41), (82), (83), (84), (86), (87), (88), (89), (90) и (91) рекомендуется вместо плотностей чистых ОВ (без подмешивания воздуха) использовать плотности \rho_{\text{эфф}i} (для первичного облака) и \rho_{\text{эфф}i}^{l} (для вторичных облаков), рассчитанные как указано в Приложении № 8 к настоящему Руководству с учетом подмешенного на начальном этапе воздуха.
Расчеты процессов распространения выбросов без подмешенного на начальном этапе воздуха или с ним идентичны. Они выполняются в порядке, рекомендованном в настоящем Руководстве, при этом при учете подмешивания воздуха для вторичных облаков рекомендуется учитывать в размере зон поражения расстояние, на которое смещается в потоке ветра облако-шлейф на начальном этапе подмешивания.
Приложение № 14 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО УЧЕТУ ТЕМПЕРАТУРНОЙ ЗАВИСИМОСТИ СВОЙСТВ ОПАСНОГО ВЕЩЕСТВА
При использовании Руководства рекомендуется также учитывать зависимость свойств ОВ от температуры, в том числе для многокомпонентных систем, при этом рекомендуется рассматривать один из двух вариантов:
Использовать термодинамические расчеты согласно Приложениям № 11 и 12 к настоящему Руководству. При этом в термодинамических расчетах рекомендуется учитывать температурную зависимость свойств ОВ от температуры.
Использовать в качестве теплоемкостей, показателей адиабат и теплот испарения осредненные значения для соответствующих температурных зависимостей. Плотность же жидкой фазы в этом случае рекомендуется брать с учетом температуры, при которой она находится.
Приложение № 15 к Руководству по безопасности «Методика моделирования распространения аварийных выбросов опасных веществ», утвержденному Приказом Федеральной службы по экологическому, технологическому и атомному надзору от 20 апреля 2015 г. № 158
РЕКОМЕНДАЦИИ ПО УЧЕТУ ИЗМЕНЕНИЯ ПЛОТНОСТИ ПРИ РАСПРОСТРАНЕНИИ ВЫБРОСА
При использовании Руководства рекомендуется также учитывать изменение плотности выброса, когда при определенных условиях распространения выброс перестает быть «тяжелым» газом и плотность облака станет меньше плотности воздуха.
В таком случае для завершения расчета допускается проводить расчеты по данному Руководству, но при этом:
формулы (109) и (126), описывающие гравитационное растекание, не используются,
величина ri (или, соответственно, b_{i}^{l} ) фиксируются и перестают изменяться;
величины эффективного радиуса (Rэффi) (или, соответственно, полуширины (B_{\text{эфф}i}^{l}) ) рассчитываются по формулам (99) и (118).
Если на стадии рассеяния происходит повторное изменение плотности — выброс снова становится «тяжелым» газом, то процедура расчета возвращается к стандартной процедуре расчета распространения «тяжелого» газа, описанной в данном Руководстве.
Данная процедура повторяется до завершения расчета.