1464
.pdf§ 2. Плоское радиальное движение по линейному закону
Исследуем горизонтальный плоский радиальный поток несжимае мой жидкости в условиях водонапорного режима. Будем считать, что жидкость притекает к гидродинамически совершенной скважине по линейному закону фильтрации, пласт однороден и его мощность постоянна и равна Ь(см. рис. 45).
На рис. 54 приведено изобра жение в плане исследуемого пото ка: Ас — горизонтальное сечение сква жины (индекс с всегда будет соответ ствовать первой букве слова «сква жина»); Ак — сечение контура об ласти питания; Ас и Ак — концен тричные окружности, радиусы кото рых равны Rc и RK. Допустим, что первоначально во всем пласте и на за бое скважины приведенное давление одинаково и равно р£. Величину р* можем назвать начальным статиче ским приведенным пластовым давле нием. Пусть в какой-то момент вре мени из скважины начали отбирать жидкость, благодаря чему должно понизиться давление на забое скважины. Конечно, давление в пласте по низится и всюду вокруг скважины,
но мы будем считать, что на границе области питания все вре мя поддерживается постоянное начальное давление р*. Приведен ное давление на забое скважины обозначим через р* и будем назы вать его приведенным динамическим забойным давлением. Итак, при ток жидкости к скважине обеспечивается за счет перепада давления в пласте (р* —р*).
Заметим, что величину р* можно и не называть именно начальным статическим приведенным давлением, а называть просто статическим приведенным пластовым давлением. Действительно, по условиям за дачи жидкость считается несжимаемой, величина р* постоянна, а по тому после остановки скважины в ней самой и во всем пласте вновь восстанавливается давление р*.
Пусть круговые сечения Ак и Ас на рис. 54 проведены на уровне нижней горизонтальной границы (подошвы) пласта. Давления вдоль
окружностей Ас и Ак пусть также будут известными и равными соот ветственно рс и рк. Ясно, что
Рк - Рс = Рк - Рс = ЛР, |
(16, IX) |
где через Ар обозначено понижение (перепад) давления в скважине. Из ранее сделанных замечаний вытекает, что вполне достаточно
изучить исследуемый поток лишь в одной плоскости, например в плос кости, изображенной на рис. 54.
Задача исследования потока состоит в том, чтобы определить дебит скважины, давление, градиент давления и скорость филь трации в любой точке пласта, а также установить закон движе ния частиц жидкости вдоль их траекторий. Все перечисленные неиз вестные должны быть выражены через заданные величины рК1 рс (или р*, р*), Ь, Дк, Дс, /*, А;, т .
Согласно линейному закону фильтрации в форме (9, VUE) и равен
ству (11, VIII), получим: |
|
+ - Й 4 ' |
(17' К) |
где р — давление в некоторой точке М пласта, радиус-вектор которой равен г, <2 — постоянный дебит скважины (расход жидкости через каж дую цилиндрическую поверхность, концентричную скважине, равен ее дебиту Q).
Для определения давления в точке М проинтегрируем уравнение (17, IX):
Р к |
Фм |
Я к |
dr |
|
[ |
f |
(18, IX) |
||
Jdp=2*bkJ |
r • |
|
||
p |
|
r |
|
|
откуда |
|
|
|
|
P PK |
QL1 I |
RK |
(19, IX) |
|
2n bk |
r ' |
|
||
где In — натуральный логарифм. |
|
проинтегрируем |
уравнение |
|
Для определения дебита |
скважины |
|||
в других пределах: |
|
Як |
|
|
|
Qf1 |
|
|
|
|
Г dr |
|
||
|
2ттbk J |
г ’ |
|
|
Рс |
|
яс |
|
|
откуда1*
Q = 2тг b k (p K - pc) |
(21, IX) |
Подставив найденное значение Q из уравнения (21, IX) в (19, IX) и (9, VIII), найдем соответственно градиент давления, давление и ско рость фильтрации в любой точке плоского сечения исследуемого пото ка:
d p |
_ Р к |
Р с |
1 |
(22, |
IX) |
d r |
~ Ь |
|
Г’ |
||
|
|
|
|
||
|
1пА |
|
|
|
|
Р = Р к |
Р к - |
Р с |
ь ф , |
(23, |
IX) |
|
1 |
-*ч< |
|
|
|
1пж
г
|
1_ |
|
1п “5" |
|
itc |
с? II |
‘З4 |
|
Ц |
1П£-, itc
1
Г’
(23', IX)
(24, IX)
Заметим, что в формуле (23, IX) отношение натуральных логариф мов можно было бы заменить отношением десятичных логарифмов.
Аналогично тому, что |
было сказано в предыдущем параграфе, |
и здесь формулы (22, IX) |
— (24, IX) определяют искомые величины |
не только в одной плоскости потока, но и во всем потоке; согласно ра венству (16, IX) во всех только что выведенных формулах разность давлений можно заменить разностью приведенных давлений.
Прежде чем перейти к анализу выведенных формул, выполним
некоторые преобразования; положим, что |
|
|
|
|||
|
Рк — *уНк, |
Рс — 7-^с, |
(25, |
IX) |
||
|
Р к - Р с = Р * к - |
Рс |
= 7 (# к - |
Н с ) = 7 5 . |
(26, |
IX) |
1 Формулу (21, |
IX) можно было |
бы |
вывести |
из формулы |
(19, IX), заметив, |
|
что р = рс при г = |
Rc . |
|
|
|
|
|
Так как в формулы типа (21, IX), (24, IX) входят разности дав лений, то под величинами рк и рс мы можем подразумевать и абсо лютные и избыточные (над атмосферным) давления. Соответственно с этим под величинами Нк и Нс следует понимать пьезометрические высоты, соответствующие либо абсолютным, либо избыточным стати ческому и динамическому давлениям на забое скважины; s — пониже ние пьезометрического уровня жидкости в скважине (s = Нк - Нс).
Если скважина насосная и жидкость в ней однородная с весом единицы объема 7 , то под Нк и Нс можно подразумевать фактические высоты динамического и статического уровней над забоем скважины, а под s — фактическое понижение динамического уровня под статиче ским.
Однако нефтяные скважины часто бывают обводнены, причем удельный вес минерализованной воды всегда больше удельного веса нефти и различен для различных пластов. Кроме этого, столбы неф ти в скважинах бывают газированы. Удельные же веса газированной и обводненной нефти в скважинах трудно определить. Поэтому часто предпочитают в формулах (25, IX) и (26, IX) считать величину 7 соот ветствующей химически чистой воде в нормальных условиях, т. е.
7 = 0,001 кг/см3 |
(27, IX) |
В таком случае величины Нк и Нс представляют собой высоты статического и динамического столбов жидкости в скважине, пересчи танные на воду.
Для точного определения пьзометрических высот статического и динамического уровней наиболее желательно определить глубинным манометром статическое и динамическое давления на забое скважины; тогда, подставляя величину 7 из уравнения (27, IX) в (25, IX) — (26, IX), легко подсчитать Нк, Яс, а значит, и 5 в сантиметрах водяного столба.
Последний метод одинаково применим для фонтанных, компрессо рах и насосных скважин.
Учитывая равенство (26, IX), формулу (21, IX) можно переписать
так: |
|
_ 27гЬ/ь7 (# к - # с) _ |
27rbkrys |
, RK |
. RK |
Используя формулу (28, IX), надо не забывать оговаривать, чему равен вес единицы объема жидкости 7 (как только что было упомянуто, он может отличаться от веса единицы объема жидкости в пластовых условиях), что и определяет физический смысл величин Як, # с,
Перейдем к анализу выведенных формул. Судя по формулам (21, IX) и (28, IX) дебит скважины оказывается прямо пропорциональ ным перепаду давления в ней или понижению пьезометрического уров ня.
График зависимости дебита сква жины от перепада давления (или от по нижения пьезометрического уровня) на зывается индикаторной диаграммой. Из сказанного ясно, что в рассматривае мом случае индикаторной линией будет прямая линия. Такая индикаторная диа грамма изображена на рис. 55, на оси абсцисс откладывается дебит скважи ны, а на оси ординат — перепад давле ния, или понижение пьезометрического уровня; ось ординат удобнее направлять вниз, ибо тогда понижению вдоль оси ординат наглядно соответствует пони жение уровня (снижение забойного дав ления) в скважине.
Пользуясь обычными в расчетах подземной гидравлики размерностями, т. е. измеряя к в д, р в сантипуазах, b в см, перепад давления Ар в am,
мы из формулы (21, IX) получим дебит скважины в см3/сек (поскольку в эту формулу входит отношение радиусов Дк и i?c, постольку радиусы могут иметь любую, но только одинаковую размерность).
Желая получить дебит в м3/сутки, измеряя мощность пласта в м и сохраняя для остальных величин прежние размерности, перепишем формулу (21, IX) так:
Q = 23,6 |
kbAp |
(29, IX) |
к
/•tig Rc
причем для удобства расчетов в последней формуле совершен переход к десятичным логарифмам.
Пример 1. Допустим, например, что к = 1 д,
р = 1 сп, Ар = 1 am, b = 10 м, RQ = 10 теле, Rc = 10 см
(диаметр скважины приблизительно равен 8 дюймам). Тогда из формулы (21, IX) найдем:
Q = 2 -3,14 -1000-1 1 = 546 см3/сек
1 In 105
Подстановка тех же величин в формулу (29, IX) дает:
Q = 23> |
= 47,2 м3/сутки |
Другие подсчеты при тех же исходных данных будут приведены дальше.
Если при исследовании скважины замерены ее дебит и перепад давления, известны мощность пласта, вязкость нефти в пластовых условиях, радиус скважины и величина RKможет быть, примерно, оце нена то по любой из формул (21, IX), (28, IX), (29, IX) можно опреде лить коэффициент проницаемости пласта к. Такой метод определения к по формуле (21, IX) носит название метода определения проницаемости по промысловым данным2.
Заметим, что те же формулы употребляются для определения про ницаемости керна в лабораторных условиях в тех случаях, когда вдоль оси керна просверлено отверстие (подобие скважины) и осуществлены условия радиальной фильтрации жидкости через керн.
Рис. 56. График зависимости градиента давления и скорости фильтрации от расстояния до центра скважины.
Как показывают формулы (22, IX) и (24, IX), градиент давления и скорость фильтрации в любой точке пласта обратно пропорциональ-
2Метод определения проницаемости по кернам весьма ценен, но он не может за менить метода определения проницаемости по промысловым данным; в последнем случае на основании формулы (21, IX) (с соответствующими поправками на гидро динамическое несовершенство скважины — см. последующие главы), определяется средняя эффективная проницаемость пласта вокруг скважины, что очень важно.
11 Подземная гидравлика
ны расстоянию от этой точки до оси скважины. На рис. 56 на оси ор динат откладываются величины скорости фильтрации либо градиента давления, а на оси абсцисс — расстояния г от рассматриваемой точ ки пласта до оси скважины; пунктирная линия А В , проходящая на расстоянии Rc от оси ординат* соответствует положению стенки сква жины. Начерченная на рис. 56 равнобочная гипербола В С , асимпто тами которой служат оси координат, представляет собой построенный по формулам (22, IX) и (24, IX) график зависимости градиента давле ния и скорости фильтрации от радиуса-вектора. Из графика ясно вид но, что при приближении к скважине градиенты давления и скорость фильтрации резко возрастают, причем своего наибольшего значения их величины достигают у стенки скважины.
Сформулированный сейчас вывод совершенно очевиден. Действи тельно, траектории частиц жидкости радиально сходятся к оси сква жины; одно и то же количество жидкости должно проходить через боковые цилиндрические поверхности, размеры которых убывают про порционально радиусу. Ясно, что в этих условиях скорости фильтрации должны изменяться именно обратно пропорционально радиусу цилин дрической поверхности, т. е. расстоянию до оси скважины:
Q |
(30, IX) |
|
2пгЪ' |
||
|
Рис. 57. Логарифмические пьезометрические линии АВ и CD.
Из формулы (23, IX) следует, что давление в какой-либо точке пла ста есть логарифмическая функция расстояния от этой точки до оси
скважины. Построенные на рис. 57 по формуле (23, IX) логарифми ческие кривые линии АВ и CD изображают пьезометрические линии; параллельно оси z откладываются отрезки, пропорциональные истин ному (удобнее избыточному, чем абсолютному) пластовому давлению3. Вращая линии АВ или CD вокруг оси скважины Oz, получим так на зываемую «пьезометрическую воронку депрессии». Точки А и С ле жат на стенке скважины, так что АС = 2Дс. Линия DEFB определя ет положение статического пьезометрического уровня в пласте, а ли ния СА — положение пьезометрического динамического уровня жид кости в скважине. Величина отрезка СЕ определяет понижение s пье зометрического уровня жидкости в скважине [см. формулу (26, IX}]. В условиях рассматриваемого потока пьезометрические (депрессионные) линии АВ и CD отнюдь не касаются линии D EFB , а пересекают ее в точках D и В под некоторым углом. Если бы в каких-либо точках пьезометрических линий касательные были горизонтальны, т. е. гради ент давления был бы равен нулю, то отсюда следовало бы, что в этих точках и скорость фильтрации была равна нулю. Однако в рассматри ваемых условиях такой случай невозможен, ибо при постоянном деби те скважины, при отсутствии инфильтрации в пласт через его кров лю и подошву и в предположении несжимаемости жидкости и самого пласта жидкости в любой точке пласта должна двигаться в направле нии к скважине; не может существовать никакой ограниченной зоны влияния (или зоны дренирования) скважины, вне которой жидкость якобы остается в покое, а внутри движется к скважине. В упомяну тых условиях предположение о существовании ограниченного «радиуса влияния» скважины является логически противоречивым и физически бессмысленным. Влияние скважины должно распространяться на весь пласт — до его естественных границ, т. е. в данном случае до границы области питания. Расстояние DB между точками пересечения пьезо метрических линий АВ и CD с линией DEFB равно удвоенному ра диусу контура области питания Дк. Пьезометрические линии (сечения пьезометрической воронки) имеют большой уклон около стенки сква жины и сильно выполаживаются по мере удаления от нее. На рис. 56 это не столь заметно, ибо Rc для наглядности взято непропорциональ но большим. Наличие уклона пьезометрических линий в точках D и В свидетельствует о притоке жидкости из области питания в пласт.
Радиусы Rc и Дк входят во все выведенные выше формулы под знаком логарифма, а потому в рассматриваемых условиях их величины
3Если ось абсцисс поместить на уровне опорной плоскости, то параллельно оси ординат можно откладывать отрезки, пропорциональные приведенному пластовому давлению, т. е. отрезки, пропорциональные или равные напорам в соответствующих точках пласта.
мало влияют на подсчеты дебита скважины, давлений в различных точках пласта и т.д.; об этом более подробно сказано в главе XIV.
Для определения семейства изо бар заметим, что, судя по форму ле (23, IX), давление одинаково в тех точках плоскости движения, в кото рых
г = const = с. |
(31, IX) |
Следовательно, уравнение (31, IX) представляет собой уравнение се мейства изобар; изобарами служат окружности, концентричные сече нию скважины Ас (рис. 58).
Ясно, что и здесь изобары орто гональны к траекториям, совпадаю щим с радиусами окружностей.
Указанное в § 2 главы VIII пра вило построения поля (карты) изо бар будет в данном случае выполне
но, если величины радиусов следующих друг за другом изобар будут изменяться в геометрической профессии. Действительно, при измене нии входящей под знаком логарифма в формулу (23, IX) величины г
вгеометрической прогрессии величина давления р будет изменяться
варифметической прогрессии, что и требуется (см. дальше табл. 5). Для вычерченных на рис. 58 изобар отношение радиусов любых
двух соседних окружностей равно 1,5.
Удовлетворяя упомянутому правилу, изобары располагаются тем теснее, чем ближе они к скважине. Зная, что разность давлений меж ду двумя любыми соседними изобарами всюду одинакова, мы, глядя на рис. 58, сразу замечаем, что давление быстрее всего изменяется в бли жайшей окрестности к скважине. Итак, рис. 58 является правильно построенным гидродинамическим полем плоско-радиального потока.
Понятно, что поверхностями равных приведенных давлений или поверхностями равных напоров будут служить боковые поверхности цилиндров, соосных скважине; окружности на рис. 58 можно рассмат ривать как горизонтальные сечения поверхностей равных напоров.
В зависимости от того, на каком уровне проведены горизонтальные сечения цилиндрических поверхностей равных напоров, одни и те же изобары будут характеризовать распределение либо истинных, либо приведенных пластовых давлений.
Допустим, что на расстоянии г от работающей (возмущающей) скважины находится бездействующая (реагирующая) скважина M NT (рис. 57); последнюю скважину можно рассматривать как естествен ный пьезометр, облегчающий наблюдения за распределением пласто вого давления. С помощью формулы (23, IX) легко определить соот ношение между понижениями пьезометрического уровня S H S P B воз мущающей и реагирующей скважинах. Действительно, из упомянутой формулы и из формулы (26, IX) следует:
Р к |
- Р |
_ £ р |
(32, IX) |
|
Р к |
- Pc |
s |
||
R, |
Подчеркнем, что sp = T N — понижение пьезометрического уровня в реагирующей скважине, вызванное отбором жидкости из единствен ной возмущающей скважины при тех условиях ее работы, какие были сформулированы в начале данного параграфа.
Для более ясного представления о распределении пластовых дав лений рассмотрим числовые примеры.
Пример 2. Пусть RK= 10 км = 106 см, Rc = 10 см. Тогда на осно вании формулы (32, IX) можно составить следующую табл. 5, в кото рой приведены значения относительного понижения пьезометрических уровней на разных расстояниях от возмущающей скважины (рис. 57).
Та б л и ц а 5
Относительные понижения пьезометрических уровней в пласте на разных расстояниях от возмущающей скважины; таблица рассчитана по формуле (32, IX)
г, м 0,1 |
1 |
5 |
10 |
50 |
100 |
500 |
1000 |
2000 |
5000 |
10000 |
Sp |
0,80 |
0,66 |
0,60 |
0,46 |
0,40 |
0,26 |
0,20 |
0,14 |
0,06 |
0 |
1 |
||||||||||
S |
|
|
|
|
|
|
|
|
|
|
Из табл. 5 видно, что при рассматриваемых условиях на расстоя нии 1 м от оси скважины теряется 20% перепада давления (т. е. потеря такая же, как на пути в 9 км от г = 10 км до г = 1 тш); на рассто янии 10 м теряется 40% и т.д. Таким образом, действительно, наи большие потери давления (а следовательно, и наибольшие градиенты
давления) имеют место вблизи скважины. Из той же таблицы видно,
5р что при изменении радиуса г в 10 раз величина отношения — меняется
в арифметической прогрессии с разностью 0,20.