Опыт обработки данных опытно-фильтрационных работ в трещиноватых породах с использованием аналитических решений и численного моделирования

Размещено на http://www.allbest.ru/

Опыт обработки данных опытно-фильтрационных работ в трещиноватых породах с использованием аналитических решений и численного моделирования

Овинников А.Е.

Аннотация

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

Ключевые слова: гидрогеология, опытно-фильтрационные работы, трещиноватые породы, численное моделирование.

Abstract

The construction of an energy facility in the north of Egypt involves the implementation of experimental filtration (EFW) works to determine the water availability of rocks of near-surface aquifers. This is due to the need to build drainage systems for safe construction, as well as further operation of the facility. The results of processing the EFW on the basis of analytical solutions are given; the advantages and disadvantages of the chosen calculation scheme are indicated. The necessity to create numerical models is described. Comparison of interpretation results with analytical solutions and numerical simulation is performed.

Keywords: hydrogeology, experimental filtration work, fissured strata, numerical modeling.

Введение

В рамках инженерно-геологических изысканий для строительства энергетического объекта на севере Египта, АО «СПб НИИИ «Энергоизыскания» был выполнен ряд полевых исследований, в том числе опытно-фильтрационные работы (ОФР), включавших в себя проведение одиночных и кустовых опытных откачек [4], [7], [8]. Целью проведения ОФР являлось определение водообильности приповерхностных водоносных горизонтов.

Гидрогеологические условия рассматриваемой территории

На рассматриваемой территории выделяются два гидравлически связанных между собой водоносных горизонта [9], [10]:

? прибрежный водоносный горизонт,

? водоносный горизонт известняков.

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

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

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

Анализ результатов ОФР. Аналитическое решение

Ввиду хорошей гидравлической связи между прибрежным водоносным горизонтом и водоносным горизонтом известняков, интерпретация результатов ОФР осуществлялась по схеме двухслойного пласта [6, С. 204], которая подразумевает наличие двух водоносных пластов: верхнего слабопроницаемого и нижнего хорошо проницаемого. На исследуемой территории такими являются песчаник и известняк-ракушечник соответственно.

В рамках выбранной схемы было использовано аналитическое решение Мироненко (1,2) [5, С. 365]:

(1)

(2)

водообильность горизонт откачка фильтрация

где Q - расход опытной скважины, м3/сут; R(t) - условный радиус влияния, м; r - расстояние от опытной скважины до наблюдательной, м; T - проводимость хорошо проницаемого пласта, м2/сут; Sy - гравитационная водоотдача слабопроницаемого пласта; S - упругая водоотдача хорошо проницаемого пласта; m` - начальная обводненная мощность слабопроницаемого пласта, м; k` - коэффициент фильтрации слабопроницаемого пласта, м/сут.

Решение (1) используется для наблюдательных скважин, находящихся на расстоянии меньше условного радиуса влияния: r ? R(t).

По данным понижения уровня в основном пласте с помощью уравнения (2) определяются: проводимость и водоотдача основного пласта (T, S), гравитационная водоотдача и коэффициент фильтрации слабопроницаемого слоя (Sy, k`).

На основании данных по понижениям в скважинах опытных кустов, было выявлено, что в нескольких центральных скважинах присутствует ярко выраженная зона кальматации. В связи с тем, что выбранная расчетная схема не позволяет учесть эту зону, обработка откачек производилась только по наблюдательным скважинам и центральным скважинам, в которых такая зона отсутствует. Всего было проинтерпретировано 11 одиночных и кустовых откачек. Коэффициенты фильтрации, полученные в результате обработки представлены в табл. 1.

Таблица 1 - Коэффициенты фильтрации, полученные по результатам обработки откачек

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

Численное моделирование

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

С помощью программного комплекса Processing MODFLOW были созданы численные фильтрационные модели опытных кустов. Для каждой отдельной модели плановая схематизация была одинаковой: откачивающая скважина расположена в центре моделируемой области площадью 900 м2; на расстоянии 3,5/7 м по одному/двум лучам располагаются наблюдательные скважины. Размеры блоков - 1м ? 1м, количество блоков по вертикали и по горизонтали - 30 штук. В районе центральной и наблюдательных скважин блоки имеют более мелкую разбивку - 0,2м ? 0,2м, а на участках между скважинами - 0,5м ? 0,5м (рис. 1).

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

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

Рис. 1 - Общая схематизация моделей и граничные условия

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

Каждый выделенный модельный слой принимался однородным в плане и (коэффициент анизотропии) = 1.

Рис. 2 - Схематизация разреза на примере скважины № 3

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

Таблица 2 - Коэффициенты фильтрации по результатам моделирования на примере скважины № 11

Рис. 3 - График понижения во времени по скважине № 11

Сравнение результатов интерпретации ОФР

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

Из табл. 3 видно, что наибольшая сходимость результатов наблюдается по опытному кусту № 10.

Таблица 3 - Сравнение значений коэффициентов фильтрации

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

Выводы

1. Интерпретация ОФР по схеме двухслойного пласта позволяет определить значения коэффициентов фильтрации отдельно для слоев песчаников и известняков.

2. Моделирование опытных кустовых откачек с помощью численных методов позволило расчленить всю толщу известняков на отдельные слои с различной проницаемостью:

1) известняк-ракушечник;

2) кавернозный известняк-ракушечник;

3) трещиноватый известняк-ракушечник;

4) дресвяно-щебенистый грунт известняка-ракушечника.

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

Список литературы

1. Гриневский А. В. Оценка параметров геофильтрационного поля в условиях неопределенности. - ТУЛ ГУ. -Науки о земле. - 2016. - выпуск 3.

2. Громов В. Е. Методика обработки данных опытно-фильтрационных работ в слоистых толщах с использованием моделирования и аналитических решений / В.Е. Громов, И.С. Пашковский, Е.Ю. Потапова // Вестник Московского университета. Серия 4, Геология. - 2009. - № 6. - С. 45-51.

3. Манукьян В. А. Оценка проницаемости современных тектонических нарушений (ослабленных зон) в глинах с помощью опытно-фильтрационных работ / В.А. Манукьян, А.А. Анненков, М.Л. Глинский, С.В. Святовец // Разведка и охрана недр - 2015. - №1. - С. 39-45.

4. Мироненко В. А. Теория и методы интерпретации опытно-фильтрационных работ / В.А. Мироненко, В.М. Шестаков // - М.: Недра, 1973.

5. Синдаловский Л. Н. Справочник аналитических решений для интерпретации данных опытно-фильтрационных работ. - СПб: Издательство СПбГУ, 2006. - 768 С.

6. Синдаловский Л. Н. ANSDIMAT программный комплекс для определения параметров водоносных пластов. - СПб: Издательство «Наука», 2011 335 С.

7. Шестаков В. М. Динамика подземных вод. - М.: Издательство Московского Университета, 1973. - 327 С.

8. Шестаков В. М. Определение коэффициента фильтрации для анизотропных пластов по данным откачек // Разведка и охрана недр. - 1955. - №6. С. 52-55.

9. Mohamed Bakhbakhi Hydrogeological framework of the Nubian Sandstone Aquifer System // Proceedings of the International Workshop Tripoli, Libya, 2-4 June 2002.

Размещено на Allbest.ur