CC BY
56
УДК 539.16.08 + 51-73
ОЦЕНКА ОБЛАСТИ ВЛИЯНИЯ ЗАГРЯЗНЕННОГО УЧАСТКА ПОЧВЫ ПРИ РЕШЕНИИ ЗАДАЧ РАДИАЦИОННОГО МОНИТОРИНГА
МЕТОДОМ «IN SITU»
12 2 1 Жуковский А.И. , Кутень С.А. , Хрущинский А.А. , Толкачев А.Н. ,
Гузов В.Д.1, Кожемякин В.А.1, Чудаков В.А.3.
1УП «АТОМТЕХ», г. Минск, Республика Беларусь 2НИИ ядерных проблем БГУ, г. Минск, Республика Беларусь 3Международный государственный экологический университет им. А.Д. Сахарова, г. Минск, Республика Беларусь e-mail: [email protected]
Рассчитана интенсивность нерассеянных гамма-квантов, поглощенных в цилиндрическом детекторе, помещенном над равномерно загрязненной почвой. Установлена зависимость радиуса области влияния для сцинтилляционного детектора в зависимости от глубины залегания радионуклидов цезия на примере почв региона Тохоку, Япония.
Ключевые слова: детектор, объемный радиоактивный источник, гамма-квант, функция отклика, критический радиус.
Введение
Для решения задач, связанных с радиационным мониторингом объектов окружающей среды, в частности почв, подвергшихся радиоактивному загрязнению, на практике используют различные инструментальные методы и средства контроля.
Стандартный метод экспресс-оценки радиоактивного загрязнения почвы беспробоотборным способом - измерение функции отклика детектора (например, сцинтиляционного) на некоторой высоте над загрязненным участком почвы. Функция отклика в виде энергетического распределения импульсов гамма-излучения зависит от характеристик детектора, почвы и геометрии измерения. На этапе калибровки устройства детектирования требуются калибровочные энергетические распределения импульсов гамма-излучения в геометрии измерения, получение которых слишком трудоемко и дорого. Для этих целей, а также для оптимизации задачи эффективно используют Монте-Карло моделирование процесса переноса гамма-излучения. При использовании метода «in situ» почва не может быть непосредственно использована в качестве математической модели, так как представляет собой объект измерения условно бесконечных размеров. Поэтому
очень важно определить область влияния, т.е. тот минимальный участок почвы, который вносит основной вклад в функцию отклика детектора. После этого задача Монте-Карло моделирования решается корректно.
Данная работа отражает результаты теоретических исследований в рамках решения задач по разработке и созданию средства измерений для контроля радиоактивного загрязнения почв региона Тохоку в Японии.
Целью настоящей работы являлась разработка алгоритма определения полуаналитическим методом параметров области влияния загрязненного участка почвы на функцию отклика детектора и получение численных данных для их использования в моделировании методом Монте-Карло процесса переноса гамма-излучения в геометрии измерения «in situ».
Аналитическая связь «источник-детектор» для определения влияния нерассеянных гамма-квантов на показания детектора над загрязненной почвой
В бесконечно однородной среде, ослабляющей одинаково во всех направлениях, интенсивность излучения от объемного цилиндрического
источника является функцией линейных размеров и физико-химических параметров источника и детектора, а также расстояния между ними [1]. Как правило, аналитическая связь между параметрами и интенсивностью нерассеянного излучения устанавливается в какой-либо точке пространства, что эквивалентно приближению точечного детектора. Важно получить такую связь для детектора конечных размеров. Для оценки радиуса области влияния достаточно рассмотреть только нерассеянное излучение, что в свою очередь значительно упростит выражение аналитической связи «источник - детектор».
Рассмотрим цилиндрический детектор (цилиндр высотой На и радиусом г^), расположенный на высоте И над загрязненной почвой (рисунок 1), которая представлена слоем цилиндрической формы (высота И,, радиус г,). Мы можем использовать такую простую модель реального устройства детектирования только для определения влияния нерассеянных гамма-квантов на чувствительность детектора и определения области влияния на детектор.
Интересующая функция отклика детектора будет определяться нерассеянным потоком гамма-квантов до и после поглощения в детекторе. Показания детектора (функция отклика) конечных размеров 1ЛеЬ полученные в результате поглощения детектором нерассеянного гамма-излучения, пропорциональны следующей величине:
(В) = { сК | суа^)ехр (-ц(Е К (Е )
т/ т/, АкК
Ус1
:(1 - ехр (-ц (ЕК ))
(1)
где Я - расстояние между произвольно выбранной точкой г1 на поверхности почвы и произвольно выбранной точкой г2 на поверхности детектора, представлющее собой сумму Я,, Яа и Ял - длин прямолинейных участков траекторий движения гамма-кванта в поишЛ итттухе и детекторе, соответственно; а (7[} - удельная активность нуклида; цДЕ), ца(Ё и ц^Ё - линейные коэффициенты ослабления гамма-излучения соответственно для почвы, воздуха и материала детектора для энергии Е; членами ц и ца представлен поток нерассеянных гамма-квантов, проходящий расстояние Я, в почве, а также расстояние Яа в воздухе соответственно. Последний член уравнения (1) представляет гамма-кванты, которые были выпущены в направ-
К
лении детектора = —, достигли его поверх-
К
ности и были поглощены в теле детектора в точке г2. Интегрирование выполняется для объема почвы dVs = ^Зг1 для объема детектора
dVd = d3Г2.
Последний множитель в подинтегральном выражении (1) описывает вероятность поглощения гамма-квантов в материале детектора. В общем случае он не равен единице и отражает возможность для гамма-квантов пролететь детектор не поглотившись в нем. Эти процессы могут быть существенными для детекторов малых размеров и для приповерхностных траекторий гамма-квантов.
б
Рисунок 1 - Детектор и источник в цилиндрической системе координат (а) с проекцией на плоскость ху (б)
Для упрощения задачи примем, что радионуклид с энергией Е равномерно распре-
1
X
делен по объему почвы, т.е. a(r{) = const, и подставим эту величину в уравнение (1).
Интегрирование в уравнении (1) будем проводить в цилиндрической системе координат (переменные zb рь ф1 и z2, р2, ф2 для почвы и детектора, соответственно).
Величины Rs, Ra и Rd определяются из геометрических соображений через расстояние R между точкой испускания гамма-кванта в почве и точкой поглощения его в детекторе:
R = n - n;
R = ЛR +(-2 - -i )2 ;
(2)
Rí =P¡ +P¡ - 2PlP2COs(9l -ф2)•
u - h R --R;
-2 - -l
-, - u
R = --- R;
u
= -lP 2 + - 2 P¡ -(-l + - 2 )Pl
Rí
P 2 X
X cos(
(ф, -ф¡)-(-¡ --l yÍD);
2
D = (rdRí)2 -(PlP¡sin(9l -ф2))2.
(6)
В случае пересечения гамма-квантом нижней поверхности детектора по траектории Я выражения для Яа и Ял имеют более простой вид:
Путь гамма-кванта в почве Rs определяется из соотношения:
R = ^ R.
(3)
R = R;
—2 - -l
R = -2- h - h + Kl ¡ R
(7)
-2 - -l
где И, - толщина загрязненного слоя почвы; И - высота расположения центра детектора над почвой.
Уравнения для траекторий гамма-кванта в воздухе и в материале детектора зависят от направления Пк и различаются в зависимости от пересечения гамма-квантом торцевой или боковой поверхности детектора. Эти расстояния можно найти путем решения системы из двух уравнений. Первое из них описывает траекторию фотона, проходящую через две точки - г и г2:
• - xl _ y - у,
- - —
Х2 X1 y2 У1 —2 -1
(4)
Второе уравнение является уравнением для торцевой и боковой поверхности детектора, соответственно, и имеет вид:
— = h + h-h^ ' 2
2 2 2 • + У = r¡ •
(5)
Окончательные выражения для путей Яа и Ял при пересечении гамма-квантом боковой поверхности детектора (по траектории Я) можно представить в виде:
Необходимо отметить, что все величины в выражении (1) зависят от азимутов ф] и ф2 только через их разность ф] - ф2. Поэтому, вводя новые угловые переменные ф = ф! - ф2 и и = ф] + ф2 и проводя интегрирование по и, выражение (1) можно привести к виду:
h' h' + h+hd /¡ r rd ¡x
e, (E)= J d— 1 J d— ¡ jP,dPi jP¡dP¡ J(2n-T)3(-„ — ¡, P„P¡, ф)dф;
G(—i, —¡, Pi, P¡, ф)= exp '(E R¡-* a (eR ) (i-e-, (-ц d (E)Rd )) nR
(8)
Выражение (8) позволяет рассчитать вклад «прямых» гамма-квантов в функцию отклика детектора после вычисления пятикратного интеграла.
Расчет параметров области влияния для загрязненного радионуклидами 134Cs и 137Cs плодородного слоя почвы региона Тохоку, Япония
Линейные коэффициенты ослабления для выражения (8) в почве, воздухе и материале детектора были рассчитаны при помощи программы XCOM (Версия 3.1, июнь 1999, M.J. Berger, S.M. Seltzer и др.) на основании широко известных массовых коэффициентов ослабления M.J. Berger и J.H. Hubbell [2].
- 2 -l
—2 -—1
Массовые доли элементов, формирующих плодородный слой почвы региона Тохоку, представлены в таблице 1. Элементный состав получен на основе данных представленных в работе [3]. В качестве органической составляющей вместо гумусовых кислот использовались данные элементного состава торфа. Рассчитанный массовый коэффициент ослабления гамма-излучения для почвы региона Тохоку в зависимости от энергии гамма-кванта показан на рисунке 2.
Приведенные выше массовые коэффициенты ослабления использовались в качестве основы для построения функции цДЕ), используемой в (8) как известная аналитическая функция. Аналогичным образом были построены функции для ца(£) и ^d(E).
Представленные ниже расчеты выполнялись для сцинтилляционного детектора NaI(Tl) высотой hd = 63 мм и радиусом rd = 31,5 мм. При этом геометрический центр детектора располагался над загрязненной почвой на высоте h = 53,5 мм. Типичное значение плотности почвы региона Тохоку составляет 1,3 г/см3.
Элементный состав
Пятикратный интеграл (8) рассчитан в пакете «МаЛвтаНеа» с использованием метода интегрирования «AdaptiveQuasiMonteCaгlo» [4]. На рисунке 3 показана типичная зависимость показаний детектора от радиуса участка почвы с равномерно распределенным радионуклидом (глубина до 10 см).
С ростом радиуса участка почвы показание детектора стремится к насыщению (истинному значению), поэтому зависимость показаний детектора от радиуса участка почвы может быть аппроксимирована кривой:
1 ^ (г ) = а + Ь{1- ехр (- сг )), (9)
где коэффициенты а, Ъ, с определяются в процессе аппроксимации.
Истинное показание детектора, которое выходит на асимптотику при росте радиуса, с учетом выражения (9) можно представить в виде:
= И = а + Ь. (10)
Таблица 1
зы региона Тохоку
Элемент H C N O Na Mg Al
Массовая доля 0,004 0,035 0,002 0,622 0,011 0,006 0,051
Элемент Si P S K Ca Fe
Массовая доля 0,243 0,001 0,001 0,012 0,003 0,010
Рисунок 2 - Массовый коэффициент ослабления для почвы региона Тохоку в зависимости от энергии гамма-кванта
Рисунок 3 - Зависимость показаний детектора от радиуса участка почвы с равномерным распределением 137Сб по глубине до 10 см
Из рисунка 3 видно, что показания детектора не изменяются для радиусов, превышающих некоторый критический радиус (эффективный радиус области влияния). Влияние на показание детектора бесконечного слоя почвы, равномерно загрязненной в слое заданной толщины, эквивалентно влиянию слоя такой же толщины с конечным радиусом, равным критическому.
Критический радиус фактически определяет область источника, влияющую на показание прибора.
Для определения критического радиуса целесообразно использовать допустимое относительное отклонение Р показаний детектора ^Ж?) для заданного радиуса участка почвы г относительно показаний детектора ^(да). По определению, величина Р записывается в виде:
использованием вместо уравнения (8) следующего выражения [5] для интенсив-ности излучения в центре рассматриваемого детектора при заданном распределенном источнике (приближение точечного детектора):
4. (Е) = к (Е >]- е2[Ц (Е > (Е )
+
Ц (Е К ]-
гх + И(ца (Е )-Ц, (Е))
(13)
Ц, (Е)(*+*,) Е ["
| 1 х йх},
-(Е> х^Ц, (Е)г )2 + х2 ]
где £'п(х) - интегральная показательная функция [5]:
Р = 1 (да)-1 * (г)
1 ае.
(да)
(11)
да -хх
Еп ( х ) = | —йХ.
1
Хп
(14)
В соответствии с (11) показания детектора для радиуса участка почвы г будут отклоняться в пределах погрешности Р от своего максимального (истинного) значения, соответствующего значению г = да.
Согласно уравнениям (9-11) для заданной погрешности Р значение критического радиуса участка почвы можно представить в виде:
1п
Ь
Р(а + Ь)
(12)
Зависимости показаний детектора от радиуса участка почвы с равномерным распределением по глубине и их аппроксимации кривой насыщения (9) на глубине 2 см, 5 см, 10 см и 15 см приведены на рисунке 4.
Соответствующие значения критических радиусов почв, рассчитанные с использованием выражений (8) и (12), представлены в таблице 2 для различных глубин залегания радионуклидов 137С8 и 134С8.
Следует отметить, что величина относительного отклонения в не может быть меньше типичной ошибки измерения (в нашем случае порядка 10 %).
Для этого значения относительного отклонения в = 10 % в таблице 2 в последнем столбце приведены также значения критических радиусов почвы, рассчитанные согласно (12) с
Как и следовало ожидать, влияние размера детектора на эффективный радиус области влияния максимально (порядка 9 %) для загрязнения с глубиной залегания радионуклидов 2 см и практически исчезает уже для загрязнения с глубиной 5 см.
Для Монте-Карло моделирования процесса измерения удельной или поверхностной активности с использованием сцинтилляционного детектора следует использовать критический радиус области влияния из предпоследнего столбца таблицы 2, для величины относительного отклонения в порядка 10 %.
Рисунок 4 - Зависимости показаний детектора от радиуса участка почвы с равномерным распределением 137Сб по глубине и их аппроксимации кривой насыщения (9) на глубине 2 см, 5 см, 10 см и 15 см
Гс =
с
Таблица 2
Критический радиус участка почвы региона Тохоку для различных глубин залегания радионуклидов 137С8 и 134С8
Относительное отклонение показаний детектора от истинного значения, в %
3 % 5 % 10 % 10 % 1
Глубина залегания, см Критический радиус почвы, см
2 115 98 76 69
5 83 71 55 52
10 71 60 47 46
15 64 55 42 44
1Значения критического радиуса для точечного детектора, рассчитанные с использованием уравнений (12)и(13)
Заключение
Аналитический метод, использованный для описания прохождения нерассеянного гамма-излучения, позволил получить и оптимизировать математическую модель равномерно загрязненной по глубине почвы, необходимую для последующего Монте-Карло моделирования. В частности, определена область влияния загрязненной почвы региона Тохоку, в зависимости от глубины залегания радионуклидов Cs и Cs.
Полученная аналитическая связь «источник - детектор» и рассчитанные на ее основе параметры, в частности, значения критического радиуса источника излучения (загрязненного участка почвы) для различных глубин залегания контролируемых радионуклидов использованы в ходе разработки и создания портативного гамма-спектрометра МКС-АТ6101ДР Проведенные натурные испытания МКС-АТ6101ДР и последующая эксплуатация спектрометра в префектурах Японии, подвергшихся радиоактивному загрязнению, подтвердили высокую степень достоверности представленных теоретических исследований.
Список использованных источников
1. Гусев, Н.Г. Защита от ионизирующих излучений. Том 1 : Физические основы защиты от излучений : учебное пособие / Н.Г. Гусев, В.А. Климанов, В.П. Машкович, А.П. Суворов. - 3-е изд. - М. : Энергоатомиздат, 1989. - 512 с.
2. Berger, M.J. XCOM: Photon Cross Section Database (v. 3.1) / M.J. Berger [et. al]. http://physics.nist.gov/PhysRefData/Xcom/Text/do wnload.html (accessed 17.09.2013)
3. International Commission on Radiation Units and Measurements, Gamma-ray spectrometry in the environment. ICRU Report 53, ICRU Publication Bethesda, Maryland, 1994.
4. Фихтенгольц, Г.М. Курс дифференциального и интегрального исчисления / Г.М. Фихтенгольц. - Том 3. - М. : Физматлит, 2001. - 662 с.
5. Гусев, Н.Г. Защита от ионизирующих излучений. Том 2: Защита от излучений ядерно-технических установок : учебное пособие / Н.Г. Гусев, Е.Е. Ковалев, В.П. Машкович, А.П. Суворов. - 3-е изд. - М. : Энергоатомиздат, 1990. - 352 с.
6. Прудников, А.П. Интегралы и ряды. Том 1: Элементарные функции / А.П. Прудников, Ю.А. Брычков, О. И. Маричев. - 2-е изд. -М. : Физматлит, 2003. - 631 с.
IMPACT RANGE ESTIMATION OF POLLUTED SOIL AREA FOR RADIATION MONITORING
BY «IN SITU» METHOD
1 2 2 11 1 3
Zhukouski A. , Kutsen S. , Khrutchinsky A. , Tolkachev A. , Guzov V. , Kojemiakin V. , Chudakov V.
'SPE "ATOMTEX", Minsk, Belarus 2Research Institute for Nuclear Problems, Minsk, Belarus international Sakharov Environmental University, Minsk, Belarus e-mail: [email protected]
Abstract. Intensity of «direct» gamma-quanta absorbed in cylindrical detector has been determined for the detector placed over soil uniformly contaminated in depth. The dependence of impact range radius for NaI(Tl) detector to cesium depth contamination is established for soil from Tohoku region, Japan.
Keywords: detector, volume radioactive source, gamma quantum, response function, critical radius.
References
1. Gusev N.G., Climanov V.A., Mashkovich V.P., Suvorov A.P. Zashchita ot ioniziruyushchikh izlucheniy. Tom 1: Fizicheskiye osnovy zashchity ot izlucheniy [Ionizing radiation's protection. V. 1: Physical basis of radiation protection], 3rd edition, Moscow, Energoatomizdat Publ., 1989, 512 p. (in Russian).
2. Berger M.J., Hubbell J.H., Seltzer S.M., Chang J., Coursey J.S., Sukumar R., Zucker D.S., Olsen K., XCOM: Photon Cross Section Database (v. 3.1). http://physics.nist.gov/PhysRefData/Xcom/Text/download.html (accessed 17.09.2013).
3. International Commission on Radiation Units and Measurements, Gamma-ray spectrometry in the environment. ICRU Report 53, ICRU Publication Bethesda, Maryland, 1994.
4. Fikhtengolts G.M. Kurs differentsialnogo i integralnogo ischisleniya. Tom 3 [Course of differential and integral calculus. V. 3], Moscow, Fizmatlit Publ., 2001, 662 p. (in Russian).
5. Gusev N.G., Kovalev E.E., Mashkovich V.P., Suvorov A.P. Zashchita ot ioniziruyushchikh izlucheniy. Tom 2: Zashchita ot izlucheniyyaderno-tekhnicheskikh ustanovok [Ionizing radiation's protection. V. 2: Radiation protection against nuclear installations] 3rd edition, Moscow, Energoatomizdat Publ., 1990. 352 p. (in Russian).
6. Prudnikov A.P., Brychkov U.A., Marichev O.I. Integraly i ryady. Tom 1: Elementarnyye funktsii. [Integrals and series. V. 1: Elementary Functions.] 2nd edition, Moscow, Fizmatlit Publ., 2003, 631 p. (in Russian).
Поступила в редакцию 14.03.2014.
