Нормативно-допустимый сброс (НДС) – это масса вредных веществ в водах, максимально допустимая к сбросу в данном водном объекте на единицу времени. Установление таких нормативов призвано обеспечить контроль за загрязнением воды в выбранном пункте.

НДС разрабатываются во исполнение Закона Российской Федерации «Об охране окружающей среды», постановления Правительства Российской Федерации № 469 «О Порядке утверждения нормативов допустимых сбросов веществ и микроорганизмов в водные объекты для водопользователей» и в соответствии с «Методикой разработки нормативов допустимых сбросов веществ и микроорганизмов в водные объекты для водопользователей», утвержденных Приказом № 333 Министерства природных ресурсов Российской Федерации.

Величины НДС определяются исходя из нормативов качества воды водного объекта. Если нормативы качества воды в водных объектах не могут быть достигнуты из-за воздействия природных факторов, не поддающихся регулированию, то величины НДС определяются исходя из условий соблюдения в контрольном пункте сформировавшегося природного фонового качества воды [8].

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

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

Величины НДС определяются для всех категорий водопользователей как произведение максимального часового расхода сточных вод – Qст (м3) на допустимую концентрацию загрязняющего вещества Сндс (г/м3). При расчете условий сброса сточных вод сначала определяется значение Сндс, обеспечивающее нормативное качество воды в контрольных створах с учетом требований Методики, а затем определяется НДС согласно формуле

НДС= Q ст C ндс (13)

Если фоновая концентрация загрязняющего вещества в водном объекте не превышает ПДК, то Сндс определяется в зависимости от типа водного объекта и используются расчетные формулы, приведенные ниже.

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

Основная расчетная формула для определения Сндс без учета неконсервативности вещества имеет вид

C ндс =n( C пдк C ф )+ C ф , (14)

где: С ПДК - предельно допустимая концентрация загрязняющего вещества в воде водотока, г/м3;

С ф - концентрация загрязняющего вещества в водотоке (г/м3) выше выпуска сточных вод, определяемая в соответствии с действующими методическими документами по проведению расчетов фоновых концентраций химических веществ в воде водотоков; n – кратность общего разбавления сточных вод в водотоке, равная произведению кратности начального разбавления, n н на кратность основного разбавления n о ( основное разбавление, возникающее при перемещении воды от места выпуска к расчетно-му створу):

n= n н n 0 . (15)

С учетом неконсервативности загрязняющего вещества расчетная формула имеет вид:

C ндс =n( C пдк e kt C ф )+ C ф , (16)

где k – коэффициент неконсервативности органических веществ, показывающий скорость по-требления кислорода, зависящий от характера органических веществ, 1/сут; t – время добегания от места выпуска сточных вод до расчетного створа, сутки.

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

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

Начальное разбавление обусловлено различием скорости истечения сточной жидкости из оголовка выпуска (Vст) и скорости движения окружающей массы воды (Vр). Разность скоростей приводит к возникновению турбулентного струйного потока, присоединяющего массы воды из окружающего пространства. Это приводит к увеличению поперечного сечения и расхода струи и к постепенному уменьшению скорости течения в струе и концентрации загрязняющих веществ в ней. Интенсивность начального разбавления зависит от отношения средней скорости потока к скорости истечения (mv):

m v = V ср V ст . (17)

Начальное разбавление рекомендуется учитывать для следующих условий: а) при соотношении скоростей потока и выпуска Vст ≥ 4Vср, т.е. 0 < mv  0,25; б) при абсолютных скоростях истечения струй из выпуска, больших 2,0 м/с. При меньших скоростях расчет начального разбавления не производится [1].

Начальное разбавление рассчитывается для участка от места выпуска до сечения, где струя не может больше присоединять новые объемы чистой воды. Это сечение следует принимать расположенным там, где скорость на оси струи на 0,1-0,15 м/с превышает среднюю скорость потока или где соседние струи, выходящие из рассеивающего выпуска, начинают соприкасаться.

Определение кратности начального разбавления (nн) может быть выполнено по зависимости, полученной Н.Н. Лапшевым:

n н = 0,248 1 m v d ̃ ст 2 ( m v 2 +8,1 1 m v d ̃ ст 2 m v )f( h ср d пт ) (18)

Значение функции f( h ср d пт ) принимается равным единице, если диаметр загрязненного пятна dпт в граничном слое начального разбавления меньше средней глубины водотока в зоне начального разбавления при минимальном среднемесячном расходе лимитирующего сезона hср. В противном случае значение f( h ср d пт )  определяется по специальной номограмме (рис. 1).

Для рассеивающего  напорного  выпуска  расчет  осуществляется следующим образом. Задаваясь числом выпускных отверстий оголовка выпуска N 0 и скоростью истечения сточных вод из них Vст ≥ 4Vср, или Vст ≥ 2 м/с, определяют диаметр отверстия или оголовка рассеивающего выпуска (d0) в метрах:

d 0 = 4 q ст π V ст N 0 (19)

где qст – суммарный расход сточных вод, м3/с; π  3,1416.

Относительный диаметр загрязненной струи в замыкающем створе зоны начального разбавления

( d ̃ ст )  может быть определена по формуле d ̃ ст = 7,46 Δ V m 2 ( 1 m v )+1,92Δ V m m v , (20)
где Δ V m = 0,125 V ст V ср . (21)

Диаметр загрязненного пятна dпт в граничном слое начального разбавления равен:

d пт = d ̃ ст d 0 . (22)

Расстояние до пограничного сечения зоны начального разбавления определяется по формуле

L н = d пт 0,48( 13,12 m v ) . (23)

Расход смеси сточных вод и воды водотока в том же сечении находится по формуле

Q см = n н q ст . (24)

Средняя концентрация вещества в пограничном сечении определяется по формуле

C ср,н = C ф + C ст C ф n н , (25)

где: Cст – концентрация нормируемого вещества в сточных водах, г/м3.

Максимальная концентрация в центре пятна примеси в этом сечении равна

C max,н = C ср,н 0,428 . (26)

Расчет основного разбавления по методу В.А. Фролова – И.Д. Родзиллера. Одним из широко употребляемых упрощенных методов расчета перемешивания сточных вод с водами речного потока является метод, предложенный в 1950 г. В.А. Фроловым и существенно дополненный и уточненный И.Д. Родзиллером. Метод может применяться для относительно прямых водотоков при соблюдении следующего неравенства [6]:

0,0025 q ст Q р 0,1 , (27)

где Qр – расчетный расход воды в реке на участке сброса сточных вод.

Кратность основного разбавления n0 определяется по формуле

n 0 = q ст +γ Q р q ст , (28)

где γ – коэффициент смешения, показывающий, какая часть речного расхода смешивается со сточными водами в максимально загрязненной струе расчетного створа.

При этом коэффициент смешения γ рассчитывается по формуле

γ= 1β 1+ β Q р q ст   , где β=exp( α L 3 ) , (29)

где L – расстояние от места выпуска сточных вод до контрольного створа; α – коэффициент, учитывающий влияние гидравлических условий смешения и определяемый по формуле

α=kφ D c / q ст 3 , (30)

где k – коэффициент, зависящий от места впадения сточных вод в реку (расcт. от берега), k = 1 расположении выпуска у берега, k = 1,5 при выпуске в стрежневой зоне русла реки; φ – коэффициент извилистости русла реки; Dc – коэффициент турбулентной диффузии. Применительно к рассматриваемому методу, коэффициент турбулентной диффузии Dc для летнего периода времени рассчитывают по формуле

D c = g V ср h ср 37 n ш C ш 2 , (31)

где g – ускорение свободного падения (g = 9,81 м/с2); Vср – средняя скорость водного потока; hср – средняя глубина потока; nш – коэффициент шероховатости ложа реки, определяемый по табл. 7п (по М.Ф. Срибному); Cш – коэффициент Шези (м0,5/с), определяемый по формуле Н.Н. Павловского (при hср ≤ 5 м и допущении равенства значений гидравлического радиуса и средней глубины):

C ш = h ср y п n ш , (32)

где уп – показатель степени, определяемый по формуле

y п =2,5 n ш 0,130,75 h ср ( n ш 0,1 ) . ; (33)

В случае проведения расчетов в зимний период (период ледостава) в формулы (31–33) вместо глубины потока hср вводится значение 0,5hср, а вместо коэффициента шероховатости ложа nш – его приведенное значение nшпр:

n шпр = n ш [ 1+ ( n л n ш ) 1,5 ] 0,67 , (34)

где nл – коэффициент шероховатости нижней поверхности льда по П.Н. Белоконю (табл. 8п).

Если также производится расчет начального разбавления по методу Н.Н. Лапшева, то расчетные формулы (28-30) имеют следующий вид:

n 0 = Q см +γ( Q р Q см + q ст ) Q см ; γ= 1β 1+ β Q р Q см ; α=kφ D c / Q см 3 , (35, 36, 37)

где Qсм – расход смеси сточных вод и воды водного объекта в пограничном сечении зоны начального разбавления, найденный по формуле (24).

Значение концентрации нормируемого вещества в максимально загрязненной струе водного потока на заданном расстоянии от места выпуска Cmax,L (г/м3) находится по формуле

C max,L = C ст q ст + C ф γ Q р q ст +γ Q р . (38)

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

C max,L =( C ст q ст + C ф γ Q р q ст +γ Q р )exp( kt ) . (39)

Метод В.А. Фролова – И.Д. Родзиллера достаточно прост в применении и позволяет получить достоверное представление о потенциально возможном разбавлении сточных вод в стационарных, максимально неблагоприятных условиях, что и определяет целесообразность его использования для расчета допустимых концентраций загрязняющих веществ в сточных водах [6].

Расчет основного разбавления по экспресс-методу ГГИ. В Государственном гидрологическом институте (ГГИ) М.А. Бесценной был разработан ряд упрощенных методов расчета разбавления на основе уравнения турбулентной диффузии. В данном методе на основе графических построений выполнен анализ связи между интенсивностью снижения показателя разбавления вдоль потока и гидравлическими характеристиками последнего. Получена аналитическая зависимость между этими величинами, которая приводит к определению кратности основного n0 разбавления и максимальной концентрации Cmax,L на любом расстоянии от места выпуска сточных вод:

n 0 = 1 C пдк C ст +μ( L ) , (40)
C max,L = C n +μ( L ) C ст (41)

где Cn – средняя концентрация нормируемого вещества в речном потоке ниже места выпуска, определяемая по формуле

C n = C ф Q р + C ст q ст Q р + q ст , (42)

μ(L) – эмпирический коэффициент неравномерности распределения примеси в расчетном створе, расположенном на расстоянии L от места выпуска:

μ( L )= 0,14 q ст N B р h ср B р Lφ( Q р + q ст ) , (43)

где Bр – ширина русла реки в расчетном створе; ϕ – коэффициент извилистости речного русла; N – безразмерное характеристическое число турбулентного потока, определяемое по формуле

N= M C ш g , (44)

где M – параметр, являющийся функцией коэффициента Шези Cш, связанный с Cш следующей зависимостью:

M={ 0,7 C ш +6            при         10 C ш 60 48                  при              C ш >60 (45)

Если учитывать начальное разбавление, найденное по методу Н.Н. Лапшева, то рас-четные формулы (40 – 43) имеют следующий вид:

n 0 = 1 C пдк C max,н +μ( L ) ; C max,L = C n +μ( L ) C max,н  ; (46, 47)
μ( L )= 0,14 Q см N B р h ср B р Lφ( Q р + q ст ) (48)

Экспресс-метод ГГИ рассматривается как приближенный и может применяться для относительно прямых водотоков при соблюдении неравенства (27).

Расчет разбавления сточных вод в реках по методу А.В. Караушева. При необходимости более подробного изучения условий распространения вещества в исследуемом водотоке или несоблюдении условий аналитических методов (экспресс-метода ГГИ, метода В.А. Фролова – И.Д. Родзиллера) рекомендуется использовать численные методы решения уравнения турбулентной диффузии. В частности, в нормативных документах имеются ссылки на предложенный А.В. Караушевым конечностно-разностный метод решения для условий плоской задачи, который в дальнейшем был модифицирован М.И. Бесценной и Л.И. Фаустовой для учета поперечной циркуляции в потоке и его кинематической неоднородности [5].

Метод А.В. Караушева позволяет получить поле концентрации нормируемого вещества в пределах всей расчетной области, начиная от источника загрязнения до створа водопользования.

Уравнение турбулентной диффузии для условий плоской задачи (при пренебрежении поперечными скоростями и стационарности процесса во времени) в конечных разностях имеет вид

ΔC Δx = D c ' V ср × Δ 2 C Δ z 2 . (49)

При расчете речной поток в плане разбивается сеткой, каждая вертикальная линия которой отвечает определенному поперечному сечению и отстоящей от предыдущего и последующго сечения на длину ∆х. Расстояние между горизонтальными линиями (по ширине реки) равно ∆z. Каждой клетке присвоен свой индекс по соответствующим осям координат: по оси х - к, по оси z - m (рис. 2).

Расчетная зависимость, позволяющая вычислить распределение концентрации загрязняющих веществ по длине и ширине потока (плоская задача), записывается в виде

C m,k+1 =0,5( C m+1,k + C m1,k ) . (50)

При расчете по уравнению (50) вся изучаемая область потока делится на прямоугольные параллелепипеды, объемы которых равны ∆х × ∆z × hср, где hср - средняя глубина в рассматриваемой области. При пользовании формулой (50) предполагается, что уже вблизи от выпуска сточные воды равномерно распространяются по всей глубине реки. Расстояние между расчетными сечениями определяется по формуле

Δx= V ср Δ z 2 2 D c , (51)

где Dc – коэффициент турбулентной диффузии с учетом поперечной циркуляции в речном потоке и его кинематической неоднородности:

D c = k 0 D c ' , (52)

где Dc – коэффициент турбулентной диффузии для условий плоской задачи; k0 – корректирующий коэффициент к значению Dc.

По А.В. Караушеву, значение Dc может быть найдено по формуле

D c ' = g V ср h ср M C ш . (53)

При этом значение корректирующего коэффициента k0 может быть найдено из выражения

log k 0 =0,25 θ h ( 1+0,54 ε пц )+0,589 ε пц 0,356 , (54)

где θh – параметр, выражающий изменчивость средних по поперечным сечениям глубин на рассматриваемом участке потока [6]:

θ h = h max,ср h ср h ср , (55)

где hmax,ср – максимальная из средних глубин в поперечных сечениях.

Для больших и средних рек расчетные значения θh не должны превышать 0,6, поэтому при θh > 0,6 принимается, что θh = 0,6.

Значение параметра εпц, учитывающего влияние поперечной циркуляции на коэффициент турбулентной диффузии, можно найти по формуле

ε пц =1+ 0,0042 h ср R изл M C ш M C ш (56)

где Rизл – средний радиус излучин русла реки (м) на рассматриваемом участке.

Когда раствор загрязняющего вещества достигает граничных поверхностей потока (берегов), для расчета диффузии следует использовать соотношение, учитывающее отсутствие переноса через стенки потока:

( ΔC Δz ) берег =0 . (57)

Учет граничных условий осуществляется путем введения в расчеты экстраполяционных значений концентрации. Расчетная сетка и поле концентрации условно распространяются за ограничивающие поток поверхности. При этом экстраполяционное значение концентрации в клетке, примыкающей к внешней поверхности стенки Cэкстр (рис. 2), и значение концентрации C1 в клетке, примыкающей к внутренней поверхности стенки на том же поперечнике, должны удовлетворять условию (57), что возможно только в случае, если Сэкстр = С1.

Учет начальных условий осуществляется путем задания места выпуска сточных вод, их расхода (qст) и их концентрации (Cст). На плане реки обозначают место сброса и через него проводят начальный поперечник. Ниже по течению поток схематизируется и делится на расчетные клетки. Скорость поступления сточных вод (Vст), сбрасываемых в водный объект, принимается равной скорости течения реки [1].

Часть площади поперечного сечения реки, занятая сточными водами fст, вычисляется по следующей зависимости:

f ст = q ст V ср . (58)

Затем определяется ширина загрязненной струи (bст) в начальном створе:

b ст = f ст h ср . (59)

В соответствии с величиной bст назначается ширина расчетной клетки ∆z. Обычно принимается, что ∆z ≥ bст/2, однако, если значения ∆z оказываются очень большими, то их уменьшают так, чтобы выполнялось неравенство

Δz0,1 B р . (60)

Клетки, попадающие в струю притока сточных вод в начальном поперечнике, заполняются цифрами, выражающими начальную концентрацию сточных вод (Cст), остальные клетки — цифрами, выражающими фоновую концентрацию загрязняющего вещества в реке (Cф) (в частном случае это может быть нулевая концентрация).

Кратность общего разбавления сточных вод в расчетном створе на расстоянии L от места их выпуска в поток находится по формуле

n= C ст C ф C max,L C ф . (61)

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

Расчет разбавления сточных вод в озере (водохранилище) по методу М.А. Руффеля. При наличии в водоеме устойчивых ветровых течений для расчета кратности общего разбавления n может быть использован метод М.А. Руффеля. В расчетах по этому методу рассматриваются два случая: 1) выпуск в мелководную часть или в верхнюю треть глубины водоема, загрязненная струя распространяется вдоль берега под воздействием прямого поверхностного течения, имеющего одинаковое с ветром направление; 2) выпуск в нижнюю треть глубины водоема, загрязненная струя распространяется к береговой полосе против выпуска под воздействием донного компенсационного течения, имеющего направление, обратное направлению ветра [8].

Метод М.А. Руффеля имеет следующие ограничения: глубина зоны смешения не превышает 10 м, расстояние от выпуска до контрольного створа вдоль берега в первом случае не превышает 20 км, расстояние от выхода сточных вод до берега против выпускного оголовка во втором случае не превышает 0,5 км.

Кратность общего разбавления определяется по формуле (15). Кратность начального разбавления вычисляется следующим образом:

  1. при выпуске в мелководье или верхнюю треть глубины водоема:
    n н = q ст +0,00215 V вт h ср.вд 2 q ст +0,000215 V вт h ср.вд 2 ; (62)
  2. при выпуске в нижнюю треть глубины водоема:
    n н = q ст +0,00158 V вт h ср.вд 2 q ст +0,000079 V вт h ср.вд 2 (63)
  3. где Vвт – скорость ветра над водой в месте выпуска сточных вод, м2/с; hср.вд – средняя глубина (м) участка водоема вблизи от места выпуска, при этом линейные размеры участка (Lз) зависят от средней глубины в озере (водоеме) h0 (табл. 9п).

Кратность основного разбавления вычисляется следующим образом:

  1. при выпуске в мелководье или в верхнюю треть глубины водоема:
    n 0 =1+0,412 ( l Δx ) 0,627+ 0,0002l Δx , (64)

    где: l – расстояние от места выпуска до контрольного створа, м;

    Δx=6,53 h ср.вд 1,167 ; (65)
  2. при выпуске в нижнюю треть глубины водоема:
    n 0 =1,85+2,32 ( l Δx ) 0,41+ 0,0064l Δx (66)
    Δx=4,41 h ср.вд 1,167 (67)

Величины НДС для выпусков сточных вод в водохранилища и озера определяются по расчетным формулам (13–16). Если не выполняются условия применимости метода М.А. Руффеля, то расчет кратности начального разбавления nн выполняется по методу Н.Н. Лапшева. Расчет кратности основного разбавления может быть выполнен численным методом А.В. Караушева.