Автореферат (1151768), страница 4
Текст из файла (страница 4)
Рис.7. Изменение значения диффузного стока при различных значениях коэффициента скорости распада.
В пятой главе для апробации модели расчета распространения примеси в водотоке была выбрана часть бассейна р. Волги – р. Ока (114-1500км от ее истока). В этой части бассейна р. Волги была разработана расчетная схема и проведено разбиение ствола р. Оки на участки (всего выбрано 14 участков). Для моделирования были зафиксированы показатели количества и качества воды по основному загрязняющему веществу: аммонийный азот. В качестве критического расчетного периода для оценки качества воды была взята летне-осенняя межень.
В рамках создания базы данных также была проведена привязка крупных и средних притоков, а также рассредоточенной боковой приточности к р. Оке. При этом отсутствующие данные по водности р. Оки и ее притоков были вычислены по корреляции стока в рассматриваемых створах со стоком в створе г. Горбатов. Для информационной полноты базы данных потребовалось восстановить ряд поперечников по данным наблюдений гидрометрических постов. Для восполнения недостающих данных и упрощения расчетной схемы был сделан ряд упрощений и допущений:
- модель разработана для ствола р. Оки, все притоки учитывались устьевыми расходами и концентрациями в них;
- боковая приточность (кроме крупных притоков) принималась равномерно распределенной по длине участка;
- загрязняющие вещества поступают непосредственно в ствол р. Оки, кроме поступлений с крупными притоками с известной концентрацией ЗВ;
- загрязняющие вещества не взаимодействуют друг с другом и другими загрязняющими веществами;
- коэффициент скорости распада ЗВ не меняется по длине каждого участка.
В соответствии с допущениями по данным отчета Гидрохимического института (ГХИ) и программы «Клязьма – голубая река» были определены показатели качества вод, поступающих в ствол р. Оки от ее притоков. В частности, для рек Зуша, Упа, Жиздра, Протва, Москва, Проня, Мокша и Теша по данным отчета ГХИ была сделана выборка наблюденных концентраций ЗВ в период ярко выраженной летне-осенней межени, то есть за июль-август. Полученные таким образом ряды наблюденных концентраций были подвергнуты тщательному анализу. При этом оказалось, что концентрации ЗВ в большинстве притоков колеблются в широком диапазоне и имеют ярко выраженные скачкообразные изменения значений средних концентраций на участках ряда. На основании проведенного анализа для каждого ЗВ на каждом притоке индивидуально определялся ряд значений, характеризующих концентрацию данного ЗВ в притоке. При этом среднее арифметическое значение полученного ряда применялось как исходное для калибровки модели на данных 1989г.
Сбросы рассматриваемых ЗВ в ствол р. Оки на участках и в крупных городах были определены по данным ежегодной госстатотчетности 2ТП (водхоз) за 1989г. По условиям наличия контрольных створов с наблюденными концентрациями ЗВ и в соответствии с изменчивостью коэффициента скорости распада ствол р. Оки был разбит на 7 участков: ниже г. Орел – ниже г. Белев; ниже г. Белев – ниже г. Калуга; ниже г. Калуга – ниже г. Алексин; ниже г. Алексин – выше г. Касимов; выше г. Касимов – ниже г. Муром; ниже г. Муром – устье. Сбросы ЗВ в городах Орел, Калуга, Алексин, Серпухов, Ступино, Кашира, Коломна, Рязань, Касимов, Выкса, Муром, Павлово, Богородск, Дзержинск, стоящих непосредственно на стволе р. Оки, принимались по соответствующей таблице отчетности 2ТП (водхоз) и вычитались из сбросов ЗВ в ствол р. Оки на соответствующих участках.
Для определения массы ЗВ на каждом участке предусматривались следующие действия: составлялась пропорция между расходами небольших рек и равномерно распределенной по длине участка боковой приточностью. Масса ЗВ (за исключением сбросов городов) распределялась согласно этой пропорции. Далее были заданы коэффициенты скорости распада ЗВ для всех рассматриваемых участков и проведена тарировка модели качества на основе данных за 1989г. Затем были проведены расчеты с применением усовершенствованного метода расчета распространения примеси. В этих расчетах использовались упомянутые выше данные по сбросам ЗВ и концентрациям ЗВ в притоках. На основе результатов проведенных расчетов был выполнен анализ полученных концентраций ЗВ в стволе р. Оки. Эти концентрации сравнивались с наблюдаемыми концентрациями ЗВ в створах, по данным отчета ГХИ.
В случаях, когда расчетные концентрации ЗВ превышали наблюдаемые, проводилась коррекция коэффициента скорости распада ЗВ в сторону увеличения; в случаях, когда наблюдаемые концентрации превышали моделируемые, корректировка данных для модели заключалась как в уменьшении коэффициента скорости распада, так и в увеличении фоновых концентраций. В случае невозможности достичь приемлемых результатов моделирования при изменении вышеуказанных факторов, проводилась корректировка данных по сбросам ЗВ со сточными водами. При этом принималось предположение о том, что данные отчетности 2ТП (водхоз) занижены, а сама эта отчетность охватывает не все предприятия-источники сбросов ЗВ.
В результате корректировок данных по концентрациям ЗВ в поступающих расходах и коэффициентов скорости распада загрязняющего вещества получены модельные концентрации ЗВ, мало отличающиеся от натурных концентраций.
Рис.8. Плановое положение контрольных створов.
Табл.3. Результаты расчета распространения N-NH3 в русле р. Оки.
№ створа | Расход массы, г/сек | Расход, м3/с | Область | |||
1 | 2 | 3 | Σ | |||
1 | 2 | 3 | 4 | 5 | 6 | 7 |
участок 1 | ||||||
0 | 0 | 0 | 0 | 0 | 3.32 | |
0 | 0 | 0 | 0 | 2.39 | Орловская | |
0 | 0 | 0 | 0 | 2.49 | Орловская | |
1 | 0 | 0 | 0 | 0 | 8.2 | |
участок 2 | ||||||
1 | 0 | 0 | 0 | 0 | 8.2 | |
0 | 18.924 | 0 | 18.924 | 0.1 | Орловская | |
0 | 0 | 0 | 0 | 0 | Орловская | |
2 | 0 | 18.924 | 0 | 18.924 | 8.3 | |
участок 3 | ||||||
2 | 18.924 | 0 | 0 | 18.924 | 8.3 | |
0 | 0 | 0 | 0 | 0 | Орловская | |
0 | 0 | 0.364 | 0.364 | 7.28 | Орловская | |
3 | 16.156 | 0 | 0.343 | 16.499 | 15.580 | |
участок 4 | ||||||
3 | 16.156 | 0 | 0.343 | 16.499 | 15.580 | |
0 | 10.187 | 0 | 10.187 | 22.480 | Орловская | |
0 | 0 | 20.822 | 20.822 | 35.920 | Тульская | |
4 | 13.834 | 8.804 | 20.197 | 42.834 | 73.980 | |
участок 5 | ||||||
4 | 13.834 | 8.804 | 20.197 | 42.834 | 73.980 | |
0 | 48.825 | 0 | 48.825 | 90.580 | Калужская | |
0 | 0 | 1.026 | 1.026 | 8.520 | Тульская | |
5 | 19.041 | 22.330 | 0.692 | 42.063 | 173.100 | |
участок 6 | ||||||
5 | 42.063 | 0 | 0 | 42.063 | 173.100 | |
0 | 8.643 | 0 | 8.643 | 43.620 | Московская | |
0 | 0 | 5.856 | 5.856 | 12.200 | Калужская | |
6 | 21.517 | 6.638 | 3.205 | 31.359 | 228.900 | |
участок 7 | ||||||
6 | 31.359 | 0 | 0 | 31.359 | 228.900 | |
0 | 300.127 | 0 | 300.127 | 91.700 | Московская | |
0 | 0.000 | 0 | 0 | 0 | Московская | |
7 | 22.121 | 258.724 | 0 | 280.846 | 320.600 | |
участок 8 | ||||||
7 | 280.846 | 0 | 0 | 280.846 | 320.600 | |
0 | 0.815 | 0 | 0.815 | 16.300 | Московская | |
0 | 0 | 0 | 0 | 0 | Московская | |
8 | 199.445 | 0.674 | 0 | 200.119 | 336.900 | |
1 | 2 | 3 | 4 | 5 | 6 | 7 |
участок 9 | ||||||
8 | 200.119 | 0 | 0 | 200.119 | 336.900 | |
0 | 0 | 0 | 0 | 0 | Рязанская | |
0 | 0 | 0 | 0 | 0 | Рязанская | |
9 | 102.081 | 0 | 0 | 102.081 | 336.900 | |
участок 10 | ||||||
9 | 102.081 | 0 | 0 | 102.081 | 336.900 | |
0 | 0 | 0 | 0 | 0 | Рязанская | |
0 | 0 | 9.215 | 9.215 | 16.300 | Рязанская | |
10 | 87.240 | 0 | 8.830 | 96.070 | 353.200 | |
участок 11 | ||||||
10 | 96.070 | 0 | 0 | 96.070 | 353.200 | |
0 | 0 | 0 | 0 | 0 | Рязанская | |
0 | 0 | 33.456 | 33.456 | 137.180 | Рязанская | |
11 | 49.040 | 0 | 23.539 | 72.579 | 490.400 | |
участок 12 | ||||||
11 | 72.579 | 0 | 0 | 72.579 | 490.400 | |
0 | 7.544 | 0 | 7.544 | 47.420 | Владимирская | |
0 | 0 | 15.758 | 15.758 | 55.720 | Нижегородская | |
12 | 36.797 | 6.529 | 11.277 | 54.602 | 593.500 | |
участок 13 | ||||||
12 | 54.602 | 0 | 0 | 54.602 | 593.500 | |
0 | 133.770 | 0 | 133.770 | 102.900 | Владимирская | |
0 | 0 | 0.290 | 0.290 | 2.900 | Нижегородская | |
13 | 33.566 | 105.594 | 0.000 | 139.161 | 699.300 | |
участок 14 | ||||||
13 | 139.161 | 0 | 0 | 139.161 | 699.300 | |
0 | 0 | 0. | 0 | 0 | Нижегородская | |
0 | 0 | 70.126 | 70.126 | 21.260 | Нижегородская | |
14 | 107.369 | 0 | 61.251 | 168.620 | 720.600 |