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

Полуограниченные системы: теория и анализ исследований

Полуограниченная система (semi-infinite system) — пласт, ограниченный непроницаемыми границами не со всех сторон: хотя бы в одном направлении он открыт и простирается далеко за пределы исследования. Простейшие такие системы — полуплоскость за одним разломом, прямой угол между двумя перпендикулярными разломами, канал между параллельными разломами и канал, закрытый с одного конца (U-образная система разломов). Псевдоустановившегося режима у них нет: поздняя производная давления выходит на новую полку или идёт с наклоном ½, а после остановки давление восстанавливается до начального. Урок собирает эти системы в одну картину — общее решение с отражениями, одна характерная длина, правило произведения для прямоугольных систем — и подробно разбирает канал, закрытый с одного конца.

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

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

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

  • Полуограниченная система открыта хотя бы в одном направлении: псевдоустановившегося режима нет, а давление после долгой остановки возвращается к начальному pip_i, а не к среднему давлению блока, как в замкнутой системе.
  • Каждая прямая граница — отражение. Если положение всех отражений задаёт одна длина LL, их вклад зависит только от ηt/L2\eta t/L^2: одной геометрии — одна эталонная кривая.
  • Правило произведения: в прямоугольной системе производная, отнесённая к полке mem_e, равна произведению множителей по двум направлениям. Открытое направление даёт 1, одна граница на расстоянии LL — 1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)}, две параллельные — множитель полосы, растущий как 4πηt/W\sqrt{4\pi\eta t}/W.
  • Поэтому торец канала действует на линейный приток так же, как разлом на радиальный: удваивает его. Наклон прямой Δp\Delta p — t\sqrt t становится 2mL2m_L; расстояние до торца надёжнее всего находить по одной точке производной: L3=ηt ln⁡[1/(r−1)]L_3 = \sqrt{\eta t\,\ln[1/(r - 1)]}.
  • Скважина у самой границы ведёт себя как скважина с удвоенным дебитом в системе без этой границы; скин при этом относится к настоящему дебиту.
  • После короткой работы торец канала оставляет на восстановлении давления след, обратно пропорциональный WL3WL_3, а не L2L^2, как разлом в радиальном пласте: поэтому он заметнее, но приходит поздно — около 2L32/η2L_3^2/\eta после остановки.
Ключевые соотношения
Δp=me[ln⁡tD+0,809+2S+∑jE1 ⁣(dj24ηt)]\Delta p = m_e\Big[\ln t_D + 0{,}809 + 2S + \sum_j E_1\!\Big(\dfrac{d_j^2}{4\eta t}\Big)\Big]

давление в скважине с отражениями после окончания влияния ствола; me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh)

Δp′me=Fx(t) Fy(t),F=1+e−L2/(ηt)\dfrac{\Delta p'}{m_e} = F_x(t)\,F_y(t),\qquad F = 1 + \e^{-L^2/(\eta t)}

правило произведения; множитель одной границы на расстоянии LL

L3=ηt ln⁡1r−1,r=Δp′12mLtL_3 = \sqrt{\eta t\,\ln\dfrac{1}{r - 1}},\qquad r = \dfrac{\Delta p'}{\tfrac12 m_L\sqrt t}

расстояние до торца канала по одной точке производной

L3=ηtx/π=0,0336k txϕμctL_3 = \sqrt{\eta t_x/\pi} = 0{,}0336\sqrt{\dfrac{k\,t_x}{\phi\mu c_t}}

по пересечению прямых Δp\Delta p — t\sqrt t: м, мД, ч, мПа·с, 1/МПа (требует очень длинной записи)

δpmax⁡=qB tp198 ϕcth WL3\delta p_{\max} = \dfrac{qB\,t_p}{198\,\phi c_t h\,WL_3}

наибольший след торца на кривой восстановления давления (КВД) после короткой работы tpt_p: м³/сут, ч, 1/МПа, м, МПа

Где чаще всего ошибаются
  • Считать ширину канала по позднему наклону t\sqrt t, когда торец уже сказывается: наклон удвоен, и ширина получается вдвое заниженной.
  • Проводить вторую прямую Δp\Delta p — t\sqrt t через последние точки и брать расстояние до торца по пересечению: поздняя прямая формируется медленно, и расстояние занижено.
  • Принимать восстановление давления, которое медленнее обычного выходит к pip_i, за истощение блока, не проверив открытое направление.
  • Удваивать вместе с дебитом и скин скважины, стоящей у самой границы.
  • Искать по исследованию одной скважины, с какой стороны закрыт канал: запись даёт расстояния, а не направления.

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

  • объяснить, чем полуограниченная система отличается от замкнутой и от неограниченного пласта, и назвать её признаки на производной и при восстановлении давления;
  • записать давление в скважине с отражениями в безразмерном виде с одной характерной длиной и объяснить, почему поздняя часть решения зависит только от ηt/L2\eta t/L^2;
  • строить производную прямоугольной системы разломов как произведение множителей по двум направлениям;
  • узнать канал, закрытый с одного конца, по удвоению наклона линейного притока и найти ширину канала и расстояние до торца по прямым Δp\Delta p — t\sqrt t и по одной точке производной;
  • оценить, увидит ли исследование торец канала при работе скважины и после остановки, и назвать, что запись различить не может.

Полуограниченные системы среди моделей границ

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

Простейшие полуограниченные системы и замкнутый блок для сравнения
Система
Один разлом (полуплоскость)
Отражения
одно, на 2L2L
Поздняя производная
полка 2me2m_e
Система
Прямой угол
Отражения
три
Поздняя производная
полка 4me4m_e
Система
Канал (две параллельные границы)
Отражения
бесконечный ряд
Поздняя производная
наклон ½
Система
Канал, закрытый с одного конца
Отражения
два ряда
Поздняя производная
наклон ½, вдвое выше
Где разобрано
этот урок
Система
Замкнутый прямоугольник
Отражения
решётка
Поздняя производная
наклон 1

LL — расстояние до разлома, me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh) — уровень полки производной при радиальном притоке (Диагностика по производной давления). Клин с произвольным углом — тоже полуограниченная система; его поздний уровень (360∘/θ) me(360^\circ/\theta)\,m_e разобран в уроке Несколько границ: параллельные разломы и клиновидные блоки.

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

Отражения скважины в канале, закрытом с одного концаL3скважинаторецстенки каналаканал открыт →отражения в стенкахзеркальный рядряд продолжается

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

СхемаКанал, закрытый с одного конца (схема, не в масштабе). Скважина посередине канала в L3L_3 от торца. Отражения в стенках образуют вертикальный ряд, как у открытого канала; торец отражает весь этот ряд — появляется второй, зеркальный ряд на расстоянии 2L32L_3. В сторону, где канал открыт, отражений нет.

Решение с отражениями и одна характерная длина

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

Давление в скважине с отражениями в безразмерном виде
pwD=12(ln⁡tD+0,809)+S+12∑jE1 ⁣(βj24τ),τ=ηtL2=tDLD2,βj=djLp_{wD} = \tfrac12\left(\ln t_D + 0{,}809\right) + S + \tfrac12\sum_j E_1\!\left(\frac{\beta_j^2}{4\tau}\right),\qquad \tau = \frac{\eta t}{L^2} = \frac{t_D}{L_D^2},\quad \beta_j = \frac{d_j}{L}
где:
  • pwDp_{wD}безразмерное снижение забойного давления 2πkh Δp/(qBμ)2\pi kh\,\Delta p/(qB\mu), —
  • tD, St_D,\ Sбезразмерное время kt/(ϕμctrw2)kt/(\phi\mu c_tr_w^2) и скин, —
  • L, LDL,\ L_Dхарактерная длина системы — например, расстояние до ближней границы — и она же в радиусах скважины, м, —
  • τ\tauбезразмерное время, отнесённое к L2L^2, —
  • dj, βjd_j,\ \beta_jрасстояние от скважины до jj-го отражения и оно же в единицах LL, м, —
  • E1E_1интегральная показательная функция (Решение для линейного источника и функция давления), —

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

Справедливо, если
  • прямые непроницаемые границы через всю толщину однородного пласта; отражения — линейные источники (dj≫rwd_j \gg r_w);
  • постоянный дебит, влияние ствола закончилось, tD>25t_D > 25;
  • для отражений, уходящих в бесконечность (каналы), сумма берётся до тех пор, пока слагаемые не станут пренебрежимо малыми.

Радиус скважины и скин входят только в первые слагаемые — они сдвигают кривую по вертикали. Все границы — в сумму, и она зависит от одной переменной τ\tau, если все расстояния пропорциональны одной длине LL. Поэтому для каждой геометрии строят одну эталонную кривую в координатах τ\tau: совмещение по времени даёт LL, по давлению — проводимость, как в обычном совмещении (Эталонные кривые). Производная — pD′=12[1+∑je−βj2/(4τ)]p_D' = \tfrac12\big[1 + \sum_j\e^{-\beta_j^2/(4\tau)}\big] (отражения отражений).

Если отражений конечное число nn, поздняя производная выходит на полку (1+n) me(1 + n)\,m_e, а наклон полулогарифмической прямой увеличивается во столько же раз: вдвое у одного разлома, вчетверо в прямом угле. Если их бесконечно много, как в канале, новой полки нет: число «почувствованных» отражений растёт как t\sqrt t, и производная идёт с наклоном ½. Отношение позднего уровня к раннему или поздний наклон — первый признак, по которому выбирают модель.

Правило произведения для прямоугольных систем

Если все границы параллельны двум взаимно перпендикулярным направлениям — назовём их xx и yy, — отражения стоят в узлах прямоугольной решётки. Сдвиг каждого отражения складывается из сдвига по xx и сдвига по yy, квадрат расстояния — из их квадратов, а экспонента суммы равна произведению экспонент. Сумма по всем отражениям распадается на произведение двух сумм — по каждому направлению отдельно (правило произведения):

Правило произведения: производная прямоугольной системы
Δp′me=Fx(t) Fy(t),Fx=∑iexp⁡ ⁣(−ξi24ηt),Fy=∑lexp⁡ ⁣(−ζl24ηt)\frac{\Delta p'}{m_e} = F_x(t)\,F_y(t),\qquad F_x = \sum_i \exp\!\left(-\frac{\xi_i^2}{4\eta t}\right),\quad F_y = \sum_l \exp\!\left(-\frac{\zeta_l^2}{4\eta t}\right)
где:
  • Δp′\Delta p'логарифмическая производная снижения давления  dΔp/ dln⁡t\dd\Delta p/\dd\ln t, Па
  • mem_eполка производной при радиальном притоке, qBμ/(4πkh)qB\mu/(4\pi kh), Па
  • Fx, FyF_x,\ F_yмножители направлений: суммы по одномерным отражениям вдоль xx и вдоль yy, включая саму скважину (сдвиг 0), —
  • ξi, ζl\xi_i,\ \zeta_lсдвиги одномерных отражений по xx и по yy, м
  • η\etaпьезопроводность k/(ϕμct)k/(\phi\mu c_t), м²/с

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

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

Каждое направление вносит свой множитель независимо от другого. Открытое направление — 1 (одна «скважина» без отражений). Одна граница на расстоянии LL — 1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)}: отражение на 2L2L. Две параллельные границы — множитель полосы: сумма по ряду отражений, которая на поздних временах растёт как 4πηt/W\sqrt{4\pi\eta t}/W (вывод — в дополнительном разделе).

Множители направлений и системы, которые из них собираются
По xx
1
По yy
1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)}
Система
один разлом
Поздняя производная
2me2m_e
По xx
1+e−LA2/(ηt)1 + \e^{-L_A^2/(\eta t)}
По yy
1+e−LB2/(ηt)1 + \e^{-L_B^2/(\eta t)}
Система
прямой угол
Поздняя производная
4me4m_e
По xx
1
По yy
полоса ширины WW
Система
канал
Поздняя производная
me4πηt/Wm_e\sqrt{4\pi\eta t}/W
По xx
1+e−L32/(ηt)1 + \e^{-L_3^2/(\eta t)}
По yy
полоса ширины WW
Система
канал, закрытый с одного конца
Поздняя производная
2me4πηt/W2m_e\sqrt{4\pi\eta t}/W
По xx
полоса длины LL
По yy
полоса ширины WW
Система
замкнутый прямоугольник
Поздняя производная
растёт как tt

Для прямого угла правило произведения даёт 1+e−LA2/(ηt)+e−LB2/(ηt)+e−(LA2+LB2)/(ηt)1 + \e^{-L_A^2/(\eta t)} + \e^{-L_B^2/(\eta t)} + \e^{-(L_A^2 + L_B^2)/(\eta t)} — ровно три отражения угла. Для замкнутого прямоугольника те же множители полос описывают выход на псевдоустановившийся режим (Исследование на ограниченность пласта).

Множитель полосы, множитель торца и их произведение для скважины в канале, закрытом с одного конца110100100010 000125102050100Время после пуска t, чМножительштрих — произведение Fx·Fyбирюзовая — полоса Fyоранжевая — торец Fx

График в двойном логарифмическом масштабе: время от 1 до 20 000 часов, множители от 0,8 до 200. Бирюзовая кривая множителя полосы равна единице до нескольких часов, затем поднимается и идёт прямой с наклоном одна вторая: 3,1 к 48 часам, 10 к 480 часам, 31,5 к 4800 часам. Оранжевая кривая множителя торца равна единице примерно до 200 часов, затем поднимается: 1,13 к 480 часам, 1,82 к 4800 часам, и подходит к двум. Тёмная кривая произведения совпадает с бирюзовой до 200 часов, а затем отходит от неё вверх и к 20 000 часам почти вдвое выше.

РасчётРасчётные множители для скважины посередине канала шириной 300 м, закрытого в 1200 м от неё (η\eta = 1480 м²/ч, данные примера ниже): множитель полосы FyF_y, множитель торца FxF_x и их произведение — производная, отнесённая к полке mem_e. Обе шкалы логарифмические.
Расчётный примерПрогноз производной по правилу произведения
Исходные данные (условные)
Пласт
kk = 100 мД; ϕ\phi = 0,20; μ\mu = 1,0 мПа·с; ct=1,2⋅10−3c_t = 1{,}2\cdot10^{-3} МПа⁻¹
Канал по геологической модели
ширина 300 м, скважина посередине; в 1200 м канал, по-видимому, закрыт поперечным разломом
Вопрос
какой будет производная через 2, 20 и 200 сут работы с постоянным дебитом
Решение
  1. Пьезопроводность: η=3,553⋅10−3⋅100/(0,20⋅1,0⋅1,2⋅10−3)=1480\eta = 3{,}553\cdot10^{-3}\cdot100/(0{,}20\cdot1{,}0\cdot1{,}2\cdot10^{-3}) = 1480 м²/ч. Скважина посередине канала, поэтому при ηt/W2>0,2\eta t/W^2 > 0{,}2 множитель полосы равен 2πηt/W22\sqrt{\pi\eta t/W^2} точнее чем на 0,1 %.
  2. 48 ч: ηt=71 040\eta t = 71\,040 м², ηt/W2=0,7893\eta t/W^2 = 0{,}7893. Полоса в линейном притоке: Fy=2π⋅0,7893=3,149F_y = 2\sqrt{\pi\cdot0{,}7893} = 3{,}149. Торец: L32/(ηt)=20,270L_3^2/(\eta t) = 20{,}270, e−20,270\e^{-20{,}270} ничтожно мало, Fx=1F_x = 1. Производная — 3,149 me3{,}149\,m_e: открытый канал.
  3. 480 ч: ηt/W2=7,8933\eta t/W^2 = 7{,}8933, Fy=9,959F_y = 9{,}959; L32/(ηt)=2,027L_3^2/(\eta t) = 2{,}027, Fx=1+e−2,027=1,1317F_x = 1 + \e^{-2{,}027} = 1{,}1317; производная — 9,959⋅1,1317=11,27 me9{,}959\cdot1{,}1317 = 11{,}27\,m_e (сумма всех отражений даёт 11,27).
  4. 4800 ч: ηt/W2=78,9333\eta t/W^2 = 78{,}9333, Fy=31,495F_y = 31{,}495; L32/(ηt)=0,203L_3^2/(\eta t) = 0{,}203, Fx=1,8163F_x = 1{,}8163; производная — 57,2 me57{,}2\,m_e.
  5. Торец поднимет производную на 10 %, когда e−L32/(ηt)=0,1\e^{-L_3^2/(\eta t)} = 0{,}1, то есть при t=L32/(ηln⁡10)=423t = L_3^2/(\eta\ln10) = 423 ч; почти полное удвоение (90 %) — при ηt≈9,5 L32\eta t \approx 9{,}5\,L_3^2, через 385 сут.
Ответ: Через 2 сут производная — около 3,149 mem_e (открытый канал), через 20 сут — 11,27 mem_e (торец поднял её на 13 %), через 200 сут — 57,2 mem_e.

Что это значит. Множитель полосы и множитель торца считаются независимо, и их произведение совпадает с суммой по всем отражениям до последнего знака. Торец в 1200 м заметен уже через 2–3 недели работы, но полностью удвоенную прямую можно увидеть лишь за год — это задача не для исследования, а для постоянных забойных манометров при эксплуатации. Данные примера условные.

Канал, закрытый с одного конца

Канал, закрытый с одного конца, — это вытянутый блок между параллельными разломами, перегороженный поперечным разломом, или русловое песчаное тело, которое выклинивается или упирается в глинистую пробку брошенного русла (Континентальные обстановки: эоловые и речные). В плане границы образуют букву П — отсюда второе название, U-образная система разломов (U-shaped fault system). Пусть скважина стоит в L1L_1 и L2L_2 от стенок (W=L1+L2W = L_1 + L_2) и в L3L_3 от торца.

По правилу произведения производная такой системы — производная открытого канала, умноженная на множитель торца 1+e−L32/(ηt)1 + \e^{-L_3^2/(\eta t)}. Если торец в несколько раз дальше, чем ширина канала, запись проходит четыре участка: полку радиального притока mem_e; подъём, пока сказываются стенки; линейный приток открытого канала — производная 12mLt\tfrac12m_L\sqrt t с наклоном ½ (линейный приток в канале); и, когда «почувствован» торец, переход на параллельную прямую вдвое выше, mLtm_L\sqrt t. Торец делает с линейным притоком то же, что разлом с радиальным: удваивает его, и по тому же закону.

Скважина в канале, закрытом с одного конца: поздний линейный приток
Δp=2mLt+qBμ2πkh (S+Sch)−qBμkhW L3,Δp′=mLt\Delta p = 2m_L\sqrt t + \frac{qB\mu}{2\pi kh}\,(S + S_{ch}) - \frac{qB\mu}{khW}\,L_3,\qquad \Delta p' = m_L\sqrt t
где:
  • Δp\Delta pснижение забойного давления pi−pwfp_i - p_{wf}, Па
  • mLm_Lнаклон прямой Δp\Delta p — t\sqrt t открытого канала той же ширины, qBhWμπkϕct\dfrac{qB}{hW}\sqrt{\dfrac{\mu}{\pi k\phi c_t}}, Па/√с
  • S, SchS,\ S_{ch}скин скважины и геометрический скин канала ln⁡W2πrw−ln⁡sin⁡πL1W\ln\frac{W}{2\pi r_w} - \ln\sin\frac{\pi L_1}{W} (геометрический скин канала), —
  • W, L3W,\ L_3ширина канала и расстояние от скважины до торца, м
  • q, B, μ, k, hq,\ B,\ \mu,\ k,\ hдебит в стандартных условиях, объёмный коэффициент, вязкость, проницаемость, толщина, м³/с, м³/м³, Па·с, м², м

Единицы: СИ

Справедливо, если
  • стенки и торец — непроницаемые плоскости через всю толщину однородного пласта; в другую сторону канал не ограничен;
  • торец дальше ширины канала (L3>WL_3 > W), и ηt≫L32\eta t \gg L_3^2 — торец давно «почувствован»; влияние ствола закончилось, дебит постоянен.

Наклон удвоен: флюид идёт к скважине только с открытой стороны, по одной «трубе» сечением WhWh вместо двух. Последнее слагаемое — перепад, который дал бы весь дебит, протекая по каналу на длине L3L_3: чем дальше торец, тем ниже прямая. Подход к этой прямой медленный — избыток над ней убывает как L32/πηtL_3^2/\sqrt{\pi\eta t}, а не как 1/t1/t, поэтому поздние точки ещё долго лежат выше неё.

Производные скважины посередине открытого канала и каналов, закрытых на разных расстояниях0,010,11101001000131030100300ηt / W²Δp′ / mₑоранжевая — L₃ = Wбирюзовая — L₃ = 3Wтёмная — L₃ = 10Wсерая — открытый канал

График в двойном логарифмическом масштабе: по горизонтали eta t делённое на W в квадрате от 0,01 до 1000, по вертикали отношение производной к полке от 0,8 до 300. Все кривые начинаются на уровне 1, около 0,05 поднимаются и после 0,1–0,2 выходят на нижнюю штриховую прямую с наклоном одна вторая. Кривая открытого канала идёт по ней до конца, к 1000 — около 112. Кривая для торца на расстоянии W отходит вверх около 0,3 и к 10 выходит на верхнюю штриховую прямую. Кривая для 3W отходит около 2 и выходит на верхнюю прямую к 100. Кривая для 10W отходит около 20 и к 1000 почти выходит на верхнюю прямую.

РасчётРасчётные производные скважины посередине канала шириной WW, без влияния ствола: открытый канал и каналы, закрытые в L3L_3 = WW, 3W3W и 10W10W от скважины. По вертикали — отношение производной к полке mem_e; штриховые прямые — линейный приток открытого канала и вдвое выше. Обе шкалы логарифмические.

Расстояние до торца находят тремя способами. Надёжнее всего — по одной точке производной: множитель торца 1+e−L32/(ηt)1 + \e^{-L_3^2/(\eta t)} точен в любой момент, и отношение производной к продолжению линейной прямой открытого канала сразу даёт L3L_3. Вторая возможность — пересечение двух прямых на графике Δp\Delta p — t\sqrt t, как пересечение двух полулогарифмических прямых у разлома (Непроницаемый разлом у скважины: метод изображений); но поздняя прямая формируется очень долго, и без поправки расстояние получается заниженным. Третья — подбор модели ко всей записи.

Расстояние до торца канала: по точке производной и по пересечению прямых Δp — √t
L3=ηt ln⁡1r−1,r=Δp′(t)12mLt;L3=η txπL_3 = \sqrt{\eta t\,\ln\frac{1}{r - 1}},\quad r = \frac{\Delta p'(t)}{\tfrac12 m_L\sqrt t};\qquad L_3 = \sqrt{\frac{\eta\,t_x}{\pi}}
где:
  • rrотношение производной в момент tt к продолжению прямой линейного притока открытого канала, 1<r<21 < r < 2, —
  • Δp′(t)\Delta p'(t)логарифмическая производная давления в момент tt, Па
  • txt_xвремя пересечения прямой открытого канала (наклон mLm_L) и поздней прямой (наклон 2mL2m_L) на графике Δp\Delta p — t\sqrt t, с
  • η\etaпьезопроводность k/(ϕμct)k/(\phi\mu c_t), м²/с

Единицы: СИ

Справедливо, если
  • в записи есть линейный приток открытого канала — по нему проводят прямую 12mLt\tfrac12m_L\sqrt t на производной и прямую с наклоном mLm_L на графике Δp\Delta p — t\sqrt t;
  • торец дальше ширины канала; проницаемость известна по полке радиального притока;
  • вторая формула верна, если поздняя прямая действительно достигнута: при прямой с наклоном 2mL2m_L через последний замер в момент tendt_{end} она даёт Lx=L3−L32/πηtendL_x = L_3 - L_3^2/\sqrt{\pi\eta t_{end}} (поправка — в примере).

Первая формула — та же, что для ступени производной у разлома, только отношение берут к прямой с наклоном ½, а не к полке. Вторая получается приравниванием двух прямых: разность их отрезков равна qBμL3/(khW)qB\mu L_3/(khW), а mL=qBμkhWη/πm_L = \tfrac{qB\mu}{khW}\sqrt{\eta/\pi}.

В других системах единиц (2)
Метрические промысловые единицы
L3=0,0596k tϕμctln⁡1r−1,L3=0,0336k txϕμctL_3 = 0{,}0596\sqrt{\frac{k\,t}{\phi\mu c_t}\ln\frac{1}{r - 1}},\qquad L_3 = 0{,}0336\sqrt{\frac{k\,t_x}{\phi\mu c_t}}

L3L_3 — м, kk — мД, tt и txt_x — ч, μ\mu — мПа·с, ctc_t — 1/МПа.

0,0596=3,553⋅10−30{,}0596 = \sqrt{3{,}553\cdot10^{-3}}, как в формуле ступени у разлома (Непроницаемый разлом у скважины: метод изображений); 0,0336=3,553⋅10−3/π0{,}0336 = \sqrt{3{,}553\cdot10^{-3}/\pi}.

Промысловые единицы (oilfield)
L3=0,01624k tϕμctln⁡1r−1,L3=0,00916k txϕμctL_3 = 0{,}01624\sqrt{\frac{k\,t}{\phi\mu c_t}\ln\frac{1}{r - 1}},\qquad L_3 = 0{,}00916\sqrt{\frac{k\,t_x}{\phi\mu c_t}}

L3L_3 — ft, kk — мД, tt и txt_x — ч, μ\mu — сП, ctc_t — 1/psi.

0,01624=2,637⋅10−40{,}01624 = \sqrt{2{,}637\cdot10^{-4}}; 0,00916=2,637⋅10−4/π0{,}00916 = \sqrt{2{,}637\cdot10^{-4}/\pi}.

Пробная эксплуатация скважины в русловом теле

Разберём условную запись. Разведочная скважина вскрыла песчаник руслового тела; по сейсмике и соседним скважинам тело узкое и вытянутое, но неизвестно, открыто ли оно в обе стороны. Скважину пустили открытием забойного клапана с постоянным дебитом 150 м³/сут и записывали давление около 26 суток.

Снижение давления в скважине руслового тела (условные данные, выборка)
tt
36 с
t\sqrt t, √ч
0,10
Δp\Delta p, МПа
0,4148
Δp′\Delta p', МПа
0,0398
tt
6 мин
t\sqrt t, √ч
0,32
Δp\Delta p, МПа
0,5031
Δp′\Delta p', МПа
0,0377
tt
1,00 ч
t\sqrt t, √ч
1,00
Δp\Delta p, МПа
0,5970
Δp′\Delta p', МПа
0,0553
tt
2,00 ч
t\sqrt t, √ч
1,41
Δp\Delta p, МПа
0,6419
Δp′\Delta p', МПа
0,0775
tt
6,31 ч
t\sqrt t, √ч
2,51
Δp\Delta p, МПа
0,7615
Δp′\Delta p', МПа
0,1380
tt
10,0 ч
t\sqrt t, √ч
3,16
Δp\Delta p, МПа
0,8324
Δp′\Delta p', МПа
0,1753
tt
20,0 ч
t\sqrt t, √ч
4,47
Δp\Delta p, МПа
0,9785
Δp′\Delta p', МПа
0,2663
tt
50,1 ч
t\sqrt t, √ч
7,08
Δp\Delta p, МПа
1,3207
Δp′\Delta p', МПа
0,5280
tt
100 ч
t\sqrt t, √ч
10,00
Δp\Delta p, МПа
1,7898
Δp′\Delta p', МПа
0,8742
tt
200 ч
t\sqrt t, √ч
14,13
Δp\Delta p, МПа
2,5474
Δp′\Delta p', МПа
1,3703
tt
398 ч
t\sqrt t, √ч
19,95
Δp\Delta p, МПа
3,7065
Δp′\Delta p', МПа
2,0535
tt
631 ч
t\sqrt t, √ч
25,12
Δp\Delta p, МПа
4,7740
Δp′\Delta p', МПа
—

Δp=pi−pwf\Delta p = p_i - p_{wf}, pip_i = 24 МПа. Данные рассчитаны по модели однородного пласта со стволом, скином и тремя непроницаемыми границами (две стенки и торец) и округлены до 0,0001 МПа; параметры модели раскрыты в примере. Производная — по всем 59 замерам (10 на цикл) с интервалом 0,2 цикла; в последних двух замерах её нет.

Диагностический график скважины в русловом теле, закрытом с одного конца0,0010,010,111010010000,1110Время после пуска t, чΔp и Δp′, МПаточки — Δpоранжевая — Δp′штрих — полка и прямые наклона ½mₑ = 0,0376 МПа

График в двойном логарифмическом масштабе: время от 0,001 до 1000 часов, давление от 0,02 до 10 МПа. Точки снижения давления растут от 0,3 МПа в начале до 0,6 МПа к часу и до 4,8 МПа к 631 часу. Производная после короткого горба влияния ствола лежит на горизонтальной штриховой полке 0,038 МПа примерно от 0,03 до 0,25 часа, затем поднимается и от 1 до 10 часов идёт вдоль нижней штриховой прямой с наклоном одна вторая. После 10 часов производная отходит от неё вверх: к 100 часам 0,87 МПа, к 400 часам 2,05 МПа, и приближается к верхней штриховой прямой, вдвое более высокой.

РасчётДиагностический график записи: снижение давления (точки) и производная по замерам (линия). Штриховые линии — полка mem_e, прямая линейного притока открытого канала 12mLt\tfrac12m_L\sqrt t и вдвое более высокая прямая mLtm_L\sqrt t.
Снижение давления против корня из времени: прямая открытого канала и прямая с удвоенным наклоном0510152025012345√t, √чСнижение давления Δp, МПа√tₓ = 10,8бирюзовая — прямая открытого каналаоранжевая — вдвое круче, через последний замерсерая — асимптота модели

График в обычных координатах: по горизонтали корень из времени от 0 до 26, по вертикали снижение давления от 0 до 5 МПа. Точки поднимаются от 0,4 МПа и с корня из времени около 1,4 идут по бирюзовой штриховой прямой с малым наклоном, затем загибаются вверх. Оранжевая штриховая прямая с вдвое большим наклоном проходит через последнюю точку 4,77 МПа при 25,1 и пересекает первую около 10,8, где отмечена точка. Серая штриховая асимптота модели параллельна оранжевой, идёт ниже неё на 0,23 МПа и пересекает первую прямую около 12,9; точки между 11 и 25 лежат выше обеих поздних прямых.

РасчётТа же запись в координатах Δp\Delta p — t\sqrt t (точки). Первая прямая проведена через замеры 2 и 6,3 ч (линейный приток открытого канала), вторая — с вдвое большим наклоном через последний замер; отмечено их пересечение. Светлая штриховая — асимптота модели, к которой запись подходит лишь при очень долгой работе: поздние точки ещё лежат выше неё, и прямая через последний замер пересекает первую раньше — при tx\sqrt{t_x} = 10,8 вместо 12,9.
Расчётный примерПроницаемость, ширина канала и расстояние до торца
Исходные данные (условные)
Скважина и пласт
qq = 150 м³/сут; BB = 1,2; μ\mu = 0,8 мПа·с; hh = 12 м; ϕ\phi = 0,22; ct=1,4⋅10−3c_t = 1{,}4\cdot10^{-3} МПа⁻¹; rwr_w = 0,1 м
Полка производной
среднее от 0,1 до 0,2 ч: mem_e = 0,03759 МПа; замер через 0,1 ч: Δp\Delta p = 0,5031 МПа
Линейный приток открытого канала
замеры 1,995 ч — 0,6419 МПа, 6,310 ч — 0,7615 МПа; производная на этом участке, делённая на t\sqrt t, в среднем 0,05489 МПа/√ч
Поздние данные
производная через 100 ч — 0,8742 МПа; последний замер через 631 ч — 4,7740 МПа
Решение
  1. Проводимость: qBμ=150⋅1,2⋅0,8=144qB\mu = 150\cdot1{,}2\cdot0{,}8 = 144; kh=0,933⋅144/0,03759=3574kh = 0{,}933\cdot144/0{,}03759 = 3574 мД·м, k=297,8k = 297{,}8 мД. Скин: ln⁡ktϕμctrw2=16,308\ln\dfrac{kt}{\phi\mu c_tr_w^2} = 16{,}308 при tt = 0,1 ч; S=12 (0,5031/0,03759−16,308+4,831)=0,95S = \tfrac12\,(0{,}5031/0{,}03759 - 16{,}308 + 4{,}831) = 0{,}95 (Диагностика по производной давления).
  2. Пьезопроводность: η=3,553⋅10−3⋅297,8/(0,22⋅0,8⋅1,4⋅10−3)=4294\eta = 3{,}553\cdot10^{-3}\cdot297{,}8/(0{,}22\cdot0{,}8\cdot1{,}4\cdot10^{-3}) = 4294 м²/ч.
  3. Первая прямая Δp\Delta p — t\sqrt t: mL=(0,7615−0,6419)/(2,5119−1,4125)=0,10879m_L = (0{,}7615 - 0{,}6419)/(2{,}5119 - 1{,}4125) = 0{,}10879 МПа/√ч. Ширина канала: μ/(kϕct)=2,9533\sqrt{\mu/(k\phi c_t)} = 2{,}9533; W=0,394⋅180/(12⋅0,10879)⋅2,9533=160W = 0{,}394\cdot180/(12\cdot0{,}10879)\cdot2{,}9533 = 160 м (ширина канала).
  4. Торец по точке производной: продолжение прямой открытого канала, проведённой по самой производной (так ошибка численного дифференцирования сокращается), через 100 ч — 0,05489⋅1000{,}05489\cdot\sqrt{100}; r=0,8742/(0,05489⋅10)=1,593r = 0{,}8742/(0{,}05489\cdot10) = 1{,}593; ln⁡[1/(1,593−1)]=0,523\ln[1/(1{,}593 - 1)] = 0{,}523; L3=4294⋅100⋅0,523=474L_3 = \sqrt{4294\cdot100\cdot0{,}523} = 474 м.
  5. По пересечению прямых: отрезок первой b1=0,6419−0,10879⋅1,4125=0,4882b_1 = 0{,}6419 - 0{,}10879\cdot1{,}4125 = 0{,}4882 МПа; прямая с наклоном 2mL2m_L через последний замер: b2=4,7740−2⋅0,10879⋅25,119=−0,6914b_2 = 4{,}7740 - 2\cdot0{,}10879\cdot25{,}119 = -0{,}6914 МПа. Пересечение: tx=(b1−b2)/mL=10,843\sqrt{t_x} = (b_1 - b_2)/m_L = 10{,}843, tx=117,6t_x = 117{,}6 ч; Lx=0,0336297,8⋅117,6/(0,22⋅0,8⋅1,4⋅10−3)=401L_x = 0{,}0336\sqrt{297{,}8\cdot117{,}6/(0{,}22\cdot0{,}8\cdot1{,}4\cdot10^{-3})} = 401 м.
  6. Поправка на медленный подход поздней прямой: X=πηtend=π⋅4294⋅631=2918X = \sqrt{\pi\eta t_{end}} = \sqrt{\pi\cdot4294\cdot631} = 2918 м; L3=X2(1−1−4Lx/X)=480L_3 = \tfrac{X}{2}\left(1 - \sqrt{1 - 4L_x/X}\right) = 480 м.
Ответ: k≈297,8k \approx 297{,}8 мД, скин 0,95, канал шириной около 160 м; торец — в 474 м по точке производной и в 480 м по пересечению прямых с поправкой (без поправки — 401 м).

Что это значит. Запись рассчитана по модели с kk = 300 мД, SS = 1, каналом шириной 160 м и торцом в 480 м. Проницаемость занижена примерно на 1 %: полка ещё чуть приподнята хвостом влияния ствола. Точка производной дала расстояние до торца с ошибкой около 1 % уже через 100 ч, а пересечение прямых без поправки занизило его на 16 %: за 26 сут поздние точки ещё не легли на асимптоту. Если бы канал считали открытым и взяли наклон последних точек, ширина вышла бы около 80 м — почти вдвое меньше. С какой стороны от скважины торец, запись не говорит: это дают сейсмика и геологическая модель. Данные примера условные.

Скважина у самой границы: удвоенный дебит

Если скважина стоит в нескольких метрах от разлома, её отражение почти совпадает с ней самой. Как только возмущение пройдёт расстояние до границы (время порядка δ2/η\delta^2/\eta — секунды или минуты), они работают как одна скважина с удвоенным дебитом: флюид подходит только с одной стороны, и сопротивление притоку вдвое больше (полурадиальный приток). Это позволяет свести «пристенные» системы к более простым без новой эталонной кривой: скважина у одной стенки канала шириной WW — это скважина с удвоенным дебитом посередине открытого канала шириной 2W2W, скважина у торца канала — скважина с удвоенным дебитом в открытом канале той же ширины, скважина в углу — четырёхкратный дебит в неограниченном пласте.

Скважина у непроницаемой границы как скважина с удвоенным дебитом
pwD=2 pD∗(tD)−ln⁡2δrw+Sp_{wD} = 2\,p_D^{*}(t_D) - \ln\frac{2\delta}{r_w} + S
где:
  • pwDp_{wD}безразмерное снижение забойного давления скважины у границы, —
  • pD∗(tD)p_D^{*}(t_D)безразмерное решение без скина для той же системы без ближней границы (все остальные границы и их отражения сохраняются), —
  • δ\deltaрасстояние от оси скважины до ближней границы, м
  • SSскин скважины — входит один раз, потому что относится к настоящему дебиту, —

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

Справедливо, если
  • граница прямая и непроницаемая; δ\delta мало по сравнению с остальными размерами системы;
  • время больше нескольких δ2/η\delta^2/\eta: отражение на расстоянии 2δ2\delta уже сливается со скважиной (вывод — в дополнительном разделе).

Удвоенный дебит удваивает весь «пластовый» перепад, а скин — нет. Слагаемое −ln⁡(2δ/rw)-\ln(2\delta/r_w) — поправка на то, что отражение стоит не в самой скважине, а в 2δ2\delta от неё: эквивалентная скважина с удвоенным дебитом имеет радиус 2δrw\sqrt{2\delta r_w} — среднее геометрическое радиуса скважины и расстояния до отражения. Это та же поправка, что в формуле полурадиального притока.

Отсюда ловушка. Скважина у торца канала сначала даёт полку 2me2m_e, затем линейный приток с удвоенным наклоном 2mL2m_L. Если принять её за обычную скважину в открытом канале, полка даст вдвое заниженную проводимость, а ширина канала по наклону выйдет в 2\sqrt2 раз меньше настоящей: наклон вдвое больше, но и 1/k\sqrt{1/k} в формуле ширины — в 2\sqrt2 раз больше. Запись при этом выглядит безупречно — её выдаёт только несогласие с керном, соседними скважинами и геологией.

Восстановление давления в полуограниченной системе

Кривая восстановления давления (КВД) — суперпозиция пуска и остановки, и каждое отражение тоже «останавливается» вместе со скважиной (Анализ кривой восстановления давления). Обобщение формулы для одного разлома (КВД у разлома) на любой набор отражений:

КВД скважины в системе с отражениями
pi−pws=me[ln⁡tp+ΔtΔt+∑j(E1 ⁣(dj24η (tp+Δt))−E1 ⁣(dj24η Δt))]p_i - p_{ws} = m_e\left[\ln\frac{t_p + \Delta t}{\Delta t} + \sum_j\left(E_1\!\left(\frac{d_j^2}{4\eta\,(t_p + \Delta t)}\right) - E_1\!\left(\frac{d_j^2}{4\eta\,\Delta t}\right)\right)\right]
где:
  • pwsp_{ws}забойное давление через Δt\Delta t после остановки, Па
  • tp, Δtt_p,\ \Delta tвремя работы с постоянным дебитом и время после остановки, с
  • djd_jрасстояния до отражений, м
  • me, ηm_e,\ \etaуровень полки производной и пьезопроводность, Па, м²/с

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

Справедливо, если
  • до пуска давление в пласте было pip_i; дебит постоянен в течение tpt_p; для самой скважины — логарифмическое приближение, влияние ствола закончилось;
  • для каналов сумма бесконечна, но сходится: дальние отражения дают одинаковые слагаемые при tp+Δtt_p + \Delta t и Δt\Delta t, и разность исчезает.

Результат зависит теперь от двух безразмерных времён, ηtp/L2\eta t_p/L^2 и ηΔt/L2\eta\Delta t/L^2, а не от одного, как у кривой падения давления (КПД): граница видна на КВД, только если возмущение успело до неё дойти за время работы или успевает за время записи. При Δt→∞\Delta t \to \infty все слагаемые обращаются в ноль, и pws→pip_{ws} \to p_i: полуограниченная система восстанавливается до начального давления. Продолжение ранней прямой Хорнера даёт p∗p^* ниже pip_i — на me∑jE1(dj2/(4ηtp))m_e\sum_j E_1\big(d_j^2/(4\eta t_p)\big).

У канала поздние точки КВД ложатся на прямую в координатах pwsp_{ws} — (tp+Δt−Δt)(\sqrt{t_p + \Delta t} - \sqrt{\Delta t}) (КВД в канале); у закрытого с одного конца — с вдвое большим наклоном 2mL2m_L, когда и tpt_p, и Δt\Delta t уже на участке удвоенного линейного притока. Производную такой КВД строят с учётом времени работы (Влияние времени работы на анализ КВД). Чаще, однако, торец видят по КВД после короткой работы — испытания пласта на трубах, — когда до него возмущение за время притока почти не дошло. Тогда след торца — небольшое отставание давления от КВД открытого канала, и его наибольшую величину можно оценить заранее:

Наибольший след торца канала на КВД после короткой работы
δpmax⁡=e−1/222π⋅qB tpϕcth WL3=0,121 qB tpϕcth WL3,Δtmax⁡≈2L32η\delta p_{\max} = \frac{\e^{-1/2}}{2\sqrt{2\pi}}\cdot\frac{qB\,t_p}{\phi c_t h\,WL_3} = 0{,}121\,\frac{qB\,t_p}{\phi c_t h\,WL_3},\qquad \Delta t_{\max} \approx \frac{2L_3^2}{\eta}
где:
  • δpmax⁡\delta p_{\max}наибольшее отставание давления КВД от КВД того же канала без торца, Па
  • qB tpqB\,t_pобъём, отобранный за время работы, в пластовых условиях, м³
  • ϕ, ct, h\phi,\ c_t,\ hпористость, общая сжимаемость, толщина, —, 1/Па, м
  • W, L3W,\ L_3ширина канала и расстояние до торца, м
  • Δtmax⁡\Delta t_{\max}время после остановки, когда отставание наибольшее, с

Единицы: СИ

Справедливо, если
  • за время работы торец почти не почувствовался: ηtp≪L32\eta t_p \ll L_3^2 (практически годится при L32>3ηtpL_3^2 > 3\eta t_p);
  • торец дальше ширины канала, и к моменту Δtmax⁡\Delta t_{\max} в канале давно линейный приток;
  • однородный пласт; ни дрейфа манометра, ни работы соседних скважин.

Сравните с одним разломом в радиальном пласте: там след убывает как 1/L21/L^2, здесь — как 1/(WL3)1/(WL_3). В канале возмущение не растекается по плоскости, и торец сказывается на давлении гораздо сильнее, чем разлом на том же расстоянии; приходит же сигнал так же поздно — время растёт как L32L_3^2 (вывод — в дополнительном разделе).

В других системах единиц (1)
Метрические промысловые единицы
δpmax⁡=qB tp198 ϕcth WL3\delta p_{\max} = \frac{qB\,t_p}{198\,\phi c_t h\,WL_3}

qq — м³/сут, tpt_p — ч, ctc_t — 1/МПа, hh, WW, L3L_3 — м; δpmax⁡\delta p_{\max} — МПа.

198=24/0,121198 = 24/0{,}121 (24 — часов в сутках).

Отставание КВД от КВД открытого канала из-за торца после притока разной длительности1101001000020406080Время после остановки Δt, чОтставание давления, кПаоранжевая — приток 24 чбирюзовая — 12 чсерая — 6 чточки — оценка по формуле

График отставания давления от 0 до 70 килопаскалей в зависимости от времени после остановки от 1 до 1000 часов на логарифмической шкале. Три кривые поднимаются, проходят максимум и снова убывают; кривая 24-часового притока начинается не с нуля, а около 10 килопаскалей: торец успел сказаться ещё во время притока. Для притока 6 часов максимум около 19 килопаскалей около 104 часов, для 12 часов — около 38 килопаскалей около 101 часа, для 24 часов — около 77 килопаскалей около 95 часов. Точки оценок по формуле лежат у вершин кривых при 107 часах.

РасчётРасчётное отставание давления КВД от КВД открытого канала после притока 6, 12 и 24 ч в русловом теле из примера выше (ширина 160 м, торец в 480 м). Точки — оценка δpmax⁡\delta p_{\max} по формуле при Δt=2L32/η\Delta t = 2L_3^2/\eta. Шкала времени логарифмическая.
Расчётный примерУвидит ли испытание на трубах торец руслового тела
Исходные данные (условные)
Испытание
приток 12 ч с дебитом 150 м³/сут (BB = 1,2), затем КВД; манометр различает 0,5 кПа
Пласт
как в примере выше: kk = 300 мД, ϕ\phi = 0,22, μ\mu = 0,8 мПа·с, ct=1,4⋅10−3c_t = 1{,}4\cdot10^{-3} МПа⁻¹, hh = 12 м
Канал
ширина 160 м, торец в 480 м от скважины
Решение
  1. Пьезопроводность η=4326\eta = 4326 м²/ч; ηtp=51 912\eta t_p = 51\,912 м² против L32L_3^2 = 230 400 м² — за время притока торец почти не почувствовался (ηtp/L32=0,23\eta t_p/L_3^2 = 0{,}23).
  2. Упругая ёмкость участка канала до торца: ϕcthWL3=0,22⋅1,4⋅10−3⋅12⋅160⋅480=283,85\phi c_thWL_3 = 0{,}22\cdot1{,}4\cdot10^{-3}\cdot12\cdot160\cdot480 = 283{,}85 м³/МПа; qBtp=150⋅1,2⋅12=2160qBt_p = 150\cdot1{,}2\cdot12 = 2160 (м³/сут·ч).
  3. Наибольший след торца: δpmax⁡=2160/(198⋅283,85)=0,0384\delta p_{\max} = 2160/(198\cdot283{,}85) = 0{,}0384 МПа, то есть 38,4 кПа, примерно через 2L32/η=1072L_3^2/\eta = 107 ч после остановки.
  4. Расчёт по полной сумме отражений даёт 38,4 кПа через 101 ч — формула верна.
  5. Для сравнения: разлом на тех же 480 м в радиальном пласте оставил бы след qBtp/(820 ϕcthL2)=3,1qBt_p/(820\,\phi c_thL^2) = 3{,}1 кПа (предел обнаружения разлома) — примерно в 12 раз меньше.
Ответ: Торец в 480 м отклонит КВД примерно на 38,4 кПа — в десятки раз больше порога манометра, — но максимум наступит лишь через 107 ч после остановки.

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

Что запись различить не может

  • Направления. Запись зависит только от расстояний до отражений. С какой стороны торец, в какую сторону вытянут канал, запись не скажет — это дают сейсмика, геологическая модель и соседние скважины.
  • Торец или ухудшение пласта. Если по одну сторону от скважины канал дальше сужается, утончается или теряет проницаемость, наклон t\sqrt t тоже растёт: на поздних временах — в 2/(1+k2/k)2/(1 + \sqrt{k_2/k}) раза, если по ту сторону проницаемость k2k_2, а пористость и сечение те же (одномерный приток с двух сторон к скважине). При проницаемости вдвое ниже наклон больше примерно на 17 %, в 10 раз ниже — в 1,5 раза, в 100 раз ниже — в 1,8 раза. Такую запись легко принять за ещё не закончившееся удвоение от торца; различить их помогают длина записи, геология и соседние скважины.
  • Торец, пропускающий флюид. Если поперечный разлом частично проводящий, производная поднимается к удвоенной прямой и затем возвращается к прямой открытого канала (Частично проводящие разломы).
  • Замкнутый блок. Вытянутый замкнутый блок до выхода на наклон 1 тоже даёт наклон ½; второй торец проявится поздним подъёмом с наклоном 1 и КВД, которая выходит не на pip_i, а на среднее давление (Исследование на ограниченность пласта).
  • Незаконченный переход. На производной торец сказывается по тому же закону, что разлом, но поздняя прямая Δp\Delta p — t\sqrt t подходит к асимптоте медленно, и расстояние по второй прямой без поправки занижено на долю около L3/πηtL_3/\sqrt{\pi\eta t}: чтобы ошибка была меньше 5 %, ηt\eta t должно превысить L32L_3^2 больше чем в сто раз. Точка производной и подбор модели надёжнее.
  • Прямые границы — модель. Русловые тела изгибаются, меняют ширину и толщину; расстояния по простой модели — оценки, которые сверяют с картой.

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

Полуограниченная система
пласт, ограниченный непроницаемыми границами не со всех сторон: полуплоскость, угол, канал, канал с одним закрытым концом; псевдоустановившегося режима нет.
Канал, закрытый с одного конца
пласт между двумя параллельными границами, перегороженный третьей с одной стороны (U-образная система разломов); поздний линейный приток вдвое сильнее, чем в открытом канале.
Характерная длина LL
размер, которому пропорциональны все расстояния до отражений; поздняя часть решения зависит только от ηt/L2\eta t/L^2.
Правило произведения
в прямоугольной системе производная, отнесённая к полке, равна произведению множителей по двум направлениям.
Множитель направления
сумма exp⁡(−ξ2/(4ηt))\exp(-\xi^2/(4\eta t)) по одномерным отражениям вдоль одного направления: 1 — открыто, 1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)} — одна граница, полоса — две.
Торец канала L3L_3
расстояние от скважины до закрытого конца канала; находят по точке производной или по пересечению прямых Δp\Delta p — t\sqrt t.
Удвоенный дебит
приём для скважины у самой границы: её отражение сливается с ней, и она работает как скважина с удвоенным дебитом в системе без этой границы; скин не удваивают.

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

Неверно

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

На самом деле

Сначала проверяют, не перешла ли производная на параллельную прямую вдвое выше. В примере урока ширина по поздним точкам вышла бы около 80 м вместо 160.

Неверно

«Вторая прямая на графике Δp\Delta p — t\sqrt t — это последние точки записи».

На самом деле

Поздняя прямая у торца формируется медленно: избыток над ней убывает как 1/t1/\sqrt t. Прямая с наклоном 2mL2m_L через последний замер в примере урока дала торец в 401 м вместо 480; поправка или точка производной это исправляют.

Неверно

«Давление после остановки не дошло до начального за неделю — значит, блок замкнут и истощается».

На самом деле

В канале, закрытом с одного конца, давление восстанавливается до pip_i, но медленно: торец задерживает его на время порядка 2L32/η2L_3^2/\eta. Замкнутость доказывает только выход на среднее давление ниже pip_i и наклон 1 при работе.

Неверно

«У скважины рядом с разломом удваиваются и дебит, и скин».

На самом деле

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

Неверно

«Торец далеко — короткое испытание его не увидит, как не увидит далёкий разлом».

На самом деле

В канале след торца убывает как 1/(WL3)1/(WL_3), а не как 1/L21/L^2: в примере урока 12-часовой приток дал бы на КВД след около 38 кПа — в 12 раз больше, чем разлом на том же расстоянии в радиальном пласте. Ограничение — время: максимум приходит через 2L32/η2L_3^2/\eta.

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

Выводы: правило произведения, поздняя прямая, удвоенный дебит и след торца

Правило произведения. Отражение в прямоугольной системе получают, отражая скважину несколько раз в границах вдоль xx и вдоль yy; его сдвиг — пара (ξi,ζl)(\xi_i, \zeta_l), где ξi\xi_i пробегает одномерные отражения по xx (включая 0), а ζl\zeta_l — по yy. Все пары встречаются ровно по разу, и ∑i,lexp⁡(−(ξi2+ζl2)/(4ηt))=∑iexp⁡(−ξi2/(4ηt))⋅∑lexp⁡(−ζl2/(4ηt))\sum_{i,l}\exp\big(-(\xi_i^2 + \zeta_l^2)/(4\eta t)\big) = \sum_i\exp(-\xi_i^2/(4\eta t))\cdot\sum_l\exp(-\zeta_l^2/(4\eta t)). Производная по ln⁡t\ln t от решения для линейного источника с отражением на расстоянии dd равна mee−d2/(4ηt)m_e\e^{-d^2/(4\eta t)}, отсюда правило. Для полосы шириной WW ряд по ζl\zeta_l (шаг 2W2W) по формуле суммирования Пуассона переходит в ряд по косинусам, главный член которого 4πηt/W\sqrt{4\pi\eta t}/W: множитель полосы растёт как корень из времени.

Поздняя прямая и пересечение. Когда ηt≫W2\eta t \gg W^2, ряд отражений в стенках действует как плоский источник, и одномерное решение для ряда на расстоянии 2L32L_3 даёт вклад qBμkhW[ηt/π e−a−L3 erfca]\tfrac{qB\mu}{khW}\big[\sqrt{\eta t/\pi}\,\e^{-a} - L_3\,\mathrm{erfc}\sqrt a\big], a=L32/(ηt)a = L_3^2/(\eta t). При малых aa это mLt−qBμkhW L3+qBμkhW⋅L32πηt+…m_L\sqrt t - \tfrac{qB\mu}{khW}\,L_3 + \tfrac{qB\mu}{khW}\cdot\tfrac{L_3^2}{\sqrt{\pi\eta t}} + \dots — асимптота из основного раздела и медленно убывающий избыток над ней. Приравнивая прямые mLt+bm_L\sqrt t + b и 2mLt+b−qBμkhWL32m_L\sqrt t + b - \tfrac{qB\mu}{khW}L_3 и учитывая mL=qBμkhWη/πm_L = \tfrac{qB\mu}{khW}\sqrt{\eta/\pi}, получаем tx=L3π/η\sqrt{t_x} = L_3\sqrt{\pi/\eta}. Прямая с наклоном 2mL2m_L через замер в момент tendt_{end} поднята на избыток, и пересечение даёт Lx=L3−L32/XL_x = L_3 - L_3^2/X, X=πηtendX = \sqrt{\pi\eta t_{end}}; решение этого квадратного уравнения — поправка из примера.

Удвоенный дебит. Отражение скважины в ближней границе стоит на расстоянии 2δ2\delta и при tD≫δD2t_D \gg \delta_D^2 даёт 12E1(δD2/tD)≈12[ln⁡(tD/δD2)−0,5772]\tfrac12E_1(\delta_D^2/t_D) \approx \tfrac12\big[\ln(t_D/\delta_D^2) - 0{,}5772\big]. Вместе со скважиной: ln⁡tD+0,809−ln⁡(2δD)\ln t_D + 0{,}809 - \ln(2\delta_D), то есть 2⋅12[ln⁡(tD/(2δD))+0,809]2\cdot\tfrac12\big[\ln\big(t_D/(2\delta_D)\big) + 0{,}809\big] — удвоенное решение для скважины радиусом 2δrw\sqrt{2\delta r_w}. Все остальные отражения тоже стоят парами в почти совпадающих точках и удваиваются.

След торца на КВД. Вклад зеркального ряда при КВД — разность g(tp+Δt)−g(Δt)g(t_p + \Delta t) - g(\Delta t) одномерного решения из второго абзаца. При tp≪Δtt_p \ll \Delta t она близка к tp  dg/ dtt_p\,\dd g/\dd t, а  dg/ dt=qBμkhW⋅12η/(πt) e−L32/(ηt)\dd g/\dd t = \tfrac{qB\mu}{khW}\cdot\tfrac12\sqrt{\eta/(\pi t)}\,\e^{-L_3^2/(\eta t)}. Максимум функции t−1/2e−L32/(ηt)t^{-1/2}\e^{-L_3^2/(\eta t)} — при t=2L32/ηt = 2L_3^2/\eta; подстановка даёт δpmax⁡=e−1/222π⋅qB tpϕcthWL3\delta p_{\max} = \tfrac{\e^{-1/2}}{2\sqrt{2\pi}}\cdot\tfrac{qB\,t_p}{\phi c_thWL_3}, потому что qBμηkh=qBϕcth\tfrac{qB\mu\eta}{kh} = \tfrac{qB}{\phi c_th}.

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

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

  1. 1

    После полки радиального притока производная около цикла идёт с наклоном ½, затем поднимается и выходит на параллельную прямую примерно вдвое выше. Какая система самая вероятная?

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

    Правильный ответ: Б. Канал, закрытый с одного конца

    Наклон ½ — линейный приток в канале. Переход на прямую вдвое выше при том же наклоне — множитель торца 1+e−L32/(ηt)1 + \e^{-L_3^2/(\eta t)}, выросший до 2. Разлом и угол дают горизонтальные полки, граница постоянного давления опускает производную.

  2. 2

    В записи скважины в канале через 200 ч производная в 1,42 раза выше продолжения прямой линейного притока открытого канала; η\eta = 2000 м²/ч. Как далеко торец?

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

    ln⁡[1/(1,42−1)]=ln⁡2,381=0,868\ln[1/(1{,}42 - 1)] = \ln 2{,}381 = 0{,}868. L3=2000⋅200⋅0,868=347 200≈589L_3 = \sqrt{2000\cdot200\cdot0{,}868} = \sqrt{347\,200} \approx 589 м. Отношение уже заметно больше единицы, поэтому оценка устойчива к шуму производной в несколько процентов.

  3. 3

    По позднему участку записи, где производная уже на удвоенной прямой, наклон t\sqrt t — 0,30 МПа/√ч. Интерпретатор посчитал ширину канала, приняв канал за открытый, и получил 110 м. Какова ширина на самом деле?

    1. А.55 м
    2. Б.Около 155 м
    3. В.110 м
    4. Г.Около 220 м
    Показать ответ

    Правильный ответ: Г. Около 220 м

    На удвоенной прямой наклон равен 2mL2m_L, а ширина обратно пропорциональна наклону. Настоящий mLm_L = 0,15 МПа/√ч, и ширина вдвое больше — около 220 м (если проницаемость определена по полке правильно).

  4. 4

    Скважину пробурили в 8 м от торца канала, и первая полка скрыта влиянием ствола. Интерпретатор обработал запись как скважину в открытом канале. Во сколько раз ошибутся проницаемость и ширина канала?

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

    Скважина у торца работает как скважина с удвоенным дебитом в открытом канале. Видимая полка — 2me2m_e, и проницаемость получится вдвое меньше настоящей. Наклон t\sqrt t тоже удвоен, а в формулу ширины входит 1/k\sqrt{1/k} — в 2\sqrt2 раз больше: Wвидимая=W⋅12⋅2≈0,71 WW_{\text{видимая}} = W\cdot\tfrac12\cdot\sqrt2 \approx 0{,}71\,W. Скин тоже исказится; вскрыть ошибку помогают керн, соседние скважины и геология.

  5. 5

    После долгой работы скважину в канале, закрытом с одного конца, остановили. К какому давлению стремится КВД?

    1. А.К начальному пластовому давлению pip_i
    2. Б.К давлению p∗p^* первой прямой Хорнера
    3. В.К среднему давлению участка канала до торца
    4. Г.К забойному давлению в момент остановки
    Показать ответ

    Правильный ответ: А. К начальному пластовому давлению pip_i

    Канал открыт в одну сторону, объём, охваченный возмущением, растёт без предела, и при Δt→∞\Delta t \to \infty все слагаемые суперпозиции исчезают: pws→pip_{ws} \to p_i. Только медленнее, чем в открытом канале: торец задерживает выравнивание. p∗p^* первой прямой ниже pip_i на вклад отражений.

  6. 6

    Испытание: приток 8 ч с дебитом 100 м³/сут (BB = 1,3), ϕ\phi = 0,18, ct=1,6⋅10−3c_t = 1{,}6\cdot10^{-3} МПа⁻¹, hh = 10 м; канал шириной 200 м закрыт в 600 м от скважины; η\eta = 1500 м²/ч. Каков наибольший след торца на КВД и когда его ждать?

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

    qBtp=100⋅1,3⋅8=1040qBt_p = 100\cdot1{,}3\cdot8 = 1040; ϕcthWL3=0,18⋅1,6⋅10−3⋅10⋅200⋅600=345,6\phi c_thWL_3 = 0{,}18\cdot1{,}6\cdot10^{-3}\cdot10\cdot200\cdot600 = 345{,}6 м³/МПа. δpmax⁡=1040/(198⋅345,6)=0,0152\delta p_{\max} = 1040/(198\cdot345{,}6) = 0{,}0152 МПа, около 15 кПа. Время: 2L32/η=2⋅360 000/1500=4802L_3^2/\eta = 2\cdot360\,000/1500 = 480 ч — 20 сут; условие ηtp=12 000≪L32=360 000\eta t_p = 12\,000 \ll L_3^2 = 360\,000 выполнено. Сигнал заметен, но КВД такой длины после 8-часового притока — редкость: торец вероятнее увидеть при пробной эксплуатации.

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