Естественнонаучные проблемы Арктического региона : восьмая региональная научная студенческая конференция, Мурманск, 15-16 мая 2007г. : труды конференции. Мурманск, 2008.

ВЛИЯНИЕ ВХОДНЫХ ПАРАМЕТРОВ ДИФФУЗИОННОГО БЛОКА МИГРАЦИОННОЙ ПОДМОДЕЛИ РАДИОГЕОЭКОЛОГИЧЕСКОЙ МОДЕЛИ ПОДЗЕМНОГО ХРАНИЛИЩА ЯДЕРНЫХ МАТЕРИАЛОВ НА ЗАГРЯЗНЕНИЕ ВМЕЩАЮЩЕГО МАССИВА А.Н. Попов1, П.В. Амосов1,2 1Кольский филиал Петрозаводского государственного университета 2Горный институт КНЦ РАН, vosoma@goi.kolasc.net.ru В работе рассмотрено влияние входных параметров диффузионного блока ми­ грационной подмодели, как составной части радиогеоэкологической модели подземно­ го хранилища РАО и ядерных материалов (Мельников, 2001), на потенциальное загряз­ нение вмещающего массива. В качестве анализируемого радионуклида выбран долго- 79 живущий радиологически опасный изотоп Se , для которого справедливо приближе­ ние модели предела элементной растворимости. В сферической системе координат дифференциальное уравнение, описывающее диффузионный перенос растворенной радиоактивности C (r ,t ) с нулевыми начальными и граничными условиями на бесконечности, а также использованием модели предела элементной растворимости в области источника активности могут быть записаны сле­ дующим образом: дС^ дС = А , 1 д -Л С , dt г З 'К С(г,0) = 0 для r > R 0; С(со,/) = 0 , C(i?0,/) = CS, где R0 - радиус источника; Cs - предел растворимости конкретного элемента; Я - по­ стоянная распада радиоактивного изотопа; Da - наблюдаемый коэффициент диффузии (Мельников, 2001). Параметр Da даже в простейшем приближении - линейной адсорб­ ции - является функцией пористости е , коэффициента распределения K d и плотности «скелета» породы p s - Da = sDr D„ £ + ( \ - e ) p sK d R . В последнем выражении D - ко­ эффициент диффузии радионуклида по обводненному норовому пространству породы, который определяется посредством выражения Dp = Dv ■G f , где Dv - молекулярный коэффициент диффузии ионов в свободной воде; Gf - геометрический фактор, опреде- g ляемый через G f = —у , где 8 и и г - сжимаемость и извилистость пор, соответственно; Т Rp - коэффициент задержки в матрице породы. В рамках указанной модели пространственно-временное распределение концен­ трации C (r,t) радионуклида (учтены механизмы диффузии, сорбции-десорбции и ра­ диоактивного распада) может быть описано выражением С(г, 0 = ^ 2 г ехр ~ w r - R°] erfc ( r -R p ) 2Ѵ д 7 - у/At + 62