РЕКОМЕНДАЦИИ

 

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

 

Р 52.24.627-2001

 

Дата введения 2002-07-01

 

Предисловие

 

1 РАЗРАБОТАНЫ Гидрохимическим институтом

 

2 РАЗРАБОТЧИКИ О.А. Клименко, канд. хим. наук; Н.М. Трунов, канд. техн. наук; В.Ф. Геков, канд. техн. наук; Р.С. Пятницына; О.В. Сергеева; В.А. Циркунов; Л.М. Федотова

 

3 УТВЕРЖДЕНЫ Федеральной службой России по гидрометеорологии и мониторингу окружающей среды 31.10.2001

 

4 ЗАРЕГИСТРИРОВАНЫ ЦКБ ГМП за номером Р 52.24.627-2001 от 09.11.2001

 

5 ВЗАМЕН "Усовершенствованных методических рекомендаций по оперативному прогнозированию распространения зон опасного аварийного загрязнения в водотоках и водоемах, а также уровней содержания в воде основных загрязняющих веществ" (СПб.: Гидрометеоиздат, 1992)

 

 

1 Область применения

 

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

Настоящие рекомендации предназначены для оперативно-производственных подразделений управлений по гидрометеорологии и мониторингу окружающей среды (УГМС) Федеральной службы России по метеорологии и мониторингу окружающей среды (Росгидромет), осуществляющих организацию и проведение наблюдений за состоянием поверхностных вод суши.

 

2 Нормативные ссылки

 

В настоящих рекомендациях использованы следующие стандарты:

ГОСТ 17.1.1.01-77 Охрана природы. Гидросфера. Использование и охрана воды. Основные термины и определения

ГОСТ 19179-73 Гидрология суши. Термины и определения

ГОСТ 19185-73 Гидротехника. Основные понятия. Термины и определения

ГОСТ 27065-86 (СТ СЭВ 5184-85) Качество вод. Термины и определения

 

3 Термины и определения

 

В настоящих рекомендациях использованы следующие термины и определения.

Аварийный сброс сточных вод - сброс сточных вод с превышением проектных или установленных предельно допустимых норм по расходу воды или по содержанию в ней одного или нескольких загрязняющих веществ

Вертикаль створа наблюдений - условная отвесная линия от поверхности воды (или льда) до дна в водоеме или водотоке, в пределах которой выполняют работы для получения данных о показателях состава и свойств воды

Водный объект - сосредоточение природных вод на поверхности суши либо в горных породах, имеющее характерные формы распространения и черты режима (ГОСТ 19179)

Водозабор - забор воды из водоема, водотока или подземного водоисточника (ГОСТ 19185)

Водопользование - использование водных объектов для удовлетворения любых нужд населения и народного хозяйства (ГОСТ 17.1.1.01)

Водопотребление - использование водных ресурсов с безвозвратным изъятием воды из водоисточника (ГОСТ 19185)

Водоток - водный объект, характеризующийся движением воды в направлении уклона в углублении земной поверхности (ГОСТ 19179)

Высокозагрязненные воды - воды с повышенным содержанием одного или нескольких загрязняющих веществ, исключающим или существенно ограничивающим водопользование на водном объекте

Горизонт створа наблюдений - место на вертикали (по глубине водотока или водоема), на котором производят работы для получения данных о показателях состава и свойств воды

Загрязняющее вещество - вещество в воде, вызывающее нарушение норм качества воды (ГОСТ 17.1.1.01)

Зона высокозагрязненных вод - участок водного объекта с высокозагрязненными водами

Качество воды - характеристика состава и свойств воды, определяющая пригодность ее для конкретных видов водопользования (ГОСТ 17.1.1.01)

Максимально загрязненная струя в створе водотока - масса воды с наиболее высоким содержанием загрязняющих веществ, занимающая определенную часть сечения водного потока

Самоочищение воды - совокупность природных процессов, направленных на восстановление экологического благополучия водных объектов (ГОСТ 27065)

Створ водотока (реки) - условное поперечное сечение водотока, используемое для оценок и прогноза качества воды

Сточные воды - воды, отводимые после использования в бытовой и производственной деятельности человека (ГОСТ 17.1.1.01)

 

4 Общие положения

 

4.1 Появление в речной сети в течение сравнительно непродолжительного времени высокозагрязненных вод в большинстве случаев связано с аварийным сбросом сточных вод.

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

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

 

Таблица 1 - Критерии высокого загрязнения речных вод

 

Вещества и показатели качества воды

Предельно допустимая концентрация (ПДК), мг/дм3

Уровень высокого загрязнения воды Свз, мг/дм3

Минерализация воды

1000

1400

1,4

Хлориды

300

450

1,5

Сульфаты

100

400

4,0

Железо общее

0,10

1,0

10

Марганец

0,01

0,10

10

Азот аммонийный

0,39

2,5

6,4

Азот нитратный

9,1

15,0

1,6

Фосфор минеральный

0,20

0,40

2,0

ХПК

30

150

5,0

БПК5

2,00

10

5,0

Ртуть

0,0001

0,001

10,0

Кадмий

0,001

0,020

20

Свинец

0,030

0,150

5,0

Мышьяк

0,050

0,250

5,0

Медь

0,001

0,030

30

Хром общий

0,030

0,200

6,6

Хром шестивалентный

0,020

0,100

5,0

Кобальт

0,010

0,100

10

Никель

0,010

0,100

10

Цинк

0,010

0,100

10

Цианиды

0,050

0,150

3

Фториды

0,75

1,5

2

СПАВ

0,10

1,5

15

Фенолы летучие

0,001

0,030

30

Нефтепродукты

0,05

1,5

30

Алюминий

0,04

1,00

25

Магний

40

120

3

 

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

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

В качестве контрольных следует выделять створы, расположенные:

- на 1 км выше мест основных водозаборов (перечень таких водозаборов заранее согласовывается);

- перед большими населенными пунктами (обычно эти створы совпадают со створами систематических гидрохимических наблюдений);

- в непосредственной близости к государственной границе;

- на расстоянии примерно суточного добегания водных масс до створов основных водозаборов (эти створы желательно совмещать со створами систематических гидрохимических наблюдений);

- в устьях рек на возможном пути перемещения зоны высокозагрязненных масс воды.

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

4.3 Для оперативного составления прогнозов в каждом УГМС по контролируемой территории следует иметь и при необходимости корректировать следующие материалы:

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

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

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

- сведения о гидрологическом режиме и морфометрических характеристиках русла водных объектов; для водотоков, где ведутся гидрологические наблюдения, в створах гидрологических постов целесообразно установить статистические зависимости значений максимальной umax и средней ucp по ширине реки скоростей течения, средних глубины H и ширины B реки, площади поперечного сечения реки F и коэффициента Шези c от уровня воды Hy или расхода речной воды Q;

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

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

4.4 В случае состоявшегося аварийного сброса сточных вод необходимо оперативно собрать информацию, включающую:

- время начала и продолжительность аварийного сброса сточных вод;

- перечень основных загрязняющих веществ в сточных водах на период аварийного сброса;

- режим аварийного сброса сточных вод (изменение расхода сточных вод и концентраций в них основных загрязняющих веществ за период аварийного сброса), а также расстояние от берега до места сброса сточных вод*;

- основные морфометрические и гидродинамические характеристики для характерных участков рек на период аварийного сброса.

_______________

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

 

4.5 К основным морфометрическим и гидродинамическим характеристикам водотоков следует относить:

- расход воды в речной сети на выделенных характерных участках ниже аварийного сброса сточных вод;

- средние и максимальные скорости течения воды, средние значения глубины и ширины речного потока на характерных участках, выделенных по морфометрическим характеристикам между контрольными створами, а также длину этих участков.

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

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

При отсутствии данных о расходе сточных вод, продолжительности аварийного сброса и концентрациях основных загрязняющих веществ в сточной воде или об одной из этих характеристик следует попытаться получить сведения о максимальных концентрациях основных загрязняющих веществ и объеме высокозагрязненных масс воды непосредственно на водном объекте в одном из нижележащих створов. В этом створе главная вертикаль должна располагаться в стрежне потока речных вод, дополнительные (одна или две) - на расстоянии примерно 1/3 от берега до стрежня потока со стороны аварийного сброса сточных вод и со стороны местоположения водопользователей ниже по течению реки (например, от берега, где расположен ближайший населенный пункт). Число горизонтов для отбора проб воды на вертикали определяется глубиной водотока: при глубине до 5 м устанавливают один горизонт (от 0,2 до 0,5 м от поверхности воды или нижней поверхности льда зимой); при глубине от 5 до 10 м - два горизонта (в 0,5 м от поверхности и в 0,5 м от дна); при глубине более 10 м - три горизонта (дополнительно берется промежуточный горизонт, расположенный посередине между поверхностью и дном). Для отбора одной пробы воды в каждой исследуемой точке сечения реки в течение 10 мин через равные промежутки времени следует отбирать от 10 до 15 примерно равных порций воды, смешивая их в одной емкости. Период наблюдения в нижележащем створе определяется временем прохождения всей зоны высокозагрязненных вод, примерное начало которого задается по результатам прогноза. При отсутствии прогностических данных отбор проб следует проводить через каждый час начиная с предполагаемого времени начала аварийной ситуации в интересуемом x-м створе tн,х. Последнее определяют по формуле

,                                                             (1)

где Lx - общая протяженность участка реки (по фарватеру) от места аварийного сброса сточных вод до заданного x-го створа (нижележащего), м;

umax - максимальная скорость течения воды в реке, м/с;

7200 - коэффициент перевода данных из секунд в часы.

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

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

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

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

В конце года или другого специально оцениваемого периода оправдываемость совокупности составленных за этот период прогнозов может быть рассчитана по формуле

,                                                             (2)

где h - оправдываемость совокупности прогнозов за рассматриваемый период, %;

nпр - число оправдавшихся прогнозов;

Nпр - общее число прогнозов.

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

4.9 По алгоритму, приведенному в приложении А, разработана программа для ПЭВМ "Зона ВЗ", которая является неотъемлемой частью настоящих рекомендаций. В приложении Б представлена инструкция для пользователей этой программы.

 

5 Теоретические основы проведения прогностических расчетов распространения по речной сети зон высокозагрязненных вод

 

5.1 Общие замечания

5.1.1 К основным характеристикам перемещающейся по речной сети зоны высокозагрязненных вод относят: время, через которое область зоны с максимальной концентрацией загрязняющих веществ достигает заданного контрольного створа реки ; максимальную концентрацию загрязняющего вещества в заданном контрольном створе реки при прохождении через него высокозагрязненных вод Cmax,x; время, через которое фронтальная и хвостовая части высокозагрязненных вод достигнут заданного контрольного створа реки tх,фр, tх,хв; продолжительность времени, в течение которого в заданном контрольном створе ожидаются концентрации загрязняющих веществ на уровне высокого загрязнения tх,з (рисунок 1).

 

 

Рисунок 1 - Общий вид распределения концентрации загрязняющего вещества в зоне высокозагрязненных вод

 

Наиболее важной из перечисленных характеристик является tх,фр.

5.1.2 Если при составлении прогноза необходимо рассматривать несколько загрязняющих веществ, то для определения времени перемещения фронтальной части зоны высокозагрязненных вод выбирают вещество, которому соответствует наименьшее значение tх,фр, а для определения времени перемещения хвостовой части зоны - вещество, которому соответствует наибольшее значение tх,хв. Значение параметра tх,з определяется в этом случае разницей рассчитанных значений tх,фр и tх,хв.

5.1.3 Особо следует выделять аварийные ситуации, связанные со сбросом в водный объект сырой нефти. Попав в водный объект, нефть не остается в неизменном виде. Из образовавшейся пленки, главным образом в течение первых пяти суток, происходит интенсивное испарение ее легких фракций, составляющее в зависимости от температуры воздуха от 14 до 30% от ее общей массы (коэффициент убыли нефтепродуктов примерно равен 0,03-0,06 1/сут) [1]. В результате быстрых процессов растворения, образования эмульсий и взвесей в водный поток из пленки первоначально переходит примерно от 10 до 20% нефти. Далее этот процесс продолжается постоянно (коэффициент перехода нефтепродуктов из пленки в воду примерно равен 0,15 1/сут, из воды в пленку - 0,04 1/сут). Под влиянием ветра и течений нефтяная пленка может перемещаться со скоростью до двух раз большей, чем верхний слой воды. Пригнанная к берегу, она оседает в прибрежной зоне или на частях растений [2]. При температуре воды 4 °С разложение нефтепродуктов в результате процессов биохимического и химического окисления практически не происходит. Начиная с 5 °С разложение нефтепродуктов ускоряется пропорционально повышению температуры [3].

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

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

 

5.2 Прогноз времени перемещения высокозагрязненных вод по речной сети между начальным и заданным створами

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

Как показали трассерные эксперименты, проведенные на разных реках сотрудниками Гидрохимического института (ГХИ), действительные значения скорости перемещения центра зоны высокозагрязненных вод лежат в пределах между средним и максимальным значениями скорости потока, причем на первых участках пути скорость перемещения центра зоны высокозагрязненных вод практически равна максимальной скорости течения воды в той части сечения реки, в которой находится основная масса этих вод. Учитывая этот факт, целесообразно представлять минимальное и максимальное время перемещения зоны высокозагрязненных вод, используя в расчетах соответственно значения максимальной umax и средней ucp скоростей течения речной воды.

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

,                                                                  (3)

где Lx - общая протяженность участка реки (по фарватеру) от места аварийного сброса сточных вод до заданного x-го створа, м;

u - соответственно рассматриваемому варианту расчета средняя или максимальная скорость течения воды на речном участке длиной Lx, м/с.

5.2.2 На речных участках с плавно изменяющимися морфометрическими характеристиками средние значения параметра  рассчитывают по формуле

,                                                             (4)

где uI и uII - соответственно варианту расчета средние или максимальные в сечении реки скорости течения воды в начале (индекс I) и конце (индекс II) рассматриваемого речного участка, м/с*.

_______________

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

 

5.2.3 На реке, где нет резко отличающихся по морфометрическим характеристикам участков, значения средней и максимальной скоростей течения воды между створами гидрологических постов (например, для створов водопользования) можно определить по интерполяционной формуле вида

,                                                       (5)

где ux - ожидаемая скорость течения речной воды в x-м створе, расположенном между первым I и вторым II гидрологическими постами или створами, где были проведены срочные измерения этого параметра, м/с;

LI/x - протяженность участка реки между створом первого гидрологического поста и заданным x-м створом, м;

LI/II - протяженность участка реки между створами гидрологических постов, м.

Если на рассматриваемом речном участке можно допустить линейное изменение значений по другим гидрологическим и морфометрическим характеристикам, то расчеты их средних значений можно также проводить по формулам (4) и (5).

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

5.2.4 При наличии выделенных участков реки (или рек), отличающихся по морфометрическим характеристикам, значение параметра tст определяют как сумму времени перемещения зоны высокозагрязненных вод на каждом из этих участков в виде

,                                                        (6)

где u1, u2, ..., un - используемые для расчета скорости течения речной воды на участках длиной L1, L2, ..., Ln.

5.2.5 Увеличение скорости течения речной воды в результате кратковременного сброса больших объемов сточных вод обычно имеет значение для расчетов только на первых километрах ниже места сброса. Чтобы обоснованно решить вопрос о необходимости учета этого фактора, следует рассчитать, насколько увеличилась средняя скорость течения речной воды в створе сброса сточных вод по формулам:

uн = (Q + qн) / Fн,                                                             (7)

uт = (Q + qт) / Fт,                                                            (8)

m = [(uт - uн) / uн] · 100 %,                                                       (9)

где uн и uт - ориентировочная средняя скорость речной воды соответственно до начала сброса и в момент максимального расхода аварийного сброса сточных вод, м/с;

Q - расход воды в реке выше места сброса сточных вод, м3/с;

qн и qт - расход сточных вод соответственно до начала и в момент максимального аварийного сброса, м3/с;

Fн и Fт - ориентировочная площадь поперечного сечения речного потока соответственно до начала и в момент максимального аварийного сброса сточных вод, м2;

m - коэффициент увеличения скорости течения речной воды, %.

Если повышение средней скорости течения речной воды составляет более 20%, то специалисты по гидрологическим прогнозам должны оценить длину участка реки, где будет происходить "распластывание" профиля распределения повышенного расхода речной воды и наблюдаться повышенная средняя скорость течения воды. В тех случаях, когда длина участка полного "распластывания" профиля распределения повышенного расхода воды составляет более 10% расстояния между местом сброса сточных вод и заданным для прогноза створом, при расчетах по формуле (6) следует вводить необходимые поправки в значения скоростей на выделенных речных участках (способы ориентировочного вычисления максимальных скоростей течения речной воды описаны в 5.2.6).

В простейшем случае, когда вода ниже места сброса сточных вод не выходит на пойму, примерную длину открытого (свободного ото льда) речного участка L* (м), где происходит полное "распластывание" профиля распределения повышенного расхода речной воды, можно определить по формуле Крицкого-Менкеля [4], решая ее относительно рассматриваемой длины участка реки:

,                                               (10)

где Qн и Qт - расход речной воды в створе сброса сточных вод соответственно до начала сброса и в период максимального расхода сбрасываемых сточных вод, м3/с;

I0 - уклон дна речного русла;

wст - объем сброшенных сточных вод, м3;

nш - коэффициент шероховатости русла, определяемый по таблицам 2, 3 или по формуле

,                                                             (11)

где dэ - эффективный диаметр донных отложений, соответствующий 50%-му значению крупности частиц по гранулометрической кривой, мм.

 

Таблица 2 - Значения коэффициента шероховатости nш для открытых русел (по М.Ф. Срибному)

 

Характеристика ложа

nш

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

0,025

Реки в благоприятных условиях течения

0,030

Реки в сравнительно благоприятных условиях, но с некоторым количеством камней и водорослей

0,035

Реки, имеющие сравнительно чистые русла, извилистые, с некоторыми неправильностями в направлении струй, или же прямые, но с неправильностями в рельефе дна (отмели, промоины, местами камни); некоторое увеличение количества водорослей

0,040

Русла больших и средних рек, значительно засоренные, извилистые и частично заросшие, каменистые, с неспокойным течением. Поймы больших и средних рек сравнительно разработанные, покрытые нормальным количеством растительности (травы, кустарник)

0,050

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

0,067

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

0,080

Горно-водопадного типа русла с крупновалунным строением ложа, перекаты ярко выражены, пенистость настолько сильна, что вода, потеряв прозрачность, имеет белый цвет; шум потока доминирует над всеми остальными звуками, делает разговор затруднительным

0,100

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

1,133

 

Таблица 3 - Коэффициент шероховатости nш для равнинных рек (по Б.В. Полякову)

 

Категория

Характеристика рек

nш

I

Реки с песчаным руслом, ровным, без растительности, с незначительным перемещением донных наносов

0,020-0,023

II

Реки с песчаным извилистым руслом, с большими перемещениями донных наносов, пойма, заросшая травой

0,023-0,033

III

Пойма, заросшая кустарником или редким лесом

0,033-0,045

IV

Пойма, заросшая лесом

0,045-0,060

 

На речном участке длиной L*, где повышен расход воды, средняя скорость течения воды

                                                            (12)

где uн и uк - соответственно варианту расчета средние или максимальные скорости течения в начале и конце участка.

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

,                                                            (13)

,                                                           (14)

где K1, K2 - переходные коэффициенты, приближенные значения которых приведены в таблице 4;

uср(пов) и umax(пов) - соответственно средняя и максимальная по ширине реки скорости поверхностного течения воды, м/с.

 

Таблица 4 - Данные о переходных коэффициентах K1, K2 [5]

 

Характеристика русла и условий протекания рек

Переходный коэффициент

Средняя глубина, м

менее 1

от 1 до 5

более 5

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

K1

0,78-0,86

0,87-0,88

0,89-0,90

K2

0,55-0,67

0,68-0,77

0,78-0,79

Реки большие и средние с менее благоприятными условиями протекания (значительно засоренные, частично заросшие, извилистые, каменистые, с неспокойным течением)

K1

0,70-0,77

0,78-0,85

0,86-0,87

K2

0,43-0,54

0,55-0,65

0,66-0,70

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

K1

-

0,70-0,79

0,80-0,84

K2

-

0,43-0,60

0,61-0,66

 

5.2.7 В случае разлива сырой нефти ориентировочный интервал наименьшего времени перемещения между заданными створами центра пятна нефтяной пленки tст(пл) может быть определен по формуле (6) с использованием значений uср(пов) и umax(пов).

 

5.3 Ориентировочный прогноз основных характеристик перемещающейся по речной сети зоны высокозагрязненных вод

5.3.1 Состав и содержание рекомендованных для прогноза математических моделей ориентированы на разные уровни полноты исходной информации об условиях аварийного сброса и характеристиках речного потока.

5.3.2 Исходной информацией для прогноза являются зафиксированные в пределах небольшого периода времени (часы, несколько суток) в одном из створов контроля за качеством воды в реке неединичные значения концентраций загрязняющего вещества (или ряда загрязняющих веществ), превышающие уровень Свз. Общая блок-схема проведения расчетов характеристик зоны высокозагрязненных вод в контрольных створах водного объекта показана на рисунке 2.

 

 

Рисунок 2 - Общая блок-схема проведения прогностических расчетов характеристик зоны высокозагрязненных вод в контрольных створах водного объекта

 

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

,                                   (15)

где  - средняя концентрация вещества в водном потоке на расстоянии Lx через время t, г/м3;

mв - начальная масса сброшенного в водоток вещества, г;

F - площадь поперечного сечения водного потока, м2;

Dx - коэффициент продольной дисперсии, м2/с;

t - время, прошедшее от начала сброса вещества в водоток, с;

K - коэффициент скорости самоочищения воды от сброшенного вещества, 1/с.

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

,                                                 (16)

где

;                              (17)

;                                                     (18)

CN,t - концентрация вещества в N-м сегменте распределения значений концентрации вещества в заданном контрольном створе с учетом формирования концентраций вещества во всех сегментах этого распределения значений концентрации, г/м3;

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

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

CN- концентрация загрязняющего вещества в N-м сегменте распределения значений концентрации в исходном створе реки, г/м3 (рисунок 3);

Dt0 - шаг по времени при разбиении распределения повышенных значений концентрации вещества на сегменты в рассматриваемой зоне высокозагрязненных вод, с;

N - номер сегмента в исходном распределении значений концентрации вещества, для которого ведется экстраполяционный расчет;

n - номер очередного сегмента с концентрацией CN, включенного в расчет CNn;

Nc - общее число сегментов, взятых для описания распределения значений концентрации вещества (число сегментов в начальном и конечном створах должно быть одинаковым).

 

 

а - исходная эпюра; б - вид эпюры при прохождении зоны через заданный x-й контрольный створ реки

t0 - период времени, в течение которого происходил сброс сточных вод с концентрацией загрязняющего вещества, превышающей Свз, или в течение которого в исходном створе водотока наблюдались концентрации загрязняющего вещества выше Свз, с

 

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

 

Совокупность точек CN,t позволяет воспроизвести все распределение значений концентрации вещества в заданном контрольном створе реки.

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

После введения такого коэффициента уравнение (17) будет иметь вид

,                     (19)

где

                                              (20)

a - коэффициент, корректирующий фронтальную часть распределения концентрации загрязняющего вещества в речном потоке;

b - коэффициент, корректирующий хвостовую часть этого распределения.

В случае небольших безвозвратных потерь загрязняющего вещества в слабопроточных зонах речного русла коэффициенты a и b являются взаимозависимыми величинами. Эту связь можно выразить уравнением вида

b = 0,31 a0,5.                                                                (21)

С целью учета индивидуальных особенностей речного потока коэффициенты a и b могут быть получены опытным путем, например в результате проведения трассерного эксперимента при разных расходах речной воды. Можно их определить и для момента конкретной аварийной ситуации, используя результаты измерения концентраций вещества в нижележащих контрольных створах реки (раздел 8).

Для реализации рассмотренной модели необходимо иметь значения , K и Dx.

Варианты расчета скоростей течения речной воды указаны в 5.2.1-5.2.7.

Коэффициенты скорости самоочищения загрязняющих веществ K желательно брать по результатам натурных наблюдений на рассматриваемом водном объекте. При отсутствии таких данных ориентировочные значения K можно брать из таблиц 5, 6. Для водотоков со скоростью течения воды более 0,02 м/с рекомендуется пользоваться ориентировочными коэффициентами K, приведенными в таблице 5. Для других загрязняющих веществ в аналогичных условиях можно использовать коэффициенты, помещенные в таблице 6, предварительно увеличив их в 3 раза. При отсутствии сведений о коэффициенте K принимают, что K = 0. Используя коэффициенты K для органических веществ, следует обращать внимание на наличие в зоне высокозагрязненных вод высоких значений концентраций тяжелых металлов или ядохимикатов. Если концентрации этих токсических веществ превышают значения уровня высокого загрязнения воды Свз, то первые 2-3 сут следует считать, что биохимического окисления органических веществ практически не происходит.

 

Таблица 5 - Ориентировочные значения коэффициентов скорости самоочищения речной воды от некоторых загрязняющих веществ K

 

Вещества и групповые показатели химического состава воды

Значения K (10-51/с)при температуре воды, °С

менее 10

от 10 до 15

более 15

Ионы аммония (NH)

1,04

2,08

3,12

Медь

0,69

1,38

2,08

БПК5

0,57

1,15

1,73

СПАВ

0,34

0,69

1,04

БПКполн

0,23

0,57

0,81

Фенолы

0,23

0,46

0,69

Цинк

0,11

0,34

0,69

ХПК

0,11

0,23

0,34

Нефтепродукты

0,11

0,23

0,34

Железо

0,11

0,23

0,34

Никель

0,11

0,23

0,34

Хром

0,11

0,23

0,34

Фосфорорганические пестициды

0,04

0,06

0,10

Хлорорганические пестициды

0,01

0,01

0,02


Таблица 6 - Ориентировочные значения коэффициентов скорости самоочищения воды водоемов от некоторых загрязняющих веществ K

 

Вещества и групповые показатели химического состава воды

Значения K (10-5 1/с)при температуре воды, °С

менее 10

от 10 до 15

более 15

Азот аммонийный

0,23

0,35

0,58

Аминофенол

0,06

0,17

0,23

Алкилсульфанат (АС) без наполнителя

0,17

0,58

0,81

Алкилсульфанат керосиновый

0,02

0,07

0,12

Ацетанилид

0,07

0,23

0,35

Бензальдегид

0,46

1,74

2,55

Бензин

0,02

0,06

0,09

БПК5

0,12

0,23

0,35

БПКполн

0,02

0,07

0,12

Гваякол

0,02

0,07

0,12

Гидрохинон

0,01

0,02

0,03

Глицерин

0,35

1,27

1,85

Дизельное топливо

0,02

0,06

0,09

Дисольван-4411

0,003

0,01

0,17

Капролактам

0,002

0,006

0,009

Карвакрол

0,009

0,03

0,05

Керосин

0,036

0,12

0,17

Кислота адипиновая

0,0323

0,12

0,12

Кислота бензинсульфоновая

0,05

0,17

0,23

Кислота бензойная

0,35

1,39

1,91

Кислота галловая

0,08

0,29

0,46

Кислота изофталевая

0,35

1,16

1,62

Кислота нафтеновая

0,07

0,23

0,35

Кислота сульфасалициловая

0,05

0,17

0,23

Кислота фталевая

0,35

1,10

1,56

м-Крезол

0,03

0,12

0,17

о-Крезол

0,03

0,12

0,17

n-Крезол

0,01

0,04

0,06

Крезол

0,03

0,12

0,17

Ксиленол

0,009

0,03

0,05

Мазут топочный

0,003

0,01

0,02

Мальтоза

0,12

0,35

0,58

Масло машинное

0,0

0,01

0,02

Масло солярное

0,02

0,06

0,09

Ментол

0,08

0,23

0,35

Метафос

0,003

0,01

0,02

Метилмеркаптан

0,009

0,03

0,05

Метол

0,003

0,01

0,02

a-Нафтол

0,02

0,06

0,09

b-Нафтол

0,17

0,06

0,80

Нефтепродукты

0,0

0,03

0,05

Нитробензальдегид

0,05

0,12

0,23

Нитробензин

0,06

0,23

0,29

Нитрофенол

0,06

0,17

0,23

Оксонол-17

0,07

0,23

0,35

Оп-10

0,001

0,003

0,006

"Ордрам"

0,003

0,012

0,18

Паста-АС-1

0,18

0,59

0,83

Паста-АС-2

0,12

0,47

0,65

Паста-ДНС-1

0,03

0,12

0,18

"Прогресс"

0,08

0,29

0,41

"Сатурн"

0,52

1,85

2,60

Синтанол ВТ-7

0,06

0,18

0,23

Синтанол ДС-10

0,03

0,12

0,18

СПАВ

0,06

0,12

0,18

Спирты (амиловый, гептиловый, изобутиловый метиловый, пропиловый, этиловый)

0,07

0,24

0,35

Сульфанол (всех марок)

0,003

0,012

0,017

Тимол

0,009

0,03

0,05

Тридекан

0,02

0,06

0,12

Уксусный альдегид

0,09

0,35

0,47

Фенол

0,07

0,24

0,35

Фенол (обобщенный показатель)

0,05

0,09

0,12

Формальдегид

0,23

0,81

1,16

Фосфорорганические пестициды (обобщенный показатель)

0,05

0,07

0,09

Фурфорол

0,09

0,29

0,41

Хлорорганические пестициды (обобщенный показатель)

0,01

0,01

0,02

ХПК

0,01

0,01

0,02

Хлорофенол

0,06

0,123

0,24

o-Хлоранилин

0,07

0,24

0,35

m-Хлоранилин

0,02

0,09

0,12

Целлюлоза

0,0006

0,002

0,002

"Ялан"

0,77

2,72

3,84

 

Коэффициенты продольной дисперсии Dx для водотоков шириной до 60-70 м рекомендуется определять в соответствии с работами [7, 8] по формуле*

,                                                  (22)

где u - соответственно варианту расчета максимальная или средняя по сечению речного потока скорость течения воды на рассматриваемом участке реки, м/с;

H, B - соответственно средние глубина и ширина речного потока, м;

c - коэффициент Шези, м0,5/с.

_______________

* В тексте приведены рабочие формулы расчета, их теоретический вид следующий:

; ,

где R - гидравлический радиус;

u - динамическая скорость.

 

Для водотоков с большей шириной русла

Dx = 43000 H u c-2,63.                                                        (23)

5.3.3 Исходной информацией для прогностических расчетов являются данные об аварийном сбросе сточных вод, который происходит в условиях одного из следующих режимов:

а) Сст = const, q = const;

б) Сст ¹ const, q = const;

в) Сст = const, q ¹ const;

г) Сст ¹ const, q ¹ const;

где Сст - концентрация загрязняющего вещества в сточных водах при аварийном сбросе, мг/дм3;

q - расход сточных вод в период аварийного сброса, м3/с.

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

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

(Cст,NСф(к)) yN,                                                          (24)

где Cст,N - концентрация загрязняющего вещества в N-м исходном достаточно малом сегменте (или периоде времени) распределения значений повышенных концентраций этого вещества в сточной воде в период аварийной ситуации, мг/дм3 (см. рисунок 3);

yN - коэффициент, характеризующий разбавление сточных вод в водном объекте (в долях единицы N-го расхода, выделяемого в пределах периода аварийного сброса).

Для учета процессов смешения и разбавления высокозагрязненных вод в реке за основу была взята математическая модель, разработанная Г.В. Филькиным [9, 10]. Эта модель представляет собой аналитическое решение уравнения, описывающего процесс турбулентной дисперсии для двумерной задачи, вида

,                                                          (25)

где Dy - коэффициент поперечной дисперсии;

C - концентрация загрязняющего вещества;

ucp - средняя скорость течения воды в водотоке;

x и y - соответственно продольная и поперечная координаты.

Рекомендуемое уравнение для нахождения yN имеет вид

         (26)

где

;                                            (27)

если , то принимается                                                                           (28)

если yN > 1, то принимается yN = 1;                                                                                                        (29)

y0 - координата места сброса сточных вод в сечении реки (отсчет координаты по ширине реки от берега, к которому ближе место сброса сточных вод), м;

qN - N-й измеренный или предложенный для прогностического расчета расход сточной воды, по времени соответствующий наличию в сточной воде концентрации, равной Сст,N, м3/с;

kв - коэффициент, ограничивающий значение ширины реки при расчетах смешения речных вод в случае кратковременного сброса. В формулах (26)-(28) используют средневзвешенные значения параметров B*, H*,  с учетом ряда речных участков, существенно отличающихся по морфометрическим характеристикам речного русла.

Для определения коэффициента поперечной дисперсии рекомендуется использовать следующие формулы:

а) при наличии сведений о радиусе кривизны русла Rизл (м), взятом как среднее значение для участка реки, расположенного ниже рассматриваемого места сброса сточных вод (или впадения притока) и включающего 1-2 излучины, и максимальной из средних по створам глубины на рассматриваемом участке Hmax*:

;                                                          (30)

_______________

* Формула (30) рекомендована в работе [11] для расчета разбавления сточных вод.

 

б) при отсутствии сведений о значениях Rизл и Hmax

.                                                           (31)

В формулах (30), (31) использованы обозначения:

g - ускорение свободного падения, равное 9,8 м/с2;

M - коэффициент, зависящий от коэффициента Шези c (при 10 < c < 60 м0,5M = 0,7c + 6; при c ³ 60 м0,5M = 48 = const);

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

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

Коэффициент Kобщ вычисляют по уравнению

lg Kобщ = 0,25 gн (1 + 0,54w) + 0,589w - 0,356,                                      (32)

где

,                                                   (33)

gн = (Hmax - H) / H                                                            (34)

5.3.4 Характеристики времени перемещения зоны высокозагрязненных вод (tх,фр, tх,хв, ) для рассмотренных двух случаев исходной информации об аварийной ситуации на водном объекте вычисляются следующим образом.

Время tх,фр, через которое фронтальная часть зоны высокозагрязненных вод достигнет заданного контрольного створа реки, можно определить по формуле

tх,фр = tст + (N0 - Nвз) Dt0.                                                  (35)

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

,                                                  (36)

,                                                  (37)

.                                                (38)

В формулах (35)-(38) использованы следующие обозначения:

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

 - ожидаемое по прогнозу время начиная с момента начала опасной аварийной ситуации, когда часть зоны с максимально высокой концентрацией вещества достигнет заданного контрольного створа реки, с;

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

Nвз - номер сегмента, в котором по результатам прогностических расчетов впервые отмечалась концентрация, равная или более Свз в продольном профиле распределения концентраций вещества в заданном контрольном створе;

 - номер сегмента, в котором по результатам прогностических расчетов в последний раз имела место концентрация вещества, равная или более Свз, в продольном профиле распределения концентраций вещества в заданном контрольном створе;

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

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

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

Cmax,x > Свз                                                                  (39)

Если условие (39) выполняется, то следует продолжать прогностические расчеты перемещения зоны высокозагрязненных вод по речной сети. Если условие (39) не выполняется, то их следует считать завершенными.

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

,                                                                (40)

,                                                                (41)

,                                                           (42)

,                                                               (43)

 = 0,                                                                    (44)

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

Вид расчета на следующем речном участке зависит от выполнения условия

                                                             (45)

Если условие (45) выполняется, то принимается, что имеет место попадание в реку условно нового аварийного сброса. Если условие (45) не выполняется, то следует считать, что происходит попадание относительно чистого притока в загрязненную (сточную) реку. Расчет в последнем случае коэффициента разбавления в максимально загрязненной струе yN(p) в заданном очередном контрольном створе реки следует выполнять по формуле

,                   (46)

где

;                                             (47)

B*, H*, , Dy - средневзвешенные значения параметров на участке реки от "узлового" створа до заданного очередного контрольного створа.

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

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

,                                       (48)

где Dx(max) - параметр Dx, полученный с использованием значений umax.

При использовании формулы (48) принимается, что в заданном контрольном створе в любой момент времени Cmax,x >> Свз.

 

6 Представление результатов прогностических расчетов

 

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

 

Таблица 7 - Результаты прогноза аварийной ситуации на водном объекте

 

р. Дон, 1 км выше г. Азова

(Местоположение контрольного створа)

 

38 км

 

ПО "Водоканал" г. Ростова-на-Дону

 

02.12.97; 9 ч 30 мин

(Расстояние по речной сети от места сброса до контрольного створа)

 

(Наименование источника аварийного сброса сточных вод)

 

(Сроки начала аварии:

число, месяц, год; часы, минуты)


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

Уровень высокого загрязнения воды Свз, мг/дм3

Учиты-

ваемая скорость течения воды в реке

Дата и время прибытия фронтальной части зоны высоко-

загрязненных вод tфр,х

Дата и время прибытия хвостовой части зоны высоко- загрязненных вод tхв,х

Продолжи-

тельность высокого загрязнения воды в створе tх,з, ч, мин

Макси-

мальная концен-

трация вещества Cmax,x мг/дм3

Дата и время начала появления максимальных концентраций вещества

БПК5

10,00

03.12.97;

3 ч 30 мин

03.12.97;

10 ч 10 мин

6 ч 40 мин

20,00

03.12.97;

4 ч 15 мин

 

 

03.12.97;

7 ч 20 мин

03.12.97;

14 ч 05 мин

6 ч 45 мин

19,00

03.12.97;

8 ч 15 мин

Нефтепродукты

1,50

03.12.97;

4 ч 10 мин

03.12.97;

9 ч 00 мин

5 ч 50 мин

2,25

03.12.97;

6 ч 20 мин

 

 

03.12.97;

6 ч 35 мин

03.12.97;

12 ч 35 мин

6 ч 00 мин

2,00

30.12.97;

9 ч 10 мин

 

7 Организация проведения трассерных экспериментов

 

7.1 Задачи и этапы трассерных экспериментов

Трассерные эксперименты рекомендуется использовать для решения следующих задач:

- оконтуривания возможного распространения зон высокого загрязнения речных вод путем имитации на водном объекте возможных аварийных ситуаций, в том числе от конкретных промышленных предприятий;

- уточнения или получения осредненных данных для конкретных участков рек морфометрических и гидродинамических характеристик, используемых в прогностических расчетах по математическим моделям;

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

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

- они должны быть безопасными в санитарно-гигиеническом и экологическом отношении;

- обладать устойчивостью в течение длительного времени в толще воды водного объекта и определяться с помощью несложных методов при наличии различных помех;

- иметь невысокую стоимость;

- обеспечивать их определение в воде при весьма малых концентрациях загрязняющих веществ;

- не должны изменять (ухудшать) физические свойства воды водного объекта.

Всем этим требованиям удовлетворяют флуоресцентные трассеры - маркеры, приготавливаемые из инертных полимерных материалов, не взаимодействующих с водной средой. По разработанной в ГХИ методике такие трассеры регистрируются в воде в условиях очень больших разбавлений (при значениях кратности разбавления до 1013).

Основными этапами выполнения трассерного эксперимента, имитирующего аварийную ситуацию, являются:

- изучение особенностей объекта исследований и разработка схемы эксперимента;

- приготовление индикаторов с заданными характеристиками;

- проведение гидрометрических работ;

- введение трассера в исследуемый водный объект по сценарию имитации аварийного сброса сточных вод;

- отбор проб воды по заданной схеме;

- количественное определение трассера в отобранных пробах и обработка полученной информации.

 

7.2 Изучение особенностей объекта исследований и разработка схемы эксперимента

Первый этап трасcерных экспериментов независимо от конкретных задач эксперимента включает детальное ознакомление с имеющимися данными, характеризующими объект исследования. Этот этап в случае необходимости включает: рекогносцировочное обследование участка проведения работ, в результате которого устанавливают ориентировочные скорости течения речных вод; определение расстояния между интересуемыми точками и створами; выбор створов, вертикалей и горизонтов отбора проб; составление схемы отбора проб и перемещений поперек и вдоль водного объекта; выявление возможных мешающих определению трассеров условий (наличие пленки нефтепродуктов, большое количество мелких водорослей или взвешенных веществ) и т.д. Чрезвычайно важной является и ориентировочная оценка (расчет) ожидаемых кратностей разбавления трассерного вещества в местах отбора проб. Последнее необходимо для того, чтобы при проведении трассерного эксперимента в отобранных пробах концентрации индикатора были достаточными для его надежного количественного определения.

 

7.3 Приготовление индикаторов с заданными характеристиками

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

_______________

* В настоящее время в ГХИ разработаны (и могут готовиться для заказчиков) два типа концентрированных препаратов, которые в расфасованном виде транспортируются на любые расстояния.

 

В практической работе хорошо зарекомендовали себя жидкие препараты в виде высококонцентрированных стабилизированных суспензий. Они дешевле, легче дозируются, разбавляются (распускаются) при приготовлении стартовых концентраций индикатора. Недостатком этих препаратов является дополнительная масса. Наиболее оптимальным является соотношение 10 л суспензии: 40 г сухого препарата трассера.

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

 

 

1 - емкость вместимостью 0,2 м3; 2 - патрубок; 3 - резиновый шланг; 4 - поплавок из пенопласта; 5 - воздушная трубка; 6, 7 - регулировочный винт; 8 - сменная шайба с калиброванным отверстием; 9 - уровень воды в емкости.

Рисунок 4 - Схема дозирующего устройства

 

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

 

7.4 Проведение гидрометрических работ

В объем гидрометрических работ входит:

- измерение значений расхода сбрасываемого в водный объект раствора трассера либо расхода потока сточных вод, в который произведено введение трассера, в течение всего периода имитации аварийного сброса;

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

 

7.5 Введение трассера в исследуемый водный объект

Режим введения трассера в исследуемый водный объект зависит от конкретной схемы эксперимента.

В случае имитации аварийного сброса сточных вод конкретного предприятия введение трассера проводят в сточные воды этого предприятия в режиме аварийной ситуации:

- выбирается несколько вариантов продолжительности имитируемого аварийного сброса (трассеры, отличающиеся по окраске, сбрасывают одновременно или по очереди, но с разной продолжительностью инжекции);

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

Имитацию аварийного сброса можно выполнять и путем инжекции трассера непосредственно в водный объект в заданной точке сечения водного потока.

Возможность использования достаточно больших стартовых концентраций трассера в сочетании с высокой чувствительностью метода его количественной регистрации в водном объекте позволяет применять несложные дозирующие устройства. Учитывая необходимость работы в различных режимах инжекции трассера и широкий диапазон пространственных масштабов водных объектов, где могут проводиться трассерные эксперименты, рекомендуется использовать универсальное дозирующее устройство (см. рисунок 4), позволяющее проводить введение трассера в заданном режиме практически в любой водный объект, в том числе непосредственно в поток сточных вод. Для работы дозирующего устройства не требуется энергопитания. Объем накопительной емкости дозирующего выпускного устройства составляет 0,2 м3, что является оптимальным для проведения трассерных исследований на самых разных водных объектах. Если по условиям моделирования требуются большие расходы сбрасываемых сточных вод, то в указанное оборудование должен быть включен электрический насос, с помощью которого можно подавать к смесителю дополнительные объемы воды из водного объекта.

В дозирующем устройстве с помощью сменной шайбы 8, имеющей калиброванное отверстие, можно регулировать напор воды, который определяется глубиной расположения водного отверстия шланга DH. При необходимости введения трассера в режиме многоточечной инжекции шланг 3 заменяют кабель-шлангом, состоящим из 10-20 тонких полихлорвиниловых трубок, по которым осуществляется независимое введение трассера в различные точки сечения потока. Устанавливая соответствующим образом шайбу 8, можно реализовать схему введения трассера с учетом расхода воды в каждой точке введения (режим взвешенной многоточечной инжекции).

Иногда по условию эксперимента требуется использование так называемого метода внедрения объема индикатора в исследуемый водный поток. В этом случае объем индикатора вводят (внедряют) в заданную точку сечения водного потока (имитация залпового сброса при транспортной аварии). Для реализации этого метода целесообразно использовать разработанное в ГХИ простое техническое приспособление, схема которого представлена на рисунке 5. С помощью двух штанг 1, одна из которых закреплена вертикально в носовой (кормовой) части лодки (катера) 9, системы растяжек 7 и груза 6 поплавок 3 с закрепленной на нем накопительной емкостью 2 устанавливают на нужной глубине водного объекта. Накопительную емкость изготавливают из полиэтиленового рукава необходимой вместимости. Верхний конец накопительной емкости закреплен на тяжелом металлическом кольце 5 и возвышается над поверхностью воды на 15-20 см. После того как система смонтирована, в накопительную емкость наливают индикатор, затем лодка (катер) вместе с оборудованием дрейфует к необходимой точке введения трассера. В нужный момент времени с помощью фиксатора 8 освобождают конец верхнего шнура, удерживающего тяжелое металлическое кольцо 5. Кольцо 5 падает на поверхность поплавка 3, увлекая за собой полиэтиленовый рукав. При этом стартовый объем трассера, имеющий первоначально форму цилиндра, утрачивает ограничивающую его боковую поверхность и становится участником того сложного движения водного потока, которое имеет место на данном участке водного объекта.

 

 

1 - разъемные штанги, удерживающие верхнее кольцо; 2 - накопительная емкость (полиэтиленовый рукав); 3 - поплавок; 4 - стартовый объем индикатора; 5 - верхнее кольцо;

6 - груз; 7 - растяжки (капроновые шнуры); 8 - фиксатор верхнего шнура; 9 - носовая (кормовая) часть лодки

Рисунок 5 - Схема введения индикатора в водный объект методом его внедрения

 

7.6 Отбор проб воды на водном объекте

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

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

На водном объекте используют те или иные разновидности двух основных схем отбора проб воды:

1) отбор проб в контрольных створах с целью наблюдения за основными характеристиками перемещающейся зоны загрязненных вод через эти створы;

2) отбор проб по сетке станций с целью наблюдения за пространственно-временной эволюцией шлейфа меченых трассером водных масс в водном объекте.

Первую схему отбора проб используют, как правило, на водотоках и проточных водоемах, а вторую - на малопроточных, непроточных водоемах или водоемах со сложной структурой течений.

Для реализации отбора проб воды по первой схеме используют технические средства одномоментного отбора проб во всех назначенных точках сечения потока (рисунок 6). Для этой цели в заданных створах поперек водотока натягивают капроновый трос, к которому крепят небольшие пенопластовые поплавки. Число поплавков и место их установки должны соответствовать числу и местоположению вертикалей отбора проб. На поплавки подвешивают гирлянды приспособлений для отбора проб, которые представляют собой отрезки капронового шнура длиной, несколько меньшей, чем глубина водотока на данной вертикали, с закрепленными на них емкостями для отбора проб воды. Верхний конец шнура гирлянды крепят на основном тросе, а к нижнему прикрепляют груз. В качестве емкостей для отбора проб используют полиэтиленовые бутыли вместимостью от 0,25 до 1 л. Число и местоположение бутылей в гирлянде должны соответствовать числу и положению точек отбора проб на данной вертикали. Бутыли для отбора закрывают специальными пробками со вставленными в них двумя стеклянными трубками 7 разной длины (см. рисунок 6). Изменяя диаметр трубки, через которую происходит забор воды, можно изменять время наполнения бутылей. Обычно это время выбирают равным не менее 2 мин.

 

 

Составные части устройства (внизу): 1 - плотик пенопластовый (500х500 мм);

2 - трос капроновый (диаметр сечения 6 мм); 3 - груз (массой 1-2 кг); 4 - емкость для пробы; 5 - кольцо по диаметру емкости; 6 - пластинка металлическая для крепления емкости на тросе;

7 - трубка для набора воды

Рисунок 6 - Схема размещения устройства для одномоментного отбора проб воды по сечению реки в контрольном створе (вверху)

 

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

По истечении времени наполнения бутылей первый отбор проб считается завершенным. Для следующего (очередного) отбора проб бутыли в гирляндах заменяются на новые. Гирлянды устанавливают на плотиках и отбор проб в заданное время повторяют.

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

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

Сроки отбора проб воды назначают по результатам прогностического расчета по математической модели. Начало отбора должно соответствовать минимальному времени прихода фронтальной части зоны высокозагрязненных вод (т.е. с учетом значений tх,фр(min)), окончание - максимальному времени прихода хвостовой части зоны загрязненных вод (с учетом значений tх,хв(max)). Для визуального ориентировочного контроля за зоной высокозагрязненных вод в месте введения трассера в водный объект в начале и конце имитируемого аварийного сброса можно бросить наборы поплавков, каждый из которых снабжен притопленным грузом для устойчивости поплавка против ветра. За время прохождения зоны высокозагрязненных вод через створ реки общее число отборов проб во времени должно составлять не менее 10-15. При наличии поплавков 3-4 отбора проб делают до их прихода и столько же после прихода поплавков, замыкающих аварийную ситуацию.

Для длительного хранения (более трех суток) пробы консервируют хлороформом (2-4 мл на 1 л воды).

 

7.7 Количественное определение трассера в отобранных пробах и обработка полученной информации

Количественное определение трассера в пробах воды основано на фильтрации измеренного объема воды через мембранный фильтр с размерами пор, меньшими, чем диаметры частиц трассера, и последующем подсчете числа частиц в поле зрения люминесцентного микроскопа (микроскопирование). Пробы фильтруют таким образом, чтобы исключить возможность загрязнения рабочей поверхности мембранного фильтра. Отобранную пробу воды тщательно взбалтывают и через мембранный фильтр с диаметром пор 0,45 мм фильтруют в колбу Бунзена определенный (фиксируемый) объем воды. Оптимальной для анализа является такая плотность числа частиц, когда на поверхности мембранного фильтра диаметром 12 мм содержится от 300 до 500 частиц. Если концентрация трассера в пробах мала, фильтрацию ведут до тех пор, пока процесс дальнейшей фильтрации не станет затруднительным вследствие кольматации фильтра. Залитый в воронку объем воды фильтруют до конца, и на этом процесс фильтрации заканчивают. Фильтровальную воронку осторожно, чтобы не загрязнить фильтр, развинчивают, и фильтр переносят на предметное стекло. Сверху фильтр плотно прижимают покровным стеклом, после чего предметное и покровное стекла фиксируют на торцах липкой лентой. Полученный препарат трассера маркируют фломастером, и в таком состоянии он может быть микроскопирован сразу или храниться неопределенное время. Фильтрат из колбы Бунзена переливают в мерный цилиндр, где измеряют его объем.

Микроскопирование препарата трассера на предметном стекле осуществляют с помощью люминесцентного микроскопа при возбуждении флуоресценции сине-фиолетовыми лучами (светоделительная пластина 360-440 нм со светофильтром ФС-1-4). В качестве запирающего светофильтра используют ЖС18+ЖС19. Как правило, подсчет частиц трассера проводят с помощью объектива 10х0,4 и окуляров 1,1х7. При необходимости возможно использование объектива 40х0,65Л, а также других окуляров. В поле зрения люминесцентного микроскопа частицы трассера легко идентифицируются даже при наличии посторонних частиц по яркости, цвету флуоресценции, а также по размерам и форме частиц.

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

На основе просчитывания полей зрения вычисляют среднее число частиц в поле зрения.

Число частиц N в единице объема (дм3) рассчитывают по формуле

,                                                             (48)

где S - площадь поверхности фильтра, см2;

V - объем профильтрованной воды, см3;

n - среднее число частиц в одном поле зрения;

a - площадь зрения, см2.

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


8 Использование результатов наблюдений на водном объекте для корректировки прогностических данных

 

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

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

tк,вtн,в ® tст(ут),                                                           (50)

,                                                              (51)

где uут - уточненная скорость перемещения загрязненных масс на рассматриваемом участке реки, м/с;

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

tст(ут) - уточненное время перемещения (добегания) загрязненных водных масс на рассматриваемом речном участке (берется как разница между значениями tк,в и tн,в в секундах);

L - длина рассматриваемого речного участка, м.

Если значения uут были получены при различных расходах речной воды, то целесообразно рассчитать уравнение регрессии для статистической связи типа uут = f (Q). Полученное в результате статистического анализа уравнение можно использовать для расчета параметра uут при прогнозе перемещения по реке высокозагрязненных вод при других гидрологических условиях.

Следующей важной характеристикой, для которой может оказаться целесообразной корректировка, является максимальная концентрация рассматриваемого загрязняющего вещества в зоне высокозагрязненных вод при прохождении ее через замыкающий интересуемый речной участок контрольный створ Cmax,x. Такую корректировку желательно провести, если полученное по заново проведенному прогнозу при ucp = umax = uут значение Cmax,x отличается от фактически зарегистрированного более чем на 5%.

Корректировка осуществляется подбором такого значения коэффициента самоочищения речной воды Kут, при котором разница между прогностическим и фактическим значениями Cmax,x будет практически несущественной. Для некоторых загрязняющих веществ коэффициент Kут может оказаться со знаком "минус".

После получения уточненных значений uут и Kут следует провести корректировку ожидаемого по прогнозу времени, когда фронтальная часть зоны высокозагрязненных вод достигнет замыкающего интересуемый участок реки контрольного створа tх,фр. Данную корректировку осуществляют путем прогонки программы расчета на ПЭВМ с подбором такого значения коэффициента aут, при котором рассчитываемое значение tх,фр при использовании полученных значений uут и Kут будет иметь минимальное отличие от фактического. Уточненное значение коэффициента aут следует подбирать в пределах 0,01-0,45. С увеличением значений aут рассчитываемые значения параметра tх,фр снижаются.

Аналогичную процедуру выполняют для корректировки ожидаемого по прогнозу времени, когда хвостовая часть зоны высокозагрязненных вод достигнет замыкающего интересуемый участок реки контрольного створа tх,хв. Корректировка проводится путем прогонки программы расчета на ПЭВМ с подбором такого значения коэффициента bут, при котором рассматриваемое значение tх,хв при использовании полученных значений uут, Kут, aут будет иметь минимальное отличие от фактического.

Уточненное значение коэффициента bут следует подбирать в пределах от 0,01 до 0,20. С увеличением значений bут рассчитываемые значения параметра tх,хв увеличиваются. Если в створе, замыкающем рассматриваемый участок реки, по результатам прогноза значения Cmax,x стали близки к уровню Свз и в распределении концентрации загрязняющего вещества в хвостовой его части наблюдается дополнительный пик (рост) значений концентрации, то для этого участка следует принять, что bут = 0.

Использование всех полученных уточненных параметров uут, Kут, aут, bут должно быть ограничено только рассмотренным конкретным участком реки в соответствующих корректировке гидрологических условиях.

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

При имитации аварийного сброса с ориентацией на конкретное, сравнительно устойчивое в водной среде загрязняющее вещество в качестве уровня Свз следует брать значение концентрации трассера Свз(тр), рассчитанное по формуле

,                                                            (52)

где

;                                                              (53)

Стр(н) - исходная концентрация трассера в сточных водах;

Сст(max) - максимальная концентрация интересуемого загрязняющего вещества в имитируемом аварийном сбросе сточных вод;

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

 

 

ПРИЛОЖЕНИЕ А

(обязательное)

 

Алгоритм решения задачи по прогнозу характеристик зон высокозагрязненных вод в заданных контрольных створах речной сети

 

А.1 Общие замечания

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

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

2) прогноз распространения зон высокозагрязненных вод по речной сети по результатам наблюдений, полученным в исходном контрольном створе реки;

3) прогноз распространения зон высокозагрязненных вод по речной сети при наличии сведений об изменениях расходов аварийного сброса сточных вод и содержании в них загрязняющих веществ.

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

 

А.2 Алгоритм решения задачи по первому варианту

А.2.1 Для решения первой задачи требуется следующая исходная информация:

- расстояние Lсб (м) от места (створа) аварийного сброса сточных вод до первого "узлового" контрольного створа на пути перемещения высокозагрязненных вод по речной сети;

- общая протяженность участка (участков) водного объекта Lx (м) от места (створа) аварийного сброса сточных вод до заданного x-го створа (створ, замыкающий интересуемый участок, не должен располагаться ниже "узлового" контрольного створа);

- время tн (число, месяц, часы, минуты), когда было зафиксировано начало аварийного сброса сточных вод;

- время tк (число, месяц, часы, минуты), когда было зафиксировано окончание аварийного сброса сточных вод (на момент начала прогноза данные о tк могут отсутствовать);

- протяженность Lj (м) выделенных в пределах расстояния Lx характерных по морфометрическим или другим характеристикам участков рек;

- средняя скорость течения воды ucp,j (м/с) на j-м выделенном характерном участке реки;

- максимальная скорость течения речной воды umax,j (м/с) на j-м выделенном характерном участке или коэффициент пересчета Ku,j значений ucp,j в umax,j;

- средняя глубина речного потока Hj (м) на j-м характерном участке;

- средняя ширина речного потока Bj (м) на j-м характерном участке;

- коэффициент шероховатости русла реки nш,j на j-м характерном участке;

- коэффициент шероховатости нижней поверхности льда nл,j (берется при наличии ледяного покрова по таблице А.1);

- число nj выделяемых характерных участков реки в пределах Lx.

А.2.2 Расчет максимальной скорости течения воды umax,j на j-м выделенном характерном участке, если в исходных данных отсутствует информация об этой характеристике, а известен коэффициент Ku,j, выполняют по формуле

.                                                            (А.1)

А.2.3 Расчет средневзвешенных значений средней и максимальной скоростей течения речной воды на заданном участке длиной Lx проводят по формулам

,                                                           (А.2)

.                                                        (А.3)

А.2.4 Расчет средневзвешенных значений всех приведенных в исходных данных морфометрических характеристик русла реки H*, B*, , , I осуществляют по формуле

,                                                           (А.4)

где  - средневзвешенное значение искомой характеристики;

Rj - значение характеристики на j-м выделенном характерном участке реки.

А.2.5 Расчет коэффициента Шези c выполняют по следующим формулам:

- при наличии в исходных данных информации о  (при наличии ледяного покрова на реке)

,                                                       (А.5)

где

,                                                  (А.6)

                                             (А.7)

- при отсутствии значений

.                                               (А.8)

 

Таблица А.1 - Значения коэффициента шероховатости нижней поверхности льда nл [12]

 

№ п/п

Число дней после ледостава

nл

1

До 10

0,150-0,050

2

10-20

0,100-0,040

3

21-60

0,050-0,030

4

61-80

0,040-0,015

5

81-110

0,025-0,010

Примечание - Для подпертых речных бьефов данные 1-й и 2-й строк, отвечающие рекам в бытовых условиях, следует уменьшить на 15%, 3-й и 4-й строк - на 35%. Меньшие значения nл характерны для гладкого ледяного покрова, большие - для ледяного покрова с торосами и шугой.

 

А.2.6 Расчет коэффициентов продольной дисперсии Dx проводят по следующим формулам:

- в случае, когда

B* > 70 м,                                                            (А.9)

,                                             (А.10)

,                                               (А.11)

где Dx(max) и Dx(min) - соответственно максимальное и минимальное значения коэффициента продольной дисперсии, м2/с;

- в противном случае

,                                      (А.12)

.                                         (А.13)

A.2.7 Расчет времени перемещения водных масс на речном участке между исходным и заданным контрольными створами tст выполняют по формулам

,                                                         (А.14)

,                                                        (А.15)

где tст(max) и tст(min) - соответственно максимальное и минимальное значения tст.

А.2.8 Расчет выходных данных проводится следующим образом.

А.2.8.1 Ожидаемое по прогнозу минимальное tх,фр(min) и максимальное tх,фр(max) время, когда фронтальная часть зоны высокозагрязненных вод достигнет заданного контрольного створа водного объекта, вычисляют по формулам

,                                  (А.16)

.                                   (А.17)

А.2.8.2 При наличии в исходных данных tк ожидаемое по прогнозу минимальное tх,хв(min) и максимальное tх,хв(max) время, когда хвостовая часть зоны высокозагрязненных вод достигнет контрольного створа водного объекта, считая от времени tк, определяют по формулам

;                                  (А.18)

.                                    (А.19)

A.2.9 Тип заданного контрольного створа используется в расчете следующим образом: если он является "узловым" и необходимо продолжить расчет для нижележащих створов, то полученные значения tх,фр и tх,хв переводят в абсолютное время и используют как новые характеристики аварийного сброса в этом "узловом" створе. Продолжение прогностического расчета выполняют в виде нового цикла по алгоритму первого варианта решения задачи.

 

А.3 Алгоритм решения задачи по второму варианту

А.3.1 Для решения второй задачи требуется следующая исходная информация:

- расстояние Lсб (м) от исходного до первого "узлового" контрольного створа на пути перемещения загрязненных вод по речной сети;

- общая протяженность участка (участков) водного объекта Lx (м) от исходного до заданного x-го контрольного створа (створ, замыкающий интересуемый участок, не должен располагаться ниже "узлового" створа), м;

- время tн (число, месяц, часы, минуты), когда в исходном створе было зафиксировано начало прохождения зоны высокозагрязненных вод (условное начало аварийного поступления загрязняющих веществ в водный объект);

- измеренные в исходном контрольном створе реки значения концентрации загрязняющего вещества CN (мг/дм3) во время прохождения через этот створ зоны высокозагрязненных вод;

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

- относительные значения времени tN (с), когда в исходном контрольном створе начиная с момента фиксирования высоких концентраций загрязняющего вещества было измерено соответствующее значение CN;

- продолжительность t0 (с) (к моменту начала прогностических расчетов) наличия в речной воде высоких концентраций загрязняющего вещества в исходном контрольном створе;

- примерная температура воды t (°C) в перемещающейся по реке зоне высокозагрязненных вод;

- используемое для прогностических расчетов значение коэффициента скорости самоочищения K (1/с) речных вод от загрязняющего вещества (в случае отсутствия данных принимается K = 0 1/с);

- уровень концентрации загрязняющего вещества Свз (мг/дм3), наличие которого в воде водного объекта расценивается как высокое загрязнение;

- ориентировочное время tз (с), в течение которого ниже исходного контрольного створа может не происходить самоочищение речной воды от рассматриваемого загрязняющего вещества (при отсутствии данных принимается tз = 0 с);

- протяженность Lj (м) выделенных в пределах расстояния Lx характерных по морфометрическим или другим характеристикам участков рек;

- средняя скорость течения воды ucp,j (м/с) на j-м выделенном характерном участке реки;

- средняя глубина речного потока Hj (м) на j-м выделенном характерном участке;

- средняя ширина речного потока Bj (м) на j-м выделенном характерном участке;

- коэффициент шероховатости русла реки nш,j на ее j-м выделенном участке;

- максимальная скорость течения речной воды umax,j (м/с) на j-м выделенном характерном участке или коэффициент пересчета Ku,j значений ucp,j в значения umax,j;

- коэффициент aj, корректирующий форму фронтальной части распределения высоких концентраций загрязняющих веществ в речной воде на j-м выделенном характерном участке реки (при отсутствии данных принимается aj = 0);

- коэффициент bj, корректирующий форму хвостовой части распределения концентраций загрязняющих веществ в речной воде на j-м выделенном характерном участке реки (при отсутствии данных принимается bj = 0,2);

- коэффициент j, характеризующий извилистость русла реки, представляющий собой отношение расстояния по фарватеру реки к расстоянию по прямой на участке реки в пределах смешения речных и сточных вод;

- коэффициент шероховатости нижней поверхности льда nл,j, представляемый в том случае, если река покрыта льдом (см. таблицу А.1);

- расход воды Qj3/с) в j-м контрольном створе реки.

А.3.2 Расчет максимальной скорости течения речной воды на j-м выделенном характерном участке, если в исходных данных отсутствует информация об этой характеристике, а приведен коэффициент Ku,j, выполняется по формуле (А.1).

А.3.3 Расчет средневзвешенного значения средней скорости течения речной воды  на заданном участке длиной Lx проводится по формуле (А.2).

А.3.4 Расчет средневзвешенного значения максимальной скорости течения речной воды  на заданном участке длиной Lx проводится по формуле (А.3).

А.3.5 Расчет средневзвешенных значений всех приведенных в исходных данных морфометрических характеристик русла реки H*, , Q*, B*, , , a*, b* на участке протяженностью Lx выполняется по формуле (А.4).

А.3.6 Расчет коэффициента Шези c рекомендуется проводить:

- при наличии ледяного покрова на реке (при наличии данных о ) по формулам (А.5)-(А.7);

- в противном случае (при отсутствии ледяного покрова на реке) по формуле (А.8).

А.3.7 Расчет коэффициентов продольной дисперсии Dx проводится с учетом условия (А.9) по формулам (А.10)-(А.13).

А.3.8 Расчет времени перемещения водных масс на речном участке между исходным и заданным контрольными створами tст осуществляют по формулам (А.14) и (А.15).

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

,                                                            (А.20)

,                                                      (А.21)

при этом если , то принимается, что ; если , то принимается, что .

После расчета числа основных сегментов Nc в зоне высоких концентраций вещества в исходном створе каждому из них присваивают порядковый номер и время, прошедшее от начала регистрации в этом створе концентраций на уровне Свз. Временные характеристики сегментов распределения концентраций вещества в исходном створе будут иметь вид NDt0, в заданном контрольном створе через время tст – (tст + NDt0), где N - порядковый номер сегмента.

А.3.10 Исходный продольный профиль распределения концентраций вещества обычно имеет малое число сегментов (по числу наблюдений за зоной высокозагрязненных вод в исходном контрольном створе), поэтому пересчет профиля распределения концентраций вещества в намеченные сегменты общим числом 50 и более выполняют методом интерполяции сплайнами второго порядка (описание использования такой интерполяции и библиографию по этому вопросу можно найти в работе [13]).

А.3.11 Сравнение времени прохождения зоны высокозагрязненных вод через исходный створ с временем добегания до заданного контрольного створа осуществляют в зависимости от варианта расчета с учетом средней или максимальной скорости течения воды по условию (А.22) или (А.23):

t0 ³ tст(min),                                                            (А.22)

t0 ³ tст(max).                                                           (А.23)

Если условие (А.22) или (А.23) выполняется, то переходят к А.3.13, в противном случае - к А.3.12.

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

А.3.12.1 Определение параметров  и :

,                                                   (А.24)

,                                                  (А.25)

где , - значения времени для вычисления изменения концентрации вещества в n-х сегментах, "участвующих" в формировании значений CN (соответственно в параллельных расчетах учитываются средняя и максимальная скорости возможного перемещения зоны загрязненных вод);

n - порядковый номер очередного используемого для расчета значений CN,n сегмента с концентрацией CN (вычисление номера такого сегмента возможно по обе "стороны" от главного сегмента, взятого для расчета CN,t: для n-х сегментов, имеющих большее время от начала аварии, чем у главного, порядковый номер берется со знаком "плюс" и, наоборот, для сегментов, имеющих меньшее время, - со знаком "минус"; номер главного сегмента в этой процедуре приравнивается к нулю).

А.3.12.2 Предварительное вычисление параметров  и :

,             (А.26)

.                 (А.27)

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

А.3.12.3 Предварительный расчет концентрации загрязняющего вещества в N-x сегментах в заданном контрольном створе реки:

,                                                  (А.28)

.                                                  (А.29)

Расчет параметров  и  по формулам (А.28) и (А.29) ведется для всех сегментов, в том числе для добавленных в контрольном створе. Значения CN,t < Свз, полученные на концах распределения концентраций вещества, в качестве результата расчета не учитывают.

А.3.12.4 Определение коэффициентов a* и b* проводят следующим образом.

А.3.12.4.1 Число основных сегментов во фронтальной части распределения высоких концентраций загрязняющего вещества Na, для которых при расчетах значений CN,n следует учитывать значения коэффициента a* при b* = 0, определяют по формуле

,                                                      (А.30)

где  - время начала появления в рассматриваемом контрольном створе концентраций загрязняющего вещества на уровне Cmax,x, с.

При Na £ 10 принимается a* = b* = 0.

А.3.12.4.2 Число сегментов в средней части распределения высоких концентраций загрязняющего вещества N0, для которых при расчетах параметра CN,n следует принимать a* = b* = 0, рассчитывают по формуле

,                                                    (А.31)

где  - время последней регистрации концентрации вещества на уровне Cmax,x в рассматриваемом контрольном створе.

А.3.12.4.3 Число сегментов в хвостовой части распределения высоких концентраций загрязняющего вещества Nb, для которых при расчетах параметра CN,n следует учитывать значения коэффициента b* при a* = 0, определяют по формуле

Nb = Nc - Na - N0.                                                        (А.32)

При Nb £ 10 принимается a* = b* = 0.

А.3.12.4.4 Определение значений коэффициента a*, необходимого для расчета CN,n, формирующих концентрацию CN в соответствующем N-м сегменте, проводят по формулам:

- для N = 1

,                                                            (А.33)

- для N в пределах от 2 до Na

,                                                    (А.34)

где  - исходное (максимальное) значение коэффициента a*.

А.3.12.4.5 Определение значений коэффициента b*, необходимого для расчета CN,n, формирующих концентрацию CN,t в соответствующем N-м сегменте хвостовой части распределения концентрации загрязняющего вещества, вычисляют по формулам:

- для N = Na + N0

b*(N) = 0;                                                              (А.35)

- для N в пределах от (Na + N0 + 1) до Nb

.                                             (А.36)

Если значение b было задано в исходных данных , то

.                                                    (А.37)

При корректировке исходных данных, чтобы сделать независимыми друг от друга вводимые значения коэффициентов a* и b*, вначале в расчет вводится значение a*, а затем значение b*.

А.3.12.5 Вычисление значений коэффициентов  и  проводят по формулам

                                (А.38)

                                   (А.39)

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

А.3.12.6 Для окончательного вычисления параметров  и  используют формулы

,     (А.40)

         (А.41)

А.3.12.7 Окончательный расчет концентраций загрязняющего вещества в N-х сегментах в заданном контрольном створе реки осуществляют по формулам (А.28) и (А.29), в которых используются параметры CN,n, полученные по формулам (А.40) и (А.41).

А.3.13 В случае, если время прохождения зоны высокозагрязненных вод через исходный створ t0 не менее времени перемещения водных масс между исходным и контрольным створами tст, расчет выполняют по формулам (А.24)-(А.41), при этом сегменты, для которых NDt0 ³ tст, в расчеты по формулам (А.24), (А.25) включать не следует.

А.3.14 Строят графики зависимостей

, ,

где tав - время, прошедшее от начала регистрации в исходном контрольном створе концентрации вещества на уровне Свз.

А.3.15 Осуществляют подготовку выходных данных:

tх,фр(min); tх,фр(max); tх,хв(min); tх,хв(max); ; ; ; ; ;

где характеристики с одним штрихом получены при использовании значений , с двумя штрихами - при использовании значений .

А.3.16 Если в очередном заданном контрольном створе оказалось, что Cmax,x £ Свз, то прогностические расчеты с использованием алгоритма второго варианта решения задачи заканчивают. Если заданный створ оказался "узловым" и в нем Cmax,x > Свз, то осуществляют преобразование данных, полученных в этом створе. При этом принимается

;                                                          (А.42)

;                                                          (А.43)

tх,з = t0;                                                               (А.44)

y0 = 0;                                                                (А.45)

Q = q.                                                                (А.46)

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

 

А.4 Алгоритм решения задачи по третьему варианту

А.4.1 Для решения третьей задачи требуется следующая исходная информация:

- расстояние Lсб (м) от "узлового" контрольного створа реки до места (створа) аварийного сброса сточных вод;

- общая протяженность участка (участков) водного объекта Lx (м) от места (створа) аварийного сброса сточных вод до заданного x-го створа (створ, замыкающий интересуемый участок, не должен располагаться ниже "узлового" створа);

- протяженность Lj (м) выделенных для расчетов на участке Lx характерных по морфометрическим или другим характеристикам участков рек;

- время tн (число, месяц, часы, минуты), когда было зафиксировано начало аварийного сброса сточных вод;

- продолжительность t0 (с) (к моменту начала прогностических расчетов) аварийного сброса сточных вод;

- максимальный расход сточных вод qmax3/с) в известный период аварийного сброса сточных вод;

- расстояние y0 (м) от места аварийного выпуска сточных вод до ближайшего берега (при сбросе с берега y0 = 0 м, при рассеивающем выпуске y0 равно расстоянию от берега до середины рассеивающего устройства);

- уровень концентрации рассматриваемого загрязняющего вещества Свз (мг/дм3), наличие которого в воде водного объекта расценивается как высокое загрязнение;

- концентрация загрязняющего вещества в заданном контрольном створе реки Сф(к) (мг/дм3) до прохождения через него зоны высокозагрязненных вод;

- параметры, характеризующие режим аварийного сброса интересуемого загрязняющего вещества (целесообразно представлять в виде таблицы А.2);

 

Таблица А.2 - Сведения о режиме аварийного сброса загрязняющего вещества

 

Порядковый номер проведенного измерения

, с

qN, м3/c

Cст,N, мг/дм3

 

 

 

 

Примечание - В таблицу включены следующие параметры: Cст,N, - N-я измеренная или предлагаемая для прогностического расчета концентрация загрязняющего вещества в аварийном сбросе сточных вод; qN  - N-й измеренный или предложенный для прогностического расчета расход воды в аварийном сбросе сточных вод;  - относительное время, когда в сточных водах, начиная с момента фиксации начала аварийного сброса, было измерено значение Cст,N или значение qN.

 

- используемый для прогностических расчетов коэффициент K (1/с) скорости самоочищения речных вод от загрязняющего вещества (в случае отсутствия данных принимается K = 0 1/с);

- примерная температура воды в перемещающейся по реке зоне высокозагрязненных вод t (°С);

- рекомендуемое для учета ориентировочное время tз (с), в течение которого ниже аварийного сброса сточных вод может не происходить самоочищение речной воды от рассматриваемого загрязняющего вещества (при отсутствии данных принимается tз = 0 с);

- средняя скорость течения воды ucp,j (м/с) на j-м выделенном характерном участке реки;

- максимальная скорость течения речной воды umax,j (м/с) на j-м выделенном характерном участке или коэффициент Ku,j пересчета значений ucp,j в umax,j;

- средняя глубина речного потока Hj (м) на j-м выделенном характерном участке;

- среднее значение максимальных глубин речного потока Hmax,j (м) на j-м выделенном характерном участке (при отсутствии данных не учитывается);

- средняя ширина речного потока Bj (м) на j-м выделенном характерном участке;

- коэффициент шероховатости nл,j нижней поверхности льда (см. таблицу А.1);

- коэффициент aj, корректирующий форму фронтальной части профиля распределения высоких концентраций загрязняющих веществ в речной воде на j-м выделенном характерном участке реки (при отсутствии данных принимается aj = 0);

- коэффициент bj, корректирующий форму хвостовой части профиля распределения концентраций загрязняющих веществ в речной воде на j-м выделенном характерном участке реки (при отсутствии данных принимается bj = 0,2);

- расход воды Qj3/с) в j-м контрольном створе;

- радиус кривизны русла Rизл (м), взятый как среднее значение для участка реки, расположенного непосредственно ниже места аварийного сброса сточных вод;

- коэффициент извилистости русла реки j, представляющий собой на участке смешения речных и сточных вод отношение расстояния, измеренного по фарватеру реки, к расстоянию, измеренному по прямой;

- коэффициент kз, используемый для выделения "активной" ширины реки в створе аварийного сброса сточных вод (ширины, где средние по вертикали скорости течения воды в речном потоке могут отличаться от максимальной не более чем на 40%; при отсутствии таких данных принимается kз = 0,7).

А.4.2 Расчет периода Dt0 и общего числа основных сегментов Nc, используемых для описания профиля распределения концентраций вещества в аварийном сбросе сточных вод выполняют по формулам (А.20), (А.21).

А.4.3 Исходный профиль распределения концентраций вещества обычно имеет малое число сегментов (по числу наблюдений за аварийным сбросом сточных вод), поэтому пересчет профиля распределения концентраций вещества в намеченные сегменты общим числом 50 и более выполняют методом интерполяции сплайнами второго порядка [13]. При этом каждому из сегментов присваивают последовательные номера N и относительное время, исходя из выражения NDt0, для заданного контрольного створа – (tст + NDt0).

А.4.4 Расчет максимальной скорости течения речной воды на j-м выделенном характерном участке umax,j, если в исходных данных отсутствует этот параметр, а приведен коэффициент Ku,j, выполняют по формуле (А.1).

А.4.5 Расчет средневзвешенного значения средней скорости течения речной воды до заданного x-го контрольного створа  проводят по формуле (А.2).

А.4.6 Расчет средневзвешенного значения максимальной скорости течения речной воды до заданного x-го створа  осуществляют по формуле (А.3).

А.4.7 Расчет средневзвешенных значений всех приведенных в исходных данных морфометрических характеристик русла реки H*, , Q*, B*, , , a*, b* на участке протяженностью Lx выполняют по формуле (А.4).

А.4.8 Расчет коэффициента Шези c проводят:

- при наличии ледяного покрова на реке (при наличии данных о ) по формулам (А.5)-(А. 7);

- при отсутствии ледяного покрова на реке по формуле (А.8)

А.4.9 Коэффициент поперечной дисперсии Dy рассчитывают по формулам:

а) при наличии данных о характеристиках Rизл и

;                                                       (А.47)

                                         (А.48)

;                                              (А.49)

;                                   (А.50)

Kобщ = 10z;                                                           (А.51)

                                                      (А.52)

где g = 9,8 м/с2;

б) при отсутствии данных о характеристиках Rизл и

                                                        (А.53)

А.4.10 Если условие

                                                          (А.54)

не выполняется, то осуществляется переход к А.4.11, в противном случае - переход к А.4.12.

А.4.11 Расчет коэффициента yN(p), учитывающего степень смешения речных и сточных вод для N-го сегмента, проводят по формулам

,                                              (А.55)

                 (А.56)

Если , то принимается

                                                        (А.57)

А.4.12 Если условие

                                                               (А.58)

выполняется, то осуществляется переход к А.4.13, в противном случае - переход к А.4.15.

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

,                                                           (А.59)

,                                                 (А.60)

                   (А.61)

(А.62)

, если (yБ - yС) £ 0,01                                           (А.63)

Условия вычисления : формулы (А.61) и (А.62) используют для подбора значения Lp, необходимого для условия (А.63). При вычислении Lp используют шаг DLp = 0,2B*. Вычисления по формулам (А.61) и (А.62) проводят до тех пор, пока не будет иметь место (yБ - yС) £ 0,01.

А.4.14 Откорректированное значение  получают с использованием формул:

если Lx £ , то

,                                                         (А.64)

если Lx > , то

,                      (А.65)

где L1 - протяженность первого ниже аварийного сброса выделенного характерного участка реки, м;

 - максимальная скорость течения воды в реке на первом ниже аварийного сброса выделенном характерном участке реки, м/с.

В дальнейшем принимают, что .

А.4.15 Расчет времени tст перемещения водных масс между створом, где произошел аварийный сброс сточных вод, и заданным контрольным створом осуществляют по формулам (А.14), (А.15).

А.4.16 Ориентировочное расстояние до створа практически полного смешения речных и сточных вод Lсм вычисляют по формуле

                         (А.66)

А.4.17 Расчет коэффициентов продольной дисперсии Dx проводят с учетом условия (А.9) по формулам (А.10)-(А.13).

А.4.18 Если выполняется условие

Lсм > t0 ,                                                           (А.67)

то осуществляют переход к А.4.19, в противном случае - переход к А.4.21 и принимают kв=1.

А.4.19 Расчет примерной длины зоны высокозагрязненных вод Lз:

.                                            (А.68)

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

,                                                 (А.69)

.                                                   (А.70)

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

,                                        (А.71)

             (А.72)

Обязательными условиями использования формулы (А.72) являются следующие: если yN > 1, то следует принимать yN = 1; если в результате расчета максимальное значение yN соответствует неравенству

                                                     (А.73)

то принимается

                                                     (А.74)

А.4.22 Сравнение продолжительности аварийного сброса сточных вод с временем добегания водных масс до заданного контрольного створа осуществляется по условиям

t0 ³ tст(min),                                                           (А.75)

t0 ³ tст(max).                                                           (А.76)

Если условие (А.75) или (А.76) выполняется, то следует перейти к А.4.24, в противном случае - к А.4.23.

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

А.4.23.1 Вычисление параметров  и  выполняют по формулам (А.24), (А. 25).

А.4.23.2 Расчет параметра Ap осуществляют по формулам:

а) при 2qN £  B* H*

Ap = (Cст,NСф(к)) yN;                                                   (А.77)

б) при 2qN >  B* H*

Ap = (Cст,NСф(к)) (1 - yN(p));                                             (А.78)

А.4.23.3 Параметры  и  предварительно рассчитывают по формулам (А.26) и (А.27), в которых выражение (CNСф(к)) заменено на параметр Ap (по А.4.23.2).

А.4.23.4 Значения параметров  и  предварительно вычисляют по формулам (А.28) и (А.29).

А.4.23.5 Коэффициенты a* и b* определяют в соответствии с указаниями в А.3.12.4 алгоритма решения задачи по второму варианту (по формулам (А.30)-(А.37)).

А.4.23.6 Коэффициенты  и  вычисляют в соответствии с указаниями в А.3.12.5 алгоритма решения задачи по второму варианту (по формулам (А.38), (А.39)).

А.4.23.7 Окончательный расчет параметров  и  выполняют по формулам (А.40), (А.41), в которых выражение (CNСф(к)) заменено на параметр Ap.

А.4.23.8 Окончательный расчет концентраций загрязняющего вещества в N-х сегментах в заданном контрольном створе реки осуществляют по формулам (А.28), (А.29), для которых используются значения параметра CN,n, полученные в соответствии с указаниями в А.4.23.7.

А.4.24 В случае, когда продолжительность аварийного сброса t0 превышает время перемещения водных масс от места сброса сточных вод до контрольного створа tст, расчет выполняют по формулам (А.24)-(А.41) с учетом указаний в А.4.23.4-А.4.23.7 и при условии, что при вычислении характеристик ,  по формулам (А.24), (А.25) сегменты, для которых NDt0 ³ tст, в расчете не участвуют.

А.4.25 Строят графики зависимостей

, ,

где tав - время, прошедшее от начала аварийного сброса сточных вод, с.

А.4.26 Осуществляют подготовку выходных данных:

tх,фр(min); tх,фр(max); tх,хв(min); tх,хв(max); ; ; ; ; ;

А.4.27 Если в очередном заданном контрольном створе оказалось, что Cmax,x £ Свз, то прогностические расчеты с использованием алгоритма третьего варианта решения задачи заканчивают. Если заданный створ оказался "узловым" и в нем Cmax,x > Свз, то осуществляют преобразование данных, полученных в этом створе. При этом принимается:

;                                                         (А.79)

;                                                        (А.80)

tх,з = t0;                                                               (А.81)

y0 = 0;                                                                (А.82)

qp,N = qN,                                                               (А.83)

где при наличии данных о yN

                                                            (А.84)

при наличии данных о yN(p)

                                                        (А.85)

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

 

 

Приложение Б

(обязательное)

 

Инструкция по использованию программы "Зона ВЗ" для проведения прогностических расчетов на ПЭВМ характеристик зон высокозагрязненных вод в водотоках

 

Б.1 Общие сведения

Программа "Зона ВЗ" предназначена для выполнения на ПЭВМ прогностических расчетов количественных характеристик перемещающейся по речной сети зоны высокозагрязненных вод в задаваемых контрольных створах.

Программа разработана в соответствии с настоящими рекомендациями в среде Delphi 5 и может работать на ПЭВМ IBM 486 и выше под управлением операционных систем Windows 95/98/NT.

 

Б.2 Выполнение программы "Зона ВЗ"

Запуск программы осуществляется в соответствии с правилами работы с приложениями в среде Windows.

Запуск сопровождается появлением на мониторе главной экранной формы программы, представленной на рисунке Б.1.

 

 

Рисунок Б.1 - Главная экранная форма программы "Зона ВЗ" с горизонтальным меню

 

В горизонтальном меню главной экранной формы предусмотрены кнопки вызова следующих видов работ:

"Створы" - ввод и корректировка необходимых для прогностических расчетов параметров, характеризующих место аварийного поступления в водоток загрязняющего вещества (исходный створ водотока), а также параметров, характеризующих речной поток на отдельных выделяемых характерных речных участках и условно относимых к замыкающим эти участки контрольным створам; открытие файлов с ранее подготовленными данными; запись текущих данных в отдельный файл;

"Параметры сброса" - ввод и корректировка данных, характеризующих условия аварийного поступления загрязняющего вещества в водоток; открытие файлов с ранее подготовленными данными; запись текущих данных в отдельный файл;

"Расчет" - расчет характеристик зоны высокозагрязненных вод при прохождении ее через заданные контрольные створы, замыкающие выделенные характерные речные участки;

"Просмотр" - просмотр полученных результатов текущего прогностического расчета или результатов ранее проведенного расчета, хранящихся в отдельном файле;

"Выход" - выход из программы.

 

Б.3 Формирование исходной информации

После нажатия в горизонтальном меню главной экранной формы кнопки "Створы" на экране появляется подменю с разделами: "Создание-правка"; "Открыть"; "Записать как..."

Активизация мышью раздела "Создание-правка" позволяет открыть экранную форму с редакционными окнами (рисунок Б.2), в которые нужно внести (или откорректировать) необходимую для расчета исходную информацию.

 

 

Рисунок Б.2 - Первая экранная форма раздела "Создание-правка" для внесения и корректировки информации

 

В данной и последующих экранных формах перемещаться между редакционными окнами можно с помощью мыши или клавиши "Tab". В окне раскрывающегося списка выбирается интересуемое загрязняющее вещество (активизируется мышью). Далее в соответствии с особенностями исходной информации переключателем на панели радиокнопок выбирается тип аварийной ситуации на водном объекте. Пункт "Аварийный сброс в водоток" выбирают в том случае, когда имеется информация о параметрах аварийного сброса сточных вод. Пункт "Обнаружена зона высокозагрязненных вод в водотоке" соответствует случаю, когда имеются данные о прохождении зоны высокозагрязненных вод через исходный створ наблюдения. Пункт "При неполной информации об аварии" выбирают при наличии информации только о времени и месте аварийного сброса сточных вод.

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

Назначение других кнопок:

"Уничтожить список створов" - очистка всех ранее заполненных редакционных окон и выход в главную экранную форму (начинать ввод новой информации по створам водотока рекомендуется после нажатия данной кнопки);

"Выход" - выход в главную экранную форму после внесения всей необходимой информации;

"След." - переход к следующей экранной форме для внесения или редактирования исходной информации;

"В конец" - переход к экранной форме с информацией по последнему из вводимых створов водотоков.

Окна "Всего створов" и "Текущий створ" являются информационными (не подлежат заполнению).

После нажатия кнопки "След." открывается новая форма с редакционными окнами, в которые заносится информация по следующему отдельному (очередному) участку водотока (рисунок Б.3).

 

 

Рисунок Б.3 - Вторая экранная форма раздела "Создание-правка" для внесения и корректировки информации

 

Занесение информации начинают с введения кода створа (по умолчанию порядковый номер створа присваивается автоматически). Затем уточняется (с помощью мыши или клавиатуры) местоположение контрольного створа: створ является "узловым", если он расположен в устье водотока или в месте впадения крупного притока, существенно меняющего ширину реки.

Окно "Максим. значен. средней глубины на участке" заполняется при наличии информации. Одно из окон "Гидравлич. уклон" или "Коэф. шероховатости" можно оставить незаполненным. При наличии других данных значения, приведенные в окнах "Коэф. извилистости", "Коэф., корректирующий фронтальную часть зоны" и "Коэф., корректирующий хвостовую часть зоны", могут быть заменены.

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

Из окон, в которые должна быть внесена необходимая для расчета информация, выйти нельзя, пока они не будут заполнены.

Назначение других кнопок:

"След." - переход к аналогичной экранной форме для корректировки или внесения данных по следующему отдельному участку (контрольному створу) водотока на пути перемещения зоны высокозагрязненных вод; если текущий створ был последним в списке введенных створов, то по нажатию этой клавиши будет открыта форма для создания нового створа в конце списка;

"Предыд." - переход к экранной форме, в которую были занесены данные по предыдущему участку (контрольному створу) водотока;

"В начало" - переход к экранной форме, с которой был начат цикл внесения информации по контрольным створам водотока;

"В конец" - переход к экранной форме, которой был закончен цикл внесения информации по контрольным створам водотока;

"Вставить створ" - выполнение вставки данных по новому (пропущенному) участку водотока в соответствующее ему место перед текущим в последовательности ранее занесенных данных по характерным речным участкам на пути перемещения зоны высокозагрязненных вод;

"Удалить створ" - удаление текущего (представленного на экранной форме) створа из списка контрольных створов;

"Уничтожить список створов" - уничтожение занесенной информации по всем контрольным створам;

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

Нажав кнопку "Створы" и затем активизировав раздел "Записать как...", можно записать занесенную перед этим информацию по контрольным створам в отдельный файл. Активизировав раздел "Открыть", можно открыть для работы ранее записанный файл с данными по участкам интересуемой системы водотоков.

Для перехода к записи данных о режиме аварийного сброса сточных вод в меню главной экранной формы предусмотрена кнопка "Параметры сброса". После ее нажатия в появившемся подменю активизируют раздел "Создание-правка". В открывшейся экранной форме (рисунок Б.4), используя редакционные окна, вносят результаты контроля за аварийным сбросом сточных вод. Вначале записывают (корректируют) время первого контрольного наблюдения за аварийным сбросом сточных вод (год, месяц, день, час, минуты) и результаты контроля: расход сточных вод и концентрацию загрязняющего вещества в них. Кнопка "Добавить" служит для занесения записанной информации в список исходных данных, характеризующих режим аварийного сброса. Далее таким же образом заносятся другие имеющиеся результаты контрольных наблюдений за аварийным сбросом сточных вод. Любой результат в списке может быть удален путем активизации нужной строки с последующим нажатием кнопки "Удалить". Вводимые в список результаты автоматически сортируются по времени контрольного наблюдения. Для очистки списка от всех записей служит кнопка "Уничтожить список".

 

 

Рисунок Б.4 - Экранная форма раздела "Параметры сброса" для внесения и корректировки информации

 

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

Для открытия файла с ранее занесенными данными о режиме аварийного сброса сточных вод в подменю вместо раздела "Записать как..." следует активизировать раздел "Открыть" и выбрать на жестком диске интересуемый файл.

 

Б.4 Прогностический расчет и просмотр результатов прогноза

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

Для рассмотрения полученных результатов в главном меню нажимается кнопка "Просмотр". В появившемся подменю активизируют раздел "Текущие результаты". В открывшейся экранной форме в окне "прокрутки" "Код, ингредиент, размерность" активизируют название загрязняющего вещества, а в окне "Номер створа" - номер интересуемого створа водотока (нулевой номер дает выход на исходный профиль распределения концентрации вещества в аварийном сбросе сточных вод). На экране появляется таблица с прогностическими значениями всех основных характеристик зоны высокозагрязненных вод при прохождении ее через заданный (запрашиваемый) створ водотока (рисунок Б.5). На графике проиллюстрированы результаты прогностических расчетов двух профилей распределений концентраций загрязняющего вещества при прохождении зоны высокозагрязненных вод через заданный створ: один из них рассчитан с учетом максимальной скорости течения речной воды, другой - с учетом средней скорости течения воды. Если с целью просмотра результатов расчета вместо раздела "Текущие результаты" активизировать раздел "Результаты из файла", то можно открыть и просмотреть интересуемые результаты другого прогностического расчета, сохраненные в отдельном файле. Кнопка "Выход" в итоговой и главной экранной формах служит для последовательного выхода из программы.

 

Рисунок Б.5 - Итоговая форма, иллюстрирующая вид представляемых результатов прогностического расчета

 

Б.5 Другие варианты прогностических расчетов

В случае если в основу расчета должна быть положена информация не о режиме аварийного сброса сточных вод, а о зоне высокозагрязненных вод при прохождении ее через исходный контрольный створ, в разделе "Создание-правка" на форме "Ввод параметров в створе аварийного сброса" выбирается пункт "Обнаружена зона высокозагрязненных вод в водотоке". Далее процедура внесения информации практически не отличается от вышерассмотренной. Особенность ее внесения заключается лишь в том, что значения в редакционных окнах "Радиус кривизны русла...", "Расстояние от берега до точки сброса...", "Коэф. доли ширины водотока..." не заполняют и не корректируют. Информация в них нужна только для условий аварийного сброса сточных вод.

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

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

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

Представленные значения времени характеризуют перемещение этой зоны с учетом максимальной и средней скоростей течения речной воды.

 

 

Приложение В

(справочное)

 

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

 

Пример В.1 - Исходные данные: тип аварийной ситуации - "При неполной информации об аварии"; tн - 07.07.2000 г. 11 ч 20 мин; основное загрязняющее вещество - СПАВ;

Lсб = 30 км; Lx = 30 км; L1 = 10 км; B1 = 40 м; H1 = 1,2 м; ucp,1 = 0,45 м/с; umax,1 = 0,60 м/с; Q1 = 21,6 м3/с; nш,1 = 0,02; j1 = 1,2; L2 = 20 км; B2 = 45 м; H2 = 1,3 м; ucp,2 = 0,50 м/с; umax,2 = 0,71 м/с; Q2 = 29,2 м3/с; nш,2 = 0,02; j2 = 1,2. Индексы "1" и "2" обозначают номер j-го участка реки.

Результаты прогностического расчета:

- минимальное ориентировочное время прихода фронтальной части зоны загрязненных вод: 07.07.2000 г. 21 ч 36 мин;

- максимальное ориентировочное время прихода фронтальной части зоны загрязненных вод: 08.07.2000 г. 1 ч 35 мин.

Пример В.2 - Исходные данные: тип аварийной ситуации "Обнаружена зона высокозагрязненных вод в водотоке"; основное загрязняющее вещество - фенолы;

Lсб = 80 км; Lx = 70 км; L1 = 40 км; B1 = 100 м; H1 = 1,82 м; ucp,1 = 0,25 м/с; umax,1 = 0,32 м/с; Q1 = 45,5 м3/с; I1 =0,012+; j1 = 1,1; L2 = 30 км; B2 = 100 м; H2 = 2,14 м; ucp,2 = 0,30 м/с; umax,2 = 0,36 м/с; Q2 = 64,2 м3/с; I2 = 0,01+; j2 = 1,1; Сф(к) = 0,001 мг/дм3; Свз = 0,03 мг/дм3.

Измеренные характеристики зоны высокозагрязненных вод в водотоке приведены в таблице В.1. Результаты прогноза представлены в таблице В.2.

 

Таблица B.1 - Измеренные характеристики зоны высокозагрязненных вод при прохождении через исходный контрольный створ

 

Номер проведенного измерения

Дата и время

Концентрация, мг/дм3

Расход воды в реке, м3/c

1

28.10.2000 г.

05 ч 30 мин

0,200

38

2

28.10.2000 г.

08 ч 00 мин

0,350

38

3

28.10.2000 г.

13 ч 00 мин

0,100

38

4

28.10.2000 г.

18 ч 00 мин

0,100

38

5

29.10.2000 г.

00 ч 00 мин

0,400

38

6

29.10.2000 г.

03 ч 00 мин

0,350

38

7

29.10.2000 г.

05 ч 00 мин

0,010

38

 

Таблица В.2 - Результаты прогностического расчета

 

Характеристика

При максимальной скорости течения воды

При средней скорости течения воды

Время появления фронтальной части зоны высокозагрязненных вод

30.10.2000 г.

14 ч 53 мин 38 с

31.10.2000 г.

5 ч 5 мин 13 с

Время появления хвостовой части зоны высокозагрязненных вод

31.10.2000 г.

15 ч 06 мин 30 с

01.11.2000 г.

5 ч 37 мин 9 с

Продолжительность времени прохождения зоны высокозагрязненных вод

1 д. 0 ч 12 мин

1 д. 0 ч 31 мин

Продолжительность времени прохождения части зоны с наибольшей концентрацией вещества

0 д. 2 ч 56 мин

0 д. 2 ч 47 мин

Время появления максимальной концентрации вещества

31.10.2000 г.

10 ч 58 мин 48 с

01.11.2000 г.

1 ч 19 мин 55 с

Значение максимальной концентрации вещества

0,39996 мг/дм3

0,39973 мг/дм3

 

Пример В.3 - Исходные данные: тип аварийной ситуации - "Аварийный сброс в водоток"; основное загрязняющее вещество - азот аммонийный;

Lсб = 30 км; Lx = 80 км; L1 = Lсб = 30 км; B1 = 20 м; H1 = 0,35 м; ucp,1 = 0,70 м/с; umax,1 = 0,90 м/с; Q1 = 5,6 м3/с; y0 = 8 м; nш,1 = 0,025; j1 = 1,3; L2 = 50 км; B2 = 50 м; H2 = 1,5 м; ucp,2 = 0,30 м/с; umax,2 = 0,38 м/с; Q2 = 22,5 м3/с; nш,2 = 0,03; j2 = 1,1; Сф(к) = 0,1 мг/дм3; Свз = 2,5 мг/дм3.

Измеренные характеристики режима сброса сточных вод приведены в таблице В.3. Результаты прогноза представлены в таблице В.4.

 

Таблица В.3 - Характеристики режима аварийного сброса сточных вод

 

Номер проведенного измерения

Дата и время

Концентрация, мг/дм3

Расход сточных вод, м3

1

08.07.2000 г.

10 ч 00 мин

240

0,40

2

08.07.2000 г.

11 ч 00 мин

150

0,35

3

08.07.2000 г.

12 ч 00 мин

200

0,27

4

08.07.2000 г.

13 ч 00 мин

200

0,15


Таблица В.4 - Результаты прогностического расчета

 

Характеристика

При максимальной скорости течения воды

При средней скорости течения воды

Время появления фронтальной части зоны высокозагрязненных вод

10.07.2000 г.

8 ч 16 мин 12 с

10.07.2000 г.

20 ч 26 мин 30 с

Время появления хвостовой части зоны высокозагрязненных вод

10.07.2000 г.

19 ч 20 мин 39 с

10.07.2000 г.

22 ч 9 мин 50 с

Продолжительность времени прохождения зоны высокозагрязненных вод

0 д. 1 ч 4 мин

0 д. 1 ч 43 мин

Продолжительность времени прохождения части зоны с наибольшей концентрацией вещества

0 д. 1 ч 4 мин

0 д. 1 ч 28 мин

Время появления максимальной концентрации вещества

10.07.2000 г.

18 ч 48 мин 25 с

10.07.2000 г.

21 ч 17 мин 37 с

Значение максимальной концентрации вещества

2,5569 мг/дм3

2,68184 мг/дм3

 

 

Приложение Г

(справочное)

 

Библиография

 

1 Голубева М.Т. Влияние сточных вод, содержащих нефть и нефтепродукты, на санитарное состояние водоемов и обоснование гигиенического нормирования их в воде водоемов // Производственные сточные воды. - 1960. - Вып.5. - С.14-33.

2 Изъюрова А.И. Поведение нефти в водоеме // Гигиена и санитария. - 1955. - № 5. - С.15-18.

3 Петру А. Промышленные сточные воды. - М.: Стройиздат, 1965. - 334 с.

4 Грушевский М.С. Неустановившееся движение воды в реках и каналах. - Л.: Гидрометеоиздат, 1982. - 341 с.

5 Быков В.Д. Гидрометрия. - Л.: Гидрометеоиздат, 1965. - 275 с.

6 Gray H.F. Sewerage in Ancient and Medieval Times // Sewage Works J. - 1940. - No 939.

7 Быстров А.В., Клименко О.А. К вопросу определения коэффициента продольной дисперсии в турбулентных потоках // Водные ресурсы. - 1990. - № 5. - С.174-177.

8 Справочник по гидрохимии. - Л.: Гидрометиоиздат, 1989. - 391 с.

9 Филькин Г.В. О подобии полей концентрации загрязняющих веществ, возникающих при сбросе сточных вод в водотоки. - Новочеркасск, 1986. - 11 с. - Деп. в ВНИИТЭИСХ, № 105 ВС-86.

10 Филькин Г.В., Медведева Л.Д. Поле концентраций, создаваемое линейным выпуском. - Новочеркасск, 1980. - 6 с. - Деп. в ВИНИТИ, № 3949-85.

11 Методические основы оценки антропогенного влияния на качество поверхностных вод / Под ред. А.В.Караушева. - Л.: Гидрометеоиздат, 1981. - 175 с.

12 Справочник проектировщика. Канализация населенных мест и промышленных предприятий. - М.: Стойиздат. - 638 с.

13 Де Бор Л. Практическое руководство по сплайнам: Пер. с англ. - М.: Радио и связь, 1985. - 304 с.