Экопрогнозы техногенных процессов
При изучении и оценке свойств геологической среды, формирования и развития техногенных процессов, необходимых для эколого-гидрогеологического обоснования мер защиты биосферы и, прежде всего, человека, выполняются экопрогнозы этих процессов различными методами.[7,9,10] Методы экопрогнозов важнейших техногенных эколого-гидрогеологических процессов по степени их изученности, делятся на три группы. К первой относят наиболее изученные техногенные процессы, для которых научные основы прогнозов в настоящее время хорошо разработаны. Это процессы вторичной консолидации осушенных рыхлых пород, депрессионного уплотнения песчано-глинистых пород при снятии пластового давления в напорных водоносных горизонтах, формирования депрессионной воронки при осушении горных разработок и на водозаборных сооружениях, вторичного засоления почв на орошаемых землях, режима уровня грунтовых вод на орошаемых массивах, подтопления территорий и сооружений и др. С этими процессами связано решение ряда прикладных гидрогеодинамических и гидрогеомеханических задач. Вторая группа охватывает менее изученные техногенные процессы, научные основы прогнозирования которых разработаны недостаточно. К ней относят группы процессов, определяющих влияние техногенеза на гидрогеохимические параметры водоносных горизонтов, подвергающихся техногенному воздействию. Это техногенное загрязнение пресных подземных вод при их эксплуатации, миграция различных загрязнений в горных породах, различных загрязнений в пластовых условиях и др. Третья группа процессов — качественно прогнозируемые — мало изучены и крайне сложные, где ведущая роль принадлежит взаимодействию разнонаправленных по тепломассообмену техногенных процессов и др. К ним относится прогнозирование процессов суффозионно-карстового выщелачивания карбонатных пород и процессов взаимодействия подземных вод и горных пород с микроорганизмами и природными газами и др. Для прогнозов первой и второй группы техногенных процессов в настоящее время используются известные гидрогеодинамические и гидрогеохимические методы (в том числе для сложных условий с применением моделирования и ЭВМ). Методика прогнозных решений геоэкологических задач с использованием названных методов подробно изложена в опубликованных работах [ 9, 10] В практике эколого-гидрогеологических исследований чаще всего приходится сталкиваться с прогнозом наиболее сложных техногенных гидрогеохимических и радиогенных процессов, негативное влияние которых проявляется в форме химического, санитарно-бактериологического загрязнения и зараженности подземных вод. Пример решения таких наиболее сложных геоэкологических задач на конкретном объекте рассмотрен в 5 разделе настоящого пособия. Как показывает многолетний опыт, при эксплуатации водозаборных сооружений питьевого назначения или системы дренажа производственного назначения прогноз техногенного изменения качества воды сводится к решению следующих задач: 1) оценка возможности проникновения загрязнения в "область захвата" водозабора; 2) определение времени подтягивания загрязненных вод к участку водозабора; 3) оценка степени изменения качества подземных вод продуктивного горизонта, происходящего под влиянием техногенного загрязнения; 4) разработке рекомендаций по защите геологической среды от негативного влияния техногенных гидрогеохимических процессов, при возможном загрязнении подземных вод. Решение этих задач для простых гидрогеологических условий можно выполнить приближенными методами расчетов; в сложных гидрогеологических условиях работы водозаборного сооружения (сложность граничных условий потока в плане, фильтрационная неоднородность продуктивной толщи и др.) целесообразно использовать метод математического моделирования. Важнейшим фактором, определяющим условия миграции загрязнений от техногенных очагов до водозабора является гидродинамический, характеризующий наличие или отсутствие на изучаемом участке естественного движения подземных вод. В связи с этим целесообразно рассматривать методы прогноза техногенного загрязнения подземных вод на водозаборных сооружениях раздельно для условий естественного потока и бассейна, где естественная скорость потока очень мала и при расчетах ею можно пренебречь. Расстояние, на котором могут находиться от водозабора существующие объекты возможного загрязнения, приближенно можно рассчитать по уравнению Tq l= -------, mn где T - время продвижения загрязненного потока; m - мощность водоносного горизонта; n - активная пористость пород; q- удельный расход потока; A - коэффициент, учитывающий сорбцию, и коэффициент распределения собрируемого компонента между водой и природой: A=(l+b)/b В случае отсутствия данных по величине коэффициента b значение l определяется приближению по уравнению l= Tq/ Amn Согласно рекомендации А.Е. Орадовской и Н.Н. Лапшина (1987) при проектировании линейного берегового (инфильтрационного) водозабора, расположенного от реки на расстоянии х0 и имеющего длину z, расходы потоков, поступающих со стороны реки Qp и от береговой Qб зоны, могут определяться по уравнениям : zkm(Hр – Hв) zkm(Hк – Hв) Qp =--------------- Qб= ------------------- , x0 xk-x0 где Hр и Hв – уровни воды в водозаборе и реке; Hк - естественный уровень подземных вод на берегу на расстоянии хк от водозабора. Те же авторы предлагают определять концентрацию любогокомпонента в водах водозабора после смешения речных и береговых вод по уравнению Qp Ср + QбСб Сi= ------------------ , Qp+ Qб где Ср, Сб – концентрации компонента в реке и береговых водах. Время подтягивания загрязненной воды к водозабору, расположенному в удалении от реки на расстоянии чф от источника загрязнения, определяется по формуле πmn (чф 2-ч0 2) Т=--------------------- , Qв
где ч0 – радиус водозабора; Qв - его расход. В случае перетока загрязненной воды из смежного водоносного горизонта максимальная концентрация загрязняющего компонента Сi в водозаборе напорных вод составляет
Сп (Сэ-Сп) (km)э Сi=---------------------- , (km)э+(km)п где Сэ и Сп - концентрации загрязняющего компонента в водоносных породах загрязненного и в водах эксплуатируемого горизонтов; (km)э и (km)п – водопроводимости эксплуатируемого и питающего загрязненного горизонтов. Для водозаборного участка, расположенного в естественном потоке, решение задачи, производится на основании учета положения в плане потока так называемой нейтральной линии тока относительно выявленного контура техногенного загрязнения подземных вод. Если очаг загрязнения располагается за пределами "области захвата" водозабора, то нет оснований ожидать техногенного загрязнения подземных вод, если в контуре "захвата", то решение таких задач возможна аналитическими и другими методами, изложенными в работах.[ ] Время достижения загрязненными водами с постоянным расходом уровня грунтовых вод определяется по уравнению: t=n*H/k *[m/H-ln(1+m/H)] , где n – пористость пород зоны аэрации, Н – высота слоя сточных вод в хвостохранилище, k – средний коэффициент фильтрации, m – мощность пород зоны аэрации. Время попадания в горные выработки вод с повышенной минерализацией вычисляется по формулам: - при отсутствии естественного потока, т.е. в бассейне для «большого колодца» T=πmnx12/Q, - для траншеи T=2πmnlnx/Q, - в неограниченном в плане водоносном горизонте при наличии естественного потока по его направлению T=n/Ki *[x1-xAln(x1/xA+1)]; - в полуограниченном пласте с контуром питания T=πnmd/3Q* [2+(x1/d)3-3x1/d]; Если граница вод повышенной минерализации совпадает с кнтуром питания водоносного горизонта, то T=2πmnd2/3Q, где Т – время подтягивания минерализованных вод; n- активная пористость; d – расстояние от выработок до контура с постоянным напором; x1 – расстояние от контура питания до границы с минерализованными водами: xA – расстояние от центра выработок до контура минерализованных вод; i – уклон естественного потока. Минерализация дренажных вод на любой момент времени после смещения вод для условий неограниченного бассейна определяется для «большого колодца» по зависимости C=C0+(C1+C0)/π*arccosT/t, где С,С1,С0 – соответственно минерализация дренажных вод на момент времени t, минерализация пресных вод и повышенная минерализация вод; T – время подтягивания первых порций воды повышенной минерализации; t – время, на которое определяется минерализация. Прогноз техногенного изменения химического состава подземных вод при оценке миграции от загрязнителя до водозабора следует рассматривать преимущественно в виде массопереноса в подземных водах. Массопереносом в подземных водах называют движение, при котором происходит одновременное изменение количества воды, ее минерализации и химического состава, сопровождаемое различными физико-химическими процессами (ионно-сорбционным обменом, растворением, осаждением и др.). Такое совместное изменение получило название гидрогеохимического процесса или гидрогеохимическая миграция подземных вод.Линии, проведенные на гидрогеохимических картах или разрезах перпендикулярно изолиниям минерализации (изоминам) или изолиниям содержания различных компонентов химического состава (изохорны) подземных вод, показывают направление массопереноса. При изучении массопереноса предполагается, что гидродинамика потока известна, т.е. определены градиенты и скорости фильтрации, а также расход потока подземных вод. Перенос массы растворенного вещества подземными водами осуществляется тремя основными способами. 1) конвекцией, 2) диффузией и 3) гидродисперсией. Конвекция (поршневое вытеснение) возникает, когда движущаяся вода одной минерализации и состава механически, как поршень, вытесняет воду, имеющую другую минерализацию и состав. Предполагается что никаких гидрогеохимических взаимодействий между этими водами, а также водами и породой не происходит. На границе фронта вытеснения происходит резкий скачок минерализаций и химического состава воды. Скорость движения этой границы полностью зависит от скорости фильтрации v потока подземных вод. Смещение Δl (в м, км) границы раздела за время t (в сут, годах) можно определить из следующей формуле: Δl= = , где k, I, n0, v- соответственно коэффициент фильтрации, градиент потока, активная пористость и скорость фильтрации на участке, где идет перемещение границы раздела (фронта вытеснения). Так, река, (Пример по И.К. Гавич) которая загрязняет напорный водоносный горизонт, представленный песками, k = 10 м/сут, nc=0,1. Средний градиент потока I в зоне влияния водозабора равен 10-4, водозабор находится на расстоянии 100 м от реки. Определить, за сколько суток загрязненная вода, имеющая М = 5 г/л, дойдет до водозаборных скважин. Используя формулу Δ l= , запишем ее относительно t. Получим: t=Δln0/kI. Подставив исходные данные, найдем, что вода с минерализацией 5 г/л дойдет до водозаборных скважин через 104 сут, т.е. через 25 лет. Процесс конвективного переноса наблюдается, например, при внедрении загрязненных речных вод в "чистые" подземные воды при их откачке водозабором, расположенном вблизи реки. Основными параметрами процесса поршневого вытеснения являются n0 и v. При конвекции можно в обобщенном виде учесть сорбцию, тогда вместо n0 используют в расчетах nэ - эффективную пористость породы. Значения пэ могут быть самыми различными; чем сильнее выявляется влияние сорбции, тем больше эта величина и может быть больше единицы. Определяется она экспериментально. Диффузия характеризует молекулярный перенос вещества и проявляется в тех случаях, когда скорость фильтрации v воды близка или равна нулю. Перенос компонентов, растворенных в подземной воде, идет под действием градиента их концентрации и определяется уравнением Фика, аналогичным уравнению Дарси: Qc=DμFIc , где Qc - количество или поток вещества, перемещающегося через площадь поперечного сечения потока F под влиянием градиента концентрации Ic; Dm - коэффициент диффузии, зависящий от типа пород, структуры порово-трещинного пространства, влажности пород и других факторов (измеряется в м 2/сут и по физическому смыслу аналогичен коэффициенту фильтрации). Обычно он равен n 10-4 м2/сут, n=1 9. Градиент концентрации определяется формулой Iс= (С1-С2)/l , где С1 С2 - соответственно концентрации изучаемого компонента, измеренные в сечениях, находящихся на расстоянии l друг от друга. Основным параметром процесса молекулярной диффузии является Dм (определяется экспериментально или по справочной литературе). Процесс гидродисперсии более сложный и соединяет в себе как бы два предыдущих: поршневое вытеснение и сопровождающее его рассеивание вещества в порово-трещинном пространстве горных пород., возникающее вследствие неравномерной скорости движения частиц подземного раствора (молекул воды, ионов, комплексов и т.п.). В результате образуется зона дисперсии шириной 2Δхn, где наблюдается закономерное изменение концентрации компонента от С° до С0. Если рассеивание вещества осуществляемое потоком воды, идет в порах и трещинах, то называется микродисперсией, если в слоистых толщах или трещиноватых породах, разбитых сетью крупных трещин на блоки разных размеров, - макродисперсией. Влияние пород на характер рассеивания оценивается коэффициентом дисперсии D,измеряемым в м2/сут, который зависит от типа породы, структуры порово-трещинного пространства, размеров и формы частиц, трещин и блоков породы, скорости и направления фильтрации. Коэффициент микродисперсииравен n 10-5 - n 10-4 м2/сут, макродисперсии - n 10-3 - n 10-2 м2/сут и более. Оба коэффициента определяют экспериментально, первый в основном по лабораторным опытам, второй - по данным полевых опытов. Решение для одномерного процесса гидродисперсии имеет следующий вид с= = 0,5 erfcξ при ζ= , где С - относительная концентрация изучаемого компонента; С, С0, С°- концентрация компонента соответственно искомая, начальная в момент времени t = 0, максимальная, заданная на границе х = 0, до которой она мгновенно повышается или уменьшается (в зависимости от того, рассматривается процесс засоления, "загрязнения" или рассоления, "очищения"); erfcξ- табулированная функция, аналогичная функции erfcλ, равная 1 - erfξ, а erfξ, - известный интеграл вероятности; х - координата точки, находящейся на расстоянии х от границы, в которой определяется относительная концентрация компонента С на момент времени t. Размер переходной зоны (зоны дисперсии) характеризуется величиной Δх, на которую передняя граница этой зоны обгоняет фронт поршневого вытеснения, находящегося на расстоянии х0 от границы, где задана концентрация С°. Половина ширины зоны дисперсии равна Δxn=2 = 2
при x0= vt/n0
Известны также аналитические расчеты распространения загрязнения в области захвата водозабора, изложенные в работе Е.Л. Минкина ( 12 ). При стационарном режиме фильтрации нейтральная граница возмущенного потока остается постоянной во времени и химическое загрязнение было обнаружено за пределами области захвата водозабора, то загрязнение практически не может оказать своего влияния на ухудшение качества воды на водозаборном участке, за исключением тех случаев, когда загрязнение обнаружено выше по потоку от водозабора. Если в процессе эксплуатации загрязнение выявлено непосредственно в пределах площади захвата водозабора, то следует произвести расчет времени возможного его распространения и появления в эксплуатационных скважинах. Методика аналитического решения таких задач подробно изложена в работах В. Д. Бабушкина, И. С. Глазунова, В. М. Гольдберга, В. М. Шестакова, Ф.М. Бочевера (3, 4, 5, 7 ). Н.И. Плотников рассматривает в этом отношении приближенное решение задач для двух наиболее характерных случаев [14].
с1 2
х1 с0 l l l l 1 Рис. Схема к расчету времени подтягивания загрязнения в бесконечном пласте 1 — линейный ряд равнодебитных равномерно размещенных скважин; 2 — граница раздела пресных вод продуктивного горизонта с минерализацией с0 и загрязненных вод с минирализацией c1
Для линейного ряда взаимодействующих равнодебитных скважин, расположенных в неограниченном водоносном пласте на равных расстояниях друг от друга (рис.) время начала подтягивания загрязнения со стороны первоначального положения границы раздела пресных вод с загрязненными приближенно можно .определить по следующему выражению: T= , где — дебит отдельной взаимодействующей скважины; h — мощность водоносного горизонта; п0 — активная пористость пород водоносного горизонта. Это уравнение действительно при Приведенное выше уравнение отвечают таким гидрогеологическим условиям, когда водозаборное сооружение расположено на площади артезианского бассейна. Для инфильтрационных водозаборов, расположенных в долинах рек (рис. ), время начала подтягивания к водозабору фронта загрязненных вод от очага их распространения на террасовой поверхности можно определить по следующей формуле: T= , где Q — дебит отдельной скважины линейного ряда, для которой производится определение; d, — расстояние линейного ряда скважин от реки; l— расстояние между скважинами; у — расстояние от реки до контура загрязненных вод. У1 l l d
Рис.. Схема к расчету времени подтягивания загрязнения в полуограниченном пласте 1 — контур постоянного напора; г — линейный ряд равнодебитных равномерно размещенных скважин инфильтрапионного водозабора; з — граница раздела пресных и загрязненных вод Применение для расчетов приведенной формулы требует выполнения следующих условий: 1) линейный ряд скважин имеет сравнительно большую протяженность; 2) дебит скважин является одинаковым; гидравлическая связь с рекой хорошая, дебит скважин обеспечивается преимущественно за счет береговой инфильтрации поверхностных вод (притоком со стороны террасовой поверхности практически можно пренебречь); 3) расстояние между скважинами в ряду одинаковое. Если дебиты скважин линейного водозабора неодинаковы, то в расчетах можно принимать средний дебит скважин; то же относится и к расстоянию между скважинами в ряду. В том случае, если загрязненные флюиды имеют по сравнению с пресными подземными водами более высокую минерализацию и распространены на значительной площади, прогнозную оценку возможного изменения минерализации вод продуктивного горизонта можно выполнить (для линейного водозабора в неограниченном пласте) с использованием следующего уравнения: С = С0+ где С0 — минерализация пресных подземных вод; С1 — минерализация загрязненных (соленых) вод; Т — время начала подсоса соленых вод; определяется по формуле; t — расчетный момент времени, для которого определяется минерализация воды в водозаборных скважинах. При этом b= Следует иметь в виду, что выражение можно использовать при условии, когда t > Т, т.е, после начала подтягивания загрязненных вод. Предельное содержание загрязнений в воде на водозаборном участке через достаточно большой промежуток времени (t ), будет Смакс= , т.е. равно полусумме минерализации пресных и загрязненных вод/ Миграционные параметры пласта, которые используются в приведенных формулах при гидрогеологических прогнозах определяются по данным опытно-миграциооных работ (ОМР), рассмотренных ниже.[9] Интерпретация результатов ОМР, проведенная по специально разработанной методике, позволяет рассчитать активную пористость пласта n и макродисперсию δ2. Определение исходных параметров: активной пористости (при наличии физико-химических реакций взаимодействия вода-пласт n рассматривается как эффективная пористость) и макродисперсии δ2- возможно только на основе полевых миграционных опытов. Наиболее информативным методом полевых опытов определения миграционных параметров является кустовой налив. Однако его применение связано с бурение одной или нескольких наблюдательных скважин, что делает проведение опыта малоэффективным с технико-экономической точки зрения. Более простой, хотя и менее информативной, является односкажинная налив-откачка, когда в опробуемую скажину наливают трассерный раствор в течение некоторого времени. Затем из этой скажины осуществляется откачка с дебитом, превышающим расход налива. В качестве трассера при опробовании водоносных продуктивных горизонтов могут быть использованы хлор-, бром-, нитрит- ионы, радиоактивные и другие элементы. Предлагаемые трассеры могут отличаться существенно от пластовой воды как по химическим, так и по физическим характеристикам. Тогда обработка данных опыта идет не по одному, а по двум индикаторам, что существенно повышает достоверность определения искомых миграционных параметров. Методика определения искомых миграционных характеристик пласта для односкважинной налив-откачки основывается на решении радиального уравнения массопереноса, которое имеет вид при пренебрежительно малой роли молекулярной диффузии
2πnmr дс дс д2c ------- * --- + --- = δ2 ------, Q дt дr дr2 Для определения исходных параметров при значительном времени налива может быть использовано приближенное решение вышеприведенного уравнения: с = 0,5 erfc(ζ), где
Qн (tн - t) ζ = ------ * ---------- , Qот 2√δ2nt которое справедливо с погрешностью не более 5% при условии tн>0,25Qот/Qн *δ2n, где Qн и Qот – дебит налива трассера и откачки соответственно; tн – время налива. с= где С и С0 – текущие и исходные значения концентрации индикатора; С0- фоновое значение концентрации индикатора в опробуемом водоносном горизонте. Для полной обработки опытных данных по формуле C = 0,5 erfc(ξ), Строится график зависимости ζ =√tξ от t, где ξ=inferf(2C), имеющий вид прямой, пересекающей оси ζ и t в точках ζ 0 и t0. По величине этих отрезков может быть определен комплексный параметр δ2n=(t0/2 ζ)2 Имея значения пористости, по данным наблюдений за индикатором непосредственно определяют параметр дисперсии. При малом времени откачки применимо другое асимптотическое решение: C = erfc(ξ), где
Qн tн ξ= ------ * ---------- , Qот 2√δ2nt при условии t=<2,5*10-4δ2n Метод определения параметров аналогичен выше рассмотренному. Обработка полученных в процессе опыта данных выполняется в следующей последовательности: 1. По каждому определения концентрации индикатора в ходе опыта его относительная концентрация в откачиваемой воде вычисляется по формуле: С(t)= С(t) – C0 / C0 – C0, где С(t) – текущая концентрация индикатора на время t, С0—фоновая концентрация индикатора в пластовой воде до закачки рабочей жидкости; С0 – концентрация индикатора в рабочей жидкости. 2. Используя литературные данные, определяем показатель ζ 2(t)=inferf(2C(t)-1) 3. Строим график зависимости ζ 2 от t, имеющий вид прямой, точки пересечения которой с осями координат есть t0от и ζ 02.
Затем из следующих соотношений определяем параметры миграции nн,nот и δ2: где nот – активная пористость при откачке, nн – активная пористость при нагнетании, Qот-дебит откачки, Qн-расход нанетания, t0от- время откачки, определяется из графика по п.3, tн – время нагнетания,δ2-макродисперсия, ζ 02-определяется из графика по п.3. Проведение полевых ОМР с использованием односкважинной налив-откачки, несмотря на меньшую информативность по сравнению с кустовым наливом, экономически целесообразно. Интерпретация результатов ОМР по рассмотренной выше методике позволяет определить миграционные параметры пласта, учет которых значительно повышает достоверность прогнозных гидродинамических расчетов продвижения контура вод некондиционного состава. Односкажинный метод исследований (налив-откачка), опробованный В.Н. Лукиным, Д.А. Манукяном, Г.И. Потаповым на одном из месторождений промышленных вод для получения мигарционных параметорв пласта, рассмотрен в работе [9] На объектах со сложными эколого-гидрогеологическими условиями, прогнозную оценку техногенных процессов и техногенеза в целом, целесообразно выполнять методом математического моделирования. Это позволяет наиболее полно учесть все детали фильтрационного объекта. Так математическое моделирование было выполнено для прогнозной оценки биосферных и ландшафтных условий, изменяющихся в стадии длительной эксплуатации группы железорудных месторождений Курской магнитной аномалии (Лебядинского, Михайловского, Стойлинского и др.). В этом районе действует группа дренажных систем с целью осушения карьерных разработок, групповых водоразборов подземных вод, используемых для водоснабжения городов и других населенных пунктов, система хвосто-, шламо- и водохранилищ. При такой тесной компоновки водонесущих систем на сравнительно небольшой площади с помощью математического моделирования изучались: 1) балансовая структура источников обводнения горных разработок и техногенного питания подземных вод за счет неизбежных инфильтрационных потерь со стороны водонесущих сооружений; 2) влияние мощных дренажных систем на производительность действующих водозаборов для предотвращения возможного снижения их производительности и др. В общей сложности оценивались изменения ландшафтных и биосферных условий района в результате промышленного освоения большой группы железорудных месторождений. Для региональных прогнозов техногенных загрязнений подземных вод изучаемого региона необходимо иметь информацию о фоновых значениях химического состава подземных вод и длительного ряда наблюдений за режимом качества. На основании результатов региональных натурных исследований и составления прогнозных эколого-гидрохимических карт в региональном плане оценивается опасность загрязнения исследуемой территории. Эта оценка носит сравнительный приближенный характер.
Популярное: Почему человек чувствует себя несчастным?: Для начала определим, что такое несчастье. Несчастьем мы будем считать психологическое состояние... Как вы ведете себя при стрессе?: Вы можете самостоятельно управлять стрессом! Каждый из нас имеет право и возможность уменьшить его воздействие на нас... Как построить свою речь (словесное оформление):
При подготовке публичного выступления перед оратором возникает вопрос, как лучше словесно оформить свою... ©2015-2024 megaobuchalka.ru Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. (635)
|
Почему 1285321 студент выбрали МегаОбучалку... Система поиска информации Мобильная версия сайта Удобная навигация Нет шокирующей рекламы |