
книги из ГПНТБ / Проблемы охраны природных и использования сточных вод межведомственный сборник
..pdfбы рандомизировать те систематически действующие факторы, которые трудно поддаются учету и контролю, с тем, чтобы можно было рассматривать их как случай
ные величины и, следовательно, учитывать |
статистиче |
||
ски. Взятие проб воды проводилось в течение |
несколь |
||
ких лет (1970— 1973 гг.) |
из различных |
водоемов в |
|
окрестностях г. Минска, |
в различное время года |
(с мар |
та по ноябрь), со случайными интервалами между опы тами, при различных метеорологических условиях. Одно временно с взятием проб регистрировалась температура воздуха и воды, определялись бихроматная и перманганатная окисляемость, содержание растворенного кисло рода и некоторые другие показатели. Определение зна чений ВПК при различных сроках экспозиции проб воды проводилось в соответствии с рекомендациями сборника «Унифицированные методы анализа вод» [11].
Регрессионный анализ зависимости ВПК от времени экспозиции проб воды. Функции регрессии, как известно, описывают статистическую зависимость изучаемой слу чайной величины от одной или нескольких других вели чии [4]. В простейшем случае статистическая или кор реляционная зависимость проявляется в том, что при изменении одной из величин (в данном случае время экспозиции i) изменяется среднее значение другой (БПК(). При проведении регрессионного анализа необ ходимо решить две основные задачи:
а) установить форму корреляционной связи; б) оценить тесноту (силу) корреляционной связи.
Так как характер зависимости между исследуемыми величинами (t и ВПК) заранее неизвестен и вообще не существуют строгих математических методов, которые «априори» позволили бы выбрать общий вид искомой функции, вид кривой определялся по характеру распо ложения точек на корреляционном поле (рис. 1). Оче видно, что при таком методе выбор общего вида эмпири ческой функции регрессии не является однозначным и возникает необходимость в последующей оценке степени достигнутого приближения эмпирической функции к «истинной».
При нелинейном характере исследуемой зависимости для упрощения вычислительных операций по определе
ний параметров |
функций регрессий ^методом наимень |
ших квадратов |
рекомендуется произвести преобразова- |
31
нпе исходных величин или системы координат таким образом, чтобы полный набор данных давал по возмож ности прямую линию.
Четыре варианта таких преобразовании для одних и тех же данных (см. рис. I) показаны на рис. 2; в табли це даны необходимые -пояснения к рис. 2, приведены со ответствующие уравнения (2) — (5) регрессии БПК на I, для облегчения логического анализа все рассматри ваемые зависимости записаны также в дифференциаль ной форме.
Данные опытных определении БПК при различных сроках экспозиции по 59 пробам воды из различных во дохранилищ в окрестностях г. Минска п в Таджикской ССР были обработаны на ЭЦВМ «Мписк-22» с целью подбора коэффициентов аппроксимирующих уравнении регрессии (21) — (5‘). Программой расчета предусматри валось также проведение дисперсионного анализа и рас чет эмпирических корреляционных отношений для конт роля соответствия опытных данных полученным уравне ниям регрессии.
Результаты расчета п анализ принятых допущений позволили прийти к следующему.
Все исследованные зависимости, кроме описанной дифференциальным ура-виеннем (3) и соответствующим ему уравнением регрессии (З1), дают вполне удовлетво рительную [5] аппроксимацию опытных данных, прн
Рис. 1. Опытные данные по зависимости БПК от времени t для раз личных проб воды: У— Заславль, 26.VI 1970; 2 — Волковнчи 15.VII 1970; 3 — Заславль, 28.IV 1970; 4 — Станки, 25.VI 1970
32
этом коэффициент корреляции для всех 59 проб нахо дится в пределах 0,880—0,995 и доля остаточной, т. е. необъяснениой линейной регрессией, дисперсии состав ляет 5— 10% при максимальном значении 22%. Провер ка гипотез о значимости выборочных коэффициентов корреляции' по критерию Стыодента показала, что все рассчитанные коэффициенты корреляции значимо отли
чаются от |
нуля, |
так как |
/Выч= (8,42 |
— 42,24) > / Прод = |
|
= (1,860 — 2,015) |
при а = 0,05. |
|
|
||
Зависимость |
(З1) следует признать |
неудовлетвори |
|||
тельной и по той причине, |
что физический смысл |
про |
|||
цесса так, |
как его описывает дифференциальное |
урав |
нение (3), явно не соответствует действительному ходу биохимического потребления кислорода в воде. Пред ставляется спорным н физический смысл уравнения (5). Уравнение (21) неудобно для практического использо вания в целях прогнозирования величины полного ВПК
Рис. 2. Графики зависимости БПК от t в различных системах коор динат (по опытным данным рис. 1)
3. Зак . 769 |
зз |
Расчетные значения
|
Преобразова |
|
Рисунок |
ние координат |
|
осипо абсцисс |
по оси |
|
|
|
ординат |
2, а t 1п(БПКп— -ВПК,)
|
Зависимость между ! и ВПК |
||
в диффе |
|
физический |
уравнение регрессии |
ренциаль |
|
смысл |
БПК на t |
ной форме |
|
|
|
dxt |
(2) |
Скорость |
In (БПКп— (2') |
dt = |
|
процесса |
—БПК,)= |
|
пропорцио =1пБПКп—М |
||
=/г(/С0 X/) |
|
||
|
нальна кон |
|
центрации
2, |
б |
t |
1п ВПК/ |
dxt |
|
(3) |
Скорость |
In БПК/= |
(3') |
|
|
|
|
|
Ц Г = Ь ' |
процесса |
:= In a-\-kt |
|
|||
|
|
|
|
пропорцио |
|
|
||||
|
|
|
|
|
|
|
нальна по |
|
|
|
|
|
|
|
|
|
|
требленному |
|
|
|
|
|
|
|
|
|
|
кислороду |
|
|
|
2, |
а |
In t |
ВПК/ |
dxt |
1 |
(4) |
Скорость |
БПК<=С+ |
(4') |
|
|
|
|
|
dt |
kt |
|
процесса |
4~— In t |
|
|
|
|
|
|
|
обратно |
|
||||
|
|
|
|
|
|
|
пропорцио |
1 k |
|
|
|
|
|
|
|
|
|
нальна вре |
|
|
|
|
|
|
|
|
|
|
мени |
|
|
|
2, г |
In t |
In БПК< |
dxt |
xt |
(5) |
Скорость |
In БПК/=1п a + |
(S') |
||
|
|
|
|
dt |
~~ kt |
|
процесса |
- f— In t |
|
|
|
|
|
|
|
пропорцио |
|
||||
|
|
|
|
|
|
|
нальна по |
k |
|
|
|
|
|
|
|
|
|
требленно |
|
|
|
|
|
|
|
|
|
|
му кислоро |
|
|
|
|
|
|
|
|
|
|
ду и обрат |
|
|
|
|
|
|
|
|
|
|
но пропор |
|
|
|
|
|
|
|
|
|
|
циональна |
|
|
|
|
|
|
|
|
|
|
времени |
|
|
|
|
|
Обозначения: |
Ка — количество |
кислорода, необходимое |
для |
|||||
окисления всего органического вещества к |
началу |
процесса; xt— по |
||||||||
требление кислорода за время |
t; а, |
с, k — константы. |
|
по результатам краткосрочных наблюдений. Уравнение
(4) [7] свободно от указанных недостатков, физический смысл его вполне согласуется с логическими представ лениями о ходе процесса, кроме того, как уже указыва лось, оно достаточно хорошо аппроксимирует опытные данные. Все сказанное позволяет сделать вывод о прак-
34
тической возможности описания зависимости величины БПК от времени экспозиции пробы воды в исследован ной области изменения переменной t от 1 до 20 сут в ви де уравнения регрессии (41)
БПК, = с И- — In t.
k
Некоторые из линий регрессии и соответствующие им уравнения регрессии приведены на рис. 3. Там же пока заны данные опытных определений БПК.
6ПК,Мгйг/ л
Рис. 3. Эмпирические линии и уравнения |
регрессии Б П К tна: 1— За- |
|
||||||||||||
славль, |
26.V I |
1970; |
Б П К ; = 1,27+7,30 |
1\lg |
2 — |
Стайки, |
1 5 .V II |
1970; |
||||||
Б П К , = |
1,05+5,30 |
lgt\ |
3 — Сурхоб, 5 .V III |
1971; |
Б П К , = |
0 ,38+ 5 ,3 0 lg |
/; |
|||||||
4 — Плпсса, |
11.I I I |
1970; |
Б П К , = 0 ,5 3 + 4 ,0 8 |
t;lg5 — Волковичн, 2 6 |
.V III |
|||||||||
1970; |
Б П К ,= 0 ,2 7 + 3 ,8 1 |
lgt; 6 — |
Свислочь, |
2 8 .IV 1971; |
Б П К ;= 0 |
,7 1 |
+ |
|||||||
+ |
2,96 |
lg U |
7 — |
Солигорск, |
19.Ш 1971; |
Б П К ,= 0 ,5 9 + 1 ,7 4 tlg |
|
Коэффициент регрессии с, как показал анализ, не имеет статистической связи ни с одним из исследованных показателей качества и состояния воды в период отбора проб и является случайной величиной, распределенной по нормальному закону. Проверка гипотезы о нормаль ном распределении проводилась по критерию согласия Пирсона.
Расчетное значение критерия %1абл = 9,038 потаблице критических точек распределения %а [2] при уровне значи-
з* |
35 |
мости а = 0 ,0 5 и числе степеней |
свободы |
v= /e— 3 = 6, где |
k = 9 —число групп выборки (рис. |
4); х2р |
(0.05; 6) = 1 2 ,6 . |
Так как х2абл < Х2Р> гипотеза о нормальном распределении генеральной совокупности случайной величины с не отвер гается.. Расчетом получены следующие статистические ха рактеристики нормального распределения величины с: сред нее значение с=0,503; среднее квадратическое отклонение о - = 0,715. Таким образом, 95% значений коэффициента
Рнс. 4. Гистограмма и нормальная кривая распределения значений коэффициента с
с лежат в интервале с+;2сгс = 0,503 К1,430. С учетом по лученных данных после перехода к десятичным логарифмам уравнение (4) может быть записано
ВПК, = 0,503- | ~ l g / |
|
(6) |
|
и в таком виде рекомендовано для оценки |
ожидаемой |
ве |
|
личины ВПК при известной зависимости |
константы |
ско |
|
рости К от показателей качества |
и состояния воды в |
пе |
|
риод отбора проб. |
факторов |
на скорость |
|
Анализ влияния различных |
биохимического потребления кислорода. Как уже указы валось, в уравнении регрессии -(6) величины ВПК на время экспозиции t коэффициент К имеет смысл кон станты скорости протекания процесса биохимического потребления кислорода,
36
' Для проведения дисперсионного анализа [2] влия ния различных факторов на величину константы произ водились выборка опытных данных и пх группировка по значениям исследуемого фактора. При этом влияние всех остальных факторов усреднялось (такой прием из вестен под названием «Метод случайного баланса»). Расчеты выполнялись на ЭЦВМ «Мпнск-22».
Результаты анализа позволили установить, что ско рость окисления является характеристикой, статистиче ски связанной с величиной pH воды. Коэффициент кор
реляции |
при |
этом равен |
/'= 0,686. |
Существенность |
ре |
|||||||
грессии явствует из сравнения’ |
вычисленного критерия |
|||||||||||
Фишера |
с его предельным |
значением |
при |
а = 0,05— |
||||||||
— /7выч= 29,'49>Кпр = 4,16. Доля |
остаточной |
дисперсии |
||||||||||
составляет 0,529, т. е. около 47% |
дисперсии константы |
|||||||||||
скорости |
окисления |
обусловлено |
изменением |
pH |
воды. |
|||||||
Скорость |
окисления |
также |
коррелпруется |
с величиной |
||||||||
бнхроматной |
окисляемости воды |
|
(г = 0,472). |
Проверка |
||||||||
по критерию Фишера подтвердила существенность |
|
ли |
||||||||||
нейной регрессии (/7Выч = 7,009Жо,о5 = 4,35). Доля |
оста |
|||||||||||
точной дисперсии в этом случае составляет |
0,778, т. е. |
|||||||||||
около 22% дисперсии константы |
скорости |
обусловлено |
||||||||||
изменением бнхроматной окисляемости. |
|
|
|
|
|
|||||||
Корреляционной зависимости |
между |
исходным |
со |
|||||||||
держанием кислорода в |
воде (KI,bI4 = 2,268<.Fo,o5 = 4,10), |
|||||||||||
а также |
перманганатной |
окисляемостыо |
|
(Квыч= 1,353< |
||||||||
< ^ 0.05 = 4,12) |
и константой скорости не обнаружено. |
|
||||||||||
Не установлено |
также |
статистически |
|
существенной |
разницы между скоростью процесса биохимического по требления кислорода по годам для каждого из обследо
ванных водохранилищ, а также по водохранилищам |
в |
|
пределах одного года. |
с |
|
Расчеты парной корреляции константы скорости |
||
величиной pH и бнхроматной окисляемостыо воды (БО) |
|
|
позволили получить следующие линейные уравнения ре |
||
грессии: |
|
|
— = — 19,466 + 2,893 pH, |
(7) |
|
К |
|
|
'- ^ = 2,554+ 0,117БО. |
(8) |
- |
Учитывая независимость действия указанных факто ров на константу скорости, взаимное влияние обоих фак
торов можно описать в виде суммы линейных функций с удельным весом каждой, пропорциональным вкладу фактора в общую дисперсию. Таким образом, получен ная эмпирическая формула для расчета константы ско рости
— = — 12,4+ 1,971 pH-!- 0.038БО, |
(9) |
К
и после подстановки (9) в уравнение (6) формула пред варительной оценки ожидаемой величины БПК по из вестным значениям pH и бпхроматной окисляемости (БО) в момент отбора проб воды
БПКг = 0,503 + (— 12,4 -f- 1,971 pH -|- 0.038БО) lg/. (10)
Расчеты и сравнение с опытными данными показали, что формула (10) позволяет оценить ожидаемую вели чину БПКго с ошибкой не более 1,5—2,0 лггОг/л непо средственно в день отбора проб воды. Довольно значи тельная величина ошибки вызвана тем, что изменение pH и БО объясняет только 69% дисперсии константы скорости, а 31% обусловлен неучтенными причинами. Тем не менее зависимость (10) может найти применение на практике.
Данные опытного определения п статистического ана лиза различных показателей качества природных вод по зволили сделать следующие выводы.
Процесс биохимического потребления кислорода в во де может быть описан линейным уравнением регрессии значения БПК на логарифм величины срока экспози ции.
Скорость процесса биохимического потребления кис лорода статистически связана с величиной pH воды и бпхроматной окисляемостыо ее. Их суммарное влияние обусловливает около 70,% дисперсии константы скорости процесса.
Для ориентировочной оценки ожидаемого значения БПКго получена эмпирическая зависимость между БПК, pH и БО.
Исследования должны быть продолжены в направ лении выявления дополнительных факторов, влияющих на скорость процесса окисления, уточнения полученных эмпирических зависимостей, дающих возможность уско ренной оценки качества природных вод,
38
Литература
1. Б а з я к и н а Н . А. Значение константы скорости потребл ния кислорода при определении биологической потребности в кисло
роде |
сточной жидкости. |
«Санитарная |
техника», |
1933, |
№ 2. |
|
|
||||||||||
|
2. |
Г м у р м а н |
В. |
Е. |
|
Теория |
вероятностей |
п |
математическая |
||||||||
статистика. М ., 1972. |
|
|
|
|
загрязнением |
рек, |
озер и |
водохра |
|||||||||
|
3. |
Д р а ч е в С. М . Борьба с |
|||||||||||||||
нилищ |
промышленными и бытовыми стоками. М .— Л ., |
1964. |
анализ. |
||||||||||||||
М , |
4. |
Д р е й п е р |
М , |
С м и Гт. Прикладной |
регрессионный |
||||||||||||
1973. |
р у ж |
и и п н |
Н . |
И . |
Л огика |
оценки |
статистических гипоте |
||||||||||
М „ |
5. |
Д |
|||||||||||||||
1973. |
у к о |
в А. И . и |
др. Канализация. М ., |
1969. |
|
|
|
||||||||||
|
6. |
Ж |
В П К |
воды |
по |
||||||||||||
|
7. |
К у д р я ш е в а А. |
|
Г. |
Расчетное |
определение |
|||||||||||
результатам краткосрочных |
наблюдений. |
В сб. «О чистка н |
использо |
||||||||||||||
вание природных и сточных |
вод». М инск, |
1973. |
|
|
1971. |
|
|
||||||||||
|
8. |
Н а л и м о в В. |
В. |
Теория эксперимента. М ., |
|
пол |
|||||||||||
|
9. |
П о п о в а Н . |
М ., |
К |
а р ю х п на |
Т. А . К |
вопросу расчета |
ной биохимической потребности в кислороде для сточных вод различ
ной степени загрязнения. |
«Водоснабжение и санитарная техника |
1963, № 10. |
потребления кислорода в загрязненных |
10. Т е р ь о Э. Скорость |
реках. В сб. «Вопросы загрязнения и самоочищения вод'оемов». М . 1937.
11. |
Унифицированные |
методы анализа вод. |
П од общей ре |
Ю . Ю . |
Лурье. М ., 1973. |
Н . Санитарные условия |
спуска сточных |
12. |
Ч е р к н н с к и й О. |
||
вод в водоемы. М .; 1962. |
|
|
П . В. О С Т А П Е Н Я , Л . Г. М И Х А Н Л И К
(ЦНИИ комплексного использования водных ресурсов)
ПРИРОДНЫЕ ФАКТОРЫ, ОБУСЛОВЛИВАЮЩИЕ ПОВЫШЕННОЕ СОДЕРЖАНИЕ
ЖЕЛЕЗА В ПОДЗЕМНЫХ ВОДАХ
y\yy/Vv^/V\y\AAAAAAAAAAAA/\AAAAAA/\AAA/VV\AAAAA/\/VV^'v\AAAAAAAAAAAAA/VVV\AAAA
Источником различных химических элементов, вхо дящих в состав подземных вод, в том числе и железа, яв ляются породы, с которыми контапктпрует вода при сво ем движении.
Со сменой геологических эпох менялись причины и условия выноса железа из пород в врду. Считается, что в давних геологических эпохах содержалось меньше кис лорода и больше углекислоты, что облегчало миграцию железа. Изменялось, также содержание соединений же леза в породах. В более древних отложениях находится обычно больше соединений железа, следовательно, и воды могли растворять большее количество этих соеди нений. Вероятно, что чем старше становилась Земля, тем меньше железа содержали воды. Подтверждением этого является то, что большинство известных место рождений железа сформировалось в давних геологиче ских эпохах (прокембрнй, кембрий, силур, карбон и юра) [1, 2]..
Вразных видах пород количество железа различно.
Впесчаных отложениях количество железа колеблется от следов до 10%, в глинисто-илистых отложениях его обычно больше, известняки содержат следы железа. Со
держание железа в базальтах довольно значительно — 8,7%, кислые гранитные породы содержат его меньше — 2,4%, в гидросфере еще меньше — 5— 10-6% [9].
Многие специалисты считают, что количество Fe в осадочных песчано-илистых отложрниях зависит от зер нистости и возрастает по мере уменьшения грануляции зерен (пород).
Комбинированное воздействие на железосодержащие породы воды, углекислоты н кислорода приводит к раз
40