Basic ecological-hydrogeological examination of european part of Russia that suffered from Chernobyl nuclear power plant emergency

Full Text

ВВЕДЕНИЕ

26 апреля 2016 г. исполнилось тридцать лет с момента самой крупной антропогенной катастрофы ХХ в. на территории СССР $―$ аварии на Чернобыльской атомной электростанции (ЧАЭС), которая произошла в основном в результате так называемого “человеческого фактора” и привела к катастрофическим последствиям. В Европейской части России от этой аварии пострадали почти все центральные регионы, а наиболее сильно Брянская, Тульская и Калужская области. Проблемы с загрязнением окружающей среды радионуклидами будут сказываться еще долгое время. Известно, что активность радиоактивных выпадений исчезает полностью после 10 периодов полураспада радионуклидов.

Авария на атомной станции “Фукусима 1” в Японии произошла в марте 2011 г., предположительно в результате ошибки в проектировании станции (станция устояла при землетрясении, но не выдержала воздействия цунами). Выбросы в атмосферу, загрязнение почв и подземных вод как после аварии на ЧАЭС, так и после аварии на Фукусиме характеризуются практически аналогичными наборами радионуклидов, так как на обеих атомных станциях влиянию аварийной ситуации подверглось топливо. В настоящее время подземные воды являются уязвимым к загрязнению радионуклидами объектом окружающей среды.

В связи с этим рассмотрим базовые принципы оценки защищенности подземных вод от загрязнения, в исследуемом случае $—$ от радионуклидов на территории Калужской области, пострадавшей от аварии на ЧАЭС.

Базовая, основополагающая методика $—$ разработанная и усовершенствованная авторами статьи методика оценки защищенности грунтовых вод от загрязнения радионуклидами [1$–$3]. На основе этой методики проведены следующие исследования:

$•$ оценены защищенность и уязвимость грунтовых вод к загрязнению радионуклидами на территории Калужской области на различные расчетные сроки с шагом, равным периоду полураспада радионуклидов, это, по существу, $—$ прогнозная, предварительная оценка изменения экологической ситуации;

$•$ разработана методика оценки защищенности напорных подземных вод от загрязнения;

$•$ проведена унификация данной методики, что позволяет упростить и ускорить соответствующие оценки.

БАЗОВАЯ СРЕДНЕМАСШТАБНАЯ (1:200 000) ОЦЕНКА ЗАЩИЩЕННОСТИ И УЯЗВИМОСТИ ГРУНТОВЫХ ВОД ОТ ЗАГРЯЗНЕНИЯ РАДИОНУКЛИДАМИ

Авторами разработана оригинальная методика оценки защищенности и уязвимости грунтовых вод к загрязнению различными загрязняющими веществами (ЗВ), включая радионуклиды для мелко- и среднемасштабных эколого-гидрогеологических исследований [1$–$3]. Рассмотрим более подробно изменение степени защищенности и уязвимости грунтовых вод к загрязнению радионуклидами на территории Калужской области в масштабе 1:200 000.

Основные положения, которые следует учитывать при оценке и картографировании защищенности грунтовых вод от загрязнения различного масштаба, представлены в [1$–$3]. Дадим авторские определения оцениваемых показателей.

Защитная зона $—$ зона, отделяющая подземные воды от поверхностного загрязнения и имеющая двухуровневое строение: почвы и породы зоны аэрации.

Защищенность $—$ способность защитной зоны препятствовать проникновению загрязнения в подземные воды в течение определенного времени.

Природный защитный потенциал $—$ способность геологической среды (почв и пород зоны аэрации) удерживать загрязнение в защитной зоне, зависящая от литологических, фильтрационных и сорбционных свойств почв и пород.

Отношение реальной техногенной нагрузки изучаемой территории к естественной защищенности подземных вод называется их уязвимостью к загрязнению.

Основные радионуклиды при аварии на ЧАЭС в выпадениях $—$ ¹³⁷Cs и ⁹⁰Sr. За прошедшие после аварии 32 года активность этих радионуклидов в поверхностных выпадениях должна снизиться вдвое против начальной после аварии, так как период полураспада их составляет ~30 лет.

Для построения карт защищенности и уязвимости грунтовых вод к загрязнению необходимо иметь комплект карт, последовательность построения которых отражена в [1$–$3].

Карта защитной зоны строится по методике [1$–$3], она является базовой для построения карт защищенности и уязвимости грунтовых вод к любым ЗВ.

Карты защищенности грунтовых вод от загрязнения радионуклидами

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

Сорбция радионуклидов оказывает наибольшее влияние на задержку их породами и почвами. Основной параметр, характеризующий сорбцию, $―$ коэффициент распределения радионуклида [1$–$3]. Почвы обладают наибольшей удерживающей способностью по отношению к радионуклидам, что и определяет их как буфер более высокого порядка, чем зона аэрации.

Самые токсичные из долгоживущих радионуклидов $—$ ⁹⁰Sr и ¹³⁷Cs, поэтому оценка защищенности должна проводиться отдельно по каждому из них.

Защищенность грунтовых вод от любого ЗВ зависит от времени достижения фронтом загрязненных инфильтрационных вод водоносного горизонта (t_з) в пределах каждого выделенного на карте защитной зоны типового участка. Время прохождения растворенным в воде радионуклидом толщи почв и пород зоны аэрации мощностью М с заполнением их сорбционной емкости и последующим достижением уровня грунтовых вод можно определить так [2, 3]:

$t_{з}=\frac{M\vartheta п}{\text{υ}}+\frac{M\vartheta \text{δ}{K}_{р}}{W},\text{}$ (1)

где t_з $–$ время достижения загрязненными водами водоносного горизонта; K_p, л/кг, $—$ коэффициент распределения сорбции; δ, кг/дм³, $—$ объемная масса скелета грунта; п $—$ полная влагоемкость (в долях единицы); υ, м/сут, $—$ скорость просачивания инфильтрационного потока [4].

$\upsilon =\frac{1}{\text{θ}}\sqrt[3]{W^2{k}_{\text{ф}}}$ (2)

где θ $—$ естественная влажность пород (в долях единицы), W $—$ инфильтрационное питание (м/ сут); k_ф $—$ коэффициент фильтрации (м/сут).

Первое слагаемое в (1) характеризует движение влаги в ненасыщенной зоне (или движение нейтрального ЗВ), второе $―$ физико-химическое взаимодействие (сорбцию) в системе порода$–$вода (или задержку ЗВ почвой или породой). В случае, когда коэффициент распределения значительно превышает единицу (как это характерно для радионуклидов), первым слагаемым в (1) можно пренебречь; вторым можно пренебречь, когда коэффициент распределения значительно меньше единицы.

Шкалу категорий естественной защищенности грунтовых вод от загрязнения целесообразно строить в зависимости от Т $—$ периода полураспада радионуклида. В этом случае выделяются следующие категории:

незащищенные грунтовые воды: t_з < Т; t_з < 30 лет;

слабозащищенные грунтовые воды: Т < t_з < 2Т; 30 лет < t_з < 60 лет;

среднезащищенные грунтовые воды: 2Т < t_з < 3Т; 60 лет < t_з < 100 лет;

условно защищенные грунтовые воды: 3Т < t_з < 10Т; 100 лет < t_з < 300 лет;

защищенные грунтовые воды: t_з >10Т; t_з >300 лет.

Выделение категории защищенности по времени продвижения ЗВ через защитную зону, по сути, $—$ приближенная прогнозная оценка загрязнения грунтовых вод, в данном случае радионуклидами.

Рассмотренный подход к построению карт естественной защищенности грунтовых вод от радиоактивного загрязнения может быть использован при составлении аналогичных карт по оценке загрязнения другими высокотоксичными ЗВ.

Карта уязвимости грунтовых вод к загрязнению ¹³⁷Cs

Оценка уязвимости грунтовых вод проводится по реальному загрязнению и осуществляется путем совмещения карты защищенности грунтовых вод от загрязнения ¹³⁷Cs и карты радиационной обстановки после аварии на ЧАЭС [8].

Для характеристики уязвимости грунтовых вод к загрязнению ¹³⁷Cs разработана категоризация степени уязвимости по отношению техногенной нагрузки по ¹³⁷Cs к защищенности грунтовых вод от загрязнения ¹³⁷Cs (табл. 1).

Степень уязвимости грунтовых вод к загрязнению ¹³⁷Cs в зоне радиоактивного следа сразу после аварии на ЧАЭС выражена на рис. 1 так: очень сильно уязвимые грунтовые воды сосредоточены на отдельных участках в среднем течении р. Рессеты, сильно уязвимые $―$ в долинах рек в центральной части следа, среднеуязвимые $―$ по долинам рек, притоков, оврагов и на востоке у р. Вытебети, слабо уязвимые $―$ по водоразделам рек, очень слабо уязвимые $—$ по периферии следа, неуязвимые сосредоточены на высоких террасах рек и частично водоразделах рек.

Для оценки уязвимости грунтовых вод спустя тридцать лет после аварии активность поверхностных выпадений ¹³⁷Cs уменьшена вдвое по сравнению с зафиксированной на карте [8], согласно истечению одного периода полураспада радионуклидов (табл. 2). На карте уязвимости грунтовых вод (рис. 2) на этот срок показано, что в зоне следа очень сильно уязвимых грунтовых вод уже нет, сильноуязвимые сохранились узкими полосами в долинах рек и вдоль русел, среднеуязвимые $—$ на низких террасах рек, значительно уменьшившись по площади; площади слабоуязвимых грунтовых вод увеличились, также увеличились площади очень слабо уязвимых грунтовых вод, площади условно уязвимых и неуязвимых не изменились. Кроме этого, построен комплект карт на сроки истечения остальных периодов полураспада.

СРЕДНЕМАСШТАБНАЯ МЕТОДИКА ОЦЕНКИ ЗАЩИЩЕННОСТИ ПОДЗЕМНЫХ, НАПОРНЫХ ВОД ОТ ЗАГРЯЗНЕНИЯ

Рассмотрим подходы к оценке защищенности и уязвимости подземных напорных вод к загрязнению, что также важно для оценки экологической ситуации с питьевыми подземными водами, которые, как правило, сосредоточены в напорных водоносных горизонтах, залегающих ниже горизонтов грунтовых вод. Исследования по данной проблеме в России начаты В.М. Гольдбергом [5$–$7]. В работе [3] использованы оценки с помощью специальных индексов проникновения.

Таблица 1. Степень уязвимости грунтовых вод к загрязнению ¹³⁷Cs (определяется суммой весовых значений защищенности (0–4) и техногенной нагрузки (0–6), которые указаны в скобках для соответствующих показателей): после аварии на ЧАЭС (2–10 — суммарные весовые значения cтепени уязвимости): неуязвимые — t₃ > 300 лет, условно уязвимые <1, очень слабо уязвимые — 2÷1, слабо уязвимые — 4÷3, среднеуязвимые ― 6÷5, сильноуязвимые ― 8÷7, очень сильно уязвимые ―10÷9)

Оценка защищенности напорных подземных вод от загрязнения имеет особенности по сравнению с аналогичной выше приведенной оценкой для грунтовых вод. Основное отличие содержится в определении защитной зоны напорных вод.

Рис. 1. Картосхема уязвимости грунтовых вод к загрязнению 137Cs после аварии на ЧАЭС на территории Калужской области.

Защитная зона напорных вод имеет также два уровня:

первый $—$ собственно грунтовые воды (или другие напорные водоносные горизонты, залегающие выше изучаемого напорного горизонта), которые имеют двойное воздействие на напорные воды: они могут быть источником загрязнения напорного водоносного горизонта, если они разгружаются (перетекают) в нижезалегающий напорный водоносный горизонт; или они являются областью защиты напорных вод, если в них разгружаются (перетекают) нижезалегающие напорные воды;

второй $—$ водоупорная толща (водоупор), отделяющий горизонт напорных вод от горизонта грунтовых вод.

Защитная зона $—$ это зона, отделяющая напорные подземные воды от загрязнения, поступающего из вышезалегающего водоносного горизонта, и имеющая двухуровневое строение: вышезалегающий водоносный горизонт и водоупор, разделяющий эти водоносные горизонты.

Таблица 2. Степень уязвимости грунтовых вод к загрязнению ¹³⁷Cs (определяется суммой весовых значений защищенности (0–4) и техногенной нагрузки (0–4), которые указаны в скобках для соответствующих показателей) спустя 30 лет после аварии на ЧАЭС (2−8 — суммарные весовые значения уязвимости)

Рис. 2. Картосхема уязвимости грунтовых вод к загрязнению 137Cs спустя 30 лет после аварии на ЧАЭС на территории Калужской области.

Защищенность $—$ способность защитной зоны препятствовать проникновению загрязнения в подземные напорные воды из вышезалегающего водоносного горизонта в течение определенного времени.

Природный защитный потенциал $—$ способность геологической среды удерживать загрязнение в защитной зоне, зависящая от литологических, фильтрационных, гидродинамических и сорбционных свойств пород.

Оценка и картирование защищенности напорных подземных вод, залегающих ниже горизонта грунтовых вод, проводилась для части территории Калужской области в зоне радиоактивного следа от аварии на ЧАЭС (масштаб карт 1:200 000).

Гидрогеологические условия на юге Калужской области: безнапорные водоносные горизонты (грунтовые воды) включают в себя воды четвертичных (аллювиальных, гляциальных, флювиогляциальных, болотных, пролювиальных горизонтов), меловых и юрских отложений, все эти водоносные горизонты связаны между собой и не имеют выдержанных водоупоров внутри этой толщи (комплекса); напорные водоносные горизонты пресных подземных вод включают в себя множество связанных между собой водоносных горизонтов каменноугольного возраста; между двумя этими толщами (водоносными комплексами) залегает выдержанный по площади верхнеюрский водоупор, разделяющий эти толщи.

Для построения карты защитной зоны необходимо исследовать первый ее уровень с точки зрения оценки степени влияния грунтового водоносного комплекса на возможность загрязнения напорного водоносного комплекса через него. Степень этого влияния оценивается величиной водопроводимости (km), или способностью профильтровать через себя загрязненные воды, и разностью уровней первого и второго водоносных комплексов (H₁ − H₂), обеспечивающей интенсивность и направленность перетекания подземных вод между этими комплексами.

Второй уровень защитной зоны характеризуется параметрами водоупора: коэффициентом фильтрации (k_ф) и мощностью (m₀). В связи с тем, что нет подробных сведений о литологическом строении водоупора (раздельного слоя), коэффициент фильтрации задавался постоянным для всей его толщи, значения его изменялись только для различных расчетных сценариев.

Таким образом, для оценки потенциала защитной зоны использовались три показателя: km, k_ф и m₀ (табл. 3).

Таблица 3. Определение потенциала защитной зоны (сумма трех показателей: весовых значений водопроводимости (km), мощности водуопора m₀ и разницы уровней водоносных горизонтов (H₁ – H₂); для km вес присвоен каждому отдельному интервалу реальных размерных значений (м²/сут) — в квадратных скобках (1–6); для m₀ вес присвоен интервалам размерных значений (м) — в квадратных скобках (1, 3, 6); для (H₁ – H₂) вес присвоен интервалам размерных значений (м) — в квадратных скобках (1–6))

Оценка проводилась следующим образом:

$–$ всем показателям присваивались весовые значения в порядке увеличения в зависимости от степени ухудшения защитных свойств (увеличение значения km приводит к ухудшению защитных свойств, т.е. загрязненная вода интенсивнее профильтруется внутри грунтового потока и быстрее достигнет водоупора, поэтому весовое значение по мере ухудшения защитных свойств увеличивается; аналогичная ситуация и с разностью уровней (H₁ − H₂); весовое значение m₀, наоборот, увеличивается при уменьшении мощности водоупора);

$–$ затем проводилось суммирование весовых значений трех показателей и устанавливался потенциал защитной зоны: высокий (весовое значение 0$–$5), средний (5$–$10), слабый (10$–$15), очень слабый (15$–$18), очень высокий $―$ эта категория вводится дополнительно для области, где грунтовые воды не являются источником загрязнения напорных вод, где их уровень залегает ниже пьезометрического уровня напорных вод.

При наличии сведений о литологическом строении водоупора и его фильтрационных свойствах последний показатель может быть включен при оценке защитного потенциала. Желательно, чтобы диапазоны изменения весовых величин отдельных показателей были близкими, чтобы не искажать суммарные значения всех показателей. Результаты оценки потенциала защитной зоны приведены в табл. 3.

На схематической карте защитной зоны (рис. 3) выделены 42 типовых участка с различным строением защитной зоны по трем показателям и по величине защитного потенциала (табл. 3), описание которых приведено в экспликации к карте (табл. 4).

Таблица 4. Экспликация к схематической карте защитной зоны напорных подземных вод (фрагмент)

(пример описания типового участка № 29 (3-4-6-13 слабый защитный потенциал): первая цифра в скобке — мощность водоупора, разделяющего грунтовые воды от напорных вод, ― составляет 0–10 м; вторая — водопроводимость напорного водоносного горизонта — 50–100 м²/сут; третья — превышение уровня грунтовых вод над пьезометрическим уровнем напорного водоносного горизонта, равное нулю; четвертая — сумма весовых значений показателей составляет 13, что соответствует слабому защитному потенциалу защитной зоны напорного водоносного горизонта (табл. 3); жирный шрифт — участки, на которых пьезометрический уровень напорного водоносного горизонта превышает уровень грунтового водоносного горизонта, что фиксирует перетекание напорных подземных вод из напорного водоносного горизонта в грунтовый и определяет защищенность напорного горизонта от загрязнения через грунтовый водоносный горизонт)

Потенциал защитной зоны

На карте (рис. 3) видно, что в долинах рек Болвы, Рессеты и Жиздры распространены территории с очень слабым защитным потенциалом, где напорные и грунтовые воды не разделены водоупором, а представляют собой единый водоносный комплекс. На незначительных территориях пьезометрический уровень напорного водоносного комплекса превышает уровень грунтового водоносного комплекса (хотя и незначительно), что определяет отсутствие перетекания загрязненных грунтовых вод в напорные. На большей части территории отмечается слабый защитный потенциал защитной зоны напорного водоносного комплекса, и только на незначительных участках распространен средний и высокий защитный потенциал.

Карта защищенности подземных вод от загрязнения строится путем наложения карты защитной зоны и карты интенсивности проникновения нейтрального, несорбируемого мигранта через защитную зону.

Защищенность напорных вод от любого ЗВ зависит от t_з. Время прохождения растворенным в воде несорбируемым ЗВ толщи почв и пород зоны аэрации мощностью М вычислялось с использованием модели поршневого вытеснения по первому члену уравнения (1).

Для определения времени t_з по (1) необходимо рассчитать скорость просачивания по формуле (2), в которой один из главных параметров $—$ инфильтрационное питание (W).

В ислледуемом случае W $—$ не прямой параметр, характеризующий поступление атмосферных осадков в напорный водоносный комплекс. Питание поступает в этот комплекс через перетекание из грунтового водоносного комплекса. Это и определяет некоторые трудности определения этой величины.

Попытки определения величины питания напорного горизонта через известный параметр перетекания [10, 11] привели к получению величин, во много раз превышающих значение инфильтрационного питания грунтового водоносного горизонта, что связано с тем, что параметр перетекания оценивается в режиме водоотбора из водоносного горизонта с привлечением значительных объемов воды к водозаборным скважинам в радиусе их влияния.

При использовании уравнения (3), приведенного в работе [9] для определения величины взаимосвязи W_вз при перетекании, также получается его завышенное значение, хотя и гораздо меньшее, чем в случае использования метода [10, 11] по аналогичной причине:

$W_{\text{вз}}=k_0\frac{H_1-H_2}{m_0}$ , (3)

где k₀, м/сут, $—$ коэффициент фильтрации водоупора; m₀, м, $—$ мощность водоупора; Н₁ $–$ Н₂, м, $—$ разность уровней напорного и грунтового водоносных горизонтов.

Для решения этой проблемы авторами выбран другой способ определения питания $―$ перетекания из грунтового водоносного горизонта в напорный. Основой данного способа явились режимные гидрогеологические наблюдения за уровнями грунтового и напорного водоносного горизонтов на соседней от изучаемой территории. По данным этих наблюдений, максимальная величина инфильтрационного питания верхнего грунтового комплекса составила 4.59 м/сут, а питание нижнего напорного каменноугольного комплекса $—$ 1.14 м/ сут. Сравнивая эти две величины, можно сделать вывод, что инфильтрационное питание напорного водоносного комплекса составляет 25% от питания грунтового водоносного комплекса.

Таким образом, для определения инфильтрационного питания напорного водоносного горизонта установлены два значения: 25% (по данным режимных наблюдений) и 10% (минимальное значение) от величины инфильтрационного питания грунтового комплекса.

За исходное инфильтрационное питание грунтового водоносного комплекса приняты величины 200, 100 и 60 мм/год, которые установлены для изучаемого региона при оценке защищенности грунтовых вод от загрязнения [1$–$3].

Значения k₀ на изучаемой территории не определялись, поэтому расчеты проводились для трех значений: 10^-4, 10^-3, 10^-2 м/сут (0.0365, 0.365, 3.65 м/год).

Мощность водоупора (m₀) определялась по карте защитной зоны (рис. 3) в интервалах 5, 10, 20, 30 м.

Вместо θп (естественная влажность пород) принято значение активной пористости, равное 0.2 (в долях единицы).

Таким образом, рассмотрены три варианта сценария оценки защищенности напорных подземных вод от загрязнения нейтральными ЗВ для разных величин инфильтрационного питания грунтовых вод: 200, 100 и 60 мм/год; в каждом варианте рассмотрены по два значения инфильтрационного питания напорных вод: 25 и 10%, три значения коэффициента фильтрации водозабора: 10^-4, 10^-3, 10^-2 м/сут. Результаты расчета времени проникновения нейтрального ЗВ (t_з) в подземные воды при W = 100 мм/год приведены в табл. 5. Соответственно, построены три схематические карты защищенности напорных подземных вод от загрязнения нейтральными ЗВ.

По результатам расчетов установлены следующие степени защищенности напорных подземных вод от загрязнения нейтральными ЗВ: t_з = 0$–$10 лет $—$ незащищенные, t_з = 10$–$30 лет $—$ слабозащищенные, t_з = 30$–$50 лет $—$ среднезащищенные, t_з > 50 лет $—$ защищенные.

Рис. 3. Картосхема защитной зоны напорных подземных вод.

Таблица 5. Расчет времени проникновения нейтрального ЗВ (t_з) в подземные воды (год) при исходном инфильтрационном питании грунтовых вод W = 100 мм/год

Рис. 4. Картосхема защищенности напорных подземных вод от загрязнения нейтральными ЗВ.

На рис. 4 (при инфильтрационном питании грунтовых вод W = 100 мм/год $—$ среднее значение) показано, что на низких террасах рек Рессеты и Жиздры и на северо-западе долины р. Болвы простирается единый водоносный комплекс, что характеризует незащищенность напорных вод от загрязнения, поступающего из грунтового потока, а на незначительных участках уровни напорного комплекса превышают уровни грунтового, что характеризует его защищенность от загрязнения сверху (рис. 3). При 25%-й инфильтрации на более высоких террасах рек Рессеты и Жиздры и на северо-западе долины Болвы, а также на водораздельном пространстве р. Рессеты на востоке изучаемой территории распространены незащищенные ПВ, уязвимые к загрязнению нейтральными ЗВ; на западе территории распространены слабозащищенные ПВ. При 10%-й инфильтрации на востоке территории распространены слабозащищенные ПВ, а на западе защищенные.

При анализе трех сценариев показано, что по мере увеличения инфильтрационного питания степень защищенности подземных вод уменьшается, увеличивается их уязвимость к загрязнению нейтральными ЗВ. К числу нейтральных ЗВ могут быть отнесены макрокомпоненты (хлориды, сульфаты), нитраты, нефтепродукты и др.

УНИФИЦИРОВАННЫЙ МЕТОД ОЦЕНКИ ЗАЩИЩЕННОСТИ ПОДЗЕМНЫХ ВОД ОТ ЗАГРЯЗНЕНИЯ

С использованием многолетнего опыта оценки защищенности подземных вод от загрязнения различными ЗВ, включая радионуклиды, проведено усовершенствование методики [1$–$3], которое позволило использовать ее для решения комплексных экологических задач при загрязнении подземных вод различными ЗВ и при изменениях климатических условий или антропогенных факторов.

Унификация методики оценки защищенности подземных вод сводится к упрощению расчетной процедуры по выражениям (1) и (2), базирующейся на выделении некоторого математического ядра, которое в расчетах не меняется, и некоторой изменяющейся величины, с которой проводятся вычисления в зависимости от решаемой задачи.

Рассмотрим вычислительные процедуры для трех решаемых задач.

Первая задача $―$ миграция различных загрязняющих веществ

Как отмечено выше, защищенность грунтовых вод от любого ЗВ зависит от t_з в пределах каждого выделенного на карте защитной зоны типового участка, которое определяется по выражениям (1) и (2).

Объединим (1) и (2):

$t_{з}=\frac{M\vartheta п}{\sqrt[3]{W^2{k}_{\text{ф}}}}+\frac{M\vartheta \text{δ}{K}_p}{W},$ (3)

перепишем (3):

$t_{з}=a+bx,$ (4)

где: $a_1=\frac{M\vartheta п}{\sqrt[3]{W^2{k}_{ф}}},$$b_1=\frac{M\vartheta \text{δ}}{W}, x={К}_{р}.$

Выражение (4) $—$ рекуррентная формула для расчета времени миграции ЗВ с различными значениями коэффициента распределения при всех других неизменяющихся величинах: величины a и b характеризуют природные особенности защитной зоны.

Вторая задача $—$ изменения инфильтрационного питания

Преобразуем (1):

$t_{з}=\frac{M\vartheta п}{\sqrt[3]{W^2}\sqrt[3]{k_{\text{ф}}}}+\frac{M\vartheta \text{δ}{K}_p}{W},$ (5)

$t_{з}=\frac{a_2}{\sqrt[3]{y^2}}+\frac{b_2}{y},$ (6)

где $a_2=\frac{M\vartheta п}{\sqrt[3]{k_{ф}}},$ $b_2=M\vartheta \text{δ}{K}_{р},$ y = W.

Выражение (6) $—$ рекуррентная формула для расчета времени миграции ЗВ при изменении инфильтрационного питания при неизменных других величинах: величина a характеризует природные особенности защитной зоны, а величина b характеризует взаимодействие $―$ сорбцию ЗВ природной средой $―$ породой.

Третья задача $—$ миграция различных загрязняющих веществ при различных величинах инфильтрационного питания

$t_{з}=\frac{a_2}{\sqrt[3]{y^2}}+\frac{cx}{y},$ (7)

где $a_2=\frac{M\vartheta п}{\sqrt[3]{k_{ф}}}, c=M\vartheta \text{δ},$ x = K_р, y = W.

Выражение (7) $—$ рекуррентная формула для расчета времени: величины a и c характеризуют природную среду в условиях двух меняющихся показателей $—$ ЗВ и инфильтрационного питания.

Предложенная усовершенствованная методика среднемасштабной оценки защищенности подземных вод от загрязнения позволяет проследить возможные вариации развития процесса загрязнения в зависимости от индивидуального ЗВ $—$ от нейтрального до высокотоксичного при различных сценариях изменения климата: от иссушения до его увлажнения.

ВЫВОДЫ

С использованием базовой методики среднемасштабной оценки защищенности и уязвимости грунтовых вод к загрязнению радионуклидами получены предварительные прогнозные оценки изменения эколого-гидрогеологической ситуации (на части территории Калужской области, пострадавшей от аварии на ЧАЭС) на различные прогнозные сроки, которые могут быть использованы для планирования научных исследований более крупного масштаба, выбора мест для расположения водозаборов грунтовых вод и ведения экологического мониторинга.

Разработана оригинальная новая методика оценки защищенности подземных, напорных вод от загрязнения, которая имеет свои особенности и отличия по сравнению с аналогичной оценкой для грунтовых вод.

Для оценки защитного потенциала защитной зоны использованы три показателя: величина водоотдачи грунтового водоносного горизонта, разность уровней первого и второго водоносных горизонтов $—$ первый уровень защитной зоны; коэффициент фильтрации водоупора и его мощность $—$ второй уровень защитной зоны. Построена схематическая карта защитной зоны напорных подземных вод.

Оценка защищенности подземных напорных вод от загрязнения нейтральными ЗВ проводилась путем определения времени их миграции по модели поршневого вытеснения. Рассмотрены три варианта-сценария оценки защищенности напорных подземных вод от загрязнения нейтральными ЗВ для разных значений инфильтрационного питания грунтовых вод (200, 100 и 60 мм/год), в каждом варианте рассмотрены по два значения инфильтрационного питания напорных вод (25 и 10%) и по три значения коэффициента фильтрации водоупора (10⁻⁴, 10⁻³, 10⁻² м/сут). Построены три схематические карты защищенности напорных подземных вод от загрязнения нейтральными ЗВ.

Методика может быть использована для выбора мест расположения водозаборов пресных напорных вод и ведения экологического мониторинга.

Проведено усовершенствование ранее разработанной методики оценки защищенности подземных вод от загрязнения различными ЗВ, включая радионуклиды. Данная методика может быть использована для решения комплексных экологических задач при загрязнении подземных вод различными ЗВ и в условиях изменяющегося климата или действия антропогенных факторов. Усовершенствование методики сводится к унификации оценки защищенности подземных вод, к упрощению расчетной процедуры, базирующейся на выделении некоторого ядра, которое в расчетах не изменяется, и некоторой изменяющейся величины, с которой проводятся вычисления в зависимости от решаемой задачи. Рассмотрены три типа задач: миграция различных ЗВ, изменение климата, а также изменение климатических условий и набора ЗВ. Унифицированный метод оценки защищенности подземных вод от загрязнения позволяет уменьшить затраты времени на проведение расчетных процедур, что в свою очередь позволяет заранее построить и иметь комплекты карт защищенности подземных вод от загрязнения любыми ЗВ, при любых природных и антропогенных условиях и на любые расчетные прогнозные сроки.

Работа выполнена по государственному заданию ИГМ СО РАН. Финансирующая организация: Министерство науки и высшего образования Российской Федерации.

About the authors

A. P. Belousova

Author for correspondence.
Email: anabel@iwp.ru
Russian Federation, Moscow

E. E. Rudenko

Email: anabel@iwp.ru
Russian Federation, Moscow

References

Supplementary files

Supplementary Files

Action

1. Fig. 1. Map of the vulnerability of groundwater to 137Cs pollution after the Chernobyl accident in the Kaluga region.

Download (296KB)

Indexing metadata

2. Fig. 2. Map of groundwater vulnerability to 137Cs pollution 30 years after the Chernobyl accident in the Kaluga Region.

Download (297KB)

Indexing metadata

3. Fig. 3. Map of the protective zone of pressure groundwater.

Download (259KB)

Indexing metadata

4. Fig. 4. Map of the protection of pressure groundwater from pollution by neutral pollutants.  

Download (300KB)

Indexing metadata