Материал: Разработка нефтяных месторождений.

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

Уравнение (11.42) можно записать также в следующем об­

щем виде:

 

div (рр)+ д(^W) ■= 0.

(Н.43)

Уравнения (11.42),

(11.43)— уравнения неразрывности мас­

сы вещества во время его движения при трехмерном измерении. Если в пласте одновременно движутся несколько веществ, на­ ходящихся как в газовой, так и в жидкой фазе, составляют уравнения неразрывности массы каждого вещества (компонен­ та) в соответствующих фазах.

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

ской энергии веществ, движущихся в пласте

со скоростью ш..

Поэтому

 

Eu= ua-{-z-{-w2/{2g).

(II. 44}

Из закона сохранения энергии или, точнее, из первого нача­ ла термодинамики следует, что изменение энергии пласта ДЕп и произведенной удельной работы равно количеству подве­ денного к пласту тепла 6QT, умноженного на механический эк­ вивалент тепла А, т. е.

AEn+6W = A6QT

(11.45)

или с учетом (П.44)

 

A/«II+ z + -g -)+ 6 W ' = 4«QT.

(11.46)

Дадим количественную оценку входящих в (11.46) величин. Удельная внутренняя энергия пласта ип при отсутствии в нем химических или ядерных превращений вещества представляет собой тепловую энергию в единице массы пласта, так что

Аип= АсАТ,

(11.47)

где с — удельная теплоемкость пласта; Т — температура. Поло­ жим, что пористый пласт насыщен водой. Тогда с=ст(1—т ) + А~свт (с-г — удельная теплоемкость пород пласта, св— удельная теплоемкость воды, m — пористость). Пусть ст= 1,046 кДж/(кг* •К), св = 4,184 кДж/(кг-К), АТ = 1 К, т = 0,2. Тогда с= 1,046• (1— —0,2) +4,184-0,2= 1,67 кДж/(кг-К), Аип= 102-1,67-1 = 170 м. Удельная потенциальная энергия z в пластах может изменяться

всоответствии с возможными изменениями уровня движущихся

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

гии. Скорость w движения в пласте насыщающих его веществ

61

изменяется в значительных

пределах — от 0 до

10 м/сут =

= 3650 м/год= 1,16-10-4 м/с.

Сравнивая удельные

потенциаль­

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

эффициент е = где рт — плотность горных пород;

рв — плотность насыщающих пласт веществ, и умножать все ви­ ды удельной энергии, кроме внутренней, на е. При рв = 103кг/м3, рт = 2,25 • 103 кг/м3 т = 0,2, е= 0,1.

Тогда для изменения удельной кинетической энергии полу­ чим

= 0,68-10-10 м.

2-9,81

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

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

W = p6AV/(pgAV),

 

(11.48)

где р — давление; AV — объем вещества, насыщающего пласт в

элементарном объеме пласта; р — плотность

этого

вещества;

g — ускорение свободного падения.

говоря,

неизмен­

Поровый объем пласта остается, вообще

ным, поскольку не изменяются геометрия пласта и его порис­ тость. Работа вещества в пласте связана всегда с его расшире­ нием. Поэтому в (11.48) и введена величина 6ДУ, характеризу­ ющая расширение вещества. При этом условно можно считать, что вещество, насыщающее пласт, расширяясь, как бы выходит за пределы элементарного объема пласта. Будем считать, что при бесконечно малом расширении вещества в элементарном объеме пласта масса вещества ДМ —рДУ остается неизменной.

62

Тогда бАМ= 6рАУ+рбАУ=0 и, следовательно,

6ДУ/ДУ = —бр/р.

 

(11.49)

Подставляя

 

(11.49)

и (11.48)

получим

6W ------4 е—

- & ( — )

(11.50)

р2g

g

I p

)

 

Оценим возможную работу вещества, насыщающего пласт. Очевидно, что наибольшую работу может производить в пласте газ. Для простоты оценки будем считать газ идеальным, для которого р/р = ро/ро, где ро, ро — давление и плотность газа при начальных условиях. Отсюда для идеального газа

е8Г------(11.51)

Р

'

'

Роs

Пусть при

снижении давления

бр = —10-105 Па,

р =

= 100♦ 105 Па, ро= 105 Па, ро=1 кг/м3, е= 0,1. Тогда

0,1-106.10-106

EW = 1-9,81-100-Ю5 = 102 м.

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

Рассмотрим

Boinpoc о

том, чему равняется входящая в

(11.45) и (11.46)

величина

бQT. Тепловыделение в элементе пла­

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

сенной

к единице

 

массы

движущегося

вещества в элементе

пласта, имеем следующее выражение:

 

 

AN

1

 

 

 

 

|XUa

'

 

 

(11.52)

pgWn

mpg

v grad p = mpgk '

 

 

Допустим,

что

 

в

 

пласте

движется

газ вязкостью

ц =

= 0,02-10-3 Па-с

со скоростью и=10-6

м/с« 8 6 ,4 -10_3

м/сут.

Проницаемость пласта

£ « 0 ,1 2

мкм, пористость т = 0,2,

плот­

ность газа р

при

давлении

р = 100 МПа

составляет 100

кг/м3.

Тогда

0,02 -10_3-10~12

 

 

 

 

 

 

\IVZ

=

1,02-10"6 м/с.

 

 

mpgk

0 ,2 -10—13-981

 

 

 

63

В сутки из килограмма движущегося в пласте газа будет выделяться 1,02 ■ 10~6- 0,864 -105= 0,088 м энергии. Это, конечно, незначительная величина. Однако, например, в призабойной зоне скважин скорость фильтрации того же газа может дости­ гать 10~4 м/с, а иногда и более. Тогда при тех же остальных условиях, что и выше, значение \LV2/{mpgk)« 10-2 м/с. В сутки из килограмма фильтрующегося в пласте газа выделится энер­ гии почти 9 кДж. Таким образом, можно заключить, что наи­ более существенное изменение энергии в элементе пласта свя­ зано с переносом тепла за счет теплопроводности и конвекции. Определенный вклад в энергетический баланс пласта, особенно при высоких скоростях движения насыщающих его веществ, вносят работа расширения-сжатия веществ и гидравлическое трение.

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

В соответствии с (11.48) и (11.49) работу движущегося ве­ щества в элементарном объеме пласта в целом можно пред­ ставить в следующем виде:

W = m S W = m p - ^ - = —m p -^ -.

(11.53)

Работу W' можно приравнять к энергии сжатия Ер, поэтому

W = -m 6 E p

рбр

(11.54)

Р2

 

 

PI

где pi и Р2 —плотности.

Рассматривая, как и при выводе уравнения неразрывности массы фильтрующегося в пласте вещества, поток внутренней энергии и=срТ и энергии сжатия Ер, а также считая, что тепло поступает в элементарный объем только за счет гидравличе­ ского трения, т. е. что j46^T = ^gradр, получим

Л

j-divo£«j = m ^ ^ 5 - + d iv £ ppo^ = t>grad р.

(11.55)

Здесь

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

в пласте

за счет теплопроводности

и конвекции, v — вектор

скорости

фильтрации. Выражение

(11.55) и есть дифференциальное

уравнение сохранения энергии в пласте, выведенное при указан­ ных выше предположениях.

Рассмотрим законы фильтрации. Основным законом под­ земной гидромеханики является закон фильтрации однородной жидкости или газа — з а к о н Да р с и . Все известные законы фильтрации базируются на этом основном законе.

При фильтрации неоднородной жидкости или смесей жид­ кости и газа справедлив закон двухфазной фильтрации. В слу-

4>4

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

н

ккп (s)

дРп

 

 

 

и*

дх

’

 

 

____ ккъ (s)

дрв

 

 

 

v ‘ ~

и.

дх

’

 

 

где

Ун — вектор

скорости фильтра­

 

ции нефти; vB— вектор фильтрации

 

воды, &н (з),

kji(s)— относитель­

 

ные

проницаемости

соответственно

Рис. 40. Графики зависимости

для

нефти и воды,

зависящие от

&и и kBот 5

водонасыщенности s; рв и рв — дав-

ления в нефти и воде. Графики относительных проницаемостей для нефти и воды имеют вид, показанный на рис. 40, на кото­ ром по оси абсцисс отмечены две характерные точки: sCBи з*: В точке з=Зев относительная проницаемость для воды равна нулю, так что kB(sen) =0. В точке s=s* относительная прони­

цаемость для нефти kB(s*)= 0,

несмотря на то что в точке s=

= з св в пласте присутствует

вода, а в точке з=з* имеется

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

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

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

насыщенности пористой среды газом зг и жидкостью зж. По­ скольку

 

(И.57)

где Sa, sm— соответственно насыщенности пласта

водой и

нефтью, можно шашписатъ следующие выражения:

 

S = —

(11.58)

5 Ю.. Ш.Ж®ашив

Источник: https://tut-files.ru/previewfile/82937