УДК 532.54

НЕСТАЦИОНАРНОЕ ТЕЧЕНИЕ ВЯЗКОПЛАСТИЧНОЙ ЖИДКОСТИ

Хусаинова Гузалия Ядкаровна
Стерлитамакский филиал ФГБОУ ВПО «Башкирский государственный университет»
кандидат физико-математических наук, доцент кафедры прикладной информатики и программирования

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

Ключевые слова: вязкопластичная жидкость, нестациональное течение, расперделение давления, температуры


NON-STATIONARY CURRENT OF VISCOUSLY PLASTIC LIQUID

Khusainova Guzaliya Yadkarovna
Sterlitamak Branch of Bashkir State University
associate professor, candidate physical and mathematical sciences

Abstract
In work non-stationary current of viscously plastic liquid is at a plainly-parallel filtration. The analytical decision of distribution of pressure and temperatures is found. The analysis is made.

Рубрика: 01.00.00 ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

Библиографическая ссылка на статью:
Хусаинова Г.Я. Нестационарное течение вязкопластичной жидкости // Современные научные исследования и инновации. 2015. № 10 [Электронный ресурс]. URL: http://web.snauka.ru/issues/2015/10/58442 (дата обращения: 30.09.2017).

Рассматрим нестационарное течение упругой вязкопластичной жидкости. При течении аномальной жидкости в пласте образуется две области фильтрации, на границе которой величина градиента давления должен равняться предельному значению градиента, а давление – пластовому [1-2].
Для однотемпературного приближения уравнение энергии имеет вид: 
(1)
В дальнейшим будем пренебрегать продольной теплопроводностью, а также будем считать, что пористость достаточно мала, так что последнее слагаемое в первой части уравнения теплопроводности (1) не существенна.
Для температуры примем следующие начальноеи граничное   условия.
В начальный момент времени  на границе пласта  начинает работать добывающая галерея, на которой поддерживается постоянное давление . При таком воздействии в пласте образуются две зоны, а именно зона фильтрации жидкости и зона застоя. Граница раздела между этими зонами перемещается со временем по закону причем. В области фильтрации давление удовлетворяет уравнению[1-2]:

, (2)

с граничными условиями 
 при  (3)
и вне области фильтрации давление равно пластовому давлению: при . (4)Необходимо найти распределение давления и температуры в области .
Скорость фильтрации будет определятся из следующего соотношения:

, при (5)
, при .

Решение уравнения (2) найдено методом интегральных соотношений [1]. Распределение давления в области  имеет вид:
, (6)

, (7)

где - перепад давления.
В силу изложенного выше допущений, а также используя (8), уравнение (1) представим в виде:
 (8)
Уравнение (8) решается модифицированным методом характеристик. Уравнение характеристики имеет вид:

, (9)

 определяется из уравнения (3) с помощью интегрального соотношения, характеризующего условие материального баланса [1]: 
. (10)
Решение уравнения (9) имеет вид:
, (11)
где ,

.

На рис.1 приведена динамика движения границы для перепада давления P=106 Па при разных значениях начального градиента давления.
Уравнение для изменения температуры системы, вдоль характеристики в возмущенной области , имеет вид:
 (12)
Отсюда получено выражение, описывающее распределение температуры:

 (13)

Хотя получено аналитическое решение, анализ затруднен, так как зависимость  явно не выражена, поэтому ниже приводится численный расчет[3].
Результаты вычисления с использованием конечно-разностной схемы для плоско-параллельной задачи представлены ниже в виде графиков. На рис. 2 приведены зависимости температуры от времени при разных значениях начального градиента давления. Где кривые 1, 2 соответствуют значениям х=0,4 м; 0,6 м при G0=0,95*105 Па/м; 3,4 – х=0,4 м; 0,6 м при G0=0,5*102 Па/м. На рис.3 иллюстрируется график зависимости температуры от расстояния в разные моменты времени, при разных значениях начального градиента давления и фиксированном перепаде давления (G0=0,95*10Па/м, 1-t=60 мин, 2-t=120 мин; G0=0,01 Па/м, 3-t=60 мин, 4-t=120 мин, P=106 Па). Из графиков следует, что при фильтрации аномальной жидкости вблизи границы температурного эффекта больше, чем при обычной жидкости, но в то же время глубина прогрева в первом случае меньше. Это обстоятельство объясняется тем, что в случае аномальной жидкости при заданном перепаде давления движение прекращается раньше, нежели в случае обычной жидкости.

 
Рисунок 1. Зависимость движения границы от времени (1-G0=0,95*10Па/м, 2-G0=0,085*10Па/м, 3-G0=0,0075*10Па/м, P=106 Па)

Рисунок 2. Зависимость температуры от времени при различных значениях х (1-G0=0,95*10Па/м, х=0,4 м; 2 – G0=0,95*10Па/м, х=0,6 м, 3 – G0=0 Па/м, х=0,4 м, 4 – G0=0 Па/м, х=0,6 м, P=106 Па)

Рисунок 3. Зависимости температуры от расстояния в различные моменты времени (G0=0,95*10Па/м, 1-t=60 мин, 2-t=120 мин; G0=0,01 Па/м, 3-t=60 мин, 4-t=120 мин, P=106 Па).

Библиографический список
  1. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидродинамика. М.: Недра, – 1993. – 416 c.
  2. Алишаев М.Г., Розенберг М.Д., Теслюк Е.В. Неизотермическая фильтрация при разработке нефтяных месторождений. – М.: Недра, 1985. – 271 с.
  3. Хусаинова Г.Я. Исследование температурных полей при фильтрации аномальных жидкостей. Автореферат диссертации на соискание ученой степени кандидата физико-математических наук / Уфа, 1998. – С. 14.


Все статьи автора «Хусаинова Гузалия Ядкаровна»


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

Связь с автором (комментарии/рецензии к статье)

Оставить комментарий

Вы должны авторизоваться, чтобы оставить комментарий.

Если Вы еще не зарегистрированы на сайте, то Вам необходимо зарегистрироваться: