Диффузионно-адвективный перенос радона в многослойных геологических средах

Abstract

Представлены решения стационарного диффузионно-адвективного уравнения переноса радона в многослойных геологических средах численными методами. Для построения однородных консервативных разностных схем использовали интегро-интерполяционный метод (метод баланса). Полученные явную и неявную конечноразностные схемы решали методом правой и левой прогонки. Произведен выбор оптимального шага сетки из сопоставления данных численных расчетов и аналитического решения, полученного в пакете Mathematica 5 для стационарного случая. С помощью полученной численной модели диффузионно-адвективного переноса радона можно рассчитывать распределения поровой активности радона и плотности потока радона по глубине и во внешнюю атмосферу. Приведены результаты численных расчетов функции поровой активности радона по глубине для случая трехслойной среды, представленной разными типами грунтов и пород с сильно различающимися физико-геологическими характеристиками.