Skip to content
Brent$117,6Газ TTF$1 192USD453,87 ₸
MunaiHub
Физика пластаВытеснение нефти Средний 19 мин чтения

Функция доли воды в потоке и соотношение подвижностей

Доля воды в потоке (fractional flow of water, fwf_w) — доля воды в суммарном расходе нефти и воды через сечение пласта в пластовых условиях: fw=qw/(qw+qo)f_w = q_w/(q_w + q_o). Она зависит от водонасыщенности, кривых относительной фазовой проницаемости (ОФП) и вязкостей фаз, а в наклонном пласте — ещё от силы тяжести, темпа закачки и капиллярного давления. Кривая fw(Sw)f_w(S_w) — главный инструмент расчёта вытеснения нефти водой; в отечественной литературе её называют также функцией Баклея–Леверетта или функцией распределения потоков фаз.

Зачем это нужно

По кривой доли воды видно, при какой насыщенности пласта вода начинает преобладать в потоке, как быстро растёт обводнённость и что меняют вязкая нефть, гидрофобная порода, наклон пласта или другой темп закачки. На этой кривой построен расчёт фронта вытеснения и прорыва воды (Теория Баклея–Леверетта и метод Уэлджа), а в гидродинамической модели она переносит насыщенность от ячейки к ячейке (Уравнения двухфазной фильтрации).

Кратко о главном

  • fwf_w — доля воды в потоке, а не в порах. При одной и той же насыщенности вода может давать и 20 %, и 85 % расхода: всё решает отношение подвижностей krw/μwk_{rw}/\mu_w и kro/μok_{ro}/\mu_o.
  • В горизонтальном пласте без капиллярных сил fw=λw/λtf_w = \lambda_w/\lambda_t и от темпа закачки не зависит. Сила тяжести и капиллярное давление добавляют поправки, обратно пропорциональные суммарной скорости фильтрации utu_t.
  • Чем вязче нефть, тем левее кривая: вода начинает преобладать в потоке при меньшей насыщенности, и большую часть нефти добывают при высокой обводнённости. Половина потока — вода там, где krw/kro=μw/μok_{rw}/k_{ro} = \mu_w/\mu_o.
  • В гидрофобной породе вода течёт по крупным порам, её ОФП высока, и кривая сдвигается влево — вытеснение хуже, чем в гидрофильной.
  • Вода тяжелее нефти: при движении вверх по восстанию сила тяжести её придерживает, и кривая сдвигается вправо, при движении вниз по падению — влево. Меру эффекта даёт гравитационное число NgN_g; оно растёт при снижении темпа.
  • Капиллярный член в гидрофильной породе увеличивает fwf_w в зоне фронта и размывает фронт на длину порядка λoΔpc/ut\lambda_o\Delta p_c/u_t: в пласте это от десятков сантиметров до метров, в опыте на керне — величина, сравнимая с длиной образца.
Ключевые соотношения
fw=qwqw+qo=1+kkroμout(∂pc∂x−Δρ gsin⁡α)1+krokrw⋅μwμof_w = \dfrac{q_w}{q_w + q_o} = \dfrac{1 + \dfrac{k k_{ro}}{\mu_o u_t}\left(\dfrac{\partial p_c}{\partial x} - \Delta\rho\,g\sin\alpha\right)}{1 + \dfrac{k_{ro}}{k_{rw}}\cdot\dfrac{\mu_w}{\mu_o}}

общая форма для одномерного потока; α>0\alpha > 0, если вода движется вверх по восстанию

fw=λwλt(1−Ng krokro′),Ng=k kro′ Δρ gsin⁡αμoutf_w = \dfrac{\lambda_w}{\lambda_t}\left(1 - N_g\,\dfrac{k_{ro}}{k_{ro}'}\right),\qquad N_g = \dfrac{k\,k_{ro}'\,\Delta\rho\,g\sin\alpha}{\mu_o u_t}

без капиллярного члена: вязкая часть и поправка на силу тяжести

fw=0,5 ⟺ krwkro=μwμof_w = 0{,}5\ \Longleftrightarrow\ \dfrac{k_{rw}}{k_{ro}} = \dfrac{\mu_w}{\mu_o}

горизонтальный пласт: воды и нефти течёт поровну

ℓc=k kro Δpcμout\ell_c = \dfrac{k\,k_{ro}\,\Delta p_c}{\mu_o u_t}

капиллярная длина — масштаб размытия фронта

Где чаще всего ошибаются
  • Путать долю воды в пластовом потоке с обводнённостью продукции на поверхности: их связывают объёмные коэффициенты нефти и воды.
  • Считать, что темп закачки меняет кривую fwf_w горизонтального пласта: без силы тяжести и капиллярных сил он сокращается.
  • Переносить вывод «сила тяжести почти не влияет» с одномерного потока вдоль падения на мощный пласт: поперёк пласта вода и нефть разделяются, и одномерная модель перестаёт работать.
  • Сравнивать кривые разных пород только по водонасыщенности, не глядя на их подвижные интервалы 1−Swi−Sor1 - S_{wi} - S_{or}.
  • Судить о том, сколько нефти осталось у скважины, по одной обводнённости: та же обводнённость при вязкой нефти означает гораздо больше оставшейся нефти.

После урока вы сможете

  • вывести уравнение доли воды в потоке из закона Дарси для двух фаз с силой тяжести и капиллярным давлением и объяснить смысл каждого члена;
  • построить кривую fw(Sw)f_w(S_w) по кривым ОФП и вязкостям и связать её форму с соотношением подвижностей;
  • оценить, как вязкость нефти, смачиваемость породы, наклон пласта и темп закачки сдвигают кривую и что это значит для вытеснения;
  • рассчитать гравитационное число и капиллярную длину и решить, можно ли пренебречь силой тяжести и капиллярными силами;
  • пользоваться кривой как инструментом: переводить долю воды в обводнённость, находить насыщенность у скважины по её обводнённости, сравнивать варианты разработки.

Что показывает доля воды в потоке

Там, где вода вытесняет нефть, в порах движутся обе фазы. Каждая течёт под действием градиента своего давления со своей подвижностью λ=k kr/μ\lambda = k\,k_r/\mu: проницаемостью породы, умноженной на относительную проницаемость фазы и делённой на её вязкость (Относительная проницаемость). Поэтому расход делится между фазами не пропорционально их насыщенностям, а пропорционально подвижностям.

Доля воды в потоке fw=qw/(qw+qo)f_w = q_w/(q_w + q_o) — это и есть результат такого деления в данной точке пласта. Расходы берут в пластовых условиях. На выходе из пласта, у добывающей скважины, fwf_w определяет состав её продукции; внутри пласта — то, как быстро вода переносится вперёд и как движется фронт вытеснения.

Три величины, которые часто путают
Величина
Водонасыщенность SwS_w
Что показывает
долю объёма пор, занятую водой
От чего зависит
истории насыщения: сколько воды пришло и сколько нефти ушло
Величина
Доля воды в потоке fwf_w
Что показывает
долю воды в пластовом расходе через сечение
От чего зависит
SwS_w, кривых ОФП, вязкостей; в наклонном пласте — ещё от темпа и наклона
Величина
Обводнённость продукции
Что показывает
долю воды в добытой жидкости в стандартных условиях
От чего зависит
fwf_w на забое и объёмных коэффициентов BoB_o, BwB_w

Обводнённость (water cut) — доля воды в продукции; её связь с fwf_w выведена в уроке об ОФП.

Уравнение доли воды с силой тяжести и капиллярным давлением

Рассмотрим одномерный поток вдоль оси xx от нагнетательной стороны к добывающей; ось наклонена к горизонту под углом α\alpha. Флюиды и порода несжимаемы, поэтому суммарная скорость фильтрации ut=uo+uw=qt/Au_t = u_o + u_w = q_t/A одинакова в любом сечении. Насыщенность считаем одинаковой по сечению — так называемый диффузный (рассредоточенный) поток. Для каждой фазы записываем закон Дарси с её давлением и весом (двухфазный закон Дарси), выражаем градиенты давлений, вычитаем одно уравнение из другого и пользуемся тем, что разность давлений фаз — капиллярное давление pc=po−pwp_c = p_o - p_w. Выкладки приведены в разделе для углублённого изучения; результат такой:

Доля воды в потоке с силой тяжести и капиллярным давлением
fw=1+k kroμo ut(∂pc∂x−Δρ gsin⁡α)1+krokrw⋅μwμof_w = \frac{1 + \dfrac{k\,k_{ro}}{\mu_o\,u_t}\left(\dfrac{\partial p_c}{\partial x} - \Delta\rho\,g\sin\alpha\right)}{1 + \dfrac{k_{ro}}{k_{rw}}\cdot\dfrac{\mu_w}{\mu_o}}
где:
  • fwf_wдоля воды в суммарном потоке в пластовых условиях, —
  • kkабсолютная проницаемость породы, м²
  • kro, krwk_{ro},\ k_{rw}ОФП нефти и воды при текущей водонасыщенности, —
  • μo, μw\mu_o,\ \mu_wвязкости нефти и воды в пластовых условиях, Па·с
  • utu_tсуммарная скорость фильтрации qt/Aq_t/A — расход обеих фаз на единицу площади сечения, м/с
  • pcp_cкапиллярное давление po−pwp_o - p_w, Па
  • xxкоордината вдоль потока, от нагнетательной стороны к добывающей, м
  • Δρ\Delta\rhoразность плотностей воды и нефти ρw−ρo\rho_w - \rho_o в пластовых условиях, кг/м³
  • ggускорение свободного падения, 9,81, м/с²
  • α\alphaугол наклона оси xx к горизонту: положительный, если вода движется вверх по восстанию, отрицательный — вниз по падению, рад

Единицы: СИ (любая согласованная система единиц)

Справедливо, если
  • одномерный линейный поток, насыщенность одинакова по сечению (диффузный поток);
  • порода и флюиды несжимаемы, поэтому utu_t не меняется вдоль потока; свободного газа нет;
  • фазы несмешивающиеся, ОФП и капиллярное давление зависят только от насыщенности и направления её изменения (при заводнении гидрофильного пласта — кривые пропитки);
  • закон Дарси выполняется для каждой фазы: течение ламинарное.

Знаменатель — вязкая часть: как поток делился бы в горизонтальном пласте без капиллярных сил. Числитель содержит две поправки. Капиллярная зависит от того, как насыщенность меняется вдоль потока; гравитационная — от наклона и разности плотностей. Обе поправки делятся на utu_t, поэтому они тем сильнее, чем медленнее закачка.

В других системах единиц (1)
Метрические промысловые единицы
fw=1+0,0853 k kro Aqt μo(∂pc∂x−0,00981 Δρsin⁡α)1+krokrw⋅μwμof_w = \frac{1 + \dfrac{0{,}0853\,k\,k_{ro}\,A}{q_t\,\mu_o}\left(\dfrac{\partial p_c}{\partial x} - 0{,}00981\,\Delta\rho\sin\alpha\right)}{1 + \dfrac{k_{ro}}{k_{rw}}\cdot\dfrac{\mu_w}{\mu_o}}

kk — мД, AA — м², qtq_t — м³/сут в пластовых условиях, μ\mu — мПа·с, ∂pc/∂x\partial p_c/\partial x — МПа/м, Δρ\Delta\rho — г/см³.

0,08530{,}0853 — постоянная закона Дарси в метрических промысловых единицах (Проницаемость горных пород); 0,00981=9,81⋅1000/1060{,}00981 = 9{,}81\cdot1000/10^6 переводит г/см³ в кг/м³ и Па в МПа, как в уроке Капиллярное давление и переходная зона.

Ту же формулу удобно записать через подвижности λo=k kro/μo\lambda_o = k\,k_{ro}/\mu_o, λw=k krw/μw\lambda_w = k\,k_{rw}/\mu_w и суммарную подвижность λt=λo+λw\lambda_t = \lambda_o + \lambda_w. Тогда видно, что вязкая часть — это отношение λw/λt\lambda_w/\lambda_t, а поправки умножаются на множитель λoλw/λt\lambda_o\lambda_w/\lambda_t, который равен нулю там, где подвижна только одна фаза:

Доля воды в потоке через подвижности фаз
fw=λwλt+λoλwλt⋅1ut(∂pc∂x−Δρ gsin⁡α)f_w = \frac{\lambda_w}{\lambda_t} + \frac{\lambda_o\lambda_w}{\lambda_t}\cdot\frac{1}{u_t}\left(\frac{\partial p_c}{\partial x} - \Delta\rho\,g\sin\alpha\right)
где:
  • λo, λw\lambda_o,\ \lambda_wподвижности нефти и воды k kro/μok\,k_{ro}/\mu_o и k krw/μwk\,k_{rw}/\mu_w, м²/(Па·с)
  • λt\lambda_tсуммарная подвижность λo+λw\lambda_o + \lambda_w, м²/(Па·с)
  • ut, pc, Δρ, g, αu_t,\ p_c,\ \Delta\rho,\ g,\ \alphaкак в общей форме выше, СИ

Единицы: СИ

Справедливо, если
  • те же, что для общей формы.

Это то же выражение, что поток воды через суммарную скорость в уроке Уравнения двухфазной фильтрации, делённый на utu_t: там направление силы тяжести задано градиентом глубины, здесь — углом наклона оси потока. В гидродинамическом симуляторе доля воды без поправок, λw/λt\lambda_w/\lambda_t, переносит насыщенность, а поправки считаются отдельными членами.

Вязкость нефти и воды: форма кривой

В горизонтальном пласте без капиллярных сил остаётся одна вязкая часть fw=λw/λtf_w = \lambda_w/\lambda_t — знаменатель общей формулы. Кривая зависит только от отношения ОФП (свойство породы) и отношения вязкостей (свойство флюидов). Вязкая нефть течёт хуже, поэтому при той же насыщенности большую часть потока даёт вода, и кривая сдвигается влево и вверх.

Удобная точка отсчёта — насыщенность, при которой воды и нефти течёт поровну. Из формулы fw=0,5f_w = 0{,}5, когда krw/kro=μw/μok_{rw}/k_{ro} = \mu_w/\mu_o: отношение ОФП должно скомпенсировать отношение вязкостей. Если нефть в 60 раз вязче воды, достаточно, чтобы ОФП воды была в 60 раз меньше ОФП нефти, — а это случается вскоре после того, как вода стала подвижной.

Доля воды в потоке при трёх вязкостях нефти0,20,30,40,50,60,70,80,00,20,40,60,81,0Водонасыщенность SwДоля воды в потоке fw30 мПа·с3 мПа·с1,2 мПа·с

График: по горизонтали водонасыщенность от 0,2 до 0,8, по вертикали доля воды в потоке от 0 до 1. Три кривые начинаются с нуля при водонасыщенности 0,22 и выходят на единицу к 0,75. Самая левая и крутая — для нефти вязкостью 30 мПа·с, она достигает 0,5 при водонасыщенности около 0,37. Средняя — для 3 мПа·с, 0,5 около 0,47. Правая — для 1,2 мПа·с, 0,5 около 0,52; в начале она дольше других остаётся близкой к нулю.

РасчётДоля воды в потоке для кривых ОФП урока Относительная проницаемость (гидрофильная порода) при вязкости воды 0,5 мПа·с и трёх вязкостях нефти. Точки — насыщенности, при которых воды и нефти течёт поровну.
Расчётный примерОдна порода, три нефти
Исходные данные (условные)
Кривые ОФП (урок Относительная проницаемость)
Swi=0,22S_{wi} = 0{,}22, Sor=0,25S_{or} = 0{,}25, kro′=0,85k_{ro}' = 0{,}85, krw′=0,25k_{rw}' = 0{,}25, no=2,5n_o = 2{,}5, nw=3n_w = 3
Вязкость воды μw\mu_w
0,5 мПа·с
Вязкость нефти μo\mu_o
1,2; 3,0 и 30 мПа·с
Водонасыщенность для сравнения
Sw=0,45S_w = 0{,}45
Решение
  1. Соотношение подвижностей в концевых точках M′=(0,25/0,5)/(0,85/μo)M' = (0{,}25/0{,}5)/(0{,}85/\mu_o): 0,71; 1,76 и 17,6.
  2. ОФП при Sw=0,45S_w = 0{,}45: Sn=(0,45−0,22)/0,53=0,4340S_n = (0{,}45 - 0{,}22)/0{,}53 = 0{,}4340; krw=0,25⋅0,43403=0,02044k_{rw} = 0{,}25\cdot0{,}4340^3 = 0{,}02044; kro=0,85⋅(1−0,4340)2,5=0,2049k_{ro} = 0{,}85\cdot(1 - 0{,}4340)^{2{,}5} = 0{,}2049; отношение kro/krw=0,2049/0,02044=10,02k_{ro}/k_{rw} = 0{,}2049/0{,}02044 = 10{,}02.
  3. Доля воды fw=1/(1+10,02⋅μw/μo)f_w = 1/(1 + 10{,}02\cdot\mu_w/\mu_o): для 1,2 мПа·с — 1/(1+10,02⋅0,4167)=0,1931/(1 + 10{,}02\cdot0{,}4167) = 0{,}193; для 3,0 — 1/(1+10,02⋅0,1667)=0,3741/(1 + 10{,}02\cdot0{,}1667) = 0{,}374; для 30 — 1/(1+10,02⋅0,0167)=0,8571/(1 + 10{,}02\cdot0{,}0167) = 0{,}857.
  4. Насыщенность, при которой fw=0,5f_w = 0{,}5, то есть krw/kro=μw/μok_{rw}/k_{ro} = \mu_w/\mu_o, находят по кривым подбором: 0,518; 0,474 и 0,373. Проверка для 3,0 мПа·с: при SwS_w = 0,474 krw=0,0275k_{rw} = 0{,}0275, kro=0,1663k_{ro} = 0{,}1663, их отношение 0,165 — с точностью округления это 0,5/3,0≈0,1670{,}5/3{,}0 \approx 0{,}167.
  5. Доля воды 0,98 — уровень, близкий к пределу рентабельной добычи, — достигается при SwS_w 0,670; 0,643 и 0,550.
Ответ: При Sw=0,45S_w = 0{,}45 доля воды в потоке — 19 %, 37 % и 86 % для нефти вязкостью 1,2, 3,0 и 30 мПа·с; половина потока становится водой при SwS_w 0,518, 0,474 и 0,373.

Что это значит. Порода одна и та же, а картина разная. При маловязкой нефти (M′<1M' < 1) вода преобладает в потоке лишь после того, как вытеснено больше половины подвижной нефти: точка fw=0,5f_w = 0{,}5 лежит на 56 % подвижного интервала. При нефти 30 мПа·с — уже на 29 %: дальше нефть идёт вместе с большим количеством воды, а доля 0,98 наступает задолго до насыщенности 1−Sor=0,751 - S_{or} = 0{,}75. Отсюда правило: чем больше M′M', тем большую часть нефти добывают после прорыва воды и тем больше воды приходится прокачать. Данные примера условные.

Поздний подъём
при очень маловязкой нефти или загущённой воде, когда M′M' много меньше единицы: доля воды почти до конца интервала остаётся малой, а затем резко растёт. Вытеснение близко к поршневому — почти вся подвижная нефть уходит до прихода воды (поршневое вытеснение).
S-образная кривая
типичный случай, к нему относятся все три кривые рисунка: сначала доля воды растёт всё быстрее, после точки перегиба — всё медленнее приближается к единице. На фронте образуется скачок насыщенности, за ним вода и нефть текут вместе.
Ранний подъём
чем вязче нефть, тем ближе перегиб к SwiS_{wi} (у кривой для 30 мПа·с — около 0,36): доля воды круто растёт почти сразу после того, как вода стала подвижной. Вода приходит к скважинам рано, и основную нефть добывают при высокой обводнённости.

Смачиваемость породы

Смачиваемость решает, какие поры занимает вода. В гидрофильной породе вода покрывает зёрна плёнкой и занимает мелкие поры, а нефть — крупные; закачанная вода долго движется по узким путям вдоль зёрен, и её ОФП остаётся низкой. В гидрофобной породе всё наоборот: вода — несмачивающая фаза, она идёт по крупным порам, её ОФП быстро растёт, а нефть остаётся в мелких порах и плёнках на зёрнах. Признаки смачиваемости по кривым ОФП сведены в правила, приведённые в уроке об ОФП.

Кривые ОФП двух пород для сравнения
Параметр
SwiS_{wi}; SorS_{or}
Гидрофильная (урок об ОФП)
0,22; 0,25
Гидрофобная (условно)
0,12; 0,20
Параметр
kro′k_{ro}'; krw′k_{rw}'
Гидрофильная (урок об ОФП)
0,85; 0,25
Гидрофобная (условно)
0,80; 0,60
Параметр
non_o; nwn_w
Гидрофильная (урок об ОФП)
2,5; 3
Гидрофобная (условно)
3,5; 1,8
Параметр
Пересечение кривых ОФП
Гидрофильная (урок об ОФП)
Sw≈0,56S_w \approx 0{,}56
Гидрофобная (условно)
Sw≈0,40S_w \approx 0{,}40

Параметры гидрофобной породы приняты условно в соответствии с правилами Крейга: малая неснижаемая водонасыщенность, пересечение кривых левее 0,5 и высокая концевая ОФП воды.

Доля воды в потоке в гидрофильной и гидрофобной породе0,10,20,30,40,50,60,70,80,90,00,20,40,60,81,0Водонасыщенность SwДоля воды в потоке fw0,160,85гидрофобнаягидрофильная

График: по горизонтали водонасыщенность от 0,1 до 0,9, по вертикали доля воды в потоке от 0 до 1. Кривая гидрофобной породы начинается при 0,12 и поднимается круто: при водонасыщенности 0,29 доля воды уже 0,5, при 0,40 — 0,85. Кривая гидрофильной породы начинается при 0,22, медленно поднимается и достигает 0,5 при 0,47; при 0,40 она равна 0,16. Обе точки при 0,40 отмечены.

РасчётДоля воды в потоке в гидрофильной и гидрофобной породе при одних и тех же флюидах (нефть 3,0 мПа·с, вода 0,5 мПа·с). Точки — значения при Sw=0,40S_w = 0{,}40 из примера ниже.
Расчётный примерОдна и та же водонасыщенность в двух породах
Исходные данные (условные)
Гидрофильная порода
кривые урока Относительная проницаемость
Гидрофобная порода (условно)
Swi=0,12S_{wi} = 0{,}12, Sor=0,20S_{or} = 0{,}20, kro′=0,80k_{ro}' = 0{,}80, krw′=0,60k_{rw}' = 0{,}60, no=3,5n_o = 3{,}5, nw=1,8n_w = 1{,}8
Флюиды
нефть 3,0 мПа·с, вода 0,5 мПа·с
Водонасыщенность
Sw=0,40S_w = 0{,}40
Решение
  1. Гидрофильная: Sn=(0,40−0,22)/0,53=0,3396S_n = (0{,}40 - 0{,}22)/0{,}53 = 0{,}3396; krw=0,25⋅0,33963=0,00979k_{rw} = 0{,}25\cdot0{,}3396^3 = 0{,}00979, kro=0,85⋅(1−0,3396)2,5=0,301k_{ro} = 0{,}85\cdot(1 - 0{,}3396)^{2{,}5} = 0{,}301. Подвижности на единицу проницаемости: λw=0,00979/0,5=0,01958\lambda_w = 0{,}00979/0{,}5 = 0{,}01958, λo=0,301/3,0=0,10033\lambda_o = 0{,}301/3{,}0 = 0{,}10033 1/(мПа·с); fw=0,01958/(0,01958+0,10033)=0,163f_w = 0{,}01958/(0{,}01958 + 0{,}10033) = 0{,}163.
  2. Гидрофобная: Sn=(0,40−0,12)/0,68=0,4118S_n = (0{,}40 - 0{,}12)/0{,}68 = 0{,}4118; krw=0,60⋅0,41181,8=0,1215k_{rw} = 0{,}60\cdot0{,}4118^{1{,}8} = 0{,}1215, kro=0,80⋅(1−0,4118)3,5=0,1249k_{ro} = 0{,}80\cdot(1 - 0{,}4118)^{3{,}5} = 0{,}1249; λw=0,1215/0,5=0,2430\lambda_w = 0{,}1215/0{,}5 = 0{,}2430, λo=0,1249/3,0=0,04163\lambda_o = 0{,}1249/3{,}0 = 0{,}04163; fw=0,2430/(0,2430+0,04163)=0,854f_w = 0{,}2430/(0{,}2430 + 0{,}04163) = 0{,}854.
  3. Половина потока — вода при SwS_w = 0,474 в гидрофильной породе и 0,289 в гидрофобной.
Ответ: При Sw=0,40S_w = 0{,}40 вода даёт 16 % потока в гидрофильной породе и 85 % — в гидрофобной.

Что это значит. В гидрофобной породе вода становится основной фазой потока почти сразу: она идёт по крупным порам, обходя нефть. Поэтому при заводнении такого пласта вода прорывается раньше, а нефть добывается дольше и с большим количеством воды. Расчёт фронта по этим кривым (Теория Баклея–Леверетта и метод Уэлджа) даёт коэффициент вытеснения к прорыву 0,36 против 0,51 в гидрофильной породе. Это не обязательно значит, что гидрофобный пласт в итоге отдаст меньше нефти: в условном примере его подвижный интервал даже шире (0,68 против 0,53), но эту нефть он отдаёт медленно, при большой прокачке воды. Большинство пластов — смешанной или промежуточной смачиваемости, и их кривые лежат между этими двумя. Данные примера условные.

Наклон пласта и сила тяжести

Вода тяжелее нефти. Когда она движется вверх по восстанию — например, при закачке в пониженную часть залежи или за контур нефтеносности, — её собственный вес тормозит её, а нефть, более лёгкая, получает добавочный «подъём». Вода отстаёт, доля воды в потоке при той же насыщенности меньше, и вытеснение лучше. При движении вниз по падению сила тяжести, наоборот, подгоняет воду. В формуле этому соответствует член −Δρ gsin⁡α-\Delta\rho\,g\sin\alpha. Если пренебречь капиллярным членом, его удобно собрать в одно безразмерное число.

Гравитационное число
Ng=k kro′ Δρ gsin⁡αμo utN_g = \frac{k\,k_{ro}'\,\Delta\rho\,g\sin\alpha}{\mu_o\,u_t}
где:
  • NgN_gгравитационное число: отношение градиента от разности весов фаз вдоль пласта, Δρ gsin⁡α\Delta\rho\,g\sin\alpha, к градиенту давления, который нужен, чтобы прокачать через породу одну нефть со скоростью utu_t, —
  • kro′k_{ro}'концевая ОФП нефти при неснижаемой водонасыщенности, —
  • k, Δρ, g, α, μo, utk,\ \Delta\rho,\ g,\ \alpha,\ \mu_o,\ u_tкак в уравнении доли воды, СИ

Единицы: СИ (безразмерное число)

Справедливо, если
  • диффузный одномерный поток вдоль падения: насыщенность одинакова по толщине;
  • знак: Ng>0N_g > 0, когда вода движется вверх по восстанию.

При NgN_g порядка сотых сила тяжести на вытеснение почти не влияет; при NgN_g порядка десятых и больше её надо учитывать. Число растёт с проницаемостью, разностью плотностей и углом и обратно пропорционально скорости фильтрации и вязкости нефти.

В других системах единиц (2)
Метрические промысловые единицы
Ng=0,0853 k kro′ A⋅0,00981 Δρsin⁡αqt μoN_g = \frac{0{,}0853\,k\,k_{ro}'\,A\cdot 0{,}00981\,\Delta\rho\sin\alpha}{q_t\,\mu_o}

kk — мД, AA — м², Δρ\Delta\rho — г/см³, qtq_t — м³/сут в пластовых условиях, μo\mu_o — мПа·с.

0,00981 Δρsin⁡α0{,}00981\,\Delta\rho\sin\alpha — градиент давления вдоль пласта от разности плотностей, МПа/м; 0,08530{,}0853 — постоянная закона Дарси (Проницаемость горных пород).

Промысловые единицы (oilfield)
Ng=4,886⋅10−4 k kro′ A Δγsin⁡αqt μoN_g = 4{,}886\cdot10^{-4}\,\frac{k\,k_{ro}'\,A\,\Delta\gamma\sin\alpha}{q_t\,\mu_o}

kk — мД, AA — ft², Δγ\Delta\gamma — разность относительных плотностей воды и нефти (пресная вода = 1), qtq_t — пластовые баррели в сутки, μo\mu_o — сП.

4,886⋅10−4=1,127⋅10−3⋅0,43354{,}886\cdot10^{-4} = 1{,}127\cdot10^{-3}\cdot0{,}4335: постоянная закона Дарси в промысловых единицах, умноженная на градиент давления столба пресной воды в psi/ft.

Доля воды в потоке в наклонном пласте без капиллярного члена
fw=λwλt(1−Ng krokro′)f_w = \frac{\lambda_w}{\lambda_t}\left(1 - N_g\,\frac{k_{ro}}{k_{ro}'}\right)
где:
  • λw/λt\lambda_w/\lambda_tвязкая часть — доля воды в горизонтальном пласте, —
  • kro/kro′k_{ro}/k_{ro}'ОФП нефти, отнесённая к концевой точке: 1 при SwiS_{wi}, 0 при 1−Sor1 - S_{or}, —

Единицы: безразмерная

Справедливо, если
  • капиллярным членом пренебрегают;
  • NgN_g — по формуле выше, с тем же знаком наклона.

Поправка пропорциональна kro/kro′k_{ro}/k_{ro}': в относительном выражении она больше всего там, где нефть ещё хорошо подвижна, — при малых и средних водонасыщенностях, то есть в зоне фронта. При 1−Sor1 - S_{or} нефть неподвижна, и fw=1f_w = 1 при любом наклоне.

Расчётный примерГравитационное число наклонного пласта
Исходные данные (условные)
Проницаемость kk
400 мД
Сечение потока
толщина 25 м, ширина 600 м
Угол падения пласта
12°, закачка в пониженную часть: вода идёт вверх по восстанию
Плотности в пласте
вода 1,04 г/см³, нефть 0,80 г/см³ (как в уроке Капиллярное давление и переходная зона)
Нефть
3,0 мПа·с; кривые ОФП урока Относительная проницаемость, kro′=0,85k_{ro}' = 0{,}85
Суммарный расход в пластовых условиях
600 м³/сут; для сравнения — 150 м³/сут
Решение
  1. Площадь сечения A=25⋅600=15 000A = 25\cdot600 = 15\,000 м². Разность плотностей Δρ=1,04−0,80=0,24\Delta\rho = 1{,}04 - 0{,}80 = 0{,}24 г/см³.
  2. Градиент от силы тяжести вдоль пласта: 0,00981⋅0,24⋅sin⁡12°=0,00981⋅0,24⋅0,2079=4,89⋅10−40{,}00981\cdot0{,}24\cdot\sin 12° = 0{,}00981\cdot0{,}24\cdot0{,}2079 = 4{,}89\cdot10^{-4} МПа/м.
  3. Гравитационное число: Ng=0,0853⋅400⋅0,85⋅15 000⋅4,89⋅10−4600⋅3,0=0,118N_g = \dfrac{0{,}0853\cdot400\cdot0{,}85\cdot15\,000\cdot4{,}89\cdot10^{-4}}{600\cdot3{,}0} = 0{,}118. При расходе 150 м³/сут оно в четыре раза больше: 0,472.
  4. При Sw=0,40S_w = 0{,}40 вязкая часть λw/λt=0,163\lambda_w/\lambda_t = 0{,}163 (из примера о смачиваемости), kro/kro′=0,301/0,85=0,354k_{ro}/k_{ro}' = 0{,}301/0{,}85 = 0{,}354.
  5. Вверх по восстанию при 600 м³/сут: fw=0,163⋅(1−0,118⋅0,354)=0,156f_w = 0{,}163\cdot(1 - 0{,}118\cdot0{,}354) = 0{,}156; при 150 м³/сут: 0,163⋅(1−0,472⋅0,354)=0,1360{,}163\cdot(1 - 0{,}472\cdot0{,}354) = 0{,}136.
  6. Если бы ту же воду закачивали в повышенную часть и она шла вниз по падению (150 м³/сут): 0,163⋅(1+0,472⋅0,354)=0,1900{,}163\cdot(1 + 0{,}472\cdot0{,}354) = 0{,}190.
Ответ: Ng≈N_g \approx 0,118 при 600 м³/сут и 0,472 при 150 м³/сут; при Sw=0,40S_w = 0{,}40 сила тяжести меняет долю воды от 0,163 до 0,136 (вверх по восстанию) или 0,190 (вниз по падению).

Что это значит. Даже в проницаемом крутом пласте при обычном темпе сила тяжести вдоль пласта меняет долю воды лишь на несколько процентов; при вчетверо меньшем темпе — уже на 17 %. Расчёт фронта показывает, что при закачке вверх по восстанию с расходом 150 м³/сут коэффициент вытеснения к прорыву растёт с 0,508 до 0,524. В тонких пластах с проницаемостью в десятки миллидарси и пологим залеганием NgN_g не превышает сотых долей, и наклоном в одномерном расчёте обычно пренебрегают. Данные примера условные.

Доля воды в потоке в горизонтальном и наклонном пласте0,20,30,40,50,60,70,80,00,20,40,60,81,0Водонасыщенность SwДоля воды в потоке fwвниз по падениювверх по восстанию

График: по горизонтали водонасыщенность от 0,2 до 0,8, по вертикали доля воды в потоке от 0 до 1. Три близкие S-образные кривые. Сплошная — горизонтальный пласт. Штриховая кривая для движения вниз по падению проходит левее и выше, для движения вверх по восстанию — правее и ниже; наибольшее расхождение, около 0,045 по доле воды в каждую сторону, — при водонасыщенности около 0,45–0,52. Выше 0,65 все три кривые сливаются у единицы.

РасчётДоля воды в потоке для нефти 3,0 мПа·с в горизонтальном пласте (сплошная линия) и в пласте с наклоном 12° при расходе 150 м³/сут из примера (NgN_g около 0,47): вода идёт вверх по восстанию или вниз по падению (штриховые).

Если Ng kro/kro′>1N_g\,k_{ro}/k_{ro}' > 1, формула даёт fw<0f_w < 0: при очень малом темпе вверх по восстанию вода стекает вниз против общего потока, а нефть всплывает быстрее него. Это противоток (counter-current flow) — признак того, что сила тяжести сильнее вязких сил и вытеснение идёт в режиме гравитационного разделения фаз.

Темп закачки

Темп закачки входит в уравнение только через utu_t в знаменателях капиллярной и гравитационной поправок. В горизонтальном пласте без капиллярных сил скорость сокращается, и кривая fw(Sw)f_w(S_w) от темпа не зависит: быстрая и медленная закачка проходят одни и те же состояния, только за разное время. Темп важен тогда, когда работают другие силы.

  • Вверх по восстанию — меньший темп усиливает стабилизирующее действие силы тяжести: кривая сдвигается вправо, вытеснение ближе к поршневому.
  • Вниз по падению — меньший темп, наоборот, позволяет воде «скатываться» вперёд: кривая сдвигается влево.
  • Капиллярные силы — при меньшем темпе фронт размывается на большую длину (раздел ниже); в пласте это обычно несущественно, в опыте на керне — существенно.
  • Неоднородный пласт — при меньшем темпе вода успевает перетекать между прослоями под действием капиллярных сил и силы тяжести; это уже за пределами одномерной кривой (Вытеснение в неоднородных пластах).

Капиллярный член: размытый фронт

Капиллярная поправка в числителе общей формулы пропорциональна градиенту капиллярного давления вдоль потока. Его удобно разложить на два множителя:

Градиент капиллярного давления вдоль потока
∂pc∂x= dpc dSw⋅∂Sw∂x\frac{\partial p_c}{\partial x} = \frac{\dd p_c}{\dd S_w}\cdot\frac{\partial S_w}{\partial x}
где:
  •  dpc/ dSw\dd p_c/\dd S_wнаклон кривой капиллярного давления (пропитки при заводнении), Па
  • ∂Sw/∂x\partial S_w/\partial xизменение водонасыщенности вдоль потока, 1/м

Единицы: СИ

Справедливо, если
  • капиллярное давление — функция только водонасыщенности.

В гидрофильной породе капиллярное давление тем выше, чем меньше воды:  dpc/ dSw<0\dd p_c/\dd S_w < 0. На фронте вытеснения водонасыщенность убывает по ходу потока: ∂Sw/∂x<0\partial S_w/\partial x < 0. Произведение положительно, и капиллярный член увеличивает долю воды в зоне фронта: порода «впитывает» воду впереди фронта (капиллярная пропитка), и фронт размывается.

Обычно этим членом пренебрегают по простой причине: градиент насыщенности вдоль потока заранее неизвестен — его и ищут при расчёте вытеснения. Насколько это оправдано, показывает капиллярная длина — расстояние, на котором капиллярная поправка сравнима с единицей.

Капиллярная длина — масштаб размытия фронта
ℓc=k kro Δpcμo ut\ell_c = \frac{k\,k_{ro}\,\Delta p_c}{\mu_o\,u_t}
где:
  • ℓc\ell_cдлина, на которой капиллярный член в уравнении доли воды порядка единицы, м
  • Δpc\Delta p_cперепад капиллярного давления на фронте: между водонасыщенностью перед фронтом и за ним, Па
  • k, kro, μo, utk,\ k_{ro},\ \mu_o,\ u_tкак в уравнении доли воды; krok_{ro} — характерная ОФП нефти в зоне фронта, СИ

Единицы: СИ

Справедливо, если
  • оценка порядка величины: если зона фронта длиннее ℓc\ell_c, капиллярный член мал по сравнению с вязкой частью;
  • с концевой ОФП kro′k_{ro}' вместо средней по фронту оценка получается завышенной.

Фронт вытеснения размыт капиллярными силами на длину порядка ℓc\ell_c; эту размытую область фронта называют стабилизированной зоной. Если ℓc\ell_c много меньше расстояния между скважинами, фронт можно считать скачком.

В других системах единиц (1)
Метрические промысловые единицы
ℓc [м]=0,0853 k kro Δpcμo ut\ell_c\,[\text{м}] = \frac{0{,}0853\,k\,k_{ro}\,\Delta p_c}{\mu_o\,u_t}

kk — мД, Δpc\Delta p_c — МПа, μo\mu_o — мПа·с, utu_t — м/сут.

0,08530{,}0853 — постоянная закона Дарси (Проницаемость горных пород): q=0,0853 kA Δp/(μL)q = 0{,}0853\,kA\,\Delta p/(\mu L), откуда L=0,0853 k Δp/(μu)L = 0{,}0853\,k\,\Delta p/(\mu u).

Расчётный примерПласт и керн: насколько размыт фронт
Исходные данные (условные)
Порода
50 мД, кривые ОФП урока Относительная проницаемость, kro′=0,85k_{ro}' = 0{,}85; нефть 3,0 мПа·с
Капиллярное давление пропитки
условная кривая урока Уравнения двухфазной фильтрации: pc=0,04 (1−Sn)2p_c = 0{,}04\,(1 - S_n)^2 МПа
Насыщенности по обе стороны фронта
Swi=0,22S_{wi} = 0{,}22 и насыщенность на фронте 0,564 (Теория Баклея–Леверетта и метод Уэлджа)
Скорость фильтрации в пласте
закачка элемента 255 м³/сут через сечение 300 × 12 м: 0,0708 м/сут
Опыт на керне
образец длиной 7 см, скорость фильтрации 3,0 м/сут
Решение
  1. Перепад капиллярного давления на фронте: pc(Swi)=0,040p_c(S_{wi}) = 0{,}040 МПа, pc(0,564)=0,0049p_c(0{,}564) = 0{,}0049 МПа, Δpc=0,0351\Delta p_c = 0{,}0351 МПа.
  2. Числитель с концевой ОФП (оценка сверху): 0,0853⋅50⋅0,85⋅0,0351=0,12720{,}0853\cdot50\cdot0{,}85\cdot0{,}0351 = 0{,}1272.
  3. В пласте: ℓc=0,1272/(3,0⋅0,0708)=0,60\ell_c = 0{,}1272/(3{,}0\cdot0{,}0708) = 0{,}60 м — 0,12 % расстояния 500 м между рядами скважин.
  4. В керне: ℓc=0,1272/(3,0⋅3,0)=1,4\ell_c = 0{,}1272/(3{,}0\cdot3{,}0) = 1{,}4 см — 20 % длины образца.
  5. Если бы керн прокачивали с пластовой скоростью 0,0708 м/сут, ℓc\ell_c был бы 0,60 м — в 8,6 раза длиннее образца.
Ответ: В пласте фронт размыт на 0,60 м из 500 м; в опыте на керне при 3 м/сут — на 1,4 см из 7 см.

Что это значит. В масштабе пласта капиллярные силы размывают фронт на десятки сантиметров (в более проницаемой породе или при более медленной закачке — на метры), и в расчёте вытеснения его спокойно считают скачком. В керне те же силы сопоставимы с длиной образца: при пластовой скорости опыт показал бы в основном капиллярную пропитку, а у выходного торца вода задерживалась бы капиллярными силами — это концевой эффект (как измеряют ОФП). Поэтому нестационарные опыты ведут при повышенных расходах или на длинных составных образцах, а результат проверяют на независимость от расхода. Оценка с концевой ОФП завышает ℓc\ell_c: в зоне фронта krok_{ro} меньше. Данные примера условные.

Кривая доли воды как инструмент

Построив кривую fw(Sw)f_w(S_w) по кривым ОФП и вязкостям, инженер получает ответы на несколько практических вопросов сразу.

Обводнённость скважины
доля воды на забое равна fwf_w при насыщенности у скважины; на поверхности её пересчитывают через объёмные коэффициенты (формула обводнённости). Обратно: по замеренной обводнённости находят fwf_w, а по кривой — насыщенность у скважины (пример ниже).
Дебит нефти при заданном отборе жидкости
в пластовых условиях нефть составляет долю 1−fw1 - f_w потока: qo=(1−fw) qt/Boq_o = (1 - f_w)\,q_t/B_o в стандартных кубометрах.
Скорость насыщенности
наклон кривой fw′= dfw/ dSwf_w' = \dd f_w/\dd S_w пропорционален скорости, с которой данная насыщенность продвигается по пласту; касательная к кривой из точки (Swi,0)(S_{wi}, 0) даёт насыщенность на фронте вытеснения, момент прорыва и среднюю насыщенность за фронтом (Теория Баклея–Леверетта и метод Уэлджа).
Сравнение вариантов
всё, что сдвигает кривую вправо (менее вязкая нефть или более вязкая вода, закачка вверх по восстанию, гидрофильная порода), улучшает вытеснение; сдвиг влево — ухудшает. Так на одном графике сравнивают полимерное заводнение, горячую воду, выбор стороны закачки.
Расчётный примерОбводнённость 90 %: сколько нефти осталось у скважины
Исходные данные (условные)
Добывающая скважина
обводнённость продукции 90 % (стандартные условия)
Флюиды
нефть 3,0 мПа·с, Bo=1,25B_o = 1{,}25; вода 0,5 мПа·с, Bw=1,02B_w = 1{,}02 (как в уроке Относительная проницаемость)
Кривые ОФП
урока Относительная проницаемость; пласт горизонтальный
Для сравнения
та же обводнённость при нефти 30 мПа·с
Решение
  1. Доля воды на забое (пластовые условия): fw=0,90⋅1,020,90⋅1,02+0,10⋅1,25=0,9181,043=0,880f_w = \dfrac{0{,}90\cdot1{,}02}{0{,}90\cdot1{,}02 + 0{,}10\cdot1{,}25} = \dfrac{0{,}918}{1{,}043} = 0{,}880.
  2. По кривой для нефти 3,0 мПа·с fw=0,880f_w = 0{,}880 при Sw≈S_w \approx 0,569. Нефтенасыщенность у скважины 1−0,569=0,4311 - 0{,}569 = 0{,}431, из них подвижной нефти 0,431−0,25=0,1810{,}431 - 0{,}25 = 0{,}181 — 34 % подвижного интервала 0,53.
  3. Для нефти 30 мПа·с та же доля воды — при Sw≈S_w \approx 0,460: подвижной нефти у скважины ещё 0,75−0,460=0,2900{,}75 - 0{,}460 = 0{,}290, то есть 55 % интервала.
  4. Дебит нефти при отборе жидкости 255 м³/сут в пластовых условиях: (1−0,880)⋅255/1,25=24,5(1 - 0{,}880)\cdot255/1{,}25 = 24{,}5 м³/сут.
Ответ: Обводнённости 90 % соответствует доля воды на забое 0,880 и водонасыщенность у скважины около 0,569 (при вязкой нефти — около 0,460).

Что это значит. Обводнённость сама по себе не говорит, насколько промыт пласт у скважины. При нефти 3,0 мПа·с 90 % воды означает, что вытеснено около двух третей подвижной нефти; при нефти 30 мПа·с — меньше половины, и большая её часть ещё впереди, хотя добывать её придётся с очень большим количеством воды. Кривая fwf_w переводит одно в другое. Данные примера условные.

Термины и обозначения

Доля воды в потоке fwf_w
доля воды в суммарном пластовом расходе нефти и воды через сечение пласта.
Вязкая часть λw/λt\lambda_w/\lambda_t
доля воды в горизонтальном пласте без капиллярных сил; зависит только от насыщенности, ОФП и вязкостей.
Соотношение подвижностей M′M'
(krw′/μw)/(kro′/μo)(k_{rw}'/\mu_w)/(k_{ro}'/\mu_o) — отношение подвижности воды за фронтом к подвижности нефти перед ним; быстрая характеристика формы кривой.
Гравитационное число NgN_g
отношение силы тяжести вдоль пласта к вязкому сопротивлению нефти при данной скорости фильтрации; положительно при движении воды вверх по восстанию.
Капиллярная длина ℓc\ell_c
масштаб, на котором капиллярные силы размывают фронт вытеснения; обратно пропорциональна скорости фильтрации.
Диффузный поток
течение, при котором насыщенность одинакова по толщине пласта; условие применимости одномерной кривой fwf_w.
Противоток
встречное движение фаз (вода вниз, нефть вверх) под действием силы тяжести или капиллярных сил.

Частые ошибки и заблуждения

Неверно

«Доля воды в потоке равна водонасыщенности: если в порах 50 % воды, то и в потоке примерно половина».

На самом деле

Расход делится по подвижностям. В примере урока при Sw=0,45S_w = 0{,}45 вода даёт от 19 до 86 % потока в зависимости от вязкости нефти, а в гидрофобной породе при Sw=0,40S_w = 0{,}40 — 85 %.

Неверно

«Если закачивать быстрее, вода прорвётся раньше и кривая fwf_w станет хуже».

На самом деле

В горизонтальном пласте без капиллярных сил кривая от темпа не зависит: быстрее закачка — быстрее проходят те же состояния. Темп влияет через силу тяжести (в наклонном пласте) и капиллярные силы, а в неоднородном пласте — через перетоки между прослоями.

Неверно

«Капиллярное давление мало по сравнению с перепадом между скважинами, значит, на вытеснение оно не влияет нигде».

На самом деле

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

Неверно

«Сила тяжести в одномерном расчёте оказалась мала — значит, в этом пласте её можно не учитывать».

На самом деле

Одномерный расчёт учитывает только силу тяжести вдоль пласта. В мощном пласте вода стекает к подошве поперёк напластования, и это влияет на вытеснение сильнее; такой случай рассчитывают иначе (Сегрегированный поток и гравитационное вытеснение).

Подробно: для продвинутых

Вывод уравнения доли воды

Закон Дарси для фаз при одномерном потоке вдоль оси xx, наклонённой под углом α\alpha к горизонту (двухфазный закон Дарси), выразим относительно градиентов давлений:

Градиенты давления фаз из закона Дарси
−∂po∂x=μo uok kro+ρogsin⁡α,−∂pw∂x=μw uwk krw+ρwgsin⁡α-\frac{\partial p_o}{\partial x} = \frac{\mu_o\,u_o}{k\,k_{ro}} + \rho_o g\sin\alpha,\qquad -\frac{\partial p_w}{\partial x} = \frac{\mu_w\,u_w}{k\,k_{rw}} + \rho_w g\sin\alpha
где:
  • uo, uwu_o,\ u_wскорости фильтрации нефти и воды, м/с
  • po, pwp_o,\ p_wдавления в нефти и в воде, Па
  • ρo, ρw\rho_o,\ \rho_wплотности нефти и воды, кг/м³

Единицы: СИ

Справедливо, если
  • те же, что для уравнения доли воды.

Вычтем из второго равенства первое. Слева получится ∂(po−pw)/∂x=∂pc/∂x\partial(p_o - p_w)/\partial x = \partial p_c/\partial x, справа — разность вязких членов и Δρ gsin⁡α\Delta\rho\,g\sin\alpha:

Разность уравнений Дарси для нефти и воды
∂pc∂x=μw uwk krw−μo uok kro+Δρ gsin⁡α\frac{\partial p_c}{\partial x} = \frac{\mu_w\,u_w}{k\,k_{rw}} - \frac{\mu_o\,u_o}{k\,k_{ro}} + \Delta\rho\,g\sin\alpha
где:
  • pcp_cкапиллярное давление po−pwp_o - p_w, Па
  • Δρ\Delta\rhoρw−ρo\rho_w - \rho_o, кг/м³

Единицы: СИ

Справедливо, если
  • те же, что для уравнения доли воды.

Подставим uo=ut−uwu_o = u_t - u_w и соберём члены с uwu_w слева:

Скорость воды через суммарную скорость
uw(μwk krw+μok kro)=μo utk kro+∂pc∂x−Δρ gsin⁡αu_w\left(\frac{\mu_w}{k\,k_{rw}} + \frac{\mu_o}{k\,k_{ro}}\right) = \frac{\mu_o\,u_t}{k\,k_{ro}} + \frac{\partial p_c}{\partial x} - \Delta\rho\,g\sin\alpha
где:
  • uw, utu_w,\ u_tскорость фильтрации воды и суммарная скорость фильтрации, м/с

Единицы: СИ

Справедливо, если
  • те же, что для уравнения доли воды.

Умножим обе части на k kro/(μout)k\,k_{ro}/(\mu_o u_t) и учтём, что fw=uw/utf_w = u_w/u_t. Слева останется fwf_w, умноженная на знаменатель общей формулы, справа — её числитель; деление на знаменатель и даёт общую форму. Если же сумму в скобке слева записать через подвижности, 1/λw+1/λo=λt/(λoλw)1/\lambda_w + 1/\lambda_o = \lambda_t/(\lambda_o\lambda_w), и разделить на неё, получится форма через подвижности.

Формула для газа, вытесняющего нефть, выводится так же: вместо воды — газ, вместо po−pwp_o - p_w — разность давлений нефти и газа, а разность плотностей ρg−ρo\rho_g - \rho_o отрицательна. Газ, закачанный в повышенную часть залежи и вытесняющий нефть вниз по падению, сила тяжести стабилизирует; вверх по восстанию — гонит вперёд (Вытеснение нефти газом).

Точка перегиба и наклон кривой

Для кривых Кори в горизонтальном пласте наклон доли воды записывается явно (он же используется в уроке Схемы IMPES и полностью неявная):

Наклон кривой доли воды для модели Кори
fw′= dfw dSw=fw (1−fw)1−Swi−Sor(nwSn+no1−Sn)f_w' = \frac{\dd f_w}{\dd S_w} = \frac{f_w\,(1 - f_w)}{1 - S_{wi} - S_{or}}\left(\frac{n_w}{S_n} + \frac{n_o}{1 - S_n}\right)
где:
  • fw′f_w'наклон кривой доли воды, —
  • SnS_nнормированная водонасыщенность (Sw−Swi)/(1−Swi−Sor)(S_w - S_{wi})/(1 - S_{wi} - S_{or}), —
  • nw, non_w,\ n_oпоказатели Кори для воды и нефти, —

Единицы: безразмерная

Справедливо, если
  • модель Кори, горизонтальный пласт, капиллярный член не учитывается;
  • на концах подвижного интервала, где течёт одна фаза, наклон равен нулю.

С ростом вязкости нефти точка перегиба смещается к меньшим насыщенностям; наибольший наклон для нефтей от 1 до 3 мПа·с почти не меняется, а для более вязких растёт: для кривых урока при нефти 1,2 мПа·с перегиб около Sw=0,52S_w = 0{,}52 с наклоном 5,2, при 3,0 мПа·с — около 0,47 (5,2), при 30 мПа·с — около 0,36 (6,7). Насыщенности вблизи перегиба продвигаются по пласту быстрее всех, и именно из-за этого на фронте образуется скачок (скачок насыщенности). Наибольший наклон определяет и предельный шаг явной схемы в модели (Схемы IMPES и полностью неявная).

Связь с другими темами

Проверьте себя

  1. 1

    Горизонтальный однородный пласт, капиллярными силами пренебрегают. Темп закачки увеличили вдвое. Что произойдёт с кривой fw(Sw)f_w(S_w)?

    1. А.Сдвинется вправо: вода будет вытеснять нефть лучше
    2. Б.Сдвинется влево: вода прорвётся раньше
    3. В.Не изменится: в этих условиях кривая от темпа не зависит
    4. Г.Станет круче в середине
    Показать ответ

    Правильный ответ: В. Не изменится: в этих условиях кривая от темпа не зависит

    Без силы тяжести и капиллярного давления fw=λw/λtf_w = \lambda_w/\lambda_t зависит только от насыщенности, ОФП и вязкостей: скорость фильтрации сокращается. Быстрая закачка проходит те же состояния за меньшее время.

  2. 2

    При некоторой водонасыщенности kro=0,12k_{ro} = 0{,}12, krw=0,05k_{rw} = 0{,}05, вязкость нефти 4,0 мПа·с, воды 0,6 мПа·с. Найдите долю воды в потоке горизонтального пласта.

    Показать ответ

    Подвижности на единицу проницаемости: λo=0,12/4,0=0,030\lambda_o = 0{,}12/4{,}0 = 0{,}030, λw=0,05/0,6≈0,0833\lambda_w = 0{,}05/0{,}6 \approx 0{,}0833 1/(мПа·с). fw=0,0833/(0,0833+0,030)≈0,74f_w = 0{,}0833/(0{,}0833 + 0{,}030) \approx 0{,}74: ОФП воды почти в два с половиной раза меньше ОФП нефти, но вода в 6,7 раза менее вязкая, и почти три четверти потока — вода.

  3. 3

    Воду можно закачивать в пониженную часть наклонного пласта (она пойдёт вверх по восстанию) или в повышенную (пойдёт вниз по падению). Где доля воды в потоке при той же насыщенности меньше и почему?

    1. А.Одинаково: сила тяжести действует на обе фазы
    2. Б.При закачке в пониженную часть: вес воды тормозит её движение вверх
    3. В.При закачке в повышенную часть: вода стекает вниз и не участвует в потоке
    Показать ответ

    Правильный ответ: Б. При закачке в пониженную часть: вес воды тормозит её движение вверх

    Вода тяжелее нефти, и при движении вверх по восстанию сила тяжести придерживает её, а нефть «всплывает»: член −Δρ gsin⁡α-\Delta\rho\,g\sin\alpha при α>0\alpha > 0 уменьшает fwf_w. Сила тяжести действует на обе фазы, но важна разность плотностей.

  4. 4

    Гравитационное число пласта при расходе 800 м³/сут равно 0,2. Каким оно станет при 200 м³/сут и во сколько раз уменьшится доля воды по сравнению с горизонтальным пластом при насыщенности, где kro/kro′=0,5k_{ro}/k_{ro}' = 0{,}5 (вода идёт вверх по восстанию)?

    Показать ответ

    NgN_g обратно пропорционально расходу: 0,2⋅800/200=0,80{,}2\cdot800/200 = 0{,}8. Множитель 1−Ng kro/kro′=1−0,8⋅0,5=0,61 - N_g\,k_{ro}/k_{ro}' = 1 - 0{,}8\cdot0{,}5 = 0{,}6: доля воды составит 60 % от доли в горизонтальном пласте, то есть уменьшится примерно в 1,7 раза. При 800 м³/сут множитель был бы 0,9.

  5. 5

    В какой породе при одинаковых флюидах и одинаковой водонасыщенности доля воды в потоке, скорее всего, выше?

    1. А.В гидрофильной
    2. Б.В гидрофобной
    3. В.Смачиваемость на долю воды не влияет
    Показать ответ

    Правильный ответ: Б. В гидрофобной

    В гидрофобной породе вода — несмачивающая фаза и занимает крупные поры, поэтому её ОФП высока, а ОФП нефти, оставшейся в мелких порах и плёнках, мала. Кривая fwf_w сдвинута влево: в примере урока при Sw=0,40S_w = 0{,}40 — 85 % воды против 16 % в гидрофильной породе.

  6. 6

    Почему капиллярным членом в уравнении доли воды обычно пренебрегают при расчёте заводнения пласта, но не при обработке опыта на керне длиной 5 см, прокачанного с пластовой скоростью?

    Показать ответ

    Капиллярные силы размывают фронт на длину порядка ℓc=k kroΔpc/(μout)\ell_c = k\,k_{ro}\Delta p_c/(\mu_o u_t). При пластовых скоростях это от десятков сантиметров до метров — мало по сравнению с сотнями метров между скважинами, но больше длины образца. В керне капиллярная пропитка и концевой эффект определили бы результат, поэтому опыты ведут при повышенных расходах и проверяют, что результат от расхода не зависит.

Обновлено: 14 октября 2026 г.