Оценка риска для здоровья персонала при эксплуатации объектов накопления нефтесодержащих отходов

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

Технологический процесс добычи и переработки нефти сопровождается образованием огромного количества нефтесодержащих отходов, которые накапливаются в шламовых амбарах, нефтешламонакопителях, ямах аварийного хранения продукта, полигонах нефтеотходов и др. Подавляющее большинство таких объектов открыто, расположены на местности и активно воздействуют на атмосферу. При выборе способа оценки опасности накопителей следует иметь в виду, что за период их длительной эксплуатации, исчисляемый не одним десятком лет, содержимое таких амбаров превращается в сложную расслоившуюся систему, в которой можно выделить 3 основных слоя: верхний слой (плавающий нефтешлам), средний слой (водную фазу) и нижний слой (донный осадок) [4]. Помимо потенциальной опасности необходимо рассмотреть плохо отражённый в литературе аспект, связанный с токсическим воздействием испаряющихся веществ на здоровье персонала, работающего в районе нефтешламонакопителя. Для прогнозирования экологических последствий в местах накопления нефтесодержащих отходов важную роль играет процесс испарения веществ со свободной поверхности пруда-накопителя и рассеивание в атмосфере. Модели испарения позволяют оценить интенсивность поступления в окружающую среду паров горючих и токсичных веществ, модели рассеивания – образующиеся концентрации испаряющихся веществ [4-6]. Для объектов накопления нефтесодержащих отходов целесообразно использовать комплексную модель оценки риска, основанную на определении весовых коэффициентов, характеризующих опасность различных углеводородов на основе факторного анализа – времени существования накопителя и времени экспозиции для промышленного персонала [4]. Наибольшими недостатками оценки экологического риска для накопителей нефтесодержащих отходов характеризуются методики оценки последствий, связанных с выбросами в окружающее пространство углеводородов, содержащихся в верхнем слое. Это связано с тем, что зачастую в методиках не учитывается движение облака опасных веществ под действием силы тяжести, в связи с чем обоснована необходимость перехода при моделировании подобных процессов от моделей гауссовского типа к моделям рассеяния тяжелого газа, а также разработки соответствующих адекватных программных средств.

Основная опасность нефтешламонакопи-теля связана с возможностью возникновения долговременного токсического выброса небольшой интенсивности. Последствия такой ситуации можно оценить с помощью расчета ущерба жизни и здоровью людей [7]. При расчетах математического ожидания ущерба и потерь принято, что наибольшую потенциальную опасность представляет открытый верхний углеводородный слой нефтешламонакопителя, с поверхности которого испаряются углеводороды и токсические вещества, приводящие к загазованности территории. Вероятность возгорания и взрыва испаряющихся веществ в данном случае очень мала. В рамках реализации вышеупомянутой методологии [7] разработаны алгоритм и компьютерная программа Risk, предназначенная для количественной оценки риска опасных площадочных производственных объектов накопления нефтесодержащих отходов. Программный комплекс оценки последствий выбросов углеводородов с открытой поверхности накопителя на основе совершенствования методики «ТОКСИ» позволяет наиболее точно, по сравнению с другими отечественными методиками в этой области, определять последствия распространения выброса опасных веществ в атмосфере [4, 5]. Необходимость разработки оригинального программного комплекса объясняется тем, что существующие программы по методике «ТОКСИ» не позволяют определить концентрации опасного вещества для площадных объектов накопления нефтесодержащих отходов с учетом их геометрических размеров, формы и времени эксплуатации. Преимущество разработанной авторами программы заключается в возможности нанесения зон распространения опасных веществ, рассчитанных для любого накопителя, на реальную картографическую основу.

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

a

C(x,t) = d dt (x,t)д C(x,t) dt        dx        dx dx

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

C ( x,t ) = 0 = 0 при х=0             (2)

^{C ( x,t} ) = C o при t=o             (3)

— C ( x,t ) = 0

dx         при x=H (4)

где С0 –концентрация в начальный момент времени, определяемая экспериментально; H – толщина верхнего углеводородного слоя накопителя.

Дифференциальное уравнение (1) при граничных условиях (2)-(4) описывает независимое испарение лёгких компонентов с поверхности нефтешламонакопителя и их поступление из глубины верхнего слоя толщины H к поверхности. Выражение (1) представляет собой нелинейное дифференциальное уравнение в частных производных параболического вида, где коэффициент зависит от температуры. Температура в свою очередь зависит от координаты и времени. Аналитическое решение для ограниченной глубины получить не удаётся, поэтому система уравнений решалась численными методами. Для этого необходимо перейти от непрерывных значений аргументов функции концентрации к дискретным. Разобьем слой на n – 1 частей и получим n разных значений функции по глубине слоя, находящихся на расстоянии равном шагу расчётной сетки Δ x . Временной отрезок приводится к дискретному заданием шага расчётной сетки по времени, равному Δ t . В соответствии с введёнными положениями можно перейти от функционального вида к индексам:

Cj,u = c(z,t),Cw,u+1 = c(z + Az,t + Az)

В соответствии выбранной системой отсчёта выражение (1) с учётом граничных условий можно представить следующей системой выражений:

Ci u = C          j e [0;n]

j,u     0 при u=0                  (6)

Ci u = ⁰ u e [0; ^ ) ^j,u     при j=0                      (7)

C j,u + 1

= C j,u +8 1 [ D j,uJ - J ) - D j _ 1,„ ( C J,

-

C J,U + 1

= C j,u -5 1^Di - 1,u ⁽ C j,u - C j - 1,u ) при j=n

при 0

где C0 – начальное значение концентрации; δ1, δ2 – коэффициенты, определяемые выражениями 81 = At/Ax², ⁸2 =At/Ax; Dj,u — значение функции коэффициента диффузии от средней координаты Dj,_u= D(J Ax + 0,5Ax,uAt). Расчёт по формулам (7)-(9) выполняется последовательно от j=0 до j=n для заданного промежутка времени. Интенсивность испарения определялась скоростью диффузионного потока на границе раздела фаз:

Ju ⁸i^Dj^-i,u ⁽Cj,u Cj-u) при о=о (10)

Основной задачей численного решения является получение зависимостей интенсивностей испарения веществ с разной молекулярной массой и толщиной слоя с течением времени с поверхности накопителя. Решения были аппроксимированы следующей функцией для веществ с молекулярной массой μi и глубины слоя H от 5 до 30 см по алгоритму минимизации для каждого i-го компонента:

J⁽C0,i, Hi^,H,t)= al C0,i ^exP^(a2 fit ^{(t)- a}3 ^f2^t(t))•£ ^exp^(-a4 ^{(2n - 1)1,66} t)

n=¹                                  (11)

C где , – начальная концентрация i-го компонента, определяемая из опыта;

ft (t ) = 2cos(t 2t)-[cos([t2 +t1 ]t)+ cos([t2 +t1 ]t)] f2t(t) = cos(T1t + 0,175)

τ1, τ2 – частоты смены суток и времён года, с^-1.

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

LADD =

C,_IM• CR • Ef ^(t1 -1₂)

BW • AT • 365

ности испарения октана при его концентрации в верхнем углеводородном слое накопителя 10 г/м3 и при различной глубине слоя 20 и 30 см. Из результатов расчёта видно, что ключевым параметром,

характеризующим потенциальную опасность накопителя, с течением времени является концен-

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

гдеC_атм– концентрация вещества в загрязненной среде, мг/м³; CR – скорость поступления воздействующей среды; Ef – частота воздействия; t1 – время начала воздействия от создания объекта; t₂ – время окончания воздействия от создания объекта; BW – масса тела человека, кг; AT – период усреднения экспозиции; 365 – число дней в году.

лении в организм нескольких опасных веществ обобщённый канцерогенный риск определяется по формуле CR _z:

f2t(t )= cos(T1t + 0,175)

где CR – риск развития рака у индивидуума; Sf – фактор наклона кг день/мг; LADD – среднесуточная доза в течение жизни мг/(кг день), определяется по формуле:

Рис. 1. Сравнение зависимостей интенсивности испарения октана от глубины верхнего углеводородного слоя в течение 15 лет

Изменение скорости испарения с поверхности накопителя (выброс) с течением времени можно определить как сумму скоростей испарения всех углеводородов:

nвещ

W(H,t) = S^ £ J(Co,i,Hi,H,t)

i=1                                 (16)

где j(Coi,Hi,H,t) определяется выражением (11); Sнак – площадь накопителя.

С учетом формулы (16) функция зависимости изменения концентрации углеводородов в атмосфере в пространстве и во времени имеет следующий вид:

nвещ

w(H,t )=S„a, £ j(Co,i, H i,H,t)

i-1                                 (17)

Для подстановки полученного значения концентрации, его нужно усреднить по временному промежутку экспозиции:

Ga„(x,y,z,H,t)-3,17098.10-8 Sf^,y,z)£j‘“j(Co,.,Hi,H,t)dt t2 t1        i=1                                             (18)

где 3,17098 10-8 – переводной коэффициент по времени;

Согласно [7] величины CR, BW, AT равны 20 м3/сут., 70 кг и 70 лет соответственно. Примем частоту воздействия (Ef), равной 120 дней/год. Подставляем значение величин CR, BW, AT, Ef и выражение (18) в формулу (15) получаем выражение для расчёта среднесуточной дозы:

Ga„ (x,y,z,H,t ) = 3,17098 .10-8 Sf ^,y,z)£ j“j(Go,., H i,H,t)dt t2 t1       i=1                                          (19)

Подставляем полученное выражение (19) и значение Sf=0,035 мг/(кг день) в формулу (14), получаем выражение для расчёта риска здоровью промышленного персонала:

nвещ t2

GR L(x,y,z,H,t„t2)-1,4831.10-,2S„fR (x,y,z)£jt1 j(Co,„ h „H,t)dt i-1

Результаты расчета риска с учётом распределительной функции анализа рассеивания веществ для временного промежутка, равного 15 лет, представлены графически на рис. 2, 3. В целом из приведенных графиков видно, что наиболее активный период воздействия накопителя с толщиной слоя 20 см равен первым пяти годам существования накопителя, при этом величина риска за первый год GR£ = 2,4 10^-5 соответствует половине общей величины риска GR£ = 4,9 10^-5.

Рис. 2. Динамика экологического риска при эксплуатации накопителя

Рис. 3. Определение полей индивидуального риска (×10^-5) при расчёте распределительной функции по методике ТОКСИ. Масштаб 1:8000

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