Анализ функции фильтрационного сопротивления для неустановившегося притока жидкости (газа)Страница 2
Как _ видим, дополнительное слагаемое R(rc , h, f0) в уравнении (6) зависит не только от геометрии пласта, но и от параметра Фурье (f0). В дальнейшем будем называть это слагаемое функцией фильтрационного сопротивления. Заметим, что при h=l (скважина совершенная по степени вскрытия) уравнение (2) представляет собой интегрально-показательную функцию
(7)
С учетом равенства (7) решение (6) запишем в виде
(8)
Разрешая уравнение (8) относительно функции сопротивления и учитывая уравнение (2), находим
(9)
и на основании равенства (7) приведем выражение (9) к виду
(10)
Численное значение R(rс,h,fo) рассчитано по уравнению (10) на ЭВМ в широком диапазоне изменения параметров rc, h, f0. Интеграл (2) вычислялся методом Гаусса, оценка его сходимости выполнена согласно работе [3]. С учетом равенства (7) вычисления дополнительно проконтролированы по значениям интегрально-показательной функции.
С целью выяснения поведения депрессии и функции сопротивления проанализируем их зависимость от значений безразмерных параметров.
1. Определим поведение Dр в зависимости от значений параметров rс, h, f0.
Результаты расчетов значений депрессии для каждого фиксированного rc сведены в таблицы, каждая из которых представляет собой матрицу размером 10х15. Элементы матрицы это значения депрессии Dp(rc) для фиксированных h и f0. Матрица построена таким образом, что каждый ее столбец есть численное значение депрессии в зависимости от h, .а каждая строка соответствует численному значению депрессии в зависимости от fo (табл. 1). Таким образом, осуществлен переход от значений безразмерной депрессии Dp(rc, h, f0) к относительной депрессии
Dр*i,j (rc).
Для удобства построения и иллюстрации графических зависимостей выполнена нормировка матрицы. С этой целью каждый элемент i-й строки матрицы поделен на максимальное значение депрессии в данной строке, что соответствует значению j==15. Тогда элементы новой матрицы определятся выражением
(11)
Условимся элементы матрицы называть значениями относительной депрессии. На рис. 1 приведен график изменения относительной депрессии при фиксированных значениях h. Характер поведения относительной депрессии позволяет описать графики уравнением пучка прямых
Рис. 1. Поведение относительной депрессии (rc=0,0200, hi=const, f0) при значениях h, равных: 1— 0,1; 2 — 0,3; 3—0,5; 4 — 0.7; 5 —0,9; 6—1,0.
где ki — угловой коэффициент прямой, который определяется h и от индекса j не зависит.
Анализ зависимости поведения депрессии Dp*i,j от f0 для всех rc >0,01 показывает, что графики этой зависимости можно описать уравнением пучка прямых для любого значения h. Для rc< 0,01 в графиках зависимости появляются начальные нелинейные участки, переходящие при дальнейшем уменьшении параметра f0 (или же при увеличении его обратной величины 1/foj) в прямые для всех значений h<l,0
(рис. 2). При h=l,0 поведение депрессии строго линейно. Кроме того, протяженность нелинейного участка для разных rc при h=const различна. И чем меньше значение безразмерного радиуса rc , тем больше протяженность нелинейного участка (рис. 2).