Линейный закон фильтрации Дарси
В гидродинамике рассматривается не движение отдельной частицы или слоя воды, а всей массы воды — фильтрационного потока — условного потока жидкости через пористую или пористо-трещинную среду по сообщающимся порам и трещинам. Фильтрационные потоки подземных вод различаются по характеру движения и подчиняются двум законам. Движение воды параллельно-струйчатого типа называют ламинарным, и оно подчиняется линейному закону Дарси.
Для простейших условий прямолинейно-параллельного потока линейный закон фильтрации Дарси имеет вид
Q = KфF*ΔH/ΔL
где Q — расход потока, м3/сут; Кф — коэффициент фильтрации, зависящий от свойств жидкости и фильтрующей среды, м/сут; F — площадь поперечного сечения потока, м2 ; ΔН — перепад напоров, м; ΔL — длина участка фильтрационного потока, м.
Согласно закону Дарси, количество воды Q, проходящее через трубку, заполненную дисперсным материалом, прямо пропорционально разности напоров Нв крайних сечениях трубки, прямо пропорционально площади поперечного сечения трубки F, обратно пропорционально длине пути фильтрации L и прямо пропорционально постоянному для данного материала коэффициенту Кф, характеризующему проницаемость материала, заполняющего трубку.
К основным параметрам фильтрационного потока относятся:
1. расход фильтрационного потока Q — количество воды, проходящее через поперечное сечение потока водоносного слоя за единицу времени, см3/с, м3/сути т.д.;
2. удельный расход потока q — количество воды Q, проходящее через поперечное сечение потока У7при ширине потока 1 м, м3:
F= lm,
где F — поперечное сечение потока, м2 ; l — ширина потока; m — мощность потока, м.
Подставим в формулу Дарси полученное значение:
q = Кфlm*(H1-H2)/L = Кфlm
Поскольку I= (H1-H2)/L, при ширине потока 1 м получим
q = Кфml,
где q — удельный расход потока, м3; L — длина пути фильтрации, м; Кф — коэффициент фильтрации, м/сут; H1-H2 — напор, или разность уровней в крайних сечениях потока, м; I — напорный градиент.
Произведение мощности потока на его водопроницаемость называется водопроводимостью Т потока:
Km = Т или Т = q/I м2/сут;
3. пьезометрический напор H подземных вод:
H = P/ρ + z или H = hp + z
где Р — гидростатическое давление в данной точке потока, МПа, ρ — плотность жидкости, кг/м3; z — гипсометрическая высота данной точки над выбранной плоскостью сравнения, м; Р/ρ или hp — пьезометрическая высота — та высота, на которую должна подняться вода над данной точкой потока под влиянием гидростатического давления Р в данной точке (энергия давления) (рисунок ниже), равная
hp = P/ρg = cP/ρ,
где с = 102 (переводной коэффициент значения, МПа).
Схема пьезометрического напора подземных вод (по А.И. Силину-Бекчурину)

Таким образом, пьезометрический напор — это сумма гипсометрической и пьезометрической высот, представляет собой меру энергии потока движущейся воды. При определении напора подземных вод в качестве плоскости сравнения может быть взята подошва потока или любая другая горизонтальная поверхность, например уровень Мирового океана или забой самой глубокой скважины;
4. напорный градиент (гидравлический уклон) — величина, характеризующая падение напора ΔH на единицу длины ΔL в направлении фильтрации:
I = ΔH/L = H1-H2 / L
где ΔН — перепад напоров, м; Н1 и Н2 — напоры в крайних точках потока; L — длина участка фильтрационного потока, м.
Применительно к основному закону фильтрации формула Дарси имеет вид
Q = КфFΔH/L = КфFI,
где Кф — коэффициент фильтрации, м/сут; F — площадь поперечного потока, м2; I—напорный градиент, м;
5. коэффициент фильтрации Кф— скорость фильтрации при гидравлическом уклоне, равном I = 1, характеризующий способность породы пропускать воду. На коэффициент фильтрации влияют вязкость и плотность жидкости, минеральный состав пород, температура и др. Коэффициент фильтрации для различных пород имеет разные значения; так, для очень хорошо проницаемых галечников с крупным песком, сильно закарстованных и трещиноватых пород 100-1000 м/сутки и более; для хорошо проницаемых галечников и гравия, крупного песка, среднезернистого песка, закарстованных, трещиноватых пород 100-10; проницаемых галечников и гравия, засоренных мелким песком и глиной, среднезернистого песка, сла- бозакарстованных, слаботрещиноватых пород 10—1; слабо проницаемых тонкозернистых песков, супесей, слаботрещиноватых пород 1-0,1; весьма слабопроницаемых суглинков, глин 0,1 —0,001 м/сут.
В нефтегазовой гидрогеологии коэффициент фильтрации заменяют коэффициентом проницаемости Kпр, м2:
Кф = Q/F * I = V/I или Кф = Kпр ρ/μ
где μ — вязкость жидкости, мПа * с; ρ — плотность жидкости, кг/м3; т.е. коэффициент фильтрации прямо пропорционален проница ёмости фильтрующей среды и обратно пропорционален вязкости фильтрующейся жидкости. Тогда закон Дарси принимает вид
![]()
Отсюда выразим скорость фильтрации V, м/сут, через коэффициент проницаемости:

где V— скорость фильтрации, м/сут; Кпр — коэффициент проницаемости, м2; ρ — плотность жидкости, кг/м3; μ — вязкость жидкости, мПа * с; ΔР — перепад давлений (напоров), МПа или м; L — длина пути фильтрации, м.
Коэффициенты фильтрации и проницаемости определяют в лабораторных условиях, прокачивая через образцы жидкость известной плотности и вязкости. Размерность этих величин м2 или мкм2 или дарси (Д);
6. скорость фильтрации V — количество воды, которое проходит в единицу времени через единицу поперечного сечения потока (м/сут, см/с). Скорость фильтрации Vможно получить, разделив расход потока на площадь поперечного сечения фильтрующей среды V=Q/F= KфFI/F, откуда
V=КфI.
Так как в практике гидрогеологических исследований вместо Кф используют коэффициент проницаемости породы, то скорость фильтрации определяют как произведение коэффициента проницаемости на гидравлический уклон:
V = Kпрρ / μ * i
По данной формуле определяется фиктивная скорость фильтрации, поскольку площадь поперечного сечения потока принята равной площади поперечного сечения породы. В действительности движение воды в породе происходит только по порам и площадь потока равна площади пор. Чтобы получить действительную скорость, необходимо расход воды разделить на площадь, занятую порами. Например, применительно к пескам и крупнообломочным породам:
V = Q / Fn
где Q — расход потока, м3/сут; F — площадь пор, м2; n — пористость (скважность), выраженная в долях единицы.