Материал: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ В ЗАДАЧАХ НЕФТЕГАЗОВОЙ ОТРАСЛИ

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

СПБГУАП| Институт 4 группа 4736

при x 0, t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

G1 x, t G1 x x, t G3 x, t G3 x x, t

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

u x, t t u x,t S lim

G4

S 2

 

 

 

t 0

 

 

t

 

 

 

 

t 0

t

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G1 x, t G3 x,t

 

 

 

 

 

 

 

 

 

 

 

(2.125)

 

 

u x,t

 

G4 x,t

 

 

 

 

 

 

 

 

 

S

 

 

 

S

 

 

 

 

 

 

Объем жидкости,

 

x

 

 

 

t

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

поступающей в сечение 1—1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

будет —.

 

G2

.Следовательно, в объеме

 

 

G2

— содержится масса

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

растворенного газа G

 

G2

, т. е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

2

 

 

 

 

 

 

 

 

 

G G

G2

.

 

 

 

 

 

 

 

 

 

 

 

(2.126)

 

 

 

 

 

 

 

 

 

 

 

 

3

4

2

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в трех уравнениях (2.122), (2.125) и (2.126) имеются неизвестные: Gx, G2, G3, G4, S2 и 1 .

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

Через сечение 1—1 в единицу времени прошли масса газа G1 и масса конденсата G2. Через сечение 2—2 отобраны масса газа G12 и

масса конденсата G22. Тогда масса конденсата, накопившегося в объеме в единицу времени,

G2 G22 G2 .

Это приращение массы конденсата приведет при условии несжимаемости конденсата ( = const) к изменению объема в единицу

времени на величину G12 .

2

За отрезок времени t изменение объема будет G12 t . Это

2

изменение объема приведет к изменению насыщенности

86

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

конденсатом на величину S2 .

Аналогично предыдущему случаю запишем

G2 x x,t G2

x,t

t S2

x,t t S2

x,t F x

 

 

 

 

 

(2.127)

2

 

 

 

 

 

 

 

Знак минус перед правой частью поставлен, исходя из следующих соображений: если G2 x x,t G2 x,t , то насыщенность конденсата в объеме

V уменьшается, т. е.

S2 x,t t S2 x,t , и наоборот.

Разделив обе части уравнения (2.127) на x

и t и перейдя к пределу при

x 0, t 0 , получим

 

 

G2 x,t

F 2

S2 x,t

 

(2.128)

 

x

t

 

 

 

 

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

G1 x,t

F

1

(2.129)

x

t

В уравнениях (2.128) и (2.129), описывающих процесс движения газоконденсатной системы, имеются четыре неизвестных.

Перейдем к примерам применения закона сохранения энергии для механических и тепловых процессов.

Как известно, для механических процессов сумма кинетической и потенциальной энергии постоянна.

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

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

С поглощением тепла в количестве dQ

температура повышается на

величину dT. Отношение

c

dQ

(где т

— масса тела) называется

mdT

 

 

 

 

теплоемкостью тела. В физике пользуются теплоемкостью при постоянном давлении ср и теплоемкостью при постоянном объеме cV .

87

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

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

Температура в какой-либо точке тела будет зависеть от расстояния точки до места подвода тепла и времени

Т = Т(х, t).

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

В приведенном примере через сечение 1—1 в единицу времени подводится количество тепла, равное Q (х, t). Через сечение 2—2 в единицу времени отводится количество тепла, равное Q x x,t . Разность

между этими количествами теплаQ Q x x,t Q x,t (2.130)

Рис. 2.8

в соответствии с первым законом термодинамики затрачивается на изменение температуры в отсеке между сечениями 1—1 и 2—2.

Если

Q 0 , то тело нагревается, в противном случае — охлаждается,

т. е. в первом случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

T

 

x,t t

 

T

 

x,t

 

0

,(2.131)

а во втором — T 0 .

При изменении количества подводимого и отводимого тепла па величинуQ температура тела (в соответствии с законом сохранения энергии и определением теплоемкости) должна изменяться на величину

T

Q

 

t

 

 

m1cV

 

(2.132)

 

 

 

 

 

m1 xS

 

 

 

 

где m1 — элементарная масса тела; S — площадь поперечного

сечения тела;

— плотность тела.

 

Так как Q 0

при T 0 и Q 0 при

T 0 , то в правой части

(2.132 ) поставлен знак минус.

Подставив в (2.130) значения Q и T соответственно из (2.130), (2.131

88

Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736

) н разделив его на x, t и

перейдя к пределу при x 0, t 0 ,

получим

 

 

 

1

Q

T

(2.133)

 

 

S c

t

 

 

x

 

 

 

V

 

 

 

 

Таким образом, в одном уравнении получаем два неизвестных Q и Т.

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

и Q2.

При изотермическом совместном движении газа и жидкости в

трубопроводе

неизвестными были и V, а уравнение — одно. Но, как

известно,

 

связано с давлением при помощи уравнения состояния

Клайперона — Менделеева

 

p

,

 

2RTg

где р — давление.

Однако при этом появляется новое неизвестное р, т. е. опять не хватает одного уравнения. Приведенное выше уравнение соответствовало закону постоянства массы. Для исследуемого явления необходимо применить еще закон сохранения энергии, например в виде уравнения Бернулли. Но в этом случае появится новое неизвестное — работа сил трения, т. е. требуется опять одно уравнение для определения силы трения. Это уравнение устанавливается экспериментально.

Для совместного движения несжимаемых жидкостей закон сохранения энергии приводит также к необходимости определения силы трения. Указанный процесс характеризуется и давлением. Поэтому и в этом случае необходимо экспериментально составить одно уравнение для сил трения. Если движение неизотермическое, то следует применить и первый закон термодинамики. Тогда появятся два неизвестных Q и Т, связанных одним уравнением типа

1

Q

 

T

,

S c

x

t

 

 

V

 

 

 

 

т. е. и в этом случае необходимо определить одну экспериментальную зависимость для Q.

При изотермическом движении газированной жидкости в трубопроводе или при фильтрации газированной жидкости в пористой среде плотность газа р связана с давлением р. Появляется еще одно неизвестное р. Вводя понятие массовой скорости для жидкости w1 и газа w2, можно через них вычислить соответственно G1 и G2. Неизвестными в этом случае являются w1, w2 ,G3 ,G4 , S2 , , p . Таким образом, необходимо экспериментально определить

зависимости для w1, w2 и G3.

89

Контакты | https://new.guap.ru/i03/contacts

2 rdr

СПБГУАП| Институт 4 группа 4736

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

уравнения состояния.

Таким образом, необходимо экспериментально определить две зависимости — для G1 и G2.

При исследовании процесса распространения тепла в теле имеем два неизвестных Q и Т и одно уравнение. Закон постоянства массы выполняется тождественно и недостающее уравнение для Q определяется экспериментально.

Пример Движения в круглой цилиндрической трубе вдоль кольца.

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

инерции 2 rl vt .

По принципу Д'Аламбера, сумма этих сил должна равняться нулю. Если при этом учесть, что первая и последняя из этих сил действуют в направлении, обратном действию остальных двух сил, то

2 r dr l d 2 rdr 2 r 2 rl

v

0

(2.134)

 

t

 

 

Пренебрегая величиной drdx по сравнению с другими членами, получаем

dr rd rdr

 

rdr

v

0

(2.135)

 

l

 

t

 

 

90

Контакты | https://new.guap.ru/i03/contacts

Источник: https://studfile.net/preview/16405273/