Повышение эффективности природоохранных мероприятий требует, в свою очередь, развитие и совершенствование методов борьбы с загрязнением окружающей среды и, в частности, подземных вод углеводородами, которое является проблемой для многих регионов мира (Logeshwaran et al., 2018). Основные техногенные источники загрязнения подземных горизонтов углеводородами формируются из-за наличия дефектов в коммуникациях при добыче, транспортировке, переработке или реализации нефти и нефтепродуктов. При внештатных ситуациях углеводороды обычно распространяются от источника загрязнения латерально и вниз по профилю грунта и в случае отсутствия непроницаемых для них барьеров достигают уровня грунтовых вод, формируя в подземном пространстве отдельную фазу – углеводородную линзу. При этом часть углеводородов линзы растворяется в воде и относительно быстро мигрирует вместе с потоком грунтовых вод (Соколов и др., 2015; Garg et al., 2017; Logeshwaran et al., 2018). Таким образом, грунтовые воды становятся основным транспортным агентом загрязняющих веществ в подземном пространстве. Даже после удаления линзы свободных углеводородов миграция растворенных углеводородов с потоком грунтовых вод может продолжаться достаточно длительное время (Carey et al., 2014; Kavanaugh et al., 2013; Logeshwaran et al., 2018; Maksimovich, Khmurchik, 2015; Suthersan et al., 2015), поэтому разработка, применение и совершенствование технологий очистки грунтовых вод от растворенных углеводородов остаются приоритетной задачей.

Ликвидация углеводородного загрязнения предусматривает выполнение комплекса задач, реализацию различных методов и использование технических средств (Оборин и др., 2008; Пряничникова, 2018; Gaur et al., 2022; Naeem, Qazi, 2020). Все методы по борьбе с нефтяным загрязнением можно разделить на следующие группы: механические, физико-химические и биологические. Выделяют два основных подхода к очистке объектов окружающей среды от углеводородного загрязнения: ex situ – на специальных полигонах и in situ – на месте загрязнения (Nwachukwu, 2014), при этом применение методов in situ может быть значительно дешевле и эффективнее по сравнению с методами, требующими извлечения и транспортировки загрязненных компонентов (Juwarkar et al., 2014; Logeshwaran et al., 2018; Maksimovich et al., 2019; McHugh et al., 2014; Ossai et al., 2020).

При реализации природоохранных мероприятий любыми методами важны своевременная оценка текущего состояния объектов окружающей среды и прогноз их состояния на будущее. Использование методов гидрогеологического моделирования считается наиболее подходящим для анализа и прогноза состояния грунтовых вод (Жернов, Шестаков, 1971; Bear, Cheng, 2010; Anderson et al., 2015). Программно-инструментальные средства численного моделирования позволяют симулировать различные варианты проведения природоохранных мероприятий, изменяя расположение очистных скважин и их параметрические характеристики, что позволяет, в свою очередь, подобрать наиболее оптимальную схему реализации очистки подземных вод до сооружения технологического комплекса.

Основной целью настоящей работы являлась оценка возможности применения методов математического моделирования для разработки стратегии и определения основных параметров очистки подземных вод от растворенных нефтепродуктов с применением биотехнологического метода, использующего дозированную подачу кислорода для стимулирования жизнедеятельности углеводородокисляющих микроорганизмов подземных вод, реализованного в виде технологического комплекса. Испытания этого комплекса были проведены в районе загрязнения подземных вод растворенными нефтепродуктами (Demenev et al., 2022).

Для достижения цели в рамках исследований решались следующие задачи:

– определение основных природных характеристик исследуемого участка в ходе натурных исследований;

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

– разработка геофильтрационной и геомиграционной моделей для исследуемого участка;

– моделирование процесса очистки подземных вод биотехнологическим методом с использованием пакета RT3D (модифицированная версия MT3DMS), в том числе по различным сценариям проведения работ;

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

Участок для проведения опытных работ находится на территории Пермского края в районе развитого промышленного узла, где более 50 лет функционируют предприятия нефтехимической промышленности, и располагается непосредственно в зоне загрязнения подземных вод растворенными нефтепродуктами.

Участок расположен примерно в 100 м от небольшой реки и имеет незначительный уклон в сторону водотока. В геолого-литологическом разрезе распространены отложения нижнего отдела пермской системы. Коренные отложения повсеместно перекрыты четвертичными породами: техногенными, аллювиальными и аллювиально-делювиальными отложениями, представленными преимущественно супесями и суглинками с прослойками супеси и песка. Мощность четвертичных отложений составляет 4,5–14 м. К числу природных особенностей участка относятся неоднородный геолого-литологический состав пород, анизотропия фильтрационных свойств водовмещающих пород, переход пористой среды в трещиноватую на относительно небольшой глубине и локальное присутствие напорных подземных вод. Глубина залегания грунтовых вод меняется в зависимости от микрорельефа, а также от сезона и составляет в среднем 0,5–2,7 м. Грунтовые воды в основном приурочены к слою текучепластичных суглинков, где встречаются прослои песчаных отложений, которые преимущественно обводнены. Движение грунтовых вод четвертичных отложений направлено к местным дренам – рекам, коэффициент фильтрации грунтов достигает 0,3–0,6 м/сут.

Химический состав грунтовых вод в целом достаточно однороден, воды HCO3–Ca-типа с минерализацией около 1 г/дм3. Отмечаются повышенные концентрации растворенных нефтепродуктов, ионов железа, магния и микроэлементов (марганца, мышьяка, свинца и др.). Во всех пробах грунтовых вод прослеживается стабильное превышение ПДК по нефтепродуктам, что свидетельствует о их наличии в растворенной форме. Содержание нефтепродуктов в подземных водах испытывает сезонные колебания. Средние концентрации нефтепродуктов в пробах воды составили 5,8 мг/дм3 с максимальным содержанием 9,0 мг/ дм3. В подобных случаях проводимые природоохранные мероприятия должны предусматривать очистку подземных вод от растворенных нефтепродуктов.

Отметим, что для многих объектов топливно-энергетического комплекса проблема распространения загрязнения водорастворенными нефтепродуктами является актуальной. В таких случаях эффективными методами очистки подземного пространства являются биотехнологические методы стимуляции природных процессов деградации нефтяных углеводородов микрофлорой подземных вод (Максимович и др., 2009; Gaur et al., 2022; Logeshwaran et al., 2018; Naeem, Qazi, 2020), например с помощью дозированной подачи кислорода для стимуляции жизнедеятельности аэробных углеводородокисляющих микроорганизмов (Baldwin et al., 2010; Robertson, Cherry, 1997; Wilson et al., 1997).

В рамках натурных исследований из исследовательских скважин диаметром 60 мм были отобраны пробы грунтов и подземных вод для определения химического состава и содержания нефтепродуктов, уточнения геолого-литологического разреза территории и свойств грунтов.

В полевых условиях производили экспресс-анализ проб воды по следующим параметрам: рН (кислотно-щелочные условия), Eh (окислительно-восстановительный потенциал), TDS (Total Dissolved Solids – общее содержание растворенных веществ), удельная электрическая проводимость, температура, – с использованием комбинированного анализатора воды HI 98129 Combo (Hanna Instruments, Германия) и ОВП-метра HI 98201 (Hanna Instruments, Германия) и анализатора растворенного кислорода МАРК-302М (Россия).

В лабораторных условиях химический анализ подземных вод проводили с использованием следующего оборудования: система капиллярного электрофореза «Капель-104-Т» (Россия), pH-метр «Эксперт-рН» (Россия), спектрофотометр UNICO-1200 (США). Содержание нефтепродуктов в воде определяли флуориметрическим методом на анализаторе жидкости «Флюорат-02» (Россия) из гексанового экстракта. Исследование проб грунта в лабораторных условиях проводили по следующим показателям: физические свойства, гранулометрический состав, коэффициент фильтрации.

Обустройство опытной площадки. В настоящем исследовании основное оборудование, определяющее в дальнейшем параметры обустройства скважин, – это эмиттеры, устройства, предназначенные для рассеивания различных веществ через проницаемые трубчатые материалы. Нами использовались эмиттеры Waterloo™ (Solinst, Канада).

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

Обработка загрязненных вод с использованием эмиттеров. Физико-химические процессы играют существенную роль в самоочищении биоценозов от углеводородов нефти и нефтепродуктов, при этом полная деструкция нефтяных углеводородов возможна лишь при участии углеводородокисляющих бактерий (Оборин и др., 2008; Logeshwaran et al., 2018). Разрушение нефтяных углеводородов в аэробных условиях идет значительно быстрее, поэтому в биотехнологических методах очистки объектов окружающей среды от углеводородного загрязнения в качестве окислителя предпочтительнее использовать кислород (Коронелли, 1996; Logeshwaran et al., 2018; Tischer et al., 2013). В связи с расходом кислорода в окислительных процессах возникает необходимость поддержания его содержания на уровне, необходимом для жизнедеятельности углеводородокисляющих бактерий – не менее 5–10 мг/л. Поэтому основным контролируемым параметром в ходе проведения работ стало содержание кислорода в воде, которое с помощью эмиттеров поддерживали на уровне не ниже оптимального.

Мониторинг экологического состояния подземных вод и эффективности их очистки от загрязнения. Оценка эффективности очистки подземных вод проводилась на 30-е, 45-е, 60-е, 75-е и 92-е сутки исследования. К основному контролируемому параметру (содержанию в воде растворенного кислорода) добавлен второй – содержание нефтепродуктов. Дополнительно велись наблюдения по следующим показателям: TDS, рН, температура, удельная электрическая проводимость, уровень грунтовых вод. В качестве контрольных отобраны пробы воды до начала очистки и выполнен их анализ по этим же показателям.

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

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

На первом этапе разработана геофильтрационная модель участка, на базе которой решалась обратная задача по восстановлению естественного положения уровня подземных вод в водоносном горизонте. Геофильтрационная модель разрабатывалась с использованием программного обеспечения Processing Modflow ver 8.3. В расчетах использовалась конечно-разностная схема с применением пакета Modflow 2005 (Harbaugh, 2005).

Размеры геофильтрационной модели составили 120 × 115 м, для выполнения расчетов бралась равномерная сеточная разбивка с размерами блока 3 × 3 м. Для дополнительной детализации в районе очистных скважин проведена детализация сеточной разбивки, размеры блока составили 1 × 1 м. На западной границе модели с использованием граничного условия третьего рода задавалась река. Абсолютная отметка значения уровня грунтовых вод соответствовала результату полевого замера и составляла 108,93 м.

Для имитации естественного распределения уровней подземных вод в водоносном горизонте на западной границе модели задалось граничное условие первого рода с постоянной отметкой уровня, равной 110,9 м. Введением такого граничного условия удалось воссоздать распределение естественного потока подземных вод, направленного в сторону реки. Фильтрационные параметры водоносного горизонта задавались в соответствии с результатами лабораторных испытаний грунтов и в дальнейшем корректировались в ходе калибровки модели.

Геофильтрационная модель состоит из одного расчетного слоя, описывающего первый от поверхности водоносный горизонт. В качестве отметки кровли расчетного слоя использовался рельеф земной поверхности, мощность слоя составила 10–15 м. Калибровка модели выполнялась по данным трех наблюдательных скважин и проводилась в стационарном режиме. В ходе калибровки подбиралось значение сопротивления ложа реки на западной границе модели. По результатам калибровки удалось добиться приемлемого совпадения модельного и фактического уровней подземных вод в наблюдательных скважинах. Пространственная дискретизация и распределение уровня подземных вод на участке исследований в соответствии с геофильтрационной моделью представлены на рис. 1.

На втором этапе на основе геофильтрационной модели разработана геомиграционная модель, главной задачей которой являлась имитация естественного процесса аэробной деградации углеводородов. Геомиграционная модель переноса загрязняющих веществ в подземных водах через пористые среды создавалась на базе пакета программ MT3DMS, реактивный перенос загрязняющих веществ моделировался в трех измерениях. Пакет MT3DMS решает уравнения, включающие адвекцию, диспергирование, диффузию и химические реакции загрязняющих веществ в системе стока подземных вод в обобщенных гидрогеологических условиях, полученных путем моделирования с помощью кода MODFLOW (Zheng et al., 2012). Для моделирования процессов аэробной деградации углеводородов использовался пакет RT3D (модифицированная версия MT3DMS), описывающий многокомпонентный перенос различных загрязнителей. В программном коде учитываются механизмы конвекции, дисперсии, сорбции, распада и биодеградации по кинетике Моно. Программный код применяется для описания естественного разложения загрязнителя и усиленной биоремедации загрязненных областей (Румынин, 2020; Поташев и др., 2011).

На первом этапе геомиграционного моделирования восстанавливался уровень загрязнения водоносного горизонта углеводородами до его обработки. Так, в ходе полевых исследований по результатам опробования 10 скважин установлено, что осредненное содержание растворенных нефтепродуктов в скважинах до начала очистки составило 4–5 мг/дм3. На восточной границе модели задавалась область с постоянным содержанием растворенных углеводородов. При эпигнозном моделировании достигнуто значение содержания нефтепродуктов в скважинах, равное 4–5 мг/дм3. Значение активной пористости водоносного горизонта задавалось по экспертным данным с учетом литологического разреза и составило 0,08. Параметры гидродисперсии и молекулярной диффузии, от значений которых зависит рассеивание нефтепродуктов вдоль и поперек фильтрационного потока, ввиду их слабой изученности на участке проведения исследований также брались из опубликованных данных для территорий со схожими гидрогеологическими условиями и уточнялись в ходе калибровки (Румынин, 2020). Распределение по модели содержания растворенных в грунтовых водах нефтепродуктов до начала обработки представлено на рис. 2.

На втором этапе моделью воспроизводилась работа очистных скважин по подаче кислорода в водоносный горизонт. Модель разбивалась на 22 стресc-периода, продолжительность которых варьировала от 2 до 5 сут. Общая продолжительность работы составила 92 сут. Калибровку модели проводили по данным, полученным в ходе натурных исследований.

В ходе наблюдений за основными параметрами подземных вод, а также за параметрами неустойчивых показателей выявлены следующие закономерности для исследуемой территории. Параметр TDS, характеризующий общее содержание растворенных частиц в воде, за время проведения экспериментальных работ изменялся незначительно и в среднем составлял 593,9–765,0 ppm. Это свидетельствует о том, что воды относятся к пресным, что характерно для обрабатываемого водоносного горизонта. Показатель pH воды практически не изменялся на всем протяжении исследований и составлял в среднем 6,6–7,0, что позволяет отнести воды к нейтральным.

Температура вод первого от поверхности водоносного горизонта составляла в среднем 7,5 °С с выявленной тенденцией к ее умеренному росту, связанному с сезонными колебаниями в летний период и незначительным общим техногенным влиянием промышленной площадки, расположенной в непосредственной близости от территории исследований. Глубина залегания грунтовых вод незначительно и постепенно снижалась в среднем на 0,4 м, что также вызвано естественными сезонными колебаниями уровня вод в летний период.

Показатели содержания растворенного кислорода в воде и степени насыщения им воды в исследуемых водах являются одними из ключевых параметров при проведении мониторинга и оценке работы системы эмиттеров в целом. Чтобы избежать возможного искажения данных при отборе проб воды и их транспортировке, замер этих параметров проводили непосредственно в полевых условиях. Динамика содержания растворенного кислорода в исследуемых водах представлена в табл. 1.

До начала работ среднее содержание растворенного кислорода в воде менялось в пределах 1,0–2,0 мг/дм3. После установки системы эмиттеров и начала подачи кислорода уже на 10-й день эксперимента среднее содержание растворенного кислорода значительно возросло, достигнув максимального значения 45 мг/дм3. За время эксперимента максимальное среднее содержание растворенного кислорода в воде составило 39,2 мг/дм3, при этом насыщение воды кислородом достигало 321,5%. В ходе исследования установлены периоды снижения и стабилизации концентрации водорастворенного кислорода, связанные, скорее всего, с температурными колебаниями, разбавлением подземных вод атмосферными осадками, а также с особенностями технологического процесса подачи кислорода через систему эмиттеров.

Изменение концентрации растворенных нефтепродуктов. До начала эксплуатации технологического комплекса проводилось опробование подземных вод территории исследований с минимально возможным временным интервалом до начала работы системы. По результатам отбора среднее содержание растворенных нефтепродуктов в подземных водах исследуемого участка составило 4–5 мг/дм3. Снижение содержания нефтепродуктов в зоне очистки начинает прослеживаться в первые 30 сут подачи кислорода, при этом формируется общая тенденция к снижению содержания нефтепродуктов в зоне очистки (Demenev et al., 2022).

Довольно ярко выраженный минимум среднего содержания нефтепродуктов зафиксирован при отборе проб на 47-е сутки эксперимента в скважинах 1, 6 и 9, он равен 0,11 мг/дм3 (средняя степень очистки для системы – 95%). Последующий рост содержания нефтепродуктов может быть объяснен поступлением новой массы углеводородов из постоянного источника загрязнения или грунтов, обладающих значительной остаточной нефтенасыщенностью, при их промывании атмосферными осадками. Однако из полученных данных видно, что система очистки потенциально способна справляться с поступлением новых порций загрязнителя. Отметим, что тенденция к снижению содержания нефтепродуктов в зоне очистки в дальнейшем сохраняется на протяжении всего эксперимента, а зафиксированные особенности в целом могут быть характерны для большинства реальных объектов, где сформирован источник загрязнения и происходит постоянное поступление нефтепродуктов в подземные воды (Demenev et al., 2022). Тот факт, что разрабатываемая система очистки продолжает эффективно функционировать в подобных условиях, является ее явным преимуществом.

Результаты моделирования позволили подтвердить тенденцию снижения содержания нефтепродуктов, зафиксированную при проведении натурных работ. Так, среднее значение содержания растворенных нефтепродуктов по наиболее репрезентативным очистным скважинам 1, 6 и 9 на время окончания исследования составило 0,62 мг/ дм3. По данным, полученным на основе геомиграционной модели, среднее значение содержания углеводородов на момент заключительного временного интервала составило 0,53 мг/дм3. Для оценки скорости снижения содержания нефтепродуктов и дополнительной калибровки модели проводилось сравнение фактической и модельной динамик снижения содержания углеводородов (рис. 3).

Из рис. 3 видно, что динамика снижения содержания углеводородов в процессе аэробной биодеградации, установленная в результате моделирования, близка к фактической. Локальные отклонения данных полевых замеров от общей линии тренда можно объяснить привносом новых порций загрязнения при промывании загрязненных грунтов зоны аэрации в процессе инфильтрации осадков, процессом, который не учитывался при моделировании.

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

На заключительном этапе с помощью откалиброванной геомиграционной модели проведены прогнозные расчеты снижения содержания углеводородов при продолжающейся подаче кислорода в очистные скважины. Параметры работы очистных скважин не менялись. Прогнозные расчеты выполнялись на следующие временные точки: 120-е, 180-е и 365-е сутки с момента подачи кислорода в скважины. В результате получено распределение содержания углеводородов исследуемого участка, которое приведено на рис. 5.

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

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

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

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

Распределение содержания нефтепродуктов в подземных водах в соответствии со вторым сценарием представлено на рис. 7. Эффективность системы из шести скважин приближается к эффективности исходной системы из десяти скважин, что позволяет провести очистку водоносного горизонта вблизи скважин в короткие сроки.

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

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

Биотехнологический метод очистки подземных вод с применением дозированной подачи кислорода показал довольно высокую результативность при проведении натурных работ – по данным мониторинга содержание нефтепродуктов снизилось более чем на 90%.

Насыщение воды кислородом наблюдалось на протяжении всего эксперимента, а среднее содержание растворенного кислорода в подземных водах в зоне очистки составило 28,5 мг/дм3.

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

В результате моделирования установлено, что динамика снижения содержания углеводородов в процессе аэробной биодеградации близка к фактической: среднее значение содержания растворенных нефтепродуктов по наиболее репрезентативным очистным скважинам на время окончания исследований составило 0,62 мг/дм3, по данным моделирования – 0,53 мг/дм3.

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