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

?????????? ?????? ??????? (7) ????? ????????, ??? c_ij =mr_ij , ??? ?????? r_ij ??????? ?????? ?? ?????????????? ????????????? ???????? ????? ? ???????? ???????? ???????? ?????????????? ?????????????; ?????????? ??????? r_ij ????? ??????????? ???????? ???????. ???????, ????????, ?????????? ??????? ???????? ????? ?????????? ??????? ????????? ????????, ????????

(8)

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

Покажем теперь, что тензор сопротивлений r_ij и тензор k_ij являются симметричными, т.е. r_ij = r_ij , k_{ij =} k_ij . В самом деле, на пористую среду со стороны фильтрующейся жидкости действует объемная сила, пропорциональная градиенту давления; безразмерный множитель пропорциональности зависит только от геометрических характеристик пористой среды. Удельная работа этой силы, т.е. работа за единицу времени на единицу объема системы жидкость - пористая среда, равная удельной диссипации энергии жидкостью в пористой среде, равна скалярному произведению

(9)

Очевидно, что удельная работа сил взаимодействия жидкости с пористой средой не должна зависеть от выбора осей координат х₁ , х₂ , х_3. Но для того чтобы квадратичная форма r _βα , пропорциональная этой удельной работе, не зависела от выбора системы координат, необходимо и достаточно, чтобы r _αβ = r _βα Аналогично можно показать, что k _αβ = k _βα

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

k₁₁ =k₂₂ =k; k₃₃ =k₀ ; k₁₂ =k₂₁ =k₃₂ =k₂₃ =k₃₁ = k ₁₃ =0. (10)

Закон Дарси в выбранной системе координат записывается в силу соотношений (10) следующим образом:

(11)

1.2. Зависимость параметров жидкости и пористой среды

от давления

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

1. Рассмотрим прежде всего влияние давления на свойства жидкости - плотность r и вязкость m.

Для капельных жидкостей - воды и нефти - изменения плотности обычно невелики. Встречающиеся в фильтрационных движениях перепады давления (десятки кгс/см² ) весьма малы по сравнению с модулями объемного сжатия К_r капельных жидкостей (5×10³ - 2×10⁴ кгс/см² ). Поэтому для приложений достаточно ограничиться линейной зависимостью

(12)

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

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

(13)

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

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

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

s+mp =q +p ; s = q +p (1-m ), (14)

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

q =s^f = s-p (1-m ). (15)

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