Радиоэкологическое моделирование динамики активности 131I в молоке в центральной части Мазовии по материалам «варшавского» сценария проекта МАГАТЭ EMRAS: реконструкция, верификация, оценки достоверности
Автор: Власов О.К., Звонова И.А., Щукина Н.В.
Рубрика: Научные статьи
Статья в выпуске: 2 т.33, 2024 года.
Бесплатный доступ
Работа посвящена исследованиям динамики транспорта радионуклидов 137Cs и 131I по трофической цепочке: «атмосфера - почва - растительность - молочный скот - молоко - организм человека» с использованием имитационной радиоэкологической модели на основе инструментальных данных «Варшавского» сценария проекта МАГАТЭ EMRAS, полученных после аварии на ЧАЭС. С использованием радиоэкологической модели выполнена реконструкция активностей 131I в молоке в районах Warsaw Area и Ostroleka Area центральной части Мазовии, где в апреле-мае 1986 г. проводилась его радиометрия. Цель исследования - верификация радиоэкологической модели, оценка достоверности и неопределённости молочного блока расчётной модели. Установлено, что последовательный учёт невязок расчётных и реконструированных инструментальных данных динамики активностей 131I в зелёном корме на фермах в округах Мазовии и недостатка запасов чистых кормов урожая 1985 г. приводят к последовательному улучшению согласия расчётных и инструментальных данных активностей 131I в молоке. Для всех вариантов учёта поправок отношения расчётных данных к инструментальным с пересчитанными активностями 131I в атмосфере района Ostroleka Area в среднем в 2,5-3 раза ближе к их идеальному значению, равному 1, чем в расчётах с инструментальными данными для облака. Прогностические свойства расчётной модели на этапе реконструкции удельных активностей 131I в молоке, оценённые в виде отношения расчётных данных к инструментальным данным, можно оценить в 1,05±2,0; величины этих отношений для корректировки оценок доз облучения щитовидной железы населения за счёт потребления загрязнённого молока оцениваются в 1,3±2,5.
Авария на чаэс, проект магатэ emras, «варшавский» сценарий, динамическая агрорадиоэкологическая модель, молочные фермы, погодные условия, обеспеченность ферм чистыми кормами, удельная активность 131i в молоке, верификация и оценки неопределённостей модели, ингаляционная и кормовая составляющая загрязнения молока, радиобиология, состояние окружающей среды
Короткий адрес: https://sciup.org/170204502
IDR: 170204502 | DOI: 10.21870/0131-3878-2024-33-2-17-33
Текст научной статьи Радиоэкологическое моделирование динамики активности 131I в молоке в центральной части Мазовии по материалам «варшавского» сценария проекта МАГАТЭ EMRAS: реконструкция, верификация, оценки достоверности
Настоящая работа является шестой в серии работ авторов [1-5], посвящённой исследованиям динамики транспорта радионуклидов 137Cs и 131I по трофической цепочке: «атмосфера – почва – растительность – молочный скот – молоко – организм человека» с использованием имитационной радиоэкологической модели [6, 7] на основе инструментальных данных «Варшавского» сценария проекта МАГАТЭ EMRAS [8-10], полученных после аварии на ЧАЭС.
В «Варшавском» сценарии для 29 населённых пунктов (НП), 27 молочных ферм, 24 молочных округов и 2 молочных районов Warsaw Area и Ostroleka Area представлена следующая база инструментальных данных за 1986 г.:
-
• динамика суточных удельных объёмных активностей 134,137Cs, аэрозольной, органической и газообразной форм 131I в атмосфере за период с 16.04.86 по 02.06.86 гг., измерения проводили в Центральной лаборатории радиологической защиты г. Варшава (молочный район Warsaw Area);
Власов О.К.* – зав. лаб., д.т.н.; Щукина Н.В. – ст. науч. сотр. МРНЦ им. А.Ф. Цыба – филиал ФГБУ «НМИЦ радиологии» Минздрава России. Звонова И.А. – гл. науч. сотр., д.т.н. ФБУН НИИРГ им. П.В. Рамзаева.
-
• метеоданные об осадках на 33 метеостанциях на время их фактических выпадений за период с 24.04.86 по 16.05.86 гг.;
-
• плотности выпадения 137Cs для НП и их распределения в молочных округах;
-
• годовой ход среднесуточных температур воздуха с 10.03.86 по 31.09.86 гг.;
-
• измерения удельной активности 131I в траве на территории метеостанции Warszawa Obserw Astr с 03.05.86 по 18.05.86 гг.;
-
• измерения удельной активности 131I в молоке с 29.04.86 по 07.06.86 гг. на 7 молочных фермах района Ostroleka Area и на 10 молочных фермах района Warsaw Area.
Созданная в предыдущих работах [1-5] база используется в настоящей работе для корректировки результатов их реконструкции, что должно обеспечить увеличение достоверности оценок динамики удельных активностей 131I в молоке, организме жителей и, в конечном счёте, при оценках доз внутреннего облучения их щитовидной железы (ЩЖ).
Цель исследования – верификация, оценки неопределённостей реконструкции динамики удельных активностей 131I в молоке в молочных округах центральной части Мазовии после аварии на ЧАЭС.
Материалы и методы
Исследования динамики удельных активностей в молоке в пунктах отбора его проб на радиометрию в молочных округах Мазовии проводятся с использованием блока КРС (крупный рогатый скот) радиоэкологической модели [7]. Блок представляет собой программный комплекс для расчётов метаболизма радионуклидов 131I в органах КРС путём численного решения системы линейных дифференциальных уравнений для многокамерной модели (рис. 1). Коэффициенты модели настроены на экспериментальные данные динамики активностей 131I в органах и каналах его выведения из организма КРС с удоем 12 л/сут [7]. Верификация модели и оценки её неопределённостей проводились на основе расчётов среднегеометрических величин и стандартных отклонений отношений расчётных данных к инструментальным.

Рис. 1. Модель организма «йодной» коровы [7].
Для демонстрации возможностей этого блока был выполнен сценарный расчёт с входными данными, представленными в табл. 1.
Пример динамики активностей 131I в органах и каналах выведения организма молочных коров в центральной части Мазовии с входными данными табл. 1 приведён на рис. 2.
Таблица 1
Параметры многокамерной модели
Параметр |
Mazovia, Warsaw Area |
Начало радиоактивных выпадений Время начала выпаса Режим содержания КРС Время перевода на чистый корм Удой молока Потребление зелёного корма |
2,4 сут 25 сут пастбищный 45 сут 11,6 л/сут 50 кг/сут |
Регион "Мазовия"

Рис. 2. Динамика активностей 131I в органах и каналах выведения КРС весной в год аварии на ЧАЭС в центральной части Мазовии [4].
Из особенностей данных на рис. 2 можно отметить, что в расчётном сценарии выпас КРС на пастбища начинается только спустя 25 сут после аварии, поэтому два первых максимума активности молока в период основных выпадений обусловлены ингаляцией, а второй, примерно в 10 раз больший, с максимумом на 28 сут, за счёт потребления загрязнённой травы пастбищ.
После перевода КРС на чистый корм на 45 сут сначала начинается быстрый спад активности 131I во всех органах и каналах его выведения из организма и через 2 сут активность молока уменьшается в 5 раз. Через 50 сут скорость уменьшения активности во всех органах остаётся постоянной в соответствии с радиоактивным распадом йода и медленной компонентой его скорости выведения из ЩЖ.
Исследования динамики активности 131I в молоке выполнялись с использованием следующих параметров созданной в работах [1-3, 5, 6] базой радиоэкологических данных для молочных ферм и округов в молочных районах Warsaw Area и Ostroleka Area, в которых производился отбор проб молока на радиометрию 131I:
-
• время начала выпаса КРС, равное 2 сут, рассчитано в агроклиматическом блоке имитационной модели [2] при выполнении двух условий: среднесуточная температура воздуха в течение нескольких суток должна быть выше 8 °С, плотность свежей биомассы пастбищной травы должна быть больше 0,4 кг/м2;
-
• величины эффективных осадков за период основных выпадений в ареале молочных ферм, согласованные с инструментальными данными динамики активности радионуклидов 137Cs в атмосфере и плотностями его выпадений [1-3];
-
• реконструированные инструментальные данные удельных активностей 1311 в траве культурных пастбищ, согласованные с базой данных эффективных осадков и активностью радионуклидов 131I в атмосфере [5];
-
• величины невязок (отношений расчётных и инструментальных данных удельных активностей 131I в траве культурных пастбищ) [5].
Все расчёты выполнялись по трём вариантам параметров атмосферы: прямой расчёт, однородное и неоднородное облако.
Прямой расчёт: для всех НП региона (часть территории, на которой были получены инструментальные данные динамики активностей радионуклидов в атмосфере) используются инструментальные данные динамики активности радионуклидов 137Cs, 13II и фаз его существования в атмосфере, метеоданные об осадках за период основных выпадений на ближайших к НП метеостанциях.
Однородное облако – неоднородные осадки: пространственная однородность удельных объёмных активностей радионуклидов в атмосфере над регионом – неоднородные эффективные осадки в НП (осадки в НП с постоянной интенсивностью в период основных выпадений, реконструированные по инструментальным данным активности 137Cs в атмосфере и плотностям его выпадения в НП).
Неоднородное облако – однородные осадки: в регион неоднородного облака входят все НП в районе ближайшей к ним метеостанции, для каждого НП региона – данные этой метеостанции об осадках за период основных выпадений, динамика удельных активностей радионуклидов в атмосфере над каждым НП нормируется на плотность выпадения 137Cs.
Результаты и обсуждение
Оценки долей загрязнённого корма в рационе КРС в округах центральной части Мазовии
По данным «Варшавского» сценария большинство частных хозяйств в округах района Ostroleka Area страдали от нехватки незагрязнённого корма прошлогоднего урожая и в рацион КРС частично входил свежезаготовленный загрязнённый зелёный корм. В среднем в Мазовии за период с 1.05.86 по 20.10.86 гг. зелёный корм входил в рацион КРС только на половине ферм.
В сценарии для округов в районах Warsaw Area и Ostroleka Area приведены данные о валовом производстве молока, структуре кормовых угодий, а также данные о плотностях выпадения 137Cs и осадках за период основных выпадений, рассчитанным по данным метеостанций.
Эта информация, а также расчётные данные зависимостей удельной активности травы культурных пастбищ от осадков за время основных выпадений, приведённые в работе [5], использовались для оценок относительных долей загрязнённого молока в округах (8^ак), в которых проводились измерения его удельной активности 131I. Оценки 8^а1к для к-го округа проводились по следующему соотношению:
S mUkj ^R r ) -
{^пЬ+^агу^т^ k -+Ikf .(я^1 '
[‘inh,j+I gr f,j(Rj )\
где j – модель расчёта: 1 – прямой расчёт, 2 – неоднородное облако, 3 – однородное облако; H ilary — доля ферм в округе, не имеющих запасов чистых кормов в мае 1986 г., принимается, что она соответствует средней доле загрязнённого корма в рационе КРС молочного округа; 8 k mllk j -доля радиоактивного молока, производимого в округе; R k - средние осадки за период основных выпадений, используемые в j -ой модели расчёта, мм; l k nh и i g r f - интегралы удельных активностей 131I в молоке за счёт ингаляции и загрязнённого зелёного корма соответственно, (кБк/кг)⋅сут.
Здесь:
I k. = I0 . J l inh,j l inh,j j
k С l k=2 k С l°
1J-1,3
,
lgrfj(RJ) Qgrfj(Rl) ^ (RO ^, где kClk - коэффициент пересчёта удельной активности 1311 в атмосфере для j-ой модели; kClJ0=2,1°lnh и Qgrfj оценивались по модели прямого расчёта для метеостанции Warszawa Obserw Astr с осадками за период основных выпадений 0,8 мм, равные 1,72, 14,4 и 696 (кБк/кг)⋅сут соответственно; Rk - осадки в к-ом молочном округе, мм.
Пример результатов расчёта радиоэкологических параметров для округа Warszawski Zachodniс с долей ферм, равной 30%, без запасов чистого корма к дате аварии приведён в табл. 2.
Радиоэкологические параметры для округа Warszawski Zachodni
Таблица 2
Модель расчёта |
kQ cloud |
R , мм |
q^, кБк/кг |
lmllk llnh , (кБк/кг)⋅сут |
tmllk l gr f , (кБк/кг)⋅сут |
8dlary о/ 8mllk , |
Прямой расчёт |
1,0 |
0,8 |
38 |
14,4 |
604 |
31,6 |
Неоднородное облако |
3,9 |
0,8 |
140,6 |
56 |
2230 |
31,7 |
Однородное облако |
1,0 |
7,8 |
33,6 |
14,4 |
534 |
31,8 |
Данные табл. 2 и аналогичные результаты для остальных округов показывают, что интеграл удельных активностей молока за счёт ингаляции в 35-40 раз меньше, чем за счёт зелёного корма. Соответственно, доли радиоактивного молока в них практически не зависят от модели расчёта и несущественно отличаются от доли ферм без запасов чистого корма к дате аварии.
Результаты расчётов для других округов, в которых производился отбор проб молока на радиометрию, приведённые в табл. 3, показывают, что в округе Ostroleka КРС был полностью обеспечен запасами чистого корма урожая прошлого года. Поэтому доля загрязнённого молока в этом округе равна доле ингаляционной составляющей в общей активности в молоке. По оценкам имитационной модели её величина равна 2%. В табл. 3 вместе с данными о величинах 8dlary и ^kmiik приведены также данные «Варшавского» сценария о структуре плотностей выпадения 137Cs в округах, метеоданные об осадках за период основных выпадений и число измерений динамики активностей 131I в пробах молока на пунктах его приёма.
Таблица 3
Радиоэкологические данные для молочных округов центральной части Мазовии
Молочный район |
Молочный округ |
Плотность выпадения 137Cs, кБк/м2 |
Осадки, мм |
^ diary Q/ Я т11к |
N изм |
|||
среднее |
минимум |
максимум |
R meteo |
R rec |
||||
Legionowski |
4,1 |
3 |
7 |
0 |
40,1 |
20/22 |
24 |
|
Minski |
4,5 |
4 |
5 |
0 |
40,2 |
30/32 |
4 |
|
Nowodworski Mazowiecki |
3,8 |
2 |
6 |
0 |
30,2 |
50/51 |
5 |
|
Warsaw Area |
Warszawski |
4,5 |
3 |
6 |
0 |
40,7 |
15/17 |
28 |
Warszawski Zachodni++ |
6,7 |
10,5 |
14 |
0 |
70,8 |
30/31 |
20+ |
|
Wolominski++ |
5,7 |
6 |
8 |
10 |
60 |
30/31 |
4 |
|
Falenty Duze diary*++ |
0,8 |
20,5 |
0,0/2,0** |
8 |
||||
Ostrolecki++ |
4,6 |
3 |
3,5 |
0 |
00,8 |
50/51 |
20+ |
|
Ostroleka |
4,6 |
3 |
3,5 |
0 |
0 |
0,0/2,0** |
20 |
|
Ostroleka |
Piski |
3,8 |
20,5 |
6 |
0 |
00,3 |
40/40 |
8 |
Area |
Przasnyski |
3,6 |
20,5 |
4 |
0 |
00,2 |
80/80 |
14 |
Szczycienski |
4,3 |
20,5 |
5 |
0 |
00,5 |
90/90,2 |
3 |
|
Zambrowski |
4,1 |
20,5 |
6,5 |
0 |
00,4 |
50/51 |
7 |
* – молочная ферма; + пересчёт от инструментальных данных фермы Ostroleka Diary округа Ostroleka по соотношению (4); ++ пересчёт от инструментальных данных фермы Warsaw Town Diary округа Warszawski по соотношению (4); R meteo и R rec – осадки за период основных выпадений 2,2-5,5 сут после аварии; ** – доля ингаляционного пути поступления 131I в организм молочных КРС; N изм – число последовательных измерений проб молока в одном месте.
Величины 8 mi l k использовали для корректировок расчётных Q^j и реконструированных инструментальных Q^ данных удельных активностей 1311 в молоке в к -ом молочном округе по инструментальным данным для округа k0 по соотношениям (4) и (5):
*'кк пСк як \ПМ№г+\ nC^(tW max[qCgr(t)]
ЛС,к ( _ л^'к Л+Л । як , як л С,к r.x
QcorJ(t) — QinhJ(t) + ^rad bgrf Q grfjktj’ (5)
где индексы M и C относятся к инструментальным и расчётным данным; j - модель расчёта; ( к и Q k cO r,j — прямые и скорректированные зависимости удельных активностей 1311 в молоке соответственно; q^ - расчётные данные динамики активностей 131I в пастбищной траве; Q k nh j и Q g r f j -расчётные зависимости ингаляционной и пищевой составляющей загрязнения молока соответственно; 8 g r f - поправки на невязки расчётных и реконструированных инструментальных данных удельных активностей 131I в пастбищной траве [6].
Отметим, что размах величин коэффициентов к _ ma^rjV]
Kco'l max^t)] (5)
для пересчёта инструментальных данных активностей 1311 в молоке K C or j по моделям прямого расчёта ( j =1) и однородного облака ( j =3) равен (0,95-1,05), для неоднородного облака ( j =2) равен (0,5-1,7).
Согласно данным «Варшавского» сценария, полное обеспечение чистыми кормами было только в округе Ostroleka района Ostroleka и на ферме Falenty Duze diary района Warsaw Area
(табл. 3). Наибольшая доля загрязнённого молока (90%) была в округе Szczycienski района Ostroleka Area. В среднем по всем округам с радиометрией 131I в молоке доля загрязнённого молока в районе Warsaw Area была в 2-3 раза меньше, чем в районе Ostroleka Area. С учётом практически совпадающих в этих районах параметров распределений плотностей выпадений 137Cs этот результат, возможно, может объяснить и особенность инструментальных данных активностей 131I в молоке, представленных на рис. 3.
Наиболее подробные многократные инструментальные данные динамики активностей 131I в молоке округов района Ostroleka Area были получены более 10 раз для фермы Ostroleka diary, Przasnyski diary, в районе Warsaw Area для фермы Warsaw town diary, округа Legionowski и фермы Falenty Duze diary (рис. 3).


б
Рис. 3. Инструментальные данные динамики удельных активностей 131I в молоке [8, 9]. а – молочный район Ostroleka; б – молочный район Warsaw.
Отметим существенную нерегулярность инструментальных данных для всех округов в области их максимальных значений и явную несогласованность между данными на рис. 2, в табл. 2 и на рис. 3 для округа Ostroleka.
Во-первых, инструментальные данные на рис. 3 для этого округа, единственного в районе Ostroleka Area со 100% обеспечением чистыми кормами (табл. 1), имеют одинаковые особенности с данными измерений во всех остальных округах этого района. Во всех этих округах обеспечение чистыми кормами было лишь частичным, и загрязнение молока в них было обусловлено в основном потреблением загрязнённой травы пастбищ. Во-вторых, активность131I в молоке на ферме Falenty Duze на время окончания основных выпадений была в 10 раз меньше, чем в молочных округах района Ostroleka Area.
Данные, полученные с помощью карты относительных долей ферм в молочных округах, не имеющих запасов чистых кормов в мае 1986 г., больше соответствуют форме инструментальных данных динамики активностей 131I в молоке для этого округа, полностью совпадающих с формой данных округа Przasnyski на рис. 3б с величиной 8d i ary , равной 80%. По картографическим данным «Варшавского» сценария в округе Ostroleka величина 8dlary составляла (0,5-0,8).
С учётом этих разногласий для округа Ostroleka проводился дополнительный тестовый расчёт с величиной доли загрязнённых кормов, равной 60%, как средней величины для территории всего района Ostroleka Area.
Результаты исследований динамики загрязнения молока и верификации радиоэкологической модели
Расчёты динамики активностей 131I в молоке проводились для округов, для которых в «Варшавском» сценарии имеются многократные инструментальные данные. На первом этапе для района Warsaw Area по всем трём вариантам входных данных об осадках и активностей 137Cs и 131I в атмосфере последовательно проводились следующие расчёты:
-
1. расчёт без учёта поправок (рис. 4а – прямой расчёт);
-
2. с учётом невязок расчётных и реконструированных инструментальных активностей 131I в зелёном корме (рис. 4б);
-
3. с учётом доли чистого зелёного корма в рационе КРС в молочных округах (расчёт с учётом поправок на загрязнение и относительных долей ферм (%) в молочных округах, не имеющих запасов чистых кормов, рис. 4в);
-
4. c учётом поправок пп. 2 и 3 (рис. 4г).
Результаты расчётов и инструментальные данные удельных активностей 131I в молоке представлены для фермы Falenty Duze и округа Legionowski без учёта поправок на рис. 4а и на рис. 4б для этого же округа с учётом поправок на загрязнение и долю загрязнённого корма в рационе КРС.
Реконструкция ингаляционной составляющей поступления 131I в молоко, выполненная на примере инструментальных данных для фермы Falenty Duze, в рацион КРС которого входил только чистый корм урожая прошлого года, показывает, что результаты расчётов по модели однородного облака вполне удовлетворительно согласуются с инструментальными данными как по статистике их отношений, так и по их форме, отражающей детали динамики объёмных активностей 131I в атмосфере.
Данные на рис. 4 показывают, что модели прямого расчёта и однородного облака дают практически одинаковые расчётные результаты и имеют в 2,5 раза лучшее согласие с инструментальными данными, чем модель неоднородного облака. Причём, отношение расчётных к инструментальным данным для всей их серии ближе к 1, чем в области их максимальных значений. Вероятно, это связано с существенной нерегулярностью инструментальных данных вблизи их максимальных значений.
Исследования динамики удельной активности 131I в молоке округа Ostroleka проводились и для района Warsaw Area с использованием трёх моделей расчёта и двух вариантов динамики активностей 137Cs и 131I в атмосфере:
-
• прямые инструментальные данные «Варшавского» сценария, полученные в г. Варшаве;
-
• реконструированные инструментальные данные, равные прямым инструментальным данным, умноженным на величину отношения расчётной плотности «сухих» выпадений 137Cs, оценённой в работе [3] по прямым инструментальным данным его активности в атмосфере, к минимальной инструментальной плотности его выпадений.
Облако = "прямые данные" County = "County legionowski"
* * * измерение
^^^^^^— измерение, интерполяция Л---А--А модель прямого расчета
■---■---■ модель неоднородного облака
----------модель однородного облака в---в---&■ Falenty. измерение
-----------Falenty. прямой расчет, однородное облако ........ Falenty. неоднородное обдако
б
Облако = "все данные прямые"
поправка на долю чистого корма
в
County = "County legionowski" Облако = "все данные прямые" поправка на загрязнение и долю чистого корма
Рис. 4. Динамика удельных активностей 131I в молоке в округе Legionowski района Warsaw Area и на ферме Falenty Duze.
а – расчёт без учёта поправок; б – расчёт с учётом поправок на загрязнение корма; в – расчёт с учётом доли чистого корма; г – расчёт с учётом поправок на загрязнение и долю чистого корма.
Инструментальные и расчётные временные зависимости удельных активностей 131I в молоке без учёта и с учётом поправок на загрязнение и долю чистого корма представлены на рис. 5, 6 для округов Ostroleka и Przasnyski соответственно, а также на рис. 7 в виде гистограмм среднегеометрических отношений расчёт/измерение для всех их серий инструментальных данных и в области их максимальных значений.
Сопоставление данных на рис. 5а и 5б показывает, что при исследованиях динамики загрязнения травы пастбищ правомерность использования варианта с пересчётом удельных объёмных активностей 131I в атмосфере над районом Ostroleka, рассмотренного нами в работе [2], при создании в работе [5] и использовании взаимно согласованной базы входных данных: удельные объёмные активности 137Cs и 131I в атмосфере – осадки – плотности выпадения 137Cs. Это подтверждается лучшим согласием расчётных и инструментальных данных активностей 131I в молоке в округах района Ostroleka, чем в варианте с их прямыми инструментальными данными.


б

в

г
Рис. 5. Динамика удельных активностей 131I в молоке в округе Ostroleka.
а – расчёт без учёта поправок, облако – прямые данные; б – расчёт без учёта поправок, облако – данные с пересчётом; в – расчёт с учётом поправок на загрязнение и долю чистого корма, облако – прямые данные; г – расчёт с учётом поправок на загрязнение и долю чистого корма, облако – данные с пересчётом.
Отметим, что инструментальные данные для округа Przasnyski на рис. 6 имеют более регулярный вид без провала в области 12-18 сут, как в округе Osrtoleka на рис. 5.


Рис. 6. Динамика удельных активностей 131I в молоке в округе Przasnyski. а – облако – прямые данные; б – облако – данные с пересчётом.

Рис. 7. Гистограммы среднегеометрических отношений расчётных к инструментальным данным активностей 131I в молоке для округов Ostroleka и Przasnyski.
* – без учёта поправок; ** – с учётом поправок.
Анализ сопоставления расчётных и инструментальных данных показал, что расчётные данные для Warsaw Area лучше воспроизводят инструментальные данные после 30 сут после аварии, чем для Osrtoleka Area, и хуже на времена, меньшие 20 сут. В диапазоне 20-30 сут для обоих вариантов имеется практически вполне удовлетворительное согласие расчётных и инструментальных данных.
Улучшение сходимости расчётных и инструментальных данных при последовательном учёте поправок характерно как для всех округов обоих районов, так и для этих районов в целом (рис. 8).

Warsaw Area, все данные

б


Рис. 8. Отношение реконструированных активностей 131I в молоке к инструментальным данным. а, б – для всей области измерений; в, г – область максимальных значений.
Отметим основные особенности величин невязок расчётных и инструментальных данных на рис. 8:
-
• во всех трёх вариантах расчёта последовательный учёт поправок приводит к улучшению сходимости величин невязок к 1;
-
• для всех вариантов расчётов отношения расчётных к инструментальным данным с пересчитанными активностями 131I в атмосфере района Ostroleka Area в среднем на 0,25-0,3 ближе к их идеальному значению, равному 1, чем в расчётах с прямыми данными облака;
-
• величины отношений для серии инструментальных данных в области их максимальных значений больше на 30%, чем для всей их серии;
-
• во всех вариантах учёта поправок невязки для района Warsaw Area на 0,2-0,25 ближе к 1, чем для района Ostroleka Area.
Стандартные среднегеометрические отклонения невязок в виде минимальных/сред-них/максимальных значений равны для района Ostroleka area 1,4/1,8/2,4, для района Warsaw Area – 1,6/2,1/2,5.
Таким образом, данные на рис. 8 подтверждают ранее сделанное в работе [2] предположение о том, что удельная активность 131I в атмосфере района Ostroleka Area больше, чем данные её измерений в г. Варшава района Warsaw Area. Так, по оценкам, сделанным на основании сопоставления минимальных плотностей выпадения 137Cs в этих районах, величина этого превышения равна 2,7 раза, а по данным рис. 8 она равна в среднем 2,5-3 раза.
Из физики переноса радионуклидов по цепочке: «растительность – организм КРС – молоко
– организм человека» следует, что максимальные активности радионуклидов в организме человека за счёт их потребления с молоком и зеленью будут прямо пропорциональны их максимальным активностям. Соответственно, пищевая составляющая доз внутреннего облучения ЩЖ и всего тела человека будет так же прямо пропорциональны этим величинам. С учётом этого невязки в виде отношения расчётных и инструментальных данных удельных активностей радионуклидов в молоке и зелени могут использоваться для корректировки расчётных данных инте- грала их поступления в организм человека по следующему соотношению:
I cor _ n . fC , л-1 , D rC . л-1
1 *gr Igr лgr + rmilk Imilk лmllk, где Icor - скорректированный суммарный интеграл поступления активности радионуклидов с зеленью и молоком в организм человека, кБк; Pgr и Pmilk - суточное потребление зелени и молока соответственно, кг/сут; Icgr и imilk - расчётные интегралы удельных активностей радионуклидов в зелени и молоке соответственно, кБк/кг⋅сут;
л = ЫВ1 л= [ (^ dm^ g [g™)^ , '
- max '-max
Лgr и Лmllk - невязки расчётных qdlr и инструментальных qlnstr данных удельных активностей радионуклидов в зелени gr и в молоке milk в области максимальных значений инструменталь ных данных -max.
В целом можно принять, что прогностические свойства расчётной модели с учётом всех поправок, оценённых в виде стандартных среднегеометрических значений и среднегеометрических отклонений отношения расчётных данных к инструментальным данным активностей молока, можно оценить как 1,05±2,0, а величины этих отношений для корректировки оценок доз облучения ЩЖ населения как Л gr=0,75±1,7; Л mllk =1,3±2,5.
Заключение
Выполнена реконструкция динамики активностей 131I в молоке в округах центральной части Мазовии, Польша. Оценена достоверность и неопределённость блока динамической радиоэкологической модели.
Правомерность пересчёта параметров загрязнения атмосферы над районом Ostroleka от инструментальных данных в Варшаве, предложенная нами при создании взаимно согласованной базы входных данных: удельные активности 137Cs и 131I в атмосфере – осадки в период основных выпадений – плотности выпадения 137Cs на местность, и при исследованиях динамики загрязнения травы пастбищ, подтверждается более лучшим согласием расчётных и инструментальных данных динамики активностей 131I в молоке, чем в расчётах без такого пересчёта.
Основные численные результаты выполненной работы заключаются в следующем:
-
• получено, что последовательный учёт невязок расчётных и реконструированных инструментальных данных динамики активностей 131I в пастбищной траве, зелёном корме КРС на фермах в округах Мазовии и учёт недостатка запасов чистых кормов урожая 1985 г. приводит к
улучшению согласия расчётных и инструментальных данных активностей 131I в молоке для всех округов с радиометрией молока;
-
• для всех вариантов учёта поправок отношения расчётных к инструментальным данным с пересчитанными активностями 131I в атмосфере района Ostroleka Area в среднем в 2,5-3 раза ближе к их идеальному значению, равному 1, чем в расчётах с прямыми данными облака;
-
• вследствие значительной нерегулярности инструментальных данных в области их максимальных значений в период 5-15 сут после аварии, величины отношений инструментальных данных за этот период области больше на 30%, чем для всей их серии;
-
• невязки для района Warsaw Area в 2-2,5 раза больше, чем для района Ostroleka Area;
-
• стандартные среднегеометрические отклонения невязок в виде минимальных/сред-них/максимальных значений равны для района Ostroleka area 1,4/1,8/2,4 для района Warsaw Area – 1,6/2,1/2,5;
-
• прогностические свойства расчётной модели на этапе реконструкции удельных активностей 131I в молоке, оценённые в виде отношения расчётных к инструментальным данным, можно оценить в 1,05±2,0, а величины этих отношений для корректировки оценок доз облучения ЩЖ населения за счёт потребления загрязнённого молока – в 1,3±2,5 .
Список литературы Радиоэкологическое моделирование динамики активности 131I в молоке в центральной части Мазовии по материалам «варшавского» сценария проекта МАГАТЭ EMRAS: реконструкция, верификация, оценки достоверности
- Власов О.К., Краевский П., Бартускова М., Малатова И., Щукина Н.В., Чекин С.Ю., Туманов К.А., Звонова И.А. Создание и верификация взаимосогласованной базы входной информации имитационной модели транспорта радионуклидов 137Cs и 131I по пищевой цепочке на инструментальных радиоэкологических данных (по материалам "Пражского" и "Варшавского" сценариев проекта МАГАТЭ EMRAS) //Радиация и риск. 2019. T. 28, № 3. С. 5-23. EDN: DZFBCA
- Власов О.К., Краевский П., Звонова И.А., Щукина Н.В., Чекин С.Ю., Туманов К.А. Реконструкция параметров радиоэкологической модели по результатам мониторинга радиационной обстановки в Мазовии после аварии на ЧАЭС (по материалам "Варшавского" сценария проекта МАГАТЭ EMRAS) //Радиационная гигиена. 2020. Т. 13, № 4. С. 26-33. EDN: TOVSYL
- Власов О.К., Краевский П., Звонова И.А., Малатова И., Бартускова М., Щукина Н.В., Чекин С.Ю., Туманов К.А. Верификация радиоэкологической модели на инструментальных данных удельной активности 131I в траве в Мазовии и Богемии после аварии на ЧАЭС (по материалам "Пражского" и "Варшавского" сценариев проекта МАГАТЭ EMRAS) //Радиационная гигиена. 2019. Т. 12, № 4. С. 6-17. EDN: YYOXSU
- Власов О.К., Звонова И.А., Краевский П., Щукина Н.В., Чекин С.Ю., Туманов К.А. Радиоэкологическое моделирование динамики активности 131I в травяной растительности разных видов в год аварии на ЧАЭС //Известия вузов. Ядерная энергетика. 2020. № 4. С. 106-117. EDN: XRGIEA
- Власов О.К., Звонова И.А., Щукина Н.В., Чекин С.Ю. Радиоэкологическое моделирование динамики активности 131I в пастбищной растительности Мазовии в год аварии на ЧАЭС: реконструкция, верификация, оценки достоверности //Известия вузов. Ядерная энергетика. 2023. № 2. С. 107-120. EDN: WRBCEF
- Власов О.К. Радиоэкологическая модель транспорта радионуклидов йода и цезия по пищевым цепочкам после радиационных аварий с выбросом в атмосферу для исследований закономерностей формирования доз внутреннего облучения населения. Часть 1. Описание, постановка и свойства агроклиматического блока модели //Радиация и риск. 2013. Т. 22, № 2. С. 16-34. EDN: QZKJBX
- Власов О.К. Радиоэкологическая модель транспорта радионуклидов йода и цезия по пищевым цепочкам после радиационных аварий с их выбросом в атмосферу для исследований закономерностей формирования доз внутреннего облучения населения. Часть 2. Описание, постановка и свойства блоков метаболизма радионуклидов 131I и 137Cs в органах крупного молочного скота //Радиация и риск. 2013. Т. 22, № 3. С. 21-39. EDN: ROTSXD
- Environmental Modelling for Radiation Safety (EMRAS): a summary report of the results of the EMRAS programme (2003-2007). IAEA-TECDOC-1678. Vienna: IAEA, 2012. 60 p.
- Krajewski P., Ammann M., Bartusková M., Duffa C., Filistovic V., Homma T., Kanyár B., Malátová I., Nedveckaite T., Simon S., Vlasov O., Webbe-Wood D., Zvonova I. Validation of environmental transfer models and assessment of the effectiveness of countermeasures using data on 131I releases from Chernobyl //Appl. Radiat. Isot. 2008. V. 66, N 11. P.1730-1735. EDN: LLDEIH
- Bartusková M., Malátová I., Berkovskyy V., Krajewski P., Ammann M., Filistovic V., Homma T., Horyna J., Kanyár B., Nedveckaite T., Vlasov O., Zvonova I. Radioecological assessments of the iodine working group of IAEA's EMRAS programme: presentation of input data and analysis of results of the Prague scenario //Radioprotection. 2009. V. 44, N 5. P. 295-300. EDN: OZLLZX