Методы, утвержденные данным документом, введены в действие с 1 января 2018 года.

X. Метод расчета долгопериодных средних концентраций ЗВ

в атмосферном воздухе

10.1. Расчет поля долгопериодных средних концентраций ЗВ от одиночного точечного источника выброса.

10.1.1. Для расчета поля долгопериодных средних концентраций от одиночного точечного источника выброса значения осредненной мощности M выброса ЗВ, осредненной скорости w0 выхода ГВС из устья источника и ее вертикальной составляющей ws, объемного расхода V1 ГВС и осредненного перегрева ГВС относительно окружающего атмосферного воздуха 00000361.wmz определяются в соответствии с методиками расчета выбросов вредных (загрязняющих) веществ в атмосферный воздух стационарными источниками <10>. В случае отсутствия в указанных методиках необходимых сведений, практические расчеты по формулам, приведенным в главе X настоящих Методов, в соответствующих отраслях не производятся.

--------------------------------

<10> Правила разработки и утверждения методик расчета выбросов вредных (загрязняющих) веществ в атмосферный воздух стационарными источниками, утвержденные постановлением Правительства Российской Федерации от 16.05.2016 N 422 (Собрание законодательства Российской Федерации, 2016, N 21, ст. 3018).

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

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

00000362.wmz (109)

где rр и 00000363.wmz - полярные координаты расчетной точки относительно источника выброса;

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

Функция 00000365.wmz вычисляется по формуле (110):

00000366.wmz, (110)

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

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

p2(u) и 00000368.wmz - соответствующие рассматриваемому периоду времени плотности вероятностей параметров u и 00000369.wmz (пункт 10.2 настоящих Методов);

00000370.wmz - подынтегральная функция, формулы для определения которой с учетом влияния рельефа местности и застройки приведены в пункте 10.5 настоящих Методов.

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

00000371.wmz, (111)

где K1 - коэффициент вертикального турбулентного обмена на уровне z1 = 1 м, м2/с;

u1 - скорость ветра на уровне z1 = 1 м, м/с.

Упрощенный метод расчета среднегодовых концентраций, не связанный с использованием указанных плотностей вероятностей, приведен в пункте 10.6 настоящих Методов.

10.1.3. Эффективная высота источника выброса Hе определяется по формуле (112):

00000372.wmz, (112)

где 00000373.wmz - начальный подъем факела, то есть эффективное изменение высоты источника выброса под влиянием начальной скорости и/или перегрева выбрасываемой из источника газовоздушной (пылегазовоздушной) смеси за счет скоростного напора и/или сил плавучести;

Для источника выброса с круглым устьем значение 00000374.wmz устанавливается в зависимости от u, 00000375.wmz и параметров выброса, по которым рассчитываются вспомогательные величины Fm, м42, и Fb, м43:

00000376.wmz, (113а)

00000377.wmz, (113б)

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

g = 9,81 м/с2.

Температура Tа определяется по формуле Tа = 273 + tа, где tа, °C - средняя температура атмосферного воздуха за период времени, использованный при определении функций 00000378.wmz, p2(u) и 00000379.wmz (подпункт 10.2.1 настоящих Методов). При расчетах допускается использовать постоянное значение Tа, равное 283 K.

При выполнении неравенства 00000380.wmz °C расчет 00000381.wmz производится для 00000382.wmz.

При 00000383.wmz начальный подъем факела определяется по формуле (114а):

00000384.wmz, (114а)

где

00000385.wmz. (114б)

При 00000386.wmz значение 00000387.wmz вычисляется как минимальное из значений 00000388.wmz и 00000389.wmz. При этом 00000390.wmz определяется на основе следующего алгоритма:

а) Величина S' принимается равной 6,7 · 10-4 с-2 при 00000391.wmz и 1,17 · 10-3 с-2 при 00000392.wmz;

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

uH = u · [0,6667 + 0,1448 · ln(H)]; (115)

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

в) Вычисляется значение 00000393.wmz по формуле (116а):

00000394.wmz; (116а)

г) При 00000395.wmz значение 00000396.wmz определяется по формуле (116б):

00000397.wmz; (116б)

д) При 00000398.wmz значение 00000399.wmz находится по формуле (116в):

00000400.wmz. (116в)

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

Для источников выбросов с прямоугольным устьем расчет 00000401.wmz проводится по приведенным выше формулам с использованием эффективного диаметра устья источника Dэ, м, и соответствующего эффективного объемного расхода ГВС V, м3/с, которые определяются по формулам (31) - (33).

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

10.1.4.1. Для источника выброса, мощность выброса которого определяется метеорологическими параметрами u и 00000402.wmz (например, для пылящей поверхности), 00000403.wmz вычисляется по формуле (117):

00000404.wmz (117)

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

10.1.4.2. Для источников, параметры выброса которых не зависят от метеорологических параметров u и 00000407.wmz, расчет проводится по формулам (109) - (116в). При этом М, V1, ws и 00000408.wmz задаются средними за рассматриваемый период времени значениями в соответствии с пунктом 10.1 настоящих Методов.

10.1.4.3. В общем случае для источника с переменными параметрами выброса рассматриваемый период времени разбивается на непересекающиеся интервалы, каждый из которых соответствует изменяющимся в пределах 00000409.wmz 10% значениям этих параметров. Для каждого из указанных интервалов продолжительностью 00000410.wmz с использованием соответствующих этому интервалу функций 00000411.wmz, p2(u) и 00000412.wmz (пункт 10.2 настоящих Методов) по формулам (109) - (117) определяется средняя концентрация 00000413.wmz ЗВ, а затем находится средняя за рассматриваемый период времени концентрация ЗВ по формуле (118):

00000414.wmz, (118)

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

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

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

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

10.2.1. При расчете долгопериодных средних концентраций ЗВ изменчивость во времени метеорологических параметров характеризуется функциями 00000417.wmz, p2(u) и 00000418.wmz.

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

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

Для приближенных оценок функций 00000419.wmz, p2(u) и 00000420.wmz допускается использование данных наблюдений продолжительностью не менее трех лет.

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

10.2.2. Функция 00000421.wmz определяется по розе ветров Pj, соответствующей рассматриваемому интервалу времени. При этом повторяемость штилей P0 исключается из рассмотрения, а повторяемости румбов Pj (j = 1, 2, 3..) нормируются по формуле (119):

00000422.wmz. (119)

Значения 00000423.wmz определяются таким образом, чтобы выполнялись условия:

00000424.wmz, (120)

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

Середине румба j соответствует отличающееся от направления ветра на 180° направление факела 00000426.wmz. Внутри каждого румба функция 00000427.wmz аппроксимируется полиномом второго порядка, коэффициенты которого находятся из формулы (120) и условий непрерывности 00000428.wmz при переходе через границы румбов.

При аппроксимации функции 00000429.wmz должно быть выполнено условие 00000430.wmz.

10.2.3. Выборочная оценка плотностей вероятности 00000431.wmz и p2(u) проводится по данным регулярных наблюдений за направлением и скоростью ветра на уровне флюгера (около 10 м над подстилающей поверхностью), выполненных на репрезентативной для рассматриваемой местности метеостанции.

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

10.2.4. Для выборочной оценки плотности вероятностей 00000432.wmz используются данные регулярных срочных наблюдений за основными метеорологическими элементами (температурой и влажностью атмосферного воздуха, скоростью ветра), производимых одновременно на нескольких высотах в приземном слое атмосферного воздуха (далее - градиентные наблюдения) на теплобалансовых станциях. При этом значения 00000433.wmz, используемые при оценке 00000434.wmz, определяются по формуле (121):

00000435.wmz, (121)

где K1 - значение коэффициента вертикального турбулентного обмена за указанные сроки наблюдений на уровне z1 = 1 м, м2/с;

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

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

10.3. Расчет долгопериодных средних концентраций ЗВ от группы точечных, линейных и площадных источников загрязнения атмосферного воздуха

10.3.1. Долгопериодная средняя концентрация ЗВ C при наличии N источников выбросов определяется как сумма долгопериодных средних концентраций ЗВ от отдельных источников:

C = C1 + C2 +... + CN, (122)

где C1, C2,..., CN - концентрации ЗВ соответственно от первого, второго, N-го источника в рассматриваемой расчетной точке.

Концентрации C1, C2,..., CN ЗВ могут соответствовать точечным, линейным и площадным источникам выбросов.

В случае, когда имеются неучтенные источники выброса того же ЗВ (глава XI настоящих Методов), в правую часть формулы (122) добавляется слагаемое C'Ф, характеризующее фоновое загрязнение атмосферного воздуха от неучтенных источников выброса.

10.3.2. Долгопериодная средняя концентрация Cl ЗВ от линейного источника выброса, расположенного вдоль отрезка l длиной L, рассчитывается по формуле (123):

00000437.wmz, (123)

где C(x, y) - концентрация ЗВ, создаваемая в расчетной точке (x, y) точечным источником выброса, находящимся в точке 00000438.wmz отрезка l, интеграл в формуле (123) вычисляется вдоль этого отрезка. Подынтегральная функция в формуле (123) вычисляется по формулам, приведенным в пункте 10.1 настоящих Методов.

При расчетах для аэрационного фонаря подынтегральная функция в формуле (123) рассчитывается с использованием суммарного выброса из фонаря и эффективного диаметра, определяемого по формуле (37).

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

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

Для линейного источника, мощность выброса от которого изменяется вдоль отрезка l, подынтегральная функция в формуле (123) умножается на функцию 00000439.wmz, характеризующую изменение мощности выброса вдоль рассматриваемого источника по отношению к ее характерному значению, применяемому при расчете C. Учет изменения мощности выброса осуществляется и конкретный вид функции 00000440.wmz задается при наличии утвержденной в порядке, установленном законодательством Российской Федерации, методики расчета выбросов ЗВ в атмосферный воздух стационарными источниками для рассматриваемого производства <11>. В случае отсутствия соответствующей методики, учет изменения мощности выброса не осуществляется.

--------------------------------

<11> Правила разработки и утверждения методик расчета выбросов вредных (загрязняющих) веществ в атмосферный воздух стационарными источниками, утвержденные постановлением Правительства Российской Федерации от 16.05.2016 N 422 (Собрание законодательства Российской Федерации, 2016, N 21, ст. 3018).

10.3.3. Долгопериодная средняя концентрация Cs (x, y) ЗВ от площадного источника выброса, занимающего область S площадью SП, рассчитывается по формуле (124):

00000441.wmz, (124)

где C(x, y) - осредненная за этот период времени концентрация ЗВ, создаваемая в расчетной точке (x, y) точечным источником выброса, находящимся в точке 00000442.wmz области SП, по которой вычисляется интеграл.

Типы источников выбросов, которые могут аппроксимироваться площадными источниками, перечислены в пункте 8.8 настоящих Методов.

Подынтегральная функция в формуле (124) вычисляется по формулам, приведенным в пункте 10.1 настоящих Методов, с использованием суммарного выброса от всего площадного источника.

Для площадного источника, мощность выброса которого изменяется от точки к точке, подынтегральную функцию в формуле (124) следует умножить на функцию 00000443.wmz, характеризующую изменение удельной мощности выброса в точках рассматриваемого источника по отношению к ее характерному значению, применяемому при расчете C. Учет изменения мощности выброса осуществляется и конкретный вид функции 00000444.wmz задается при наличии утвержденной в порядке, установленном законодательством Российской Федерации, методики расчета выбросов ЗВ в атмосферный воздух стационарными источниками для рассматриваемого производства <12>. В случае отсутствия соответствующей методики, учет изменения мощности выброса не осуществляется.

--------------------------------

<12> Правила разработки и утверждения методик расчета выбросов вредных (загрязняющих) веществ в атмосферный воздух стационарными источниками, утвержденные постановлением Правительства Российской Федерации от 16.05.2016 N 422 (Собрание законодательства Российской Федерации, 2016, N 21, ст. 3018).

Таким же образом с заменой в формуле (124) двойного интеграла на тройной (аналогично формуле (64) проводится расчет среднегодовых концентраций от объемного источника выброса.

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

10.3.4. Требования к шагам расчетной сетки и общему количеству ее узлов принимаются в соответствии с пунктом 8.10 настоящих Методов.

10.4. Расчет максимальных значений долгопериодных средних концентраций ЗВ в атмосферном воздухе.

10.4.1. Максимальные значения долгопериодных средних концентраций CMAX ЗВ вычисляются по формуле (125):

CMAX = (1 + VC) · C, (125)

где C - осредненные концентрации ЗВ, относящиеся к средним за рассматриваемый период времени метеорологическим условиям согласно пункту 10.2 настоящих Методов;

VC - соответствующий коэффициент вариации.

10.4.2 Значение VC определяется по формуле (126):

00000445.wmz, (126)

где 00000446.wmz - стандартное отклонение долгопериодных средних концентраций ЗВ.

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

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

При расчете среднегодовых концентраций ЗВ в случае отсутствия необходимой исходной информации для расчета стандартного отклонения в формуле (126) и данных мониторинга загрязнения атмосферного воздуха допускается определять CMAX по формуле (125) при VC = 0,5.

10.5. Расчетные формулы для определения подынтегральных функций

10.5.1. Подынтегральные функции 00000448.wmz в формулах (110) и (117) вычисляются с использованием вспомогательной функции G:

00000449.wmz, (127)

где 00000450.wmz;

00000451.wmz - определяемый по формуле (41) безразмерный коэффициент, учитывающий влияние рельефа местности;

Hе - определяется согласно подпункту 10.1.3 настоящих Методов;

00000452.wmz, (128а)

00000453.wmz, (128б)

00000454.wmz, (129)

00000455.wmz (130)

00000456.wmz, (131)

00000457.wmz, (132а)

00000458.wmz. (132б)

Значение h в формуле (129) находится из соотношения:

00000459.wmz, (133а)

00000460.wmz, (133б)

где tзв = 530 с.

10.5.2. Долгопериодная средняя концентрация нерастворимых ЗВ или концентрация растворимых ЗВ, осредненная за период времени, в который отсутствовали осадки, согласно подпункту 10.5.3 настоящих Методов, определяется с использованием подынтегральной функции q0 (то есть qi при i = 0):

q0 = G + G1 + G2 + G3 + G4 при Hе 00000461.wmz 10 h, (134а)

q0 = 0 при Hе > 10 h, (134б)

где вспомогательная функция G определена по формуле (127), а параметры G1 - G4 вычисляются по формуле (127) при замене Hе на 20 · h - Hе, 20 · h + Hе, 40 · h - Hе и 40 · h + Hе, соответственно.

10.5.3. Для периодов времени с осадками подынтегральная функция q1 (то есть qi при i = 1) вычисляется по формуле (135):

00000462.wmz. (135)

Значение Rр вычисляется по формуле (136):

00000463.wmz, (136)

где I - средняя интенсивность осадков за рассматриваемый период времени, мм/ч;

00000464.wmz, с-1 - коэффициент вымывания рассматриваемого ЗВ твердыми либо жидкими осадками интенсивностью I0 = 1 мм/ч.

Коэффициент вымывания 00000465.wmz зависит от физико-химических свойств и дисперсного состава вымываемого ЗВ. Для мелкодисперсного аэрозоля диаметром не более 10 мкм и для диоксида серы в расчетах коэффициент вымывания следует принимать равным 1,3 · 10-5 с-1. По остальным ЗВ 00000466.wmz устанавливаются с использованием данных натурных и лабораторных экспериментов по определению характеристик вымывания этих ЗВ в атмосфере.

В случае, когда значения 00000467.wmz не установлены, для растворимых ЗВ допускается вычисление q1 по формулам (137а) - (137б):

q1 = q0 при I 00000468.wmz Iв, (137а)

q1 = 0 при I > Iв, (137б)

где Iв = 0,17 мм/ч.

При этом для растворимых ЗВ функции распределения 00000469.wmz, p2(u) и 00000470.wmz оцениваются для части интервала времени, соответствующей отсутствию осадков, а долгопериодные средние концентрации ЗВ вычисляются по формуле (118), в которой при наличии осадков принимается 00000471.wmz.

10.5.4. Влияние рельефа местности на поле долгопериодных средних концентраций ЗВ учитывается с помощью безразмерного коэффициента влияния рельефа 00000472.wmz. Этот коэффициент определяется согласно главе VII настоящих Методов для сечений рельефа местности, соответствующих 8 направлениям оси факела вдоль середины румбов, которые использованы при определении функции 00000473.wmz по формулам, приведенным в пункте 10.2 настоящих Методов. Для промежуточных направлений значение 00000474.wmz устанавливается с помощью линейной интерполяции.

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

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

10.5.5. Влияние застройки учитывается в соответствии с главой IX настоящих Методов. При этом в качестве масштаба длины xм используется коэффициент rм0, вычисленный по формуле (131), а построение объединенных ветровых теней для группы зданий проводится с учетом не более одного определяющего здания с наветренной и подветренной стороны.

10.5.6. При необходимости проведения расчета распространения ЗВ от автомагистрали вспомогательную функцию G в формулах (127), (134а) следует умножить на коэффициент sН:

00000477.wmz, (138а)

00000478.wmz, (138б)

где

00000479.wmz. (139)

Использование формул (138а), (138б) и (139) соответствует аппроксимации автомагистрали в виде объемного источника выброса высотой 2 м.

10.5.7. Расчет по формулам (127) - (132б) при 00000480.wmz проводится для расстояний rр, удовлетворяющих условию (140а - 140б):

00000481.wmz, (140а)

где

00000482.wmz. (140б)

Для расстояний, не удовлетворяющих этому условию, в формулах (127) и (10.24а), (10.24б) принимается 00000483.wmz = 1.

10.5.8. При расчете долгопериодных средних концентраций мелкодисперсных аэрозолей, для которых в соответствии с пунктом 5.6 настоящих Методов принято F > 1, правую часть формулы (127) следует умножить на безразмерный коэффициент 00000484.wmz, который вычисляется по формуле (141):

00000485.wmz, (141)

где 00000486.wmz (с/м) - вспомогательный параметр, определяемый по таблице 1 в зависимости от параметра F, установленного согласно пункту 5.6 настоящих Методов;