Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

DYNAMICS OF THE RELAXATION OF PRESSURE IN THE CAVITY WITH PLAINLY-PARALLEL WALLS AFTER IT PRESSURE TESTS

Khusainov I.G. 1
1 Sterlitamak branch "Bashkir state university"
In paper dynamics of a relaxation of pressure in a cavity of the form of the crack surrounded by the saturated liquid by the porous and nontight environment, after it pressure tests is researched. Mathematical statement of a problem includes the equation of preservation of weight of a liquid inside of a cavity. The filtration of gas in an environment is absent, it is considered, that gas is in the special container. For definition of speed of a filtration of a liquid from a cavity in surrounding porous environment use Darcy’s law. The field of pressure around of a cavity is described by means of the equation piezoconduction. Compressibility of the liquid which are being a cavity and in the porous environment, is considered in acoustic approximation. Gas is considered calorifically perfect and its behaviour is described by means of the polytropical law. The integrated equation describing a relaxation of pressure in a cavity is received. The integrated equation is solved a numerical method. Dependence of dynamics of restoration of pressure on parameters of the porous environment is researched. It is installed, that dependence of time of a relaxation of pressure on factors of permeability and porosity is back-proportional, and from the initial volumetric contents of gas is a straight line square-law.
acoustic sounding
perforated wellbore
perforating tunnels

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

В [10; 14] метод опрессовки используется для оценки коллекторских характеристик пористой среды, насыщенной газом. Использование волн давления с целью нагрева и очистки призабойной зоны, а также исследование нефтяных скважин рассмотрено в [4; 6-8; 16]. Процесс распространения импульса давления в пористой среде, насыщенной газом, исследуется в [2; 9; 11; 12]. Колебательные процессы, происходящие в неньютоновских жидкостях, исследуются в [1; 5; 15].

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

1. Основные уравнения. Пусть в исходном состоянии (t < 0) давление жидкости во всем пористом пласте вокруг полости постоянно и равно , а сама полость частично заполнена жидкостью и частично газом (рис. 1). В момент времени t = 0 давление в полости мгновенно увеличивается до значения , например, введением некоторого количества газа. Далее, за счет фильтрации жидкости в окружающее пористое пространство, давление в полости будет снижаться до значения .

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

Рис. 1. Схематическое изображение полости, окруженной насыщенной жидкостью пористой средой.

В рамках вышеизложенных допущений уравнение сохранения массы жидкости внутри полости запишем в виде [3; 13]:

, (1)

где – средняя плотность жидкости внутри полости, которая определяется по формуле ; – плотность жидкости; – объемная доля газа в полости; – полутолщина полости с плоскопараллельными стенками; – скорость фильтрации жидкости через стенки полости.

Для определения скорости фильтрации жидкости из полости в окружающую пористую среду используем закон Дарси:

, (2)

где – динамический коэффициент вязкости жидкости; – коэффициент проницаемости; – давление и скорость фильтрации жидкости вокруг полости.

Поле давления вокруг полости описывается с помощью уравнения пьезопроводности:

. (3)

Здесь – коэффициент пьезопроводности; – коэффициент пористости; – скорость звука в жидкости.

Сжимаемость жидкости, находящейся в полости и в пористой среде, будем учитывать в акустическом приближении:

. (4)

Газ будем считать калорически совершенным, тогда для его поведения примем политропический закон:

, (5)

где – показатель политропы, – начальное давление в полости. Здесь и далее индекс (0) внизу соответствует начальному значению величины.

Начальное условие для уравнения (3) запишем в виде:

. (6)

Граничные условия на стенке полости (r = a) для уравнений (2), (3) могут быть записаны в виде:

. (7)

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

, (8)

где .

2. Результаты численных расчетов. Уравнение (8) решается численным методом. В работе для изучения динамики релаксации давления в полости предлагается использовать период полувосстановления давления. Периодом полувосстановления давления будем называть промежуток времени, в течение которого разница между значениями давлений в полости и пористой среде снижается в два раза от начальной разницы.

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

Рис. 2. Эволюция давления в полости с плоскопараллельными стенками при различных значениях коэффициента пористости : 1 – ; 2 –.

На рис. 3 представлены зависимости периода полувосстановления давления от начального объемного содержания газа при различных значениях коэффициента проницаемости . Анализ графиков показывает, что зависимость периода полувосстановления давления от начального объемного содержания газа квадратичная, т.е. .

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

Рис. 3. Зависимости периода полувосстановления давления от начальной объемной доли газа в полости при различных значениях коэффициента проницаемости : 1 – , 2 – .

Выводы. Разработана математическая модель и получено нелинейное интегральное уравнение, описывающее процесс релаксации давления в полости формы трещины, окруженной насыщенной жидкостью пористой средой и опрессованной введением газа. В результате численного анализа интегрального уравнения установлено: зависимость периода полувосстановления от коэффициентов проницаемости и пористости является обратно пропорциональной, а от начального объемного содержания газа является прямой квадратичной.

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

Работа выполнена при поддержке гранта СФ БашГУ № В14-19.

Рецензенты:

Мустафина С.А., д.ф.-м.н., профессор кафедры математического моделирования, Стерлитамакский филиал ФГБОУ ВПО «Башкирский государственный университет», Республика Башкортостан, г. Стерлитамак;

Михайлов П.Н., д.ф.-м.н., профессор кафедры алгебры, геометрии и методики обучения математике, Стерлитамакский филиал ФГБОУ ВПО «Башкирский государственный университет», Республика Башкортостан, г. Стерлитамак.