Уравнения фильтрации в пористой среде

Видео:Двухфазная фильтрация в пористой средеСкачать

Двухфазная фильтрация в пористой среде

Лекции № 14, 15, 16. Установившаяся фильтрация газа в пористой среде

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

1. Дифференциальное уравнение установившейся фильтрации идеального газа.

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

Уравнение неразрывности потока:

Уравнения фильтрации в пористой среде(1)

Уравнения фильтрации в пористой среде(2)

Уравнение состояния идеального газа:

Уравнения фильтрации в пористой среде(3)

При установившейся фильтрации идеального газа

Уравнения фильтрации в пористой среде= 0 (4)

С учетом (2) – (4) уравнение (1) принимает вид:

Уравнения фильтрации в пористой среде= 0

или Уравнения фильтрации в пористой среде= 0 (5)

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

2. Прямолинейно-параллельный поток идеального газа.

Дифференциальное уравнение (5) в этом случае будет иметь вид:

Уравнения фильтрации в пористой среде= 0 (6)

Примем граничные условия:

Уравнения фильтрации в пористой средепри x=0; Уравнения фильтрации в пористой средепри x=L (7)

По аналогии с установившимся движением несжимаемой жидкости решение уравнения (6) при условиях (7) дает закон распределения давления в виде:

Уравнения фильтрации в пористой среде

или Уравнения фильтрации в пористой среде(8)

Уравнения фильтрации в пористой среде(9)

Дебит галереи, приведенный к атмосферному давлению

Уравнения фильтрации в пористой среде(10)

Средневзвешенное по объему пластовое давление

Уравнения фильтрации в пористой среде(11)

где Уравнения фильтрации в пористой среде, Уравнения фильтрации в пористой среде, Уравнения фильтрации в пористой среде.

Тогда Уравнения фильтрации в пористой среде

После интегрирования получим

Уравнения фильтрации в пористой среде(12)

3. Плоскорадиальный фильтрационный поток идеального газа.

Дифференциальное уравнение (5) будет иметь вид:

Уравнения фильтрации в пористой среде= 0 (13)

которое решается при следующих граничных условиях

Уравнения фильтрации в пористой средепри Уравнения фильтрации в пористой среде Уравнения фильтрации в пористой средепри Уравнения фильтрации в пористой среде(14)

Решение уравнения (13) имеет вид:

Уравнения фильтрации в пористой среде(15)

Уравнение (15) представляет собой закон распределения давления в пласте.

Уравнения фильтрации в пористой среде(16)

Дебит газовой скважины, приведенный к атмосферному давлению,

Уравнения фильтрации в пористой среде(17)

Средневзвешенное по объему пластовое давление определяется по (11) при Уравнения фильтрации в пористой среде. После интегрирования (11) получим, пренебрегая величиной Уравнения фильтрации в пористой среде,

Уравнения фильтрации в пористой среде(18)

4. Плоскорадиальный фильтрационный поток реального газа по закону Дарси.

Если пластовое давление выше 10МПа и депрессия не слишком мала ( Уравнения фильтрации в пористой среде0,9), то уравнение состояния газа значительно отличается от уравнения состояния идеального газа и плотность газа определяется по формуле:

Уравнения фильтрации в пористой среде(19)

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

Уравнения фильтрации в пористой среде(20)

или при малых изменениях давления

Уравнения фильтрации в пористой среде(21)

где µ0 — вязкость при фиксированном давлении;

αµ — коэффициент, определяемый экспериментально и зависящий от состава газа.

Проницаемость пласта принимается постоянной.

Дебит газовой скважины, приведенный к атмосферному давлению,

Уравнения фильтрации в пористой среде(22)

Имеется несколько способов вычисления интеграла в формуле (22), наиболее употребляем из которых следующий: по графикам зависимостей Уравнения фильтрации в пористой средеи Уравнения фильтрации в пористой средеопределяются значения Уравнения фильтрации в пористой средепеременные Уравнения фильтрации в пористой средеи z под знаком интеграла заменяются постоянными, равными Уравнения фильтрации в пористой среде, Уравнения фильтрации в пористой среде.

Тогда интеграл в формуле (22) вычисляется и (22) принимает следующий вид:

Уравнения фильтрации в пористой среде.

Основная литература: 2 117

Дополнительная литература: 4 71

Контрольные вопросы:

1. Какой газ называется идеальным?

2. Уравнения состояния идеального и реального газов.

3. Дифференциальное уравнение установившейся фильтрации идеального газа.

4. Закон распределения давления при фильтрации идеального газа.

5. Дебит газовой скважины.

6. Средневзвешенное по объему пластовое давление.

7. Коэффициент сверхжимаемости.

8. Зависимость вязкости реального газа от давления.

9. Дебит газовой скважины при фильтрации реального газа.

Дата добавления: 2016-03-10 ; просмотров: 1692 ; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ

Видео:15.10.20 | Фильтрация жидкости в однородной среде со смешанным граничным условиемСкачать

15.10.20 | Фильтрация жидкости в однородной среде со смешанным граничным условием

Модели пористых сред. Закон Дарси

Вы будете перенаправлены на Автор24

Видео:Раздел 4 Основные модели фильтрацииСкачать

Раздел 4   Основные модели фильтрации

Закон Дарси. Техническая реализация закона. Способы записи закона Дарси

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

Закон Дарси – это закономерность фильтрации газов и жидкостей в пористой среде.

Закон Дарси устанавливает линейную зависимость между гидравлическим градиентом (перепад давления, уклон) и объемным расходом газа или жидкости в таких пористых структурах, как:

  1. Мелкозернистые грунты.
  2. Глинистые грунты.
  3. Песчаные грунты.

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

Рисунок 1. Формула. Автор24 — интернет-биржа студенческих работ

где, Q — объемный расход жидкости; v — скорость фильтрации жидкости или газа; F — эффективная площадь рассматриваемого объема пористой среды (площадь поперечного сечения образца); k — коэффициент проницаемости пористой среды; р1, р2 — давления, которые созданы на концах испытуемого образца; u — абсолютная вязкость (для жидкости); L — длина фильтрующей части породы.

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

Готовые работы на аналогичную тему

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

Рисунок 2. Формула. Автор24 — интернет-биржа студенческих работ

где, P — внешнее давление; p — плотность флюида; n — динамическая вязкость флюида; g — ускорение свободного падения; z — вертикальная координата; K — коэффициент проницаемости.

Закон Дарси также возможно представить в виде уравнения баланса сил:

Рисунок 3. Формула. Автор24 — интернет-биржа студенческих работ

где, f — поле внешних сил; n — динамическая вязкость газа или жидкости; K = nk / pg — коэффициент проницаемости.

Таким образом, полную систему уравнений фильтрации, включающую условие не сжимаемости, можно записать следующим образом:

Рисунок 4. Система уравнений. Автор24 — интернет-биржа студенческих работ

При условии, что поле скорости фильтрации имеет скалярный потенциал (то есть коэффициент проницаемости постоянен), то закон Дарси можно записать в форме уравнения Лапласа:

Рисунок 5. Уравнение Лапласа. Автор24 — интернет-биржа студенческих работ

Из условия не сжимаемости получается уравнение Лапласа с граничным условием:

Рисунок 6. Уравнение. Автор24 — интернет-биржа студенческих работ

где n — вектор нормали к поверхности.

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

Рисунок 7. Схема установки. Автор24 — интернет-биржа студенческих работ

1 — насос с жидкостью или баллон с газом; 2 — вентиль; 3 — регулятор давления; 4 — входной манометр; 5 — модель пористой среды; 6 — выходной манометр; 7 — расходомер.

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

Видео:Лекция_Проницаемость_Закон_Дарси_09_10_2018Скачать

Лекция_Проницаемость_Закон_Дарси_09_10_2018

Модели пористых сред

Пористая среда – это твердое тело, пронизанное системой пустот, которые имеют нерегулярный характер и делают среду проницаемой для газов и жидкостей.

Первые исследования порового пространства осуществлялись при помощи идеализированных моделей грунта, называемых:

  1. Идеальным грунтом.
  2. Фиктивным грунтом.

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

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

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

Рисунок 8. Формула. Автор24 — интернет-биржа студенческих работ

где, di — средний диаметр i фракции; n — число частиц фракции.

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

Видео:Особенности фильтрации в трещиноватых и трещиновато-пористых средах {подземная гидрогазодинамика}Скачать

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

Реферат: Физическое описание явления фильтрации жидкости

1. Физическое описание явления фильтрации жидкости

1.1. Закон фильтрации однородной жидкости

Фильтрация представляет собой движение жидкости в пористой среде под действием перепада давления. Основной характеристикой фильтрационного движения является вектор скорости фильтрации u определяемый следующим образом. Выберем точку М пористой среды и проведем через нее элементарную площадку S . Через выделенную площадку в единицу времени протекает масса жидкости Q . Тогда проекция вектора u на нормаль к выделенной площадке равна lim Δ Q/(p S) , где p – плотность жидкости. Подчеркнем, что масса жидкости делится на полную площадь S , а не на ее часть, занятую порами.

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

Окружим точку пористой среды некоторой малой окрестностью; поле скоростей фильтрации в этой окрестности можно считать непрерывным, а все параметры пористой среды и насыщающей ее жидкости — постоянным. Нельзя пренебречь лишь изменением давления, как бы мало оно не было, поскольку при постоянном по пространству давлении движение полностью отсутствует (по существу это утверждение является основной гипотезой). Поскольку изменение давления в окрестности данной точки определяется градиентом давления, основное предположение при установлении вида закона фильтрации состоит в том, что вектор скорости фильтрации в данной точке пористой среды определяется свойствами жидкости и пористой среды и градиентом давления grad p . Пористая среда характеризуется геометрическими параметрами — характерным размером d и некоторыми безразмерными характеристиками: пористостью m , безразмерными параметрами кривой распределения и др. Закон фильтрации должен являться следствием уравнений количества движения жидкости в поровом пространстве, поэтому в систему определяющих величин следует включить также те характеристики жидкости, которые входят в эти уравнения, т.е. плотность p и вязкостьm. Таким образом, предполагается, что существует зависимость градиента давления grad p от вектора скорости фильтрации u , геометрических характеристик пористой среды m, d и т.д. и характеристик жидкости r и m. Среди величин, от которых зависит grad p , только скорость фильтрации u является вектором. В силу изотропии среды (т.е. независимости ее свойств от вращений и отражений системы отсчета) зависимость grad p от u должна быть инвариантной относительно вращения вокруг направления вектора u .

Поэтому вектор grad p должен быть направлен по одной прямой с вектором u . В самом деле, предположим обратное, т.е. пусть вектор grad p составляет некоторый угол с направлением вектора u . Если повернуть выбранную произвольную систему координат относительно направления вектора u на некоторый угол, то ни вектор u , ни какой-либо другой из определяющих параметров не изменится. Следовательно, не должен измениться и вектор grad p , зависящий только от этих параметров. Но если grad p составляет некоторый угол с направлением вектора u , то при повороте его направление относительно координатных осей обязательно изменится. Отсюда вытекает, что вектор grad p может обращен только по направлению вектора u , так что

где с — некоторая скалярная величина, зависящая от модуля вектора скорости u , а также величин d, m, p, m .

Рассмотрим сначала такие фильтрационные движения, для которых несущественны силы инерции. К числу подобных безынерционных движений принадлежит, в силу их крайней медленности, большинство фильтрационных движений, встречающихся на практике. При этом плотность р , характеризующая инерционные свойства жидкости, несущественна и исключается из числа определяющих параметров. Таким образом, при безынерционных движениях величина c зависит только от u, d, m и m. Выпишем размерности интересующих нас величин:

Уравнения фильтрации в пористой среде

Из пяти величин (2) можно выбрать три с независимыми размерностями (например, u, m, и d). Тогда, согласно p — теореме, анализа размерностей искомая зависимость будет связывать две безразмерные комбинации указанных величин. В качестве одной из безразмерных величин удобно взять пористость m, в качестве другой выберем cd 2 /m. Таким образом, имеем

cd 2 /m= f (m), c=md -2 f(m). (3)

После этого уравнение (1) может быть представлено в виде:

Уравнения фильтрации в пористой среде

Это соотношение называется законом фильтрации Дарси (по имени французского ученого, установившего его экспериментально в 1856г.). Величина k=d 2 /f(m), вводимая уравнением (4), носит название проницаемости. Проницаемость имеет разномерность площади; она не зависит от свойств жидкости и является чисто геометрической характеристикой пористой среды.

В физической системе единиц проницаемость измеряется в см 2 . Однако проницаемость большинства горных пород выражается при этом весьма малыми числами. Так, проницаемость крупнозернистых песчаников составляет 10 -8 -10 -9 см 2 ; проницаемость плотных песчаников — около 10 -10 см 2 . Ввиду этого в нефтепромысловой практике получила распространение единица проницаемости 1Д (дарси)= 1,02× 10 -8 см 2 .

В практике гидротехнических расчетов вместо давления обычно используется напор H = p/rg, и закон Дарси записывается в виде:

Уравнения фильтрации в пористой среде

Величина C, имеющая размерность скорости, называется коэффициентом фильтрации.

Функция f в выражении (3) зависит не только от пористости, но и от других безразмерных характеристик геометрии порового пространства. Были сделаны многочисленные попытки представить в качестве функции пористости и характерного размера для типичных пористых сред как путем рассмотрения простейших моделей, так и путем обработки опытных данных. Все полученные результаты носят частный характер и имеют узкую область применимости. Наибольшей известностью из формул этого рода пользуется уравнение Козени — Кармана, полученное на основе аналогии между пористой средой и системой параллельных трубок, выражающее проницаемость через удельную поверхность å и пористость m:

Уравнения фильтрации в пористой среде

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

Закон Дарси является следствием предположения о безинерционности движения жидкости. Фильтрационное течение, следующее закону Дарси, является частным случаем ползущего течения (широко известным примером ползущего течения является стоксовское обтекание сферы). Течения такого типа характеризуются преобладанием вязких сил над инерционными, т. е. очень малыми числами Рейнольдса (Re 2 ) весьма малы по сравнению с модулями объемного сжатия Кr капельных жидкостей (5×10 3 — 2×10 4 кгс/см 2 ). Поэтому для приложений достаточно ограничиться линейной зависимостью

Уравнения фильтрации в пористой среде

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

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

Уравнения фильтрации в пористой среде

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

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

2. Рассмотрим теперь вопрос, как зависят от давления жидкости свойства пористой среды — ее пористость m и проницаемость k. Обе эти величины характеризуют структуру порового пространства, и их изменение в любой точке определяется давлением жидкости и тензором напряжений, действующих в скелете пористой среды. При этом следует отметить, что в опытах определяется их зависимость не от истинных напряжений, действующих в скелете, а от некоторой их части, которую мы назовем фиктивными напряжениями. Для выяснения этого обстоятельства Уравнения фильтрации в пористой средеразберем следующую элементарную схему опыта. Пусть в цилиндрическом сосуде с площадью поперечного сечения, равной единице, находится некоторый объем пористой среды, в котором содержится жидкость под давлением p. На верхней грани этого объема лежит непроницаемый поршень, по другую сторону которого находится жидкость под тем же давлением p. В силу известного принципа гидростатики — принципа отвердевания — эта система находится в состоянии равновесия. Для выяснения зависимости пористости от нагрузки приложим к поршню дополнительную нагрузку q. Вычислим сжимающее нормальное напряжение, действующее в сечении объема пористой среды плоскостью, параллельной поршню; для этого составим уравнение равновесия части рассматриваемого объёма, ограниченной поршнем и плоскостью сечения. Пренебрегая силами трения о стенки вмещающего сосуда и собственным весом среды и жидкости, получаем

где s — истинное напряжение, действующее в пористой среде (в расчете на единицу площади общего сечения) и, очевидно, не равное приложенной нагрузке q. Изменение пористости в зависимости от давления при фиксированной нагрузке в целом мало существенное, учитывается отдельно (это изменение обусловливается сжимаемостью материала зерен, составляющих пористую среду, которая мала сравнительно со сжимаемостью пористой среды в целом, так как изменение пористости происходит в основном за счет более плотной упаковки зерен и лишь в очень небольшой мере — за счет их сжатия; если вообще не учитывать сжимаемость материала зерен, составляющих пористую среду, то пористость при фиксированной нагрузке не будет зависеть от давления жидкости). Можно показать также, что при фиксированных напряжениях s изменение давления жидкости вообще не будет приводить к изменению объема скелета, независимо от того, какова сжимаемость его материала. Таким образом, рассматриваемый опыт дает нам зависимость пористости от нагрузки q, составляющей лишь часть истинных напряжений, действующих в скелете пористой среды:

Уравнения фильтрации в пористой средеВеличину s f будем в дальнейшем называть фиктивным напряжением.

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

Аналогичные соображения применимы и в более общих случаях. Таким образом, опыт, поставленный в условиях произвольного нагружения, даст нам зависимость пористости не от тензора истинных напряжений, действующих в скелете пористой среды, а от тензора фиктивных напряжений. Ввиду того что при действии на пористую среду одного гидростатического касательные напряжения в пористой среде не возникают, касательные компоненты тензора истинных напряжений и тензора фиктивных напряжений и тензора совпадают, а нормальные компоненты отличаются на величину р(1-m), имеем

Уравнения фильтрации в пористой среде

Название: Физическое описание явления фильтрации жидкости
Раздел: Рефераты по физике
Тип: реферат Добавлен 17:10:55 23 июля 2005 Похожие работы
Просмотров: 3226 Комментариев: 21 Оценило: 6 человек Средний балл: 3.7 Оценка: 4 Скачать

(i , j =1, 2, 3,…) (16)

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

где компоненты тензора фиктивных напряжений; компоненты тензора истинных напряжений; δij =1 при j =i , δij =0 при i ?j ,

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

Уравнения фильтрации в пористой среде

где ==- главные нормальные фиктивные напряжения, аq среднее напряжение.

Величину q можно связать с давлением р, если рассматривать напряженные состояние в пласте. Пусть Н — глубина залегания пласта, h-его мощность, а r0 — средняя плотность горных пород. Обыкновенно нефтяные пласты располагаются на значительной глубине под дневной поверхностью и их мощность мала сравнительно с глубиной залегания, т. е. h 3 ), равно по порядку величины

Уравнения фильтрации в пористой среде(21)

Значение dрмакс /r0 gH обычно не превышает одной-двух десятых; величина h/Н исчезающе мала, так что изменение напряжения во всем вышележащем массиве и, в частности, на его границах мало сравнительно с исходным напряжением. Поэтому можно считать, что при изменении давления жидкости в пласте напряжения, действующие на кровле и подошве пласта, остаются постоянным.

Предыдущее рассуждение существенно основано на том, что модуль Юнга системы жидкость — пористая среда Е и модуль вышележащего массива горных пород Е1 имеют одинаковый порядок величины (что обычно имеет место в действительности). Если бы эти модули Юнга сильно отличались между собой, то выражение (21) содержало бы дополнительный множитель Е1 /Е и при Е1 >> Е отношение напряжений могло бы и не быть малым. Физически это означает, что в случае, когда вышележащая толща сложена из очень жестких пород, могут образоваться своды, и при изменении давления жидкости напряжения на кровле и подошве пласта будут меняться.

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

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Откуда (22)

Свертывая уравнения (22) (т. е. полагая i, j=1, 2, 3 и суммируя получающие уравнения), имеем

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

откуда вытекает важное соотношение

2. Основные задачи нестационарной фильтрации

2.1. Уравнение неразрывности

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

Приравнивая выражения (24) и (25) и используя формулу преобразования поверхностного интеграла в объёмный

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

откуда в силу произвольности элемента V и вытекает уравнение неразрывности

2.2. Упругий режим фильтрации

Уравнения фильтрации в пористой среде

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

Для того чтобы получить замкнутую систему уравнений, нужно воспользоваться тем, что свойства жидкости (плотность r и вязкость m), так же как и пористость и проницаемость пористой среды, являются функциями давления (мы предполагаем движение изотермическим).

Уравнения фильтрации в пористой среде

В силу (23) имеем

Уравнения фильтрации в пористой среде

исходя из предположения о слабой сжимаемости жидкости и пористой среды, можно считать относительные изменения величин r и m малыми и коэффициенты при dp/dt в предыдущих формулах постоянными:

Опытные данные показывают, что в реальных случаях

(p-p0 )/Кm 2 , а второй (dp) 2 /L 2 К. Отсюда следует, что вторым членом в принятом приближении также следует пренебречь. Таким образом, имеем

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

где коэффициент

носит название коэффициента пьезопроводности. Уравнение (29) обычно называется уравнением упругого режима или, по предложению В.Н.Щелкачева, уравнением пьезопроводности. Оно совпадает с хорошо известным классическим уравнением теплопроводности.

2. Рассмотрим постановку основных задач теории упругого режима. Определим распределение давления р в некоторой замкнутой области пространства D на протяжении промежутка времени 0 £ t£ T. Из теории уравнения теплопроводности известно, что если задать на границе Г области D линейную комбинацию давления и его производной по нормали к границе области

Уравнения фильтрации в пористой среде

и задать начальное распределение давления в области D

то существует распределение давления p(x, y, z, t), и при том единственное, удовлетворяющее уравнению (29), непрерывное в замкнутой области D, включая границу, и удовлетворяющее условия (31) и (32).

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

Рассмотрим подробнее физический смысл тех или иных дополнительных условий.

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

Уравнения фильтрации в пористой среде

откуда, используя закон Дарси, получаем (33)

На участках границы с областями, в которых перераспределения давления практически не происходит (“области питания”), давление можно считать постоянным и известным, так что

Такое условие справедливо, если, например, рассматриваемый пласт граничит с высокопроницаемой областью,

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

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

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

Уравнения фильтрации в пористой среде

Условия этого типа выполняются на тех участках границы, через которые может происходить обмен жидкости с соседними пластами через сравнительно слабопроницаемые перемычки. Если толщина перемычки D мала, а давление р¢ за ней можно считать постоянным, то расход вытекающей жидкости через участок перемычки площадью ds составит . Это количество жидкости должно быть равно

Уравнения фильтрации в пористой среде

где un — нормальная проекция скорости фильтрации на рассматриваемом участке границы. Отсюда имеем

Уравнения фильтрации в пористой среде

т.е. условия третьего рода.

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

2.3. Уравнения безнапорной фильтрации

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

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

Закон Дарси в рассматриваемом случае можно записать в виде:

Уравнения фильтрации в пористой среде

Величина С, имеющая размерность скорости, называется коэффициентом фильтрации, =- напором, а функция Сh- фильтрационным потенциалом. Заметим, что для безнапорного движения изменения давления обычно настолько малы, что пористую среду можно считать недеформируемой, а жидкость несжимаемой, так что С =const, rg = const.

В точной постановке исследование безнапорного фильтрационного движения представляет исключительные трудности математического характера; относящиеся сюда постановки задач и результаты можно найти в книге П. Я. Полубариновой-Кочиной [94]. Поэтому приходится обращаться к некоторым упрощенным постановкам.

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

Это неравенство можно переписать еще так:

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой средеВыведем важное для дальнейших рассуждений соотношение. Рассмотрим объем V, ограниченный свободной поверхностью жидкости и некоторой цилиндрической поверхностью с вертикальными образующими. Обозначим через h расстояние от свободной поверхности жидкости до водоупора, а через H- расстояние от свободной поверхности до горизонтальной плоскости z = 0; очевидно, dh/dt = dH/dt. Объем жидкости, заключенной в объеме V, равен

Уравнения фильтрации в пористой среде

где площадка S представляет собой проекцию объема на горизонтальную плоскость. Изменение количества жидкости в объеме V за бесконечно малый промежуток времени dt равно поэтому

Уравнения фильтрации в пористой среде

Вместе с тем это изменение равно притоку жидкости в объем V извне за время dt, равному

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

где g — замкнутый контур, ограничивающий площадку S, а un — нормальная компонента вектора потока , определяемого соотношением

Уравнения фильтрации в пористой среде

Приравнивая (40) и (41) и используя формулу преобразования контурного интеграла в интеграл по площади

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

Согласно закону Дарси, скорость фильтрации определяется соотношением (36)

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

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

(x , y , z , t ) = H (x , y , z , t ) + O (u z /C ); =C grad H + O (uz ).

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Таким образом, скорость можно, пренебрегая малыми величинами, вынести из-под знака интегрирования по вертикали в соотношении (42), определяющем вектор Тогда получаем

Уравнения фильтрации в пористой среде

= – Сh grad H . (45)

Представляя (45) в (44), имеем

Уравнения фильтрации в пористой среде

В это уравнение следует подставить соотношение

определяющее вертикальную координату свободной поверхности Н через ее расстояние h до водоупора и расстояние h0 от водоупора до плоскости отсчета z = 0; получим окончательное уравнение для определения h . В частности, если поверхность водоупора представляет собой горизонтальную плоскость, то ее можно принять за плоскость отсчета и, следовательно, h 0 (x ,y ) можно считать равным нулю. Тогда Н = h , и уравнение (46) принимает вид:

Уравнения фильтрации в пористой среде

Уравнения (46) и (47) были даны Буссинеском.

2.4. Основные уравнения фильтрации газа

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

Уравнения фильтрации в пористой среде

изменением пористости m во времени можно пренебречь, так что получим

Уравнения фильтрации в пористой среде

Для того чтобы получить замкнутую систему уравнений, снова нужно использовать связь плотности газа r с его давлением р и температурой Т:

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

где Т0 — постоянная температура.

Присоединяя к уравнениям (49) и (51) уравнение закона фильтрации (предполагаемого линейным)

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

В уравнении (53) r — известная функция давления. Аналогично и вязкость газа, зависящая в общем случае от давления и температуры, может быть представлена в виде:

Таким образом, и вязкость может считаться известной функцией одного лишь давления.

Введем теперь функции

Уравнения фильтрации в пористой среде

Уравнение (53) принимает при этом вид:

Уравнения фильтрации в пористой среде

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

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

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

и уравнение (56) преобразуется к виду:

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Уравнения (59) и (60) выведены в предположении постоянства температуры газа Т0 . Поэтому их обычно называют уравнениями изотермической фильтрации газа.

Уравнение (60) — основное для теории фильтрации газа — получено впервые Л. С. Лейбензоном, а затем, несколько позднее, в работе Маскета и Ботсета. Преобразования (55) также берет свое начало от работ Л. С. Лейбензона. Далее уравнение (60) совпадает с уравнением Буссинеска (47) для напора при пологих безнапорных фильтрационных движениях. Эта аналогия, впервые обнаруженная Л. С. Лейбензоном, позволяет рассматривать исследование изотермической фильтрации газа и пологих безнапорных движений несжимаемой жидкости как одну задачу.

3. НЕЛИНЕЙНЫЕ ИНВАРИАНТНЫЕ ЗАДАЧИ

ЖИДКОСТЕЙ И ГАЗОВ

3.1. Общая характеристика инвариантных задач

теории нестационарной фильтрации.

Автомодельные пологие безнапорные движения

при нулевом начальном уровне жидкости

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

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

Важным свойством рассматриваемых ниже решений является их инвариантность: для одних из этих решений — “автомодельных” — распределение давлений, напоров, плотностей и т. п. оказывается все время подобным самому себе, для других — перемещается как твердое тело с постоянной скоростью и т. д. Это свойство связано с особым характером задач, приводящих к таким решениям. Выполнение определенных преобразований зависимых и независимых переменных оставляет уравнения, граничные и начальные условия задачи неизменными. Как говорят в математике, эти задачи инвариантны относительно некоторой группы непрерывных преобразований. Такие задачи называются инвариантными, они рассматриваются ниже.

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

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

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

Пусть начальный напор жидкости в пласте равен нулю, а напор на вертикальной границе пласта изменяется по степенному закону, начиная с исходного момента t =t 0 :

где s > 0, а a — некоторая константа, которую будем выбирать в пределах –Ѕ -1 L 2 T -1 ; [t — t0 ] = T; [x] = L; [s] = [h] T — a , (64)

где через [h], L и Т обозначены соответственно размерности напора, длины и времени; константа a безразмерна. Из аргументов, от которых зависит напор жидкости, можно составить только две независимые безразмерные комбинации:

Уравнения фильтрации в пористой среде

В силу p- теоремы анализа размерностей выражение для напора можно представить в виде произведения комбинации определяющих параметров, имеющей размерность напора (в качестве нее можно взять s (t — t0 ) a ), на безразмерную функцию от безразмерных комбинаций (65). Имеем таким образом

где f — безразмерная функция, а параметр l введен вместо параметра a для удобства последующего изложения. Очевидно, что l лежит в интервале -1 2 /dx.

При непрерывной функции f(x) и f ¹ 0 требование непрерывности функции df 2 /dx = 2fdf/dx совпадает с требованием непрерывности производной df/dx. Однако при f = 0 из непрерывности df 2 /dx непрерывность df/dx не вытекает. Напротив, как увидим далее, искомая функция f(x, l) имеет в точке, где f обращается в нуль, разрыв первой производной.

Условие (69) удобнее привести к другому виду. Умножим обе части основного уравнения (62) на х и проинтегрируем по х от нуля до бесконечности. В результате получим

Уравнения фильтрации в пористой средеt5555555555

Очевидно, что dh 2 /dx стремится к нулю при х ®¥, быстрее, чем х -1 , в противном случае h не стремилось бы к нулю при х ®¥. Используя это обстоятельство и условие на бесконечности (63), получаем

Уравнения фильтрации в пористой среде

Интегрируя это соотношение в пределах от t = t 0 до t и используя граничное условие (61) и представление решения (66), имеем

Уравнения фильтрации в пористой среде

(напомним, что считаем a удовлетворяющим неравенству -1/2

Уравнения фильтрации в пористой среде

В интересующей нас области изменения a и l правая часть (71) конечна и положительна.

3.2. ПОЛОГИЕ БЕЗНАПОРНЫЕ ДВИЖЕНИЯ

С НУЛЕВЫМ НАЧАЛЬНЫМ НАПОРОМ:

ПРЕДЕЛЬНЫЕ АВТОМОДЕЛЬНЫЕ ДВИЖЕНИЯ,

ОСЕСИММЕТРИЧНЫЕ АВТОМОДЕЛЬНЫЕ ДВИЖЕНИЯ

3.2.1. Предельные автомодельные движения. Рассмотрим теперь для того же полубесконечного пласта несколько иную задачу. Будем исследовать движение на полубесконечном интервале времени (-¥, t ), поэтому начальное распределение напора по пласту несущественно.

Предположим, что на больших расстояниях от границы пласта, т.е. при х ® ¥, напор жидкости равен нулю; следовательно,

Пусть, далее, напор жидкости на границе пласта возрастет со временем по экспоненциальному закону:

h(0, t) = h0 e h t . (73)

Напор жидкости внутри пласта h (x , t) по-прежнему удовлетворяет уравнению

Уравнения фильтрации в пористой среде

Составим полный список аргументов, от которых зависит это решение. Помимо координаты х и времени t, в этот список войдут также величины h0 , == и a. Тогда размерности всех определяющих параметров решения представляются в виде:

[x]=L; [t]=T; [a]=[h] -1 L 2 T -1 ; [h0 ]=[h]; [c]=T -1 , (75)

где по-прежнему символы L, T и [h] означают соответственно размерности длины, времени и напора. Из пяти аргументов (75) с тремя независимыми размерностями можно составить две независимые комбинации, которые удобно взять в виде:

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

отсюда на основе p- теоремы решение рассматриваемой задачи будет

где j — безразмерная функция.

Положим теперь t = t¢ + t , где t — произвольная константа. При этом условие (72) и уравнение (74), как нетрудно проверить, записываются через новую переменную t¢, так же как и через прежнюю переменную, а условие (73) принимает вид:

Уравнения фильтрации в пористой среде

Таким образом, сдвиг во времени влияет лишь на некоторое преобразование величины h0 , и постановка задачи оказывается по отношению к группе преобразований переноса по времени; для определения h в переменных х, t¢, a, c, h¢0 получается та же задача, что и для определения h в переменных (75). Стало быть, на основе соотношений (76) и (77) имеем

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Отсюда следует, что при любом t имеет место тождество

Положим теперь t = t и получим

Уравнения фильтрации в пористой среде

Итак, функция h , зависящая от пяти аргументов (75), представляется через функцию одного аргумента:

Уравнения фильтрации в пористой среде

Подставляя (81) в основное уравнение (74), получаем для функции f (x) обыкновенное дифференциальное уравнение

Уравнения фильтрации в пористой среде

Подставляя выражение (81) в условие на бесконечности (72) и граничное условие (73), имеем граничные условия для функции f (x):

В силу непрерывности напора жидкости и потока жидкости функция f(x) по-прежнему должна быть непрерывной и иметь непрерывную производную от квадрата df 2 /dx. Мы получили, таким образом, для определения функции f(x) граничную задачу того же типа, что и граничные задачи для автомодельных решений, рассмотренных в предыдущем параграфе, и соответствующую значению параметра a, равному бесконечности, т. е. l = 1. Функция f(x) = f(x, 1) тождественно равна нулю при x ³ x0 = 1,810; передний фронт х0 (t) перемещается, таким образом, по закону

Уравнения фильтрации в пористой среде

а скорость его перемещения равна

Уравнения фильтрации в пористой среде

Полученное решение является в некотором смысле предельным для автомодельных решений, рассмотренных выше. В самом деле, положим в формуле (66)

где h0 — константа размерности напора; t — константа размерности времени, причем, очевидно, эти константы выбираются с точностью до некоторого постоянного множителя. Решение (66) принимает вид

Уравнения фильтрации в пористой среде

Будем неограниченно увеличивать в этом решении a при начальном моменте t0 ® — ¥ по закону

Раскрывая неопределенность, получаем, что при a ® ¥

Уравнения фильтрации в пористой среде

Уравнение (67) в пределе при a ® ¥ переходит в уравнение (82), а условия (68) и (69) совпадают с условиями (83); f (x, l) ® f (x, 1) = f (x).

Обозначая t через 1/c, получаем, что при a ® ¥ решение (87) стремится к решению (81). Поэтому решение (81) было названо предельным автомодельным решением. Это решение было получено в работе Г. И. Баренблатта. предельные автомодельные решения представляют и принципиальный интерес в том отношении, что для доказательства автомодельности этих решений уже недостаточно соображений анализа размерности, т. е. недостаточно инвариантности постановки задачи относительно группы преобразования подобия величин с независимыми размерностями, как это было в ранее рассмотренных автомодельных задачах, а требуется дополнительно воспользоваться инвариантностью постановки задачи относительно еще одной группы — группы преобразований переноса по времени.

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

Уравнения фильтрации в пористой среде

Задача. На границе х=0 полубесконечного пласта с непроницаемым горизонтальным водоупором задается поток (расход) жидкости как степенная функция времени

Начальный напор во всем пласте равен нулю.

Решение задачи представляется в виде:

Уравнения фильтрации в пористой среде

где м (l) =-df 2 (0, l)/dx , а координата переднего фронта жидкости х0 (t) — в виде:

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

где r — расстояние рассматриваемой точки пласта от оси симметрии.

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

Предположим далее, что расход закачиваемой жидкости изменяется со временем по степенному закону. Выражение для полного расхода жидкости, закачиваемой через скважину радиусом R, имеет вид:

Уравнения фильтрации в пористой среде

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

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

где t > 0 и b > -1. В частности, случай b = 0 соответствует закачке жидкости в пласт с постоянным расходом. Таким образом, решение задачи удовлетворяет уравнению (93) и условиям (94) и (96). По-прежнему, используя p- теорему анализа размерности, можно показать, что это решение является автомодельным и представляется в виде:

представляют собой две независимые безразмерные комбинации определяющих параметров решения; других независимых комбинаций этих параметров не существует. Постоянный множитель снова введен в формулу для x с целью удобства последующего изложения. Как и прежде, искомая функция должна быть непрерывной и иметь непрерывную производную от квадрата. Подставляя выражение (77) в уравнение (93) и условия (94) и (96), находим, что функция f1 (x, l) удовлетворяет уравнению

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

Исследование этой граничной задачи проводится аналогично предыдущему; также единственным образом строится функция f1 (x, l), отличающаяся от нуля лишь при 0 £ x £ x1 (l), где x1 (l) — некоторая функция x, а при x ³ x1 (l) тождественно равная нулю. Функция f1 (x, l) при x®0 имеет особенность, как нетрудно видеть из первого условия (100):

Уравнения фильтрации в пористой среде

Второе условие (100) может быть приведено к другой форме: умножая уравнение (99) на x и интегрируя в пределах от x = 0 до x = ¥, получаем, используя оба условия (100) и условия

Уравнения фильтрации в пористой среде

Уравнения фильтрации в пористой среде

следующее интегральное соотношение:

Первое условие (102) непосредственно следует из условия, которому функция f1 (x, l) на бесконечности, так как если бы предел == при x ® ¥ не был равен нулю, то функция f1 (x, l) не стремилась бы к нулю при x ® ¥. Второе условие (102) непосредственно следует из (101).

Эффективное вычисление функции f1 (x, l) удобно проводить следующим образом. Строим решение задачи Коши Ф1 (x, l) для уравнения (99). обращающееся в нуль при x = 1 и имеющее в этой точке конечную первую производную. Исследование, в точности аналогичное приведенному в п. 3 x1, показывает, что эта производная равна -1/4. Строить решение задачи Коши удобно так: вблизи x = 1 можно представить решение в виде ряда, при помощи которого находится надлежащее число начальных значений, после чего применяется метод численного интегрирования Адамса — Штермера. Далее численно вычисляется величина

📺 Видео

МОДЕЛИРОВАНИЕ МНОГОФАЗНЫХ ФИЛЬТРАЦИОННЫХ ТЕЧЕНИЙСкачать

МОДЕЛИРОВАНИЕ МНОГОФАЗНЫХ ФИЛЬТРАЦИОННЫХ ТЕЧЕНИЙ

Урок 132. Основные понятия гидродинамики. Уравнение непрерывностиСкачать

Урок 132. Основные понятия гидродинамики. Уравнение непрерывности

Определение коэффициента фильтрации неуплотненного песчаного грунтаСкачать

Определение коэффициента фильтрации неуплотненного песчаного грунта

Гидродинамика, теплопередача и механика грунтов в COMSOL Multiphysics®. Версия 5.6Скачать

Гидродинамика, теплопередача и механика грунтов в COMSOL Multiphysics®. Версия 5.6

К.Д.Толстыгин «Сравнение метода поро-сетевых моделей и уравнения ван Генухтена-Муалема»Скачать

К.Д.Толстыгин «Сравнение метода поро-сетевых моделей и уравнения ван Генухтена-Муалема»

Механика и термодинамика пористой средыСкачать

Механика и термодинамика пористой среды

Процесс фильтрования | Лекция по ПАППСкачать

Процесс фильтрования | Лекция по ПАПП

Закон БернуллиСкачать

Закон Бернулли

Эффект Вентури и трубка Пито (видео 16) | Жидкости | ФизикаСкачать

Эффект Вентури и трубка Пито (видео 16) | Жидкости  | Физика

Подземная гидромеханика Лекция 9 Ильдар ИбрагимовСкачать

Подземная гидромеханика Лекция 9 Ильдар Ибрагимов

Супер ЖЕСТЬ для продвинутых: sqrt(2+sqrt(2-sqrt(2+x)))=xСкачать

Супер ЖЕСТЬ для продвинутых: sqrt(2+sqrt(2-sqrt(2+x)))=x

10 класс, 43 урок, Уравнение касательной к графику функцииСкачать

10 класс, 43 урок, Уравнение касательной к графику функции

Вязкость и течение Пуазёйля (видео 14) | Жидкости | ФизикаСкачать

Вязкость и течение Пуазёйля (видео 14) | Жидкости  | Физика
Поделиться или сохранить к себе: