Основы теории фильтрации доклад

Обновлено: 19.05.2024

Фильтрация представляет собой движение жидкости в пористой среде под действием перепада давления. Основной характеристикой фильтрационного движения является вектор скорости фильтрации u определяемый следующим образом. Выберем точку М пористой среды и проведем через нее элементарную площадку S . Через выделенную площадку в единицу времени протекает масса жидкости Q . Тогда проекция вектора u на нормаль к выделенной площадке равна lim Δ Q/(p S) , где p – плотность жидкости. Подчеркнем, что масса жидкости делится на полную площадь S , а не на ее часть, занятую порами.

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

Окружим точку пористой среды некоторой малой окрестностью; поле скоростей фильтрации в этой окрестности можно считать непрерывным, а все параметры пористой среды и насыщающей ее жидкости - постоянным. Нельзя пренебречь лишь изменением давления, как бы мало оно не было, поскольку при постоянном по пространству давлении движение полностью отсутствует (по существу это утверждение является основной гипотезой). Поскольку изменение давления в окрестности данной точки определяется градиентом давления, основное предположение при установлении вида закона фильтрации состоит в том, что вектор скорости фильтрации в данной точке пористой среды определяется свойствами жидкости и пористой среды и градиентом давления grad p . Пористая среда характеризуется геометрическими параметрами - характерным размером d и некоторыми безразмерными характеристиками: пористостью m , безразмерными параметрами кривой распределения и др. Закон фильтрации должен являться следствием уравнений количества движения жидкости в поровом пространстве, поэтому в систему определяющих величин следует включить также те характеристики жидкости, которые входят в эти уравнения, т.е. плотность p и вязкостьm. Таким образом, предполагается, что существует зависимость градиента давления grad p от вектора скорости фильтрации u , геометрических характеристик пористой среды m, d и т.д. и характеристик жидкости r и m. Среди величин, от которых зависит grad p , только скорость фильтрации u является вектором. В силу изотропии среды (т.е. независимости ее свойств от вращений и отражений системы отсчета) зависимость grad p от u должна быть инвариантной относительно вращения вокруг направления вектора u .

Поэтому вектор grad p должен быть направлен по одной прямой с вектором u . В самом деле, предположим обратное, т.е. пусть вектор grad p составляет некоторый угол с направлением вектора u . Если повернуть выбранную произвольную систему координат относительно направления вектора u на некоторый угол, то ни вектор u , ни какой-либо другой из определяющих параметров не изменится. Следовательно, не должен измениться и вектор grad p , зависящий только от этих параметров. Но если grad p составляет некоторый угол с направлением вектора u , то при повороте его направление относительно координатных осей обязательно изменится. Отсюда вытекает, что вектор grad p может обращен только по направлению вектора u , так что

grad p= - си , (1)

где с - некоторая скалярная величина, зависящая от модуля вектора скорости u , а также величин d, m, p, m .

Рассмотрим сначала такие фильтрационные движения, для которых несущественны силы инерции. К числу подобных безынерционных движений принадлежит, в силу их крайней медленности, большинство фильтрационных движений, встречающихся на практике. При этом плотность р , характеризующая инерционные свойства жидкости, несущественна и исключается из числа определяющих параметров. Таким образом, при безынерционных движениях величина c зависит только от u, d, m и m. Выпишем размерности интересующих нас величин:

Из пяти величин (2) можно выбрать три с независимыми размерностями (например, u, m, и d). Тогда, согласно p - теореме, анализа размерностей искомая зависимость будет связывать две безразмерные комбинации указанных величин. В качестве одной из безразмерных величин удобно взять пористость m, в качестве другой выберем cd 2 /m. Таким образом, имеем

cd 2 /m= f (m), c=md -2 f(m). (3)

После этого уравнение (1) может быть представлено в виде:

Это соотношение называется законом фильтрации Дарси (по имени французского ученого, установившего его экспериментально в 1856г.). Величина k=d 2 /f(m), вводимая уравнением (4), носит название проницаемости. Проницаемость имеет разномерность площади; она не зависит от свойств жидкости и является чисто геометрической характеристикой пористой среды.

В физической системе единиц проницаемость измеряется в см 2 . Однако проницаемость большинства горных пород выражается при этом весьма малыми числами. Так, проницаемость крупнозернистых песчаников составляет 10 -8 -10 -9 см 2 ; проницаемость плотных песчаников - около 10 -10 см 2 . Ввиду этого в нефтепромысловой практике получила распространение единица проницаемости 1Д (дарси)= 1,02× 10 -8 см 2 .

В практике гидротехнических расчетов вместо давления обычно используется напор H = p/rg, и закон Дарси записывается в виде:

Величина C, имеющая размерность скорости, называется коэффициентом фильтрации.

Функция f в выражении (3) зависит не только от пористости, но и от других безразмерных характеристик геометрии порового пространства. Были сделаны многочисленные попытки представить в качестве функции пористости и характерного размера для типичных пористых сред как путем рассмотрения простейших моделей, так и путем обработки опытных данных. Все полученные результаты носят частный характер и имеют узкую область применимости. Наибольшей известностью из формул этого рода пользуется уравнение Козени - Кармана, полученное на основе аналогии между пористой средой и системой параллельных трубок, выражающее проницаемость через удельную поверхность å и пористость m:

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

Закон Дарси является следствием предположения о безинерционности движения жидкости. Фильтрационное течение, следующее закону Дарси, является частным случаем ползущего течения (широко известным примером ползущего течения является стоксовское обтекание сферы). Течения такого типа характеризуются преобладанием вязких сил над инерционными, т. е. очень малыми числами Рейнольдса (Re 2 ) весьма малы по сравнению с модулями объемного сжатия К_r капельных жидкостей (5×10 3 - 2×10 4 кгс/см 2 ). Поэтому для приложений достаточно ограничиться линейной зависимостью

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

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

Преимущества такой записи связаны с тем, что для функции 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 будем в дальнейшем называть фиктивным напряжением.

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

Аналогичные соображения применимы и в более общих случаях. Таким образом, опыт, поставленный в условиях произвольного нагружения, даст нам зависимость пористости не от тензора истинных напряжений, действующих в скелете пористой среды, а от тензора фиктивных напряжений. Ввиду того что при действии на пористую среду одного гидростатического касательные напряжения в пористой среде не возникают, касательные компоненты тензора истинных напряжений и тензора фиктивных напряжений и тензора совпадают, а нормальные компоненты отличаются на величину р(1-m), имеем

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

где ==- главные нормальные фиктивные напряжения, аq среднее напряжение.

Величину q можно связать с давлением р, если рассматривать напряженные состояние в пласте. Пусть Н - глубина залегания пласта, h-его мощность, а r₀ - средняя плотность горных пород. Обыкновенно нефтяные пласты располагаются на значительной глубине под дневной поверхностью и их мощность мала сравнительно с глубиной залегания, т. е. h 3 ), равно по порядку величины

(21)

Значение dр_макс /r₀ gH обычно не превышает одной-двух десятых; величина h/Н исчезающе мала, так что изменение напряжения во всем вышележащем массиве и, в частности, на его границах мало сравнительно с исходным напряжением. Поэтому можно считать, что при изменении давления жидкости в пласте напряжения, действующие на кровле и подошве пласта, остаются постоянным.

Предыдущее рассуждение существенно основано на том, что модуль Юнга системы жидкость - пористая среда Е и модуль вышележащего массива горных пород Е₁ имеют одинаковый порядок величины (что обычно имеет место в действительности). Если бы эти модули Юнга сильно отличались между собой, то выражение (21) содержало бы дополнительный множитель Е₁ /Е и при Е₁ >> Е отношение напряжений могло бы и не быть малым. Физически это означает, что в случае, когда вышележащая толща сложена из очень жестких пород, могут образоваться своды, и при изменении давления жидкости напряжения на кровле и подошве пласта будут меняться.

Есть теперь пренебречь влиянием таких границ области фильтрации, как стенки скважин (эти границы имеют сравнительно очень малую протяжность; их влияние будет оценено ниже), то из независимости от времени уравнений равновесия системы жидкость - пористая среда (20) и напряжений на кровле и подошве пласта следует важный вывод о независимости суммарного напряженного состояния в системе жидкость - пористая среда от времени, так что

Свертывая уравнения (22) (т. е. полагая i, j=1, 2, 3 и суммируя получающие уравнения), имеем

Читайте также: