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

Частично проводящие разломы

Частично проводящий разлом (partially communicating fault, leaky fault) — разлом, через который флюид перетекает, но с добавочным сопротивлением: зона разлома тонкая и плохо проницаемая. В исследовании скважины такой разлом сначала ведёт себя как непроницаемый — производная давления поднимается к удвоенному уровню, — но не доходит до него: по мере того как растёт перепад давления на самом разломе, переток через него усиливается, и производная возвращается к полке радиального притока. Высоту и время этого «горба» задаёт безразмерная проводимость разлома (dimensionless fault transmissibility) α=kfL/(klf)\alpha = k_fL/(kl_f); по ним оценивают проводимость зоны разлома kf/lfk_f/l_f.

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

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

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

  • Модель разлома — тонкая вертикальная зона шириной lfl_f и проницаемостью kfk_f через всю толщину пласта, без собственной упругоёмкости: поток через неё пропорционален скачку давления на ней.
  • Безразмерная проводимость α=(kf/lf)/(k/L)\alpha = (k_f/l_f)/(k/L) сравнивает проводимость зоны разлома с проводимостью пласта на пути от скважины до разлома. α→0\alpha \to 0 — непроницаемый разлом, α→∞\alpha \to \infty — разлома нет.
  • Отражённое поле — то же отражение, что у непроницаемого разлома, минус «утечка»: за отражением стоит линия воображаемых нагнетательных скважин с общим расходом qq, рассеянных на длину порядка L/(2α)L/(2\alpha).
  • Производная поднимается, как у непроницаемого разлома, достигает максимума ниже 2me2m_e (тем ниже, чем больше α\alpha) и возвращается к mem_e за время порядка (L/α)2/η(L/\alpha)^2/\eta. Поздняя полулогарифмическая прямая параллельна первой и сдвинута на 2meSf2m_eS_f, Sf=e4αE1(4α)S_f = \e^{4\alpha}E_1(4\alpha).
  • По высоте максимума находят α\alpha, по началу подъёма — расстояние LL, затем kf/lf=αk/Lk_f/l_f = \alpha k/L. Ширину и проницаемость зоны по отдельности исследование не даёт.
  • Окно разрешения ограничено: при α\alpha меньше примерно 0,01 разлом в исследовании обычной длительности почти неотличим от непроницаемого, при α\alpha больше 2–3 его горб тонет в шуме производной.
Ключевые соотношения
α=kf/lfk/L=kfLk lf\alpha = \dfrac{k_f/l_f}{k/L} = \dfrac{k_fL}{k\,l_f}

безразмерная проводимость разлома (одинаковая толщина по обе стороны)

Δp′me=1+e−1/tDL[1−2π αtDL  erfcx(2αtDL+1tDL)]\dfrac{\Delta p'}{m_e} = 1 + \e^{-1/t_{DL}}\Big[1 - 2\sqrt{\pi}\,\alpha\sqrt{t_{DL}}\;\mathrm{erfcx}\Big(2\alpha\sqrt{t_{DL}} + \dfrac{1}{\sqrt{t_{DL}}}\Big)\Big]

производная после окончания влияния ствола; tDL=ηt/L2t_{DL} = \eta t/L^2

Sf=e4αE1(4α)S_f = \e^{4\alpha}E_1(4\alpha)

кажущийся скин разлома: сдвиг поздней прямой 2meSf2m_eS_f

kflf=α kL\dfrac{k_f}{l_f} = \dfrac{\alpha\,k}{L}

проводимость зоны разлома по результату исследования, мД/м

F=α Δxα Δx+LF = \dfrac{\alpha\,\Delta x}{\alpha\,\Delta x + L}

множитель проводимости разлома для сетки с ячейками Δx\Delta x поперёк разлома

Где чаще всего ошибаются
  • Принимать максимум производной ниже 2me2m_e за признак того, что разлом далеко: расстояние даёт начало подъёма, а высоту максимума — утечка.
  • Считать, что возврат производной к полке означает отсутствие разлома: поздняя прямая сдвинута на кажущийся скин разлома.
  • Приписывать исследованию проницаемость зоны разлома: оно даёт только kf/lfk_f/l_f.
  • Переносить множитель проводимости разлома в модель с другой сеткой без пересчёта.
  • Принимать спад производной в конце кривой восстановления давления (КВД) за утечку через разлом, не проверив влияние времени работы скважины перед остановкой.

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

  • описать модель тонкой полупроницаемой зоны разлома и объяснить смысл безразмерной проводимости α\alpha как отношения двух сопротивлений;
  • объяснить отклик давления и производной — подъём, максимум ниже 2me2m_e, возврат к полке и сдвиг поздней прямой на кажущийся скин разлома — через отражение с «утечкой»;
  • оценить α\alpha по высоте максимума производной, расстояние до разлома — по началу подъёма, и пересчитать их в проводимость зоны разлома kf/lfk_f/l_f, переток через разлом и множитель проводимости модели;
  • определить, в каком диапазоне α\alpha исследование заданной длительности отличает частично проводящий разлом от непроницаемого и от отсутствия разлома;
  • назвать модели, дающие похожую картину, и данные, которые помогают их различить.

Разлом, который пропускает флюид

Непроницаемым разлом бывает, когда напротив коллектора оказывается непроницаемая порода или когда зона разлома запечатана глинистой примазкой, раздроблёнными зёрнами или цементом (разлом как граница). Но часто смещение меньше толщины пласта: проницаемые пропластки по разные стороны остаются в контакте, а зона разлома лишь ухудшена — зёрна частично раздроблены, глина размазана тонкой плёнкой, поры частично залечены. Тогда флюид проходит через разлом, но с добавочным сопротивлением. Геология разломов — в уроке Разломы.

Для разработки важна именно эта промежуточная проводимость. Через частично проводящий разлом нагнетание в одном блоке со временем поддерживает давление в соседнем, но с запаздыванием и потерями; блок за разломом «отдаёт» нефть скважине медленнее, чем без разлома. В гидродинамической модели такой разлом задают множителем проводимости между ячейками от 0 до 1 (множитель проводимости), и исследование скважины — один из немногих способов оценить его по измерениям, а не по аналогии.

Простейшая модель — полупроницаемая перегородка (semipermeable barrier): вертикальная зона шириной lfl_f и проницаемостью kfk_f через всю толщину пласта; по обе стороны — один и тот же однородный пласт. Зона тонкая, и её собственной упругоёмкостью пренебрегают: сколько флюида вошло в неё с одной стороны, столько вышло с другой, а на ней самой давление скачком меняется. Поток через единицу её площади подчиняется закону Дарси для слоя толщиной lfl_f:

Поток через зону разлома
u=kfμ⋅p1−p2lfu = \frac{k_f}{\mu}\cdot\frac{p_1 - p_2}{l_f}
где:
  • uuскорость фильтрации через разлом — расход в пластовых условиях на единицу площади зоны, м/с
  • kf, lfk_f,\ l_fпроницаемость и ширина зоны разлома (поперёк разлома), м², м
  • p1, p2p_1,\ p_2давление непосредственно по обе стороны зоны, Па
  • μ\muвязкость флюида, Па·с

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

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

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

В других системах единиц (1)
Метрические промысловые единицы
qr=0,0853 kflf⋅A (p1−p2)μq_r = 0{,}0853\,\frac{k_f}{l_f}\cdot\frac{A\,(p_1 - p_2)}{\mu}

qrq_r — м³/сут в пластовых условиях, kfk_f — мД, lfl_f — м, AA — площадь зоны (длина участка разлома на толщину пласта), м²; давления — МПа; μ\mu — мПа·с.

0,08530{,}0853 — постоянная закона Дарси в этих единицах (Проницаемость горных пород).

Безразмерная проводимость разлома

Насколько разлом заметен для скважины, зависит не от его проводимости самой по себе, а от сравнения с проводимостью пласта. Сопротивление зоны разлома на единицу площади пропорционально lf/kfl_f/k_f, сопротивление пласта на пути от скважины до разлома — L/kL/k. Их отношение и есть главный параметр задачи:

Безразмерная проводимость разлома
α=kf/lfk/L=kf Lk lf\alpha = \frac{k_f/l_f}{k/L} = \frac{k_f\,L}{k\,l_f}
где:
  • α\alphaбезразмерная проводимость разлома, —
  • kf, lfk_f,\ l_fпроницаемость и ширина зоны разлома, мД, м
  • kkпроницаемость пласта, мД
  • LLрасстояние от скважины до разлома, м

Единицы: безразмерная; kfk_f и kk, lfl_f и LL — в одинаковых единицах

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

α=0\alpha = 0 — непроницаемый разлом, α→∞\alpha \to \infty — разлома для потока нет. Параметр растёт с расстоянием: один и тот же разлом для близкой скважины выглядит более непроницаемым, а для далёкой — более «прозрачным», потому что сопротивление пласта до него растёт, а сопротивление зоны — нет.

Как выглядит разлом при разных α\alpha
α\alpha
меньше 0,01
Максимум Δp′/me\Delta p'/m_e
больше 1,84
Что видно в исследовании
почти как непроницаемый; возврат к полке — через десятки тысяч L2/ηL^2/\eta
α\alpha
0,01–0,3
Максимум Δp′/me\Delta p'/m_e
1,84–1,32
Что видно в исследовании
явный горб и возврат к полке; α\alpha определяется лучше всего
α\alpha
0,3–2
Максимум Δp′/me\Delta p'/m_e
1,32–1,08
Что видно в исследовании
небольшой горб, заметный при хорошем манометре
α\alpha
больше 2–3
Максимум Δp′/me\Delta p'/m_e
меньше 1,08–1,06
Что видно в исследовании
разлом теряется в шуме производной

Значения максимума — по таблице производной в разделе ниже.

Расчётный примерДве гипотезы о зоне разлома
Исходные данные (условные)
Пласт
kk = 80 мД, ϕ\phi = 0,18, μ\mu = 1,2 мПа·с, ct=1,5⋅10−3c_t = 1{,}5\cdot10^{-3} МПа⁻¹; разлом по сейсмике в 200 м от скважины
Гипотеза 1: зона дробления
ширина lfl_f = 3 м, kfk_f = 0,6 мД
Гипотеза 2: глинистая примазка
ширина lfl_f = 0,05 м, kfk_f = 0,002 мД
Решение
  1. Гипотеза 1: α=0,6⋅200/(80⋅3)=0,5\alpha = 0{,}6\cdot200/(80\cdot3) = 0{,}5. Гипотеза 2: α=0,002⋅200/(80⋅0,05)=0,1\alpha = 0{,}002\cdot200/(80\cdot0{,}05) = 0{,}1.
  2. Проводимость зоны kf/lfk_f/l_f: 0,20 и 0,04 мД/м — вторая зона в 5 раз хуже пропускает флюид, хотя она в 60 раз тоньше.
  3. Масштаб времени: η=3,553⋅10−3⋅80/(0,18⋅1,2⋅1,5⋅10−3)=877\eta = 3{,}553\cdot10^{-3}\cdot80/(0{,}18\cdot1{,}2\cdot1{,}5\cdot10^{-3}) = 877 м²/ч, L2/η=45,6L^2/\eta = 45{,}6 ч.
  4. По таблице раздела о производной: при α\alpha = 0,5 максимум 1,238 mem_e при tDLt_{DL} = 1,85, то есть через 84 ч; при α\alpha = 0,1 — 1,527 mem_e при tDLt_{DL} = 3,94, через 180 ч.
Ответ: Зона дробления дала бы небольшой горб (1,238 mem_e примерно через 3,5 сут), глинистая примазка — заметный (1,527 mem_e примерно через 7,5 сут). Исследование длительностью в неделю различает эти гипотезы: к четвёртым суткам первая уже прошла максимум, а вторая поднялась выше 1,4 mem_e и продолжает расти.

Что это значит. Тонкая плёнка с ничтожной проницаемостью может сопротивляться сильнее, чем широкая зона дробления: значение имеет только отношение kf/lfk_f/l_f. Поэтому оценки «по керну из зоны разлома» без знания её ширины и сплошности мало что говорят о проводимости — исследование измеряет её напрямую, хотя и в виде одного числа. Данные примера условные.

Отражение с утечкой

У непроницаемого разлома давление со стороны скважины описывает одно отражение с тем же дебитом (Непроницаемый разлом у скважины: метод изображений). Частично проводящий разлом отражает возмущение не полностью: часть его проходит в соседний блок. Строгое решение (вывод — в дополнительном разделе) оказывается наглядным. Отражённое поле — то же отражение, что у непроницаемого разлома, минус «утечка»: за отражением, на продолжении линии «скважина — отражение», стоит линия воображаемых нагнетательных скважин. Их общий расход равен qq, а плотность убывает с удалением от отражения по экспоненте с длиной L/(2α)L/(2\alpha).

Частично проводящий разлом: зона разлома, отражение и утечкаскважиназона разлома: kf, lfпласт kпласт kLравноценная картина для стороны скважины:отражениетот же дебит+++утечка: нагнетательные, всего q2Lдлина около L/(2α)

Схема из двух частей. Вверху: слева скважина, правее вертикальная заштрихованная полоса — зона разлома с подписями k f и l f; по обе стороны пласт с проницаемостью k; стрелки потока идут к скважине, часть из них пересекает полосу справа налево. Размерная линия L от скважины до разлома. Внизу: та же скважина слева, штриховая линия разлома, справа на расстоянии 2L пустой кружок отражения с подписью «тот же дебит», за ним вправо ряд кружков со знаком плюс, уменьшающихся по размеру, с подписью «утечка: нагнетательные, всего q».

СхемаЧастично проводящий разлом (схема, не в масштабе). Вверху — пласт в плане: скважина, тонкая зона разлома шириной lfl_f и проницаемостью kfk_f, переток через неё. Внизу — равноценная картина для стороны скважины: отражение с тем же дебитом, что у непроницаемого разлома, и за ним линия воображаемых нагнетательных скважин («утечка») с общим расходом qq; кружки тем мельче, чем меньше доля расхода.

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

Скважина у частично проводящего разлома: безразмерное давление
pD=12(ln⁡tD+0,809)+S+12E1 ⁣(1tDL)⏟отражение−α∫0∞e−2αξ E1 ⁣((2+ξ)24tDL) dξ⏟утечкаp_D = \tfrac12\left(\ln t_D + 0{,}809\right) + S + \underbrace{\tfrac12 E_1\!\left(\frac{1}{t_{DL}}\right)}_{\text{отражение}} - \underbrace{\alpha\int_0^\infty \e^{-2\alpha\xi}\,E_1\!\left(\frac{(2 + \xi)^2}{4t_{DL}}\right)\dd\xi}_{\text{утечка}}
где:
  • pDp_Dбезразмерное снижение забойного давления 2πkh Δp/(qBμ)2\pi kh\,\Delta p/(qB\mu), со скином, —
  • tD, tDLt_D,\ t_{DL}безразмерное время kt/(ϕμctrw2)kt/(\phi\mu c_tr_w^2) и оно же, отнесённое к L2L^2: tDL=kt/(ϕμctL2)=ηt/L2t_{DL} = kt/(\phi\mu c_tL^2) = \eta t/L^2, —
  • ξ\xiрасстояние за отражением в единицах LL (переменная интегрирования), —
  • α\alphaбезразмерная проводимость разлома, —
  • E1E_1интегральная показательная функция (функция E₁), —

Единицы: безразмерная; логарифм натуральный

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

При α=0\alpha = 0 интеграл исчезает, и остаётся формула непроницаемого разлома. При больших α\alpha утечка почти полностью гасит отражение. Последнее слагаемое можно записать и как интеграл по времени с дополнительной функцией ошибок — это та же величина в другом виде (дополнительный раздел).

Стадии отклика скважины у частично проводящего разлома
tDL=ηt/L2t_{DL} = \eta t/L^2
меньше 0,2
Что происходит
возмущение не дошло до разлома
Производная
полка mem_e
tDL=ηt/L2t_{DL} = \eta t/L^2
0,2 — время максимума
Что происходит
действует отражение; утечка ещё мала
Производная
подъём почти как у непроницаемого
tDL=ηt/L2t_{DL} = \eta t/L^2
около времени максимума
Что происходит
утечка набирает силу
Производная
максимум ниже 2me2m_e
tDL=ηt/L2t_{DL} = \eta t/L^2
много больше 1/α21/\alpha^2
Что происходит
отражение и утечка почти гасят друг друга
Производная
возврат к mem_e; давление сдвинуто на 2meSf2m_eS_f

Время максимума — от 1,2 до 16 единиц L2/ηL^2/\eta при α\alpha от 2 до 0,01 (таблица в следующем разделе) и около 70 при α\alpha = 0,001.

Давление и производная

Производная вклада отражения по ln⁡t\ln t — знакомая ступень 12e−1/tDL\tfrac12\e^{-1/t_{DL}}, производная утечки выражается через масштабированную дополнительную функцию ошибок erfcx(x)=ex2erfc(x)\mathrm{erfcx}(x) = \e^{x^2}\mathrm{erfc}(x). Вместе:

Производная давления у частично проводящего разлома
Δp′me=1+e−1/tDL[1−2π αtDL  erfcx(2αtDL+1tDL)]\frac{\Delta p'}{m_e} = 1 + \e^{-1/t_{DL}}\left[1 - 2\sqrt{\pi}\,\alpha\sqrt{t_{DL}}\;\mathrm{erfcx}\left(2\alpha\sqrt{t_{DL}} + \frac{1}{\sqrt{t_{DL}}}\right)\right]
где:
  • Δp′\Delta p'логарифмическая производная снижения давления, Па
  • mem_eуровень полки радиального притока, qBμ/(4πkh)qB\mu/(4\pi kh), Па
  • tDLt_{DL}безразмерное время ηt/L2\eta t/L^2, —
  • erfcx\mathrm{erfcx}масштабированная дополнительная функция ошибок, ex2erfc(x)\e^{x^2}\mathrm{erfc}(x); при больших xx близка к 1/(xπ)1/(x\sqrt\pi), —

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

Справедливо, если
  • те же, что для давления; влияние ствола закончилось.

При α=0\alpha = 0 квадратная скобка равна 1, и получается ступень непроницаемого разлома 1+e−1/tDL1 + \e^{-1/t_{DL}}. При любом α>0\alpha > 0 скобка сначала близка к 1, потом убывает к нулю: производная проходит максимум и возвращается к полке.

В других системах единиц (1)
Метрические промысловые единицы
tDL=3,553⋅10−3 k tϕμctL2,me=0,933 qBμkht_{DL} = \frac{3{,}553\cdot10^{-3}\,k\,t}{\phi\mu c_tL^2},\qquad m_e = 0{,}933\,\frac{qB\mu}{kh}

kk — мД, tt — ч, μ\mu — мПа·с, ctc_t — 1/МПа, LL и hh — м, qq — м³/сут; mem_e — МПа.

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

Производная скважины у частично проводящего разлома при разных безразмерных проводимостях0,1110100100010 0001,001,251,501,752,00ηt / L²Δp′ / mₑα = 0 (непроницаемый)сверху вниз: α = 0,01; 0,03; 0,1; 0,3; 1

График: по горизонтали безразмерное время eta t делённое на L в квадрате от 0,1 до 10 000 на логарифмической шкале, по вертикали отношение производной к полке от 0,8 до 2,1. Все кривые до 0,2 лежат на уровне 1 и вначале поднимаются вместе. Кривая непроницаемого разлома выходит на 2 и остаётся на нём. Кривая альфа 0,01 доходит до 1,84 около 16 и к 10 000 опускается до 1,1. Кривая альфа 0,03 — до 1,72 около 8, затем опускается. Кривая альфа 0,1 — до 1,53 около 4 и к 1000 почти возвращается к 1. Кривая альфа 0,3 — до 1,32 около 2. Кривая альфа 1 — небольшой горб 1,14 около 1,5.

РасчётРасчётная производная скважины у частично проводящего разлома, отнесённая к полке, без влияния ствола, при разных α\alpha; шкала времени логарифмическая. Верхняя кривая — непроницаемый разлом (α=0\alpha = 0).
Максимум производной и возврат к полке
α\alpha
0,01
Максимум Δp′/me\Delta p'/m_e
1,844
tDLt_{DL} максимума
15,57
tDLt_{DL} возврата до 1,1
9797
SfS_f
2,79
α\alpha
0,02
Максимум Δp′/me\Delta p'/m_e
1,774
tDLt_{DL} максимума
10,10
tDLt_{DL} возврата до 1,1
2568
SfS_f
2,20
α\alpha
0,05
Максимум Δp′/me\Delta p'/m_e
1,646
tDLt_{DL} максимума
5,83
tDLt_{DL} возврата до 1,1
467
SfS_f
1,49
α\alpha
0,1
Максимум Δp′/me\Delta p'/m_e
1,527
tDLt_{DL} максимума
3,94
tDLt_{DL} возврата до 1,1
139
SfS_f
1,05
α\alpha
0,2
Максимум Δp′/me\Delta p'/m_e
1,397
tDLt_{DL} максимума
2,76
tDLt_{DL} возврата до 1,1
45,5
SfS_f
0,69
α\alpha
0,5
Максимум Δp′/me\Delta p'/m_e
1,238
tDLt_{DL} максимума
1,85
tDLt_{DL} возврата до 1,1
11,7
SfS_f
0,36
α\alpha
1
Максимум Δp′/me\Delta p'/m_e
1,144
tDLt_{DL} максимума
1,46
tDLt_{DL} возврата до 1,1
4,2
SfS_f
0,21
α\alpha
2
Максимум Δp′/me\Delta p'/m_e
1,081
tDLt_{DL} максимума
1,24
tDLt_{DL} возврата до 1,1
—
SfS_f
0,11

Расчёт по формуле производной; «—» — максимум ниже 1,1. Время в единицах L2/ηL^2/\eta; SfS_f — кажущийся скин разлома (формула ниже).

Поздняя часть записи тоже несёт сведения о разломе. Когда tDLt_{DL} много больше 1/α21/\alpha^2, производная приближается к полке медленно — обратно пропорционально времени, а давление выходит на прямую, параллельную первой, но сдвинутую вверх: разлом действует как добавочный скин.

Поздний участок: остаток производной и кажущийся скин разлома
Δp′me−1≈1+4α8α2 tDL,Δp=me[ln⁡tD+0,809+2S+2Sf],Sf=e4αE1(4α)\frac{\Delta p'}{m_e} - 1 \approx \frac{1 + 4\alpha}{8\alpha^2\,t_{DL}},\qquad \Delta p = m_e\left[\ln t_D + 0{,}809 + 2S + 2S_f\right],\quad S_f = \e^{4\alpha}E_1(4\alpha)
где:
  • SfS_fкажущийся скин разлома: добавочный перепад на поздней прямой в единицах 2me2m_e, —
  • E1E_1интегральная показательная функция, —
  • tD, tDL, St_D,\ t_{DL},\ Sбезразмерное время по радиусу скважины и по расстоянию до разлома, скин скважины, —

Единицы: остаток и SfS_f безразмерны; Δp\Delta p и mem_e — в одних единицах давления; логарифм натуральный

Справедливо, если
  • tDL≫1/α2t_{DL} \gg 1/\alpha^2 и tDL≫1t_{DL} \gg 1: остаток производной мал; формула остатка приближённая (завышает его на несколько процентов);
  • других границ возмущение не достигло.

SfS_f убывает с α\alpha: около 2,8 при α\alpha = 0,01, около 1 при 0,1, 0,2 при 1. Если в записи видны и первая, и поздняя параллельная прямая, расстояние между ними 2meSf2m_eS_f даёт α\alpha без чтения максимума. Но поздней прямой надо ждать, пока производная почти вернётся к полке: десятки единиц L2/ηL^2/\eta при α\alpha около 0,3 и сотни — при 0,1 и меньше. Это удаётся только для довольно проводящих разломов.

Вклад частично проводящего разлома в безразмерное давление скважины и его поздние уровни0,1110100100010 0000,00,51,01,52,02,53,0ηt / L²Вклад разлома в pDнепроницаемыйα = 0,1: 1,05α = 0,3: 0,53α = 1: 0,21

График: по горизонтали eta t делённое на L в квадрате от 0,1 до 100 000 на логарифмической шкале, по вертикали вклад разлома в безразмерное давление от 0 до 3. Кривая непроницаемого разлома начинает расти около 0,3 и дальше поднимается прямой, доходя до 3 около 700. Кривая альфа 0,1 растёт медленнее и выходит на штриховой уровень 1,05; кривая альфа 0,3 — на 0,53; кривая альфа 1 — на 0,21.

РасчётРасчётный вклад разлома в безразмерное давление скважины pDp_D при разных α\alpha; шкала времени логарифмическая. У непроницаемого разлома вклад растёт без конца (наклон удваивается), у частично проводящего выходит на постоянную SfS_f (штриховые линии).

Проводимость разлома по исследованию

Порядок интерпретации
  1. 1
    Радиальный приток

    По полке производной — khkh и kk, по отношению давления к полке — скин (Диагностика по производной давления).

  2. 2
    Расстояние

    По самому началу подъёма (производная выше полки на 10–20 %) — расстояние LL по формуле одного разлома (ступень производной). Утечка сказывается и там, и LL получается завышенным: по точке 10 % — примерно на 2 % при α\alpha = 0,1, на 5 % при 0,3 и на 9 % при 0,5; более поздние точки завышают сильнее. Окончательное LL уточняют по времени максимума и подбором модели.

  3. 3
    Проводимость

    По высоте максимума производной — α\alpha из таблицы предыдущего раздела; по времени максимума — проверка LL.

  4. 4
    Проводимость зоны

    kf/lf=αk/Lk_f/l_f = \alpha k/L; при оценке ширины зоны по керну или аналогам — kfk_f.

  5. 5
    Подбор модели

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

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

Снижение давления у разлома (условные данные, выборка)
tt
36 с
Δp\Delta p, МПа
0,5245
Δp′\Delta p', МПа
0,07730
tt
6 мин
Δp\Delta p, МПа
0,6711
Δp′\Delta p', МПа
0,05841
tt
30 мин
Δp\Delta p, МПа
0,7638
Δp′\Delta p', МПа
0,05700
tt
1,00 ч
Δp\Delta p, МПа
0,8031
Δp′\Delta p', МПа
0,05678
tt
3,16 ч
Δp\Delta p, МПа
0,8686
Δp′\Delta p', МПа
0,05776
tt
6,31 ч
Δp\Delta p, МПа
0,9096
Δp′\Delta p', МПа
0,06286
tt
10,0 ч
Δp\Delta p, МПа
0,9398
Δp′\Delta p', МПа
0,06894
tt
20,0 ч
Δp\Delta p, МПа
0,9908
Δp′\Delta p', МПа
0,07774
tt
31,6 ч
Δp\Delta p, МПа
1,0277
Δp′\Delta p', МПа
0,08110
tt
50,1 ч
Δp\Delta p, МПа
1,0655
Δp′\Delta p', МПа
0,08176
tt
100 ч
Δp\Delta p, МПа
1,1213
Δp′\Delta p', МПа
0,07872
tt
200 ч
Δp\Delta p, МПа
1,1738
Δp′\Delta p', МПа
0,07329
tt
316 ч
Δp\Delta p, МПа
1,2067
Δp′\Delta p', МПа
—

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

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

График в двойном логарифмическом масштабе: время от 0,001 до 1000 часов, давление от 0,03 до 3 МПа. Точки снижения давления круто растут в первые секунды и дальше медленно доходят до 1,21 МПа к 316 часам. Оранжевая производная после горба влияния ствола ложится на полку 0,057 МПа от 0,3 до 2 часов, затем поднимается, около 50 часов достигает максимума 0,082 МПа и к 200 часам опускается до 0,073 МПа. Серая штриховая производная непроницаемого разлома поднимается вместе с ней до 10 часов, а потом продолжает расти почти до 0,11 МПа.

РасчётДиагностический график записи: снижение давления (точки) и производная по замерам (оранжевая линия). Серая штриховая линия — производная той же скважины у непроницаемого разлома на том же расстоянии; горизонтальная — полка mem_e.
Расчётный примерБезразмерная проводимость разлома и проводимость его зоны
Исходные данные (условные)
Скважина и пласт
qq = 150 м³/сут; BB = 1,3; μ\mu = 0,7 мПа·с; hh = 15 м; ϕ\phi = 0,2; ct=1,6⋅10−3c_t = 1{,}6\cdot10^{-3} МПа⁻¹; rwr_w = 0,1 м
Полка производной
среднее от 0,32 до 1,6 ч: mem_e = 0,05695 МПа; замер через 1 ч: Δp\Delta p = 0,8031 МПа
Начало подъёма
через 6,31 ч производная 0,06286 МПа
Максимум производной
0,08176 МПа через 50,1 ч
Решение
  1. Проводимость: qBμ=150⋅1,3⋅0,7=136,5qB\mu = 150\cdot1{,}3\cdot0{,}7 = 136{,}5; kh=0,933⋅136,5/0,05695=2236kh = 0{,}933\cdot136{,}5/0{,}05695 = 2236 мД·м, k=149,1k = 149{,}1 мД. Скин: ln⁡ktϕμctrw2=18,014\ln\dfrac{kt}{\phi\mu c_tr_w^2} = 18{,}014 при 1 ч; S=12 (0,8031/0,05695−18,014+4,831)=0,46S = \tfrac12\,(0{,}8031/0{,}05695 - 18{,}014 + 4{,}831) = 0{,}46.
  2. Расстояние по началу подъёма: η=3,553⋅10−3⋅149,1/(0,2⋅0,7⋅1,6⋅10−3)=2365\eta = 3{,}553\cdot10^{-3}\cdot149{,}1/(0{,}2\cdot0{,}7\cdot1{,}6\cdot10^{-3}) = 2365 м²/ч; Δp′/me=1,104\Delta p'/m_e = 1{,}104, ln⁡[1/(1,104−1)]=2,263\ln[1/(1{,}104 - 1)] = 2{,}263; L=2365⋅6,31⋅2,263=184L = \sqrt{2365\cdot6{,}31\cdot2{,}263} = 184 м.
  3. Максимум: 0,08176/0,05695=1,4360{,}08176/0{,}05695 = 1{,}436 — заметно меньше 2. Непроницаемый разлом к этому времени дал бы почти 2me2m_e: разлом пропускает флюид.
  4. По таблице максимум 1,436 лежит между строками α\alpha = 0,1 (1,527) и α\alpha = 0,2 (1,397). Доля интервала: (1,436−1,527)/(1,397−1,527)=0,700(1{,}436 - 1{,}527)/(1{,}397 - 1{,}527) = 0{,}700; интерполируя по lg⁡α\lg\alpha: lg⁡α=−1+0,700⋅0,301=−0,789\lg\alpha = -1 + 0{,}700\cdot0{,}301 = -0{,}789, α≈0,163\alpha \approx 0{,}163.
  5. Проверка по времени максимума: тем же способом lg⁡tDL=0,5955+0,700⋅(0,4409−0,5955)=0,4873\lg t_{DL} = 0{,}5955 + 0{,}700\cdot(0{,}4409 - 0{,}5955) = 0{,}4873, tDLt_{DL} = 3,07; L=ηt/tDL=2365⋅50,1/3,07=196L = \sqrt{\eta t/t_{DL}} = \sqrt{2365\cdot50{,}1/3{,}07} = 196 м — в пределах точности чтения пологого максимума.
  6. Проводимость зоны: kf/lf=αk/L=0,163⋅149,1/184=0,132k_f/l_f = \alpha k/L = 0{,}163\cdot149{,}1/184 = 0{,}132 мД/м. Если ширина зоны около 2 м, kf≈k_f \approx 0,26 мД; если это плёнка шириной 0,2 м — около 0,026 мД.
Ответ: k≈149,1k \approx 149{,}1 мД, скин 0,46; разлом в 184 м от скважины частично проводящий: α≈0,163\alpha \approx 0{,}163, kf/lf≈0,132k_f/l_f \approx 0{,}132 мД/м.

Что это значит. Запись рассчитана по модели с kk = 150 мД, SS = 0,5, разломом в 180 м и α\alpha = 0,15 (kf/lfk_f/l_f = 0,125 мД/м). Расстояние найдено с ошибкой около 2 %, α\alpha — около 9 %. Ошибка α\alpha — от высоты максимума: отношение по замерам примерно на 1 % ниже точного максимума 1,451 (производная по замерам сглаживает пологую вершину, а полка чуть приподнята хвостом влияния ствола), а точный обратный пересчёт этого отношения по формуле даёт 0,163 — интерполяция по таблице здесь ни при чём. Чувствительность велика: 1 % в отношении — около 8 % в α\alpha. Время максимума читается грубо — вершина пологая, а замеров 10 на цикл, — поэтому расстояние по нему (196 м) хуже, чем по началу подъёма. Запись кончилась при tDL≈22t_{DL} \approx 22: возврат к полке ещё впереди, но отличие от непроницаемого разлома уже бесспорно. Данные примера условные.

От безразмерной проводимости к разлому в модели

Исследование даёт α\alpha, а вместе с kk и LL — отношение kf/lfk_f/l_f. Это всё, что можно получить: ширину зоны и её проницаемость по отдельности давление в скважине не различает. Ширину зон разломов оценивают по керну скважин, пересёкших разлом, по обнажениям и сейсмическим атрибутам; она бывает от сантиметров (плёнка глинистой примазки) до метров и десятков метров (зона дробления). Поэтому kfk_f «по исследованию» всегда условна и зависит от принятой ширины.

Для расчётов перетока между блоками ширина и не нужна. Переток через участок разлома площадью AA при скачке давления на нём определяется тем же отношением kf/lfk_f/l_f:

Проводимость участка разлома
Tf=kflf A,qr=Tfμ (p1−p2)T_f = \frac{k_f}{l_f}\,A,\qquad q_r = \frac{T_f}{\mu}\,(p_1 - p_2)
где:
  • TfT_fпроводимость участка разлома — геометрическая часть, без вязкости, как проводимость между ячейками TijT_{ij} в уроке Уравнение давления для однофазной фильтрации, м³
  • AAплощадь участка: длина по простиранию на толщину пласта, контактирующего по разлому, м²
  • qrq_rпереток через участок в пластовых условиях, м³/с
  • p1−p2p_1 - p_2скачок давления на самой зоне разлома, Па
  • μ\muвязкость флюида, Па·с

Единицы: СИ

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

Отношение Tf/μT_f/\mu — межблочный «коэффициент продуктивности» разлома: переток на единицу скачка давления. В модели разлом задают через множитель к проводимости между ячейками по обе стороны разлома.

В других системах единиц (1)
Метрические промысловые единицы
qr=0,0853 Tfμ (p1−p2),Tf=kflf Aq_r = 0{,}0853\,\frac{T_f}{\mu}\,(p_1 - p_2),\qquad T_f = \frac{k_f}{l_f}\,A

TfT_f — мД·м (kf/lfk_f/l_f — мД/м, AA — м²), μ\mu — мПа·с, давления — МПа; qrq_r — м³/сут в пластовых условиях.

0,08530{,}0853 — постоянная закона Дарси в этих единицах; 0,0853 Tf/μ0{,}0853\,T_f/\mu — переток на 1 МПа скачка, м³/(сут·МПа), как коэффициент связи aija_{ij} ячеек модели.

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

Множитель проводимости разлома
F=Δx/kΔx/k+lf/kf=α Δxα Δx+LF = \frac{\Delta x/k}{\Delta x/k + l_f/k_f} = \frac{\alpha\,\Delta x}{\alpha\,\Delta x + L}
где:
  • FFмножитель проводимости между ячейками по разные стороны разлома, от 0 до 1, —
  • Δx\Delta xразмер ячеек поперёк разлома (одинаковых), м
  • kkпроницаемость ячеек (как у пласта в исследовании), мД
  • α, L\alpha,\ Lбезразмерная проводимость разлома и расстояние, к которому она отнесена, —, м

Единицы: безразмерная; Δx\Delta x и LL — в одних единицах

Справедливо, если
  • соседние ячейки одинакового размера и проницаемости; для разных ячеек половины сопротивлений Δxi/(2ki)\Delta x_i/(2k_i) складывают отдельно;
  • ширина зоны мала по сравнению с ячейкой.

Множитель зависит от сетки: при измельчении ячеек сопротивление пласта между центрами ячеек уменьшается, а сопротивление зоны — нет, и тот же разлом требует меньшего множителя. Переносить множитель из модели в модель без пересчёта нельзя; переносить можно kf/lfk_f/l_f.

Расчётный примерМножитель проводимости и переток между блоками
Исходные данные (условные)
Результат исследования
α\alpha ≈ 0,163 на расстоянии LL = 184 м; kk = 149,1 мД; kf/lfk_f/l_f ≈ 0,132 мД/м
Модель
ячейки поперёк разлома 100 м или 50 м
Участок разлома
длина 1500 м, толщина контакта пластов 15 м; вязкость нефти 0,7 мПа·с
Решение
  1. Ячейки 100 м: F=0,163⋅100/(0,163⋅100+184)=0,081F = 0{,}163\cdot100/(0{,}163\cdot100 + 184) = 0{,}081. Ячейки 50 м: F=0,163⋅50/(0,163⋅50+184)=0,042F = 0{,}163\cdot50/(0{,}163\cdot50 + 184) = 0{,}042.
  2. Площадь участка: A=1500⋅15=22 500A = 1500\cdot15 = 22\,500 м². Проводимость: Tf=0,132⋅22 500=2970T_f = 0{,}132\cdot22\,500 = 2970 мД·м; переток на 1 МПа скачка: 0,0853⋅2970/0,7=3620{,}0853\cdot2970/0{,}7 = 362 м³/(сут·МПа).
  3. При скачке давления на разломе 0,5 МПа переток qr=362⋅0,5=181q_r = 362\cdot0{,}5 = 181 м³/сут в пластовых условиях — около 139 м³/сут нефти в стандартных условиях при BB = 1,3.
Ответ: Множитель 0,081 для ячеек 100 м и 0,042 для ячеек 50 м; разлом пропускает около 362 м³/(сут·МПа) через участок 1500 м.

Что это значит. Разлом, который в двухнедельном исследовании заметно «держит» давление, в масштабе месторождения пропускает сотни кубометров в сутки на каждый мегапаскаль скачка. Для разработки это значит: через годы давление по обе стороны выровняется, и нагнетание с одной стороны будет поддерживать другую — но с запаздыванием и потерями на скачке. Отсюда практический вывод: исследование отвечает на вопрос «насколько», а не «да или нет», и результат передают в модель как kf/lfk_f/l_f, пересчитывая множитель для каждой сетки. Данные примера условные.

Чего исследование не различит

Частично проводящий разлом виден, только если запись достаточно длинная, чтобы производная заметно отстала от непроницаемого случая, и если горб выше шума. Это задаёт окно α\alpha, которое конкретное исследование способно оценить.

Расчётный примерКакие разломы покажет исследование
Исходные данные (условные)
Пласт и разлом
из примера выше: η\eta ≈ 2365 м²/ч, LL ≈ 184 м
Длительность
13 суток (316 ч) или 3 суток (72 ч)
Критерии
к концу записи производная ниже непроницаемого случая хотя бы на 0,1 me0{,}1\,m_e; горб выше полки хотя бы на 5 %
Решение
  1. Конец записи в единицах L2/ηL^2/\eta: tDL=2365⋅316/1842=22,1t_{DL} = 2365\cdot316/184^2 = 22{,}1 и 2365⋅72/1842=5,032365\cdot72/184^2 = 5{,}03.
  2. Наименьшее α\alpha, при котором производная к концу записи отстаёт от непроницаемого случая на 0,1 me0{,}1\,m_e (расчёт по формуле производной): 0,008 для 13 суток и 0,026 для 3 суток.
  3. Наибольшее α\alpha, при котором максимум ещё выше полки на 5 % (по формуле производной): около 3,4 — не зависит от длительности, если запись дошла до максимума.
Ответ: Тринадцатисуточное исследование различает разломы с α\alpha примерно от 0,008 до 3,4, трёхсуточное — от 0,026 до 3,4.

Что это значит. Нижняя граница определяется длительностью: чем короче запись, тем больше разломов выглядят непроницаемыми, и для них можно указать только верхнюю границу α\alpha. Верхнюю границу задаёт шум производной: если различим горб в 3 %, она поднимается примерно до 6, если нужен горб в 10 % — опускается до 1,6. Сами критерии условны: в каждом исследовании их выбирают по фактическому шуму. Данные примера условные.

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

Различить эти модели по одной записи часто нельзя; помогают геология, сейсмика и данные других скважин.

  • Только отношение. Исследование даёт α\alpha, то есть kf/lfk_f/l_f при известных kk и LL; ширину и проницаемость зоны по отдельности — нет.
  • Связь с расстоянием. Ошибка LL переходит в kf/lfk_f/l_f напрямую, ошибка ϕct\phi c_t — через LL как корень. Если начало подъёма не прочитано (скрыто влиянием ствола), LL и α\alpha находят только вместе, подбором, и неоднозначно.
  • Прямые измерения сообщаемости. Гидропрослушивание через разлом (interference test) — возмущающая скважина по одну сторону, реагирующая по другую — показывает сам переток. Замеры давления по пропласткам в скважинах соседних блоков во времени (Замеры давления в разрабатываемых пластах: барьеры и точность) и история разработки в модели (Адаптация модели к истории разработки) проверяют результат в масштабе месторождения.

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

Частично проводящий разлом
разлом, через который флюид перетекает с добавочным сопротивлением; в модели — тонкая полупроницаемая зона.
Безразмерная проводимость разлома α\alpha
kfL/(klf)k_fL/(kl_f): отношение проводимости зоны разлома к проводимости пласта на пути от скважины до разлома.
Проводимость зоны разлома kf/lfk_f/l_f
единственная характеристика зоны, которую даёт исследование; мД/м.
Утечка
в методе отражений — линия воображаемых нагнетательных скважин за отражением с общим расходом qq, описывающая переток через разлом.
Кажущийся скин разлома SfS_f
e4αE1(4α)\e^{4\alpha}E_1(4\alpha): постоянный добавочный перепад на поздней прямой после возврата производной к полке, в единицах 2me2m_e.
Множитель проводимости разлома FF
отношение проводимости между ячейками модели с разломом и без него; зависит от размера ячеек.
Окно разрешения
диапазон α\alpha, в котором исследование данной длительности и качества отличает частично проводящий разлом и от непроницаемого, и от отсутствия разлома.

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

Неверно

«Производная не дошла до 2me2m_e — значит, разлом далеко, и запись просто коротка».

На самом деле

Расстояние показывает начало подъёма. Если подъём начался вовремя, а максимум пройден ниже 2me2m_e и производная пошла вниз, разлом пропускает флюид; высота максимума даёт α\alpha.

Неверно

«Производная вернулась к полке — разлома нет».

На самом деле

Поздняя прямая параллельна первой, но сдвинута на 2meSf2m_eS_f: разлом остаётся добавочным сопротивлением, как скин. Отсутствие разлома означало бы отсутствие горба вообще.

Неверно

«По исследованию определили проницаемость зоны разлома».

На самом деле

Определили kf/lfk_f/l_f. Проницаемость зоны получают, только приняв её ширину, — и она меняется вместе с этим допущением.

Неверно

«Множитель проводимости разлома — его свойство, его можно взять из другой модели».

На самом деле

Множитель зависит от размеров ячеек: F=αΔx/(αΔx+L)F = \alpha\Delta x/(\alpha\Delta x + L). Переносить между моделями нужно kf/lfk_f/l_f, а множитель пересчитывать для каждой сетки.

Неверно

«Спад производной в конце КВД — признак утечки через разлом».

На самом деле

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

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

Выводы: отражение с утечкой, производная и поздние асимптоты

Коэффициент отражения. Разложим возмущение давления вдоль разлома на гармоники eiωy\e^{\mathrm i\omega y}; в изображениях Лапласа каждая из них затухает поперёк разлома как e−β∣x∣\e^{-\beta|x|}, β=s+ω2\beta = \sqrt{s + \omega^2} (длины — в единицах LL). Со стороны скважины к падающей волне добавляется отражённая с коэффициентом RR, за разломом идёт прошедшая с коэффициентом TT. Непрерывность потока через разлом даёт T=1−RT = 1 - R, условие на зоне — поток пропорционален скачку давления, ∂p/∂x=α (p+−p−)\partial p/\partial x = \alpha\,(p_+ - p_-), — даёт β(1−R)=2αR\beta(1 - R) = 2\alpha R, то есть R=β/(β+2α)R = \beta/(\beta + 2\alpha). При α=0\alpha = 0 получаем R=1R = 1 — полное отражение (непроницаемый разлом), при α→∞\alpha \to \infty — R=0R = 0.

Утечка как линия скважин. Запишем R=1−2α/(β+2α)R = 1 - 2\alpha/(\beta + 2\alpha) и представим 1/(β+2α)1/(\beta + 2\alpha) интегралом ∫0∞e−(β+2α)ξ dξ\int_0^\infty \e^{-(\beta + 2\alpha)\xi}\dd\xi. Множитель e−βξ\e^{-\beta\xi} — это сдвиг источника ещё на ξ\xi дальше. Значит, отражённое поле — отражение минус источники, распределённые за ним с плотностью 2αe−2αξ2\alpha\e^{-2\alpha\xi}, общий расход которых равен qq. В скважине вклад каждого источника на расстоянии 2+ξ2 + \xi равен 12E1((2+ξ)2/(4tDL))\tfrac12E_1\big((2 + \xi)^2/(4t_{DL})\big) — отсюда формула давления урока. В изображениях Лапласа то же самое: spˉ=∫0∞e−2B/(B+2α)  dus\bar p = \int_0^\infty \e^{-2B}/(B + 2\alpha)\,\dd u, B=s+u2B = \sqrt{s + u^2}; при α=0\alpha = 0 это K0(2s)K_0(2\sqrt s) — изображение одного отражения.

Производная. Продифференцируем вклад утечки по ln⁡t\ln t, пользуясь тем, что t  dE1(a/t)/ dt=e−a/tt\,\dd E_1(a/t)/\dd t = \e^{-a/t}. Подстановка w=(2+ξ)/(2tDL)w = (2 + \xi)/(2\sqrt{t_{DL}}) и выделение полного квадрата в показателе дают выражение с дополнительной функцией ошибок erfc\mathrm{erfc}; переход к erfcx\mathrm{erfcx} превращает его в формулу производной урока. Проинтегрировав производную по времени, получаем вклад утечки в давление в том виде, в каком точное решение обычно и записывают:

Вклад утечки в давление: форма с дополнительной функцией ошибок
α∫0∞e−2αξE1 ⁣((2+ξ)24tDL) dξ=π α∫0tDLe4α2u+4α erfc(2αu+1u) duu\alpha\int_0^\infty \e^{-2\alpha\xi}E_1\!\left(\frac{(2+\xi)^2}{4t_{DL}}\right)\dd\xi = \sqrt{\pi}\,\alpha\int_0^{t_{DL}} \e^{4\alpha^2u + 4\alpha}\,\mathrm{erfc}\left(2\alpha\sqrt u + \frac{1}{\sqrt u}\right)\frac{\dd u}{\sqrt u}
где:
  • uuбезразмерное время tDLt_{DL} под знаком интеграла, —
  • erfc\mathrm{erfc}дополнительная функция ошибок, —
  • α, ξ, tDL\alpha,\ \xi,\ t_{DL}как в формуле давления урока, —

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

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

Левая часть — утечка как линия источников, правая — та же величина как интеграл по времени. Подынтегральное выражение справа — скорость роста вклада утечки во времени.

Поздние асимптоты. При больших аргументах erfcx(x)≈1xπ(1−12x2)\mathrm{erfcx}(x) \approx \tfrac{1}{x\sqrt\pi}\big(1 - \tfrac{1}{2x^2}\big); подстановка в квадратную скобку даёт остаток производной (1+4α)/(8α2tDL)(1 + 4\alpha)/(8\alpha^2t_{DL}). Предел вклада разлома в давление — по теореме о конечном значении: lim⁡spˉ\lim s\bar p при s→0s \to 0, то есть интеграл ∫0∞e−2u du/(u+2α)\int_0^\infty \e^{-2u}\dd u/(u + 2\alpha), равный e4αE1(4α)\e^{4\alpha}E_1(4\alpha). Это и есть SfS_f. Все три формы решения — интеграл по утечке, интеграл с erfc и обратное преобразование Лапласа — совпадают численно до шестого знака и лучше.

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

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

  1. 1

    Полка производной — 0,050 МПа, максимум после неё — 0,076 МПа, затем производная пошла вниз. Оцените безразмерную проводимость разлома по таблице урока.

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

    Отношение 0,076/0,050=1,520{,}076/0{,}050 = 1{,}52 — почти строка α\alpha = 0,1 (1,527). Значит, α≈0,1\alpha \approx 0{,}1: интерполяция между строками 0,1 и 0,2 даёт около 0,10. Расстояние берут по началу подъёма, после чего kf/lf=αk/Lk_f/l_f = \alpha k/L.

  2. 2

    Зона разлома шириной 0,5 м с проницаемостью 0,05 мД, пласт 40 мД, скважина в 300 м от разлома. Найдите α\alpha и оцените, заметит ли исследование этот разлом.

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

    α=0,05⋅300/(40⋅0,5)=0,75\alpha = 0{,}05\cdot300/(40\cdot0{,}5) = 0{,}75. По таблице максимум производной между 1,24 (α\alpha = 0,5) и 1,14 (α\alpha = 1), то есть около 1,18 полки. Горб в 15–20 % заметен при хорошем манометре и чистой записи, но α\alpha по нему будет неточной: в этой области максимум меняется медленно.

  3. 3

    Один и тот же разлом: скважина A стоит в 100 м от него, скважина B — в 400 м. У какой исследование покажет более высокий горб производной?

    1. А.У A
    2. Б.У B
    3. В.Одинаковый: разлом один
    4. Г.Горба не будет ни у одной
    Показать ответ

    Правильный ответ: А. У A

    α=kfL/(klf)\alpha = k_fL/(kl_f) пропорционально расстоянию: у A она вчетверо меньше, и разлом для неё ближе к непроницаемому — горб выше. У B подъём начинается в 16 раз позже (L2L^2 больше в 16 раз), а горб ниже. Сопротивление пласта на пути до разлома растёт с расстоянием, а сопротивление зоны — нет.

  4. 4

    По исследованию α\alpha = 0,2 на расстоянии 250 м. Какой множитель проводимости задать для ячеек 200 м поперёк разлома, а какой для ячеек 50 м?

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

    F=αΔx/(αΔx+L)F = \alpha\Delta x/(\alpha\Delta x + L). Для 200 м: 40/(40+250)≈0,1440/(40 + 250) \approx 0{,}14. Для 50 м: 10/(10+250)≈0,03810/(10 + 250) \approx 0{,}038. Мелкие ячейки требуют меньшего множителя: сопротивление пласта между их центрами меньше, а сопротивление зоны разлома прежнее.

  5. 5

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

    1. А.Да, если подъём был больше 50 %
    2. Б.Да, спад после подъёма — признак утечки
    3. В.Нет: сначала надо исключить влияние короткой работы перед остановкой
    4. Г.Нет: спад производной означает канал
    Показать ответ

    Правильный ответ: В. Нет: сначала надо исключить влияние короткой работы перед остановкой

    После короткой работы производная КВД опускается в конце сама по себе. Её строят по эквивалентному времени или функции суперпозиции (Влияние времени работы на анализ КВД); если спад исчезает, разлом может быть и непроницаемым. Утечку обсуждают, только если спад остаётся.

  6. 6

    Почему исследование скважины даёт только отношение kf/lfk_f/l_f, а не проницаемость зоны разлома?

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

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

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