Основные уравнения фильтрации жидкости и газа

Фильтрация жидкости и газа

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

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

Здесь -вектор скорости фильтрации; — плотность жидкости. Для того, чтобы получить замкнутую систему уравнений, нужно воспользоваться тем, что свойства жидкости (плотность ρ и вязкость μ), так же как и пористость и проницаемость пористой среды, являются функциями давления (мы предположим движение изотермическим). Исходя, из предположения, слабой сжимаемости жидкости и пористой среды получается уравнение 14:

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

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

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

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

,

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

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

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

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

.

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

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

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

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

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

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

Задача 3. Одинарное прямолинейно-параллельное движение. Пусть скорость фильтрации жидкости параллельны оси х и не зависит от координат y и z. Давление при этом удовлетворяет уравнению Предположим, что в плоскости х = L давление сохраняет постоянное значение, равное контурному.

При х = 0 задан расход жидкости q(t). При заданных значениях μ, k, x и L определить давление пласта.

Решение. В работе [5], предполагая, что получено решение поставленной задачи. Она описывается в виде:

Полученная формула имеет двоякий смысл. С одной стороны, он описывает распределение давления в пласте конечной длины L при малом времени. С другой стороны он дает распределение давления в произвольный момент времени. В пласте «бесконечной» протяженностью L → ∞. Дело в том, что конечное (не бесконечно малое) изменения давления распространяется за заданное время лишь на конечное расстояние и, если рассматриваются малые времена, можно считать пласт бесконечным. Решение задачи для бесконечного пласта автомодельное: независимые переменные х и t входят в решение не порознь, а лишь в комбинации

Если , то [5]

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

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

Оценим это время для систем различных режимов. При этом для χ характерное значение х = 10 4 см 2 /сек. В результате имеем: при L = 1м (переходный процесс в одном пуске-блике породы) t = 0,1сек; при L = 300м (порядка расстояния между скважинами) t = 10 4 сек »3г. L = 1км (порядка размеров месторождения) t = 10 7 сек £ 100 суткам; L = 100км (порядка размеров крупной водонапорной системы) t = 10 9 сек = 30 лет.

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

Задания для лабораторной работы.

№ п/п	Известные параметры	Определить	Число разбивания
μ, k, χ, q(t) = α0 t, t £ T, P0 (x) μ, k, χ, q(t) = α0tsin , t £ T, P0 (x) μ, k, χ, q(t) = α0ttg ,t £ T, P0 (x) μ, k, χ, q(t) = α0t 2 /(t 0 + 2) t £ T, P0 (x) μ, х, q(t) = α0t, t £ T, P (t1), P1 (t1+ Δt) k, х, q(t) = α0tsin P (t1), P (t1+ Δt) k, х, q(t) = α0t, ρ (t1), P (t1+ Δt) μ, х, q(t) = α0tsin , ρP(t1), P (t1+ Δt) μ, k, q(t) = α0t, ρ (t1), P(t1+ Δt) μ, k, q(t) = α0tsin , P(t1), P(t1+ Δt)	P (x, t) P (x, t) P (x, t) P (x, t) k, ΔP (x, t) μ, ΔP (x, t) μ, ΔP (x, t) k, ΔP (x, t) x, ΔP (x, t) x, ΔP (x, t)	n=200, 400 n=300, 600 n=400, 800 n=500, 1000 n=250, 500 n=350, 700 n=400, 800 n=500, 1000 n=500, 1000 n=600, 1200

Задача 4.Рассмотрим теперь одномерное осесимметричное (плоско-радиальное) движение при упругом режиме. Распределение давления определяется при этом как решение уравнения теплопроводности в полярных координатах (r, j):

удовлетворяющие начальному условию

а граничным условиям при r = ρ и r = R.

По прежнему основной интерес представляют решения, отвечающие стационарному начальному распределению давления

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

где – постоянная размерность давления.

Предполагая, что в работе [5] получена решение поставленной задачи в виде

(28)

Отметим важное обстоятельство: в соотношения для P (r, t) не входит радиус скважины . Это означает, что в области применяемости условия распределения не зависит от радиуса скважины.

(28) – является несобственным интегралом, с помощью замены переменной приведем его собственному интегралу. Для этого сделаем замену:

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

Задания для лабораторной работы

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

(60)

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

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

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

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

Поэтому в уравнении неразрывности (2.32) пренебрегают изменением пористости во времени и представляют это уравнение в виде

.

(2.57)

К этому уравнению необходимо присоединить уравнение состояния газа

и закон фильтрации, который при небольшой скорости фильтрации имеет вид закона Дарси

( (2.58)

где в общем случае ; — температура.

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

,

(2.59)

где — постоянные величины.

Подстановка (2.58) и (2.59) в (2.57) дает основное нелинейное уравнение теории фильтрации газа

,

(2.60)

которое впервые было получено Л. С. Лейбензоном в 1930г.

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

Тогда вместо (2.60) необходимо решить линейное уравнение

,

(2.61)

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

Нелинейный закон фильтрации

Экспериментально установлено, что иногда линейный закон фильтрации жидкости (2.58) нарушается и зависимость между и принимает вид выпуклой или вогнутой кривой, как показано на рис. 11.

Рис. 14.Возможные виды нелинейного закона фильтрации

Основные причины проявления нелинейных эффектов следующие:

1. высокая скорость фильтрации, когда параметр Рейнольдса превышает критическое значение (зависимость изображена кривой 1 на рис. 11);

2. ламинарная фильтрация жидкостей с неньютоновскими свойствами (кривая 2);

3. малая скорость фильтрации в слабопроницаемых и неоднородных пластах (кривая 2).

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

,

(2.62)

а кривая 2 – законом фильтрации с предельным градиентом

(2.63)

где, по данным Е. М. Минского, , а, по данным Б. И. Султанова, ; — эффективный диаметр пор; — предельное напряжение сдвига.

В общем случае к обоим типам кривых применимы степенная и кусочно-линейная аппроксимации

,

(2.64)

,

(2.65)

которыми удобно пользоваться при расчетах. Здесь — параметры модели; — характерное значение градиента давления; — безразмерная функция, описывающая ломаную линию (см. рис. 14).

Вопросы к 3-ему разделу

4. Силы обуславливающие фильтрацию

5. Поровое пространство осадочных горных пород

6. Объект изучения в теории фильтрации

7. Вектор скорости фильтрации

8. Компоненты скорости фильтрации

9. Соотношение скорости фильтрации и истинной скорости движения жидкости

10.Виды пород в зависимости от степени влияния трещин на фильтрацию жидкости

12. Раскрытие трещин

13.Объемная плотность трещиноватости

14. Поверхностная плотность трещиноватости

18.Основные факторы влияющие на характер закона фильтрации

19.Модель ламинарной фильтрации ньютоновских однородных жидкостей в изотропной среде

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

21.Классическое уравнение теории фильтарции

22. Коэффициент пьезопроводности

23. Уравнение Лапласа

24. Модель ламинарной фильтрации ньютоновских однородных жидкостей для анизотропной среды

25.Обобщенный закон Дарси

26.Закономерности фильтрации жидкости в трещиновато-пористых пластах для однородной и изотропной среды

27. Распределение давления в трещинах

28. Основное значение при изучении фильтрации газа

29. Основное нелинейное уравнение теории фильтрации газа

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