ФИЗИКА
УДК 550.835.2
ПОЛЕВОЙ МЕТОД ИЗМЕРЕНИЯ КОЭФФИЦИЕНТА ДИФФУЗИИ РАДОНА
И ТОРОНА В ГРУНТЕ
В.С. Яковлева
Национальный исследовательский Томский политехнический университет,
634050, г.Томск, пр. Ленина, 30
E-mail: [email protected]
Разработан простой и достоверный полевой метод измерения эффективного коэффициента диффузии радона и торона в грунте и других пористых материалах. Анализ результатов численных исследований переноса радона и торона в поверхностных слоях грунта показал, что плотность потока торона с поверхности грунта практически не зависит от скорости адвекции почвенных газов и изменяется только с изменением коэффициента диффузии. Этот результат показал преимущества использования торона по сравнению с радоном в предлагаемом методе. Приведено сравнение нового метода с существующими, разработанными ранее. Метод может быть полезен при решении задач массопереноса радона в пористых средах и газообмена между почвой и атмосферой.
Ключевые слова: радон, торон, коэффициент диффузии, грунт, полевой метод, плотность потока.
(с) Яковлева В.С., 2014
PHYSICS
MSC 76R50
IN-SITU MEASURING METHOD OF RADON AND THORON DIFFUSION
COEFFICIENT IN SOIL
V.S. Yakovleva
National Research Tomsk Polytechnic University, 634050, Tomsk c., Lenin st., 30,
Russia
E-mail: [email protected]
A simple and valid in-situ measurement method of effective diffusion coefficient of radon and thoron in soil and other porous materials was designed. The analysis of numerical investigation of radon and thoron transport in upper layers of soil revealed that thoron flux density from the earth surface does not depend on soil gas advective velocity and varies only with diffusion coefficient changes. This result showed the advantages of thoron using versus radon using in the suggested method. The comparison of the new method with existing ones previously developed. The method could be helpful for solving of problems of radon mass-transport in porous media and gaseous exchange between soil and atmosphere.
Key words: radon, thoron, diffusion coefficient, soil, in-situ method, flux density
(c) Yakovleva V.S., 2014
Введение
Актуальность радоновой тематики в различных областях науки и практики до сих пор продолжает расти [1]-[3]. Это связано с замечательными индикаторными свойствами радиоактивного газа радона, а также его изотопа торона.
Результаты моделирования процессов переноса радона и торона в грунте имеют большое значение не только в области радиоэкологии, в основном для построения карт радоноопасности территорий, но и в строительстве, геологоразведке и геофизике, например, при изучении литосферно-атмосферных связей, или поиске предвестников землетрясений [2]-[6]. При решении подобного рода задач коэффициент диффузии радона в грунте необходимо знать непосредственно для той точки, или на той территории, которая представляет интерес.
Коэффициент диффузии радона в пористой среде
Согласно исследованиям, произведенным в работе [7], коэффициент диффузии определяется содержанием воды в порах грунта и его пористостью
Бе = п ■ Б0 ■ е-6тп -6т14п (1)
где £о - коэффициент молекулярной диффузии радона (торона) в воздухе (00=0,11 см2/с); т - коэффициент водонасыщения, т.е. степень влажности или заполнения объема пор водой, п - пористость грунта, отн. ед.
Выражение (1) было получено на основе анализа результатов измерений коэффициента диффузии в 1073 образцах грунта разного типа (от песчано-гравийных до глинистых) с различной влажностью (от полностью водонасыщенных до сухих), пористостью и плотностью.
Известно, что коэффициент диффузии радона/торона для рыхлых поверхностных отложениях изменяется незначительно, со средним значением 0,03 см2/с. Содержание влаги в поверхностных грунтах при нормальных условиях находится, в основном, в пределах от 10 до 30%, что соответствует диапазону изменения коэффициента диффузии от 0,02 до 0,04 см2/с. При полном заполнении пор водой коэффициент диффузии снижается до значения 10-5 см2/с, а в сухом грунте - достигает 0,05 см2/с.
Однако, если учитывать геолого-географические различия исследуемых территорий, это выражение не всегда может давать точный результат. Например, в работе [8] получено существенное расхождение экспериментально определенных значений коэффициента диффузии с рассчитанными из выражения (1) при т >19%. Поэтому, для получения более надежных результатов желательно параллельно использовать полевые методы измерения коэффициента диффузии.
Метод измерения коэффициента диффузии
К настоящему времени известно не так много полевых методов измерения эффективного коэффициента (молекулярной) диффузии Бе радона в грунте [9]-[10], и они имеют ряд существенных недостатков. Например, метод [9] заключается во введении мощного мгновенного источника радона на определенную глубину с последующим измерением изменения во времени активности радона в точке детектирования. Затем
фиксируют момент времени, в который измерена максимальная активность радона и рассчитывают эффективный коэффициент диффузии радона. При этом имеются следующие недостатки:
1) необходимость использования достаточно мощного источника радона;
2) использование специального устройства для мгновенного введения радона на заданную глубину грунта;
3) большая продолжительность одного измерения (не менее суток);
4) большая погрешность определения коэффициента диффузии, связанная с конечными размерами используемых в опыте источника радона и детектора.
Другой метод [10], основанный на диффузионной модели переноса радона в пористых средах, заключается в том, что одновременно измеряют объемную активность (ОА) радона в двух точках, расположенных на расстоянии от 0,2 м до 1 м друг от друга. Измерение в первой точке производят на глубине ^1 от 0,2 м до 0,5 м, а во второй - на удвоенной глубине. Этот метод обладает существенным ограничением
- он применим только при скорости адвективного переноса радона и=0 см/с, поскольку используется упрощенная диффузионная модель переноса и игнорируется адвективный процесс переноса, хотя в реальных условиях и может достигать 4-10-4 см/с и более [11]. При наличии адвективных потоков в грунте со скоростью, например, и = 4-10-4 см/с оценки _Ое получаются заниженными в ~80 раз по сравнению с истинным значением, что указывает на один из основных недостатков метода [10]
- огромная погрешность оценок эффективного коэффициента диффузии радона.
Задачей настоящей работы являлась разработка простого и достоверного полевого метода измерения эффективного коэффициента диффузии радона и торона в грунте. Поскольку физико-химические свойства радиоактивных газов радона и то-рона практически не отличаются [12], коэффициенты диффузии радона и торона являются при одних и тех же физико-геологических условиях одинаковыми. Из-за различий в ядерно-физических характеристиках радона и торона, а именно, в периодах полураспада, торон обладает большим преимуществом для использования его в предлагаемом методе, по сравнению с радоном. Согласно исследованиям [13] было получено, что плотность потока торона (ППТ) с поверхности грунта практически не зависит от скорости адвекции и почвенных газов (не более 2,6% в широком диапазоне значений 0<и <10-3 см/с) и изменяется только с изменением коэффициента диффузии.
Таким образом, для определения _Ое можно использовать диффузионную модель переноса изотопов радона в пористых средах и рассчитывать эффективный коэффициент диффузии радона и торона в грунте из простого аналитического соотношения
где дтп - ППТ с поверхности грунта, Бк м-2 с-1; Атп - ОА торона в почвенном воздухе, Бк м-3; _Ое - эффективный (объемный) коэффициент диффузии радона и торона, м2 с-1; п - пористость грунта, отн. ед.; Xтп - постоянная распада торона, с-1.
В предлагаемом методе (2) существует ограничение на минимальную глубину измерения: она должна быть не менее 15 см. Это ограничение обусловлено тем, что ОА торона в почвенном воздухе на глубине ~10 см достигает своего равновесного значения, при котором торон находится в радиоактивном равновесии с материнским
(2)
радионуклидом 224Ra, образующимся в природной цепочке 232Th, и далее не изменяется с глубиной и не зависит от и. При измерениях на меньших глубинах увеличивается погрешность оценки эффективного коэффициента диффузии. Таким образом, разработанный метод является:
1) простым, постольку не требует использования мощного источника радона и специального устройства для мгновенного ведения радона на заданную глубину грунта;
2) достоверным, поскольку не содержит погрешностей, связанных с влиянием скорости адвективного переноса почвенных газов на результат измерения.
Метод был апробирован на экспериментальной площадке ТПУ-ИМКЭС [3]. Плотность потока торона определяли способом измерения плотности потока радона и торона с поверхности грунта по а-излучению, описанным в работе [14].
Объемную активность торона в почвенном воздухе на глубине 15 см измеряли с помощью радиометра радона и торона RTM 2200 (SARAD GmbH, Германия) в соответствии с инструкцией по измерению радона и торона в грунте. На площадке также произвели отбор пробы поверхностного грунта для последующего измерения пористости грунта в лабораторных условиях в соответствии со стандартными методиками. Значение эффективного коэффициента диффузии радона и торона в грунте De составило, в среднем, 3-10 6 м2 с \ что хорошо согласуется с литературными данными для подобных типов поверхностных грунтов.
Заключение
Разработан простой и достоверный полевой метод измерения эффективного коэффициента диффузии радона и торона в грунте и других пористых материалах, который может быть полезен при решении задач массопереноса радона в пористых средах, газообмена между почвой и атмосферой, и пр. При небольшой доработке его можно модифицировать в метод неразрушающего контроля для измерения коэффициента диффузии радона и торона в (монолитных) строительных конструкциях.
Библиографический список
1. Паровик Р.И., Макаров Е.О., Фирстов П.П. Математическое моделирование фрактальной размерности геосреды и сейсмическая активность Южной Камчатки // Вестник КРАУНЦ. Физикоматематические науки. 2011. №2(3). С. 42-49.
2. Паровик Р.И. Фирстов П.П. Фазовый анализ временных рядов геофизических полей // Вестник КРАУНЦ. Физ.-мат. науки. 2013. №1(6). С. 23-29.
3. Яковлева В.С., Каратаев В.Д., Вуколов А.В., Ипполитов И.И., Кабанов М.В., Нагорский П.М., Смирнов С.В., Фирстов П.П., Паровик Р.И. Скоординированный многофакторный эксперимент по анализу процессов поступления почвенного радона в приземный слой атмосферы // АНРИ. 2009. № 4. C. 55-60.
4. Методы измерения плотности потока радона и торона с поверхности пористых материалов: монография / В.С. Яковлева; Томский политехнический университет. Томск: Изд-во ТПУ, 2011. 174 с.
5. Firstov P.P., Yakovleva V.S., Shirokov V.A., Rulenko O.P., Filippov Yu.A., Malysheva O.P. The nexus of soil radon and hydrogen dynamics and seismicity of the northern flank of the Kuril-Kamchatka
subduction zone // Annals of Geophysics. 2007. V. 50. № 4. Р. 547-556.
6. Паровик Р.И., Фирстов П.П. Апробация новой методики расчета плотности потока радона с поверхности (на примере Петропавловск-Камчатского геодинамического полигона) // АНРИ. 2009. № 3. С. 52-57.
7. Rogers V.C., Nielson K.K. Multiphase radon generation and transport in porous materials // Health Phys. 1991. V. 60. P. 807-815.
8. Hosoda M., Shimo M., Sugino M., Furakawa M., Fukushi M. Effect of Soil Moisture Content on Radon and Thoron Exhalation // J. Nucl. Sci. Technol. 2007. V. 44. № 4. P. 664-672.
9. Булашевич Ю.П., Карташов Н.П. Определение коэффициента диффузии радона в горных породах
методом мгновенного источника // Изв. АН СССР. Сер. физика земли. 1967, №10, С.71-76.
10. Патент РФ №2332687. G01T1/178. опубл. 27.08.2008.
11. Яковлева В.С., Рыжакова Н.К. Оценка скорости конвекции радона в грунтах по измеренным значениям поровой активности // Геоэкология. Инженерная геология. Гидрогеология. Геокриология. 2003. №5. С. 466-469.
12. Сердюкова А.С., Капитанов Ю.Т. Изотопы радона и короткоживущие продукты их распада в
природе. М.: Атомиздат, 1979. 294 c.
13. Яковлева В.С. Моделирование влияния состояния атмосферы и литосферы на динамику плотности потока радона и торона с поверхности земли // Известия ТПУ. 2010. Т. 317. №2. С.162-166.
14. Яковлева В.С. Методы и приборы контроля полей а-,в-,Y излучений и радона в системе "грунт-атмосфера"[Электронный ресурс]: автореферат диссертации на соискание ученой степени доктора технических наук; спец. 05.11.13 / В. С. Яковлева; НИ ТПУ; науч. конс. П. М. Нагорский. Томск, 2013. URL: http://www.lib.tpu.ru/fulltext/a/2013/71.pdf.
Поступила в редакцию / Original article submitted: 20.05.2014