Взаимодействие скважин с потоком подземных вод тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

На правах рукописи

СУ ЮЧЕВ А Диляра Таировна

ВЗАИМОДЕЙСТВИЕ СКВАЖИН С ПОТОКОМ ПОДЗЕМНЫХ ВОД

01 02 05 - механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

КАЗАНЬ-2008

003169925

Работа выполнена кафедре моделирования экологических систем Казанского государственного университета им В И Ульянова-Ленина

Научный руководитель доктор физико-математических наук,

профессор

Скворцов Эдуард Викторович

Официальные оппоненты доктор физико-математических наук;'

профессор

Маклаков Дмитрий Владимирович

кандидат физико-математических наук, доцент

Шевченко Денис Вячеславович

Ведущая организация Институт механики и машиностроения

Каз НЦ РАН, г Казань

Защита состоится 19 июня 2008 г в 14 часов 30 минут на заседании диссертационного совета Д 212 081 11 при Казанском государственном университете по адресу 420008, г Казань, ул Кремлевская, 18

С диссертацией можно ознакомиться в научной библиотеке им Н И Лобачевского Казанского государственного университета

Автореферат разослан "_" мая 2008 г

Ученый секретарь диссертационное совета кандидат физ -мат наук, доцент

А А Саченков

Общая характеристика работы

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

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

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

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

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

Научная новизна результатов.

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

\

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

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

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

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

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

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

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

- при оценке риска попадания загрязнителя в нагнетательную скважину, находящуюся в загрязненном потоке

Достоверность результатов следует из того, что в рамках принятых в работе физических допущений с применением классической модели подземной гидродинамики результаты получены либо строгими аналитическими методами, либо численно-аналитически и подтверждаются их совпадением с известными для частных случаев теоретическими данными Апробация работы. По мере их получения основные результаты докладывались и обсуждались на Международной научной конференции «Геометрическая теория функций, краевые задачи и их приложения» (Казань, 18-24 марта, 2002г), на Международном семинаре «Нелинейное моделирование и управление» (г Самара, 2004), на VI республиканской научной конференции «Актуальные экологические проблемы Республики Татарстан» (г Казань, 2004), на Всероссийской конференции «Современные аспекты экологии и экологического образования» (г Казань, 2005 г), на итоговых научных конференциях КГУ (2005, 2006, 2008 г г). Работа в целом докладывалась на научном семинаре кафедры аэрогидромеханики КГУ и отделения механики НИИ математики и механики им Н Г Чеботарева КГУ (г Казань, 2007г)

Публикации. Основные результаты опубликованы в работах [1-8], список которых приведен в конце автореферата, работы [6], [8] опубликованы в изданиях, рекомендованных экспертным советом ВАК по математике и механике

Структура и объем работы. Диссертация состоит из введения, постановки задачи, шести параграфов, примечания, заключения и списка литературы Работа изтожена на 80 страницах, содержит 28 рисунков Список литературы насчитывает 35 наименований

Содержание работы

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

В этих условиях ряд авторов использовал теорию функций комплексного переменного для развития аналитических и полуаналитических методов при оценке воздействия скважин на поток Дан обзор результатов, полученных такими авторами, как Crist J А , Goltz М N , Huang J , Shan С, Javandel I, Tsang С -F , Strack О , Erdmann J В , Скворцов Э В , Кузьмин Р С

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

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

Во введении кратко описываются основные результаты работы Систематическое изложение результатов предваряет краткая общая постановка задачи.

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

Пусть X, У - координаты точки плоскости, у - проекции скорости течения, причем - V, V J = 0 - проекции скорости в бесконечно

удаленной точке При указанных выше условиях течение жидкости потенциально и описывается потенциалом скорости фильтрации Ф и функцией тока Вводится комплексный потенциал течения = Ф + ¡4', где

2 = Х + /У, / - мнимая единица При известной величине №(2) производная с= -V, + IVу определяет значение скорости течения в произвольной точке рассматриваемой области Уравнение сШ/ /¿2 = 0 позволяет найти значения координат точек, в которых скорость течения обращается в нуль (точек нулевой скорости) Соотношения IV = \¥ (2) и сНУ/сИ = 0 служат основой для анализа течения

В §1 исследуется фильтрация под действием пары источников одинакового расхода, произвольно расположенных в потоке загрязненной жидкости. В безразмерном виде приводится комплексный потенциал течения (далее <2 - безразмерный расход)

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

6

ческого расхода Q = Q0 источника и углом а, характеризующим расположение источников относительно направления потока

cos 2а - 2у,2 & '

где

1 1г -1V2

я =w L6"4"2Q«2 cos2a+1] ~ +cos2a

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

Введением параметрического переменного / = соз2а - 2у*

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

0</<1 & =

-1</<0 & =

0-У

мог

0-'!Г

1

а = — arccos

2

\V2

,,, а = — arccos

2

h-t2)vi

где

Я(г) = агс1ё О-^Г

Показано, что кривая зависимости 2 = (2о(а)близка к расположенному между осями и а отрезку прямой с угловым коэффициентом -2/я- (рис 1) Поэтому практически критический расход допустимо находить по формуле 20 = 1-2а) ж

а г 1

0 5 I 15 а

Рис 1 График зависимости критического расхода течения от угла а

В конце параграфа приводятся рисунки (см, например, рис 2), на которых изображены границы для загрязненного потока при различных значениях угла а и расхода источника

7

О

2

,1

-2

О

2

Рис 2 Границы для загрязненного потока при а = я/3 Кривые 1 - 4 соответствуют значениям параметра £3=0 7, 0 5, 0, =0 336, 0 2

В §2 рассматривается взаимодействие загрязненного потока с парой источник-сток (с равными по модулю расходами Q), также произвольно расположенными относительно направления потока Такая задача представляет интерес в связи с используемым практически способом извлечения загрязненной воды через откачивающую скважину на поверхность, ее очисткой и закачкой в пласт через нагнетательную скважину При подобной рециркуляции не требуется подводить к нагнетательной скважине воду из независимого поверхностного источника

При определении координат точек нулевой скорости оказывается, что абсциссы этих точек могут быть как положительными, так и отрицательными Найдено условие <2 = -ъ\па, при котором абсцисса меняет знак Получены уравнения, позволяющие по заданным ординатам построить проходящие через точки нулевой скорости характерные линии тока течения, которые являются барьерами для загрязненного потока При этом возможны три схемы течения, соответствующие значениям расхода в > во> 0 = во и в <во В последнем случае переток жидкости между источником и стоком прекращается (рис 3)

Получено трансцендентное уравнение для определения характерного расхода ()а как функции угла а

А=__й

где

У, = (а2 + Оа эт « + -^У2 - Оо« + ^ -81п2«

Рис 3 Схема характерных линий тока, проходящих через точки нулевой скорости а - при расходе &>£?0,б- при расходе () = £>„, в - при расходе (2<£>0

Для варианта ()<-$№ а вводится параметрическое переменное

во 2

позволяющее выразить через него величины и а

&=-

а ■

-агат-

эт.? + ¿собя

■Г + БтЯСОБЛ Л + 31П.УС03.У

для варианта £>>-Б1 па вводится параметр и такой, что

я-_м = А 0<«<7Г/2

а

В этом случае получены следующие зависимости

Б1ПМ (л--м)с03й-51ПМ

-, а = агсБШ---

л'-м-бшмсобм к — и — ъх'аиьаъи

Найдено значение угла а, при котором характерный расход достигает максимума Построена кривая зависимости 0^=0,о(°0 (Рис 4) Область, заключенная между отрезками осей = 0, а = О и кривой Яс, = О» {а) I соответствует случаю отсутствия перетока жидкости между источником и стоком

& = -

Рис 4 График зависимости характерного расхода стока 20 от угла а при взаимодействии потока с системой источник-сток

На рис 5 показаны границы областей, защищенных от загрязненного потока, при различных значениях угла а и расхода 2

Рис 5 Границы для загрязненного потока слева- при а = /г/3, кривые 1 - 4 соответствуют значениям параметра Q = 1 8,0 7, Qn= 0 188,0 1,справа-при а = -ж!3,кривые 1 -3соответствуют значениям параметра <2 = 1, Q„ =0 613,0 45

Далее определяется концентрация загрязнителя в жидкости как доля загрязненной жидкости в общем расходе стока, извлекаемой на поверхность через сток Анализируется поведение величины концентрации как функции параметров Q и а В частности, на рис 6 видны диапазоны величины Q, при которых для фиксированных значений угла а концентрация загрязнителя равна единице При увеличении параметра Q значение концентрации падает, стремясь к нулю

в жидкости, извлекаемой на поверхность через сток Кривым 1-16 соответствуют значения угла а=л/2-як/20 (¿ = 1-13), -arcsin2/;r, -9я72, -я72

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

Модуль расхода источника равен сумме расходов стоков Получены алгебраические уравнения четвертой степени, описывающие зависимость абсцисс границ-барьеров для загрязненного потока от их заданных ординат, а также трансцендентное уравнение для определения величины характерного безразмерного расхода Q0 (его смысл тот же, что и в §2) Численно найдены координаты границ для загрязненного потока при различных значениях Q, величина Q() и концентрация жидкости, извлекаемой на поверхность через стоки

В §4 изучается взаимодействие потока незагрязненных подземных вод с водозаборной скважиной, в окрестности которой расположен источник загрязнения Требуется ответить на вопрос, попадет загрязнитель в отбираемую воду или нет В соответствии с этим рассматривается течение под действием произвольно расположенных в потоке стока и источника с произвольными безразмерными расходами Q1 и -Qi соответственно

На основе выражения для комплексного потенциала течения выписаны формулы для определения безразмерных координат двух точек х,, j/, и х2,у2 нулевой скорости течения

Пусть значение расхода источника фиксировано, а расход стока достаточно мал Тогда между стоком и источником нет перетока, и в сток загрязненная жидкость не попадает Если расход стока увеличить, то при его некотором критическом значении Q2' такой переток возникает, и сущест-

вует линия тока, проходящая через обе точки нулевой скорости Функция тока в этих точках имеет одно и то же значение Отсюда следует нелинейное уравнение в комплексных величинах для определения критического расхода стока Qг при заданных значениях расхода источника загрязнения <2\ и угла а, определяющего расположение стока и источника в потоке

Если угол а = л/2, загрязненная жидкость из источника попадает в сток при любых значениях Qí и <2г При а = -7г/2 возможность загрязнения стока исследована аналитически При -п/! <а <я/2 критический расход стока находился численно В конце параграфа представлена фазовая диаграмма, демонстрирующая зависимость величины критического расхода стока (водозабора) <92' от заданного расхода источника ^ ПРИ фиксированных значениях угла а Эта диаграмма определяет диапазон изменения расхода водозабора, в котором гарантируется его защита от загрязнения (см рис 7)

Рис 7. Графики зависимости критического расхода водозабора Q1' от расхода источника загрязнения для различных углов а кривые 1 - 6 соответствуют значениям а = я/3, я/6, 0, -агс51п2/я-,-1 24,-я-/2

В §5 рассматривается задача о взаимодействии загрязненного потока с батареей п источников одинакового расхода, расположенных на прямой поперек направления потока на равных расстояниях один от другого Эта задача привлекала внимание исследователей, но решения для произвольного числа п не получила Достаточно полно изучены лишь случаи п = 2, 3, 4

При п > 4 возможны два и более прорыва загрязненного потока между источниками, и представляет интерес выяснить, где и в какой последовательности при уменьшении безразмерного расхода Q источника возникают такие прорывы Течения при четном и нечетном числе источников имеют свои особенности и потому рассматриваются отдельно

В п. 5 1 исследуется случай четного числа источников Выписано соответствующее уравнение для определения координат точек нулевой скорости течения

Оно является алгебраическим степени п относительно комплексной переменной 2 Корни этого уравнения, которые находятся численно, соответствуют координатам точек нулевой скорости При достаточно большом значении О, две такие точки располагаются на оси х, а остальные -попарно симметрично относительно оси х, и существует единственный гидродинамический барьер Г для загрязненного потока Ввиду симметрии течения рассматривается область у>0 Значение критического расхода Ql соответствует первому прорыву потока между источниками и находится из условия слияния двух действительных корней уравнения Функция тока представлена рекуррентной формулой, позволяющей вычислить эту функцию для произвольного п и затем получить алгебраическое уравнение степени п для определения абсциссы х кривой Г при заданном значении ординаты у

Дальнейший анализ течения иллюстрируется случаем п = 6 Вычисления показывают, что при = точки нулевой скорости, не находящиеся на оси х, расположены внутри барьера Г При уменьшении расхода от Ql до значения второго критического расхода ()г возникает струйка загрязненной жидкости, примыкающая к оси х, а уравнение барьера Г несколько видоизменяется Где будет осуществляться второй прорыв барьера, - между ближайшим к оси х первым и вторым либо между вторым и третьим источниками, заранее неизвестно Из анализа геометрической схемы течения следует, что при реализации второго прорыва между первым и вторым источниками должен существовать действительный корень Q системы двух нелинейных уравнений, свидетельствующий о выходе на границу Г еще одной точки нулевой скорости Вычисления показывают, что такой корень £>2 < 0 существует, так что при 2з - б < 62 возникает вторая струйка загрязненной жидкости, прорывающая барьер Значение <2 = <2г отвечает третьему прорыву между вторым и третьим источниками и находится из соответствующей системы двух нелинейных уравнений При 2 < появляется третья струйка загрязненной жидкости Для всех вариантов схем течения построены кривые - барьеры для загрязненного потока и найдены значения 0? (см рис 8)

Описанный подход к определению критических расходов источников в случае п = 6 может быть аналогичным образом применен и к обще-

му случаю п=2т Расчеты показывают, что при п = 8 первый - четвертый прорывы появляются сначала вдоль оси х, затем между первым и вторым, вторым и третьим, третьим и четвертым источниками Таким образом, прослеживается тенденция последовательного осуществления прорывов от середины батареи источников к ее периферии, причем соответствующие величины критических расходов, уменьшаясь, постепенно сближаются Отмечено, что это значение стремится к известному значению критического расхода 0 = 1/л для случая предельной схемы, когда

И-»со

Рис 8 Границы для потока в случае батареи шести источников Кривые 1-5 соответствуют значениям параметра 0=0 5, 0 383, а =0 351, д 3=0 346, 0 3

В п. 5.2 изучается случай нечетного числа источников В отличие от предыдущего случая первый прорыв барьера при у> 0 происходит между первым и вторым источниками Дальнейший анализ от описанного выше принципиально не отличается и иллюстрируется случаем и = 5 Численные расчеты с определением критических расходов проведены для я = 5, 7, 9, для случаев п- 5, 7 построены кривые-барьеры для всех возможных схем течения (см в качестве примера рис 9) Как и в случае четного числа источников, вычисления подтверждают предположение о том, что при последовательном уменьшении расхода реализуется такой «сценарий» возникновения каждого очередного прорыва исходного барьера, при котором прорывы постепенно удаляются от оси симметрии течения При этом величины критических расходов с ростом п также сближаются, стремясь к значению (2 = 1/я , соответствующему п -» оо

о

х

Рис 9 Границы для потока в случае батареи пяти источников кривые 1-5 соответствуют значениям 2 = 0 43,0 39, 2, =0 369, 02 =0 352, 0 25

В §6 рассматривается взаимодействие загрязненного потока с нагнетательной скважиной заданного радиуса, на контуре которой давление считается постоянным

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

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

Далее вводится зависящее от координат х,у и безразмерного расхода через окружность 0 параметрическое переменное С.

2_ 26

которое при 2 > 1 позволяет выразить ординату у и абсциссу х искомой границы-барьера Г через величины I и 2 следующим образом

У = Я' + (в1'2 - яш2 г)'/2, х = уО%{к - (),0 < I < л Показано, что при снижении расхода 2 до критической величины 20 = 1 загрязненный поток достигает контура скважины При этом касательная к границе Г в точке ее соприкосновения с окружностью составляет с осью х угол, равный 2л/3

При 0<2<1 загрязненный поток втекает через часть окружности (контура скважины) внутрь, а из остальной части окружности жидкость вытекает. Получены соответствующие параметрические уравнения линий, ограничивающих втекающий и вытекающий потоки При допущении, что внутри скважины потоки чистой и загрязненной жидкости равномерно перемешиваются, определена концентрация загрязнителя в потоке, вытекающем из скважины Результаты иллюстрируются графиками (см рис 10) В частности, они показывают, что при значениях безразмерного расхода 2 в диапазоне 0<<2<1 + £, £ «1 течение в окрестности скважины существенно неодномерно

0 2 0 4 0 6 0 8 Я Рис 10 Границы течения от скважины в потоке (а) для 2 = 2, 1 5, 20 = 1, 0 5, 0 001 (линии / - 5), зависимость концентрации С жидкости, вытекающей из скважины, от величины 2 (б)

Известное значение безразмерного критического расхода (50 = 1 позволяет оценить критический объемный расход скважины при котором загрязненная жидкость достигает ее контура Показано, что при достаточно малых расходах загрязнение может попасть внутрь скважины, следовательно, и в защищаемую ею область течения

Полученная оценка сопоставлена с результатами §5 Установлено, что прорывы между скважинами в батарее возможны при значениях расхода скважины, значительно больших тех, для которых загрязненный поток может достичь ее контура

В примечании к §6 в качестве возможной интерпретации исходной постановки задачи рассматривается натекание потенциального потока на контур в виде окружности с заданным на ней условием постоянства потенциала Анализируется картина характерных линий тока внутри окружности, в центре которой в соответствии с видом комплексного потенциала расположены источник и диполь

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

Основные результаты работы, выносимые на защиту, перечислены в пункте научная новизна

Список опубликованных работ по теме диссертации

1 Скворцов Э В Управление потоком подземных вод через скважины / Э В Скворцов, Д Т Суючева // Труды математического центра им Н И Лобачевского Геометрическая теория функций, краевые задачи и их приложения Материалы международной научной конференции (Казань, 18 - 24 марта 2002г) - Казань Казан ма-тем об-во, 2002 -Т 14 -С 254-259

2 Скворцов Э В Моделирование управления потоком загрязненных подземных вод / Э В Скворцов, Д Т Суючева // В сб На рубеже веков НИИ математики и механики им НГ Чеботарева КГУ 1998 — 2002 гт - Казань Изд-во Казан матем об-ва, 2003 - С 329 - 337

3 Скворцов Э В Гидробарьеры для загрязненных подземных вод / Э В Скворцов, Д Т Суючева // Нелинейное моделирование и управление Тезисы докладов международного семинара - Самара, 2004г - с 48-49

4. Скворцов Э В Моделирование гидробарьеров для потока загрязненных подземных вод / Э В Скворцов, Д Т Суючева // Актуальные экологические проблемы Республики Татарстан Тезисы докладов VI республиканской научной конференции - Казань Отечество, 2004г -С 207-208

5 Скворцов Э В Управление потоками загрязненных подземных вод через скважины / Э В Скворцов, Д Т Суючева // Ученые записки Казанского государственного университета Естественные науки -

2005 -Т 147 -Кн 1-С 21-31

6 Скворцов Э В Взаимодействие скважин с потоком подземных вод / Э В Скворцов, Д Т Суючева // Изв РАН Механика жидкости и газа -2005 -№4 - С 86-96

7 Скворцов Э В Оценка опасности загрязнения водозабора в потоке подземных вод / Э.В Скворцов, Д Т Суючева // Ученые записки Казанского государственного университета Естественные науки -

2006 -Т 148 - Кн 4 - С 100-104

8 Скворцов Э В Взаимодействие батареи скважин с потоком подземных вод / Э В Скворцов, Д Т Суючева // Экологический вестник научных центров ЧЭС -2007 -№4 -С 49-53

Отпечатано с готового оригинал-макета в типографии Издательства Казанского государственного университета Тираж 100 экз Заказ 9/5

420008, ул Профессора Нужина, 1/37 тел 231-53-59,292-65-60

Введение

Постановка задачи

§ 1. Два источника с одинаковыми расходами в потоке

§2. Источник и сток при одинаковых по модулю расходах в потоке

§3. Источник и два стока в потоке

§4. Источник и сток при произвольных расходах в потоке

§5. Батарея источников в потоке

5.1 Четное число источников в потоке

5.2 Нечетное число источников в потоке

§6. Нагнетательная скважина в потоке 67 Примечание 73 Заключение 75 Литература

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

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

Для определения влияния множества откачивающих и нагнетательных скважин на поле регионального потока широко используется численное моделирование [3-8]. Численные методы позволяют рассмотреть сложные системы, которые, в частности, включают в себя неоднородность пористой среды и нестационарность течения. Особый интерес представляет определение огибающей зоны захвата, что важно при использовании системы скважин для защиты подземных вод от загрязнения. При двумерном течении огибающая зоны захвата определяется как линия, которая отделяет воду, текущую мимо скважины, от воды, текущей к скважине. К сожалению, численные методы не позволяют осуществить прямое вычисление огибающей зоны захвата, и существует необходимость численно воспроизводить непрерывно увеличивающееся количество линий тока для более точного определения зоны захвата [9]. Показательным в этом отношении является рис. 1 из работы [2], где изображены найденные численно линии тока течения под действием пары нагнетательной и откачивающей скважин, расположенной поперёк потока: среди этих линий огибающих зон захвата нет.

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

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

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

В [18] исследовалась рециркуляция между нагнетательной и эксплуатационной скважинами, размещенными поперек направления набегающего потока. Для развития дальнейших исследований важную роль сыграла работа [19], где были определены огибающие зоны захвата в случаях двух, трех и четырех скважин в прямолинейных батареях, расположенных поперек потока. При этом варианты с прорывом течения между скважинами не рассматривались. В статье [2] в связи с активным развитием технологий локальной очистки подземных вод дан анализ течения при размещении в потоке пар источников и стоков. При этом основное внимание обращено на определение таких параметров системы, как отношение доли жидкости, поступающей из источника в сток, к полному расходу источника, и ширина защищаемой зоны. Формы границ этих зон не определялись. В работе [9] представлено аналитическое решение задачи об определении границы зоны захвата в случае пары эксплуатационных скважин, расположенных произвольно в потоке. Прорыв потока между скважинами с разделением границы зоны захвата на две не исследовался. В [20] определены те значения безразмерного критического расхода скважины в батарее скважин (до ста), при которых происходит первый в числе возможных прорыв потока между скважинами. В работе [21] предложен полуаналитический метод для определения критического расхода в случае расположения в потоке трех скважин с одинаковыми расходами в вершинах равнобедренного треугольника.

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

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

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

В зависимости от постановки конкретной задачи определение характерных расходов отвечает на следующие вопросы:

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

- в случае использования способа рециркуляции жидкости между откачивающей и нагнетательной скважинами в загрязненном потоке с целью очистки жидкости на поверхности и повторной закачки ее в пласт какова концентрация загрязнителя в откачивающей скважине и при каких расходах она достигает максимума;

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

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

В соответствии с принятой в работе схематизацией исследуемого фильтрационного течения это течение потенциально и описывается комплексным потенциалом [24, 25]. Функция тока и потенциал такого течения удовлетворяют линейному уравнению Лапласа, и при нахождении комплексного потенциала для конкретного течения может быть использован принцип суперпозиции.

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

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

Заключение

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

Дан параметрический анализ следующих фильтрационных течений:

- под действием пары источников одинакового расхода, произвольно расположенных в потоке загрязненной жидкости;

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

- при работе стока (водозабора) и действия источника загрязнения произвольного расхода в условиях природного потока;

- под действием батарей п источников {п = 5- 9), расположенных на прямой поперек загрязненного потока на равных расстояниях один от другого;

- при взаимодействии окружности с постоянным на ней потенциалом (нагнетательной скважины) и загрязненного потока.

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

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

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

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

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

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

- при оценке риска попадания загрязнителя в нагнетательную скважину, находящуюся в загрязненном потоке.

1. Мироненко В.А. Проблемы гидрогеологии. Т. 3 ( ки. 1). Прикладные исследования / В.А. Мироненко, В.Г. Румынии . - М.: Изд-во Московского государственного горного университета, 1999. -311 с.

2. Christ J.A. Development and application of an analytical model to aid design and implementation of in situ remediation technologies / J.A. Christ, M.N. Goltz, J. Huang // Journal of Contaminant Hydrology. -1999.- №37.-P. 295-317.

3. Shafer J.M. Reverse path line calculation of time-related capture zones in nonuniform flow / J.M.Shafer // Ground Water. 1987. - № 25 (3). -P. 283-289.

4. Ahlfeld D.P. Well location in capture zone design using simulation and optimization techniques / D.P. Ahlfeld, C.S. Sawyer // Ground Water. -1990. -№28 (4).-P. 507-512.

5. Bair E.S. Comparison of flow models used to delineate capture zones of wells. 1: leaky-confined fractured-carbonate aquifer / E.S. Bair, G.S. Roadcap // Ground Water. 1992. - № 30 (2). - P. 199 - 211.

6. Ratzlaff S.A. Optimal design of ground-water capture systems using segmental velocity direction constraints / S.A. Ratzlaff, M.M. Aral, F. Al-Khayyal // Ground Water. 1990. - № 30 (4). - P. 606 - 612.

7. Springer A.E. Comparison of methods used to delineate capture zones of wells: 2. Stratified-drift buried-valley aquifer / A.E. Springer, E.S. Bair // Ground Water. 1992. - № 30 (6). - P. 908 - 917.

8. Gaily R.M. Design of optimal, reliable plume capture schemes: application to the Gloucester Landfill ground-water contamination / R.M. Gaily, S.M. Gorelick // Ground Water. 1993. - № 31 (1). - P. 107-114.

9. Shan С. An analytical solution for the capture zone of two arbitrarily located wells / C. Shan // Journal of Hydrology. 1999. - №222. - P. 123 - 128.

10. Jacob C.E. Flow of groundwater / C.E. Jacob // Engineering Hydraulics, ed. Rose H. N.Y.: John Wiley and Sons. 1950. - P. 321 - 386.

11. Полубаринова-Кочина П.Я. Теория движения грунтовых вод / П.Я. Полубаринова-Кочина. М: Наука, 1977. - 664 с.

12. Бэр Я. Физико-математические основы фильтрации воды / Я. Бэр, Д. Заславски, С. Ирмей. М: Мир, 1971.-452 с.

13. Bear J. Dynamics of fluids in porous media / J. Bear. N.Y.: Amer. Els., 1972.-764 p.

14. Бочевер Ф.Н. Защита подземных вод от загрязнения / Ф.Н. Бочевер, Н.Н. Лапшин, А.Е. Орадовская. М.: Недра, 1979. - 254 с.

15. Гольдберг В.М. Гидрогеологические прогнозы качества подземных вод на водозаборах / В.М. Гольдберг. — М.: Недра, 1976. 152с.

16. Гольдберг В.М. Гидрогеологические основы охраны подземных вод от загрязнения / В.М. Гольдберг, С. Газда. М.: Недра, 1984. -262 с.

17. Strack О. Groundwater Mechanics / О. Strack. NJ: Prentice-Hall, Englewood cliffs, 1989. - P. 227 - 240.

18. Javandel I. Capture zone type curves: a tool for aquifer cleanup / I. Javandel, C.-F. Tsang // Ground Water. 1986. - №24 (5). - P. 616 -625.

19. Erdmann J.B. On capture width and capture zone gaps in multiple-well systems / J.B. Erdmann // Ground Water. 2000. - № 38 (4). - P. 497 -504.

20. Christ J.A. Hydraulic contaminent: analytical and semi-analytical models for capture zone curve delineation / J.A. Christ, M.N. Goltz // Journal of Hydrology. 2002. - № 262. - P. 224 - 244.

21. Кузьмин P.C. Влияние скважин и галерей на поток подземных вод / Р.С. Кузьмин, Э.В. Скворцов // Труды математического центра им. Н.И. Лобачевского. Краевые задачи и их приложения. Казань: Унипресс, 1999. - Т.З. - С. 78 - 84.

22. Кузьмин Р.С. Управление потоком подземных вод через галереи / Р.С. Кузьмин, Э.В. Скворцов // Труды математического центра им. Н.И. Лобачевского. Краевые задачи аэрогидромеханики и их применения. Казань: ДАС, 2000. - Т.7. - С. 196 - 262.

23. Лойцянский Л.Г. Механика жидкости и газа / Л.Г. Лойцянский. М: Наука, 1978.-736 с.

24. Басниев К.С. Подземная гидравлика / К.С. Басниев, A.M. Власов, И.М. Кочина, В.М. Максимов. -М.: Недра, 1986. 303 с.

25. Лаврентьев М.А. Методы теории функций комплексного переменного. / М.А. Лаврентьев, Б.В. Шабат. — М: Наука, 1973. -736 с.

26. Bakker М. Capture zone delineation in two-dimentional groundwater flow models / M. Bakker, O.L.D. Strack // Water Resources Research -№32 (5).- 1996.-P. 1309- 1315.

27. Скворцов Э.В. Гидробарьеры для загрязненных подземных вод / Э.В. Скворцов, Д.Т. Суючева // Нелинейное моделирование и управление: Тезисы докладов международного семинара. Самара, 2004г. - с. 48 - 49.

28. Скворцов Э.В. Управление потоками загрязненных подземных вод через скважины / Э.В. Скворцов, Д.Т. Суючева // Ученые записки Казанского государственного университета. Естественные науки.2005.-Т. 147. Кн. 1 -С. 21 -31.

29. Скворцов Э.В. Взаимодействие скважин с потоком подземных вод / Э.В. Скворцов, Д.Т. Суючева // Изв. РАН. Механика жидкости и газа. 2005. - №4. - С. 86 - 96.

30. Скворцов Э.В. Оценка опасности загрязнения водозабора в потоке подземных вод / Э.В. Скворцов, Д.Т. Суючева // Ученые записки Казанского государственного университета. Естественные науки.2006.-Т. 148.-Кн. 4.-С. 100-104.

31. Скворцов Э.В. Взаимодействие батареи скважин с потоком подземных вод / Э.В.Скворцов, Д.Т.Суючева // Экологический вестник научных центров ЧЭС. 2007. - №4. - С. 49 - 53.