14873 (596911), страница 2
Текст из файла (страница 2)
(4)
і також допускає аналітичний розв'язок.
Нехай постійне надходження забруднювача є відсутнім. Тоді для оперативного прогнозу якості води використовують рівняння поздовжньої дисперсії для несталих умов:
при С(х,0)=0. (5)
У загальному вигляді це рівняння розв'язується чисельними методами.
Двокомпонентні моделі якості води
Двокомпонентні моделі якості води одержали широке поширення при прогнозних розрахунках вмісту органічних речовин і розчиненого кисню у воді. Наявність розчиненого кисню є важливим фактором, що визначає якість води. Кисень відіграє важливу роль як для підтримання життєдіяльності водних організмів, так і для деструкції органічних речовин – забруднювачів. Вміст кисню в поверхневих водах визначається співвідношенням між його надходженням і витратою на окислення органічних речовин.
БСК (біохімічна потреба кисню) є мірою кількості органічної речовини у воді та інтенсивності її деструкції. БСК дорівнює кількості кисню, яка затрачається на окислення органічних речовин, що містяться в одиниці об’єму води за одиницю часу. Найчастіше БСК оцінюють за 5 діб – БСК5 і за 20 діб – БСК20 – повне БСК, ознакою якого є початок процесів нітрифікації в пробі води.
Прогноз величини БСК і вмісту розчиненого кисню в поверхневих водах, як правило, проводиться на основі математичної моделі Стріттера-Фелпса (1925). У загальному випадку система рівнянь Стріттера-Фелпса для турбулентного потоку записується у виді:
(6)
Тут c1 – концентрація органічної речовини;
c2 – концентрація розчиненого кисню;
k1 – константа швидкості біохімічного окислення;
k2 – коефіцієнт аерації, що залежить від температури;
c2* - концентрація насиченого кисню.
Систему можна розв’язати одним з числових методів.
Аналіз розв'язку системи показує, що при відсутності алохтонного надходження органічних речовин в результаті процесу самоочищення вміст кисню у воді виходить на рівень насичення. Якщо у водний об'єкт надходять органічні речовини, то дефіцит кисню збільшується. При великому надходженні органічних речовин розчинений кисень цілком витрачається на біохімічне окислення органічних речовин і вміст його у воді наближається до нуля.
Багатокомпонентні моделі якості води
Багатокомпонентні моделі якості води використовують для прогнозу вмісту у водоймах речовин, пов'язаних між собою процесами взаємної трансформації. Зокрема, в основі кругообігу азоту, фосфору, вуглецю лежать процеси взаємної трансформації речовин, і з цим явищем тісно пов'язаний процес евтрофування водних об'єктів.
Одним з основних циклів взаємної трансформації речовин у водному об'єкті є кругообіг азоту - одного з головних біогенних елементів. В основі трансформації з'єднань азоту у водному середовищі лежить процес нітрифікації. Нітрифікація являє собою процес окислення мінеральних форм азоту і здійснюється в дві стадії. На першій стадії відбувається окислення амонійних форм азоту до нітритних, на другій стадії нітритні форми окислюються до нітратних.
Короткостроковий прогноз (до 4 діб) вмісту у водному об'єкті мінеральних форм азоту з урахуванням процесу нітрифікації в умовах турбулентного потоку провадиться на основі математичної моделі Харлемана. Ця модель має вид:
(7)
де D - коефіцієнт поздовжньої дисперсії, м2/с; v - швидкість потоку, м/с; N1, N2, N3 - концентрації азоту амонійного, нітритного і нітратного відповідно, г/м3; k1, k2 - коефіцієнти трансформації амонійної і нітритної форм азоту, 1/с.
Аналіз розв’язку системи показує, що по завершенні процесу нітрифікації весь азот переходить у нітратну форму. При цьому вміст азоту амонійного у воді постійно зменшується, а азоту нітратного збільшується. Інша картина спостерігається для проміжної, нітритної форми азоту. У процесі нітрифікації вміст азоту нітритного спочатку збільшується, а потім знижується до нуля. Максимум вмісту нітритного азоту може істотно перевищувати його гранично допустиму концентрацію, навіть у випадку відсутності його у водному об'єкті до початку процесу нітрифікації.
Для проведення середньострокових і довгострокових прогнозів використовують моделі Міхаеліса-Ментен-Моно. В процесі нітрифікації мінеральні форми азоту виконують роль субстратів для нітрифікуючих бактерій. Система рівнянь має наступнийвигляд:
(8)
Тут X1, X2 – густина біомаси бактерій 1-го і 2-го виду; max1, max2 – максимальні швидкості росту бактерій; y1,y2 - коефіцієнти урожайності (збільшення біомаси бактерій за рахунок споживання одиниці субстрату); ks1,ks2 - константи напівнасичення (концентрація субстрату, при якій швидкість росту біомаси дорівнює половині максимальної); kd1,kd2 - коефіцієнти смертності.
Ця система диференціальних рівнянь не має аналітичного розв'язку і розв'язується чисельними методами. Дана модель використовується для прогнозу на період до 10 діб. При прогнозуванні на більш тривалий період використовують більш складні моделі, які враховують процеси мінералізації органічних форм азоту, роль фіто- і зоопланктону в колообігу азоту у водній екосистемі.
Модель управлiння водосховищем
Для рiзних господарських потреб необхiдно утворити певний запас прiсної води i управляти цим запасом так, щоб найкращим чином задовольняти потреби у прiснiй водi. Цей запас створюється у водосховищі. Перш за все необхідно визначити оптимальний об’єм водосховища W. При малому вмісті водосховище не зможе оптимально управляти запасом води. Видiлимо фактори, якi впливають на величину запасу води Y.
Приток води по рiчцi (включаючи боковий) - Р.
Випадання опадiв - O.
Випаровування з поверхнi водосховища - V.
Фiльтрацiя води у нижньому створi водосховища - F.
Частина води витрачається на фактори антропогенного походження – сільске господарство комунальне господарство - S.
При повному заповненні басейну водосховища відбувається сток води через греблю – H.
Запас води у водосховищі не може бути меншим деякої мiнiмальної величини хmin, і не може перевищувати об’єму водосховища W. Розподiл води мiж споживачами здiйснюеться пропорцiйно до їх запитів Sсх i Sкг (Sсх + Sкг = S). Необхiдно з’ясувати, як змiнюються з часом фактори, що впливають на величину запасу води х. Цi залежностi виясняють на основi природнiх спостережень.
Залежностi притокiв i опадiв від часу носять частково періодичний а частково випадковий (стохастичний) характер. Це означає, що функцiї їх змiни можна представити у виглядi
(12)
тут індекс p позначає періодичну складову, індекс s – стохастичну складову.
Об’єм фільтрації прямо пропорційний кількості води у сховищі F = ax, де а - деякий коефiцiєнт. Витрати води на греблi залежать вiд запасiв води у водосховищi H = X – W. Величини споживань Scx i Sкв для простоти можна вважати постiйними. Тоді різницеве рівняння для запасу води в басейні матиме вигляд:
(13)
Рiвняння (13) є рiвнянням балансу запасу води у водосховищi в залежностi вiд природнiх ресурсiв i стратегiї розподiлу цього запасу мiж можливостями. Графічна ілюстрація часової залежності деяких величин зображена на рисунку:
Це рiвняння дозволяе дати вiдповiдь на запитання: “Чому дорiвнює запас води у водосховищi у кожний момент часу”. Використовуючи статистичні дані для даної місцевості ми можемо змоделювати функції P(t), O(t), V(t). При цьому періодична складова базуватиметься на даних для середньозабезпечених вологою років, а стохастична складова враховуватиме можливі відхилення в межах спостережуваних за минулі роки. Побудувавши таку математико – імітаційну модель можна змоделювати роботу водосховища в різних умовах, визначити оптимальний об’єм водосховища, визначити можливі обсяги водогосподарського споживання в різні пори року. Все це допоможе побудувати оптимальну стратегію управління водосховищем.
Література
Основна література
-
Ковальчук П.И., Лахно Е.С. Прогнозирование и оптимизация санитарного состояния окружающей среды. - К.:Вища школа, 1988. - 187с.
-
Лаврик В.І. Методи математичного моделювання в екології. Київ: Фітосоціоцентр, 1998. - 132с.
-
Горєв Л.М. Основи моделювання в гідроекології. К.:Либідь, 1996. - 336с.
-
Кузнецов Ю.Н. и др. Математическое программирование. - М.: Высшая школа, 1980.
-
Стрижак Т.Г., Коновалова Н.Р. Диференціальні рівняння. - К.: Либідь,1994.
-
Практикум по рациональному природопользованию / Н.С. Демченко и др.; Под ред. И.Г. Черванева. - К.: УМК ВО, 1991.
Додаткова література
-
Робертс Ф.С. Дискретные математические модели с приложениями к социальным, биологическим и экологическим задачам. - М.: Наука,1986. - 496с.
-
Cемевский Ф.Н., Семенов С.М. Математическое моделирование экологических процессов. - Л.: Гидрометеоиздат, 1982. - 280 с.
-
Петросян Л.А., Захаров В.В. Введение в математическую экологию. - Л.: Изд. Ленингр. ун-та, 1986.
-
Руденко В.Д., Макарчук О.М., Патланжоглу М.О. Практичний курс інформатики. К.:Фенікс, 1997. - 304с.
Навчально - методичні посібники
-
Методические указания по изучению дисциплины “ВТ и программирование: ЧМРЗ на ЭВМ”. Тема: Численное интегрирование функций. 044 -62.
-
Методичні вказівки по вивченню дисципліни “ОТ та програмування; ЧМРЗ на ЕОМ”. Тема: Чисельне розв’язування рівнянь. 044 - 61.
-
Методичні вказівки по вивченню дисципліни “ОТ та програмування; ЧМРЗ на ЕОМ”. Тема: Безумовна оптимізація функції багатьох змінних. 044 - 63.
-
Методичні вказівки по вивченню дисципліни “ОТ та програмування; ЧМРЗ на ЕОМ”. Тема: Одновимірна оптимізація функції. 044 - 67.
-
Методичні вказівки по вивченню дисципліни “ОТ та програмування; ЧМРЗ на ЕОМ”. Тема: Числове розв’язування звичайних дифрівнянь. 044 - 72.
-
Рабочая программа, методические указания и задания для выполнения контрольной работы по курсу "ВТ и программирование: ЧМРЗ на ЭВМ”. 044-56.
-
Робоча програма, методичні вказівки і завдання для виконання контрольної роботи з курсу “Обчислювальна техніка та економіко-математичні методи у землевпорядкуванні” для студ. II курсу заочної форми навчання спеціальності 6.070900. 100 – 08.
-
“Геоінформаційні системи і технології”. Робоча програма, методичні вказівки та завдання для виконання лабораторних і контрольних робіт для студентів денної та заочної форм навчання спеціальності 6.070900. 100 – 17
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
