Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Основы разработки нефтяных и газовых месторождений

..pdf
Скачиваний:
15
Добавлен:
12.11.2023
Размер:
21.27 Mб
Скачать

4.2. ЗАКОН ДАРСИ. ПОТЕНЦИАЛЬНАЯ ЭНЕРГИЯ ФЛЮИДОВ

В каждой отрасли науки и техники есть свои выдающиеся лично­ сти. Можно вспомнить такие прославленные имена, как Ньютон и Эйнштейн в физике или Дарвин в биологии. В области разработки ме­ сторождений личностью такого масштаба является французский ин­ женер Анри Дарси (XIX век), который, сам того не осознавая, вошел в историю как первый экспериментатор в данной отрасли. В 1856 г. Дар­ си опубликовал подробный отчет о своей работе2 по усовершенствова­ нию водопровода г. Дижона. Входе этих исследований он, в частности, разработал фильтр, рассчитанный на суточный расход воды для горо­ да. И, хотя в это время гидродинамика была уже достаточно хорошо развита, публикаций по движению жидкостей в пористой среде еще не было. Будучи специалистом-практиком, Дарси создал для исследо­ вания этого предмета фильтр, показанный на рис. 4.1.

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

Ь. - Ь,

АЬ

( а л \

и = К - 1---- = К

--------I

(4Л)

I

 

Здесь:

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

1= полная длина песчаной набивки,

К= постоянная.

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

Рис. 4.1. Схема экспериментальной установки Дарси

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

Позже другие исследователи повторили опыты Дарси, изменяя при этом другие параметры. Прежде всего они стали располагать песча­ ную набивку под различными углами к вертикали, как это показано на рис. 4.2. Однако было выявлено, что независимо от угла наклона пес­ чаной набивки разность уровней жидкости в пьезометрах АЬ остава­ лась неизменной для данного расхода. Таким образом, оказалось, что установленный по экспериментальным данным закон Дарси действует независимо от направления потока в гравитационном поле Земли.

плоскость отсчета; 2 = 0, р = 0,1 атм

Рис. 4 .2 . Ориентация установки Дарси в гравитационном поле Земли

Следует отметить значение члена АЬ в уравнении Дарси. Давление

влюбой точке потока, имеющей высоту г относительно плоскости от­ счета (рис. 4.2), будет превышать атмосферное давление на величину

р= Р§ (Ь - г).

Вэтом уравнении Ь представляет собой высоту подъема жидкости

вверхнем пьезометре по отношению к отметке г - 0, а р - плотность жидкости (воды). Это уравнение можно записать в другом виде:

Ь§ =

(4.2)

Если записать уравнение (4.1) в дифференциальной форме

 

аь

(4.3)

и = К -

<11

 

а затем продифференцировать уравнение (4.2) и подставить его в (4.3), то получим

К а(ь§)

(4.4)

8

а!

8 а!

 

Выражение (р / р + §г) в последнем уравнении имеет такую же раз­ мерность, что и Ь§ - произведение силы на расстояние, отнесенное к единице массы, то есть это потенциальная энергия, отнесенная к еди­ нице массы. Обычно эта удельная потенциальная энергия обознача­ ется символом Ф. Она определяется как работа, требуемая для пере­ мещения без трения единицы массы флюида из точки, где давление равно атмосферному и высота равна нулю, в данную точку. Таким об­ разом,

Ф= (•

— +§г.

(4.5)

}

р

 

1 - а1ш

 

 

Несмотря на данное определение, удельную потенциальную энер­ гию флюида не всегда измеряют по отношению к атмосферному давлению и нулевой высоте. Иногда ее измеряют относительно про-

извольно взятых базовых давления и высоты (рь, 2Ь). С учетом изло­ женного, можно записать уравнение (4.5) в другом виде:

р

Рь

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

ра

Рв

Ра

г ф

г ф

Г Ф

Ф А - Ф В = ] ~

+ 8 ( ^ - 2 , , ) - ] ~ +

8 ( 2 в - 2ь ) = ^ — + 8 ( 2 а - 2в).

Рь

Рь

Рв

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

Ф = - + § 2 ,

(4.7)

Р

 

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

Постоянная К / § применима только для потока воды, которая была единственной жидкостью в опытах Дарси. Исследования с различны­ ми жидкостями показали, что закон можно записать в общем виде

кр

(4.8)

и =

 

|! 31

В таком случае зависимость скорости фильтрации от плотности жидкости р и вязкости |Л вполне очевидна. Новая постоянная к здесь зависит только от типа песка, и она характеризует проницаемость. Фактически, это абсолютная проницаемость песка, если он полно­ стью насыщен флюидом. Исходя из того, каким образом определяет­ ся абсолютная проницаемость, она всегда будет иметь одно и то же значение независимо от типа флюида.

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

Благодаря очень низкой вязкости реального газа скорость его дви­ жения в пласте намного превышает скорость движения нефти или воды. В ограниченной области вокруг скважины, где велика депрес­ сия, скорость газа может стать настолько высокой, что течение не подчиняется закону Дарси5. Это явление и его количественное вы­ ражение в уравнениях фильтрации газа будет подробно описано в главе 8 (раздел 6).

4.3. ПРИСВОЕНИЕ ЗНАКОВ

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

Если же уравнение (4.8) применяется совместно с другими мате­ матическими зависимостями, например с выражением для изотер­ мической сжимаемости (см. главу 1, раздел 4), то необходимо при­ своить знаки.

Прямолинейно-параллельный фильтрационный поток

Если расстояние увеличивается в направлении движения, то гра­ диент потенциальной энергии 6Ф / <11 должен быть отрицательным, так как флюиды движутся из зоны, где они обладают более высокой потенциальной энергией, в зону, где они обладают более низкой по­ тенциальной энергией. Поэтому закон Дарси запишется в следую­ щем виде:

ц = _ к р а Ф

(49)

|д (II

 

Плоскорадиальный фильтрационный поток

В этой книге принято, что поток из пласта к скважине является положительным. Расстояние по радиусу считается положительным в направлении, противоположном потоку. Таким образом, <1Ф / бг будет положительным, и закон Дарси имеет вид

кр^ао

(410)

(I <1г

 

4.4. ЕДИНИЦЫ ИЗМЕРЕНИЯ. ПЕРЕХОД ОТ ОДНОЙ СИСТЕМЫ ЕДИНИЦ К ДРУГОЙ

В любой абсолютной системе единиц измерения уравнение Дарси для прямолинейно-параллельного потока записывается в виде (4.9)

кр 6Ф

и = -------- ;—

р 61

Здесь параметры имеют следующую размерность:

и = Ь / Т; р = М / V; ц = М / ЬТ; I = Ь, и Ф (удельная потенциальная энергия) = Ь2 / Т2.

Анализ размерностей в уравнении (4.9)

1 гк1 [м/ь3] [ь2 / т2]

Т

1 ] [М /Ь Т ]

[Ь]

Таким образом, единица проницаемости должна быть выражена в м2 (в системе СИ). Для большинства горных пород эта единица слиш­ ком велика и неудобна для практического применения. Поэтому была разработана система единиц, в которой проницаемость имеет более удобное численное выражение. Это так называемая «система Дарси» (см. табл. 4.1), в которой единицей проницаемости является дарси (Д). Эту единицу можно определить из закона Дарси для горизонтального прямолинейно-параллельного потока несжимаемой жидкости

поэтому проницаемость к =1 Д, когда и = 1 см / с, ц = 1 сП, и ф / <11 = 1 атм / см.

Как следует из табл. 4.1, «система Дарси» - это смешанная система, основанная на системе СГС. Единственное отличие заключается в том, что давление здесь выражается в атмосферах, вязкость в санти­ пуазах и, как следствие, проницаемость в дарси. При создании этой4

 

Обозна­

 

Абсолютные системы

Параметр

Размерность

СГС

СИ

 

чение

 

Длина

1

Ь

см

м

Масса

ш

м

г

КГ

Время

г

т

с

с

Скорость

и

ь / т

см / с

м /с

Расход

Ч

V / Т

см3 / с

м3/ с

Давление

Р

(М1/Г2)

дин / см2

Н / м 2

/V

(паскаль)

 

 

Плотность

Р

м /и

г / см3

КГ / м3

Вязкость

 

м/ьт

г / см с

кг / м с

 

(пуаз)

 

 

 

 

Проницае­

к

V

см2

м2

мость

 

 

 

 

Смешанные системы

«Система

«Промысловая

Дарси»

система»

см

фут

г

фунт

с

час

см / с

фут / с

 

ст. барр. / сут

см3 / с

(жидкость)

тыс. ст. фут3

 

 

/ сут (газ)

атм

фунт / дюйм2

(абс)

 

г / см3

фунт / фут3

сП

сП

(сантипуаз)

 

д

мД

(дарси)

(миллидарси)

Таблица 4 .1 . Абсолютные и смешанные системы единиц измерения,

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

4 Приведенная в правом столбце таблицы «промысловая система единиц» распространена в американской нефтегазовой отрасли. - Прим. ред.

системы единиц преимущества ее виделись не только в том, что еди­ ница проницаемости будет иметь удобное для использования чис­ ленное значение, но и в том, что уравнения, записанные в единицах этой системы, будут иметь такую же форму, что и уравнения, запи­ санные в единицах абсолютной системы. Иначе говоря, уравнения не будут содержать громоздких констант, кроме числа я, введенного для учета геометрии системы. К сожалению, последнее предположе­ ние оправдывается не всегда, поскольку дарси определяется исходя из соотношения (4.11), не выражающего закон Дарси в полной фор­ мулировке. Разумеется, соотношение (4.11) будет иметь одну и ту же форму и в единицах абсолютной системы, и в единицах «системы Дарси», но при использовании закона фильтрации в полной форму­ лировке (4.9) применительно к несжимаемой жидкости

(4.12)

в единицах абсолютной системы, в то время как

(4.13)

в единицах «системы Дарси». Постоянная 1,0133 х 106 представляет собой количество дин / см2, равное одной атм. Она необходима, по­ скольку р и § выражаются в одних и тех же единицах и в системе СГС и в «системе Дарси», и в то же время второй член выражения в скоб­ ках в формуле (4.13) должен быть выражен в тех же единицах, что и первый, то есть в атм / см.

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

Однако при рассмотрении практических аспектов разработки ме­ сторождений, например при анализе результатов исследования сква­ жин, описанном в главах 7 и 8, обычно переходят на так называемую «практическую», или «промысловую», систему единиц/ Эту систе-*

* В настоящее время во многих странах, в том числе в России, принята система СИ, поэтому при переводе этой книги все практические расчеты за пределами данного раздела воспроизводились в системе СИ (в некоторых случаях для расчетов использовались и кратные единицы, такие как МПа, м3/сут, мкм2). - Прим. ред.

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

Так как в отрасли применяются различные системы единиц, очень важно, чтобы инженеры умели уверенно переводить уравнения, за­ писанные в единицах «системы Дарси», в эквивалентную форму в единицах «промысловой» системы или любой другой системы еди­ ниц. Существует системный подход к выполнению таких преобразо­ ваний, который, при правильном применении, исключает возмож­ ность ошибки. Рассмотрим, например, перевод выражения (4.11) из «системы Дарси» в «промысловую» систему. Выражение

^ = и (см / с) х А (см2)

можно записать в форме, более подходящей для практического при­ менения, используя «систему Дарси»:

к (Д) А (см2)

<1р

 

(4.14)

ц(см3/с) = - ------ -----------— (атм/см)

р (сП)

61

 

 

или «промысловую» систему:

 

 

 

к (мД) А (фут2) бр

(фуиТф/д,’йм2/

(4.15)

^ (ст. барр. / сут) = - (соп81ап!)

61

р (сП)

 

 

В обоих выражениях использованы одни и те же обозначения па­ раметров.

Преобразование сводится к определению постоянной в уравнении (4.15). Для этого нужно помнить, что численные значения одних и тех же параметров в обоих уравнениях должны быть равны. То есть, если д имеет в уравнении (4.14) значение, например 200 см3/с, то ле­ вая часть уравнения (4.15) также должна иметь численное значение 200, хотя в последнем случае дебит выражен в ст. барр. / сут

^ (ст. барр. / сут) х [пересчетный коэффициент] = ^ (см3 / с)

или

^ (ст. барр. / сут) х [см3 / с / ст. барр. / сут] = ^ (см3 / с).

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

можно развернуть пересчетный коэффициент

 

Г, см3 / с

г, см3 / с

пл. барр. / сут

ст. барр. / сут

пл. барр. / сут

ст. барр. / сут

Применяя этот метод ко всем членам выражения, получаем (4.16)

ц (ст. барр. / сут)

п см3 / с

пл. барр. / сут

 

 

 

 

. пл. барр. / сут .

. ст. барр. / сут .

 

к , мД [Д / мД] А, фут2 [см2 / фут2]

 

 

 

атм

 

 

фр фунт / дюйм2

дина / см2

------------------------------------------х

41

фут

[см / фут]

Ц, сП

 

 

Поскольку [ д

]

1

см

= 30,48 и

атм

 

1

" 1000;

фут

фунт / дюйм

14,7’

мД

 

 

 

уравнение (4. "6) можно записать в виде

 

 

 

 

 

 

 

кА с!р

 

 

(4.17)

 

 

^ = - 1,127 х 10‘3 ------ - (ст. барр. / сут).

 

 

 

 

цВ„ 41

 

 

 

УПРАЖНЕНИЕ 4.1. ПЕРЕХОД ОТ ОДНОЙ СИСТЕМЫ ЕДИНИЦ К ДРУГОЙ

Требуется:

1)Определить пересчетный коэффициент для проницаемости при переходе от Д к см2 и м2.

2)Перевести полное уравнение прямолинейно-параллельного по­ тока несжимаемой жидкости, выраженное в единицах «системы Дарси» в «промысловую» систему.

кА /4р

р§

Чр 41 + 1,0133 х 106 41

УПРАЖНЕНИЕ 4.1. РЕШЕНИЕ

Для прямолинейно-параллельного горизонтального потока не­ сжимаемой жидкости