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

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

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

. (2.4)

В проекциях на оси координат:

- дифференциальные уравнения движения

Если ось z направлена вертикально вверх, то дифференциальные уравнения движения примут вид:

. (2.5)

Определим из системы уравнений (2.5) производные , , :

После подстановки полученных значений в уравнение неразрывности (2.3’) получим:

,

откуда

. (2.6)

Уравнение (2.6) является дифференциальным уравнением установившейся фильтрации несжимаемой жидкости по закону Дарси в недеформируемой пористой среде и называется уравнением Лапласа.

Основное уравнение фильтрации выводится путем объединения уравнения неразрывности, уравнений движения и, если флюид сжимаем, уравнений состояния. В наиболее общем виде оно записывается в виде:

. (2.7)

Формула (2.4) справедлива для изотропной среды, для которой характерно постоянство проницаемости К по всем направлениям в окрестности рассматриваемой точки.

Пористые среды, в которых коэффициент проницаемости зависит от направления потока, называются анизотропными. Большинство пород-коллекторов имеют отчетливую слоистую структуру. Если плоскость xy совместить с плоскостью слоя, а координатную ось z направить перпендикулярно, то закон Дарси можно записать в виде:

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

-

Здесь К – проницаемость параллельно напластованию; К_z – проницаемость перпендикулярно напластованию.

Уравнения состояния

Дифференциальные уравнения (2.3) и (2.5) содержат коэффициенты плотности r и вязкости m флюида, а также коэффициенты пористости m и проницаемости К породы. Зависимости этих параметров от давления называются уравнениями состояния флюида и пористой среды.

Параметр	Экспоненциальная зависимость	Линейная зависимость
Плотность
Вязкость
Пористость
Проницаемость

r₀, m₀, m₀, К₀ – значения параметров при фиксированном (начальном давлении Р₀;

r, m, m, К – значения параметров при текущем давлении Р;

b_ж – коэффициент объемного сжатия жидкости, Па^-1;

b_с – коэффициент объемной упругости среды, Па^-1;

a_m, a_m, a_к – коэффициенты, определяемые экспериментально и зависящие от свойств жидкости и породы, Па^-1.

Уравнения состояния пласта и насыщающих его флюидов замыкают систему дифференциальных уравнений.