Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Радіобіологія фул вершин (передмовалесс, вступлесс етс. едишн)

.pdf
Скачиваний:
101
Добавлен:
07.03.2016
Размер:
5.97 Mб
Скачать

330

Сьогодні дуже багато методів та способів математичного моделювання. В даній роботі був використаний метод камерних моделей. Цей метод передбачає: розбиття ланцюга подій системи «середовище» на сукупність ланок (камер); встановлення найбільш суттєвих зв’язків, котрі описують перенесення радіонуклідів окремими камерами об’єднують компоненти в одну систему; створення математичного апарату, який моделює динаміку поводження радіонуклідів в системі [2,3].

В якості програмного забезпечення використовувався математичний продукт – MAPLE V.

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

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

значення коефіцієнтів, які отримані з натурних досліджень та літературних даних):

aij Значення коефіцієнта

a12 0,01 a 23 0,3 a34 0,15 а 45 0,1 а56 0,05 а67 0,2 а78 0,4 а79 0,2 а69 0,1 а87 0,15 а89 0,01

331

 

Вершина

x(t)

1

 

 

гори

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а12

 

 

 

 

Льодовик

 

2

 

 

 

 

 

y(t)

 

 

 

 

 

 

 

 

 

 

 

a23

 

 

 

Альпійський луг

3

 

 

 

 

z(t)

 

 

 

 

 

 

 

 

 

 

 

a34

 

 

 

 

Чагарники

4

 

 

 

 

k(t)

 

 

 

 

 

 

 

 

 

 

 

a45

 

 

 

 

Ліс

 

l(t)

5

 

 

 

 

 

 

 

 

 

 

 

 

 

a56

 

 

 

 

 

Пасовище

6

 

 

 

 

 

m(t)

 

 

 

 

 

 

 

 

 

 

 

a67

a69

 

Агрозона

 

7

 

 

 

 

 

 

 

n(t)

a79

 

 

 

 

 

 

 

 

a87

 

 

a78

 

Ріка p(t)

 

Людина

 

a89

r(t)

8

 

9

Таблиця11.3.

Значення

Рис.11.9. Блок-камерна модель

коефіцієнтів

переходу

гірської

екосистеми,

де

aij

радіонуклідів із камери в камеру

швидкості переходу радіонукліда.

 

 

Для дослідження

міграції радіонукліда

137Cs складемо

систему

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

332

dx

 

 

 

0,04x(t),

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dy

 

 

 

0,01x(t) 0,33y(t),

 

dt

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

0,3y(t) 0,18z(t),

 

 

 

 

dt

 

 

dk

 

 

 

0,15z(t) 0,13k (t),

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

dl

 

 

 

 

0,1k(t) 0,08l(t),

 

 

 

dt

 

 

 

 

 

dm

0,05l(t) 0,33m(t),

 

 

 

 

 

 

 

 

dt

 

 

 

 

dn

 

 

 

0,2m(t) 0,15 p(t) 0,63n(t),

 

dt

 

 

 

 

 

 

 

 

 

dp

 

 

 

0,4n(t) 0,28 p(t),

dt

 

 

 

dr

 

 

0,01p(t) 0,2n(t) 0,1m(t). (11.11)

 

 

 

 

 

 

 

 

dt

 

 

 

 

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

1)забруднення радіонуклідом відбулося безпосередньо на вершині гори;

2)радіонуклідне забруднення сталося одразу на кількох складових: вершині гори, льодовику та альпійському лузі.

У першому випадку початкові умови системи диференціальних рівнянь мають вигляд:

x(0) 100, y(0) 0, z(0) 0, k(0) 0, l(0) 0, m(0) 0, n(0) 0, p(0) 0, r(0) 0.

А у другому, коли забруднення почалося одночасно із забруднення вершини гори, льодовика та лугу, мають наступний вигляд:

x(0) 0,33, y(0) 0,33, z(0) 0,33, k(0) 0, l(0) 0, m(0) 0, n(0) 0, p(0) 0, r(0) 0.

333

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

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

Таблиця 11.4 Накопичення радіонуклідів у камерах для двох випадків викиду . 1-випадок коли викид забруднив тільки вершину

гори, 2 випадок -

коли викид

забруднив більшу частину гори

(результати в лапках)

 

 

 

 

 

Камери

 

Максимальна

Час (рік)

 

 

активність

 

 

 

 

радіонуклідів

в

 

 

 

камері (%)

 

 

Трава

 

3.2

(3)

12

(36)

 

 

Чагарники

3.2

(10)

24

(21)

 

 

Ліс

 

3.2

(16)

36

(20)

 

 

Пасовище

 

0.5

(2.4)

40

(24)

 

 

Агрозона

 

0.22

(1.05)

46

(28)

 

 

Ріка

 

0.32

(1.42)

50

(32)

 

 

Людина

 

5.4

(17)

108

(100)

 

 

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

забруднення, тобто забруднення великої частини гори: це підвищує величину прогнозованої колективної дози для населення.

Якщо порівняти результати моделювання розподілу радіонукліда Cs 137 у схиловій екосистемі (для середніх коефіцієнтів переходу) (рис.11.6,11.7), де для людини максимальна активність радіонукліда становить 10% на 50 рік, то все ж таки найбільш небезпечним для людини буде рівномірне забруднення гірської екосистеми (рис.11.10,

рис.11.11):

334

r

 

 

 

r

 

 

 

 

 

20

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

15

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

3

 

 

 

 

 

5

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

0

20 40 60

80 t

 

 

 

 

 

 

 

 

 

 

0

50

100

150

t

 

Рис.11.10.Накопичення

Рис.11.11.Накопичення

 

колективної

дози

від

колективної

 

дози

від

радіонукліда для населення

радіонукліда

для

населення у

у схиловій екосистемі(викид

гірській

екосистемі

(викид

відбувся на 3-х верхніх

радіонукліда

відбувся

на

частинах схилу), де r

- доля

вершині гори), де r - доля (%) від

(%)

від запасу радіонукліда

запасу

 

радіонукліда

в

векосистемі, яка екосистемі, яка реалізується

реалізується

через

через колективну дозу, t - час (

колективну дозу, t

- час

роки).

(роки).

 

 

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

r 18 16 14 12 10 8 6 4 2

0

0 20 40 60 80 100 120 t

335

Рис.11.12. Накопичення колективної дози від радіонукліда для населення у гірській екосистемі (рівномірне забруднення), де r - доля (%) від запасу радіонукліда в екосистемі, яка реалізується через колективну дозу, t - час ( роки).

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

1.Запропонована в даній роботі модель працює і придатна для моделювання природних процесів у гірських екосистемах.

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

3.Аналіз результатів моделювання показує, що відбувається інтенсивний перерозподіл радіонукліда в даній гірський екосистемі

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

5.Положення пікових значень вмісту радіонукліда по камерам зсувається по часу. Чим нижче розташована камера, тим дальше формується пікове значення.

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

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

значення aij .

8. Показано, що відбувається накопичення радіонукліда в популяції населення, що використовує дану гірську екосистему. Розрахунки по моделі вказують на те, що частина радіонукліда, яка формує колективну дозу у людей зростає у часі і досягає помітних значень 5,4% у першому випадку та 17% у другому. Це означає зростання прогнозованої колективної дози у разі рівномірного забруднення значної частини гори.

Моделювання радіоекологічних процесів методом камерних моделей при оцінці екологічного стану на прикладі села Галузія Маневицького району Волинській області. В даному розділі використані, дані натурних експедиційних досліджень вмісту 137Cs у ґрунті, рослинах, молоці, коровах (прижиттєво) та у людей (вимірювання на лічильниках випромінювань людини – ЛВЛ), опитування населення на територіях вибраних сіл Волинської області,

336

віднесених до другої та третьої зон; здійснене математичне моделювання з використанням камерних моделей. Численними дослідженнями було встановлено, що на території с. Галузія формуються значні дозові навантаження у людей внаслідок великих значень коефіцієнтів переходу в системі ґрунт-рослини і високих рівнів забруднення сіна, молока, м`яса, лісових продуктів – грибів та ягід.Наші попередні оцінки показали, що для даної території є три основні напрями надходження радіонуклідів 137Cs до людини, які більш докладно розглянуті нижче. Перший з них – через урочища (пасовища та сіножаті), що є кормовою базою для молочної і м`ясної худоби. Цей напрям формує в середньому до 70 % усієї дози. Другий напрям – надходження через лісові продукти – гриби та лісові ягоди, що складає в середньому 20 % загальної дози. Третій – присадибна ділянка (город) забезпечує в середньому 10 % від загальної дози для жителів цих населених пунктів.

Блок-схема камерної моделі села представлена на рис.11.13, а схема розрахунку камерної моделі– на рис. 11.14. Було розглянуто камерну модель потоків радіонуклідів у селі Галузія Маневицького району Волинської області та проведений розрахунок результатів даної моделі. За отриманим рішенням побудовані графіки виносу радіонуклідів з усіх полігонів (пасовищ, лісу і городу) та проведено аналіз отриманих результатів.

337

1 пасовище:

S = 80 га; врожай – 0,1 кг/м2 Густина забруднення ґрунту – 0,7 Кі/км2

КП = 100 – 300

Корів – 40

2 пасовище:

S = 130 га; урожай – 0,3 кг/м2

Густина забруднення ґрунту – 1,3 Кі/км2

КП = 50 – 200

Корів - 70

Село Галузія Населення: 1450, в тому числі: діти – 380,

пенсіонери – 610, робітники – 460 Густина забруднення городніх культур – 0,6 Кі/км2

3пасовище (сіянка):IV S = 45 га; врожай – 1,2 кг/м2 Густина забруднення ґрунту – 1,4 Кі/км2 КП = 120 – 270,Корів – 28

4 пасовище: V

S = 64 га; врожай – 0,4 кг/м2 Густина забруднення ґрунту – 1,7 Кі/км2 КП = 80 – 340 Корів – 44

Ліс

Площа = 240 га Врожай грибів – 7 кг/км2;

густина забруднення – 2,1 Кі/кг; КП = 150 – 480. Врожай ягід (чорниці) – 30 кг/км2 (чорничник); густина забруднення –700-1100 Бк/кг; КП = 50-70.

Деревина, густина забруднення – 10-30 Бк/кг,КП = 0,5 – 2

Рис. 11.13. Принципова блок-схема основних складових екосистеми с. Галузія Маневицького району Волинської області.

338

Ґрунт

а12

Трава

а23

Корова

а34

Молоко

а47

 

Діти

1

x(t)

2

y(t)

3

z(t)

4

k(t)

 

7 l(t)

 

 

 

 

 

а12

 

 

а57

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М’ясо

а48

 

Пенсіонери

 

 

 

 

 

 

 

5

 

8

 

 

 

 

 

а36

 

m(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o(t)

 

 

а59

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Відходи

 

 

Робітники

 

 

 

 

 

 

 

6

 

 

9

 

 

 

 

 

 

 

 

а410

 

 

 

 

 

 

 

 

 

 

r(t)

 

 

 

 

 

 

 

 

 

 

 

 

Експорт

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

Рис. 11.14. Структурована блок-схема екосистеми, що досліджується.

Дані, що отримані за моделями, добре корелюються з даними натурних досліджень.В даному розділі роботи зроблено такі висновки: 1. Розроблено та проаналізовано камерну модель реальної екосистеми села Галузія Маневицького району Волинської області, що попало під вплив аварії на Чорнобильській АЕС.

2. В моделі враховані всі основні потоки радіонуклідів -137Cs.

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

4.Отримано розрахункові дані щодо динаміки потоків радіонуклідів по камерах досліджуваної екосистеми.

5.Проведено оцінку та спрогнозовано очікуване забруднення радіонуклідами основних продуктів харчування – молока, м’яса, картоплі та овочів, а також дозових навантажень людей протягом певного часу після аварії на ЧАЕС.

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

339

забруднених пасовищах. Доза від вживання молока становить 4060 % від усієї дози для людей.

7. Показано, що помітні рівні забруднення молока сформувались на цій території не відразу після аварії, а тільки у 1992-1994 роках. 8.Підтверджено, що важливим компонентом дозоутворення у цьому регіоні є використання лісових продуктів – грибів та ягід.

9.Показано, що значна частина колективної дози не формується локально у даному селі, а експортується на інші території через вивіз молока, м’яса та лісових продуктів.

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

Моделювання радіоекологічних процесів методом камерних моделей при оцінці екологічного стану на прикладі села Коцюбинчики Чортківського району Тернопільської області.

Аналогічно до с. Галузія, здійснено розрахунки і по с. Коцюбинчики Тернопільської області. Це дозволило зробити порівняльну характеристику двох досліджуваних полігонів.

 

 

 

 

2 пасовище:

 

 

 

 

 

 

S = 29 га;

 

 

 

 

 

 

Врожай – 0,2 кг/м2

 

 

 

1 пасовище (сіянка):

 

 

3 пасовище:

 

 

Густина забруднення

 

 

S = 84 га;

 

 

S = 90 га;

 

 

ґрунту – 1,29 Кі/км2

 

 

Врожай – 1,4 кг/м2

 

 

Врожай – 0,5 кг/м2

 

 

КП = 3 – 8,Корів – 20

 

 

Густина забруднення

 

 

Густина забруднення

 

 

 

 

 

 

 

ґрунту – 1,27 Кі/км2

 

ІІІ

 

ґрунту – 1,34 Кі/км2

 

КП = 1 – 5,Корів – 30

 

 

 

КП = 5 – 11, Корів – 35

 

 

Село Коцюбинчики

 

 

 

 

 

 

 

 

 

 

 

 

ІV

 

 

 

 

Населення – 850,

 

 

 

 

 

в тому числі діти – 61,

 

 

 

V

робітники – 330,

 

 

 

пенсіонери – 459

 

 

 

 

 

 

 

 

 

 

 

 

Густина забруднення

 

 

 

4 пасовище:

 

ґрунту – 1,2 Кі/км2

 

5 пасовище (сіянка):

 

S = 40 га; врожай – 0,3

 

 

 

S = 32,5 га; врожай – 1,3

 

кг/м2

 

 

 

кг/м2 Густина

 

Густина забруднення

 

 

 

забруднення;ґрунту –

 

ґрунту – 1,19 Кі/км2

 

 

 

1,30 Кі/км2; КП = 4 – 9

 

КП = 0,5 – 3,Корів – 28

 

 

 

Корів – 22