книги из ГПНТБ / Булыгин В.Я. Гидромеханика нефтяного пласта
.pdfИспользуя этп разложения, урав ѳнпе семейства изобар представим в виде
рт-\- В (Ѳ) рт~1-)-. . . = С,
П
ігде т = 2 2 = —2L по условию теоремы отлична от пуля. £=і
Дальнейшие рассуждения аналогичны приведенным в предыдущем дока зательстве.
«Точки-скважпны» для системы (II.6. Г)' являются особыми точками типа «узел». Интегральные кривые системы (II.6.1) — траектории п совпадающие •с нпмн линии тока в этих точках мало отличаются от пучка прямых с центром
в точке хі, у,-. Справедливость этого |
заключения вытекает пз того, |
что лпнип |
||||||||||||
'тока исходят пз источника п |
перпендикулярны изобарам, а последние вблизи |
|||||||||||||
|
|
|
|
скважины |
мало |
отличаются |
от окружностей. |
|||||||
|
|
|
|
Можно |
показать |
также, |
что точки покоя, где |
|||||||
|
|
|
|
др/дх = |
0, др/ду = 0, для системы (II.6.1) явля |
|||||||||
|
|
|
|
ются особыми точками типа «седло» и |
через этп |
|||||||||
|
|
|
|
точки |
проходят |
две |
или несколько |
эквипотен |
||||||
|
|
|
|
циален и интегральных кривых. |
|
|
||||||||
|
|
|
|
Линиями поворота первой и второй пазовем |
||||||||||
|
|
|
|
кривые, представляемые |
соответственно уравне |
|||||||||
|
|
|
|
ниями |
Fx (х, у) = |
О, |
|
(х, у) = |
0, где Fx и F„ |
|||||
|
|
|
|
заданы системой (II.6.3). |
первой |
(второй) линией |
||||||||
|
|
|
|
При пересечеішп |
с |
|||||||||
|
|
|
|
поворота |
линии |
тока, |
определяемые системой |
|||||||
|
|
|
|
(II.6.1), имеют касательную, параллельную осп |
||||||||||
Тис. 21. Линии тока |
при |
оу {ох), если только точка пересечения не яв |
||||||||||||
ляется особой точкой. В точке пересечения имеем: |
||||||||||||||
работе двух равнодебнтных |
||||||||||||||
эксплуатационных |
сква |
|
|
dx |
= 0; |
|
dy |
|
= 0. |
|
|
|||
жин. |
|
|
|
|
dt |
|
dt |
|
|
|
||||
-Линия поворота: |
1 — первая; |
|
|
|
|
|
|
|
|
|||||
Следовательно, dy/dx = |
сю. |
|
|
|
|
|||||||||
2 — вторая. |
|
|
|
|
|
|
||||||||
|
|
|
|
Затем через все «точки-скважпны» проходит |
||||||||||
|
|
|
|
первая |
(вторая) |
линия |
поворота и в этих точ |
|||||||
ках имеет касательную, параллельную |
осп оу {ох). |
Поместим начало коорди |
||||||||||||
нат в первую |
скважину (хг, у±) п уравнение F1 (х, у) = |
0 представим в сле |
||||||||||||
дующем виде: |
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
Fx {x, I/) = М |
(x?+ yf —2Ж£Х —2УіУ+ х2 + у2) + |
|
|
у2) В (х, у), |
||||||||||
2 |
(х2 + |
|||||||||||||
|
|
1=2 |
|
|
|
|
|
|
|
|
|
|
|
|
где В {х, у) — полином от х, |
|
|
|
|
|
|
|
|
|
|
||||
Отсюда имеем: |
|
|
дРг (0, 0) |
|
dFx (0, |
0) |
|
|
|
|
||||
|
Т і(0, |
0 )= 0 ; |
:0; |
= |
0, |
|
|
|||||||
ято и доказывает утверждение. |
дх |
|
|
ду |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||
Высказанные предложения могут помочь построить картину течения жидко |
сти в пластах горных пород, вскрытых скважинам. Для построения линий тока
.можно опираться на следующие положения.
1. Чтобы построить направление поля системы (II.6.1), данную точку М (х, у) «соединим со скважиной А х {х1Уу г) и по этой прямой отложим отрезок
|
|
_________ ____________ |
^«х |
|
|
|
Ѵ { х - Х 1)2 + { у - у 1)2 |
R l ’ |
|
где |
В 1 — расстояние |
от точки М до скважины, по направлению к скважине |
||
при |
стоке и от нее |
при источнике. Геометрическая |
сумма векторов М А г, |
|
М А 2, . . ., М А п дает |
направление поля, определяемого |
системой (II.6.1), что |
непосредственно следует из принципа суперпозиции.
60
2. Будем называть сепаратрпссами лпшш тока, проходящие через точки покоя, среди них особый интерес представляют те, которые ограничивают область притяжения к данной скважине. Через точку покоя, представляющую простое «седло», проходят две сепаратрпссы. Бслп примем точку покоя за начало, а направление вектора скорости за ось ох, то уравнение сепаратрпссы вблизи точки покоя может быть найдено в виде
у = Cjx2-]- С2х3,
причем коэффициенты определяются из дифференциального уравнения (II.6.4). Вдали от точки покоя уравнение строится методом, изложенным в (II.6.1).
3. Линии тока исходят из нагнетательных скважин и кончаются в эксплу атационных. Бесконечно удаленную точку принимают также за скважину с ха рактеристикой
п
£ = - 2 ^ . (=і
При пересечении линий поворота (первой п второй) линии тока меняют направление и имеют касательную, параллельную оси оу (ох).
Липпи тока не пересекают друг друга (и сепаратрисе), кроме как в особых точках. Приведем примеры.
Пример 1. Скважины расположены в точках (а, о), (— а, о) и имеют харак теристики Х± = Х2 = —1. Система (II.6-1) запишется в следующем виде:
|
dx |
|
X—а |
|
|
27-f- а |
|
|
|
|
dt |
|
(x -a )2 -\ -iß |
(х-)-а)2-|-у2 ’ |
|||||
|
d y |
= |
_______________У_________________________ У |
|
|
||||
откуда |
dt |
|
(х —а)'2-)-і/2 |
(x-j-a)a-|-y2 ’ |
|
||||
|
|
dy |
У(.r2 -fy 2 -)-a2 ) |
|
|
||||
|
|
|
|
|
|||||
|
|
|
dx |
X (х2 -{- 1/2— 0,2)' |
|
|
|||
Первая линия поворота состоит из оси оу и окружности |
|||||||||
|
|
|
|
а-2 + |
у2 = |
а 2 . |
|
|
|
Вторая линия поворота — из осп ох. Линии тока |
показаны на рис. 21. |
||||||||
Общий интеграл уравнения будет |
|
|
|
|
|
|
|||
|
|
|
х2 +У 3 + Сху = |
а2. |
|
|
|||
Пример 2. |
Скважина |
(а, о) — эксплуатационная |
с характеристикой — Аі, |
||||||
(—а, о) — нагнетательная с характеристикой Х2 (—а, о). |
В бесконечно удален |
||||||||
ной точке сток обильности Х2 — X.x. Система (а) имеет вид: |
|
||||||||
|
dx |
|
Хг (х — а) |
|
X2 (x + a ) |
|
|
||
|
dt |
|
(x — a)2+ |
A2 |
1 |
(x-j-a)2 + |
y2 |
> |
|
|
dy |
|
hV |
|
1 |
X2y |
2/2 |
|
|
откуда имеем |
dt |
|
(x —a)2+!/2 |
1 |
(X + O) 2 + |
> |
|||
|
|
|
|
|
|
|
|
|
|
dx |
—Xi (x — а) [(х + а)2 + |
і/2 ]4 |
-А2 У [(x —a)2- |
||||||
dt |
|
|
[(* -a )2 + p 2 ][(* + a )2 + 0 ® ] |
|
|
||||
dy |
ЯхУ[(а;+ а)2 + |
у2] + Я2Р[(а: —a)2 + y2]_ |
|||||||
dt |
|
[(я —а)2 + р2] [(® +0)2+у2] |
|
|
|||||
dy |
________ Хгу [(х — а) 2 |
+ р2 |
] _ А 2г/((х — д) 2 4 -У2]______ |
||||||
dx |
А,! (а —а) [(х + |
я.)2 |
-)-і/2 ] _ Х 2 (х —а) [(а —а)2 + і/2 ] |
61
Первая лшшя поворота будет: |
|
|
||
*■1 (*г —а) [(-г*+ |
«)3 + І/3] — ^2 (®+<0 [(а; —о)2 + |
1/2] = 0 , |
||
откуда |
а) |
|
у2 [Ѵ? (,г’ + а) — |
|
д2) [Хо (# |
Я-і (а;+ о)] + |
(а; —о)] = О, |
||
плп |
|
Ъ - х |
|
|
|
|
|
|
|
где |
|
6 + .г- (х'і- аі) , |
|
|
|
Ло “{- |
|
|
|
|
|
а = Ь. |
|
|
|
|
Х>2 — ^ |
|
|
Кривая состоит из |
двух ветвей: одна — замкнутая, пересекающая ось оаг |
|||
в точках (а, о) и (Ь, о), |
другая пересекает ось ох в точке х = —а и удаляется |
|||
в бесконечность при х = |
•—Ь. |
|
|
Рис. 23. Линии тока при работе двух равнодебптиых нагнетательных сква
скважины. жин. Обозначение см. па рис. 21.
Вторая линия поворота будет
h y [О® + а)2 + уЦ — Ъ2У [(а- — я)2 + т/2] = 0.
Она состоит из оси ох и кривой
[(а- + о)2 + 2/2]_Х27/ [(я —а)2 + і/2]=0.
Упрощая, получим:
|
|
(я — b)2 + J /2 = b 2 _ a2. |
Это |
окружность |
с центром в точке Ъ и радиуса І^Ь2 — а2. Изображение |
течения дано на рис. |
22. |
|
При |
= Х2 уравнение интегрируется в квадратурах, общий интеграл его |
*2 + (1/— С)2 = а2 + С2.
Интегральные кривые представлены на рис. 23.
Пример 3. Пусть три скважины расположены в точках (я, о), (—а, о), (о, о) и имеют равные характеристики — единицу. Первая линия поворота предста вляется следующим уравнением:
Я [Зу4 — 6я2і/2_|_3^4 — 4о2я2 -]- я4] — Q
62
II состоит из осп оу и двух замкнутых кривых, пересекающих ось ох в точках {± а , о), (± а /] /з , о). Сепаратрпссы, проходящие через точки покоя ( ± а / у 2 , о),
в первом приближении имеют следующее уравнение:
яЗ
1Г ’
Изображение течения дано па рис. 24.
Рис. 24. Линии тока при работе |
Рис. 25. Линли тока при работе двух |
||
трех эксплуатационных |
скважин. |
нагнетательных и одной эксплуатацион |
|
Обозначения см. на |
рис. |
21. |
ной скважины. |
|
|
|
Обозначения см. на рис. 21. |
Две эксплуатационные |
скважины |
с характеристиками ^ расположены |
в точках (о, а), (о, —а) и одна эксплуатационная скважина расположена в точке (о, о) и имеет характеристику —Х0. Первая линия поворота будет:
(2Х]. - |
%о) г/1 + [2h (2*2 + |
я2) - 2h (xi - Я2)] учд |
2h ^ (*2 + я2) - |
h |
(я2 + |
я2)2 = |
о. |
|||
При 2Л.Х— %0 |
0 — кривая замкнутая, |
пересекающая |
ось |
ох |
в точке |
|||||
(± а , |
о), а ось оу — в точке ( ±ß, |
о), где |
|
|
|
|
|
|||
|
„ |
а |
. |
во |
„о —(^о4-^і) + 1^?+4Ао?ч |
|
|
|
||
|
а - V ü ^ h ' Р |
|
2 h = T o |
• |
|
|
||||
Вторая линия поворота также замкнута и пересекает ось ох в точках ( ± у , |
о), |
|||||||||
где |
• |
|
|
|
_________ |
|
|
|
|
|
|
|
|
о _.о |
^-о + |
^і + 1^^і-Ь4Х0Х1 |
|
|
|
|
|
|
|
|
Ѵ ' |
|
2?ы-Хо |
|
|
|
|
|
Изображение течения дано на рис. 25.
Г л а в а III
СОВМЕСТНАЯ ФИЛЬТРАЦИЯ ДВУХ НЕСМЕШИВАЮЩИХСЯ ЖИДКОСТЕЙ
ВПЛАСТАХ ГОРНЫХ ПОРОД
Впроцессе разработки нефтяного месторождения извлечение нефти обычно сопровождается замещением нефти краевыми и нагне таемыми водами или газом. При этом по разрезу пласта нефть и вода (газ) продвигаются неравномерно, т. е. не образуют сплошного фронта вытеснения.
Например, на рис. 10, а наблюдается замещение нефти водой по всему разрезу пласта с образованием обводнеппых зон с высокой
Н | / |
2 |
Рис. 26. Расчленение фильтрации жидкости по мощности пласта.
1 — нефть; 2 — вода.
степенью вытеснения нефти, причем здесь не видим подъема перво начального положения ВНК, находящегося в нижней части пласта. На рис. 10, б, в происходит постепенное замещение нефти в нижних горизонтах, которые можно рассматривать как подъем ВНК. На рис. 10, г наблюдается послойное движение нефти и водонефтяной смеси. Еще более расчлененная по мощности картина движения жидкостей наблюдается на рис. 26 для екв. 183, 184, 54 Покровского нефтяного месторождения.
Такая же картина наблюдается и при экспериментальном изуче нии фильтрации жидкости на моделях нефтяных пластов, составлен ных из одного или нескольких пропластков, набранных средами различной проницаемости. Опыты по вытеснению жидкости, модели рующей нефть с водой, проводились многими исследователями,
64
в том числе сотрудниками лаборатории, Казанского университета. Приведем результаты опытов Б. И. Плещинского и др. [80, 81].
Для установления структуры течения был взят небольшой лоток с внутренними размерами 20 X 50 см, позволяющий набрать модель пласта мощностью до 2,2 см. Лоток изготовлялся из органического стекла и закрывался стеклянной крышкой. В качестве жидкости, моделирующей нефть, применялась смесь керосина и а-монобром- нафталина или масел с коэффициентом преломления света, равным коэффициенту преломления стекла. В качестве жидкости, модели рующей воду, бралась подкрашенная вода. Жидкость подводилась и отбиралась с помощью галерей, представляющих 10-мм трубки, рассверленные по стороне, обращенной в сторону рабочего поля модели, и прикрытые сеткой. Путем укладки стеклянной крошки разных размеров послойно создавалась одно-, двух- и трехслойная модель. Жидкостью модель заполнялась под вакуумом.
При проведении экспериментов без остаточной нефти модель целиком заполнялась жидкостью, моделирующей нефть. Определя лась проницаемость модели и проводился опыт по вытеснению. Для проведения экспериментов при наличии остаточной воды модель первоначально заполнялась дистиллированной водой. Замерялась проницаемость. Затем прокачивалась нефть до прекращения выхода воды и проводился эксперимент по вытеснению нефти. Во всех экспе риментах регистрировалось количество вышедшей жидкости, фото графировалась картина вытеснения.
На рис. 27 изображены положения границы раздела воды и нефти. Для трехслойных моделей (см. рис. 27, б, в) нижний слой виден только на правых частях рисунка, а верхний — на левых.
Сделаны следующие выводы.
1.Распределение скоростей движения частично на границе жидко стей близко к нормальному (гауссовскому).
2.При отсутствии остаточной воды и отношениях вязкостей
р,н/р,в = |
1 |
-г 2 (см. |
рис. 27, а) ВНК движется сплошным фронтом. |
||
3. При Цн/р-в = |
3-^-4 |
фронт вытеснения рвется (теряет устой |
|||
чивость), |
|
образуя |
целые |
необводненные участки (см. рис. 27, б). |
|
При рн/ра |
> 4 |
вода |
движется отдельными струйками (см. рис. 27, в— |
||
е). Для |
рн/рв |
= 2 — 3 (в |
промежутках интервала) опыт мог быть |
отнесен к тому или к другому случаю.
4. При отсутствии остаточной воды вытесняющая вода быстро продвигается по более проницаемому слою и из него проникает в менее проницаемые слои, образуя локальные, изолированные участки заводнения, которые постепенно растут, сливаются и обра зуют зону заводнения.
5. При наличии остаточной воды движение происходит более равномерно, даже при значительном отношении (в опытах до 6) вязкостей (см. рис. 27, г). В слоистой модели фронт вытеснения движется также более равномерно и его движение по малопроница емым слоям мало отстает от движения по более проницаемому (27,9). При небольших скоростях вытеснения и наличии остаточной воды
5 Заказ 322 |
65 |
|
скорости движ ения воды по слоям размой проницаемости вы равни
ваются (см. |
рис. 27, е). |
П еред фронтом |
движ ения |
нагнетаемой |
п од |
||
краш енной |
воды |
часто |
наблю дается |
движ ение |
вала |
погребенной |
|
прозрачной |
воды. |
|
|
|
|
|
|
Перечисленные |
обстоятельства движ ения границы |
раздела |
воды |
и нефти ' необходим о учитывать при построении математической модели вы теснения нефти краевой водой.
Рис. 27. Результаты экспериментов на прозрачных моделях по вытеснению сверху, нижние — вид
Модели: а, г, 0, е — однослойные; б,
66
Результаты промысловых и опытных данных говорят о том, что картина фильтрации меняется при изменении условий. Для расчетов движения нефти и воды в пластах горных пород в настоящее время
пользуются несколькими схемами или моделями, |
описывающими |
|||||||||||||||||
это движение. |
В зависимости |
от |
структуры потоков, характерных |
|||||||||||||||
особенностей строения |
нефтяных |
коллекторов, |
свойств жидкостей и |
|||||||||||||||
целей |
расчетов |
эти модели должны |
меняться |
иногда |
в |
пределах |
||||||||||||
одной расчетной |
площади. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Построение |
универсальной |
модели, |
описывающей все случаи |
|||||||||||||||
|
|
|
|
|
|
|
движения, |
едва ли целесооб |
||||||||||
|
|
|
|
|
|
|
разно ввиду усложнений, ко |
|||||||||||
|
|
|
|
|
|
|
торые вносятся |
вместе с уче |
||||||||||
|
|
|
|
|
|
|
том |
дополнительных |
факто |
|||||||||
|
|
|
|
|
|
|
ров. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Можно пойти другим пу |
||||||||||
|
|
|
|
|
|
|
тем: разработать унифициро |
|||||||||||
|
|
|
|
|
|
|
ванный ряд моделей, в ко |
|||||||||||
|
|
|
|
|
|
|
тором можно вести счет по |
|||||||||||
|
|
|
|
|
|
|
одной |
из отработанных схем, |
||||||||||
|
|
|
|
|
|
|
в |
пределах |
целого |
пласта |
||||||||
|
|
|
|
|
|
|
(месторождения) |
|
или |
его ча |
||||||||
|
|
|
|
|
|
|
сти, |
и переходить |
к |
другой |
||||||||
|
|
|
|
|
|
|
при |
изменении контрольных |
||||||||||
|
|
|
|
|
|
|
условий. |
|
Переход |
|
должен |
|||||||
|
|
|
|
|
|
|
быть |
|
выполним |
|
без |
суще |
||||||
|
|
|
|
|
|
|
ственного |
нарушения |
всей |
|||||||||
|
|
|
|
|
|
|
модели и только путем за |
|||||||||||
|
|
|
|
|
|
|
мены |
|
ее |
части |
(блока), |
при |
||||||
|
|
|
|
|
|
|
этом |
следует |
иметь |
возмож |
||||||||
|
|
|
|
|
|
|
ность использовать |
основной |
||||||||||
|
|
|
|
|
|
|
массив параметров |
и |
исход |
|||||||||
|
|
|
|
|
|
|
ных |
данных. |
Для |
выполне |
||||||||
|
|
|
|
|
|
|
ния |
последнего |
условия |
не |
||||||||
|
|
|
|
|
|
|
обходимо, |
чтобы методы из |
||||||||||
|
|
|
|
|
|
|
мерения |
параметров и исход |
||||||||||
|
|
|
|
|
|
|
ных |
данных |
были |
унифици |
||||||||
|
|
|
|
|
|
|
рованы |
в |
данном |
|
случае, |
|||||||
|
|
|
|
|
|
|
чтобы |
входящие |
величины |
|||||||||
|
|
|
|
|
|
|
измерялись |
одинаковым |
ме |
|||||||||
|
|
|
|
|
|
|
тодом. |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
Такой ряд вместе с зада |
||||||||||
|
|
|
|
|
|
|
нием |
|
признаков |
|
перехода |
|||||||
|
|
|
|
|
|
|
может играть |
роль |
универ |
|||||||||
|
|
|
|
|
|
|
сальной |
схемы |
|
счета |
при |
|||||||
|
|
|
|
|
|
|
большой конкретизации и уп |
|||||||||||
керосина, |
водой |
(левые |
рисунки — впд |
рощениях внутри каждой мо |
||||||||||||||
снизу). |
|
|
|
|
|
|
дели ряда. Все это существен |
|||||||||||
я — трехелойиые. |
|
|
|
|
|
но |
|
упрощает |
составление |
67
программ для ЭВМ, давая возможность придать им блочную структуру и комплектовать программы из отработанных блоков.
В настоящей главе предлагается ряд из 12 моделей, который удовлетворяет сформулированным требованиям. Классификация проведена в зависимости от структуры (слоистости) пластов и от вида функции Баклея — Леверетта (см. § 2).
Некоторые моделщ ряда уже использовались ранее, однако здесь они приведены к единым параметрам и значительно расширены условия их применимости. Не будем касаться истории развития этой
темы, |
отсылая читателя |
к работе |
[126]. |
Для |
написания |
§ 1, 5, 6 |
||||||
использованы работы |
[26, |
53, |
78, |
142, |
143, 147], § 3 — [137, |
138, |
||||||
145], |
§ 4 |
- [1 2 |
-1 4 , |
43, |
44, |
105, |
112], |
§ |
10 - |
[37, 38, 41], |
§ |
12 - |
[1, 2, |
26, |
86, 87, |
90, 120, |
149] |
и др. |
|
|
|
|
|
Желающих ознакомиться с более общими в определенном смысле моделями отсылаем к работам [70, 50]. Отметим, что при исполь зовании этих моделей объем вычислений весьма велик, что заставляет искать пути его уменьшения, например упрощая схемы.
Мы не считаем, что следует ограничиться только предложенными в этой главе моделями. Предлагаемые модели являются простей шими, дающими возможность в некоторой степени учесть особенности движения двухжидкостных систем в пластах горных пород, при нимая во внимание геологические особенности их строения и условия разработки. Схемы упрощаются от § 1 до § 8, а затем усложняются снова вследствие учета дополнительных факторов.
§ 1. Схема трехмерной фильтрации с пропиткой
Рассмотрим движение двухжидкостных систем.
Считается, что в каждой точке пласта существует два давления — нефти рі и воды рв. Разность этих давлений рн — pt = рк объяс няется наличием капиллярных сил и называется капиллярным скачком.
Запишем выражения для скоростей фильтрации:
Vix = |
— |
kxk*i |
dp\ |
.* |
|
|
|
|
Pi |
|
dx |
|
|
|
|||
Vl y = |
- |
kyk* |
dPi |
(і = н, |
B) |
(III.1.1) |
||
Pi |
dy |
’ |
|
|||||
Viz = |
— |
k7k* |
( |
dp* |
|
■ |
|
|
Pi |
1 |
dz |
|
f p i g ) |
|
Равенство і = н, в означает, что і принимает значения н или в. Введем приведенное пластовое давление pt для нефти и воды:
. . |
Рі = Р*і + Pigz. |
(III.1.2) |
68
Тогда, подставляя р* в (III.1.1), получаем
|
dpi |
|
I |
|
Рі |
дх |
|
|
|
k y k t |
dpi |
А-тА-г* V pp, |
(IIU .3) |
|
Pi |
Оу |
|||
Рі |
|
|||
A'zA'i |
др; |
|
|
|
Pi |
dz |
|
|
Для капиллярного скачка будем иметь следующее выражение:
|
|
Рк = Р н ~ Р„ — gz(pH— рв). |
|
|
|
|
(III.1.4) |
||||
Уравнение неразрывности и состояния запишутся: |
|
|
|||||||||
~ |
(РPu) + ^ |
{ ? Ы « ) + 4 г |
|
+ |
aiPötm) |
= °; |
|
(ІП -1-5) |
|||
|
р = |
р(р); т — т (р); |
к* = к*(р)-, |
|
|
|
(III.1.6) |
||||
|
К = кЛр)'> ку=ку(р); кг =Ісг(рУ, |
|
|
(III.1.7) |
|||||||
|
|
|
s„ + sB= |
i. |
|
|
|
|
|
(III.1.8) |
|
Подставим теперь |
в |
(III.1.5) |
значения |
скоростей фильтрации |
|||||||
из (III.1.3) и продифференцируем последний член. Тогда |
получим: |
||||||||||
д |
|
|
к у к І _ д р і _ _ \ , ± _ |
( |
_ k ^ l |
dPi \ |
|||||
дх |
Рі |
|
Pi |
|
ду ) ~ г |
dz |
\ |
vt |
Ц,. |
dz |
J |
|
|
|
dm |
1 |
dpi |
|
|
dsi |
|
|
(III.1.9) |
|
|
|
dpi |
J |
dt |
9im |
|
dt |
|
|
Разложим функции, определенные формулами (III.1.6)—(III.1.8),
в ряд Тейлора и оценим члены, входящие в (III.1.5), |
аналогично |
|||||||||||
гл. 2, § 2. Тогда получим: |
|
|
|
|
|
|
|
|
||||
д_ / |
кхк? |
dpi |
\ |
д |
/ |
кук* |
dpt \ |
. д |
/ |
кгЩ |
dpi |
\ _ |
дх \ |
р(- |
дх |
) |
‘ ду |
\ |
р,- |
diI ) |
‘ dz |
\ |
ці |
dz |
) |
|
|
= |
Si (nzß* + |
ßc) |
|
4 r |
|
• |
|
(III.1.10) |
Будем считать, что капиллярный скачок зависит только от водо насыщенности. Тогда может быть задана функция
PK= PK(S); S = sB. |
(III.1.11) |
Известными функциями водонасыщенности являются также вели чины относительных проницаемостей. Используя выражение
А* = **(«), |
(III.1.12) |
из (III. 1.4) получаем: |
|
Рн — Рв = Рк (s) — gz (р„ — р„). |
(III.1.13) |
69