книги / Основы разработки нефтяных и газовых месторождений
..pdf4.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). Таким образом, видно, что параметр Ь в уравнении Дарси пря мо пропорционален разности потенциальных энергий на концах песчаной набивки.
Постоянная К / § применима только для потока воды, которая была единственной жидкостью в опытах Дарси. Исследования с различны ми жидкостями показали, что закон можно записать в общем виде
кр |
6Ф |
(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р |
р§ |
4г |
Чр 41 + 1,0133 х 106 41
УПРАЖНЕНИЕ 4.1. РЕШЕНИЕ
Для прямолинейно-параллельного горизонтального потока не сжимаемой жидкости