.

Фильтрация газов(баротермический эффект)

Язык: русский
Формат: дипломна
Тип документа: Word Doc
0 1119
Скачать документ

Фильтрация газов(баротермический эффект)

Язык: русский
Формат: реферат
Тип документа: Word Doc
0 1038
Скачать документ

министерство общего и профессионального образования российской федерации

стерлитамакский государственный педагогический

институт

Кафедра теоретической физики

МАРИО1980mail.ru

Исследование влияния сжимаемости

на величину баротермического эффекта при фильтрации газа

Дипломная работа

Научный руководитель:

д.т.н., проф. Филиппов А.И.

ст. пр. Миколайчук Н.П.

Стерлитамак 2002

содержание

TOC \o “1-5” введение 4

глава 1. постановка задачи о температурном поле в нефтегазовом пласте 8

1.1. Уравнения состояния реального газа 8

1.2. Основные уравнения, описывающие процесс фильтрации газа в пористой
среде………………………………..…………………….11

1.3. Описание задачи 13

1.4. Математическая постановка задачи 14

1.4.1. Математическая постановка температурной задачи 14

1.4.2. Математическая постановка гидродинамической задачи 15

1.5. Основные идеи метода характеристик………………………….…15

1.6. Выводы………………………………………………………..……..22

Глава 2. Аналитическое решение задачи о баротерми-ческом эффекте для
реальных уравнений состояния 23

2.1. Решение гидродинамической задачи 23

2.2. Решение температурной задачи 25

2.3. Выводы 27

Глава 3. Получение Аналитических выражений решения задачи о
баротермическом эффекте с учетом барической сжимаемости 27

3.1. Решение гидродинамической задачи для линеаризованного уравнения
состояния 27

3.2. Температурная задача в линеаризованном случае 28

3.3. Выводы 30

Глава 4. анализ результатов расчетов и Исследование температурных полей,
возникающих при фильтрации газа 30

4.1. Анализ результатов расчетов температурных полей 31

4.2. Изучение вклада сжимаемости в величину баротермического эффекта 38

4.3. Выводы 40

заключение 41

СПИСОК ЛИТЕРАТУРЫ 42

ПРИЛОЖЕНИЯ 43

введение

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

Для решения практических задач необходимо знать зависимость температуры
от расстояния; температуры от времени при различных параметрах пластов.

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

Задачами исследования являются

– разработка математической модели термодинамических эффектов в
прискважинной зоне газовых пластов;

– постановка задачи о баротермическом эффекте в прискважинной зоне,
построение аналитического решения;

– проведение расчетов и анализ вклада различных физических процессов в
температурное поле в прискважинной зоне;

– изучение влияния сжимаемости на величину баротермического эффекта.

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

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

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

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

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

Краткая характеристика содержания работы: Дипломная работа состоит из
введения, трех глав и заключения.

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

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

Во второй главе найдено решение гидродинамической задачи методом
разделения переменных, методом характеристик построено решение
температурной задачи и осуществлен анализ полученного аналитического
решения на частных случаях.

В третьей главе осуществлены численные расчеты тепловых полей с помощью
программного пакета Mathcad. Описан анализ вклада различных физических
процессов.

В заключении подводятся итоги проведенного исследования.

При выполнении работы оказали большую помощь д.т.н., проф. Филиппов
А.И., ст. пр. Миколайчук Н.П., ст. лаб. Скворцова О.В. В связи с этим,
хочу выразить им большую благодарность за оказанную помощь в выполнении
дипломной работы, указание путей решения возникающих трудностей, советы
по рациональной организации труда.

СПИСОК ОБОЗНАЧЕНИЙ

;

;

;

;

;

?=с ?/cpl;

m – пористость;

;

;

;

;

;

;

;

;

;

;

глава 1. постановка задачи о температурном поле в нефтегазовом пласте

1.1. Уравнения состояния реального газа

Модель идеального газа хорошо описывает свойства газообразного состояния
вещества при средних и высоких температурах (от комнатной и выше) и
небольших давлениях (около атмосферного). Расчет свойств газов в широком
интервале экспериментальных условий требует использования уравнения
состояния реального газа[1].

Реальным газом называется газ, между молекулами которого существуют
заметные силы межмолекулярного взаимодействия. Оно имеет
электромагнитную и квантовую природу и осуществляется посредством сил
межмолекулярного притяжения и отталкивания.

Силы притяжения, проявляющиеся на расстояниях r между центрами молекул
порядка 10-7 см, называются ван-дер-ваальсовыми силами. Они убывают с
расстоянием ( r –7, что соответствует изменению потенциальной энергии по
закону r –6.

Различают три вида ван-дер-ваальсовых сил [7]:

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

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

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

Полная потенциальная энергия ван-дер-ваальсовых сил описывается суммой:

U = Uор + Uинд + Uдисп. (I.1.1)

Для полярных молекул основную роль играют ориентационные силы
притяжения, для остальных молекул – дисперсионные силы. Энергия
ван-дер-ваальсового притяжения составляет (0,1 – 1) ккал/моль [7]. В
большинстве случаев ван-дер-ваальсовы силы притяжения перекрываются
значительно превосходящими их химическими валентными силами притяжения с
энергиями порядка (10 – 100) ккал/моль.

Согласно упрощенной модели ван-дер-ваальсовых сил, молекулы газа –
абсолютно упругие шары – притягиваются с силами, достигающими
наибольшего значения при непосредственном их соприкосновении. Силы
отталкивания проявляют себя на значительно меньших расстояниях.

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

Уравнение Ван-дер-Ваальса состояния реального газа имеет следующий вид
[7]:

– внутреннее давление, обусловленное силами притяжения между
молекулами, b – поправка за собственный объем молекул, учитывающая
действие сил отталкивания между молекулами и равная учетверенному объему
молекул в 1 моле газа:

, (I.1.3)

. (I.1.4)

Здесь NA – число Авогадро, d – диаметр молекулы, U(r) – потенциальная
энергия притяжения двух молекул.

Уравнение состояния Бертло (1900г.):

. (I.1.5)

Здесь а и b связаны с параметрами критического состояния (в критической
точке) соотношениями [8]:

. (I.1.6)

Уравнение состояния Вукаловича и Новикова [7]:

. (I.1.7)

Здесь B1, B2 и т.д. – так называемые вириальные коэффициенты весьма
сложного вида. Их вычисление производится с учетом ассоциации молекул –
объединения под влиянием ван-дер-ваальсовых сил притяжения.

Уравнение состояния Майера [7]:

d(i=dqi1*…dqin.

Здесь Uпij – взаимная потенциальная энергия i-й и j-й молекул,
взаимодействующих по закону центральных сил, qi1,……,qin – обобщенные
координаты i-той молекулы, обладающей n степенями свободы.

Уравнение Камерлинг-Оннеса (1901) [8]:

.

Уравнение Редлиха-Квонга (1949 г.) [8]:

0,42748·R2·T2,5k/Pk, b = 0,08664·R·Tk/Pk. Уравнение Редлиха-Квонга
считается наилучшим двухконстантным уравнением. При его выводе авторы не
руководствовались какими-то определенными теоретическими обоснованиями
[8]. Это уравнение следует рассматривать как произвольную, но удачную
эмпирическую модификацию предшествующих уравнений состояния.

Уравнение Мартина (1967 г.) [8]:

27·R2·T2k/(64Pk), b = R·Tk/(8Pk).

1.2. Основные уравнения, описывающие процесс фильтрации газа в пористой
среде

В последнее время наблюдается рост интереса к различным
термодинамическим эффектам в пористых средах. Это связано с их
многообразными практическими приложениями[4,5].

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

с учетом термодинамических эффектов высокого порядка [9]

(I.2.1)

где первое слагаемое в левой части уравнения (I.2.1) описывает изменение
температуры в пласте со временем, второе – за счет конвекции
(перемещения больших объемов газа). Первое слагаемое в правой части
ответственно за теплопроводность, второе – за межфракционный теплообмен,
третье описывает адиабатический эффект, четвертое – эффект
Джоуля-Томсона и пятое – влияние поля тяготения Земли.

Вторым уравнением системы является уравнение неразрывности, которое
записывается в виде:

. (I.2.2)

Фильтрация газа подчиняется закону Дарси

. (I.2.3)

К системе добавляется уравнение состояния

. (I.2.4)

Система (I.2.1)-(I.2.4) является нелинейной, кроме того, уравнения
(I.2.1)-(I.2.2) являются взаимосвязанными.

1.3. Описание задачи

Рис. 1. постановка задачи

При описании температурной задачи примем следующие допущения:

пористый пласт считается однородным и изотропным по гидродинамическим и
теплофизическим свойствам;

давления в скважине и на контуре питания остаются неизменными;

породы, окружающие пласт предполагаются непроницаемыми и однородными по
своим теплофизическим свойствам;

температуры газа и скелета пористой среды в каждой точке совпадают;

естественное тепловое поле Земли считается стационарным;

пласт расположен на глубине порядка 1 – 2 км, поэтому суточные и
сезонные колебания температуры не достигают пласта;

адиабатическим эффектом, обусловленным гравитационным полем
пренебрегаем.

1.4. Математическая постановка задачи

Математическая постановка задачи включает температурную задачу,
гидродинамическую задачу, уравнение состояния и соотношение для поля
скорости конвективного переноса тепла. Ниже рассматриваются
соответствующие постановки задач.

1.4.1. Математическая постановка температурной задачи

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

. (I.4.1.1)

Будем рассматривать задачу при следующих условиях температуры:

начальном

, (I.4.1.2)

и граничном

. (I.4.1.3)

1.4.2. Математическая постановка гидродинамической задачи

Математическая постановка гидродинамической задачи в полярной системе
координат примет следующий вид. Учитывая, что для осесимметричного
течения поле давления является функцией координаты r уравнение можно
представить в виде:

, (1.4.2.1)

Будем рассматривать задачу при следующих условиях. Пусть PC – давление
на границе контура питания. При значении радиуса, равном радиусу контура
питания

, (1.4.2.2)

давление поддерживается равным Рс:

, (1.4.2.3)

Pс – давление на контуре питания.

При значении радиуса, равном радиусу скважины

, (1.4.1.3)

давление поддерживается равным PW:

, (1.4.1.4)

где PW – давление в скважине.

1.4. Основные идеи метода характеристик[6]

В данном разделе рассмотрим метод характеристик. Любое линейное
дифференциальное уравнение второго порядка (при двух независимых
переменных) может быть записано в следующем виде:

(1.4.1)

где а, b, с, d, e, f, g — заданные непрерывные функции от x и y (или в
частном случае, постоянные).

Попытаемся упростить это уравнение с помощью замены независимых
переменных:

) нигде в области G не обращался в нуль.

Для того чтобы сделать требуемую замену переменных, выразим частные
производные от функции u по х и у через производные от и по ( и (:

(1.4.31)

, через производные по ( и (, учтем формулу (1.4.31) и применим снова
правило дифференцирования сложной функции:

Следовательно,

(1.4.41)

Аналогично найдем:

(1.4.42)

Подставляя u’x, u’y, u’xx,… из этих формул в уравнение (1), мы
получим снова линейное уравнение второго порядка с неизвестной функцией
и и независимыми переменными( и (:

(1.4.5)

где

— функция, линейная относительно и’( , u’( , u .

Уравнение (1.4.5) становится особенно простым, если в нем коэффициенты а
и с окажутся равными нулю. Для того чтобы первоначально заданное
уравнение (1.4.1) можно было привести к такому простому виду, надо в нем
сделать замену переменных

подобрав функции ( и ( так, чтобы они являлись решениями уравнения:

(1.4.6)

Это уравнение является нелинейным уравнением в частных производных
первого порядка. Следующая теорема покажет, как связаны его решения с
общим решением некоторого обыкновенного уравнения.

Теорема. Для того чтобы функция z = f(x, у) во всех точках области G
удовлетворяла уравнению (6), необходимо и достаточно, чтобы, семейство

(1.4.7)

было общим интегралом уравнения

(1.4.8)

в той же области G.

Доказательство. Необходимость. Пусть z = f(x, у)— решение уравнения
(1.4.6). Рассмотрим семейство кривых f(x, у) — k и докажем, что любая
кривая этого семейства удовлетворяет уравнению (1.4.7).

В любой точке, лежащей на кривой f(x, у) = k (где k — фиксировано),
выполняется следующее равенство:

действительно вдоль данной кривой функция f(x, у) постоянна, и поэтому
ее полный дифференциал равен нулю.

Следовательно, всюду на кривой имеет место равенство:

обозначим каждое из этих отношений через (; тогда

Подставляя эти выражения для dx и dy в левую часть уравнения (1.4.8),
получим:

Выражение, стоящее в квадратных скобках, равно нулю, так как, по
условию, функция f(x, у) есть решение уравнения (1.4.6). Следовательно,
во всех точках нашей кривой имеет место равенство

откуда вытекает, что она является интегральной кривой уравнения
(1.4.8).

Итак, любая кривая вида f(x, у) = k является интегральной кривой
уравнения (1.4.8); с другой стороны, через каждую точку области G
проходит кривая такого вида; это вытекает из того, что функция f(x, у)
определена всюду в области G и поэтому, например, через точку (х0, у0)
проходит кривая f(x,y)=f(x0,y0).

Отсюда следует, что семейство f(x, у) = k является общим интегралом
уравнения (1.4.8).

Достаточность. Пусть семейство f(х, у)= k будет общим интегралом
уравнения (1.4.8). Возьмем произвольную точку (х0, у0) из G и выделим ту
кривую семейства, которая проходит через эту точку:

f(x, у) = k0.

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

откуда

(1.4.10)

Так как кривая является интегральной кривой уравнения (1.4.8), то при
подстановке в это уравнение dx и dy из (1.4.10), получим тождество:

или, после сокращения на (2:

В частности, в точке (х0, у0) имеет место:

Но последнее равенство означает, что функция двух переменных f(x, у)
удовлетворяет в точке (х0, у0) уравнению (1.4.7). Так как точка (х0, y0)
была взята произвольно в области G, то функция f(x, у) удовлетворяет
уравнению (1.4.7) во всех точках этой области, т. е. эта функция
является одним из решений уравнения (1.4.7).

Таким образом, теорема доказана.

Рассмотренная теорема открывает путь для упрощения исходного уравнения
(1.4.1). Для этого сначала составляем вспомогательное уравнение (1.4.8);
оно называется характеристическим уравнением для данного уравнения
(1.4.1). Характеристическое уравнение есть обыкновенное дифференциальное
уравнение первого порядка, но второй степени. Разрешая его относительно
y’x (предварительно разделив все члены уравнения на dx2), получим два
уравнения:

(1.4.101)

(1.4.102)

(предполагается, что ас — b2<0, b2 —ас>0 всюду в области G). Пусть общий
интеграл уравнения (1.4.101) имеет вид

((х, у)= k , (1.4.111)

а общий интеграл уравнения (1.4.102)

((х, у)= k. (1.4.112)

Интегральные кривые характеристического уравнения (т. е. все кривые,
входящие в семейства (1.4.111) и (1.4.112)) называются характеристиками
заданного дифференциального уравнения (1.4.1). В связи с этим
рассматриваемый метод упрощения уравнения (1.4.1) называется методом
характеристик.

Семейства (1.4.111) и (1.4.112) можно рассматривать, как общие интегралы
уравнения (1.4.8) (это уравнение распадается на два уравнения (1.4.101)
и (1.4.102)).

Следовательно, согласно доказанной теореме, функции

z=((х, у) и z=((х, у)

являются решениями уравнения в частных производных (1.4.6).

Функции ((х, у) и ((х, у) независимы друг от друга (можно доказать,
что их якобиан отличен от нуля, если ас- b2<0). Поэтому, возвращаясь к уравнению (1.4.1), мы можем в нем сделать замену переменных:. Следовательно, уравнение (1.4.1) преобразуется к виду:или, после деления на 2b и переноса в другую часть равенства:– функция, линейная относительно и’( , u’( , u (см. выше, формула (1.4.5)).Полученное уравнение имеет более простой вид, чем исходное уравнение (1.4.1); если мы его сможем решить (т. е. найти и как функцию от ( и (), то для того, чтобы найти решение исходного уравнения, достаточно вернуться к старым переменным (выразив ( и ( через х и у).1.5. ВыводыВ данной главе представлены основные уравнения состояния реального газа, уравнения, описывающие процесс фильтрации газа в пористой среде. Дано описание задачи. Сформулированы физическая и математическая постановки температурной и гидродинамической задач. Описан метод характеристик, который использован для получения решения задачи в главе 2.Глава 2. Аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состоянияВ данной главе приведены аналитические решения гидродинамической и температурной задач для произвольного уравнения состояния.2.1. Решение гидродинамической задачиРешение уравнения (I.4.1.1) приводит к необходимости решения дополнительной гидродинамической задачи для отыскания поля давления. Для описания движения газа воспользуемся квазистационарным уравнением неразрывности:. (2.1.1)Скорость фильтрации газа сквозь пористую среду определяется законом Дарси:не меняются при движении газа к скважине, и что плотность газа зависит только от давления (баротропное приближение), перепишем уравнение неразрывности в следующем виде:(2.1.3)Функцию Лейбензона представим в виде:и А задаются граничными условиями. Подставляя функцию Лейбензона в уравнение (2.1.3), получим:, уравнение (2.1.5) можно представить в виде:(2.1.6)Решение этого уравнения представим в виде:), тогда согласно (2.1.4) и (2.1.7):(2.1.10)(2.1.11)Подставив эти выражения в уравнение (2.1.7), получим:Используя (2.1.12), получим выражение для градиента давления в виде:(2.1.13)Найденное выражение для градиента давления позволяет преобразовать температурную задачу к виду удобному для аналитического решения. Решение температурной задачи рассматривается в следующем параграфе.2.2. Решение температурной задачиС учетом (2.1.13) и закона фильтрации Дарси (2.1.2) задача (I.4.1.1) – (I.4.1.3) преобразуется к виду:(2.2.1)Условия (I.4.1.1) – (I.4.1.3) остаются неизменными. Вводя обозначение(2.2.2)представим уравнение Чекалюка в виде:(2.2.3)Будем рассматривать задачу при следующих условиях для температуры:начальном(2.2.4)и граничном::(2.2.7)определяет область применимости нестационарного решения(2.2.8)Уравнение (2.2.3) с учетом (2.2.6) можно представить в виде:(2.2.9)откуда(2.2.10)где(2.2.11)Исключив константу С из (2.2.10) и (2.2.11), окончательно получаем нестационарное решение:. Для этой области(2.2.13)и стационарное решение, удовлетворяющее условию (2.2.5), имеет вид:(2.2.14)Выражения (2.2.12) и (2.2.14) полностью решают поставленную задачу для любого уравнения состояния2.3. ВыводыВ данной главе представлено аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состояния, которая включает в себя температурную и гидродинамическую задачи.Глава 3. Получение Аналитических выражений решения задачи о баротермическом эффекте с учетом барической сжимаемости3.1. Решение гидродинамической задачи для линеаризованного уравнения состоянияВыпишем полученные решения для линеаризованного баротропного уравнения состояния(3.1.1)Вычислив интеграл, входящий в (2.2.3):(3.1.2)Представим зависимость между давлением и радиальной координатой r в виде:(3.1.3)Введем обозначение(3.1.4)Тогда уравнение (3.1.3) преобразуется к виду:(3.1.5)Откуда найдем(3.1.6)Физический смысл имеет только значение полученного выражения со знаком плюс перед квадратным корнем. Введем обозначения(3.1.7)и упростить запись выражения (3.1.6)(3.1.9)Подставив (3.1.9) в (3.1.1), получим зависимость плотности от радиальной координаты r:(3.1.10)Полученные в данном разделе выражения позволяют построить решения задачи о баротермическом эффекте в случае линеаризованного уравнения состояния.3.2. Температурная задача в линеаризованном случаеВ этом случае нестационарное решение для температуры (3.1.5) записывается в виде:(3.2.1)Интеграл в (3.2.1) легко вычисляется; окончательно нестационарное решение представляется в виде(3.2.2)Выражения для G и H представляются формулами (3.2.5) и (3.2.6), а - для V представляется в виде, следующем из(2.2.8)(3.2.3)В пределе при ??0 из (3.2.2)-(3.2.3) следует известное решение для несжимаемой жидкости[4]:(3.2.4)Аналогично в стационарном случае из (2.2.14) получим:(3.2.5)В пределе при ??0 из (3.2.5) и (3.2.3) следует известное решение для несжимаемой жидкости[4]:(3.2.6)Выражения (3.2.2), (3.2.4) решают поставленную задачу о баротермическом эффекте при фильтрации газа в прискважинной зоне реальных газовых пластов. Такое решение поставленной задачи получено впервые. Поэтому представляет значительный и практический интерес анализ результатов расчетов на основе полученных решений, что и приведено в четвертой главе.3.3. ВыводыВ данной главе получено аналитическое решение задачи о баротермическом эффекте с учетом барической сжимаемости, которая включает в себя решение гидродинамической задачи для линеаризованного уравнения состояния и температурную задачу в линеаризованном случае.Глава 4. анализ результатов расчетов и Исследование температурных полей, возникающих при фильтрации газаВ данной главе приведен анализ результатов расчетов баротермического эффекта в прискважинной зоне газовых пластов применительно к реальным месторождениям газа.4.1. Анализ результатов расчетов температурных полей.Из рисунка видно, что изменение температуры подчиняется следующим закономерностям. Линейное нарастание температурного эффекта при малых временах сменяется логарифмической стабилизацией - при больших временах. Время, при котором происходит смена линейного нарастания на логарифмическую стабилизацию, зависит от барической сжимаемости; с увеличением сжимаемости это время уменьшается.Величина температурного эффекта также сильно зависит от сжимаемости. С увеличением сжимаемости величина температурного эффекта возрастает. Коэффициент барической сжимаемости приблизительно обратно пропорционален давлению. Реальные значения этого коэффициента в условиях газовых пластов лежат в пределах от 3 10-8 Па-1 до 10-5; поэтому величина эффекта лежит в пределах до –10 ( –15 К.. Это хорошо согласуется с величиной измеряемых в скважинных условиях температурных эффектов.Рис.1. Зависимость величины баротермического эффекта от времени при различных барических сжимаемостях. Обозначения: 1 - ( = 3 10-4 Па-1, 2 – 10-5, 3 – 10-6, 4 – 10-7, 5 – 5 10-8Важно отметить, что согласно разработанной нами теории время установления температурного эффекта при ( ( 10-8 Па-1, что часто встречается на практике, составляет около суток. Этот факт чрезвычайно важен при практическом использовании баротермического эффекта....Рис 4. Зависимость нестационарной температуры от времени при различных радиусах скважины. Обозначения: 1- r =0.1 м; 2 -0.2 ; 3 – 0.3; 4 -0.5..Рис 5. Зависимость нестационарной температуры от радиуса скважины при различных временах. Обозначения: 1- t =10 000 с; 2 -100 000 ; 3 – 1 000 000...Из рисунка видно, что величина температурного эффекта возрастает при увеличении теплоемкости.Рис 6. Зависимость нестационарной температуры от времени при различных радиусах контура питания. Обозначения: 1- R =25 м; 2 -50; 3 – 100; 4 -200; 5 - 250.Рис 7. Зависимость нестационарной температуры от теплоёмкости при различных временах. Обозначения: 1- t =100 000 c; 2 -1 000 000; 3 – 10 000 000..Рис 8. Зависимость нестационарной температуры от относительной вязкости при различных временах. Обозначения: 1- t =100 000 c; 2 -1 000 000; 3 – 1 500 000.На рис. 9. показана зависимость баротермического эффекта от времени при различных коэффициентах барической сжимаемости. Из рисунка видно,.4.2. Изучение вклада сжимаемости в величину баротермического эффекта.Рис. 10. Зависимость нестационарной температуры от коэффициента барической сжимаемости при различных временах. Обозначения: 1- t = 100 c; 2 -1000 ; 3 – 10 000; 4 -100 000.На рис. 11. показана зависимость стационарной температуры от коэффициента барической сжимаемости. Из рисунка видно что величина температурного эффекта в стационарном случае не зависит от коэффициента барической сжимаемости. В расчетах принято: ?=-0.5?10-5; rW=0.1; с=850; k=10-15; сPL=84000000; µ=10-5; R=100; ?=150; P=100?105; P0=150?105; PC=200?105; PW=150?105.Рис. 11 Зависимость стационарной температуры от коэффициента барической сжимаемости.На рис. 12. приведена зависимость времени установления температуры от коэффициента барической сжимаемости.Рис 12. Зависимость времени установления температуры от коэффициента барической сжимаемости.Итак, изучение вклада сжимаемости в величину баротермического эффекта показывает, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.4. 3. ВыводыВ данной главе сделан анализ результатов расчетов и исследованы температурные поля, возникающих при фильтрации газа. Показано, что величина температурного эффекта составляет около 20 К. Время установления температурного эффекта сильно зависит от проницаемости и для реальных значений проницаемости составляет приблизительно сутки. Это важно учитывать при интерпритации результатов термических исследований скважин.Изучен вклад сжимаемости в величину баротермического эффекта. Показано, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.Показано, что время установления баротермического эффекта зависит от барической сжимаемости и лежит в пределах до 109 с при ?~10-8 Па-1. При ?~10-8 Па-1 время полного установления составляет (приблизительно) три года. Значит температурные поля в газовом пласте практически всегда нестационарны.ЗаключениеВ ходе проделанной работы были получены следующие результаты:Описаны основные уравнения состояния реального газа, уравнения, описывающие процесс фильтрации газа в пористой среде.Представлено аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состояния.Получено аналитическое решение задачи о баротермическом эффекте с учетом барической сжимаемости.Сделан анализ результатов расчетов и исследование температурных полей, возникающих при фильтрации газа.Исследованы температурные поля и изучен вклад сжимаемости в величину баротермического эффекта. Показано, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.При ?~10-8 Па-1 время полного установления температуры составляет (приблизительно) три года. Это означает, что температурные поля в газовом пласте практически всегда нестационарные. Следует отметить при этом что логарифмическая стабилизация достигается при времени около суток.Список использованной литературыЛандау Л. Д., Лившиц Е. М. Статистическая физика// М.,1964.Карслоу Г., Егер Д. Теплопроводность твердых тел// М: Наука. 1964. 487с.Варгафтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей// М., Наука. 1972.Филиппов А. И., Фридман А. А., Девяткин Е. М. Баротермический эффект при фильтрации газированной жидкости: Монография. - Стерлитамак: Стерлитамак. гос. пед. ин-т; Стерлитамакский филиал Академии наук Республики Башкортостан, 2000. – 175с.Филиппов А. И. Скважинная термометрия переходных процессов. - Саратов: Изд-во Сарат. ун-та, 1989. – 116с.Очан Ю. С. Методы математической физики// М: Высшая школа. 1965. 383с.Яворский Б. М., Детлаф А. А. Справочник по физике. – М.: Наука, 1971. – 940с.Морачевский А. Г., Сладков И. Б. Физико – химические свойства молекулярных неорганических соединений. – С. Пб.: Химия, 1996. – 312с.Баскаков А. П., Гуревич М. И., Решетин Н. И. и др. Общая теплотехника. – М.-Л.: Государственное энергетическое издательство, 1963. – 392с.PAGE 1PAGE 42-25.15?10-4?-24.9-25.0?T,K431210-38?10-46?10-41234?10-42?10-48?1046?10444?1042?104-20-16-12-8-4t,c001?1055?T,K0-20-15-10-5??T,K0t,c?T,K7?105?-4-3-2-10-5-60-741231?1055?1056?1052?1053?1054?10500?T,K1234t,c-2-10-8-6-42?1054?1056?1058?1051?1060?T, Kt-4-24321-60-85?1051?1060r123-5-10-15-200?T, K0.511.50?T, K-2-4-6-8-101325?1050541?106t, c0-5321с10005000-10?T, K0-5321-10-15-20-2502?10-64?10-66?10-68?10-61?10-7?T, Kµt, c-1510-30-25-2023451·10102·10103·10104·1010?T, K?, с6?1084?1082?1081?10-31?10-85?10-4?0?

Нашли опечатку? Выделите и нажмите CTRL+Enter

министерство общего и профессионального образования российской федерации

стерлитамакский государственный педагогический

институт

Кафедра теоретической физики

МАРИО1980mail.ru

Исследование влияния сжимаемости

на величину баротермического эффекта при фильтрации газа

Дипломная работа

Научный руководитель:

д.т.н., проф. Филиппов А.И.

ст. пр. Миколайчук Н.П.

Стерлитамак 2002

содержание

TOC \o “1-5” введение 4

глава 1. постановка задачи о температурном поле в нефтегазовом пласте 8

1.1. Уравнения состояния реального газа 8

1.2. Основные уравнения, описывающие процесс фильтрации газа в пористой
среде………………………………..…………………….11

1.3. Описание задачи 13

1.4. Математическая постановка задачи 14

1.4.1. Математическая постановка температурной задачи 14

1.4.2. Математическая постановка гидродинамической задачи 15

1.5. Основные идеи метода характеристик………………………….…15

1.6. Выводы………………………………………………………..……..22

Глава 2. Аналитическое решение задачи о баротерми-ческом эффекте для
реальных уравнений состояния 23

2.1. Решение гидродинамической задачи 23

2.2. Решение температурной задачи 25

2.3. Выводы 27

Глава 3. Получение Аналитических выражений решения задачи о
баротермическом эффекте с учетом барической сжимаемости 27

3.1. Решение гидродинамической задачи для линеаризованного уравнения
состояния 27

3.2. Температурная задача в линеаризованном случае 28

3.3. Выводы 30

Глава 4. анализ результатов расчетов и Исследование температурных полей,
возникающих при фильтрации газа 30

4.1. Анализ результатов расчетов температурных полей 31

4.2. Изучение вклада сжимаемости в величину баротермического эффекта 38

4.3. Выводы 40

заключение 41

СПИСОК ЛИТЕРАТУРЫ 42

ПРИЛОЖЕНИЯ 43

введение

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

Для решения практических задач необходимо знать зависимость температуры
от расстояния; температуры от времени при различных параметрах пластов.

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

Задачами исследования являются

– разработка математической модели термодинамических эффектов в
прискважинной зоне газовых пластов;

– постановка задачи о баротермическом эффекте в прискважинной зоне,
построение аналитического решения;

– проведение расчетов и анализ вклада различных физических процессов в
температурное поле в прискважинной зоне;

– изучение влияния сжимаемости на величину баротермического эффекта.

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

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

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

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

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

Краткая характеристика содержания работы: Дипломная работа состоит из
введения, трех глав и заключения.

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

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

Во второй главе найдено решение гидродинамической задачи методом
разделения переменных, методом характеристик построено решение
температурной задачи и осуществлен анализ полученного аналитического
решения на частных случаях.

В третьей главе осуществлены численные расчеты тепловых полей с помощью
программного пакета Mathcad. Описан анализ вклада различных физических
процессов.

В заключении подводятся итоги проведенного исследования.

При выполнении работы оказали большую помощь д.т.н., проф. Филиппов
А.И., ст. пр. Миколайчук Н.П., ст. лаб. Скворцова О.В. В связи с этим,
хочу выразить им большую благодарность за оказанную помощь в выполнении
дипломной работы, указание путей решения возникающих трудностей, советы
по рациональной организации труда.

СПИСОК ОБОЗНАЧЕНИЙ

;

;

;

;

;

?=с ?/cpl;

m – пористость;

;

;

;

;

;

;

;

;

;

;

глава 1. постановка задачи о температурном поле в нефтегазовом пласте

1.1. Уравнения состояния реального газа

Модель идеального газа хорошо описывает свойства газообразного состояния
вещества при средних и высоких температурах (от комнатной и выше) и
небольших давлениях (около атмосферного). Расчет свойств газов в широком
интервале экспериментальных условий требует использования уравнения
состояния реального газа[1].

Реальным газом называется газ, между молекулами которого существуют
заметные силы межмолекулярного взаимодействия. Оно имеет
электромагнитную и квантовую природу и осуществляется посредством сил
межмолекулярного притяжения и отталкивания.

Силы притяжения, проявляющиеся на расстояниях r между центрами молекул
порядка 10-7 см, называются ван-дер-ваальсовыми силами. Они убывают с
расстоянием ( r –7, что соответствует изменению потенциальной энергии по
закону r –6.

Различают три вида ван-дер-ваальсовых сил [7]:

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

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

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

Полная потенциальная энергия ван-дер-ваальсовых сил описывается суммой:

U = Uор + Uинд + Uдисп. (I.1.1)

Для полярных молекул основную роль играют ориентационные силы
притяжения, для остальных молекул – дисперсионные силы. Энергия
ван-дер-ваальсового притяжения составляет (0,1 – 1) ккал/моль [7]. В
большинстве случаев ван-дер-ваальсовы силы притяжения перекрываются
значительно превосходящими их химическими валентными силами притяжения с
энергиями порядка (10 – 100) ккал/моль.

Согласно упрощенной модели ван-дер-ваальсовых сил, молекулы газа –
абсолютно упругие шары – притягиваются с силами, достигающими
наибольшего значения при непосредственном их соприкосновении. Силы
отталкивания проявляют себя на значительно меньших расстояниях.

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

Уравнение Ван-дер-Ваальса состояния реального газа имеет следующий вид
[7]:

– внутреннее давление, обусловленное силами притяжения между
молекулами, b – поправка за собственный объем молекул, учитывающая
действие сил отталкивания между молекулами и равная учетверенному объему
молекул в 1 моле газа:

, (I.1.3)

. (I.1.4)

Здесь NA – число Авогадро, d – диаметр молекулы, U(r) – потенциальная
энергия притяжения двух молекул.

Уравнение состояния Бертло (1900г.):

. (I.1.5)

Здесь а и b связаны с параметрами критического состояния (в критической
точке) соотношениями [8]:

. (I.1.6)

Уравнение состояния Вукаловича и Новикова [7]:

. (I.1.7)

Здесь B1, B2 и т.д. – так называемые вириальные коэффициенты весьма
сложного вида. Их вычисление производится с учетом ассоциации молекул –
объединения под влиянием ван-дер-ваальсовых сил притяжения.

Уравнение состояния Майера [7]:

d(i=dqi1*…dqin.

Здесь Uпij – взаимная потенциальная энергия i-й и j-й молекул,
взаимодействующих по закону центральных сил, qi1,……,qin – обобщенные
координаты i-той молекулы, обладающей n степенями свободы.

Уравнение Камерлинг-Оннеса (1901) [8]:

.

Уравнение Редлиха-Квонга (1949 г.) [8]:

0,42748·R2·T2,5k/Pk, b = 0,08664·R·Tk/Pk. Уравнение Редлиха-Квонга
считается наилучшим двухконстантным уравнением. При его выводе авторы не
руководствовались какими-то определенными теоретическими обоснованиями
[8]. Это уравнение следует рассматривать как произвольную, но удачную
эмпирическую модификацию предшествующих уравнений состояния.

Уравнение Мартина (1967 г.) [8]:

27·R2·T2k/(64Pk), b = R·Tk/(8Pk).

1.2. Основные уравнения, описывающие процесс фильтрации газа в пористой
среде

В последнее время наблюдается рост интереса к различным
термодинамическим эффектам в пористых средах. Это связано с их
многообразными практическими приложениями[4,5].

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

с учетом термодинамических эффектов высокого порядка [9]

(I.2.1)

где первое слагаемое в левой части уравнения (I.2.1) описывает изменение
температуры в пласте со временем, второе – за счет конвекции
(перемещения больших объемов газа). Первое слагаемое в правой части
ответственно за теплопроводность, второе – за межфракционный теплообмен,
третье описывает адиабатический эффект, четвертое – эффект
Джоуля-Томсона и пятое – влияние поля тяготения Земли.

Вторым уравнением системы является уравнение неразрывности, которое
записывается в виде:

. (I.2.2)

Фильтрация газа подчиняется закону Дарси

. (I.2.3)

К системе добавляется уравнение состояния

. (I.2.4)

Система (I.2.1)-(I.2.4) является нелинейной, кроме того, уравнения
(I.2.1)-(I.2.2) являются взаимосвязанными.

1.3. Описание задачи

Рис. 1. постановка задачи

При описании температурной задачи примем следующие допущения:

пористый пласт считается однородным и изотропным по гидродинамическим и
теплофизическим свойствам;

давления в скважине и на контуре питания остаются неизменными;

породы, окружающие пласт предполагаются непроницаемыми и однородными по
своим теплофизическим свойствам;

температуры газа и скелета пористой среды в каждой точке совпадают;

естественное тепловое поле Земли считается стационарным;

пласт расположен на глубине порядка 1 – 2 км, поэтому суточные и
сезонные колебания температуры не достигают пласта;

адиабатическим эффектом, обусловленным гравитационным полем
пренебрегаем.

1.4. Математическая постановка задачи

Математическая постановка задачи включает температурную задачу,
гидродинамическую задачу, уравнение состояния и соотношение для поля
скорости конвективного переноса тепла. Ниже рассматриваются
соответствующие постановки задач.

1.4.1. Математическая постановка температурной задачи

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

. (I.4.1.1)

Будем рассматривать задачу при следующих условиях температуры:

начальном

, (I.4.1.2)

и граничном

. (I.4.1.3)

1.4.2. Математическая постановка гидродинамической задачи

Математическая постановка гидродинамической задачи в полярной системе
координат примет следующий вид. Учитывая, что для осесимметричного
течения поле давления является функцией координаты r уравнение можно
представить в виде:

, (1.4.2.1)

Будем рассматривать задачу при следующих условиях. Пусть PC – давление
на границе контура питания. При значении радиуса, равном радиусу контура
питания

, (1.4.2.2)

давление поддерживается равным Рс:

, (1.4.2.3)

Pс – давление на контуре питания.

При значении радиуса, равном радиусу скважины

, (1.4.1.3)

давление поддерживается равным PW:

, (1.4.1.4)

где PW – давление в скважине.

1.4. Основные идеи метода характеристик[6]

В данном разделе рассмотрим метод характеристик. Любое линейное
дифференциальное уравнение второго порядка (при двух независимых
переменных) может быть записано в следующем виде:

(1.4.1)

где а, b, с, d, e, f, g — заданные непрерывные функции от x и y (или в
частном случае, постоянные).

Попытаемся упростить это уравнение с помощью замены независимых
переменных:

) нигде в области G не обращался в нуль.

Для того чтобы сделать требуемую замену переменных, выразим частные
производные от функции u по х и у через производные от и по ( и (:

(1.4.31)

, через производные по ( и (, учтем формулу (1.4.31) и применим снова
правило дифференцирования сложной функции:

Следовательно,

(1.4.41)

Аналогично найдем:

(1.4.42)

Подставляя u’x, u’y, u’xx,… из этих формул в уравнение (1), мы
получим снова линейное уравнение второго порядка с неизвестной функцией
и и независимыми переменными( и (:

(1.4.5)

где

— функция, линейная относительно и’( , u’( , u .

Уравнение (1.4.5) становится особенно простым, если в нем коэффициенты а
и с окажутся равными нулю. Для того чтобы первоначально заданное
уравнение (1.4.1) можно было привести к такому простому виду, надо в нем
сделать замену переменных

подобрав функции ( и ( так, чтобы они являлись решениями уравнения:

(1.4.6)

Это уравнение является нелинейным уравнением в частных производных
первого порядка. Следующая теорема покажет, как связаны его решения с
общим решением некоторого обыкновенного уравнения.

Теорема. Для того чтобы функция z = f(x, у) во всех точках области G
удовлетворяла уравнению (6), необходимо и достаточно, чтобы, семейство

(1.4.7)

было общим интегралом уравнения

(1.4.8)

в той же области G.

Доказательство. Необходимость. Пусть z = f(x, у)— решение уравнения
(1.4.6). Рассмотрим семейство кривых f(x, у) — k и докажем, что любая
кривая этого семейства удовлетворяет уравнению (1.4.7).

В любой точке, лежащей на кривой f(x, у) = k (где k — фиксировано),
выполняется следующее равенство:

действительно вдоль данной кривой функция f(x, у) постоянна, и поэтому
ее полный дифференциал равен нулю.

Следовательно, всюду на кривой имеет место равенство:

обозначим каждое из этих отношений через (; тогда

Подставляя эти выражения для dx и dy в левую часть уравнения (1.4.8),
получим:

Выражение, стоящее в квадратных скобках, равно нулю, так как, по
условию, функция f(x, у) есть решение уравнения (1.4.6). Следовательно,
во всех точках нашей кривой имеет место равенство

откуда вытекает, что она является интегральной кривой уравнения
(1.4.8).

Итак, любая кривая вида f(x, у) = k является интегральной кривой
уравнения (1.4.8); с другой стороны, через каждую точку области G
проходит кривая такого вида; это вытекает из того, что функция f(x, у)
определена всюду в области G и поэтому, например, через точку (х0, у0)
проходит кривая f(x,y)=f(x0,y0).

Отсюда следует, что семейство f(x, у) = k является общим интегралом
уравнения (1.4.8).

Достаточность. Пусть семейство f(х, у)= k будет общим интегралом
уравнения (1.4.8). Возьмем произвольную точку (х0, у0) из G и выделим ту
кривую семейства, которая проходит через эту точку:

f(x, у) = k0.

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

откуда

(1.4.10)

Так как кривая является интегральной кривой уравнения (1.4.8), то при
подстановке в это уравнение dx и dy из (1.4.10), получим тождество:

или, после сокращения на (2:

В частности, в точке (х0, у0) имеет место:

Но последнее равенство означает, что функция двух переменных f(x, у)
удовлетворяет в точке (х0, у0) уравнению (1.4.7). Так как точка (х0, y0)
была взята произвольно в области G, то функция f(x, у) удовлетворяет
уравнению (1.4.7) во всех точках этой области, т. е. эта функция
является одним из решений уравнения (1.4.7).

Таким образом, теорема доказана.

Рассмотренная теорема открывает путь для упрощения исходного уравнения
(1.4.1). Для этого сначала составляем вспомогательное уравнение (1.4.8);
оно называется характеристическим уравнением для данного уравнения
(1.4.1). Характеристическое уравнение есть обыкновенное дифференциальное
уравнение первого порядка, но второй степени. Разрешая его относительно
y’x (предварительно разделив все члены уравнения на dx2), получим два
уравнения:

(1.4.101)

(1.4.102)

(предполагается, что ас — b2<0, b2 —ас>0 всюду в области G). Пусть общий
интеграл уравнения (1.4.101) имеет вид

((х, у)= k , (1.4.111)

а общий интеграл уравнения (1.4.102)

((х, у)= k. (1.4.112)

Интегральные кривые характеристического уравнения (т. е. все кривые,
входящие в семейства (1.4.111) и (1.4.112)) называются характеристиками
заданного дифференциального уравнения (1.4.1). В связи с этим
рассматриваемый метод упрощения уравнения (1.4.1) называется методом
характеристик.

Семейства (1.4.111) и (1.4.112) можно рассматривать, как общие интегралы
уравнения (1.4.8) (это уравнение распадается на два уравнения (1.4.101)
и (1.4.102)).

Следовательно, согласно доказанной теореме, функции

z=((х, у) и z=((х, у)

являются решениями уравнения в частных производных (1.4.6).

Функции ((х, у) и ((х, у) независимы друг от друга (можно доказать,
что их якобиан отличен от нуля, если ас- b2<0). Поэтому, возвращаясь к уравнению (1.4.1), мы можем в нем сделать замену переменных:. Следовательно, уравнение (1.4.1) преобразуется к виду:или, после деления на 2b и переноса в другую часть равенства:– функция, линейная относительно и’( , u’( , u (см. выше, формула (1.4.5)).Полученное уравнение имеет более простой вид, чем исходное уравнение (1.4.1); если мы его сможем решить (т. е. найти и как функцию от ( и (), то для того, чтобы найти решение исходного уравнения, достаточно вернуться к старым переменным (выразив ( и ( через х и у).1.5. ВыводыВ данной главе представлены основные уравнения состояния реального газа, уравнения, описывающие процесс фильтрации газа в пористой среде. Дано описание задачи. Сформулированы физическая и математическая постановки температурной и гидродинамической задач. Описан метод характеристик, который использован для получения решения задачи в главе 2.Глава 2. Аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состоянияВ данной главе приведены аналитические решения гидродинамической и температурной задач для произвольного уравнения состояния.2.1. Решение гидродинамической задачиРешение уравнения (I.4.1.1) приводит к необходимости решения дополнительной гидродинамической задачи для отыскания поля давления. Для описания движения газа воспользуемся квазистационарным уравнением неразрывности:. (2.1.1)Скорость фильтрации газа сквозь пористую среду определяется законом Дарси:не меняются при движении газа к скважине, и что плотность газа зависит только от давления (баротропное приближение), перепишем уравнение неразрывности в следующем виде:(2.1.3)Функцию Лейбензона представим в виде:и А задаются граничными условиями. Подставляя функцию Лейбензона в уравнение (2.1.3), получим:, уравнение (2.1.5) можно представить в виде:(2.1.6)Решение этого уравнения представим в виде:), тогда согласно (2.1.4) и (2.1.7):(2.1.10)(2.1.11)Подставив эти выражения в уравнение (2.1.7), получим:Используя (2.1.12), получим выражение для градиента давления в виде:(2.1.13)Найденное выражение для градиента давления позволяет преобразовать температурную задачу к виду удобному для аналитического решения. Решение температурной задачи рассматривается в следующем параграфе.2.2. Решение температурной задачиС учетом (2.1.13) и закона фильтрации Дарси (2.1.2) задача (I.4.1.1) – (I.4.1.3) преобразуется к виду:(2.2.1)Условия (I.4.1.1) – (I.4.1.3) остаются неизменными. Вводя обозначение(2.2.2)представим уравнение Чекалюка в виде:(2.2.3)Будем рассматривать задачу при следующих условиях для температуры:начальном(2.2.4)и граничном::(2.2.7)определяет область применимости нестационарного решения(2.2.8)Уравнение (2.2.3) с учетом (2.2.6) можно представить в виде:(2.2.9)откуда(2.2.10)где(2.2.11)Исключив константу С из (2.2.10) и (2.2.11), окончательно получаем нестационарное решение:. Для этой области(2.2.13)и стационарное решение, удовлетворяющее условию (2.2.5), имеет вид:(2.2.14)Выражения (2.2.12) и (2.2.14) полностью решают поставленную задачу для любого уравнения состояния2.3. ВыводыВ данной главе представлено аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состояния, которая включает в себя температурную и гидродинамическую задачи.Глава 3. Получение Аналитических выражений решения задачи о баротермическом эффекте с учетом барической сжимаемости3.1. Решение гидродинамической задачи для линеаризованного уравнения состоянияВыпишем полученные решения для линеаризованного баротропного уравнения состояния(3.1.1)Вычислив интеграл, входящий в (2.2.3):(3.1.2)Представим зависимость между давлением и радиальной координатой r в виде:(3.1.3)Введем обозначение(3.1.4)Тогда уравнение (3.1.3) преобразуется к виду:(3.1.5)Откуда найдем(3.1.6)Физический смысл имеет только значение полученного выражения со знаком плюс перед квадратным корнем. Введем обозначения(3.1.7)и упростить запись выражения (3.1.6)(3.1.9)Подставив (3.1.9) в (3.1.1), получим зависимость плотности от радиальной координаты r:(3.1.10)Полученные в данном разделе выражения позволяют построить решения задачи о баротермическом эффекте в случае линеаризованного уравнения состояния.3.2. Температурная задача в линеаризованном случаеВ этом случае нестационарное решение для температуры (3.1.5) записывается в виде:(3.2.1)Интеграл в (3.2.1) легко вычисляется; окончательно нестационарное решение представляется в виде(3.2.2)Выражения для G и H представляются формулами (3.2.5) и (3.2.6), а - для V представляется в виде, следующем из(2.2.8)(3.2.3)В пределе при ??0 из (3.2.2)-(3.2.3) следует известное решение для несжимаемой жидкости[4]:(3.2.4)Аналогично в стационарном случае из (2.2.14) получим:(3.2.5)В пределе при ??0 из (3.2.5) и (3.2.3) следует известное решение для несжимаемой жидкости[4]:(3.2.6)Выражения (3.2.2), (3.2.4) решают поставленную задачу о баротермическом эффекте при фильтрации газа в прискважинной зоне реальных газовых пластов. Такое решение поставленной задачи получено впервые. Поэтому представляет значительный и практический интерес анализ результатов расчетов на основе полученных решений, что и приведено в четвертой главе.3.3. ВыводыВ данной главе получено аналитическое решение задачи о баротермическом эффекте с учетом барической сжимаемости, которая включает в себя решение гидродинамической задачи для линеаризованного уравнения состояния и температурную задачу в линеаризованном случае.Глава 4. анализ результатов расчетов и Исследование температурных полей, возникающих при фильтрации газаВ данной главе приведен анализ результатов расчетов баротермического эффекта в прискважинной зоне газовых пластов применительно к реальным месторождениям газа.4.1. Анализ результатов расчетов температурных полей.Из рисунка видно, что изменение температуры подчиняется следующим закономерностям. Линейное нарастание температурного эффекта при малых временах сменяется логарифмической стабилизацией - при больших временах. Время, при котором происходит смена линейного нарастания на логарифмическую стабилизацию, зависит от барической сжимаемости; с увеличением сжимаемости это время уменьшается.Величина температурного эффекта также сильно зависит от сжимаемости. С увеличением сжимаемости величина температурного эффекта возрастает. Коэффициент барической сжимаемости приблизительно обратно пропорционален давлению. Реальные значения этого коэффициента в условиях газовых пластов лежат в пределах от 3 10-8 Па-1 до 10-5; поэтому величина эффекта лежит в пределах до –10 ( –15 К.. Это хорошо согласуется с величиной измеряемых в скважинных условиях температурных эффектов.Рис.1. Зависимость величины баротермического эффекта от времени при различных барических сжимаемостях. Обозначения: 1 - ( = 3 10-4 Па-1, 2 – 10-5, 3 – 10-6, 4 – 10-7, 5 – 5 10-8Важно отметить, что согласно разработанной нами теории время установления температурного эффекта при ( ( 10-8 Па-1, что часто встречается на практике, составляет около суток. Этот факт чрезвычайно важен при практическом использовании баротермического эффекта....Рис 4. Зависимость нестационарной температуры от времени при различных радиусах скважины. Обозначения: 1- r =0.1 м; 2 -0.2 ; 3 – 0.3; 4 -0.5..Рис 5. Зависимость нестационарной температуры от радиуса скважины при различных временах. Обозначения: 1- t =10 000 с; 2 -100 000 ; 3 – 1 000 000...Из рисунка видно, что величина температурного эффекта возрастает при увеличении теплоемкости.Рис 6. Зависимость нестационарной температуры от времени при различных радиусах контура питания. Обозначения: 1- R =25 м; 2 -50; 3 – 100; 4 -200; 5 - 250.Рис 7. Зависимость нестационарной температуры от теплоёмкости при различных временах. Обозначения: 1- t =100 000 c; 2 -1 000 000; 3 – 10 000 000..Рис 8. Зависимость нестационарной температуры от относительной вязкости при различных временах. Обозначения: 1- t =100 000 c; 2 -1 000 000; 3 – 1 500 000.На рис. 9. показана зависимость баротермического эффекта от времени при различных коэффициентах барической сжимаемости. Из рисунка видно,.4.2. Изучение вклада сжимаемости в величину баротермического эффекта.Рис. 10. Зависимость нестационарной температуры от коэффициента барической сжимаемости при различных временах. Обозначения: 1- t = 100 c; 2 -1000 ; 3 – 10 000; 4 -100 000.На рис. 11. показана зависимость стационарной температуры от коэффициента барической сжимаемости. Из рисунка видно что величина температурного эффекта в стационарном случае не зависит от коэффициента барической сжимаемости. В расчетах принято: ?=-0.5?10-5; rW=0.1; с=850; k=10-15; сPL=84000000; µ=10-5; R=100; ?=150; P=100?105; P0=150?105; PC=200?105; PW=150?105.Рис. 11 Зависимость стационарной температуры от коэффициента барической сжимаемости.На рис. 12. приведена зависимость времени установления температуры от коэффициента барической сжимаемости.Рис 12. Зависимость времени установления температуры от коэффициента барической сжимаемости.Итак, изучение вклада сжимаемости в величину баротермического эффекта показывает, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.4. 3. ВыводыВ данной главе сделан анализ результатов расчетов и исследованы температурные поля, возникающих при фильтрации газа. Показано, что величина температурного эффекта составляет около 20 К. Время установления температурного эффекта сильно зависит от проницаемости и для реальных значений проницаемости составляет приблизительно сутки. Это важно учитывать при интерпритации результатов термических исследований скважин.Изучен вклад сжимаемости в величину баротермического эффекта. Показано, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.Показано, что время установления баротермического эффекта зависит от барической сжимаемости и лежит в пределах до 109 с при ?~10-8 Па-1. При ?~10-8 Па-1 время полного установления составляет (приблизительно) три года. Значит температурные поля в газовом пласте практически всегда нестационарны.ЗаключениеВ ходе проделанной работы были получены следующие результаты:Описаны основные уравнения состояния реального газа, уравнения, описывающие процесс фильтрации газа в пористой среде.Представлено аналитическое решение задачи о баротермическом эффекте с учетом реального уравнения состояния.Получено аналитическое решение задачи о баротермическом эффекте с учетом барической сжимаемости.Сделан анализ результатов расчетов и исследование температурных полей, возникающих при фильтрации газа.Исследованы температурные поля и изучен вклад сжимаемости в величину баротермического эффекта. Показано, что в нестационарных полях величина температурного эффекта сильно зависит от сжимаемости, а после установления температуры не зависит от сжимаемости.При ?~10-8 Па-1 время полного установления температуры составляет (приблизительно) три года. Это означает, что температурные поля в газовом пласте практически всегда нестационарные. Следует отметить при этом что логарифмическая стабилизация достигается при времени около суток.Список использованной литературыЛандау Л. Д., Лившиц Е. М. Статистическая физика// М.,1964.Карслоу Г., Егер Д. Теплопроводность твердых тел// М: Наука. 1964. 487с.Варгафтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей// М., Наука. 1972.Филиппов А. И., Фридман А. А., Девяткин Е. М. Баротермический эффект при фильтрации газированной жидкости: Монография. - Стерлитамак: Стерлитамак. гос. пед. ин-т; Стерлитамакский филиал Академии наук Республики Башкортостан, 2000. – 175с.Филиппов А. И. Скважинная термометрия переходных процессов. - Саратов: Изд-во Сарат. ун-та, 1989. – 116с.Очан Ю. С. Методы математической физики// М: Высшая школа. 1965. 383с.Яворский Б. М., Детлаф А. А. Справочник по физике. – М.: Наука, 1971. – 940с.Морачевский А. Г., Сладков И. Б. Физико – химические свойства молекулярных неорганических соединений. – С. Пб.: Химия, 1996. – 312с.Баскаков А. П., Гуревич М. И., Решетин Н. И. и др. Общая теплотехника. – М.-Л.: Государственное энергетическое издательство, 1963. – 392с.PAGE 1PAGE 42-25.15?10-4?-24.9-25.0?T,K431210-38?10-46?10-41234?10-42?10-48?1046?10444?1042?104-20-16-12-8-4t,c001?1055?T,K0-20-15-10-5??T,K0t,c?T,K7?105?-4-3-2-10-5-60-741231?1055?1056?1052?1053?1054?10500?T,K1234t,c-2-10-8-6-42?1054?1056?1058?1051?1060?T, Kt-4-24321-60-85?1051?1060r123-5-10-15-200?T, K0.511.50?T, K-2-4-6-8-101325?1050541?106t, c0-5321с10005000-10?T, K0-5321-10-15-20-2502?10-64?10-66?10-68?10-61?10-7?T, Kµt, c-1510-30-25-2023451·10102·10103·10104·1010?T, K?, с6?1084?1082?1081?10-31?10-85?10-4?0?

Нашли опечатку? Выделите и нажмите CTRL+Enter

Похожие документы
Обсуждение

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

avatar
  Подписаться  
Уведомление о
Заказать реферат
UkrReferat.com. Всі права захищені. 2000-2019