Расчет поля долгопериодных средних концентраций загрязняющих веществ от одиночного точечного источника

10.1.1. Для расчета поля концентраций ЗВ, осредненных за длительный (сезон, год) период времени (далее -долгопериодных средних концентраций) ЗВ от одиночного точечного источника выброса значения осредненной мощности M выброса ЗВ, осредненной скорости w₀ выхода ГВС из устья источника и ее вертикальной составляющей w_s, объемного расхода V₁ ГВС и осредненного перегрева ГВС относительно окружающего атмосферного воздуха ΔT определяются в соответствии с отраслевыми нормативно-техническими документами при их наличии. В случае отсутствия таких документов, практические расчеты по формулам раздела 10 в соответствующих отраслях могут быть произведены после их разработки.

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

10.1.2. Для источников выбросов с постоянными в течение рассматриваемого периода времени параметрами выброса долгопериодные средние приземные концентрации C ЗВ определяются по формуле

(10.1)

где r_p и φ_p- полярные координаты расчетной точки относительно источника выброса;

p₁(φ) - функция, характеризующая угловое распределение концентрации ЗВ, которая выражается через розу ветров для рассматриваемого периода времени (см. 10.2.2).

Функция C¢(r_p, φ_p) вычисляется по формуле

, (10.2)

где u - скорость ветра на уровне флюгера (z_ф = 10 м), м/с;

λ- безразмерный параметр, характеризующий условия турбулентного перемешивания, рассчитываемый по формуле (10.3);

p₂(u) и p₃(λ) - соответствующие рассматриваемому периоду времени плотности вероятностей параметров u и λ (см. 10.2);

q_i (r_p, u, j_p, l, H_e) - подынтегральная функция, формулы для определения которой с учетом влияния рельефа местности и застройки приведены в 10.5.

H_e - эффективная высота источника выброса, определяемая согласно 10.1.3, м.

, (10.3)

где K₁ - коэффициент вертикального турбулентного обмена на уровне z₁=1 м, м²/с, определяемый согласно «Руководству по теплобалансовым наблюдениям» Л., Гидрометеоиздат, 1977;

u₁ - скорость ветра на уровне z₁=1 м, м/с.

10.1.3 Эффективная высота источника H_e определяется по формуле

, (10.4)

где ΔH - начальный подъем факела.

Для источника выброса с круглым устьем значение ΔH устанавливается в зависимости от u, λ и параметров выброса, по которым рассчитываются вспомогательные величины F_m, м⁴/с², и F_b, м⁴/с³:

, (10.5а)

, (10.5б)

где T_a - средняя абсолютная температура атмосферного воздуха в Кельвинах за рассматриваемый период времени;

g=9,81 м/с².

Температура Т_а определяется по формуле T_a =273+t_a, где t_a, °С – средняя температура атмосферного воздуха за период времени, использованный при определении функций p₁(j), p₂(u) и p₃(l) (см. 10.2.1). При расчетах допускается использовать постоянное значение Т_а, равное 283 К.

При выполнении неравенства 0> ΔT³ -5 °С расчет ΔH производится для ΔT=0. При ΔT<-5 °С за рекомендациями по расчету следует обращаться в территориальные органы Росгидромета.

При λ ³0,02 начальный подъем факела определяется по формуле

, (10.6а)

где

. (10.6б)

При λ <0,02 значение ΔH вычисляется как минимальное из значений ΔH₁ и ΔH₂. При этом ΔH₂ определяется на основе следующего алгоритма:

а) Величина S¢ принимается равной 6,7·10^-4 c^-2 при 0,01£ λ<0,02 и 1,17·10^-3 c^-2при λ<0,01;

б) Для источников выбросов высотой H>10 м определяется скорость ветра u_H:

; (10.7)

при H<10 м полагается u_H = u;

в) Вычисляется значение ΔT_c по формуле

; (10.8а)

г) При ΔT> ΔT_c значение ΔH₂ определяется по формуле

; (10.8б)

д) При ΔT£ ΔT_c значение ΔH₂ находится по формуле

. (10.8в)

Если источник оборудован специальным укрытием или выброс из источника осуществляется в горизонтальном направлении, то в формулах (10.6б) и (10.8в) F_m = 0.

Для источников выбросов с прямоугольным устьем расчет ΔH проводится по приведенным выше формулам с использованием эффективного диаметра устья источника D_э, м, и соответствующего эффективного объемного расхода ГВС V_1э, м³/с, которые определяются по формулам (5.36-5.38).

10.1.4. Для источников выбросов, параметры выброса которых изменяются в течение рассматриваемого периода времени более, чем на 10%, расчет долгопериодных средних концентраций ЗВ производится по одному из алгоритмов, устанавливаемых в 10.1.4.1- 10.1.4.3.

10.1.4.1. Для источника, мощность выброса которого определяется метеорологическими параметрами u и λ (например, для пылящей поверхности), C¢(r_p,j) вычисляется по формуле

(10.9)

где F(u, λ) - функция, характеризующая изменчивость мощности выбросов в зависимости от метеорологических параметров u и λ по отношению к ее характерному (например, среднему) значению M.

10.1.4.2 Для источников, параметры выброса которых не зависят от метеорологических параметров u и λ, расчет проводится по формулам (10.1 - 10.8). При этом M, V₁, w_s и DT задаются средними за рассматриваемый период времени значениями (см. пункт 10.1).

10.1.4.3. В общем случае для источника с переменными параметрами выброса рассматриваемый период времени разбивается на непересекающиеся интервалы, каждый из которых соответствует изменяющимся в пределах ±10 % значениям этих параметров. Для каждого из указанных интервалов продолжительностью t_j с использованием соответствующих этому интервалу функций p₁(φ), p₂(u) и p₃(λ) (см. 10.2) по формулам (10.1 - 10.9) определяется средняя концентрация C_j(r_p, φ_p) ЗВ, а затем находится средняя за рассматриваемый период времени концентрация ЗВ по формуле

, (10.10)

где t_Σ - общая продолжительность рассматриваемого периода времени.

П р и м е ч а н и е - Формула (10.10) позволяет рассчитать экспозицию для рассматриваемого ЗВ как произведение C(r_p, φ_p) на общую продолжительность рассматриваемого периода времени.

10.1.5. Вычисление интегралов в формулах (10.2) и (10.9) проводится по области изменения u и λ в данной местности за рассматриваемый период времени. При этом допускается использование квадратурных формул, обеспечивающих во всех узлах расчетной сетки относительную погрешность не более 3 %.

10.1.6. Расчет долгопериодных средних концентраций ЗВ проводится с учетом их частичной трансформации (см. пункт 4.4 настоящих Методов). Пример расчета долгопериодных средних концентраций для окислов азота приведен в Приложении Д.

П р и м е ч а н и е – Для конкретной рассматриваемой территории и конкретных трансформирующихся ЗВ по данным систематических специализированных натурных измерений концентраций этих примесей в атмосферном воздухе, могут определяться экспериментальные значения соответствующих коэффициентов трансформации согласно разъяснениям Росгидромета.

10.2. Учет распределения метеорологических параметров при расчете долгопериодных средних концентраций загрязняющих веществ в атмосферном воздухе

10.2.1. При расчете долгопериодных средних концентраций ЗВ изменчивость во времени метеорологических параметров характеризуется функциями p₁(j), p₂(u) и p₃(λ).

При вычислении средних концентраций ЗВ для конкретного интервала времени используются выборочные оценки указанных функций по данным измерений в течение этого интервала времени.

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

Для приближенных оценок функций p₁(j), p₂(u) и p₃(l) допускается использование рядов наблюдений продолжительностью не менее трех лет.

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

10.2.2. Функция p₁(j) определяется по розе ветров P_j, соответствующей рассматриваемому интервалу времени. При этом повторяемость штилей P₀ исключается из рассмотрения, а повторяемости румбов P_j (j = 1,2,3 ..) нормируются по формуле

. (10.11)

Значения p₁(j) определяются таким образом, чтобы выполнялись условия

(10.12)

где d - полуширина румба.

Середине румба j соответствует отличающееся от направления ветра на 180^о направление факела j. Внутри каждого румба функция p₁(j) аппроксимируется полиномом второго порядка, коэффициенты которого находятся из формулы (10.12) и условий непрерывности p₁(j) при переходе через границы румбов.

При аппроксимации функции p₁(j) должно быть выполнено условие p₁(j) >0.

10.2.3. Выборочная оценка плотностей вероятности p₁(j) и p₂(u) проводится по данным регулярных наблюдений за направлением и скоростью ветра на уровне флюгера (около 10 м над подстилающей поверхностью), выполненных на репрезентативной для рассматриваемой местности метеостанции в соответствии с нормативно-методическими документами Росгидромета [3].

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

При расчетах долгопериодных средних концентраций ЗВ для условий сложного рельефа местности с перепадом высот более 250 м на 1 км значения функций p₁(j) и p₂(u) могут быть получены в установленном порядке в подведомственных учреждениях Росгидромета при предоставлении соответствующего требованиям пункта 7.1 настоящих Методов, картографического материала, характеризующего физико-географические условия (рельеф и др.) района расположения метеостанции, предприятия и расчетной области.

10.2.4. Для выборочной оценки плотности вероятностей p₃(l) используются данные регулярных градиентных наблюдений на теплобалансовых станциях, условия размещения которых соответствуют требованиям, изложенным в 10.2.3. Сроки этих наблюдений должны соответствовать требованиям нормативно-методических документов Росгидромета. Значенияl за указанные сроки наблюдений на теплобалансовых станциях, используемые при оценке p₃(l), определяются по формуле

, (10.13)

где K₁ - значение коэффициента вертикального турбулентного обмена за указанные сроки наблюдений на уровне z₁=1 м, определенное согласно Руководству по теплобалансовым наблюдениям [4], м²/с;

u₂ - измеренная в рассматриваемый срок наблюдения скорость ветра на уровне z=2 м, м/с.

При отсутствии данных градиентных наблюдений значения функции p₃(l) могут быть получены в подведомственных учреждениях Росгидромета в установленном порядке.