CC BY
23
Оценка влияния обжига на радиоактивность глинистых горных пород и
материалов
П.Э. Соколов1, С.А. Сентенберг2
'Волгоградский государственный технический университет, 2Волгоградский технологический колледж
Аннотация: Проведены исследования по оценке влияния обжига на удельные активности и эффективные удельные активности глинистых горных пород. Корреляционный анализ позволил установить зависимости между исследованными характеристиками. Это позволило построить модель полиномиальной регрессии для удельных активностей и эффективной удельной активности естественных радионуклидов. Результаты исследований позволяют, используя полученные уравнения регрессии и технологические режимы производства уменьшить эффективную удельную активность и, следовательно, дозы облучения населения от этих материалов.
Ключевые слова: естественная радиоактивность, обжиг, глинистые горные породы, радионуклид, корреляционный анализ, полиномиальная регрессия.
В настоящее время признана гипотеза об отсутствии порога действия радиации. Считается, что любые, сколь угодно малые, дополнительные дозы радиации оказывают негативное воздействие на живой организм [1].
Значительный вклад в дозу облучения населения вносят природные источники ионизирующего излучения, «привносимые» в жилые, общественные и производственные здания строительными материалами [2].
Концентрация естественных радионуклидов (ЕРН) в горных породах обуславливает радиоактивность готовых строительных материалов.
Литературные источники содержат значительное количество данных об удельных активностях ЕРН как для горных пород, применяемых для производства строительных материалов, так и для готовых материалов [3-5]. Наряду с этим имеют место исследования влияния тепловой обработки на коэффициент эманирования радона [6-8].
Исследованию связи между различными материалами, средами и радионуклидами посвящен ряд публикаций [9, 10]. Однако вопрос влияния тепловой обработки на удельные активности ЕРН и эффективную удельную активность глинистых горных пород и материалов остается малоизученным.
Поэтому изучение влияния обжига на радиоактивность глинистых горных пород, а также поиск путей и возможностей ее снижения при производстве строительных материалов представляет интерес. С этой целью были проведены исследования влияния тепловой обработки на радиоактивность глинистых горных пород и материалов.
Исследуемые материалы представлены тремя видами глин из различных месторождений Волгоградской области. Все исследованные глины применяются для производства различных видов строительных материалов (керамического кирпича, портландцемента, керамзита). Применение этих строительных материалов обязывает соблюдать все нормы и требования безопасности. Одной из нормируемых величин является эффективная удельная активность ЕРН, а как следствие удельные активности нормируемых радионуклидов - 40К, 226Яа и 232ТЪ (ГОСТ 30108-94, НРБ-99/2009).
Отобранные пробы исследованных материалов отбирались по стандартной методике [11]. Определение удельных активностей ЕРН производилось на сцинтилляционном гамма-спектрометре. Полученные значения удельных активностей использовались для расчета эффективной удельной активности ЕРН по формуле (ГОСТ 30108-94):
Афф = 0,09Ак + ЛКа +1,3Лп , (1)
где АК - удельная активность 40К, Бк/кг; АКа - удельная активность 226Яа, Бк/кг; А-^ - удельная активность 232Т^ Бк/кг.
Далее пробы подвергались постадийной тепловой обработке, с шагом 150°С, в интервале температур от 20 до 1500°С. После каждой стадии обжига для каждой пробы определялись удельные активности и рассчитывалась эффективная удельная активность ЕРН. Полученные результаты исследований приведены в таблице № 1.
Таблица № 1
Удельные активности и эффективные удельные активности ЕРН (Бк/кг) при
различных температурах обжига (°С)
Показатель Температура обжига
20 150 300 450 600 750
Глина 1
40К 1912,0 1208,0 1502,0 1817,0 1999,0 1164,0
226Яа 60,4 39,2 28,6 47,8 89,4 42,7
232ТЪ 75,5 74,9 111,2 64,1 103,8 65,6
Аэфф 330,6 245,3 308,3 294,6 404,6 232,7
т, г 1000 948 919 913 884 884
Глина 2
40К 320,3 201,8 249,8 304,3 333,1 192,2
226Яа 13,7 8,9 6,4 10,8 20,3 9,6
232ТИ 16,4 16,2 24,1 13,8 22,5 14,1
Аэфф 63,8 48,2 60,3 56,1 79,5 45,2
т, г 1000 961 937 924 869 857
Глина 3
40К 972,4 583,4 729,3 904,3 1035,6 573,7
226Яа 35,2 22,2 17,2 28,5 51,0 25,3
232ТИ 61,1 57,4 87,4 55,6 86,1 54,9
Аэфф 202,1 149,3 196,5 182,2 256,2 148,5
т, г 1000 952 926 909 891 883
Продолжение таблицы № 1
Показатель Температура обжига
900 1050 1200 1350 1500
Глина 1
40К 1099,0 879,0 487,0 219,0 21,3
226Яа 33,8 31,7 20,1 18,8 9,3
232ТИ 125,7 94,3 60,6 36,0 4,2
Аэфф 296,1 233,3 142,6 85,4 16,6
т, г 795 716 644 579 521
Глина 2
40К 182,6 144,1 80,1 35,2 3,2
226Яа 7,5 7,1 4,5 4,2 2,1
232ТИ 27,2 20,3 13,1 7,7 0,8
Аэфф 59,4 46,5 28,8 17,4 3,4
т, г 801 721 652 586 499
Глина 3
40К 593,2 457,0 252,8 87,5 19,5
Яа 20,4 17,9 10,2 11,3 4,9
232ТЪ 97,1 78,7 47,7 91,0 6,7
Аэфф 200,1 160,8 94,9 137,5 15,4
т, г 786 709 647 568 514
Визуализация результатов исследований (см. рис. 1) приведенных в таблице №1 осуществлялась с помощью пакета ggplot2, который представляет собой систему для построения графиков на языке Я для статистических вычислений и графики [12].
Температура
- Глина 1 1-0
■ Глина 2
■ Глина 3 ^
■ Глина 1 ш
■ Глина 2 Ц
■ Глина 3 5
500 1000
Температура
Температура
■ Глина 1
■ Глина 2
■ Глина 3
■ Глина 1
■ Глина 2
■ Глина 3
Рис. 1. - Изменения удельных и эффективных удельных активностей ЕРН
при различных температурах обжига
Как видно из рис. 1 удельные активности и эффективная удельная активность ЕРН в процессе обжига изменяются различным образом.
Удельная активность 40K увеличивается в интервале температур от 20°С до 600°С. Повышение температуры обжига до 1500°С приводит к значительному уменьшению удельной активности. Однако такие высокие температуры обжига имеют место достаточно редко (например, при производстве портландцемента). В остальных случаях (керамический кирпич, керамзит) температуры обжига колеблются от 900°С до 1100°С, что обусловлено различным минералогическим и химическим составом глинистого сырья.
Несколько иначе ведет себя при обжиге 226Ra. В интервале температур 20 - 250°С его удельная активность уменьшается, далее наблюдается рост вплоть до 650°С. В дальнейшем при температурах более 650°С удельная
226т-)
активность Ra снижается по пологой кривой, достигая минимальных значений.
Удельная активность 232Th в отличие от первых двух радионуклидов возрастает вплоть до температур 900°С, после чего снижается вплоть до минимальных активностей.
Полученные экспериментальные данные позволяют получить зависимости удельных активностей и эффективной удельной активности ЕРН от температуры обжига.
Оценим связь и зависимость удельных активностей и эффективной удельной активности ЕРН от температуры обжига с помощью полиномиальной модели регрессии. Осуществим это на примере глины 1.
Полученные при исследовании данные (см. таблицу №1) имеют стохастический (вероятностный) характер, определим, имеет ли место связь между ними. Для этого воспользуемся функцией cor () входящей в базовый
пакет base языка R. Результаты корреляционного анализа представлены в таблице №2.
Таблица № 2
Величина коэффициента корреляции между исследованными показателями
Т, °С 40K, Бк/кг 226Ra, Бк/кг 232Th, Бк/кг Аэфф, Бк/кг m, г
Т, °С 1,000 -0,862 -0,611 -0,508 -0,780 -0,967
40K, Бк/кг -0,862 1,000 0,859 0,634 0,951 0,918
226Ra, Бк/кг -0,611 0,859 1,000 0,482 0,844 0,696
232Th, Бк/кг -0,508 0,634 0,482 1,000 0,828 0,598
Аэфф, Бк/кг -0,780 0,951 0,844 0,828 1,000 0,862
m, г -0,967 0,918 0,696 0,598 0,862 1,000
Представим полученные результаты с помощью матрицы диаграмм рассеяния значений зависимой и независимых переменных. Для ее создания воспользуемся функцией ера1гБ() из пакета gclus языка Я (см. рис. 2). Пары переменных, диаграммы которых расположены ближе к главной (диагональной) оси, скоррелированы сильнее.
Как видно из таблицы №2 и рис. 2 между исследованными показателями имеет место корреляция. Так между температурой обжига и остальными показателями она отрицательная, т.е. зависимость обратно пропорциональная. Причем между температурой и удельной активностью 40К, а также массой - зависимость весьма сильная. Между температурой и эффективной удельной активностью ЕРН - зависимость сильная. Умеренная зависимость наблюдается между температурой и удельной активностью
226^ 232г-р-|
Яа, а также Тп.
Весьма сильная зависимость просматривается между удельной
активностью 40К и удельной активностью 226Яа, эффективной удельной
активностью ЕРН и массой пробы. Зависимость между удельной
40т^ 232гг.,
активностью К и Тп носит умеренный характер.
Удельная активность 226Яа весьма сильно коррелирует с эффективной удельной активностью ЕРН, сильно с массой пробы и умеренно с удельной активностью 232Т^
Сильная зависимость имеет место между удельной активностью 232^ и эффективной удельной активностью ЕРН, а умеренная между удельной активностью 232^ и массой пробы. Эффективная удельная активность ЕРН имеет весьма сильную зависимость от массы пробы.
0 500 1БОО
1 I I I
0 БОО 1500
1 I I I I
20 40 60 00 I I I I
i i i i i i
О 40 ВО 120
i i i i i
600 ВОО 1000
I I I I
100 300
Рис. 2. - Матрица диаграмм рассеяния, упорядоченные и раскрашенные в соответствии с коэффициентами корреляции переменные: V! - температура,
°С; v2 - удельная активность 40К, Бк/кг; v3 - удельная активность 226Яа, Бк/кг; v4 - удельная активность 232ТП, Бк/кг; v5 - эффективная удельная активность ЕРН, Бк/кг; v6 - масса пробы, г Воспользуемся функцией 1т() из программы Я для нахождения модели полиномиальной регрессии, а для визуализации результатов пакетом ggp1ot2. Результат визуализации представлен на рис. 3.
Рис. 3. - Графическое представление связи и зависимости удельных активностей и эффективной удельной активности ЕРН для глины 1: 1 -
кривая экспериментально полученных данных; 2 - кривая уравнения
полиномиальной регрессии Уравнение полиномиальной регрессий для удельной активности 40К представлено формулой (2), при величине достоверной аппроксимации равной Я2 = 0,9463:
у = -1Е - 11х5 + 6Е - 08х4 - 8Е - 05х3 + 0,0495х2 - 10,743х + 2064 , (2)
для 226Яа формулой (3), величина достоверной аппроксимации составляет Я2 = 0,7622:
у = -3Е-16х6 + 7Е- 13х5 + 7Е- 10х4 -3Е- 06х3 + 0,0021х2 -0,5352х + 72,394 , (3) для 232^ формулой (4), при величине достоверной аппроксимации Я2 = 0,7527:
у = 7Е -16х6 - 3Е -12х5 + 3Е - 09х4 - 1Е - 06х3 - 0,0001х2 + 0,1528х + 69,72 , (4)
для Аэфф формула полиномиальной регрессии примет вид (5), а
2
величина достоверной аппроксимации составит Я = 0,8965:
у = 3Е - 15х6 - 1Е - 11х5 + 3Е - 08х4 - 3Е - 05х3 + 0,0118х2 - 2,0509х + 366,7 . (5) Из приведенных выше формул и величин достоверной аппроксимации видно, что расчетные параметры модели наиболее близки к фактическим данным для 40К и Аэфф, более низкие показатели 226Яа и 232Т^
Рассмотрим использование полученных зависимостей на следующем примере. Предположим, что на заводе по производству керамического кирпича обжиг сырца происходит при температуре 900°С. При этой температуре готовый керамический кирпич будет иметь эффективную удельную активность Аэфф - 296,1 Бк/кг. Используя полученные уравнения регрессии, рассчитаем, каковым будет этот показатель для температуры обжига - 950°С. В этом случае эффективная удельная активность ЕРН составит Аэфф - 272,85 Бк/кг. Следовательно с точки зрения радиационной безопасности увеличение температуры обжига на 50°С приводит к
уменьшению эффективной удельной активности ЕРН на 23,75 Бк/кг, что в свою очередь приведет к уменьшению величины облучения населения которое будет проживать в зданиях возведенных из этого материала [13].
Применение полученных регрессионных моделей на практике позволит рассчитывать удельные активности и эффективную удельную активность ЕРН при различных температурах обжига, а также в случаях, если будут наблюдаться изменения удельных активностей (например, при переходе на разработку следующего пласта сырья). Наряду с регрессионным возможно применение и других видов анализа, например карт Кохонена [14].
Таким образом, из изложенного выше следует, что увеличение температуры обжига глинистого сырья на определенном этапе приводит к уменьшению удельных и эффективной удельной активности ЕРН. Это способствует уменьшению дозы облучения населения. Однако при этом необходимо учитывать, что температуру обжига необходимо контролировать и по возможности уменьшать, снижая тем самым расходы на топливо. Следовательно, одним из перспективных путей направленных на уменьшение дозы облучения будет поиск оптимальной температуры обжига, при которой получение качественной керамической продукции будет сопровождаться наименьшими показателями эффективной удельной активности ЕРН.
Литература
1. Effects of Radiation on the Environment: Report to the General Assembly with Scientific Annex. United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR 2000). New York: UN, 2000. 842 p.
2. Соколов П.Э. Радиационные аспекты производства строительных материалов // Форум. 2016. №2 (8). С. 115-120.
3. Соколов П.Э., Дрик В.Ю. Естественная радиоактивность глинистых горных пород применяемых для производства строительных материалов //
Международный научно-исследовательский журнал. 2017. №12(66). Ч.5. С. 148-151.
4. Szabo Z., Völgyesi P, Nagy H., Szabo C., Kis Z., Csorba O. Radioactivity of natural and artificial building materials a comparative study. Journal of Environmental Radioactivity. 2013. №118. pp. 64-74.
5. Соколов П.Э., Чернышев П.Д. Влияние обжига на радиоактивность керамических материалов // Научное обозрение. 2017. №8. С. 34-40.
6. Kovacs T., Shahrokhi A., Sas Z., Vigh T., Somlai J. Radon exhalation study of manganese clay residue and usability in brick production. Journal of Environmental Radioactivity. 2017. №168. pp.15-20.
7. Sas Z., Szanto J., Kovacs J., Somlai J., Kovacs T. Influencing effect of heat-treatment on radon emanation and exhalation characteristic of red mud. Journal of Environmental Radioactivity. 2015. №148. pp. 27-32.
8. Li Y., Lu X., Zhang X. Determination of natural radioactivity, 222Rn and 220Rn exhalation rates and radiation hazards of fly ash and fly ash brick used in Baotou, China. Nuclear Technology and Radiation Protection. 2016. №31(3). pp. 282-290.
9. Стасов В.В., Зорина Л.В., Морозов А.Н. и др. Исследование связи радионуклидного состава почвенных фракций и атмосферных аэрозолей // Инженерный вестник Дона, 2007, №2. URL: ivdon.ru/ru/magazine/archive/n2y2007/29/.
10. Зорина Л.В., Стасов В.В., Бураева Е.А. Оценка техногенной составляющей загрязнения приземного слоя атмосферы свинцом-210 (на примере Юго-Востока г. Ростова-на-Дону) // Инженерный вестник Дона, 2007, №2. URL: ivdon.ru/ru/magazine/archive/n2y2007/18/.
11. Измерение активности гамма-излучающих радионуклидов на сцинтилляционном спектрометре с использованием пакетов программ SM и EXPRESS: Методические рекомендации / ВНИИФТРИ. М. 1993. 31с.
12. Кабаков Р.И. R в действии. Анализ и визуализация данных в программе R. Пер. с англ. М.: ДМК Пресс. 2014. 588с.
13. Ленников М.С., Соколов П.Э. Расчет эффективной эквивалентной дозы облучения населения с применением электронных таблиц // Международная научно-практическая конференция «Научные исследования высшей школы в области строительства и архитектуры». Новосибирск: НИЦ АЭТЕРНА, 2018. С. 87-92.
14. Соколов П.Э. Использование карт Кохонена для оценки радиоактивности строительного сырья // Вестник науки и образования Северо-Запада России. 2017. №1(3). С. 70-78.
References
1. Effects of Radiation on the Environment: Report to the General Assembly with Scientific Annex. United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR 2000). New York: UN, 2000. 842 p.
2. Sokolov P.E. Forum. 2016. №2 (8). pp. 115-120.
3. Sokolov P.E., Drik V.Y. Mezhdunarodnyj nauchno-issledovatel'skij zhurnal. 2017. №12(66). C.5. pp. 148-151.
4. Szabo Z., Völgyesi P, Nagy H., Szabo C., Kis Z., Csorba O. Journal of Environmental Radioactivity. 2013. №118. pp. 64-74.
5. Sokolov P.E., Tchernyshov P.D. Nauchnoe obozrenie. 2017. №8. pp. 34-40.
6. Kovacs T., Shahrokhi A., Sas Z., Vigh T., Somlai J. Journal of Environmental Radioactivity. 2017. №168. pp.15-20.
7. Sas Z., Szanto J., Kovacs J., Somlai J., Kovacs T. Journal of Environmental Radioactivity. 2015. №148. pp. 27-32.
8. Li Y., Lu X., Zhang X. Nuclear Technology and Radiation Protection. 2016. №31(3). pp. 282-290.
9. Stasov V.V., Zorina L.V., Morozov A.N. Inzenernyj vestnik Dona (Rus), 2007, №2. URL: ivdon.ru/ru/magazine/archive/n2y2007/29/.
10. Zorina L.V., Stasov V.V., Buraeva E.A. Inzenernyj vestnik Dona (Rus), 2007, №2. URL: ivdon.ru/ru/magazine/archive/n2y2007/18/.
11. Izmerenie aktivnosti gamma-izluchajushhih radionuklidov na scintil-ljacionnom spektrometre s ispol'zovaniem paketov programm SM i EXPRESS: Metodicheskie rekomendacii [Measurement of activity of gamma-emitting radionuclides in scintil-translational spectrometer using software packages such as SM and EXPRESS: guidelines]. VNIIFTRI. M., 1993. 31 p.
12. Kabakov R.I. R v dejstvii. Analiz I vizualizaciya dannyh v R. [R in action. Data analysis and graphics with R.] М.: DMK Press. 2014. 588p.
13. Lennikov M.C., Sokolov P.E. Mezhdunarodnaya nauchno-prakticheskaya konferenciya «Nauchnye issledovaniya vysshej shkoly v oblasti stroitel'stva i arhitektury». (International scientific and practical conference «Research of higher school in the field of construction and architecture») Novosibirsk: NIZ AETERNA, 2018. pp. 87-92.
14. Sokolov P.E. Ispol'zovanie kart Kohonena dlya ocenki radioaktivnosti stroitel'nogo syr'ya Vestnik nauki i obrazovaniya Severo-Zapada Rossii. 2017. №1(3). pp. 70-78.
