В настоящее время водородсодержащий газ (ВСГ) широко используется в различных отраслях промышленности, причем с увеличением производственных нужд увеличивается объемы его производства [1]. По прогнозам компании по анализу данных и консалтингу GlobalData, глобальные мощности по производству «низкоуглеродного» водорода с большой вероятностью достигнут 14 млн т в год к 2030 г., что в 20 раз больше, чем на сегодняшний день (рис. 1) [2].С развитием производственной водородной инфраструктуры наблюдается параллельное развитие автоматизированных систем управления технологическим процессом (АСУ ТП), в связи с чем происходит последовательное усложнение задач, стоящих перед системами управления от управления отдельными установками и параметрами, к автоматизации процессом в целом. Применение современных АСУ ТП дает возможностьне только с высокой эффективностью реализовывать контроль и управление на производстве, но и отчасти избежать человеческих ошибок в оперативном управлении. К актуальным вопросам в данном случае относят вопросы развития интеллектуальных технологий АСУ ТПв направлении разработки алгоритмов реагирования на возникающие критические ситуации в динамическом режиме.Актуальность исследования обусловлена тем, что применяемая на технологических производствах АСУ ТП не учитывает всех существующих показателей пожарной опасности обращаемых в технологических аппаратах веществ, в данном случае водородсодержащего газа, что не позволяет полноценно оценивать риски возникновения аварийных ситуаций. Внедрение расчетной оценки уровня безопасности технологического оборудования, использующего водород, позволит обеспечить существенное снижение аварийности и прочих производственных рисков. Работа направлена на решение данной проблемы.В странах с развитой водородной инфраструктурой данная проблема частично решена. Так, в работе [3] представлен обзор критериев опасности, которые могут быть использованы при количественной оценке пожарного уровня пожарной опасности водородных установок.В исследовании [4] предлагается модель оценки рисков на основе индекса, оценивающая водородную инфраструктуру, используя относительный рейтинг рисков, таких как: производство, хранение и транспортировка водорода.Разработан программный инструментарий HyRAM для оценки безопасности инфраструктуры заправки топливом и хранения водорода. Его также используют для количественной оценки пожарного риска с интегрированным анализом последствий аварии [5].Целью работы является разработка и развитие комплексного подхода в обеспечении пожарной безопасности для объектов водородной энергетики, позволяющей оценивать состояние взрывопожароопасного технологического оборудования в целом, в том числев динамике.Поставленные задачи представлены в виде разработки критериальной модели, а также программы для ЭВМ, содержащих в своем составе: описание и построение параметров критериальной модели, а также разработку методики отнесения технологического участка, на которой размещено водородное оборудование, к определенной категории уровня пожарной опасности (риска).Методы исследованияНа сегодняшний день существует множество теоретических и экспериментальных методов определения показателей пожарной опасности газообразных веществ [6], однако некоторые из них невозможно применить для оценки пожарной опасности водородсодержащего газа из-за существенных отличий физико-химических и пожаровзрывоопасных свойств с углеводородными газами.С точки зрения достижения целей: моделирования безопасности технологического процесса обращения водорода и создания системы обеспечения его безопасности в заданном технологическом процессе следует выделить два вида параметров – управляемые (внутренние (А) и управляющие (внешние (В). Первые влияют на безопасность технологического процесса, улучшая или ухудшая ее, вторые воздействуют на состояние защищенности, в котором управляемые параметры находятся в зоне запрещенных значений. Таким образом, управляемые параметры зависят от управляющих. Следовательно, внутренние параметры в технологическом процессе можно рассматривать как набор показателей пожарной опасности. В соответствии с табл. 1 приложения к Федеральному закону от 22 июля 2008 г. № 123-ФЗ «Технический регламент о требованиях пожарной безопасности» к перечню показателей, необходимых для оценки пожарной опасности газообразных веществ, относятся 19 показателей [7].Любое состояние защищенности технологического процесса можно охарактеризовать следующим конечным множеством параметров: X = (x1, x2, x3, …, xi).При xi≥xikдоп  действующее значение параметра источника опасности превышает допустимое значение, что означает опасное воздействие этого источника опасностина защищаемый объект [8].Естественно, что необходимым условием обеспечения пожарной безопасности объекта является выполнение следующего неравенства:где xi – действующее значение параметра источника опасности; xiдоп – допустимое значение параметра источника опасности.Следовательно, необходимым условием при оценке уровня безопасности технологического оборудования, в котором обращается ВСГ, будет выглядеть следующим образом:Состояние защищенности технологического процесса обращения ВСГ можно охарактеризовать следующим конечным множеством параметров:где Wmin – минимальная энергия зажигания, Дж; φв – верхний концентрационный предел распространения пламени по смеси горючего вещества с воздухом, % (об.); φн – нижний концентрационный предел распространения пламени по смеси газообразного горючего вещества с воздухом, % (об.); φф – минимальная флегматизирующая концентрация флегматизатора, % (об.); ΔP – максимальное избыточное давление взрыва газа, кПа; dр/dτmax – максимальная скорость нарастания давления взрыва водородо-воздушной смеси, кПа/с; Su – нормальная скорость распространения пламени, м/с; φО2 (МВСК) – минимальное взрывоопасное содержание кислорода в горючей смеси, % (об.); φПДГ – концентрационный предел диффузионного горения газовых смесей в воздухе, % (об.); tсам – температура самовоспламенения, К; dкр – критический диаметр длинного цилиндрического огнегасящего канала, мм; δбэмз (БЭМЗ) – максимальный экспериментальный безопасный зазор, мм; u0 – предельная скорость срыва диффузионного горения, м/с; Q – низшая теплота сгорания вещества, МДж/кг; q – удельная теплота сгорания вещества, МДж/кг; Tапп – температура технологического аппарата, К; Т0 – температура вещества при нормальных условиях, 273 К; Pапп – давление в аппарате, кПа; Vапп – объем реакционного сосуда, м3; P0 – атмосферное давление, 101,3 кПа; ρв – плотность воздуха до взрыва при начальной температуре Т0, кг∙м-3; А – кратность воздухообмена, создаваемого аварийной вентиляцией, с–1; τпост – продолжительность поступления ВСГ в объем помещения, с (120 с, если вероятность отказа системы автоматики превышает 0,000001 в год и не обеспечено резервирование ее элементов; 300 с при ручном отключении трубопроводов).Результаты исследования и их обсуждениеВ табл. 1 представлены основные параметры, влияющие на безопасность технологического оборудования, в котором обращается ВСГ, их допустимые значения, методы расчета, а также исходные данные, задаваемые оператором АСУ ТП для оценки уровня пожарной опасности.Допустимые значения основных параметров, влияющих на безопасность технологического оборудования, в котором обращается ВСГ ПараметрПараметры, задаваемые операторомДопустимое (предельное) значениеЕдиницы измеренияРасчетная формулаφнTапп (задается, в К), jн0 (4,0 % (об.),Тг (2318 К)4,0 [9]% (об.)φн=φн01-Tапп-Т0Тг-Т0  [14]φвTапп (задается, в К), jв0 (75,0 % (об.),Тг (2318 К)75,0 [9]% (об.)φв=φв01+Tапп-Т0Тг-Т0  [14]δбэмзdкp (м)0,00029 [10]мδбэмз=dкр3  [17]∆Рb (м), l (м), h (м),Qн (кДж/кг),Сp (кДж/(кг∙К),ρв (кг/м3), Z (1), Тв (К), P0 (101,3 кПа),Pапп (кПа), Pтруб (кПа), Vапп (м3), τпост (120 с – автоматическое отключение, 300 с – ручное отключение), qвсг (м3/с), А (0,0125 с-1), rподв (м), Lподв (м),rотвод (м), Lотвод (м)16 [11]кПа∆Р=m∙Hт∙Z∙P0Vсв∙ρв∙Cр∙T0∙1Kн  [19]φф:1) состав(85 % (масс.) CO2 и 15 % (масс.) хладон 114В2)2) CO23) N2Выбирается одиниз флегматизаторов, использующихсяв АУПТ 38,474,491,2 [12] % (об.)φф=i=1nχiφi  [15]SuPапп (кПа), Tапп (К), Su0 (3,06 м/с), m (2,23),n (-0,12), P0 (101,3 кПа), Т0 (273 К)2,1 [13]м/сSu=Su0PаппP0nTаппT0m  [16]dр/dτVапп (м3), Su (м/с),Рапп (кПа)dPdtmax(доп)=81,55∙Pапп(доп)Su(доп)Vапп13кПа/сdPdtmax=81,55∙PнSuV13  [17]tсамTапп (К)783 [12]КQсамб=Tапп≤tсам Wminа (0,5), qг (Дж/м3),lк (м)6800 [12]ДжWmin = a qг lk3  [17]φО2:1) CO22) N2 φф (соответствует концентрации выбранного флегматизатора, % (об.)  4,21,8 [12] % (об.)φО2=100-φф4,844  [15]ПараметрПараметры, задаваемые операторомДопустимое (предельное) значениеЕдиницы измеренияРасчетная формулаu0K (0,75);µ (2,016 кг/моль); Cp (20,80 кДж/(моль·K); tапп (K); P0 (101,3 кПа); Pапп (кПа), γ (1,410)600 [13]м/сu0=K2CpTаппM∙1-P0P0+Pапп1-1γ  [18] Примечание: ТГ – температура горения смеси в воздухе, К; b, l, h – ширина, длина и высота помещения, м; Сp – теплоемкость воздуха, кДж∙кг–1 ∙К–1; Pтруб – давление в трубопроводе, кПа;Z – коэффициент участия ВСГ в горении; qвсг – расход ВСГ, определяемый в соответствиис технологическим регламентом в зависимости от давления в трубопроводе, его диаметра, температуры газовой среды и т.д., м3 ∙с-1; rподв(отвод) – внутренний радиус трубопроводов, м;Lподв(отвод) – длина трубопроводов от аварийного аппарата до задвижек, м; Su0 – известное значение нормальной скорости при давлении P0 и температуре Т0, равное 3,06 м/с; n и m – соответственно барический и температурный показатели; qг – удельное объемное количество тепла, необходимое для нагрева горючей смеси от t0 до tсам, Дж/м3; а – коэффициент пропорциональности; lк – критический зазор зажигания, мм; K – коэффициент расхода форсунки, равный 0,75; Cp – теплоемкость исследуемой газовой смеси, 20,80 кДж/(моль·K); γ – показатель адиабаты (коэффициент Пауссона, равный 1,410)На базе высокоуровневого языка программирования Python был разработан программный код, позволяющий упростить определение параметров критериальной модели.В диалоговом окне вводятся исходные значения параметров, которые подставляютсяв расчетные формулы и сравниваются с допустимыми значениями параметров пожарной опасности. При превышении того или иного допустимого значения переменной присваивается 1, если допустимое значение не превышено переменной – присваивается 0. По результатам сравнения критериев отнесения технологического участка, на котором размещено водородное оборудование, к определенной категории уровня пожарной опасности (риска), технологическому участку присваивается категория(табл. 2). Таблица 2 Критерии отнесения технологического участка, на котором размещено водородное оборудование, к определенной категории уровня пожарной опасности Категория уровня пожарной опасностиКритерииВысокийПревышаются более семи допустимых значений основных параметровСреднийПревышаются не более семи допустимых значений основных параметровУмеренныйПревышаются не более пяти допустимых значений основных параметровНизкийПревышаются не более трех допустимых значений основных параметров Фрагмент критериальной модели оценки уровня безопасности технологического оборудования, где применяется водород, в виде исходного кода разработанной программы для электронной вычислительной машины (ЭВМ) на базе языка программирования Python представлен на рис. 2.Заключение Результаты исследования реализованы в виде программы для ЭВМ, которая прошла государственную регистрацию в Федеральной службе по интеллектуальной собственности [20]. В рамках проведенного исследования была достигнута цель работы путем решения следующих задач: представлены существующие модели оценки уровня пожарной опасности (риска)на объектах водородной энергетики; отмечена актуальность использования данных моделейв системах безопасности объектов водородной энергетики; разработана критериальная модель оценки уровня пожарной безопасности при эксплуатации технологического оборудованияна объектах водородной энергетики, позволяющая оценивать состояние взрывопожароопасного технологического оборудования в целом, в том числе в динамике. К отличиям разработанной модели от других можно отнести простоту использования, а также оценку уровня пожарной опасности объекта с помощью метода расчета показателей, необходимых для определения пожарной опасности газообразных веществ, установленных Федеральным законом от 22 июля 2008 г. № 123-ФЗ, и в дальнейшем – отнесения объекта к определенной категории уровня пожарной опасности. Рассматриваемая модель апробирована и внедрена на нескольких объектах водородной энергетики, в их числе: электролизный цех Каргалинской ТЭЦ и электролизный цех Балаковской ТЭЦ. По результатам апробации данным цехам присвоена средняя категория уровня пожарной опасности. В будущем планируется апробировать и внедрить данную модель на технологических участках предприятий топливно-энергетического комплекса Санкт-Петербурга в составе систем безопасности АСУ ТП, использующих водородные системы охлаждения турбогенераторов. Также актуальным применением критериальной модели будет являться ее внедрение в систему безопасности АСУ ТП установок паровой конверсии природного газа нефтеперерабатывающих, нефтехимических, химических и металлургических предприятий.