Skip to content
Brent$117,6Газ TTF$1 192USD453,87 ₸
MunaiHub
ГДИСРаспределённые замеры давления Продвинутый 20 мин чтения

Перезаряд пласта и оценка проницаемости по опробователю

Перезаряд пласта (суперзарядка, supercharging) — превышение давления у стенки скважины над пластовым из-за фильтрата бурового раствора, который продолжает уходить через глинистую корку в породу. В плотном пропластке фильтрат не успевает растечься, и опробователь регистрирует пластовое давление плюс это превышение. Проницаемость по претесту — оценки проницаемости по записи предварительного теста: по установившейся депрессии, по восстановлению давления при сферическом притоке, по отклику наблюдательного зонда и по испытанию интервала двухпакерным модулем. Каждая из них относится к своему объёму породы и к своему сочетанию горизонтальной и вертикальной проницаемостей.

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

В плотных пропластках, плотных водоносных частях и карбонатах перезаряд сдвигает давление на единицы, десятки и сотни килопаскалей. Этого хватает, чтобы ошибиться в плотности флюида, сместить зеркало свободной воды (ЗСВ) на метры и принять перезаряженный пропласток за отдельную систему давления. Распознать перезаряд помогает подвижность, а её, как и проницаемость, приходится оценивать по тем же претестам — до спуска колонны и без притока на поверхность. Понимание того, что измеряет каждый способ, нужно и для контроля замеров давления, и для использования проницаемости по претестам в модели пласта.

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

  • Скважина с репрессией работает как медленная нагнетательная: фильтрат идёт через корку со скоростью, которую задаёт в основном корка. Пласт должен отвести этот поток, и давление у стенки поднимается тем выше, чем больше его сопротивление.
  • Превышение пропорционально скорости фильтрации и вязкости фильтрата и обратно пропорционально проницаемости, с поправкой на логарифм времени: в пропластке в 10 раз менее проницаемом оно почти в 10 раз больше. Сверху его ограничивает репрессия.
  • Перезаряженная точка выглядит кондиционной: давление стабилизировалось. Выдают её низкая подвижность, сдвиг вправо от линии, растущий с обратной подвижностью, и повышенный разброс точек.
  • Подвижность по депрессии отражает сантиметры породы у зонда; в слоистой породе она лежит между сферической и горизонтальной проницаемостями. Восстановление при сферическом притоке даёт сферическую проницаемость (kh2kv)1/3(k_h^2k_v)^{1/3} на десятках сантиметров — метрах.
  • Ёмкость прибора искажает восстановление ещё долго после того, как запись выглядит стабильной: наклон по слишком ранним точкам занижает проницаемость в разы.
  • Наблюдательный зонд над зондом-стоком даёт по уровню отклика khk_h, по его наклону во времени — ksk_s, а вместе — kvk_v; двухпакерный модуль даёт радиальный приток к интервалу.
Ключевые соотношения
Δps=Δpреп RплRк+Rпл\Delta p_s = \Delta p_{\text{реп}}\,\dfrac{R_{\text{пл}}}{R_{\text{к}} + R_{\text{пл}}}

перезаряд как делитель давления: корка и пласт соединены последовательно

Δps≈u2rwμf2kf[u1u2ln⁡Td+tsts+ln⁡kftsϕμfctrw2+0,809]\Delta p_s \approx \dfrac{u_2 r_w \mu_f}{2k_f}\left[\dfrac{u_1}{u_2}\ln\dfrac{T_d + t_s}{t_s} + \ln\dfrac{k_f t_s}{\phi\mu_f c_t r_w^2} + 0{,}809\right]

модель нагнетания фильтрата, СИ

p(Δt)=p∗−ms(1Δt−1tp+Δt)p(\Delta t) = p^* - m_s\left(\dfrac{1}{\sqrt{\Delta t}} - \dfrac{1}{\sqrt{t_p + \Delta t}}\right)

сферическое восстановление после претеста длительностью tpt_p

ms=qμϕμct4π3/2ks3/2,ks=(kh2kv)1/3m_s = \dfrac{q\mu\sqrt{\phi\mu c_t}}{4\pi^{3/2}k_s^{3/2}},\quad k_s = (k_h^2k_v)^{1/3}

наклон сферической прямой и сферическая проницаемость, СИ

Где чаще всего ошибаются
  • Отбраковывать только нестабилизированные записи: перезаряженная запись стабилизирована.
  • Считать, что экстраполированное давление p∗p^* снимает перезаряд: оно равно давлению у стенки вместе с превышением.
  • Брать наклон сферического восстановления по первым минутам: там ещё работает ёмкость прибора, и проницаемость занижена в разы.
  • Сравнивать подвижность по депрессии с проницаемостью по керну или по исследованию скважины без поправок на фильтрат, анизотропию и масштаб.

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

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

Почему давление у стенки выше пластового

Пока ствол открыт и давление раствора выше пластового, скважина работает как очень медленная нагнетательная. Фильтрат уходит в проницаемые породы в три этапа (динамика фильтрации): мгновенная начальная фильтрация (spurt loss), пока не сложилась корка; динамическая фильтрация при циркуляции, когда поток раствора смывает внешний слой корки и её толщина, а с ней и скорость фильтрации, выходят на постоянное значение; статическая фильтрация после остановки циркуляции — скорость сразу заметно падает, а потом продолжает медленно уменьшаться по мере роста корки. Каротаж и опробование идут именно в статический период.

Перепад между давлением раствора и пластовым давлением делят два последовательных сопротивления: корка и порода за ней. Это делитель давления, как в электрической цепи. Обычно сопротивление корки на порядки больше, почти весь перепад приходится на неё, и давление у стенки почти не отличается от пластового. Но сопротивление породы растёт с уменьшением проницаемости, и в плотном пропластке его доля становится заметной.

Перезаряд как делитель давления
Δps=ps−pпл=Δpреп RплRк+Rпл,Rк=μf hmckmc\Delta p_s = p_s - p_{\text{пл}} = \Delta p_{\text{реп}}\,\frac{R_{\text{пл}}}{R_{\text{к}} + R_{\text{пл}}},\qquad R_{\text{к}} = \frac{\mu_f\,h_{mc}}{k_{mc}}
где:
  • Δps\Delta p_sперезаряд — превышение давления у стенки psp_s над пластовым pплp_{\text{пл}}, Па
  • Δpреп\Delta p_{\text{реп}}репрессия: давление раствора минус пластовое давление, Па
  • RкR_{\text{к}}сопротивление корки на единицу площади стенки, Па·с/м
  • RплR_{\text{пл}}сопротивление породы на единицу площади стенки: превышение давления, делённое на скорость фильтрации uu (следующий раздел), Па·с/м
  • μf, hmc, kmc\mu_f,\ h_{mc},\ k_{mc}вязкость фильтрата, толщина и проницаемость корки, Па·с, м, м²

Единицы: СИ

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

Пока Rпл≪RкR_{\text{пл}} \ll R_{\text{к}}, скорость фильтрации задаёт корка, u≈Δpреп/Rкu \approx \Delta p_{\text{реп}}/R_{\text{к}}, и перезаряд равен этой скорости, умноженной на сопротивление породы, — он обратно пропорционален проницаемости. Когда сопротивления сравнимы, поток уменьшается, и перезаряд растёт медленнее; больше репрессии он быть не может.

Превышение давления над пластовым в зависимости от расстояния от стенки0,010,11100,1110100Расстояние от стенки, мПерезаряд, кПа0,3 мД30 мД

График в логарифмических осях: по горизонтали расстояние от стенки скважины от 0,01 до 50 метров, по вертикали превышение давления от 0,1 до 200 кПа. Верхняя кривая для 0,3 мД: около 82 кПа у стенки, 54 кПа в 1 метре, 27 кПа в 10 метрах, дальше она убывает быстрее и к 50 метрам опускается до 3 кПа. Нижняя кривая для 30 мД: 1,1 кПа у стенки, 0,8 кПа в 1 метре, 0,55 кПа в 10 метрах и 0,36 кПа в 50 метрах. У стенки кривые различаются примерно в 75 раз.

РасчётРасчётное превышение давления над пластовым в зависимости от расстояния от стенки для пропластков 30 и 0,3 мД при одной и той же скорости фильтрации (модель следующего раздела, параметры примера). Почти вся репрессия, около 1,9 МПа, теряется на корке и на графике не показана.

Повышенное давление охватывает не тонкую оболочку у стенки, а область в десятки метров: с расстоянием оно убывает лишь логарифмически, как давление вокруг нагнетательной скважины. Поэтому и восстановление после претеста, которое охватывает десятки сантиметров — метры породы, целиком лежит внутри этой области и видит то же превышение.

Сколько: модель нагнетания фильтрата

Сопротивление породы находят, рассматривая скважину как нагнетательную. Через каждый квадратный метр стенки уходит поток uu; на метр толщины пласта это 2πrwu2\pi r_w u. Сначала, при циркуляции, скорость u1u_1 держится в течение TdT_d, потом, в статический период, она меньше — u2u_2. Каждый этап — радиальная закачка с постоянным расходом, и решение для линейного источника (логарифмическое приближение) вместе с суперпозицией двух расходов дают давление на стенке.

Перезаряд по модели нагнетания фильтрата
Δps=u2 rw μf2kf[u1u2ln⁡Td+tsts+ln⁡kf tsϕ μf ct rw2+0,809]\Delta p_s = \frac{u_2\,r_w\,\mu_f}{2k_f}\left[\frac{u_1}{u_2}\ln\frac{T_d + t_s}{t_s} + \ln\frac{k_f\,t_s}{\phi\,\mu_f\,c_t\,r_w^2} + 0{,}809\right]
где:
  • u1, u2u_1,\ u_2скорость фильтрации через единицу площади стенки при циркуляции и после её остановки (объём фильтрата за единицу времени на единицу площади), м/с
  • TdT_dпродолжительность динамической фильтрации перед остановкой циркуляции, с
  • tst_sвремя статической фильтрации до замера, с
  • rwr_wрадиус ствола, м
  • kfk_fпроницаемость породы для фильтрата: в водоносном пропластке — проницаемость, в нефтенасыщенном — эффективная проницаемость промытой зоны по фильтрату, м²
  • μf, ϕ, ct\mu_f,\ \phi,\ c_tвязкость фильтрата, пористость и общая сжимаемость породы у ствола, Па·с, —, 1/Па
  • 0,8090{,}809ln⁡(4/γ)\ln(4/\gamma), где γ=1,781\gamma = 1{,}781, — как в логарифмическом приближении решения для линейного источника, —

Единицы: СИ; логарифм натуральный

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

Множитель перед скобкой — то, что задаёт масштаб: скорость фильтрации, вязкость фильтрата и сопротивление породы 1/kf1/k_f. Скобка меняется медленно: первое слагаемое — память о динамической фильтрации — убывает со временем покоя, второе растёт как логарифм времени. Поэтому перезаряд почти обратно пропорционален проницаемости и мало меняется в течение каротажа.

В других системах единиц (1)
Метрические единицы
Δps [кПа]=140,7 u2 rw μfkf[u1u2ln⁡Td+tsts+ln⁡kf tsϕ μf ct rw2−4,831]\Delta p_s\,[\text{кПа}] = 140{,}7\,\frac{u_2\,r_w\,\mu_f}{k_f}\left[\frac{u_1}{u_2}\ln\frac{T_d + t_s}{t_s} + \ln\frac{k_f\,t_s}{\phi\,\mu_f\,c_t\,r_w^2} - 4{,}831\right]

uu — мм/ч, rwr_w — м, μf\mu_f — мПа·с, kfk_f — мД, TdT_d и tst_s — ч, ctc_t — 1/МПа.

140,7=(10−3/3600)⋅10−3/(2⋅0,9869⋅10−15)/1000140{,}7 = (10^{-3}/3600)\cdot 10^{-3}/(2\cdot 0{,}9869\cdot 10^{-15})/1000; −4,831=ln⁡(0,9869⋅10−15⋅3600/10−9)+0,809-4{,}831 = \ln(0{,}9869\cdot10^{-15}\cdot 3600/10^{-9}) + 0{,}809. Скорость 0,3 мм/ч — это 0,05 см³ фильтрата в минуту через 100 см² стенки.

Расчётный примерПерезаряд водоносных пропластков 30, 3 и 0,3 мД
Исходные данные (условные)
Фильтрация
0,72 мм/ч при циркуляции в течение 30 ч, затем 0,30 мм/ч; замер через 10 ч после остановки циркуляции
Ствол и фильтрат
rwr_w = 0,108 м (долото 215,9 мм); фильтрат на водной основе, 0,40 мПа·с
Порода у ствола
ϕ\phi = 0,18; ct=1,0⋅10−3c_t = 1{,}0\cdot10^{-3} МПа⁻¹
Репрессия
1,9 МПа
Решение
  1. Память о циркуляции: 0,720,30ln⁡30+1010=2,4⋅1,3863=3,327\dfrac{0{,}72}{0{,}30}\ln\dfrac{30 + 10}{10} = 2{,}4\cdot 1{,}3863 = 3{,}327.
  2. Знаменатель под логарифмом: 0,18⋅0,40⋅1,0⋅10−3⋅0,1082=8,40⋅10−70{,}18\cdot 0{,}40\cdot 1{,}0\cdot10^{-3}\cdot 0{,}108^2 = 8{,}40\cdot10^{-7}. Для 30 мД: ln⁡(30⋅10/(8,40⋅10−7))−4,831=14,86\ln(30\cdot 10/(8{,}40\cdot10^{-7})) - 4{,}831 = 14{,}86; скобка 18,19.
  3. Множитель: 140,7⋅0,30⋅0,108⋅0,40/30=0,0608140{,}7\cdot 0{,}30\cdot 0{,}108\cdot 0{,}40/30 = 0{,}0608 кПа; перезаряд 0,0608⋅18,19=1,110{,}0608\cdot 18{,}19 = 1{,}11 кПа.
  4. Для 3 мД: множитель 0,608 кПа, скобка 3,327+12,56=15,893{,}327 + 12{,}56 = 15{,}89, перезаряд 9,7 кПа. Для 0,3 мД: 6,08 кПа, скобка 13,58, перезаряд 82,6 кПа.
  5. Проверка допущения: даже для 0,3 мД перезаряд — 4,3 % репрессии; по делителю давления он 82,6/(1+82,6/1900)=79,282{,}6/(1 + 82{,}6/1900) = 79{,}2 кПа — поправка мала.
Ответ: Перезаряд — около 1,11 кПа в пропластке 30 мД, 9,7 кПа в 3 мД и 82,6 кПа в 0,3 мД.

Что это значит. При разбросе кондиционных точек около 1 кПа перезаряд пропластка 30 мД почти незаметен, 3-миллидарсиевого — уже в десять раз больше разброса, а 0,3-миллидарсиевого — десятки килопаскалей, то есть метры ошибки контакта. Ориентиры «перезаряд важен ниже такой-то проницаемости» зависят от раствора: при впятеро меньшей скорости фильтрации все числа впятеро меньше. Скорость статической фильтрации в скважине почти никогда не измеряют, поэтому модель служит для понимания и порядков величин, а сам перезаряд распознают по данным съёмки (ниже). Данные примера условные.

Перезаряд в зависимости от проницаемости для двух растворов0,11101000,11101001000Проницаемость, мДПерезаряд, кПарепрессия 1,9 МПаразброс ≈ 1 кПа0,30 мм/ч0,06 мм/ч

График в логарифмических осях: по горизонтали проницаемость от 0,03 до 100 мД, по вертикали перезаряд от 0,1 до 3000 кПа. Две почти прямые нисходящие линии: для раствора с фильтрацией 0,30 мм/ч перезаряд около 1 кПа при 30 мД, 10 кПа при 3 мД и 80 кПа при 0,3 мД; для 0,06 мм/ч — в пять раз меньше. При самой низкой проницаемости сплошная линия для 0,30 мм/ч загибается: при 0,03 мД перезаряд около 500 кПа против 690 кПа по штриховой линии постоянной скорости фильтрации. Горизонтальная штриховая линия на 1900 кПа — репрессия, предел перезаряда; пунктир на 1 кПа отмечает обычный разброс кондиционных точек.

РасчётРасчётный перезаряд в зависимости от проницаемости для двух растворов: статическая фильтрация 0,30 и 0,06 мм/ч (динамическая в 2,4 раза больше), остальное — как в примере. Штрих — модель постоянной скорости фильтрации, сплошные линии — с учётом делителя давления при репрессии 1,9 МПа.

Как перезаряд искажает давления, градиенты и контакты

Каждая перезаряженная точка сдвинута вправо на свой Δps\Delta p_s. Если перезаряд одинаков у всех точек линии, линия сдвигается параллельно: градиент и плотность флюида правдоподобны, а уровень неверен. Если перезаряд меняется по разрезу — например, порода внизу плотнее, — меняется и наклон. Сдвиг уровня переходит в сдвиг контакта так же, как ошибка давления в теме Градиенты давления, контакты флюидов и капиллярное равновесие.

Искажение градиента и контакта перезарядом
Gкаж=G+Δ(Δps)Δz,δz∗=Δps,o−Δps,wGw−GoG_{\text{каж}} = G + \frac{\Delta(\Delta p_s)}{\Delta z},\qquad \delta z^* = \frac{\Delta p_{s,o} - \Delta p_{s,w}}{G_w - G_o}
где:
  • Gкаж, GG_{\text{каж}},\ Gкажущийся и истинный градиенты по точкам одной зоны, Па/м
  • Δ(Δps)/Δz\Delta(\Delta p_s)/\Delta zизменение перезаряда по глубине в интервале точек (средний наклон), Па/м
  • δz∗\delta z^*сдвиг найденного ЗСВ (плюс — вниз), м
  • Δps,o, Δps,w\Delta p_{s,o},\ \Delta p_{s,w}перезаряд линий нефти и воды на глубине пересечения, Па

Единицы: СИ (на практике кПа, кПа/м и м)

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

Перезаряд водоносных точек поднимает ЗСВ, перезаряд точек нефтяной части опускает его. Водоносная часть ниже контакта нередко плотнее нефтяной — цементация в ней могла продолжаться и после того, как нефть заполнила ловушку, — поэтому опаснее первый случай: ЗСВ получается слишком высоким, а запасы — заниженными.

Претесты водоносной части скважины 2
Глубина, м
2547,5
Давление, МПа
28,0495
Подвижность, мД/(мПа·с)
4
Глубина, м
2553,0
Давление, МПа
28,0952
Подвижность, мД/(мПа·с)
12
Глубина, м
2559,5
Давление, МПа
28,1700
Подвижность, мД/(мПа·с)
6
Глубина, м
2566,0
Давление, МПа
28,2289
Подвижность, мД/(мПа·с)
25
Глубина, м
2573,5
Давление, МПа
28,3143
Подвижность, мД/(мПа·с)
8

Скважина 2 из темы Градиенты давления, контакты флюидов и капиллярное равновесие: в нефтяной части шесть кондиционных претестов, линия нефти GoG_o = 7,068 кПа/м. Скважину бурили на растворе из примера выше — с большей водоотдачей, чем новую скважину, — а водоносная часть здесь плотнее. Все записи стабилизированы. Данные условные.

Расчётный примерПлотный водоносный пласт: правдоподобная линия и неверный контакт
Исходные данные (условные)
Претесты водоносной части
таблица выше
Линия нефти скважины 2
GoG_o = 7,068 кПа/м; ЗСВ по региональной линии воды — 2540,3 м ± 5 м
Региональная линия воды
градиент 10,50 кПа/м (вода по анализу, 1,0706 г/см³ в пласте); 28,163 МПа на 2560,0 м
Решение
  1. Линия по пяти точкам: Gw=10,22G_w = 10{,}22 кПа/м, плотность 1,042 г/см³ — правдоподобная пластовая вода; разброс s=5,8s = 5{,}8 кПа, в несколько раз больше обычного.
  2. Пересечение этой линии с линией нефти скважины 2: 2535,5 м — на 4,8 м выше, чем по региональной линии.
  3. Отклонения точек от региональной линии воды: 17,8, 5,7, 12,3, 2,9, 9,5 кПа при обратной подвижности 0,250, 0,083, 0,167, 0,040, 0,125 мПа·с/мД: чем ниже подвижность, тем больше отклонение.
  4. Прямая d=c+a/(k/μ)dd = c + a/(k/\mu)_d методом наименьших квадратов: a=72±2a = 72 \pm 2 кПа·мД/(мПа·с), c=0,1±0,4c = 0{,}1 \pm 0{,}4 кПа; остатки 0,4 кПа — отклонения целиком объясняются перезарядом.
  5. Уровень воды без перезаряда — региональная линия плюс cc; её пересечение с линией нефти: ЗСВ = 2540,3 ± 0,5 м. Ошибка уже не содержит сравнения двух скважин: уровень воды взят по замерам той же скважины.
Ответ: Перезаряд поднял найденное ЗСВ на 4,8 м; с поправкой ЗСВ — 2540,3 ± 0,5 м, в согласии с новой скважиной.

Что это значит. Каждая точка в отдельности выглядела кондиционной, а линия — правдоподобной; выдали перезаряд только низкие подвижности, повышенный разброс и связь отклонений с обратной подвижностью. Поправка работает здесь, потому что подвижности разные и скорость фильтрации у пропластков, по-видимому, одинакова. Погрешности подвижности по депрессии ослабляют найденную связь: наклон aa выходит меньше, cc — больше, и поправка остаётся неполной — ЗСВ немного выше истинного; признак этого — остатки заметно больше разброса замеров. Если все точки одинаково плотные, отклонения от подвижности не зависят, перезаряд неотличим от сдвига уровня, и единственная защита — отбраковка по подвижности и региональная линия воды. Приём с известным градиентом воды — частный случай метода вынужденного градиента. Данные примера условные.

Как распознать и учесть перезаряд

Признаки перезаряда
  1. 1
    Подвижность

    перезаряд заметен прежде всего у точек с подвижностью меньше нескольких мД/(мПа·с); при растворе с большой водоотдачей — и при больших подвижностях.

  2. 2
    Положение точки

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

  3. 3
    Разброс и плотность

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

  4. 4
    Повторные претесты

    повторный тест на той же точке даёт другое стабилизированное давление; повторная съёмка через несколько часов — меньшее превышение.

  5. 5
    Время стабилизации

    в плотном пропластке давление дольше подходит к стабилизации: и время ёмкости прибора, и время сферического «хвоста» обратно пропорциональны подвижности (ёмкость прибора). Долгое восстановление — признак риска перезаряда.

Как учесть. Первое средство — отбраковка: в градиент и контакты не берут точки с низкой подвижностью, если их давление не подтверждено иначе. Если зону можно оценить только по таким точкам, их используют как границы: стабилизированное давление перезаряженной точки — верхняя граница пластового давления, а незаконченное восстановление — нижняя. Приближённую поправку дают отклонения, связанные с обратной подвижностью, как в примере выше; она годится как оценка порядка величины, потому что скорость фильтрации у разных пропластков может различаться в разы. Так, в съёмке темы об опробователе точка 2526,5 м перезаряжена на 46 кПа, а у точки 2520,0 м с вшестеро большей подвижностью превышения не видно, хотя при одинаковой скорости фильтрации оно было бы около 9 кПа. Наконец, перезаряд уменьшают при планировании: меньшая репрессия, раствор с низкой водоотдачей, повторный замер на подозрительных глубинах, линия воды по региональным данным для плотного водоносного пласта (региональная линия воды).

Проницаемость по претесту: депрессия и восстановление

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

Восстановление после претеста несёт сведения о более далёкой породе. Возмущение от отбора расходится от зонда сферически, и вдали от него давление описывает точечный сток в бесконечной среде: установившаяся часть плюс «хвост», убывающий как 1/t1/\sqrt t. Восстановление — суперпозиция отбора с расходом qq в течение tpt_p и его прекращения (остановка как отрицательный расход): установившиеся части сокращаются, и остаются два «хвоста». Так получается функция сферического времени fsf_s (spherical time function), и в координатах «давление — fsf_s» запись ложится на прямую.

Восстановление давления после претеста при сферическом притоке
p(Δt)=p∗−ms fs(Δt),fs=1Δt−1tp+Δt,ms=q μϕ μ ct4π3/2 ks3/2p(\Delta t) = p^* - m_s\,f_s(\Delta t),\qquad f_s = \frac{1}{\sqrt{\Delta t}} - \frac{1}{\sqrt{t_p + \Delta t}},\qquad m_s = \frac{q\,\mu\sqrt{\phi\,\mu\,c_t}}{4\pi^{3/2}\,k_s^{3/2}}
где:
  • p(Δt)p(\Delta t)давление у зонда через время Δt\Delta t после остановки поршня, Па
  • p∗p^*давление, к которому стремится восстановление (fs→0f_s \to 0): давление у стенки, Па
  • tpt_pпродолжительность отбора, с
  • qqскорость отбора поршнем, м³/с
  • msm_sнаклон сферической прямой, Па·с^½
  • ksk_sсферическая проницаемость (kh2kv)1/3(k_h^2 k_v)^{1/3}, м²
  • μ, ϕ, ct\mu,\ \phi,\ c_tвязкость притекающего флюида (обычно фильтрата), пористость и общая сжимаемость породы у ствола, Па·с, —, 1/Па

Единицы: СИ

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

Наклон зависит от ks3/2k_s^{3/2} и лишь как от корня — от пористости и сжимаемости, поэтому их неточность влияет на ksk_s слабо: в степени 1/3. По отрезку на оси fs=0f_s = 0 находят p∗p^*, даже если восстановление не закончилось.

В других системах единиц (1)
Метрические единицы
ms [МПа⋅с1/2]=45,8 q [см3/с]  μ [мПа⋅с]ϕ μ ct [МПа−1]ks3/2 [мД]m_s\,[\text{МПа·с}^{1/2}] = 45{,}8\,\frac{q\,[\text{см}^3/\text{с}]\;\mu\,[\text{мПа·с}]\sqrt{\phi\,\mu\,c_t\,[\text{МПа}^{-1}]}}{k_s^{3/2}\,[\text{мД}]}

qq — см³/с, μ\mu — мПа·с, ctc_t — 1/МПа, ksk_s — мД, Δt\Delta t и tpt_p — с, давление — МПа.

45,8=10−6⋅10−310−3⋅10−64π3/2(0,9869⋅10−15)3/2⋅10−645{,}8 = \dfrac{10^{-6}\cdot10^{-3}\sqrt{10^{-3}\cdot10^{-6}}}{4\pi^{3/2}(0{,}9869\cdot10^{-15})^{3/2}}\cdot 10^{-6}; отсюда ks=(45,8 qμϕμct/ms)2/3k_s = (45{,}8\,q\mu\sqrt{\phi\mu c_t}/m_s)^{2/3}.

Расчётный примерСферическая проницаемость по восстановлению на 2520,0 м
Исходные данные (условные)
Претест
повторный на точке 2520,0 м из темы Замеры давления опробователем пластов на кабеле: 5 см³ со скоростью 0,2 см³/с (tpt_p = 25 с); прибор держали на точке 6 минут
Конец отбора
24,1254 МПа; пластовое давление на этой точке 27,8141 МПа
Восстановление
через 40, 80, 160 и 320 с: 27,81139; 27,81342; 27,81386; 27,81401 МПа (разрешение 0,01 кПа)
Порода и фильтрат
ϕ\phi = 0,18; ct=1,0⋅10−3c_t = 1{,}0\cdot10^{-3} МПа⁻¹; μf\mu_f = 0,40 мПа·с; ёмкость прибора 0,25 см³/МПа
Решение
  1. Подвижность по депрессии: 39,02⋅0,2/(27,8141−24,1254)=2,1239{,}02\cdot 0{,}2/(27{,}8141 - 24{,}1254) = 2{,}12 мД/(мПа·с); постоянная времени ёмкости 0,25⋅39,02/2,1=4,650{,}25\cdot 39{,}02/2{,}1 = 4{,}65 с.
  2. Функция сферического времени: 1/40−1/65=0,034081/\sqrt{40} - 1/\sqrt{65} = 0{,}03408; для 80, 160 и 320 с — 0,01421, 0,00554 и 0,00206 с^−½.
  3. Постоянная: 45,8⋅0,2⋅0,400,18⋅0,40⋅1,0⋅10−3=0,0310945{,}8\cdot 0{,}2\cdot 0{,}40\sqrt{0{,}18\cdot 0{,}40\cdot 1{,}0\cdot10^{-3}} = 0{,}03109.
  4. Ранние точки (40 и 80 с): ms=(27,81342−27,81139)/(0,03408−0,01421)=0,1022m_s = (27{,}81342 - 27{,}81139)/(0{,}03408 - 0{,}01421) = 0{,}1022 МПа·с^½, ks=(0,03109/0,1022)2/3=0,45k_s = (0{,}03109/0{,}1022)^{2/3} = 0{,}45 мД — подвижность 1,1.
  5. Поздние точки (160 и 320 с): ms=0,00015/(0,00554−0,00206)=0,0431m_s = 0{,}00015/(0{,}00554 - 0{,}00206) = 0{,}0431 МПа·с^½, ks=0,80k_s = 0{,}80 мД — подвижность 2,0; p∗=27,81401+0,0431⋅0,00206=27,8141p^* = 27{,}81401 + 0{,}0431\cdot 0{,}00206 = 27{,}8141 МПа.
Ответ: По поздним точкам сферическая подвижность — около 2,0 мД/(мПа·с), по депрессии — 2,12; по ранним точкам вышло бы 1,1, почти вдвое меньше.

Что это значит. Порода в этой модели однородна и изотропна, у зонда нет повреждения, поэтому депрессия и восстановление должны дать одно и то же, — и дают, если брать восстановление достаточно поздно. Ранние точки лежат через 9–17 постоянных времени ёмкости после остановки поршня, запись там уже кажется стабилизированной, но приток, который дожимает жидкость в линии, ещё искажает наклон. Поздние же изменения давления — всего 0,15 кПа, полтора десятка шагов разрешения манометра: окно пригодных точек узкое и с ранней, и с поздней стороны. Экстраполированное давление совпало с пластовым, потому что перезаряда на этой точке нет. При депрессии давление в зонде опускается примерно до 24,1 МПа — на 3,3 МПа ниже давления насыщения нефти этой залежи (27,4 МПа); модель однофазная и не учитывает газ, который может выделиться из остаточной нефти у зонда: он увеличил бы сжимаемость и исказил начало восстановления. Данные примера условные.

Восстановление после претеста в координатах давление — функция сферического времени0,000,010,020,030,04012345fs, с−1/2p − 27,810 МПа, кПабез ёмкостипо точкам 40 и 80 сзапись

График: по горизонтали функция сферического времени от 0 до 0,04 с в степени минус одна вторая, по вертикали давление минус 27,810 МПа от 0 до 5 кПа. Кривая восстановления поднимается справа налево: около 0,1 кПа при значении функции 0,039 и 4,1 кПа у оси, где она приходит к пластовому давлению. Серая штриховая сферическая прямая идёт от 4,1 кПа на оси вниз вправо; кривая сначала лежит заметно ниже неё, а при малых значениях функции сливается с ней. Оранжевая штриховая прямая через две правые точки (40 и 80 секунд) круче и пересекает ось около 4,9 кПа.

РасчётВосстановление после претеста на 2520,0 м в координатах «давление — функция сферического времени» (модель с ёмкостью прибора 0,25 см³/МПа); по вертикали — превышение давления над 27,810 МПа. Точки — замеры из примера; серая штриховая прямая — сферическая прямая той же породы без ёмкости, оранжевая — прямая по двум ранним точкам.

Сферический и радиальный приток, анизотропия

В слоистой породе вертикальная проницаемость меньше горизонтальной, и сферический приток превращается в эллипсоидальный. Если растянуть вертикальную координату в ks/kv\sqrt{k_s/k_v} раз, а горизонтальные сжать в kh/ks\sqrt{k_h/k_s} раз, порода становится изотропной с проницаемостью ks=(kh2kv)1/3k_s = (k_h^2k_v)^{1/3}, а объёмы и расходы не меняются. Поэтому вдали от зонда, где форма самого зонда не важна, восстановление «видит» именно сферическую проницаемость. Депрессия же определяется течением в первых сантиметрах у круглого отверстия зонда на стенке, и там анизотропия сказывается иначе: подвижность по депрессии лежит между ksk_s и khk_h и с ростом анизотропии отходит от khk_h медленнее, чем ksk_s (дополнительный раздел).

Депрессионная и сферическая проницаемость в анизотропной породе (доли khk_h)
A=kh/kvA = k_h/k_v
1
kd/khk_d/k_h
1,00
ks/khk_s/k_h
1,00
kd/ksk_d/k_s
1,00
A=kh/kvA = k_h/k_v
3
kd/khk_d/k_h
0,77
ks/khk_s/k_h
0,69
kd/ksk_d/k_s
1,12
A=kh/kvA = k_h/k_v
10
kd/khk_d/k_h
0,61
ks/khk_s/k_h
0,46
kd/ksk_d/k_s
1,31
A=kh/kvA = k_h/k_v
30
kd/khk_d/k_h
0,51
ks/khk_s/k_h
0,32
kd/ksk_d/k_s
1,57
A=kh/kvA = k_h/k_v
100
kd/khk_d/k_h
0,43
ks/khk_s/k_h
0,22
kd/ksk_d/k_s
1,97

ks/kh=A−1/3k_s/k_h = A^{-1/3}; kd/khk_d/k_h — модель круглого зонда на плоской непроницаемой стенке (дополнительный раздел); решения с эллипсоидальным стоком вместо отверстия дают значения kd/khk_d/k_h примерно на 10 % меньше.

Отсюда обычный порядок оценок в слоистой однородной породе: kv<ks<kd<khk_v < k_s < k_d < k_h. Отношение kd/ksk_d/k_s растёт с анизотропией медленно, поэтому по паре «депрессия — восстановление» анизотропию оценить можно лишь грубо: ошибка в 10 % в одной из оценок при AA около 10 меняет найденную анизотропию примерно вдвое. Надёжнее измерять отклик породы на расстоянии.

Радиальный приток. Если зонд стоит в тонком пропластке между непроницаемыми слоями, возмущение доходит до них, и сферический приток сменяется радиальным (цилиндрическим) по всей толщине пропластка: на диагностическом графике производная выходит на полку (Диагностика по производной давления), а по полулогарифмической прямой находят khhk_h h пропластка. Двухпакерный модуль (dual packer) изолирует участок стенки длиной около метра и сразу создаёт радиальный приток к этому интервалу; дальше следуют сферический и, если пропласток ограничен, поздний радиальный приток — те же режимы, что у частично вскрытой скважины, только в масштабе метра (режимы течения). Большой объём между пакерами даёт и большую ёмкость, поэтому такие испытания ведут насосом в течение десятков минут.

Наблюдательный зонд. В многозондовом приборе над зондом-стоком на расстоянии Δz\Delta z (обычно около 0,7 м) по той же образующей стоит наблюдательный зонд (observation probe), иногда ещё один — на противоположной стороне ствола. Насос отбирает флюид через сток с постоянным расходом, а наблюдательный зонд записывает, как до него доходит возмущение. Это вертикальное гидропрослушивание (vertical interference test): на отклик не влияют ни повреждение у стока, ни его ёмкость.

Отклик наблюдательного зонда при постоянном отборе
Δpн(t)=qμ4πkh Δz erfc⁡Δz4ηvt  ≈  qμ4πkh Δz−mst,ηv=kvϕμct\Delta p_{\text{н}}(t) = \frac{q\mu}{4\pi k_h\,\Delta z}\,\operatorname{erfc}\frac{\Delta z}{\sqrt{4\eta_v t}} \;\approx\; \frac{q\mu}{4\pi k_h\,\Delta z} - \frac{m_s}{\sqrt t},\qquad \eta_v = \frac{k_v}{\phi\mu c_t}
где:
  • Δpн\Delta p_{\text{н}}снижение давления у наблюдательного зонда, Па
  • Δz\Delta zрасстояние между зондами по вертикали, м
  • ηv\eta_vпьезопроводность по вертикали, м²/с
  • msm_sтот же наклон, что у сферического восстановления: qμϕμct/(4π3/2ks3/2)q\mu\sqrt{\phi\mu c_t}/(4\pi^{3/2}k_s^{3/2}), Па·с^½
  • q, μ, ϕ, ct, kh, kvq,\ \mu,\ \phi,\ c_t,\ k_h,\ k_vрасход насоса, вязкость, пористость, сжимаемость, горизонтальная и вертикальная проницаемости, м³/с, Па·с, —, 1/Па, м², м²

Единицы: СИ; приближённая форма — при Δz/4ηvt<0,15\Delta z/\sqrt{4\eta_v t} < 0{,}15

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

Уровень, к которому стремится отклик, зависит только от горизонтальной проницаемости: вертикальная растягивает путь, но во столько же раз облегчает его. Скорость, с которой отклик к нему подходит, задаёт вертикальная пьезопроводность. На поздних временах отклик — прямая в координатах «Δpн\Delta p_{\text{н}} — 1/t1/\sqrt t»: отрезок на оси даёт khk_h, наклон — ksk_s, а вместе они дают kv=ks3/kh2k_v = k_s^3/k_h^2.

В других системах единиц (1)
Метрические единицы
Δpн,∞ [кПа]=80,6 q [см3/с]  μ [мПа⋅с]kh [мД]  Δz [м]\Delta p_{\text{н},\infty}\,[\text{кПа}] = 80{,}6\,\frac{q\,[\text{см}^3/\text{с}]\;\mu\,[\text{мПа·с}]}{k_h\,[\text{мД}]\;\Delta z\,[\text{м}]}

qq — см³/с, μ\mu — мПа·с, khk_h — мД, Δz\Delta z — м; наклон msm_s — по метрической форме сферического восстановления (в МПа·с^½).

80,6=10−6⋅10−3/(4π⋅0,9869⋅10−15)/100080{,}6 = 10^{-6}\cdot10^{-3}/(4\pi\cdot 0{,}9869\cdot10^{-15})/1000.

Отклик наблюдательного зонда при постоянном отборе10100100001234Время после пуска насоса, сСнижение давления, кПапредел 3,84 кПа1 мД по вертикали0,25 мД по вертикали

График: по горизонтали время от 10 до 5000 секунд в логарифмическом масштабе, по вертикали снижение давления у наблюдательного зонда от 0 до 4,2 кПа. Горизонтальная штриховая линия на 3,84 кПа — предельный уровень. Кривая для вертикальной проницаемости 1 мД: 0,6 кПа через 10 секунд, 2,5 кПа через 100 секунд, 3,2 кПа через 450 секунд и 3,5 кПа через 1800 секунд; три точки отмечают замеры примера. Кривая для 0,25 мД сдвинута вправо: почти ноль через 10 секунд, 1,4 кПа через 100 секунд, 3,0 кПа через 1000 секунд и 3,45 кПа к 5000 секундам.

РасчётРасчётный отклик наблюдательного зонда в 0,70 м над стоком при отборе 1 см³/с для двух значений вертикальной проницаемости при khk_h = 12 мД. Уровень, к которому стремятся кривые, одинаков; вертикальная проницаемость определяет, как быстро к нему подходит отклик.
Расчётный примерГоризонтальная и вертикальная проницаемость по наблюдательному зонду
Исходные данные (условные)
Испытание
водоносный пропласток; насос отбирает 1,0 см³/с через зонд-сток; наблюдательный зонд на 0,70 м выше
Флюид и порода
вода, 0,40 мПа·с; ϕ\phi = 0,20; ct=1,0⋅10−3c_t = 1{,}0\cdot10^{-3} МПа⁻¹
Отклик наблюдательного зонда
через 450, 900 и 1800 с после пуска: 3,20, 3,39 и 3,52 кПа
Решение
  1. Абсциссы 1/t1/\sqrt t: 0,04714, 0,03333 и 0,02357 с^−½; проверка условия: Δz/4ηvt\Delta z/\sqrt{4\eta_v t} при kvk_v около 1 мД и 450 с — около 0,15.
  2. Прямая по крайним точкам: наклон (3,52−3,20)/(0,04714−0,02357)=13,58(3{,}52 - 3{,}20)/(0{,}04714 - 0{,}02357) = 13{,}58 кПа·с^½; уровень 3,52+13,58⋅0,02357=3,8403{,}52 + 13{,}58\cdot 0{,}02357 = 3{,}840 кПа; средняя точка на прямой — 3,39 кПа, замер 3,39.
  3. Горизонтальная проницаемость: kh=80,6⋅1,0⋅0,40/(0,70⋅3,840)=12,0k_h = 80{,}6\cdot 1{,}0\cdot 0{,}40/(0{,}70\cdot 3{,}840) = 12{,}0 мД.
  4. Сферическая: ks=(45,8⋅1,0⋅0,400,20⋅0,40⋅10−3/0,01358)2/3=5,26k_s = \left(45{,}8\cdot 1{,}0\cdot 0{,}40\sqrt{0{,}20\cdot 0{,}40\cdot 10^{-3}}/0{,}01358\right)^{2/3} = 5{,}26 мД.
  5. Вертикальная: kv=ks3/kh2=1,01k_v = k_s^3/k_h^2 = 1{,}01 мД; анизотропия kh/kv=11,9k_h/k_v = 11{,}9.
Ответ: kh≈12,0k_h \approx 12{,}0 мД, kv≈1,01k_v \approx 1{,}01 мД; анизотропия около 12.

Что это значит. Зонд-сток в этом пропластке дал бы по депрессии около 7,1 мД, а по восстановлению — около 5,26: обе оценки лежат между kvk_v и khk_h, и только наблюдательный зонд разделяет их. Числа относятся к объёму в метр-полтора вокруг зондов. Если между зондами есть тонкий плотный прослой, отклик и «вертикальная проницаемость» отражают именно его, и это ценно для прогноза конусообразования и вертикального вытеснения (Конусообразование воды и газа). Данные примера условные.

Ёмкость прибора и пределы метода

Между зондом и манометром — линия и камера с жидкостью, и её сжатие и расширение — ёмкость прибора (tool storage) — отделяет то, что делает поршень, от того, что происходит в породе (модель ёмкости). На депрессии это постоянная времени τ=SK/(k/μ)\tau = SK/(k/\mu): если поршень работает меньше нескольких τ\tau, депрессия не успевает установиться, и подвижность завышена. На восстановлении действие ёмкости тоньше и длится дольше. После остановки поршня порода ещё дожимает жидкость в линии, и часть объёма она отдаёт с запозданием. Сферическая формула же считает, что весь объём ушёл из породы за время tpt_p. Поэтому прямая на сферическом графике появляется лишь через десятки τ\tau, а наклон по более ранним точкам слишком крут, и проницаемость получается заниженной.

Подвижность по наклону сферического графика в разных окнах (модель с ёмкостью прибора)
Точка и ёмкость
2520,0 м; 0,25 см³/МПа
Окно после остановки
40–160 с
Подвижность по наклону
1,26
Истинная
2,10
Точка и ёмкость
2520,0 м; 0,25 см³/МПа
Окно после остановки
160–640 с
Подвижность по наклону
1,97
Истинная
2,10
Точка и ёмкость
2520,0 м; 0,05 см³/МПа
Окно после остановки
40–160 с
Подвижность по наклону
2,01
Истинная
2,10
Точка и ёмкость
2526,5 м; 0,25 см³/МПа
Окно после остановки
60–240 с
Подвижность по наклону
0,013
Истинная
0,35
Точка и ёмкость
2526,5 м; 0,25 см³/МПа
Окно после остановки
960–3840 с
Подвижность по наклону
0,327
Истинная
0,35

Подвижности — мД/(мПа·с). Точка 2520,0 м — претест примера; 2526,5 м — плотный перезаряженный пропласток съёмки темы об опробователе: отбор 2 см³ со скоростью 0,05 см³/с. Модель: сферический сток, дающий ту же депрессию, что стандартный зонд, с ёмкостью линии.

Отсюда и одна из причин, почему на практике проницаемость по восстановлению часто оказывается заметно ниже, чем по депрессии: в плотных пропластках прибор снимают раньше, чем ёмкость перестаёт искажать наклон. Ёмкость уменьшают, сокращая объём линии между зондом и манометром; анализ ведут по поздним точкам или по модели, в которой ёмкость учтена явно. Верхний предел метода задаёт разрешение манометра: сигнал восстановления убывает как k−3/2k^{-3/2}, и при претесте 5 см³ со скоростью 0,5 см³/с изменение давления от Δt=tp\Delta t = t_p до 10 tp10\,t_p превышает двадцать шагов кварцевого манометра (0,01 кПа) лишь до 11 мД, а при шаге 0,7 кПа — до 0,6 мД. Депрессия работает в гораздо более проницаемых породах.

Что ограничивает оценки проницаемости по претесту
Фактор
Объём исследования
Депрессия
сантиметры у зонда
Восстановление
десятки сантиметров — метры; растёт как η Δt\sqrt{\eta\,\Delta t}
Фактор
Повреждение, уплотнение или раскол породы зондом
Депрессия
сильно: занижает, раскол — завышает
Восстановление
почти не влияет
Фактор
Анизотропия
Депрессия
между ksk_s и khk_h
Восстановление
ks=(kh2kv)1/3k_s = (k_h^2k_v)^{1/3}
Фактор
Зона проникновения
Депрессия
фильтрат при остаточной нефти
Восстановление
может выйти за зону проникновения: насыщенности смешанные, вязкость неясна
Фактор
Ёмкость прибора
Депрессия
завышает, если отбор короче нескольких τ\tau
Восстановление
занижает по ранним точкам, иногда в разы
Фактор
Разрешение манометра
Депрессия
в проницаемых породах депрессия мала, но измерима; мешают скорее помехи самого прибора
Восстановление
предел — единицы — десятки мД
Фактор
Давление насыщения
Депрессия
глубокая депрессия может выделить газ в линии
Восстановление
искажённое восстановление

Направления влияния — для обычных случаев. Сопоставление оценок по претестам с керном и исследованиями скважин разобрано в теме Сопоставление опробователя, керна и гидродинамических исследований.

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

Перезаряд Δps\Delta p_s
превышение давления у стенки над пластовым из-за фильтрата, который пласт не успевает отвести; растёт с уменьшением проницаемости.
Динамическая и статическая фильтрация
уход фильтрата при циркуляции (корка в равновесии, скорость постоянна) и после её остановки (скорость меньше и медленно убывает).
Функция сферического времени fsf_s
1/Δt−1/tp+Δt1/\sqrt{\Delta t} - 1/\sqrt{t_p + \Delta t}; в координатах «давление — fsf_s» сферическое восстановление — прямая.
Сферическая проницаемость ksk_s
(kh2kv)1/3(k_h^2k_v)^{1/3} — проницаемость, которую определяет восстановление при сферическом притоке.
p∗p^*
давление, к которому стремится восстановление: давление у стенки, вместе с перезарядом.
Вертикальное гидропрослушивание
испытание многозондовым прибором: отбор через один зонд, запись отклика другим на некотором расстоянии по вертикали; даёт khk_h и kvk_v.

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

Неверно

«Давление стабилизировалось — значит, оно пластовое».

На самом деле

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

Неверно

«Перезаряд бывает только при подвижности меньше 1 мД/(мПа·с)».

На самом деле

Перезаряд пропорционален скорости фильтрации через корку: при растворе с большой водоотдачей он заметен и при единицах мД/(мПа·с). Граница зависит от раствора, репрессии и того, какая точность нужна для контактов.

Неверно

«Восстановление охватывает больший объём, значит, его проницаемость всегда точнее депрессионной».

На самом деле

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

Неверно

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

На самом деле

В однородной породе уровень отклика задаёт горизонтальная проницаемость; вертикальная определяет, как быстро он устанавливается. Маленький уровень — высокая khk_h или малый расход; медленный рост — низкая kvk_v. Отклика нет совсем — вероятен сплошной непроницаемый прослой между зондами.

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

Откуда сферическая прямая

Для точечного стока с постоянным расходом qq в бесконечной изотропной среде снижение давления на расстоянии rr равно qμ4πkrerfc⁡r2ηt\dfrac{q\mu}{4\pi k r}\operatorname{erfc}\dfrac{r}{2\sqrt{\eta t}}, где η=k/(ϕμct)\eta = k/(\phi\mu c_t). При малом аргументе erfc⁡x≈1−2x/π\operatorname{erfc}x \approx 1 - 2x/\sqrt\pi, и у стока, где rr мал, снижение — это установившаяся часть qμ/(4πkr)q\mu/(4\pi k r) минус «хвост» qμ4πk1πηt\dfrac{q\mu}{4\pi k}\dfrac{1}{\sqrt{\pi\eta t}}, который от rr уже не зависит.

При восстановлении к отбору с t=0t = 0 добавляют закачку того же расхода с момента tpt_p. Установившиеся части сокращаются, а «хвосты» дают p∗−p=qμ4π3/2kη(1Δt−1tp+Δt)p^* - p = \dfrac{q\mu}{4\pi^{3/2} k\sqrt\eta}\left(\dfrac{1}{\sqrt{\Delta t}} - \dfrac{1}{\sqrt{t_p + \Delta t}}\right). Подставив η=k/(ϕμct)\sqrt\eta = \sqrt{k/(\phi\mu c_t)}, получаем наклон ms=qμϕμct/(4π3/2k3/2)m_s = q\mu\sqrt{\phi\mu c_t}/(4\pi^{3/2}k^{3/2}). В анизотропной среде после растяжения координат kk заменяют на ksk_s.

Почему полная сфера, а не полусфера у стенки? Сток стоит на стенке ствола, и при малых rr стенка действительно служит плоским экраном. Но восстановление анализируют, когда возмущение ушло на десятки сантиметров, и ствол радиусом около 0,1 м для него — тонкая непроницаемая трубка, почти не загораживающая пространство. Поэтому поздний «хвост» — сферический; принять полусферу значило бы удвоить наклон в формуле и завысить ksk_s в 22/3≈1,62^{2/3} \approx 1{,}6 раза. Для одиночного отбора с tp≪Δtt_p \ll \Delta t функция fs≈tp/(2Δt3/2)f_s \approx t_p/(2\Delta t^{3/2}): восстановление тем быстрее, чем выше ksk_s, а время, за которое давление подходит к p∗p^* на заданную малую величину, обратно пропорционально подвижности.

Депрессия у круглого зонда в анизотропной породе

Пусть отверстие зонда — круг радиуса rpr_p на плоской непроницаемой стенке, давление на нём одинаково. В изотропной породе такой диск пропускает установившийся поток q=4rp (k/μ) Δpq = 4 r_p\,(k/\mu)\,\Delta p — классический результат теории потенциала (вдвое меньше, чем у такого же диска в сплошной среде, потому что приток идёт только с одной стороны).

В анизотропной породе растянем координаты так, чтобы среда стала изотропной с проницаемостью ksk_s. Стенка остаётся плоскостью, а круг превращается в эллипс с полуосями rpA1/3r_pA^{1/3} по вертикали и rpA−1/6r_pA^{-1/6} по горизонтали, где A=kh/kvA = k_h/k_v. Проводимость эллиптического диска пропорциональна его большой полуоси и обратно пропорциональна полному эллиптическому интегралу K(m)K(m) с параметром m=1−(b/a)2=1−1/Am = 1 - (b/a)^2 = 1 - 1/A. После обратной подстановки получается kd/kh=π/(2K(m))k_d/k_h = \pi/(2K(m)): при A=1A = 1 — единица, при A=10A = 10 — 0,61, при A=100A = 100 — 0,43. Депрессия теряет горизонтальную проницаемость лишь логарифмически медленно, тогда как сферическая падает как A−1/3A^{-1/3}.

Модель идеализирована: стенка плоская, давление по всему отверстию одинаково. Решения, в которых отверстие заменено эллипсоидальным стоком, дают при той же анизотропии значения kd/khk_d/k_h примерно на 10 % меньше. Направление и порядок эффекта от этого не меняются: ks<kd<khk_s < k_d < k_h.

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

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

  1. 1

    Два водоносных пропластка в одной скважине: 40 и 2 мД, раствор и время замера одинаковы. Во сколько примерно раз различаются их перезаряды?

    1. А.Примерно в 4,5 раза (20\sqrt{20})
    2. Б.Примерно в 15–20 раз, у плотного больше
    3. В.Одинаковы
    4. Г.Примерно в 400 раз
    Показать ответ

    Правильный ответ: Б. Примерно в 15–20 раз, у плотного больше

    Множитель перед скобкой обратно пропорционален проницаемости — он различается в 20 раз, а скобка меняется лишь логарифмически: у плотного пропластка она немного меньше. С параметрами примера урока скобки около 15,5 и 18,5, и отношение 20⋅15,5/18,5≈1720\cdot 15{,}5/18{,}5 \approx 17.

  2. 2

    Найдите перезаряд пропластка 1,0 мД по модели урока: множитель 140,7⋅0,30⋅0,108⋅0,40/1,0140{,}7\cdot 0{,}30\cdot 0{,}108\cdot 0{,}40/1{,}0 кПа и скобка 14,8.

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

    Множитель 140,7⋅0,30⋅0,108⋅0,40=1,82140{,}7\cdot 0{,}30\cdot 0{,}108\cdot 0{,}40 = 1{,}82 кПа; перезаряд 1,82⋅14,8≈271{,}82\cdot 14{,}8 \approx 27 кПа. При разности градиентов воды и нефти 3,4 кПа/м такой перезаряд водоносной точки сдвинул бы ЗСВ вверх почти на 8 м.

  3. 3

    Стабилизированные давления пяти водоносных претестов отклоняются от линии воды по анализу на 3, 12, 6, 20 и 9 кПа; их подвижности — 30, 8, 15, 5 и 10 мД/(мПа·с). Что это скорее всего?

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

    Правильный ответ: Б. Перезаряд плотных пропластков

    Отклонение растёт с обратной подвижностью почти пропорционально: 3⋅30≈12⋅8≈6⋅15≈20⋅5≈9⋅10≈903\cdot 30 \approx 12\cdot 8 \approx 6\cdot 15 \approx 20\cdot 5 \approx 9\cdot 10 \approx 90–100100 кПа·мД/(мПа·с). Ошибка плотности дала бы отклонения, растущие с глубиной, а не с подвижностью; негерметичность — давление раствора.

  4. 4

    На сферическом графике прямая по точкам через 30–60 с после остановки поршня дала подвижность 0,6 мД/(мПа·с), а по депрессии — 2,4. Постоянная времени ёмкости — около 4 с. Какой вывод правильнее?

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

    Окно 30–60 с — всего 8–15 постоянных времени ёмкости: приток, дожимающий жидкость в линии, ещё искажает наклон, и проницаемость по нему занижена. Сначала проверяют поздние точки (через несколько минут) или анализируют восстановление по модели с ёмкостью; только если расхождение сохранится, ищут анизотропию, неоднородность или повреждение у зонда — причём повреждение занизило бы как раз депрессионную оценку.

  5. 5

    Наблюдательный зонд в 0,7 м над стоком при отборе 2 см³/с показал предельное снижение 5,0 кПа; вязкость 0,5 мПа·с. Найдите горизонтальную проницаемость в идеализированной модели и скажите, в какую сторону её смещает неучтённый ствол.

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

    kh=80,6⋅2⋅0,5/(0,7⋅5,0)=23k_h = 80{,}6\cdot 2\cdot 0{,}5/(0{,}7\cdot 5{,}0) = 23 мД. Ствол — непроницаемое препятствие; он увеличивает отклик по сравнению с бесконечной средой, поэтому идеализированная формула занижает khk_h (не больше чем вдвое — предел плоской стенки).

  6. 6

    Почему экстраполированное давление p∗p^* сферического восстановления не годится для поправки на перезаряд?

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

    Правильный ответ: Б. Оно равно давлению у стенки, а перезаряд и есть превышение давления у стенки над пластовым

    Восстановление после претеста — отклик породы в сантиметрах–метрах от зонда, внутри области, где фильтрат уже поднял давление. Прямая стремится к давлению у стенки, и перезаряд входит в p∗p^* целиком. p∗p^* помогает, когда восстановление не успело закончиться, но не отделяет перезаряд от пластового давления.

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