Математичне моделювання впливу техногенних навантажень на екологічні системи

Размещено на http://www.allbest.ru/

НАЦІОНАЛЬНА АКАДЕМІЯ НАУК УКРАЇНИ

ІНСТИТУТ ПРОБЛЕМ МОДЕЛЮВАННЯ В ЕНЕРГЕТИЦІ

ІМ. Г.Є. ПУХОВА

УДК 577.34:612.014.482; 614.876+005

Автореферат

дисертації на здобуття наукового ступеня доктора технічних наук

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ ВПЛИВУ ТЕХНОГЕННИХ НАВАНТАЖЕНЬ НА ЕКОЛОГІЧНІ СИСТЕМИ

01.05.02 - математичне моделювання та обчислювальні методи

Сердюцька Людмила Федорівна

Київ - 2004

Дисертацією є рукопис

Робота виконана в Інституті проблем моделювання в енергетиці ім. Г.Є. Пухова Національної Академії наук України

Науковий консультант доктор технічних наук, професор,член-кореспондент НАН України

Євдокимов Віктор Федорович,

Інститут проблем моделювання в енергетиціім. Г.Є. Пухова НАН України, директор

Офіційні опоненти:

доктор технічних наук, професорБідюк Петро Іванович,

Національний технічний університет України „КПІ”,кафедра математичних методів системного аналізу, професор

доктор технічних наук, старший науковий співробітникЛисиченко Георгій Віталійович,

Інститут геохімії навколишнього середовища НАН та МНС України,відділ проблем екологічної безпеки, завідувач

доктор технічних наук, професорТарапон Олександр Григорович,

Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України, відділ моделювання задач екології, завідувач

Провідна установа - Інститут кібернетики ім. В.М. Глушкова НАН України Захист відбудеться “ 22 “ червня 2004 р. о 14 годині на засіданніспеціалізованої вченої ради Інституту проблем моделювання венергетиці ім. Г.Є. Пухова НАН України за адресою: 03164, м. Київ, вул. Генерала Наумова, 15. З дисертацією можна ознайомитись у бібліотеці Інституту проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України(03164, м. Київ, вул. Генерала Наумова, 15)

Автореферат розісланий “ 19 “ травня 2004 р.

Вчений секретар

спеціалізованої вченої ради __________ Семагіна Е.П.

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

багаторівневий візуалізація екологічний моделювання

Актуальність проблеми. Зміни навколишнього середовища під впливом техногенних навантажень набули таких масштабів, що стають реальною загрозою існуванню як самої природи в цілому, так і благополуччю людства як її невід'ємного елемента. Чорнобильська катастрофа зробила свій страшний внесок в погіршення екологічної обстановки в Україні, створивши високі ризики для життя і здоров'я населення за радіаційним фактором.

Величезні масиви різнорідної інформації, що отримані в результаті радіоекологічних, токсикологічних і медико-екологічних досліджень, проведених після Чорнобильської катастрофи, без відповідної обробки не дозволяють виявити цілісну картину взаємозв'язку і взаємозумовленості факторів, під впливом яких формуються негативні наслідки для навколишнього середовища, біоти і людини.

Прагнення до найбільш точного опису складного системного об'єкта призводить до збільшення кількості врахованих факторів і процесів, що утрудняє використання математичних моделей і доказ їхньої адекватності. Аналіз і прогнозування за такими моделями виконуються з помилкою, зумовленою як похибкою застосованих обчислювальних методів, так і неточністю вихідних даних, і ці помилки можуть призводити до катастрофічних неточностей у модельних прогнозах. Наявність навіть у простих екологічних моделях таких неточностей може призвести до прояву динамічного хаосу.

Розробка спеціальних методів і засобів математичного моделювання систем з такими складними багатовимірними характеристиками об'єктів, як екологічні, є актуальною для глибокого дослідження структури системи як єдності компонентів і зв'язків, здійснення контролю над складною екологічною обстановкою при одночасному обліку величезної кількості різнорідних параметрів, оцінки ступеня ризику для здоров'я населення і прийняття обґрунтованих рішень для його мінімізації.

Зв'язок роботи з державними програмами, планами та науковими темами. Робота виконувалась у межах наступних науково-дослідних робіт: 1) ”Разработать принципы и методы эколого-экономического моделирования процессов формирования качества воды в пресных водоемах для оценки антропогенного воздействия на водные объекты, эффективности водоохранных мероприятий и оптимизаций управления”, завдання ДКНТ СРСР № 94 (Інститут гідробіології АН УРСР, 1982); 2) ”Разработка методов моделирования задач экологии и создание имитационных экологических моделей для задач энергетики и гидротехники”, № ГР 01.89.0005284 (Ін-т проблем моделювання в энергетиці АН УРСР, 1988; 3) ”Создание динамических моделей переноса радионуклидов в экосистеме и прогнозирования коллективных доз для оптимизации решений при острых выбросах АЭС”, № ГР 0193U017084 (ІПМЕ НАН України, 1993); 4) ”Моделирование влияния ЧАЭС на экологическое состояние окружающей среды и радиационные нагрузки населения. Разработка системы эколого-эпидемиологического мониторинга в регионах влияния АЭС”, № ГР 0197U019094 (ІПМЕ НАН України, 1997); 5) ”Довгострокове прогнозування стану здоров'я і профілактика захворюваності груп населення екологічно неблагополучних та чистих регіонів України на основі оцінки адаптаційних реакцій до комплексу техногенних факторів (в першу чергу - іонізуючого опромінення) ”, № ДР 0198U004097 (ІПМЕ НАН України, 1998); 6) ”Систематизацiя результатiв вивчення радiоекологiчних та медико-бiологiчних наслiдкiв аварii на ЧАЕС щодо практичного використання та розробки заходiв по зниженню негативних ефектiв цiєї катастрофи для екосистем та населення”, № ДР 0197U007176 (ІПМЕ НАН України, 1998); 7) ”Створення математичних моделей міграції радіонуклідів в екосистемах зони відчуження (в тому числі, трофічними ланцюгами, що ведуть до людини) з метою короткострокового та довгострокового прогнозування радіологічної ситуації та випереджаючої оцінки ризику для здоров'я людини від забрудненого радіонуклідами навколишнього середовища”, № ДР 0199U002714 (ІПМЕ ім. Г.Є. Пухова НАН України, 2001).

Мета і задачі роботи. Мета роботи - розробка спеціальних методів математичного моделювання складних системних об'єктів для виявлення інтегральних властивостей, що характеризують систему як єдине ціле; створення апарату дослідження внутрішньої структури багатовимірних об'єктів за допомогою математичних методів системного аналізу за даними моніторингових спостережень (бази натурних даних) на прикладі екологічних систем; створення методики контролю над складною екологічною обстановкою, що використовує математичні моделі і великі масиви різнорідних параметрів.

Для досягнення зазначеної мети були поставлені і вирішені такі основні задачі:

Розроблено багаторівневі моделі візуалізації структури великих масивів різнотипної натурної інформації і технологію перетворення вхідних потоків даних на кожнім рівні візуалізації.

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

Розроблено методику й автоматизовані засоби, що забезпечують збереження та багаторівневу обробку інформації при візуалізації складних взаємозв'язків параметрів, отриманих з багатовимірних медико-екологічних баз даних.

Об'єктом дослідження є складні багатокомпонентні системи (природні об'єкти, забруднена екологічна система ”навколишнє середовище - людина”).

Предметом дослідження є способи обробки великих масивів зашумленої натурної інформації, якісні моделі систем великої розмірності.

Наукова новизна одержаних результатів. На основі об'єднання можливостей методів багатовимірного статистичного аналізу і різноманіття сучасних комп'ютерних способів візуального відображення даних уперше синтезована нова багаторівнева стратегія моделювання впливу техногенних навантажень (спільно хімічних і радіаційних) на екологічні системи.

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

Розроблено методику оцінки забруднення території у формі багатовимірних інтегральних екологічних індексів забруднення з урахуванням динаміки основних токсикантів, у тому числі радіонуклідів, у різних ланках екосистеми.

Розроблено методи багаторівневої стратегії упорядкування великих масивів даних (натурних і модельних) для візуального аналізу тенденцій і прогнозування розвитку складної екологічної ситуації. За допомогою запропонованих індексів забруднення розроблені оцінки екологічного і медичного ризику.

Уперше розроблено моделі для картографічного відображення просторово-розподілених кореляційних залежностей між інтегральними індексами забруднення і характеристиками стану здоров'я населення, що проживає на даних територіях.

Розроблено багаторівневий підхід до моделювання, який полягає в наступному:

візуалізація 1-го рівня (структура даних, способи одержання, методи візуалізації). Результатом візуалізації 1-го рівня є формування багатовимірного вхідного потоку даних для 2-го рівня. До 1-го рівня належать бази даних моніторингових спостережень, а також математичні моделі екологічних процесів;

візуалізація 2-го рівня (структура даних, способи одержання, методи візуалізації). Завданням візуалізації 2-го рівня є системний аналіз даних вхідного потоку інформації від 1-го рівня з установленим порядком використання методів та законами нормування;

візуалізація 3-го рівня - формування візуального образу на основі математико-картографічного моделювання з просторовою інтерпретацією результатів послідовного моделювання у виді:

інформаційних карт узагальнених факторів техногенного забруднення територій (з використанням ГІС технологій);

ізоліній і контурних поверхонь узагальнених факторів;

екологічних ”портретів” місцевості;

просторово-розподілених кореляцій;

карт ризиків проживання на забруднених територіях і об'єктах;

узагальнених фазових ”портретів” математичних моделей.

За допомогою розробленого апарату якісного дослідження багаторозмірних систем створені та проаналізовані:

моделі кругообігу біогенних елементів у прісноводних і морських екосистемах;

концептуальні динамічні моделі переносу для 18 радіонуклідів;

динамічні моделі опромінення від усього спектра радіонуклідів для дослідження аварійних ситуацій природного і техногенного характеру та розробки заходів з ліквідації їх наслідків.

7. Розроблено комп'ютерні програмні засоби для реалізації методів багаторівневої стратегії упорядкування великих масивів даних (у тому числі - засоби для створення і корекції спеціальних екологічних карт просторових кореляцій і ризиків проживання на забруднених територіях).

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

Практичне значення отриманих результатів. Запропоновані методологічні підходи, методи, результати моніторингових і модельних досліджень дають змогу виявити цілісну картину взаємозв'язку і взаємозумовленості техногенних факторів, під впливом яких формуються негативні наслідки для навколишнього середовища і людини, дати обґрунтовані оцінки впливу техногенних факторів на стан здоров'я населення, оцінити ступінь екологічного ризику.

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

Реалізація результатів роботи здійснювалась під час виконання окремих складових частин роботи. Отримані результати увійшли у понад 25 наукових звітів. Основні теоретичні результати роботи впроваджено у наступних практичних розробках:

За розробленою методикою побудовано карти ризику проживання на території Житомирської області, розраховані за інтегральним індексом забруднення з урахуванням хімічних і радіаційних показників (1991-1993 рр.).

Проведено кореляційний аналіз і побудовано картографічні матеріали результатів для узагальнених факторів забруднення і репродуктивних функцій жіночого населення сільської місцевості Житомирської області (Лугинського, Народицького, Овруцького і Олевського районів). Використання багаторівневого аналізу для оцінки клініко-лабораторних даних дає змогу прогнозувати кінець вагітності та пологів у жінок з резус-конфліктом.

Розроблено структури баз даних радіоекологічного забруднення для регіонів впливу ЧАЕС - Житомирської області (1985-1993 рр.), Зони відчуження (1986-1989 рр.).

Розроблено структуру бази даних захворювань, пов'язаних з погіршенням екологічного стану (до и після аварії) в регіонах впливу ЧАЕС ”Показники здоров'я и надання медичної допомоги потерпілим від наслідків аварії на Чорнобильській АЕС”; база містить понад 100 форм, передбачає подальше поповнення цифрових показників, графічну і картографічну візуалізацію.

Розроблено базу даних щодо забруднення повітряного басейну низки промислових міст України (1986-1995 рр.); з використанням запропонованих в роботі методик отримано картографічний матеріал розподілу узагальнених індексів забруднення, кореляційні залежності очікуваної тривалості життя і атмосферного забруднення.

Розроблено базу даних щодо захворюваності та дитячої смертності Хмельницької області, отримано картографічний матеріал рівня смертності дітей до 1 року (1989-1999 рр.), впливу забруднення атмосфери на захворюваність (2000 р.).

Створено і проаналізовано за допомогою розробленого нового апарату якісного дослідження багаторозмірних систем наступні математичні моделі аварійних ситуацій природного і техногенного характеру з метою ліквідації їх наслідків:

процесів кругообігу біогенних елементів в прісноводних водоймах для прогнозування ефективності запланованих водоохоронних заходів за показниками концентрації неорганічних і органічних забруднень. Отримано оцінки максимально можливих навантажень на Куйбишевське водосховище зі збереженням існуючої структури екосистеми;

моделі дослідження анаеробних зон і виявлення найбільш чутливих параметрів Балтійського моря з визначенням параметрів біфуркацій;

моделі для перевірки антропогенного характеру виникнення та швидкого поширення сірководневої зони Чорного моря;

концептуальні динамічні моделі переносу для 18 радіонуклідів (⁸⁹Sr, ⁹⁰Sr, ⁹⁵Zr, ⁹⁵Nb, ¹⁰³Ru, ¹⁰⁶Ru, ¹³¹I, ¹³³I, ¹³²Te, ¹³⁴Cs, ¹³⁶Cs, ¹³⁷Cs, ¹⁴⁰Ba, ¹⁴¹Ce, ¹⁴⁴Ce, ²³⁹Pu, ²⁴⁰Pu, ²⁴¹Am) для екологічних і дозових прогнозів внутрішнього опромінення від заковтування для всіх вікових категорій населення за умов багатоцільового сільськогосподарського виробництва при аварійних викидах на АЕС. Моделі ідентифіковані за умовами аварії на ЧАЕС;

динамічні моделі опромінення (зовнішнє, заковтування, інгаляція) від усього спектра водночас діючих радіонуклідів.

Розроблено програмний пакет, який реалізує розроблені математичні методи, дозволяє моделювати в діалоговому режимі та відображувати в графічному виді динаміку переносу радіонуклідів (⁸⁹Sr, ⁹⁰Sr, ⁹⁵Zr, ⁹⁵Nb, ¹⁰³Ru, ¹⁰⁶Ru, ¹³¹I, ¹³³I, ¹³²Te, ¹³⁴Cs, ¹³⁶Cs, ¹³⁷Cs, ¹⁴⁰Ba, ¹⁴¹Ce, ¹⁴⁴Ce, ²³⁹Pu, ²⁴⁰Pu, ²⁴¹Am), дози (зовнішні, заковтування, інгаляція), які формуються за певний період від кожного радіонукліда окремо і по їх сумішах залежно від заданого аварійного радіонуклідного випадіння. Даний пакет має практичне використання в задачах прийняття управлінських рішень в аварійних ситуаціях на енергетичних об'єктах.

Отримані акти впровадження результатів роботи наведено в Додатку до дисертації.

Особистий внесок. У роботах, написаних у співавторстві, автору належать: [1, 3, 15, 16, 18, 19] - системний аналіз і реалізація методики багаторівневих досліджень, методи візуалізації отриманих результатів; [2] - якісний аналіз і реалізація моделі; [4] - постановка задач і розробка моделей колективних доз населення; [9, 10, 12, 13] - методика обробки і система візуалізації багатовимірних медико-екологічних даних; [17] - методика побудови радіоекологічних карт розподілу інтегрального індексу забруднення; [20] - методика розробки карт екологічного ризику проживання на забруднених територіях; [22] - розробка концептуальної схеми автоматизованої системи керування, пов'язаної з радіаційним впливом на екосистеми; [23] - розробка та верифікація математичних моделей екологічних процесів.

Апробація роботи. Основні положення дисертаційної роботи доповідалися та обговорювались на наступних наукових конференціях і семінарах: Міжнародна конференція ”InterCarto9” - “ГІС для сталого розвитку територій” (Новоросійськ - Севастополь, 2003); III Міжнародний з'їзд з радіаційних досліджень (радіобіологія і радіоекологія) (Київ, 2003); 2-га Міжнародна конференція ”Віддалені медичні наслідки Чорнобильської катастрофи” (Київ, 1998); Міжнародна конференція ”Антропогенно змінене середовище України: ризики для здоров'я та екологічних систем” (Київ, 2003); Міжнародна конференція ”TESEC” - ”П'ятнадцять років чорнобильської катастрофи. Досвід подолання” (Київ, 2001), Міжнародні конференції ”ACS 98-99” - ”Advanced Computer Systems” (Szczecin-Poland, 1998, 1999); Наукові конференції Інституту проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України (Київ, 1988, 1993-2004); Всесоюзні наукові конференції в Інституті біології південних морів АН УРСР ”Совершенствование управления развитием рекреационных систем” - (Севастополь, 1985, 1988, 1989), XII Школа-семінар “Математическое моделирование в проблемах рационального природопользования” (Абрау-Дюрсо, 1988); Всесоюзні наукові семінари з класифікацій (Міас, 1985, Пущино, 1986); Наукова конференція з Міжнародного проекту “Балтика” (Ленінград, 1987).

Публікації. За темою дисертації опубліковано 48 наукових праць в міжнародних і вітчизняних виданнях, у тому числі 1 монографія, 24 статті в спеціальних наукових виданнях, 10 публікацій у збірниках матеріалів і праць конференцій.

Структура й обсяг роботи. Дисертація складається з 368 стор. тексту, 124 рис., 46 табл.; містить вступ, сім спеціальних розділів, висновки, список використаних джерел, чотири додатки.

ОСНОВНИЙ ЗМІСТ РОБОТИ

У вступі обґрунтовано актуальність напрямку досліджень, відбито зв'язок дисертаційної роботи з науковими програмами, планами, темами, сформульовано мету і задачі досліджень, показано наукову новизну отриманих результатів і практичну цінність роботи, наведено відомості про апробацію результатів і їхнє впровадження.

У розділі 1 розглянуто основні проблеми математичного моделювання техногенного впливу на навколишнє середовище (включаючи людину) і наведено широкий огляд літературних джерел, пов'язаних з цією темою.

При розробці методологічної частини роботи і теоретичних узагальнень автором використано результати опублікованих досліджень вітчизняних і закордонних учених з проблем математичного моделювання й обчислювальних методів, які застосовувались у розв'язанні медико-екологічних задач, а саме: з України - Бар'яхтара В.Г., Бєляєва В.І., Георгієвського В.Б., Гродзинського Д.М., Євдокимова В.Ф., Івахненка О.Г., Лаврика В.І., Лисиченка Г.В., Тарапона О.Г.; з Росії - Айвазяна С.А., Ізраеля Ю.А., Логофета Д.О., Марчука Г.І., Свирежева Ю.М., Терьохіної А.Ю., Тікунова В.С., Черв'якова В.А.; із США - Джефферсона Дж., Ерроусміт Д., Краскел Дж.; зі Швеції - С'еберга С., Янсона Дж. та ін.

Опис динаміки природних об'єктів спирається на уявлення про їхню системну організацію. Системний підхід до вирішення екологічних проблем передбачає комплексне вивчення процесів, що відбуваються у ландшафтно-географічному середовищі.

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

Проаналізовано сучасні підходи до оцінки екологічного ризику за умов техногенного впливу, де екологічний ризик - це кількісна міра небезпеки виникнення негативних змін у природному середовищі та погіршення стану здоров'я людей. Пропонується інтегральний підхід до оцінки екологічного ризику проживання населення на забруднених територіях. При цьому ризик розраховується за формулою:

де R_i - екологічний ризик проживання на i-му об'єкті (0?R_i?1);

- евклідова відстань між ознаками І у багатовимірному просторі для натурних об'єктів;

I_i - інтегральний індекс забруднення i-го об'єкта (методику визначення I_i наведено у розділі 4);

I_h - еталон, що має найгірший показник;

Запропоновано наступні 9 градацій інтегрального екологічного ризику R_i:

0<0,1	- немає ризику	0,3<0,4	- нижче середнього	0,6<0,7	- істотний
0,1<0,2	- дуже слабкий	0,4<0,5	- середній	0,7<0,8	- високий
0,2<0,3	- слабкий	0,5<0,6	- вище за середній	0,8	- дуже високий

Сформульовано основні задачі практичного використання методів математичного моделювання: розробка алгоритмів інтегральної оцінки техногенного впливу на екосистеми; розробка методів системного аналізу зашумленої натурної інформації; удосконалювання методики оцінки впливу на навколишнє середовище; розробка методики дослідження й оцінки комбінованого впливу (хімічних, радіаційних) техногенних навантажень на здоров'я людини.

У розділі 2 наведено результати дослідження основних етапів моделювання екологічних наслідків техногенних навантажень.

У даному розділі розглядаються математичні методи аналізу великих масивівданих з позицій системного підходу для виявлення інтегральних характеристик системи і їхніх візуалізацій.

При дослідженні структури даних визначено, що варто обережно підходити до вибору методів багатовимірного аналізу, оскільки, як правило, в екологічних задачах відсутня апріорна інформація про інформативність того чи іншого параметра. У зв'язку з цим необхідно вибрати такий підхід, коли будується нова сукупність ознак, кожна з яких за якимось правилом пов'язана з вихідними. До сформованих нових ознак варто поставити наступні додаткові вимоги: взаємна некорельорованість; найменші викривлення образів у просторі у порівнянні з вихідним простором R^M; найбільша надійність правильної розбивки вихідної вибірки X на природні групи і т.д. Ці вимоги є визначальними для вибору методу зниження розмірності простору опису.

Існуючі математичні підходи до зниження розмірності простору можна розбити на дві групи:

Лінійні методи зниження, у яких у тій чи іншій формі використовується ортогональне проектування пучка векторів X_n з R^M у R^M*.

Нелінійні методи зниження розмірності, поєднувані тим, що оцінка якості перетворення R^M>R^M* здійснюється по величинах перекручувань взаємних відстаней між елементами вибірки Y, які представляють елементи вибірки X у просторі зниженої розмірності.

До лінійних методів багатовимірного аналізу даних належать метод головних компонент (МГК) і факторний аналіз (ФА).

Головним об'єктом перетворення у ФА є кореляційна матриця з коефіцієнтів кореляцій Пірсона. У роботі використовується перетворення кереляційної матриці R-техніка, коли коефіцієнти кореляції обчислюються між змінними , і вихідна матриця X стискається по стовпцях, тобто число ознак зменшується з m до p.

Одним з найбільш поширених прийомів пошуку факторів є МГК. Його головна відмінність від ФА полягає в тім, що головні компоненти (фактори) пов'язані зі змінними лінійними функціями перетворення.

Основну модель методу головних компонент можна представити в матричному виді:

Очевидно, що невідомими є матриці A і F.

До відомих труднощів використання лінійних методів належать:

а) необхідність перетворення вихідних даних, наведених у різних одиницях виміру, у порівнювальні одиниці. Традиційно використовується нормування за стандартними відхиленнями.

б) оцінка того, скільки головних компонент (узагальнених факторів) необхідно побудувати для представлення аналізованих вихідних змінних.

Рішення про те, коли варто зупинити процедуру виділення головних компонент, залежить головним чином від того, що вважати малою часткою дисперсії. Це рішення приймається досить довільно, хоча є два критерії: критерій Кайзера (Kaiser) і критерій ”кам'янистого осипу” Кеттелла (Cattell), що у більшості випадків дозволяють раціонально вибрати кількість головних компонент. При цьому враховується те, що аналіз складових з малими величинами власних значень навряд чи доцільний тому, що вони можуть виявитися статистично недостовірними через помилки різного походження. Через те, що ілюстративною метою ФА часто є одержання факторного відображення в графічному виді, звичайно обмежують p = 2 (зі значною втратою інформації) і дають зображення факторного простору у двовимірному зрізі, оскільки виконати це для трьох і більше виділених головних компонент (факторів) проблематично.

У даній роботі лінійні методи ФА використовуються для одержання головних компонент (факторів), кількість яких залежить тільки від накопиченої дисперсії (пропонується не менше 80%), і відсутні обмеження на розмірність простору (p ? 2).

Щоб зрозуміти, яка екологічна реальність прихована в знайдених факторах, проведено аналіз кореляцій факторних навантажень з вихідними змінними. Для підвищення інтерпретованості факторів використовується метод варимаксного обертання VARIMAX, що дозволяє домогтися більшої ”виразності” матриці факторних навантажень. Суть методу полягає в зміні координатних осей, утворених факторами, з метою отримати більш контрастні навантаження так званої простої факторної структури.

Слід зазначити, що лінійні методи поряд з такою перевагою, як представлення в просторі зниженої розмірності структури вихідних даних з мінімальною помилкою, мають недоліки: неможливість вибору для візуалізації тільки 2-3-х головних факторів, залежність значень елементів обчислюваної коваріаційної матриці від розміру аналізованої вибірки, тотожність коваріаційних матриць для деяких різних за характером розподілу вихідних даних, залежність результатів відображення від способу нормування вихідних ознак.

В основі нелінійних методів відображення векторів з багатовимірного простору опису R^L у простір R^L* (L^*?L), які використовуються в роботі, лежить ідея багатовимірного шкалювання (БШ) даних (Дж. Краскел). Відповідно до методу аналізуються матриці близькості, сформовані на вихідній множині елементів . Під близькістю звичайно розуміють подібність розходження між елементами, що оцінюється різними показниками їхньої близькості (кореляцією, ступенем перетинання, відстанню тощо). БШ - це не стільки визначена математична процедура, скільки спосіб найбільш ефективного розміщення об'єктів, що приблизно зберігає відстань між ними у новому просторі ознак, розмірність якого істотно менша за вихідний (Терьохіна А.Ю., Дж. Краскел).

Методи нелінійного шкалювання даних (МНШ) з R^L у R^L*, на відміну від лінійних методів, за своєю суттю орієнтовані на збереження структури даних при зниженні розмірності простору описів і не пов'язані з жодними обмеженнями статистичного розподілу багатовимірних даних. Крім того відображення в порівнянні з лінійними дозволяють зберегти структуру багатовимірних даних у R^L*. Однак визначено, що основним недоліком алгоритмів нелінійного відображення є односпрямованість перетворення даних, а саме - якщо зроблено перехід з вихідного багатовимірного простору R^L у R^L*, то зворотне відображення неможливе.

У даному розділі проаналізовані переваги й основні труднощі традиційного використання математичних методів зниження ознакового простору при обробці багатовимірних екологічних масивів даних.

У результаті проведеного аналізу встановлено, що розв'язання задач, пов'язаних з виявленням внутрішньої структури багатовимірних масивів натурних (і модельних) даних не забезпечується вибором одного з розглянутих методів аналізу, і потрібна розробка нових підходів і способів, що дозволять ефективно згорнути багатовимірний простір вихідних ознак і візуалізувати структуру вихідних даних у просторі 2-3 головних факторів з найменшою втратою інформації.

У розділі 3 формулюються основні методичні підходи для розробки структур баз даних як моделей візуалізації 1 рівня (збір, збереження, первинна обробка з візуальним аналізом даних) для розв'язання задач медико-екологічного моніторингу.

Бази натурних даних (БД) для моделі візуалізації першого рівня є основою для проведення системного аналізу даних медико-екологічного моніторингу з метою виявлення комбінованих ефектів техногенного впливу.

У розділі 3 визначені вимоги до екологічних баз даних:

Структура БД повинна забезпечувати можливість використання даних для вирішення різнорідних задач.

Структура БД повинна виключати надмірності даних.

Необхідне додавання нових даних і коректування старих повинні здійснюватися без зміни існуючих програм, тобто БД повинна бути незалежною.

БД повинна бути захищена від збоїв і випадків некомпетентного чи зловмисного звертання до них осіб, що можуть їх змінити. Дані про навколишнє середовище накопичуються протягом тривалого часу, багато з них є невідтворюваними. Екосистеми описуються великою кількістю показників різної природи, між якими існують багаторівневі зв'язки і відносини типу „багато до багатьох”.

Виходячи з усього вищезазначеного найбільш вдалою моделлю БД, що підходить для екологічних даних, які описують наземні і водні екосистеми, варто визнати реляційну модель БД. Реляційна модель представляє дані у виді двовимірних таблиць із змінюваними взаємозв'язками між різними стовпцями.

Екологічні показники включають: радіаційне забруднення (щільність випадіння ¹³⁷Сs на ґрунт, вміст у продуктах харчування, дозові навантаження), мінеральні добрива і пестициди (витрата на гектар ріллі і на душу населення), важкі метали (концентрація в ґрунті). Медичні дані містять 13 показників.

Оброблялися дані по: Лугинському району - для 50 сіл; Народицькому району - для 47 сіл; Олевському району - для 45 сіл; Овруцькому району - для 47 сіл.

У розділі 4 пропонується математичний метод багаторівневого аналізу багатовимірної інформації, досліджуються проблеми, пов'язані з аналізом і обробкою великих масивів натурних даних про наслідки аварії на ЧАЕС.

Запропоновано нелінійний метод зниження розмірності вихідного простору ознак (багаторівневий алгоритм візуалізації багатовимірних даних), що сполучить особливості лінійних (МГК, ФА) і нелінійних (МНШ) методів відображення.

Суть запропонованого методу полягає в наступному:

Блок 1. З бази натурних даних формується запит (вибірка) - матриця вихідних даних розмірності nЧm. Для забезпечення зіставлення отриманих результатів, усі показники нормуються за формулою:

Блок 2. Пронормовані дані вводяться в блок “лінійні методи”, де при використанні МГК вибирається стільки головних компонент, щоб кумулятивна дисперсія була не менше 75%, і здійснюється перехід до виконання блоку ФА.

Блок 3. Проводиться змістова інтерпретація нових узагальнених факторів (головних компонент) і візуалізація даних у просторі головних факторів, після чого здійснюється перехід до аналізу інформації в новому факторному просторі за допомогою блоку “нелінійні методи (МНШ).

Блок 4. МНШ застосовується не до вихідного опису багатовимірних даних, а до вторинного - до нового факторного простору ознак. Складаємо матрицю евклідових відстаней у новому факторному просторі:

Блок 5. Виконується візуалізація структури даних у новому факторному просторі у виді діаграм, гістограм, чи графіків у вигляді картографічних зображень.

На наступному малюнку наведено приклад багаторівневої візуалізації радіоекологічної обстановки в Зоні відчуження за 1995 рік. Обробці підлягали дані, що містять агрохімічні і радіаційні характеристики ґрунтів у Зоні відчуження. З вірогідністю 89,9 % вихідна таблиця даних була ”стиснута” до трьох факторів. Перший фактор можна розглядати як узагальнений агрохімічний показник. Другий фактор природно інтерпретувати як показник кислотності. У третій фактор увійшли всі показники радіоактивного забруднення, і він розглядається як узагальнений радіаційний показник.

Наводиться проекція тих же даних на осі першого і третього факторів. У цьому випадку утворилося три кластери. Найбільший з них відрізняється високою щільністю точок, що говорить про їхню однорідність за агрохімічними і радіаційними характеристиками. В окремий кластер потрапила точка 18 (ґрунт дерено-підзолистий супіщаний). Третій кластер (точки 8, 9, 10) об'єднав точки з торф'яним типом ґрунту, причому точки 8 і 9 у цьому випадку практично не розрізняються.

Представлено результати візуалізації за допомогою запропонованого багаторівневого алгоритму інтегральної характеристики системи в цілому у вигляді кругової діаграми екологічних станів ряду населених пунктів Зони відчуження. Кожен об'єкт, описаний 15 характеристиками, представлений у вигляді кола на діаграмі. Цю діаграму можна назвати своєрідним радіаційним “портретом” Зони за 1995 рік, де кожен об'єкт характеризується “близькістю” (відстанню) до “найгіршого” еталона, розташованого в центрі “портрета” (чим ближче, тим гірше). Одержання таких екологічних портретів дає змогу ранжувати самі об'єкти за ступенем сприятливості умов навколишнього середовища для проживання на даних територіях. Якщо ці портрети розглядати в динаміці, то буде отримана повна картина адаптації (чи розвитку) екосистеми в цілому за умов радіаційного впливу.

Проведені дослідження показали, що на базі методів багатовимірного аналізуданих (МГК, ФА, МНШ) можна представити складну інформацію про природні системи у виді одного чи декількох узагальнених індексів. Ці індекси дають досить адекватнеуявлення про внутрішню структуру досліджуваної множини радіоекологічних об'єктів.

У розділі 5 розроблено системний концептуальний підхід до дослідження складних динамічних систем, до яких належать екологічні системи.

Системний підхід до математичного моделювання в екології на сьогодні в основному асоціюється з побудовою імітаційної моделі окремої екосистеми чи імітаційно-балансової моделі регіону (Свирежев Ю.М., Логофет Д.О.). Проте для застосування у практичних задачах аналізу екологічної інформації необхідно визначити адекватність цих моделей, що досить складно реалізувати через необхідність знати екологічні механізми дії кожного фактора даної системи для створення адекватної моделі. Додаткові дослідження кожного об'єкта потребують розв'язку істотно складнішої задачі, ніж просто інтерпретація отриманої інформації.

Для побудови та якісного дослідження складних математичних моделей екологічних процесів розроблено методику верифікації (зв'язності), яка включає системи диференціальних рівнянь зі змінними параметрами, і методику верифікації за узагальненими фазовими портретами, яка забезпечує якісне дослідження моделей за розробленим багаторівневим методом візуалізації.

Методика зв'язності полягає в послідовному розширенні та ускладненні структури взаємозв'язків фазових змінних початкової (концептуальної) моделі з невеликою кількістю рівнянь, що легко піддається теоретичному аналізу (зокрема, повинні виконуватись умови існування позитивної та стійкої стаціонарної точки) і в якій закладена суть модельованого процесу. Наступні етапи моделювання полягають у введенні фазових змінних (мінімальна кількість нових зв'язків) зі збереженням попередніх властивостей, закладених у вихідну модель.

Продемонструємо методику на конкретному прикладі математичного моделювання процесів кругообігу азоту в прісноводній екосистемі.

Об'єктом моделювання було Куйбишевське водосховище. Ставилася задача виявити головні причини цвітіння водойми і зробити прогноз продукційної здатності водосховища, з огляду на існуючі техногенні навантаження.

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

Як динамічні змінні використовуються такі величини: P - кількість фітопланктону у водоймищі; B - кількість бактеріопланктону; Z - кількість зоопланктону; N₀ - кількість розчиненого і детритного органічного азоту; N₁ - кількість мінерального азоту у водоймищі. Останній у свою чергу підрозділяється на такі компоненти: NH₄- кількість амонійного азоту (N₁); NO₃ - кількість нітратного азоту у водоймищі (N₂). Усі перераховані вище величини розраховані в тоннах азоту.

Модельовані процеси.

Фізичні процеси: Скидання у водойму неорганічного С₁ і органічного С₂ азоту за добу (де С₁, С₂ - константи). Стік азоту N через греблю для мінерального N₁ і органічного N₂ азоту за добу. N = f (N₁, N₂), де f - лінійна функція.

Біологічні процеси: Фіксація вільного азоту. Запас зв'язаного азоту S у водоймі може поповнитися шляхом засвоєння молекулярного азоту повітря фіксаторами азоту В (різними видами Azotobacter, Clostridium pasteurianum) і деякими синьо-зеленими водоростями P. S = f₁ (B, P), де f₁ - лінійна функція. Для модельних умов записується:

- азотфіксація фітопланктоном за добу дорівнює k₀P;

- азотфіксація бактеріопланктоном за добу дорівнює k₁₃B.

Засвоєння мінерального азоту. Ріст біомаси фітопланктону залежить від кількості поживних речовин (мінеральних сполук). Запишемо, користуючись законом взаємодії мас, що споживання фітопланктоном P неорганічного азоту N₁ за добу дорівнює k₁PN₁.

При моделюванні цей процес розбивається на дві стадії:

Денітрифікація й інші процеси, пов'язані з втратою азоту у водоймі. Перехід азоту в газоподібну форму з нітратів, тобто , називається денітрифікацією.

Цей процес поділяється на дві стадії:

втрати неорганічного азоту N₁ у водоймі за рахунок денітрифікації за добу дорівнюють k₁₆BN₁ (за законом взаємодіючих мас);

перехід азоту в газоподібну форму (вихід із системи) за добу дорівнює k₁₄B (лінійна залежність).

Виділення зоопланктону Z в моделі за добу дорівнюють k₇Z (лінійна залежність).

Смертність. При моделюванні враховується також і смертність живих компонентів (P, Z, B), тобто перехід у мертву органіку (N₀):

На перших етапах розглядаються моделі тільки з постійними коефіцієнтами. Блок-схеми кругообігу азоту в Куйбишевському водосховищі з послідовним розширенням й ускладненням структури взаємозв'язків фазових змінних показані на рис. 6 (а-г).

Вираз для стаціонарної точки моделі (3) через співвідношення параметрів моделі:

Реальна стаціонарна точка відома з натурних даних:

Лінеаризуємо систему (3), зробивши підстановку:

.

Характеристичний багаточлен системи (4) має такий вигляд:

Необхідні умови стійкості:

або, підставивши вираз для стаціонарної точки, отримаємо:

Достатня умова стійкості для багаточлена третього порядку

виглядає як нерівність тобто

Умови (5) для числових значень коефіцієнтів початкової моделі (3) виконуються.

Розширимо модель (3), додавши нову фазову координату Z - кількість зоопланктону у водоймі. Вибір нової змінної заснований на тім, що введення саме її в систему пов'язаний з додаванням мінімальної кількості нових зв'язків.

Числові значення тих коефіцієнтів розширеної моделі, що відбивають ті самі процеси, що й у попередній, залишаються колишніми, а нові шукаються з умови існування відомої стаціонарної точки (натурні дані) і перевірки на стійкість.

Розширимо модель (6), додавши фазову координату B - кількість бактеріопланктону у водоймі, при цьому числові значення коефіцієнтів нової моделі, що відбивають ті ж процеси, що й у попередній, залишаються колишніми.

Розширимо модель (7), розбивши компоненту N₁ (мінеральний азот) на дві складові: N₁ - кількості амонійного азоту у водоймі і N₂ - кількість нітратного азоту у водосховищі. Числові значення коефіцієнтів нової моделі, що відбивають ті ж процеси, що й у попередній, залишаються колишніми, а інші знаходяться з умов існування позитивної стаціонарної точки, значення якої визначаються натурними даними.

На наступному етапі здійснюється перехід до змінних коефіцієнтів.

За своєю структурою ці моделі не відрізняються від моделей з постійними коефіцієнтами: вони описують ті ж процеси і використовують ті ж числові значення коефіцієнтів. Істотною відмінністю є тільки введення в деякі параметри в явному виді часу як аргументу.

Як змінні коефіцієнти розглядаються параметри: температура води, освітленість і скидання у водойму азоту. Всі інші коефіцієнти такі ж, як і у відповідних моделях з постійними коефіцієнтами.

За допомогою розробленого підходу створені і досліджені модель Балтійського моря для пошуку точок біфуркації екосистеми при скиданні токсикантів, модель кругообігу фосфору і кисню в Чорному морі для розгляду однієї з гіпотез поширення сірководневої зони в Чорному морі, пов'язаної з посиленням антропогенного забруднення моря, а також вивчення чутливості деяких параметрів моделі.

За умов посиленого розвитку атомної енергетики на перший план висуваються радіоекологічні питання, пов'язані з оцінкою відпрацьовування екосистемою і людиною забруднень від викиду в навколишнє середовище радіонуклідів. Основний принцип вибору рішень базується на оцінці альтернатив при гіпотетичних і реальних радіаційних аваріях. У межах цієї концепції представлені динамічні моделі колективних доз населення (внутрішня складова опромінення) у залежності від постачань того чи іншого продукту харчування, що входить до складу раціону певних вікових груп населення. В основу радіаційних динамічних моделей, розглянутих у даному розділі, покладена концептуальна модель переносу радіоактивних елементів з додаванням обставин виробництва, обробки, постачань харчових продуктів і розподілу населення по районах.

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

У розділі 6 розглядається автоматизація процесів багаторівневої візуалізації для виявлення цілісної картини взаємозв'язку і взаємозумовленості факторів, під впливом яких формуються негативні наслідки для навколишнього середовища, біоти і людини, спричинені аварією на ЧАЕС, а також іншим техногенним впливом.

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

Звідси випливає необхідність розробки проблемно-орієнтованих програмних продуктів для змоги доступу до постійно обновлюваних баз даних; виявлення системних зв'язків між радіохімічними факторами забруднення територій, візуалізації цих залежностей; розрахунку ризиків проживання на цих територіях, ризиків захворювання тією чи іншою хворобою і прийняття відповідних рішень.

Динамічні моделі процесів переносу радіонуклідів екосистемою, включаючи трофічні ланцюги, що ведуть до людини, реалізовані в рамках математичного пакета радіаційних прогнозів за допомогою об'єктно-орієнтованого підходу на Borland Delphi 5.0 (client/server suite).

Концептуальна схема є відображенням компонентів моделі і взаємозв'язків між ними. Оскільки в основу побудови схеми покладена камерна модель, то під компонентами розуміються камери, а під спрямованими взаємозв'язками - потоки переносу того чи іншого хімічного елемента. Кожен компонент разом зі зв'язками відбиває процес переносу в певне середовище (молоко, м'ясо і т.д.).

Проблема полягає в розв'язанні задачі надходження радіонуклідів в організм людини, тобто в одержанні числового значення дозового статусу.

Відповідно до розглянутої схеми у формуванні результуючого значення беруть участь усі компоненти (камери). Результуючі значення формуються в компоненти, що відбивають середовище шлунково-кишкового тракту (ШКТ) дорослої людини і дитини окремо.

Наступним кроком є математична реалізація даної концептуальної схеми моделі екосистеми. Модель є точковою, і тому процеси переносу потоків реалізовані за допомогою звичайних диференціальних рівнянь 1-го порядку (для одиничного забруднення 1 Кі/км²).

Програмна реалізація даної проблеми на сьогодні стосується тільки міграції радіонуклідів по окремих камерах моделі і знаходження значень дозового статусу людини.

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

Створений інтерфейс користувача дозволяє легко вносити виправлення в отримані дані, що автоматично призводить до відповідних змін на зображенні компонентів.

Підсумком досліджень, представлених у даному підрозділі, стали такі результати:

Розроблено початкову версію інтелектуального інтерфейсу системи візуалізаціїмедико-екологічних баз даних (на прикладі Житомирської області). Створено діаграму станів інтерфейсу користувача, що відбиває підмножини станів інтерфейсу для кожного етапу аналізу багатовимірної інформації.

Візуалізовані у виді графіків і діаграм дослідження, пов'язані з аналізом наявних екологічних баз даних Житомирського регіону по найбільш постраждалих від аварії на ЧАЕС  чотирьох районах - Лугинському, Народицькому, Овруцькому, Олевському (сільська місцевість, 183 села). Отримано близько сотні графічних інтерпретацій інформації, що відбиває системний підхід до оцінки екологічної ситуації в цілому.

Оброблено 57 карт по узагальнених факторах екологічного забруднення сільської місцевості Житомирської області за 1986-1995 рр., що подають розподіл показників як у виді діаграм, так і у виді ізоліній розподілу екологічного забруднення ґрунтів. Графічні картинки можна також представити у 2d- і 3d-вимірному зображенні.

У розділі 7 розглядаються математико-картографічні моделі візуалізації інтегральних індексів техногенного забруднення. Поняття математико-картографічного моделювання згідно з Тікуновим В.С. пов'язано з процесами створення, аналізу і перетворення картографічних моделей з метою одержання нових знань. Ізолінії в основному розглядаються як засіб картографічного зображення, а не як інструмент пізнання. У цьому розділі використовуються ізолінійні карти як інструмент візуалізації, виявлення і зіставлення різних сторін і властивостей складних просторово-розподілених об'єктів.

Характерними рисами ізолінійного моделювання є: охоплення процесів побудови і використання ізолінійних карт; використання апарату математичної статистики за допомогою обчислення і картографування показників середнього рівня, max, min і кореляційної залежності досліджуваних ознак.

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

Один з перспективних напрямків сучасного картування - розробка карт статистичних поверхонь, які можна віднести до області тривимірної картографії. Головною задачею цього напрямку є візуальне відтворення рельєфу такої статистичної поверхні за числовим значенням визначеної ознаки (у нашому випадку - рельєфу узагальненого фактора) з використанням зорових ілюзій. Найбільш відоме зображення статистичних поверхонь за допомогою ізоліній. Карти статистичних поверхонь, у більш широкому представленні (з урахуванням можливості застосування різних способів зображення), зможуть допомогти у розв'язанні важливої задачі візуалізації натурних даних.

Методи багатовимірного аналізу даних виявилися найбільш адекватним засобом для переходу від набору розрізнених показників забруднення до узагальнених екологічних індексів та їх картографічного зображення.