Характеристическая функция, потенциал и функция тока

Представим себе, что имеем плоский фильтрационный поток любой жидкости или газа, подчиняющийся закону Дарси. При рассмотрении одномерных течений было показано, что если фильтрация протекает по закону Дарси, существует потенциальная функция, удовлетворяющая уравнению Лапласа. Но если существует потенциальная функция, то наряду с ней существует функция, также удовлетворяющая уравнению Лапласа. Зная функцию, всегда можно определить функцию путем интегрирования уравнения (8.4).

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

Потенциальную функцию течения и функцию тока всегда можно принять за действительную и мнимую части некоторой функции F (z) комплексного переменного z (8.1).

Функция F (z) называется характеристической функцией течения (комплексным потенциалом).

Исследование любого плоского течения жидкости или газа в пористой среде должно начинаться с определения характеристической функции, соответствующей данной задаче. Найдя ее, мы можем считать задачу решенной. В самом деле, отделив в характеристической функции действительную часть от мнимой, т. е. представив ее в виде, показанном формулой (8.1), можно определить потенциальную функцию (х, у) и функцию тока (х, у). В результате можно представить полную картину потока: принимая различные значения функции, получим уравнения семейства эквипотенциальных линий (х, у) = С, а придавая различные значения, найдем уравнения семейства линий тока (х, у) = С*. По эквипотенциальным линиям определяется распределение давлений в пласте, по линиям тока — направление движения и характер поля скоростей фильтрации.

Проекции вектора массовой скорости фильтрации на оси координат можно записать в виде:

(8.5).

Примечание. Функции тока может быть дан следующий смысл. Фиксируем некоторую линию тока (х, у) = 0 и вообразим канал, ограниченный цилиндрическими поверхностями с образующими, перпендикулярными плоскости течения, проведенными через линию тока = 0 и другую линию тока (х, у) = С* и двумя плоскостями — плоскостью движения и ей параллельной, отстоящей от первой плоскости на расстояние, равное единице (рис. 8.2).

Рис. 8.2. Распределение потока между двумя параллельными плоскостями 1 и 2

При рассмотрении двух произвольных поперечных сечения канала щ1 и щ2 видно, что количество массы жидкости, протекающей через эти сечения в единицу времени (расход) будет одно и то же; внутри такого канала количество массы жидкости при установившемся движении измениться не может; через боковые стенки канала, образованные линиями тока = 0 и (х, у) = С*1, и через плоскости движения жидкость не протекает, следовательно, втекает жидкости в единицу времени через щ1 столько, сколько вытекает через щ2.

Функцией тока можно назвать функцию, принимающую на линии тока (х, у) = С* значение, равное массе жидкости (газа), протекающей в единицу времени через поперечное сечение канала, построенного на линиях = 0 и (х, у) = С*1. Функция тока определена с точностью до произвольной постоянной, зависящей от выбора начальной линии тока = 0.

Массовую скорость фильтрации можно очень просто определить в любой точке пласта, найдя производную от характеристической функции по комплексному аргументу z. Чтобы это показать, составим полный дифференциал от характеристической функции F (z):

(8.6).

Вынося во второй скобке множитель i за знак скобки и воспользовавшись затем уравнениями Коши — Римана (8.4) получим:

т.е.. (8.7).

Учитывая (8.5), перепишем (8.7) в виде.

. (8.8).

Из (8.7) и (8.8) следует, что производная dF/dz есть комплексное число, модуль которого равен модулю массовой скорости фильтрации:

. (8.9).

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

Для однородной несжимаемой жидкости функция тока будет иметь значение объемного (а не массового) расхода жидкости через поперечное сечение канала, построенного на линиях тока =0 и =С*. Модуль же производной от характеристической функции течения будет равен скорости (а не массовой скорости) фильтрации жидкости u.