n2 = i0d/(\iw), (2.29)
так что соотношение (2.29) записывается для таких жидкостей в виде
grad р = — Ф2 , m'j. (2.30)
k \ ]XW )
Переходя в соотношении (2.30) к абсолютным величинам, находим
gradр| = -^-Ф2 , tn\. (2.31)
k V /
Заметим теперь, что если мы будем устремлять скорость фильтрации к нулю, то в пределе должна получиться величина градиента давления, не равная нулю, как в случае ньютоновской жидкости, а конечная. Эта величина называется предельным градиентом давления у. Поскольку предельный градиент давления от
скорости w не зависит, ясно, что при w -*■ 0, т. е. при больших значениях параметра т0d/[iw, функция Ф2 должна быть пропорциональна этому параметру:
Фа«б(т)^-. (2.32)
IXW
Отсюда предельный градиент давления представляется в виде у = _|=-Фз (ш), (2.33)
а закон фильтрации рассматриваемой неньютоновской жидкости (такие жидкости называются вязкопластичными) принимает вид
—*■
gradp=—w—у — . ш>0; k w
(2.34)
| grad p | < y, w = 0.
В нефтяной подземной - гидравлике этот закон был применен А. X. Мирзаджанзаде.
В последующих главах (см. гл. 6, 7) теория размерностей используется при выводе законов распределения давления для не- установившейся фильтрации упругой жидкости и газа.
Глава 3
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ФИЛЬТРАЦИИ ФЛЮИДОВ В НЕФТЕГАЗОНОСНЫХ ПЛАСТАХ
§ 1. ВВОДНЫЕ ЗАМЕЧАНИЯ
Для процессов, происходящих в нефтяных и газовых пластах, зависимость от времени существенна. Такие процессы называются неустановившимися (нестационарными). Характеристики движения — давление, скорость фильтрации и т. д.— изменяются в пласте от точки к точке; говорят поэтому, что они образуют поле. Для не- установившихся процессов элементы поля изменяются с течением времени. Задачи неустановившегося движения жидкости и газа в пласте решаются методами математической физики. Для этого составляются и интегрируются дифференциальные уравнения.
Чтобы вывести дифференциальные уравнения, в пористой среде, заключающей в себе движущийся флюид (жидкость, газ), выделяется бесконечно малый элемент и затем рассматриваются изменения массы, энергии и т. д., происходящие в нем за бесконечно малый промежуток времени. При этом используются законы сохранения массы, энергии и т. д., а также дополнительные результаты экспериментального изучения свойств и поведения флюидов в пористой среде и свойств пористой среды.
Число уравнений в системе (дифференциальных и конечных) должно равняться числу неизвестных функций, характеризующих рассматриваемый фильтрационный процесс и подлежащих определению. Такая система называется замкнутой.
В этой главе ограничимся рассмотрением процессов, для которых температура флюида равна температуре среды и неизменна: Т = const. Действительно, вследствие того что фильтрация представляет собой очень медленный процесс, изменение температуры, возникающее в ходе движения вследствие наличия сопротивления и расширения вещества, успевает компенсироваться теплообменом с окружающими горными породами. Для таких изотермических процессов уравнения энергии рассматривать уже не нужно. Однако в некоторых случаях при разработке нефтяных и газовых месторождений неизотермичность фильтрации проявляется локально в призабойной зоне скважин вследствие значительных перепадов давления. Изучение неизотермических процессов имеет особо важное значение в связи с повышением нефтеотдачи при закачке в пласт теплоносителя (горячей воды, пара). Особенности этих процессов будут рассмотрены в гл. 10.
В число дифференциальных уравнений фильтрации обязательно входят уравнение баланса массы в элементе пористой среды — уравнение неразрывности, а также дифференциальные уравнения движения. Для замыкания системы дополнительно вводятся уравнения состояния рассматриваемого флюида и пористой среды. Для получения решения системы уравнений надо еще задать условия на границах пласта и в начальный момент времени.
В результате интегрирования прежде всего определяется распределение давления и скорости фильтрации по всему пласту в любой момент времени, т. е.
р = р(х, у, z, 0, wx = wx{x, у, 2, t),
Wy = Wy{x, у, z, t), wz = wz{x, у, z, t).
Если рассматривается несжимаемая жидкость (р = const) в не- деформируемой пористой среде (т — const, k = const), то число искомых функций ограничивается этими четырьмя функциями (р, wx, wy, w2). Для фильтрации сжимаемого флюида в сжимаемой пористой среде, кроме упомянутых функций, нужно определить плотность р, вязкость (х, пористость т, проницаемость k как функции координат и времени. В этом случае нужно иметь восемь уравнений — дифференциальных и конечных — для определения восьми характеристик фильтрационного потока, жидкости и пористой среды.
Аналитическое (в виде формул) решение системы дифференциальных уравнений удается получить лишь в ограниченном числе простейших очень сильно идеализированных случаев, например в задаче о притоке упругой жидкости к скважине в пласте бесконечной протяженности с постоянным дебитом.
В более сложных случаях система уравнений решается численными методами с применением ЭВМ. В настоящее время хорошо разработаны численные методы решения самых разнообразных и очень сложных задач подземной гидравлики. При этом упомянутые аналитические решения играют очень важную роль: на них апробируются численные методы.
Систему дифференциальных уравнений можно использовать также для качественного исследования процесса. Если полученные уравнения привести к безразмерному виду, то в качестве коэффициентов будут фигурировать безразмерные параметры подобия. Анализируя их строение и численные значения, можно судить, какие силы играют решающую роль в процессе, какие члены уравнения можно отбросить и т. д.
Do'stlaringiz bilan baham: |