Skip to content
Brent$117,6Газ TTF$1 192USD453,87 ₸
MunaiHub

Кітапхана әзірге орыс тілінде. Ағылшын және қазақ тілдеріндегі нұсқалары кейінірек қосылады.

ГДИСВлияние ствола и эталонные кривые Тереңдетілген 22 мин оқу

Влияние времени работы на анализ КВД

Влияние времени работы (producing-time effect) — искажение кривой восстановления давления (КВД), которое появляется потому, что перед остановкой скважина работала ограниченное время tpt_p. После остановки воронка депрессии, созданная за время работы, продолжает расходиться по пласту, поэтому рост давления через время Δt\Delta t после остановки меньше, чем снижение давления через то же время после пуска, а наклон КВД со временем уменьшается — тем сильнее, чем короче была работа. Эквивалентное время Δte=tpΔt/(tp+Δt)\Delta t_e = t_p\Delta t/(t_p + \Delta t) (Agarwal equivalent time) пересчитывает время после остановки так, что КВД можно обрабатывать как кривую падения давления после пуска.

Не үшін қажет

Короткая работа перед остановкой — обычное дело: при испытании пласта на трубах основной приток длится часы, а заключительное восстановление — в полтора раза дольше; новую скважину часто останавливают на исследование вскоре после освоения. Методы, построенные для кривой падения давления (КПД), — метод MDH (Миллера — Дайеса — Хатчинсона), эталонные кривые, логарифмическая производная по ln⁡Δt\ln\Delta t — на такой КВД завышают проницаемость и рисуют ложные «границы постоянного давления». Урок показывает, как распознать влияние времени работы и как его учесть — графиком Хорнера, эквивалентным временем, производной по эквивалентному времени — и когда ни один из этих приёмов не помогает.

Қысқаша негізгісі

  • КВД — сумма двух откликов: продолжающегося снижения давления от «работающей» скважины и роста от воображаемой нагнетательной скважины. Поэтому из отклика на пуск вычитается то, на сколько ещё снизилось бы давление за время Δt\Delta t, продолжай скважина работать.
  • Рост давления после остановки не может превысить pi−pwf(tp)p_i - p_{wf}(t_p). При радиальном притоке наклон КВД по ln⁡Δt\ln\Delta t — он же производная Δp′\Delta p' — равен me tp/(tp+Δt)m_e\,t_p/(t_p + \Delta t): на 9 % меньше полки при Δt=0,1 tp\Delta t = 0{,}1\,t_p и вдвое меньше при Δt=tp\Delta t = t_p.
  • Функция Хорнера учитывает время работы точно; метод MDH верен только при Δt≪tp\Delta t \ll t_p.
  • Эквивалентное время Δte\Delta t_e превращает КВД в КПД: рост давления через Δt\Delta t после остановки равен снижению давления через Δte\Delta t_e после пуска. Оно всегда меньше tpt_p: КВД воспроизводит КПД лишь до времени, равного времени работы.
  • Производная по ln⁡Δte\ln\Delta t_e — то же, что производная по функции Хорнера, — при радиальном притоке равна полке mem_e при любом времени работы. Её откладывают и против Δte\Delta t_e, и против настоящего времени Δt\Delta t.
  • Пересчёт времени требует, чтобы до остановки приток вышел на радиальный режим; при меняющемся дебите, границах, в газовых скважинах нужны суперпозиция всей истории дебита и модель.
Негізгі қатынастар
Δte=tp Δttp+Δt\Delta t_e = \dfrac{t_p\,\Delta t}{t_p + \Delta t}

эквивалентное время; tpt_p и Δt\Delta t в одних единицах

ΔpBU(Δt)=ΔpDD(Δte)\Delta p_{BU}(\Delta t) = \Delta p_{DD}(\Delta t_e)

КВД как КПД: радиальный приток в бесконечном пласте

Δp′=me tptp+Δt\Delta p' = m_e\,\dfrac{t_p}{t_p + \Delta t}

местный наклон КВД по ln⁡Δt\ln\Delta t — то, что видят метод MDH и производная по Δt\Delta t

Δpe′=tp+Δttp Δp′=me\Delta p'_e = \dfrac{t_p + \Delta t}{t_p}\,\Delta p' = m_e

производная по эквивалентному времени на участке радиального притока

Δt=tp Ttp−T\Delta t = \dfrac{t_p\,T}{t_p - T}

сколько ждать после остановки, чтобы эквивалентное время дошло до TT

Жиі қателіктер
  • Принимать спад производной в конце КВД после короткой работы за границу постоянного давления или подпитку.
  • Обрабатывать методом MDH КВД длиннее десятой доли времени работы.
  • Совмещать КВД с эталонными кривыми КПД против Δt\Delta t, а не против Δte\Delta t_e.
  • Брать за tpt_p длительность короткого последнего притока после долгой работы: время сожмётся слишком сильно.
  • Применять эквивалентное время, когда до остановки приток не вышел на радиальный режим.

Сабақтан кейін сіз

  • объяснить через суперпозицию, почему рост давления после короткой работы отстаёт от снижения давления после пуска и не может превысить pi−pwf(tp)p_i - p_{wf}(t_p);
  • оценить, насколько метод MDH и производная по ln⁡Δt\ln\Delta t занижают наклон при заданном отношении Δt/tp\Delta t/t_p, и выбрать между MDH, методом Хорнера и эквивалентным временем;
  • рассчитать эквивалентное время и обработать КВД как КПД: найти проводимость, скин по формуле КПД и начальное давление на прямой при Δte=tp\Delta t_e = t_p;
  • вычислить по замерам производную по эквивалентному времени и найти проводимость по её полке;
  • спланировать длительность притока и остановки, чтобы КВД воспроизвела нужный участок КПД;
  • назвать условия, при которых эквивалентное время неприменимо, и что делать вместо него.

Почему КВД «помнит» время работы

Остановку скважины удобно представить по принципу суперпозиции: скважина как бы продолжает работать с дебитом qq, а с момента остановки в той же точке начинает работать воображаемая нагнетательная скважина с тем же расходом (почему давление восстанавливается). Значит, после остановки в пласте одновременно идут два процесса. У самой скважины давление растёт. А воронка депрессии, созданная за время работы, продолжает расходиться: флюид перетекает к стволу из всё более далёких частей пласта, и там давление ещё какое-то время снижается. Чем короче была работа, тем меньше воронка и тем быстрее она заполняется — восстанавливать почти нечего.

Обозначим ΔpDD(t)=pi−pwf(t)\Delta p_{DD}(t) = p_i - p_{wf}(t) — снижение забойного давления через время tt после пуска с постоянным дебитом из пласта с выровненным давлением (кривая падения давления, КПД; drawdown) и ΔpBU(Δt)\Delta p_{BU}(\Delta t) — рост давления через время Δt\Delta t после остановки (buildup), то есть pws(Δt)−pwf(Δt=0)p_{ws}(\Delta t) - p_{wf}(\Delta t = 0). Суперпозиция связывает их одной строкой:

КВД через отклик на пуск
ΔpBU(Δt)=ΔpDD(Δt)−[ΔpDD(tp+Δt)−ΔpDD(tp)]\Delta p_{BU}(\Delta t) = \Delta p_{DD}(\Delta t) - \Big[\Delta p_{DD}(t_p + \Delta t) - \Delta p_{DD}(t_p)\Big]
мұндағы:
  • ΔpBU(Δt)\Delta p_{BU}(\Delta t)рост забойного давления через время Δt\Delta t после остановки: pws(Δt)−pwf(Δt=0)p_{ws}(\Delta t) - p_{wf}(\Delta t = 0), Па
  • ΔpDD(t)\Delta p_{DD}(t)снижение забойного давления через время tt после пуска с постоянным дебитом qq из пласта с давлением pip_i: pi−pwf(t)p_i - p_{wf}(t), Па
  • tpt_pвремя работы с дебитом qq перед остановкой, с
  • Δt\Delta tвремя после остановки, с

Бірліктер: любая согласованная система единиц

Мына жағдайда орындалады
  • задача линейна: однофазный слабосжимаемый флюид, свойства пласта и флюида не зависят от давления;
  • до пуска давление в пласте выровнено, перед остановкой дебит на устье был постоянным;
  • ΔpDD\Delta p_{DD} — отклик той же скважины с тем же скином и влиянием ствола при закрытии в той же точке; тогда равенство точное, а не приближённое.

Квадратная скобка — на сколько ещё снизилось бы забойное давление за время Δt\Delta t, продолжай скважина работать. Пока Δt≪tp\Delta t \ll t_p, она мала: после долгой работы снижение давления за короткое время почти не растёт, и КВД повторяет КПД. Когда Δt\Delta t сравнимо с tpt_p, скобка сравнима с самим ΔpDD(Δt)\Delta p_{DD}(\Delta t), и восстановление заметно отстаёт от падения давления после пуска.

Второе следствие — предел. В бесконечном пласте давление после очень долгой остановки возвращается к начальному, поэтому рост давления не может превысить pi−pwf(tp)=ΔpDD(tp)p_i - p_{wf}(t_p) = \Delta p_{DD}(t_p) — снижения, накопленного за время работы. КПД растёт неограниченно, как логарифм времени, а КВД упирается в «потолок» — тем ниже, чем короче была работа. Поэтому КВД после короткой работы раньше выполаживается и её логарифмическая производная падает.

Логарифмическая производная КВД по ln Δt после разного времени работы и производная КПД0,010,11100,010,11Время после остановки Δt, чПроизводная Δp′, МПа0,5 ч2 ч6 ч24 чtp:полка meКПД

График в двойном логарифмическом масштабе: время после остановки от 0,01 до 30 часов, производная от 0,003 до 1 МПа. Все кривые начинаются около 0,75 МПа, спускаются и к 0,1 часа сливаются на полке около 0,24 МПа; штриховая кривая падения давления остаётся на полке до конца. Четыре сплошные кривые уходят вниз: кривая для получаса работы — после 0,05 часа, достигая 0,004 МПа к 30 часам; для 2 часов — после 0,2 часа; для 6 часов (оранжевая) — после 0,6 часа, до 0,04 МПа; для 24 часов — после 2 часов, до 0,1 МПа. На каждой кривой отмечена точка при времени после остановки, равном времени работы, на уровне около 0,12 МПа.

ЕсептеуРасчётная логарифмическая производная КВД по ln⁡Δt\ln\Delta t для скважины из примеров урока после работы с одним и тем же дебитом в течение 0,5; 2; 6 и 24 ч (сплошные линии; 6 ч — оранжевая) и производная КПД (штрих). В начале все кривые одинаковы: спад после влияния ствола и полка mem_e. Затем производная каждой КВД уходит вниз — тем раньше, чем короче работа. Точка на каждой кривой — момент Δt=tp\Delta t = t_p: там производная вдвое ниже полки.

Местный наклон: метод MDH и график Хорнера

Насколько быстро КВД отстаёт, видно из логарифмического приближения. Пока и пуск, и остановка проходят в режиме радиального притока в бесконечном пласте, ΔpDD(t)=me(ln⁡tD+0,809+2S)\Delta p_{DD}(t) = m_e(\ln t_D + 0{,}809 + 2S), где tDt_D — безразмерное время, а me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh) — наклон полулогарифмической прямой на единицу натурального логарифма, он же уровень полки логарифмической производной (полка производной). Подставим это в формулу предыдущего раздела и продифференцируем по ln⁡Δt\ln\Delta t:

Местный наклон КВД по логарифму времени остановки
Δp′= d ΔpBU dln⁡Δt=me tptp+Δt\Delta p' = \frac{\dd\,\Delta p_{BU}}{\dd\ln\Delta t} = m_e\,\frac{t_p}{t_p + \Delta t}
мұндағы:
  • Δp′\Delta p'логарифмическая производная КВД по ln⁡Δt\ln\Delta t — наклон графика MDH на единицу натурального логарифма, Па
  • mem_eуровень полки производной: qBμ/(4πkh)qB\mu/(4\pi kh), Па
  • tp, Δtt_p,\ \Delta tвремя работы и время после остановки, с

Бірліктер: любая согласованная система единиц

Мына жағдайда орындалады
  • радиальный приток в бесконечном пласте в моменты Δt\Delta t и tp+Δtt_p + \Delta t: влияние ствола закончилось, границы не сказываются;
  • постоянный дебит перед остановкой.

Множитель tp/(tp+Δt)t_p/(t_p + \Delta t) не зависит ни от пласта, ни от скважины — только от отношения времён. Это и есть влияние времени работы на наклон. На графике MDH (pwsp_{ws} от lg⁡Δt\lg\Delta t) точки поэтому изгибаются вниз от прямой, а производная по ln⁡Δt\ln\Delta t уходит с полки, хотя у пласта никаких границ нет.

Местный наклон КВД по ln⁡Δt\ln\Delta t в долях полки при радиальном притоке
Δt/tp\Delta t/t_p
0,01
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,990
Что это значит
метод MDH практически точен
Δt/tp\Delta t/t_p
0,05
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,952
Что это значит
наклон меньше на 5 %
Δt/tp\Delta t/t_p
0,1
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,909
Что это значит
на 9 %
Δt/tp\Delta t/t_p
0,5
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,667
Что это значит
на треть
Δt/tp\Delta t/t_p
1
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,500
Что это значит
вдвое
Δt/tp\Delta t/t_p
1,5
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,400
Что это значит
типичный конец КВД при испытании на трубах
Δt/tp\Delta t/t_p
3
tp/(tp+Δt)t_p/(t_p + \Delta t)
0,250
Что это значит
вчетверо

Это мгновенный наклон в момент Δt\Delta t. Наклон хорды, проведённой через участок КВД, — среднее этого множителя по участку (в логарифмах времени), поэтому он искажён меньше, чем наклон в правом конце участка.

Уравнение Хорнера (уравнение и график Хорнера) такого искажения не имеет: функция lg⁡[(tp+Δt)/Δt]\lg[(t_p + \Delta t)/\Delta t] содержит время работы, и при радиальном притоке точки ложатся на прямую с наклоном m=2,303 mem = 2{,}303\,m_e при любом tpt_p. Метод MDH отбрасывает слагаемое lg⁡[(tp+Δt)/tp]\lg[(t_p + \Delta t)/t_p] и потому верен только при Δt≪tp\Delta t \ll t_p: к моменту Δt\Delta t точки графика MDH лежат ниже продолжения его начальной прямой на mlg⁡[(tp+Δt)/tp]m\lg[(t_p + \Delta t)/t_p].

Эквивалентное время

Раз КВД отстаёт от КПД, потому что время после остановки «весит» меньше, можно найти такое время, при котором они совпадут. Подставим логарифмическое приближение во все три слагаемых суперпозиции. Из трёх постоянных 0,809+2S0{,}809 + 2S две взаимно сократятся, логарифмы времён сложатся в логарифм одной дроби, и останется формула КПД — только с другим временем:

Эквивалентное время
ΔpBU(Δt)=ΔpDD(Δte),Δte=tp Δttp+Δt\Delta p_{BU}(\Delta t) = \Delta p_{DD}(\Delta t_e),\qquad \Delta t_e = \frac{t_p\,\Delta t}{t_p + \Delta t}
мұндағы:
  • Δte\Delta t_eэквивалентное время: время после пуска, через которое снижение давления на КПД равнялось бы росту давления на КВД через Δt\Delta t после остановки, с
  • tp, Δtt_p,\ \Delta tвремя работы и время после остановки, с
  • ΔpBU, ΔpDD\Delta p_{BU},\ \Delta p_{DD}рост давления после остановки и снижение давления после пуска, Па

Бірліктер: любая согласованная система; Δte\Delta t_e — в тех же единицах, что tpt_p и Δt\Delta t

Мына жағдайда орындалады
  • радиальный приток в бесконечном пласте в моменты tpt_p и tp+Δtt_p + \Delta t: до остановки приток вышел на радиальный режим, а до границ возмущение не дошло;
  • в момент Δt\Delta t — тоже радиальный приток, либо Δt≪tp\Delta t \ll t_p: тогда Δte≈Δt\Delta t_e \approx \Delta t, и равенство верно и на участке влияния ствола;
  • постоянный дебит перед остановкой.

КВД, построенная против Δte\Delta t_e, обрабатывается всеми средствами анализа КПД (Анализ кривой падения давления): полулогарифмической прямой, формулой скина КПД, эталонными кривыми пуска. Приём опубликовал Р. Агарвал в 1980 году; ту же функцию времени давно содержит уравнение Хорнера, и эквивалентное время лишь записывает её в виде времени.

Басқа бірліктер жүйелерінде (1)
Через отношение Хорнера (десятичные логарифмы)
lg⁡Δte=lg⁡tp−lg⁡tp+ΔtΔt\lg\Delta t_e = \lg t_p - \lg\frac{t_p + \Delta t}{\Delta t}

tpt_p и Δt\Delta t в одних единицах; Δte\Delta t_e — в тех же.

Постоянных нет: полулогарифмический график против lg⁡Δte\lg\Delta t_e — это график Хорнера, отражённый слева направо и сдвинутый на lg⁡tp\lg t_p. Отношению Хорнера, равному 1, отвечает Δte=tp\Delta t_e = t_p.

  • В начале КВД Δte≈Δt\Delta t_e \approx \Delta t: при Δt=0,1 tp\Delta t = 0{,}1\,t_p эквивалентное время меньше настоящего на 9 %.
  • Δte\Delta t_e всегда меньше tpt_p и стремится к нему при бесконечной остановке; при Δt=tp\Delta t = t_p оно равно tp/2t_p/2.
  • Точка прямой при Δte=tp\Delta t_e = t_p отвечает бесконечной остановке: там читают p∗p^*, в новом пласте равное начальному давлению.
  • Чем длиннее остановка, тем теснее точки на оси Δte\Delta t_e: последние часы записи сжимаются в доли часа эквивалентного времени — то же сжатие, что и на оси Хорнера.

Раз КВД против Δte\Delta t_e — это КПД, скин считают по формуле КПД (скин по КПД) с ростом давления на прямой при Δte\Delta t_e = 1 ч. Поправка на время работы, которую добавляют к формуле скина по Хорнеру, здесь уже учтена самим пересчётом времени:

Скин по прямой против эквивалентного времени
S=1,151[Δp1чm−lg⁡kϕ μ ct rw2+2,098]S = 1{,}151\left[\frac{\Delta p_{1\text{ч}}}{m} - \lg\frac{k}{\phi\,\mu\,c_t\,r_w^2} + 2{,}098\right]
мұндағы:
  • Δp1ч\Delta p_{1\text{ч}}рост давления на прямой (или её продолжении) при Δte\Delta t_e = 1 ч, отсчитанный от pwf(Δt=0)p_{wf}(\Delta t = 0), МПа
  • mmнаклон прямой против lg⁡Δte\lg\Delta t_e на десятичный цикл — тот же, что у прямой Хорнера, МПа
  • k, ϕ, μ, ct, rwk,\ \phi,\ \mu,\ c_t,\ r_wпроницаемость (по наклону), пористость, вязкость, общая сжимаемость, радиус скважины, мД, —, мПа·с, 1/МПа, м

Бірліктер: метрические промысловые: давление — МПа, kk — мД, μ\mu — мПа·с, ctc_t — 1/МПа, rwr_w — м, время — часы

Мына жағдайда орындалады
  • условия эквивалентного времени: до остановки приток вышел на радиальный режим;
  • Δp1ч\Delta p_{1\text{ч}} читают на прямой, а не по замеру; при tpt_p меньше часа точка Δte\Delta t_e = 1 ч лежит на продолжении прямой за tpt_p.

Формула та же, что для КПД (Анализ кривой падения давления), с Δp1ч\Delta p_{1\text{ч}} вместо pi−p1чp_i - p_{1\text{ч}}. С формулой Хорнера она совпадает, если в ту добавить слагаемое lg⁡[(tp+1)/tp]\lg[(t_p + 1)/t_p]; без него скин по Хорнеру после короткой работы занижен.

Басқа бірліктер жүйелерінде (1)
Промысловые единицы (oilfield)
S=1,151[Δp1hm−lg⁡kϕ μ ct rw2+3,23]S = 1{,}151\left[\frac{\Delta p_{1\text{h}}}{m} - \lg\frac{k}{\phi\,\mu\,c_t\,r_w^2} + 3{,}23\right]

давление и mm — psi, kk — мД, μ\mu — сП, ctc_t — 1/psi, rwr_w — ft, время — часы.

3,23=−lg⁡(2,637⋅10−4)−0,809/ln⁡103{,}23 = -\lg(2{,}637\cdot10^{-4}) - 0{,}809/\ln 10, как в формуле скина по КПД; в метрических единицах 2,098=−lg⁡(3,553⋅10−3)−0,3512{,}098 = -\lg(3{,}553\cdot10^{-3}) - 0{,}351.

Заключительное восстановление давления при испытании пласта на трубах (условные данные)
Δt\Delta t, ч
0,01
Δte\Delta t_e, ч
0,00998
pwsp_{ws}, МПа
21,4875
Δp\Delta p, МПа
2,6969
Δt\Delta t, ч
0,02
Δte\Delta t_e, ч
0,0199
pwsp_{ws}, МПа
21,8868
Δp\Delta p, МПа
3,0962
Δt\Delta t, ч
0,05
Δte\Delta t_e, ч
0,0496
pwsp_{ws}, МПа
22,1990
Δp\Delta p, МПа
3,4084
Δt\Delta t, ч
0,1
Δte\Delta t_e, ч
0,0984
pwsp_{ws}, МПа
22,3875
Δp\Delta p, МПа
3,5969
Δt\Delta t, ч
0,2
Δte\Delta t_e, ч
0,194
pwsp_{ws}, МПа
22,5617
Δp\Delta p, МПа
3,7711
Δt\Delta t, ч
0,5
Δte\Delta t_e, ч
0,462
pwsp_{ws}, МПа
22,7777
Δp\Delta p, МПа
3,9871
Δt\Delta t, ч
1
Δte\Delta t_e, ч
0,857
pwsp_{ws}, МПа
22,9291
Δp\Delta p, МПа
4,1385
Δt\Delta t, ч
1,5
Δte\Delta t_e, ч
1,20
pwsp_{ws}, МПа
23,0110
Δp\Delta p, МПа
4,2204
Δt\Delta t, ч
2
Δte\Delta t_e, ч
1,50
pwsp_{ws}, МПа
23,0651
Δp\Delta p, МПа
4,2745
Δt\Delta t, ч
3
Δte\Delta t_e, ч
2,00
pwsp_{ws}, МПа
23,1348
Δp\Delta p, МПа
4,3442
Δt\Delta t, ч
4,5
Δte\Delta t_e, ч
2,57
pwsp_{ws}, МПа
23,1956
Δp\Delta p, МПа
4,4050
Δt\Delta t, ч
6
Δte\Delta t_e, ч
3,00
pwsp_{ws}, МПа
23,2328
Δp\Delta p, МПа
4,4422
Δt\Delta t, ч
9
Δte\Delta t_e, ч
3,60
pwsp_{ws}, МПа
23,2768
Δp\Delta p, МПа
4,4862

Основной приток tpt_p = 6 ч с дебитом 64 м³/сут нефти; забойное давление в момент закрытия клапана pwf(Δt=0)p_{wf}(\Delta t = 0) = 18,7906 МПа; Δp=pws−pwf(Δt=0)\Delta p = p_{ws} - p_{wf}(\Delta t = 0). Данные рассчитаны по модели однородного бесконечного пласта со скином и влиянием ствола под клапаном у пакера и округлены до 0,0001 МПа — разрешения кварцевого манометра.

Давления КВД против времени после остановки и против эквивалентного времени0,111022,422,622,823,023,223,4Время, чДавление pws, МПапротив Δtпротив ΔteΔte = tp

Полулогарифмический график: время от 0,1 до 10 часов, давление от 22,3 до 23,5 МПа. Тёмные точки — те же давления против эквивалентного времени — лежат на оранжевой штриховой прямой, которая поднимается от 22,4 МПа при 0,1 часа до 23,4 МПа при 6 часах; конец прямой подписан: эквивалентное время равно времени работы. Светлые кружки против настоящего времени до 1 часа почти совпадают с тёмными точками, а после 1–2 часов отходят вправо и вниз от прямой; последний кружок, при 9 часах, на 23,28 МПа.

ЕсептеуДавления КВД из таблицы на полулогарифмическом графике дважды: против времени после остановки Δt\Delta t (светлые кружки — график MDH) и против эквивалентного времени Δte\Delta t_e (тёмные точки). Против Δt\Delta t точки после 1–2 ч изгибаются вниз; против Δte\Delta t_e они ложатся на прямую радиального притока (штрих), проведённую в примере ниже через замеры 1,5 и 6 ч. Её продолжение до Δte=tp\Delta t_e = t_p = 6 ч даёт начальное давление.
Есептеу мысалыТри способа обработать одну КВД
Бастапқы деректер (шартты)
Основной приток
qq = 64 м³/сут нефти, tpt_p = 6 ч; заключительное восстановление — 9 ч
Пласт и флюид
BB = 1,21; μ\mu = 1,35 мПа·с; hh = 9 м; ϕ\phi = 0,18; ct=1,5⋅10−3c_t = 1{,}5\cdot10^{-3} МПа⁻¹; rwr_w = 0,108 м
Замеры
таблица выше; прямую проводим через замеры 1,5 ч (23,0110 МПа) и 6 ч (23,2328 МПа) — оба после окончания влияния ствола
Шешім
  1. Разность давлений: 23,2328−23,0110=0,221823{,}2328 - 23{,}0110 = 0{,}2218 МПа; qBμ=64⋅1,21⋅1,35=104,544qB\mu = 64\cdot1{,}21\cdot1{,}35 = 104{,}544.
  2. Метод MDH: m=0,2218/lg⁡(6/1,5)=0,2218/0,60206=0,3684m = 0{,}2218/\lg(6/1{,}5) = 0{,}2218/0{,}60206 = 0{,}3684 МПа на цикл; k=2,149⋅104,544/(0,3684⋅9)=67,8k = 2{,}149\cdot104{,}544/(0{,}3684\cdot9) = 67{,}8 мД.
  3. Метод Хорнера: абсциссы lg⁡(7,5/1,5)=0,69897\lg(7{,}5/1{,}5) = 0{,}69897 и lg⁡(12/6)=0,30103\lg(12/6) = 0{,}30103; m=0,2218/(0,69897−0,30103)=0,5574m = 0{,}2218/(0{,}69897 - 0{,}30103) = 0{,}5574 МПа на цикл; k=2,149⋅104,544/(0,5574⋅9)=44,8k = 2{,}149\cdot104{,}544/(0{,}5574\cdot9) = 44{,}8 мД.
  4. Эквивалентное время: Δte=6⋅1,5/7,5=1,2\Delta t_e = 6\cdot1{,}5/7{,}5 = 1{,}2 ч и 6⋅6/12=36\cdot6/12 = 3 ч; m=0,2218/lg⁡(3/1,2)=0,2218/0,39794=0,5574m = 0{,}2218/\lg(3/1{,}2) = 0{,}2218/0{,}39794 = 0{,}5574 МПа — тот же наклон, что у Хорнера, потому что lg⁡2,5=lg⁡5−lg⁡2\lg 2{,}5 = \lg 5 - \lg 2.
  5. Скин по формуле КПД. Рост давления в замере 1,5 ч (Δte\Delta t_e = 1,2 ч): 23,0110−18,7906=4,220423{,}0110 - 18{,}7906 = 4{,}2204 МПа; на прямой при Δte\Delta t_e = 1 ч — Δp1ч=4,2204−0,5574⋅lg⁡1,2=4,2204−0,5574⋅0,07918=4,1763\Delta p_{1\text{ч}} = 4{,}2204 - 0{,}5574\cdot\lg 1{,}2 = 4{,}2204 - 0{,}5574\cdot0{,}07918 = 4{,}1763 МПа; lg⁡kϕμctrw2=lg⁡44,80,18⋅1,35⋅1,5⋅10−3⋅0,1082=7,0227\lg\dfrac{k}{\phi\mu c_tr_w^2} = \lg\dfrac{44{,}8}{0{,}18\cdot1{,}35\cdot1{,}5\cdot10^{-3}\cdot0{,}108^2} = 7{,}0227; S=1,151 (4,1763/0,5574−7,0227+2,098)=2,96S = 1{,}151\,(4{,}1763/0{,}5574 - 7{,}0227 + 2{,}098) = 2{,}96.
  6. Начальное давление — на прямой при Δte=tp\Delta t_e = t_p = 6 ч. Рост давления в замере 6 ч (Δte\Delta t_e = 3 ч): 23,2328−18,7906=4,442223{,}2328 - 18{,}7906 = 4{,}4422 МПа; на прямой при 6 ч — 4,4422+0,5574⋅lg⁡(6/3)=4,61004{,}4422 + 0{,}5574\cdot\lg(6/3) = 4{,}6100 МПа; pi=18,7906+4,6100=23,4006p_i = 18{,}7906 + 4{,}6100 = 23{,}4006 МПа.
Жауап: По Хорнеру и по эквивалентному времени k≈44,8k \approx 44{,}8 мД, скин 2,96, pi≈23,40p_i \approx 23{,}40 МПа; метод MDH завысил проницаемость до 67,8 мД — на 51 %.

Бұл нені білдіреді. Данные рассчитаны по модели с проницаемостью 45 мД, скином 3 и начальным давлением 23,4 МПа: Хорнер и эквивалентное время вернули их с погрешностью меньше процента по проницаемости и 0,04 по скину — это след не вполне закончившегося влияния ствола в замере 1,5 ч и округления замеров. MDH ошибся не из-за данных: к 6 ч местный наклон по lg⁡Δt\lg\Delta t составляет лишь 6/(6+6)=0,56/(6 + 6) = 0{,}5 полного, и хорда через 1,5 и 6 ч в полтора раза положе прямой. Скин по формуле КПД при Δte\Delta t_e = 1 ч совпадает со скином по формуле Хорнера с поправкой lg⁡[(tp+1)/tp]\lg[(t_p + 1)/t_p]; без поправки формула Хорнера дала бы на 0,08 меньше. Данные примера условные.

Производная по эквивалентному времени

Для диагностического графика КВД после короткой работы производную берут не по ln⁡Δt\ln\Delta t, а по ln⁡Δte\ln\Delta t_e. Поскольку ln⁡Δte=ln⁡tp−ln⁡[(tp+Δt)/Δt]\ln\Delta t_e = \ln t_p - \ln[(t_p + \Delta t)/\Delta t], это та же производная, что по функции Хорнера, — производная по функции суперпозиции времени для одного режима. Из обычной производной её получают умножением:

Производная по эквивалентному времени
Δpe′= d ΔpBU dln⁡Δte=tp+Δttp⋅ d ΔpBU dln⁡Δt\Delta p'_e = \frac{\dd\,\Delta p_{BU}}{\dd\ln\Delta t_e} = \frac{t_p + \Delta t}{t_p}\cdot\frac{\dd\,\Delta p_{BU}}{\dd\ln\Delta t}
мұндағы:
  • Δpe′\Delta p'_eлогарифмическая производная КВД по эквивалентному времени, Па
  •  d ΔpBU/ dln⁡Δt\dd\,\Delta p_{BU}/\dd\ln\Delta tобычная логарифмическая производная КВД по времени после остановки, Δp′\Delta p', Па
  • tp, Δtt_p,\ \Delta tвремя работы и время после остановки, с

Бірліктер: любая согласованная система; производные — в единицах давления

Мына жағдайда орындалады
  • tpt_p — время работы с постоянным дебитом перед остановкой; при меняющемся дебите производную берут по функции суперпозиции всей истории (Эквивалентное и синтетическое время работы);
  • логарифм натуральный; связь следует из  dln⁡Δte/ dln⁡Δt=tp/(tp+Δt)\dd\ln\Delta t_e/\dd\ln\Delta t = t_p/(t_p + \Delta t).

Множитель (tp+Δt)/tp(t_p + \Delta t)/t_p ровно гасит множитель tp/(tp+Δt)t_p/(t_p + \Delta t) местного наклона: при радиальном притоке Δpe′=me\Delta p'_e = m_e при любом времени работы. Пока Δt≪tp\Delta t \ll t_p, две производные совпадают; в конце длинной КВД первая больше второй во столько раз, во сколько tp+Δtt_p + \Delta t больше tpt_p. Во столько же раз при тех же соседних замерах усиливается и шум.

Строят её двумя способами. Можно отложить и ΔpBU\Delta p_{BU}, и Δpe′\Delta p'_e против Δte\Delta t_e — так, как Агарвал предложил для давления: тогда КВД накладывается на эталонные кривые пуска. В программах интерпретации обычно откладывают против настоящего времени Δt\Delta t: поздние участки записи не сжимаются к tpt_p, а радиальный приток по-прежнему виден как полка. Сами значения производной в обоих случаях одни и те же — различается только ось времени.

По замерам её считают той же формулой трёх точек, что и обычную производную (как считают производную по замерам), только вместо ln⁡Δt\ln\Delta t берут ln⁡Δte\ln\Delta t_e. Интервал дифференцирования тоже отсчитывают по ln⁡Δte\ln\Delta t_e: в конце записи замеры на этой оси сгущаются, и при том же интервале соседние точки отстоят друг от друга по настоящему времени всё дальше.

Диагностический график КВД примера: рост давления и две производные0,010,11100,1110Время после остановки Δt, чΔp и Δp′, МПаΔpпо ln Δteпо ln Δtполка me = 0,241 МПа

График в двойном логарифмическом масштабе: время после остановки от 0,01 до 10 часов, давление и производная от 0,05 до 10 МПа. Точки роста давления медленно поднимаются от 2,7 до 4,5 МПа. Обе производные начинаются около 0,75 МПа и спускаются к полке 0,241 МПа, показанной штрихом. Оранжевая производная по эквивалентному времени остаётся на полке до 9 часов; серая производная по настоящему времени после 0,5 часа уходит вниз и к 9 часам опускается до 0,1 МПа.

ЕсептеуДиагностический график КВД из таблицы: рост давления (точки) и две производные, рассчитанные по модели урока и отложенные против настоящего времени Δt\Delta t: по ln⁡Δt\ln\Delta t (серая) и по ln⁡Δte\ln\Delta t_e (оранжевая); штрих — полка mem_e. После 0,5 ч серая производная уходит вниз, оранжевая остаётся на полке до конца записи.
Есептеу мысалыПроизводная по замерам: по Δt и по Δt_e
Бастапқы деректер (шартты)
Замеры
из таблицы: рост давления 4,2204 МПа через 1,5 ч, 4,3442 МПа через 3 ч, 4,4422 МПа через 6 ч, 4,4862 МПа через 9 ч
Испытание
tpt_p = 6 ч; qBμqB\mu = 104,544; hh = 9 м
Шешім
  1. По ln⁡Δt\ln\Delta t в точке 3 ч: соседи 1,5 и 6 ч отстоят от неё на ln⁡2=0,6931\ln 2 = 0{,}6931 каждый, производная — центральная разность: Δp′=(4,4422−4,2204)/(2⋅0,6931)=0,1600\Delta p' = (4{,}4422 - 4{,}2204)/(2\cdot0{,}6931) = 0{,}1600 МПа.
  2. По ln⁡Δte\ln\Delta t_e: эквивалентные времена тех же замеров 1,2; 2 и 3 ч; расстояния ln⁡(2/1,2)=0,5108\ln(2/1{,}2) = 0{,}5108 и ln⁡(3/2)=0,4055\ln(3/2) = 0{,}4055; наклоны хорд sl=(4,3442−4,2204)/0,5108=0,2424s_l = (4{,}3442 - 4{,}2204)/0{,}5108 = 0{,}2424 МПа и sr=(4,4422−4,3442)/0,4055=0,2417s_r = (4{,}4422 - 4{,}3442)/0{,}4055 = 0{,}2417 МПа.
  3. Взвешенное среднее хорд: Δpe′=(0,2424⋅0,4055+0,2417⋅0,5108)/(0,5108+0,4055)=0,2420\Delta p'_e = (0{,}2424\cdot0{,}4055 + 0{,}2417\cdot0{,}5108)/(0{,}5108 + 0{,}4055) = 0{,}2420 МПа. Проверка множителем: 0,1600⋅(6+3)/6=0,24000{,}1600\cdot(6 + 3)/6 = 0{,}2400 МПа.
  4. В точке 6 ч (соседи 3 и 9 ч) та же формула даёт Δp′=0,1206\Delta p' = 0{,}1206 МПа по ln⁡Δt\ln\Delta t и Δpe′=0,2414\Delta p'_e = 0{,}2414 МПа по ln⁡Δte\ln\Delta t_e.
  5. Проводимость по производной по эквивалентному времени: kh=0,933⋅104,544/0,2420=403kh = 0{,}933\cdot104{,}544/0{,}2420 = 403 мД·м, k=44,8k = 44{,}8 мД. По производной по ln⁡Δt\ln\Delta t в той же точке вышло бы 610 мД·м — в полтора раза больше.
Жауап: Производная по ln⁡Δt\ln\Delta t падает — 0,1600 МПа в 3 ч и 0,1206 МПа в 6 ч, а производная по ln⁡Δte\ln\Delta t_e стоит на полке: 0,2420 и 0,2414 МПа; по ней k≈44,8k \approx 44{,}8 мД.

Бұл нені білдіреді. Спад производной по ln⁡Δt\ln\Delta t — в точности множитель tp/(tp+Δt)t_p/(t_p + \Delta t): 6/9 = 0,67 в 3 ч и 6/12 = 0,5 в 6 ч. Никакой границы в пласте нет, и производная по эквивалентному времени это показывает. Её небольшое превышение над полкой модели (0,2409 МПа) — хвост влияния ствола и округление замеров. Если бы спад остался и на производной по Δte\Delta t_e, тогда стоило бы искать его причину в пласте. Данные примера условные.

Эталонные кривые, сжатие времени и планирование притока

Эталонные кривые скважины со стволом и скином (Эталонные кривые) рассчитаны для пуска с постоянным дебитом. Для КВД пришлось бы строить своё семейство на каждое безразмерное время работы tpD=ktp/(ϕμctrw2)t_{pD} = kt_p/(\phi\mu c_tr_w^2) — и знать проницаемость ещё до совмещения, чтобы выбрать нужное. Эквивалентное время снимает эту задачу: рост давления и производную Δpe′\Delta p'_e откладывают против Δte\Delta t_e и совмещают с теми же кривыми, что и КПД. Совмещение по давлению даёт khkh, по времени — коэффициент влияния ствола CC (в формулу подставляют Δte\Delta t_e вместо Δt\Delta t), параметр кривой — скин.

Цена — сжатие времени. Ось Δte\Delta t_e кончается на tpt_p: сколько бы ни длилась остановка, КВД воспроизводит КПД только до времени, равного времени работы, а поздние часы записи сжимаются к этой границе. В примере урока 9 ч восстановления дали эквивалентное время всего 3,6 ч. Отсюда простое правило планирования: чтобы КВД показала эквивалентное время TT, остановка должна длиться

Длительность остановки для заданного эквивалентного времени
Δt=tp Ttp−T,T<tp\Delta t = \frac{t_p\,T}{t_p - T},\qquad T < t_p
мұндағы:
  • TTэквивалентное время, до которого КВД должна воспроизвести КПД, ч
  • tpt_pвремя работы перед остановкой, ч
  • Δt\Delta tнужная длительность остановки, ч

Бірліктер: любые одинаковые единицы времени

Мына жағдайда орындалады
  • бесконечный пласт, радиальный приток; формула — обращение Δte(Δt)\Delta t_e(\Delta t);
  • TT выбирают по тому, что нужно увидеть: например, цикл полки производной после окончания влияния ствола.

Чем ближе TT к tpt_p, тем длиннее остановка: при T=tp/2T = t_p/2 она равна tpt_p, при T=0,8 tpT = 0{,}8\,t_p — уже 4tp4t_p. Всё испытание, tp+Δtt_p + \Delta t, короче всего, когда приток длится 2T2T и остановка — тоже 2T2T: производная суммы по tpt_p равна 1−T2/(tp−T)21 - T^2/(t_p - T)^2 и обращается в ноль при tp=2Tt_p = 2T.

Есептеу мысалыКакой приток выбрать
Бастапқы деректер (шартты)
Цель
увидеть по КВД полку производной до эквивалентного времени TT = 5 ч (влияние ствола кончается за полчаса)
Варианты
основной приток 6, 10 или 15 ч
Шешім
  1. tpt_p = 6 ч: Δt=6⋅5/(6−5)=30\Delta t = 6\cdot5/(6 - 5) = 30 ч; всё испытание 36 ч.
  2. tpt_p = 10 ч (=2T= 2T): Δt=10⋅5/(10−5)=10\Delta t = 10\cdot5/(10 - 5) = 10 ч; всего 20 ч.
  3. tpt_p = 15 ч: Δt=15⋅5/(15−5)=7,5\Delta t = 15\cdot5/(15 - 5) = 7{,}5 ч; всего 22,5 ч.
Жауап: Короче всего — приток 10 ч и остановка 10 ч: 20 ч вместо 36 ч при шестичасовом притоке.

Бұл нені білдіреді. Удлинять приток выгодно, пока остановка сокращается быстрее, чем он растёт. Но час притока обходится иначе, чем час остановки: продукцию при испытании обычно сжигают на факеле, растёт риск выноса песка, программа испытания ограничивает отбор. Расчёт показывает, чего стоит каждый лишний час эквивалентного времени; окончательное решение — компромисс. Данные примера условные.

Когда эквивалентное время не помогает

Эквивалентное время выведено из логарифмического приближения радиального притока в бесконечном пласте. Где оно не годится хотя бы для одного из трёх слагаемых суперпозиции, пересчёт времени перестаёт быть точным.

  • Приток не вышел на радиальный режим. Если скважину закрыли на переходном участке влияния ствола, забойное давление в момент остановки ещё не легло на полулогарифмическую прямую. Тогда КВД против Δte\Delta t_e не совпадает с кривой КПД: поздние точки сжимаются к Δte≈tp\Delta t_e \approx t_p, где эталонная кривая ещё на горбе, а скин по давлению в момент остановки занижен. Производная по ln⁡Δte\ln\Delta t_e в конце длинной КВД всё же выходит на полку — её откладывают против Δt\Delta t, а запись целиком описывают моделью со всей историей дебита (пример ниже).
  • Границы. Если до границ дошло возмущение от пуска или от остановки, логарифмическое приближение неверно для ΔpDD(tp+Δt)\Delta p_{DD}(t_p + \Delta t). В замкнутой области КВД выходит на среднее давление (Восстановление давления в ограниченной области и среднее пластовое давление), а КПД — на прямую единичного наклона (Исследование на ограниченность пласта): никакое время одну в другую не превратит. У непроницаемого разлома производная по ln⁡Δte\ln\Delta t_e со временем поднимается к удвоенному уровню, но против Δte\Delta t_e этот подъём сжат к tpt_p — поэтому границы ищут на графике против Δt\Delta t.
  • Меняющийся дебит. Формула выведена для одного режима. При нескольких режимах берут функцию суперпозиции всей истории или подбирают одно условное время работы (Эквивалентное и синтетическое время работы). Опасная ошибка — после долгой работы и короткой промежуточной остановки взять за tpt_p только последний приток: время сожмётся слишком сильно, и производная исказится.
  • Газовые скважины. Уравнение для газа нелинейно, поэтому суперпозицию и эквивалентное время применяют к псевдодавлению (псевдодавление газа), а при больших изменениях давления — ещё и к псевдовремени, учитывающему изменение μct\mu c_t с давлением.
  • Меняющиеся скин и влияние ствола. Очистка призабойной зоны во время притока или перераспределение газа и жидкости в стволе после остановки нарушают линейность задачи (Неидеальное влияние ствола, очистка скважины и меняющийся скин).
Есептеу мысалыПриток короче влияния ствола
Бастапқы деректер (шартты)
Скважина
тот же пласт и тот же скин, но после освоения: отработка 1 ч с закрытием на устье, в трубах свободный газ, CC = 0,15 м³/МПа
Замеры
pwf(Δt=0)p_{wf}(\Delta t = 0) = 19,3643 МПа; через 10 ч после остановки — 23,3761 МПа, через 40 ч — 23,3940 МПа
Шешім
  1. Начало полулогарифмической прямой по оценке из урока о влиянии ствола (длительность влияния ствола): t≈44,8 (60+3,5⋅3) 1,35⋅0,1545⋅9=1,58t \approx 44{,}8\,(60 + 3{,}5\cdot3)\,\dfrac{1{,}35\cdot0{,}15}{45\cdot9} = 1{,}58 ч — дольше всего притока.
  2. Прямая Хорнера через 10 и 40 ч: абсциссы lg⁡(11/10)=0,04139\lg(11/10) = 0{,}04139 и lg⁡(41/40)=0,01072\lg(41/40) = 0{,}01072; m=(23,3940−23,3761)/(0,04139−0,01072)=0,5836m = (23{,}3940 - 23{,}3761)/(0{,}04139 - 0{,}01072) = 0{,}5836 МПа на цикл; k=2,149⋅104,544/(0,5836⋅9)=42,8k = 2{,}149\cdot104{,}544/(0{,}5836\cdot9) = 42{,}8 мД.
  3. Продолжение прямой: p∗=23,3940+0,5836⋅0,01072=23,4003p^* = 23{,}3940 + 0{,}5836\cdot0{,}01072 = 23{,}4003 МПа.
  4. Скин по формуле Хорнера с поправкой на время работы lg⁡[(1+1)/1]=0,30103\lg[(1 + 1)/1] = 0{,}30103: p1ч=23,4003−0,5836⋅0,30103=23,2246p_{1\text{ч}} = 23{,}4003 - 0{,}5836\cdot0{,}30103 = 23{,}2246 МПа, lg⁡kϕμctrw2=7,0029\lg\dfrac{k}{\phi\mu c_tr_w^2} = 7{,}0029; S=1,151 [(23,2246−19,3643)/0,5836−7,0029+2,098+0,30103]=2,31S = 1{,}151\,[(23{,}2246 - 19{,}3643)/0{,}5836 - 7{,}0029 + 2{,}098 + 0{,}30103] = 2{,}31.
  5. Эквивалентное время замеров: 1⋅10/11=0,9091\cdot10/11 = 0{,}909 ч и 1⋅40/41=0,9761\cdot40/41 = 0{,}976 ч — оба у самого предела tpt_p = 1 ч.
Жауап: Проницаемость 42,8 мД и p∗p^* = 23,40 МПа близки к правильным (45 мД и 23,4 МПа), а скин 2,31 вместо 3.

Бұл нені білдіреді. Поздние замеры уже лежат в радиальном притоке, поэтому наклон и p∗p^* почти верны: Хорнеру не нужно, чтобы на прямую успел выйти приток. Скину нужно: формула предполагает, что к моменту остановки забойное давление лежало на полулогарифмической прямой КПД, а оно было ещё выше её — пласт не успел «почувствовать» полную депрессию. На эталонных кривых против Δte\Delta t_e эти замеры легли бы у Δte≈1\Delta t_e \approx 1 ч, где кривая КПД скважины ещё на горбе (её производная там в 1,9 раза выше полки), а производная КВД по ln⁡Δte\ln\Delta t_e уже почти на полке — совмещение было бы бессмысленным. Правильный путь — модель со стволом, скином и историей дебита, подобранная ко всей записи. Данные примера условные.

Терминдер мен белгілер

Время работы tpt_p
длительность работы с постоянным дебитом перед остановкой; при меняющемся дебите — условное время, заменяющее историю (Эквивалентное и синтетическое время работы).
Эквивалентное время Δte\Delta t_e
tpΔt/(tp+Δt)t_p\Delta t/(t_p + \Delta t): время после пуска, через которое КПД дала бы тот же перепад, что КВД через Δt\Delta t после остановки; всегда меньше tpt_p.
Рост давления после остановки ΔpBU\Delta p_{BU}
pws(Δt)−pwf(Δt=0)p_{ws}(\Delta t) - p_{wf}(\Delta t = 0); ограничен сверху снижением давления за время работы pi−pwf(tp)p_i - p_{wf}(t_p).
Производная по эквивалентному времени Δpe′\Delta p'_e
 dΔpBU/ dln⁡Δte\dd\Delta p_{BU}/\dd\ln\Delta t_e; равна производной по функции Хорнера, при радиальном притоке — полке mem_e при любом tpt_p.
Сжатие времени
свойство функции Хорнера и эквивалентного времени: поздние часы КВД умещаются в узкий интервал у Δte=tp\Delta t_e = t_p.
Десуперпозиция
восстановление КПД по КВД прибавлением продолжения снижения давления, которое дала бы работающая скважина (дополнительный раздел).

Жиі қателіктер

Қате

«Производная в конце КВД упала вдвое — у пласта граница постоянного давления».

Шын мәнінде

Сначала сравните длительность записи со временем работы: после короткой работы производная по ln⁡Δt\ln\Delta t падает в (tp+Δt)/tp(t_p + \Delta t)/t_p раз и без всяких границ. Покажет, есть ли спад на самом деле, производная по эквивалентному времени.

Қате

«Метод MDH проще, а результат тот же».

Шын мәнінде

Только при Δt≪tp\Delta t \ll t_p. В примере урока после шестичасового притока MDH по замерам 1,5 и 6 ч завысил проницаемость в полтора раза.

Қате

«Чем дольше КВД, тем дальше в пласт заглядывает исследование, сколько бы ни длился приток».

Шын мәнінде

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

Қате

«Эквивалентное время годится для любой КВД».

Шын мәнінде

Только если до остановки приток вышел на радиальный режим и границы за время исследования не сказались. После часовой отработки с газом в трубах скин по КВД в примере урока получился 2,3 вместо 3.

Қате

«Эквивалентное время и эквивалентное время работы — одно и то же».

Шын мәнінде

Δte\Delta t_e пересчитывает время после остановки, Np/qN_p/q заменяет время работы при меняющемся дебите. Первое подставляют вместо Δt\Delta t, второе — вместо tpt_p.

Толық: тереңдетілген

Вывод эквивалентного времени и его погрешность

При радиальном притоке в бесконечном пласте ΔpDD(t)=me [ln⁡(at)+0,809+2S]\Delta p_{DD}(t) = m_e\,[\ln(at) + 0{,}809 + 2S], где a=k/(ϕμctrw2)a = k/(\phi\mu c_tr_w^2) переводит время в безразмерное. Подставим это в ΔpBU=ΔpDD(Δt)+ΔpDD(tp)−ΔpDD(tp+Δt)\Delta p_{BU} = \Delta p_{DD}(\Delta t) + \Delta p_{DD}(t_p) - \Delta p_{DD}(t_p + \Delta t): из трёх слагаемых 0,809+2S0{,}809 + 2S остаётся одно, а логарифмы складываются в ln⁡a+ln⁡[tpΔt/(tp+Δt)]\ln a + \ln[t_p\Delta t/(t_p + \Delta t)]. Получается ΔpBU=me[ln⁡(aΔte)+0,809+2S]=ΔpDD(Δte)\Delta p_{BU} = m_e[\ln(a\Delta t_e) + 0{,}809 + 2S] = \Delta p_{DD}(\Delta t_e).

Где приближение не точно, полезно знать, насколько. Обозначим ε(t)\varepsilon(t) отклонение КПД от её полулогарифмической прямой: на участке влияния ствола оно отрицательно, при радиальном притоке равно нулю. Из той же суперпозиции:

Погрешность эквивалентного времени
ΔpBU(Δt)−ΔpDD(Δte)=[ε(Δt)−ε(Δte)]+[ε(tp)−ε(tp+Δt)]\Delta p_{BU}(\Delta t) - \Delta p_{DD}(\Delta t_e) = \big[\varepsilon(\Delta t) - \varepsilon(\Delta t_e)\big] + \big[\varepsilon(t_p) - \varepsilon(t_p + \Delta t)\big]
мұндағы:
  • ε(t)\varepsilon(t)отклонение КПД от её полулогарифмической прямой: ΔpDD(t)−me[ln⁡(at)+0,809+2S]\Delta p_{DD}(t) - m_e[\ln(at) + 0{,}809 + 2S], Па
  • Δt, Δte, tp\Delta t,\ \Delta t_e,\ t_pвремя после остановки, эквивалентное время, время работы, с

Бірліктер: любая согласованная система

Мына жағдайда орындалады
  • бесконечный пласт; влияние ствола и скин постоянны, задача линейна.

Первая скобка исчезает, когда Δt≪tp\Delta t \ll t_p (тогда Δte≈Δt\Delta t_e \approx \Delta t) или когда оба времени уже в радиальном притоке; вторая — когда на радиальный режим вышел сам приток. Поэтому главное условие — приток в несколько раз длиннее влияния ствола. В испытании урока (tpt_p = 6 ч, влияние ствола — десятые доли часа) разность по всей записи не больше 0,00085 МПа; после часовой отработки с газом в трубах (пример в разделе о пределах) велики обе скобки: и время работы, и эквивалентное время приходятся на переход от влияния ствола к радиальному притоку. Через 10 ч после остановки первая скобка равна 0,154 МПа, вторая −0,132 МПа, и КВД лежит выше кривой КПД на их сумму, 0,022 МПа.

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

Ту же поправку можно внести не во время, а в давление. Перепишем суперпозицию как ΔpDD(Δt)=ΔpBU(Δt)+[ΔpDD(tp+Δt)−ΔpDD(tp)]\Delta p_{DD}(\Delta t) = \Delta p_{BU}(\Delta t) + [\Delta p_{DD}(t_p + \Delta t) - \Delta p_{DD}(t_p)]: к росту давления после остановки прибавляют то, на сколько снизилось бы давление в работающей скважине за время Δt\Delta t. Получается настоящая КПД, растянутая до полной длительности остановки. Этот приём называют десуперпозицией (desuperposition).

Десуперпозиция при радиальном притоке
ΔpDD(Δt)≈ΔpBU(Δt)+m lg⁡tp+Δttp\Delta p_{DD}(\Delta t) \approx \Delta p_{BU}(\Delta t) + m\,\lg\frac{t_p + \Delta t}{t_p}
мұндағы:
  • ΔpDD(Δt)\Delta p_{DD}(\Delta t)восстановленное снижение давления, как после пуска, МПа
  • mmнаклон полулогарифмической прямой на десятичный цикл, МПа
  • ΔpBU, tp, Δt\Delta p_{BU},\ t_p,\ \Delta tрост давления после остановки, время работы и время после остановки, МПа, ч, ч

Бірліктер: давление — в любых единицах (МПа), время — в любых одинаковых единицах

Мына жағдайда орындалады
  • в момент остановки и после него скважина «продолжила бы» работать в режиме радиального притока; иначе поправку берут из записанной КПД или из модели.

Наклон mm заранее неизвестен, поэтому расчёт ведут итерациями: первый проход без поправки (это метод MDH), затем поправка с найденным наклоном, и так далее. Для пары точек ошибка наклона на каждом проходе умножается на долю, которую разность поправок составляет от разности логарифмов времени, — в примере урока (замеры 1,5 и 6 ч) на 0,34, и итерации сходятся к наклону Хорнера. Эквивалентное время делает то же сразу, без итераций, поэтому им пользуются чаще. Десуперпозиция удобна, когда продолжение КПД берут не из формулы, а из самой записи — например, из тренда давления в разбуренной залежи (Исследования в разбуренной залежи: метод Слайдера).

Басқа тақырыптармен байланыс

Өзіңізді тексеріңіз

  1. 1

    Скважина работала 10 ч, КВД записана за 15 ч, границ в пласте нет. Во сколько раз производная по ln⁡Δt\ln\Delta t в конце записи ниже полки радиального притока?

    1. А.Не ниже: при радиальном притоке производная постоянна
    2. Б.В 15 раз
    3. В.В 2,5 раза
    4. Г.В 1,5 раза
    Жауабын көрсету

    Дұрыс жауап: В. В 2,5 раза

    Местный наклон по ln⁡Δt\ln\Delta t равен me tp/(tp+Δt)=me⋅10/25=0,4 mem_e\,t_p/(t_p + \Delta t) = m_e\cdot10/25 = 0{,}4\,m_e — в 2,5 раза ниже полки. Производная по эквивалентному времени в этой точке осталась бы на полке.

  2. 2

    Скважина работала 8 ч. Вычислите эквивалентное время замеров через 2, 8 и 40 ч после остановки. Что это значит для анализа поздних замеров?

    Жауабын көрсету

    Δte=8⋅2/10=1,6\Delta t_e = 8\cdot2/10 = 1{,}6 ч, 8⋅8/16=48\cdot8/16 = 4 ч, 8⋅40/48≈6,678\cdot40/48 \approx 6{,}67 ч. Последние 32 ч записи (от 8 до 40 ч) уместились в интервал эквивалентного времени от 4 до 6,67 ч — меньше четверти логарифмического цикла. Поздние замеры сжаты к tpt_p = 8 ч: по ним КВД воспроизводит КПД лишь до 6,67 ч, а для графика производной их лучше откладывать против Δt\Delta t.

  3. 3

    Какой способ обработки КВД после короткой работы даёт правильный наклон радиального притока?

    1. А.Любой: наклон от времени работы не зависит
    2. Б.Метод MDH по точкам после окончания влияния ствола
    3. В.Совмещение с эталонными кривыми пуска против Δt\Delta t
    4. Г.График Хорнера или график против lg⁡Δte\lg\Delta t_e
    Жауабын көрсету

    Дұрыс жауап: Г. График Хорнера или график против lg⁡Δte\lg\Delta t_e

    Функция Хорнера и эквивалентное время содержат время работы, поэтому при радиальном притоке точки ложатся на прямую с правильным наклоном. MDH и совмещение против Δt\Delta t его не учитывают: наклон по lg⁡Δt\lg\Delta t занижен в (tp+Δt)/tp(t_p + \Delta t)/t_p раз, проницаемость завышена.

  4. 4

    После 4-часового притока по КВД построили прямую против lg⁡Δte\lg\Delta t_e: наклон 0,42 МПа на цикл, рост давления на прямой при Δte\Delta t_e = 1 ч — 3,70 МПа, забойное давление в момент остановки — 17,30 МПа. Пласт новый, границ не видно. Найдите начальное давление и объясните, почему его читают при Δte\Delta t_e = 4 ч.

    Жауабын көрсету

    На прямой при Δte=tp\Delta t_e = t_p = 4 ч рост давления 3,70+0,42lg⁡4=3,70+0,253≈3,953{,}70 + 0{,}42\lg 4 = 3{,}70 + 0{,}253 \approx 3{,}95 МПа, поэтому pi≈17,30+3,95=21,25p_i \approx 17{,}30 + 3{,}95 = 21{,}25 МПа. Значение Δte=tp\Delta t_e = t_p отвечает бесконечно долгой остановке (отношению Хорнера, равному 1): в бесконечном пласте давление вернулось бы к начальному, а рост давления — к pi−pwf(tp)p_i - p_{wf}(t_p).

  5. 5

    Почему производную по эквивалентному времени в программах интерпретации обычно откладывают против настоящего времени Δt\Delta t, а не против Δte\Delta t_e?

    Жауабын көрсету

    Значения производной от выбора оси не зависят. Против Δte\Delta t_e поздняя часть длинной КВД сжата у tpt_p: признаки границ, наступающие при Δt\Delta t больше tpt_p, умещаются в узкую полосу, и их трудно разглядеть. Против Δt\Delta t поздние замеры остаются на своих местах, а радиальный приток по-прежнему виден как полка.

  6. 6

    Отработка новой скважины длилась 40 минут с закрытием на устье; ствол частично заполнен газом, и по оценке полулогарифмическая прямая на КПД началась бы только через 2 ч. Что можно надёжно получить из последующей КВД длиной 30 ч?

    1. А.Проводимость и p∗p^* по поздним точкам; скин по давлению в момент остановки будет занижен
    2. Б.Скин точно, проводимость — нет
    3. В.Всё, если обработать запись против Δte\Delta t_e
    4. Г.Ничего
    Жауабын көрсету

    Дұрыс жауап: А. Проводимость и p∗p^* по поздним точкам; скин по давлению в момент остановки будет занижен

    Поздние замеры КВД лежат в радиальном притоке, и по ним наклон и p∗p^* находятся. Скин требует, чтобы забойное давление в момент остановки лежало на полулогарифмической прямой КПД, а приток до неё не дошёл: в примере урока после часовой отработки скин получился 2,3 вместо 3. Пересчёт против Δte\Delta t_e здесь не помогает — нужна модель со всей историей дебита.

Жаңартылған: 19 октября 2026 г.