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

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

ГДИСОграниченная область дренирования Тереңдетілген 22 мин оқу

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

Среднее пластовое давление области дренирования (average reservoir pressure, pˉ\bar p) — давление, которое установилось бы в замкнутой области, если бы отбор из неё прекратили и дождались, пока давление выровняется; при постоянной сжимаемости оно равно среднему по объёму пор. Кривая восстановления давления (КВД) в такой области со временем выходит именно на pˉ\bar p, а не на начальное давление. Обычно столько не ждут: продолжают прямую Хорнера до давления p∗p^* и вычитают из него поправку, которая зависит от размера и формы области, — это метод MBH (по фамилиям авторов — Мэттьюза, Бронса и Хазебрука) и его упрощения.

Не үшін қажет

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

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

  • В замкнутой области КВД со временем выполаживается и выходит на среднее давление pˉ\bar p — то, что было в области в момент остановки. Ждать для этого приходится столько же, сколько скважине после пуска нужно для выхода на псевдоустановившийся режим, — часто недели.
  • Продолжение прямой Хорнера p∗p^* здесь не равно ни pip_i, ни pˉ\bar p. Обычно p∗>pˉp^* > \bar p, но у скважины рядом с границей бывает и наоборот.
  • Метод MBH: pˉ=p∗−(m/2,303) pDMBH(tpDA)\bar p = p^* - (m/2{,}303)\,p_{DMBH}(t_{pDA}). Функция pDMBHp_{DMBH} зависит только от формы области и положения скважины; её рассчитывают заранее (графики MBH).
  • Если до остановки скважина вышла на псевдоустановившийся режим, pDMBH=ln⁡(CAtpDA)p_{DMBH} = \ln(C_At_{pDA}) и pˉ=p∗−mlg⁡(CAtpDA)\bar p = p^* - m\lg(C_At_{pDA}). То же дают формулы Дитца и Рэми — Кобба.
  • Если скважина до остановки вышла на режим, то какое время работы подставить в функцию Хорнера, на pˉ\bar p не влияет — лишь бы p∗p^* и tpDAt_{pDA} были посчитаны с одним и тем же временем; прямая же сохраняется намного дольше, если взять время, при котором tpDA≈0,08t_{pDA} \approx 0{,}08.
  • Без площади и формы области поправку не сделать. Их устанавливают по первому длительному исследованию, геологии и расстановке скважин, а дальше контролируют давление короткими КВД.
Негізгі қатынастар
pDMBH=4πkh (p∗−pˉ)qBμ=2,303 (p∗−pˉ)mp_{DMBH} = \dfrac{4\pi kh\,(p^* - \bar p)}{qB\mu} = \dfrac{2{,}303\,(p^* - \bar p)}{m}

безразмерная поправка MBH

pˉ=p∗−m lg⁡(CA tpDA)\bar p = p^* - m\,\lg\left(C_A\,t_{pDA}\right)

скважина до остановки работала в псевдоустановившемся режиме

tpDA=3,553⋅10−3 k tpϕμctAt_{pDA} = \dfrac{3{,}553\cdot10^{-3}\,k\,t_p}{\phi\mu c_t A}

метрические единицы: мД, ч, мПа·с, 1/МПа, м²

Δtpˉ=ϕμctA3,553⋅10−3 k CA\Delta t_{\bar p} = \dfrac{\phi\mu c_t A}{3{,}553\cdot10^{-3}\,k\,C_A}

время Дитца, ч: на прямой MDH (давление от lg⁡Δt\lg\Delta t) в этот момент давление равно pˉ\bar p

pˉ=pi−B NpctVp\bar p = p_i - \dfrac{B\,N_p}{c_t V_p}

материальный баланс замкнутой области (однофазная нефть)

Жиі қателіктер
  • Принимать p∗p^* за среднее давление в разрабатываемой залежи.
  • Брать для поправки коэффициент формы круга 31,6, когда скважина стоит у разлома или в вытянутом блоке.
  • Подставлять ln⁡(CAtpDA)\ln(C_At_{pDA}), когда скважина до остановки ещё не вышла на псевдоустановившийся режим, — тогда нужен график MBH.
  • Считать, что давление «почти выровнялось», раз за последние часы оно росло медленно: на полулогарифмической шкале рост может продолжаться ещё недели.
  • Применять длинную КВД к скважине среди работающих соседей без учёта их влияния.

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

  • объяснить, почему в замкнутой области давление после остановки выравнивается к среднему пластовому, а не к начальному, и сколько времени это занимает;
  • показать, почему продолжение прямой Хорнера p∗p^* в ограниченной области не равно ни начальному, ни среднему давлению;
  • рассчитать среднее давление методом MBH по безразмерному времени tpDAt_{pDA} и форме области, а для скважины, вышедшей на псевдоустановившийся режим, — по формулам Дитца и Рэми — Кобба;
  • выбрать время работы для функции Хорнера и оценить, как оно влияет на длину прямой и на результат;
  • найти среднее давление по поздним точкам длинной КВД и построить по средним давлениям повторных исследований материальный баланс области.

Какое давление называют средним

Пока скважина работает, давление в области дренирования неодинаково: у ствола — воронка депрессии, вдали давление выше. Среднее пластовое давление pˉ\bar p — среднее по объёму пор. Если область замкнута и отбор прекратить, флюид перетечёт из мест с высоким давлением в места с низким, и давление выровняется. При постоянной сжимаемости выровнявшееся значение равно именно pˉ\bar p: сколько флюида расширилось там, где давление упало, столько же сжалось там, где оно выросло, и общее количество флюида в области не изменилось.

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

В разбуренной залежи остановка одной скважины не останавливает соседей: давление в её области сначала растёт, а потом начинает снижаться вместе со всей залежью. Поэтому КВД относят к области дренирования скважины такой, какой она была в момент остановки. Метод MBH пользуется только прямой радиального притока — ранней частью КВД, которую соседи ещё не исказили, — и потому применим и в разбуренной залежи; длинные КВД среди работающих скважин требуют учёта соседей (Исследования в разбуренной залежи: метод Слайдера). Примеры урока проще: скважина одна в замкнутом блоке.

Давление вдоль длинной оси блока в момент остановки скважины и через 24, 200 и 1000 часов после неё02004006008001000161718192021Расстояние вдоль оси блока, мДавление в пласте, МПау стенки скважины 16,6 МПаостановка и 24 ч200 ч1000 чпунктир — p̄ = 20,14 МПаштрих — через 24 ч:воронка заполнилась

График давления от 16 до 21,5 МПа вдоль оси блока от 0 до 1000 метров. Кривая момента остановки опускается в узкой воронке у скважины на 250-м метре до 16,6 МПа и поднимается к дальнему концу до 20,9 МПа. Штриховая кривая через 24 часа совпадает с ней вдали от скважины, но воронки уже нет: давление вдоль оси от 19,5 до 20,8 МПа. Через 200 часов — от 19,7 до 20,6 МПа, через 1000 часов почти ровная линия у 20,1–20,2 МПа рядом с пунктиром среднего давления 20,14 МПа.

ЕсептеуРасчётное давление в пласте вдоль длинной оси блока из примеров урока (1000 × 500 м, скважина в 250 м от короткой стороны) в момент остановки после 4000 ч работы и через 24, 200 и 1000 ч после неё; пунктир — среднее давление области. Перепад на скине не показан: в стволе в момент остановки давление ещё на 0,8 МПа ниже — 15,7 МПа. Воронка у скважины заполняется за сутки, а дальний конец блока выравнивается с остальным неделями.

Восстановление давления в замкнутой области

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

КВД в замкнутой области как сумма двух дебитов
pi−pws(Δt)=qBμ2πkh[pD(tp+Δt)−pD(Δt)]p_i - p_{ws}(\Delta t) = \frac{qB\mu}{2\pi kh}\Big[p_D(t_p + \Delta t) - p_D(\Delta t)\Big]
мұндағы:
  • pip_iначальное пластовое давление, Па
  • pws(Δt)p_{ws}(\Delta t)забойное давление через время Δt\Delta t после остановки, Па
  • pD(t)p_D(t)безразмерное падение давления на стенке скважины, работающей с постоянным дебитом время tt в замкнутой области, без скина (скин входит в оба слагаемых одинаково и сокращается), —
  • tpt_pвремя работы с дебитом qq до остановки, с
  • q, B, μq,\ B,\ \muдебит в стандартных условиях, объёмный коэффициент, вязкость, м³/с, м³/м³, Па·с
  • k, hk,\ hпроницаемость и эффективная толщина, м², м

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

Мына жағдайда орындалады
  • однофазный слабосжимаемый флюид, однородный пласт, область замкнута;
  • постоянный дебит до остановки (при меняющемся — эквивалентное время Np/qN_p/q);
  • влияние ствола не учитываем: скважину закрывают у забоя или ранние точки не используют.

От формулы урока Анализ кривой восстановления давления она отличается только функцией pDp_D. Пока возмущения не дошли до границ, это логарифм, и получается уравнение Хорнера. Потом pDp_D начинает расти линейно во времени — область истощается.

Когда и время работы, и время остановки больше времени выхода на псевдоустановившийся режим, оба значения pDp_D равны 2πtDA+12ln⁡(4A/(γCArw2))2\pi t_{DA} + \tfrac12\ln\big(4A/(\gamma C_Ar_w^2)\big) (забойное давление через начальное). Постоянные сокращаются, и остаётся только истощение:

Давление, на которое выходит КВД в замкнутой области
pws(Δt→∞)=pi−2πtpDA qBμ2πkh=pi−qB tpctVp=pˉp_{ws}(\Delta t \to \infty) = p_i - 2\pi t_{pDA}\,\frac{qB\mu}{2\pi kh} = p_i - \frac{qB\,t_p}{c_t V_p} = \bar p
мұндағы:
  • tpDAt_{pDA}время работы, отнесённое к площади: ktp/(ϕμctA)kt_p/(\phi\mu c_tA), —
  • AAплощадь области дренирования, м²
  • Vp=ϕhAV_p = \phi hAпоровый объём области, м³
  • ctc_tобщая сжимаемость пород и флюидов, 1/Па
  • pˉ\bar pсреднее давление области в момент остановки, Па

Бірліктер: СИ; в метрических единицах снижение давления равно qBtp/(24 ctVp)qBt_p/(24\,c_tV_p) при qq в м³/сут, tpt_p в часах и ctc_t в 1/МПа

Мына жағдайда орындалады
  • ΔtDA=kΔt/(ϕμctA)\Delta t_{DA} = k\Delta t/(\phi\mu c_tA) больше времени выхода на режим для данной формы области;
  • область замкнута, соседние скважины не работают, флюид однофазный.

КВД выходит на давление, которое материальный баланс приписывает области в момент остановки, то есть на среднее пластовое давление. Ждать приходится столько же, сколько скважине после пуска нужно для выхода на псевдоустановившийся режим: в равенстве участвует то же pD(Δt)p_D(\Delta t). Для компактной области со скважиной в центре это tDAt_{DA} около 0,1–0,15, для вытянутой или со скважиной у края — около единицы.

КВД из примеров урока: как давление подходит к среднему (расчёт модели)
Δt\Delta t, ч
24
Δt\Delta t, сут
1,0
pwsp_{ws}, МПа
19,515
pˉ−pws\bar p - p_{ws}, МПа
0,625
Δt\Delta t, ч
72
Δt\Delta t, сут
3,0
pwsp_{ws}, МПа
19,683
pˉ−pws\bar p - p_{ws}, МПа
0,458
Δt\Delta t, ч
200
Δt\Delta t, сут
8,3
pwsp_{ws}, МПа
19,831
pˉ−pws\bar p - p_{ws}, МПа
0,310
Δt\Delta t, ч
400
Δt\Delta t, сут
17
pwsp_{ws}, МПа
19,962
pˉ−pws\bar p - p_{ws}, МПа
0,179
Δt\Delta t, ч
800
Δt\Delta t, сут
33
pwsp_{ws}, МПа
20,081
pˉ−pws\bar p - p_{ws}, МПа
0,060
Δt\Delta t, ч
1600
Δt\Delta t, сут
67
pwsp_{ws}, МПа
20,134
pˉ−pws\bar p - p_{ws}, МПа
0,007

Блок 1000 × 500 м, скважина на длинной оси в 250 м от короткой стороны, 4000 ч работы с дебитом 45 м³/сут; среднее давление области в момент остановки — 20,141 МПа, начальное — 26,0 МПа. Данные условные.

На графике Хорнера такая КВД выглядит так: после влияния ствола — прямая радиального притока, затем точки отходят от неё вниз и выполаживаются к pˉ\bar p. Продолжение прямой до отношения Хорнера, равного 1, даёт p∗p^*, и здесь оно не равно ни начальному давлению (его в области давно нет), ни среднему. Уравнение Хорнера построено для бесконечного пласта, где бесконечно долгая остановка вернула бы давление к pip_i; в замкнутой области это не так.

График Хорнера КВД в замкнутом блоке: замеры, продолжение по модели, прямая радиального притока и среднее давление100 00010 000100010010118,519,019,520,020,5Отношение Хорнера (tp + Δt)/ΔtДавление pws, МПаp* = 20,61замеры 0,25–72 ч — точки;продолжение КВД по модели — линияp̄ = 20,14 МПа

График: по горизонтали отношение Хорнера от ста тысяч до единицы на логарифмической шкале, по вертикали давление от 18,4 до 20,8 МПа. Десять точек замеров поднимаются от 18,6 до 19,7 МПа; точки до суток лежат на штриховой прямой, а при 48 и 72 часах уже немного ниже неё. Прямая продолжена до отношения 1, где давление p* около 20,61 МПа. Кривая модели уходит под прямую и выполаживается к пунктиру среднего давления 20,14 МПа при отношениях Хорнера меньше 3.

ЕсептеуГрафик Хорнера КВД из примеров урока: точки — замеры за первые 72 ч, тонкая линия — продолжение той же КВД по модели до 3000 ч, штриховая прямая — прямая радиального притока и её продолжение до p∗p^*, пунктир — среднее давление. Через 72 ч давление ещё почти на 0,5 МПа ниже среднего.

Метод MBH: поправка от p* к среднему давлению

Прямая радиального притока на КВД в замкнутой области сохраняется: в первые часы после остановки давление растёт почти так же, как в бесконечном пласте (оговорка — в разделе о времени работы). Значит, по ней по-прежнему находят khkh, скин и p∗p^*. Остаётся связать p∗p^* со средним давлением. Разность между ними зависит только от формы области, положения в ней скважины и безразмерного времени работы — её удобно выразить безразмерно и рассчитать один раз для каждой геометрии. Это функция MBH (метод MBH):

Безразмерная поправка MBH
pDMBH(tpDA)=4πkhqBμ (p∗−pˉ)=2,303 (p∗−pˉ)mp_{DMBH}(t_{pDA}) = \frac{4\pi kh}{qB\mu}\,(p^* - \bar p) = \frac{2{,}303\,(p^* - \bar p)}{m}
мұндағы:
  • p∗p^*давление на продолжении прямой Хорнера, построенной с фактическим временем работы tpt_p, при отношении Хорнера 1, Па
  • pˉ\bar pсреднее давление области в момент остановки, Па
  • mmнаклон прямой Хорнера на десятичный цикл (положительный), Па
  • tpDAt_{pDA}время работы, отнесённое к площади: ktp/(ϕμctA)kt_p/(\phi\mu c_tA), —
  • q, B, μ, k, hq,\ B,\ \mu,\ k,\ hдебит до остановки, объёмный коэффициент, вязкость, проницаемость, толщина, м³/с, м³/м³, Па·с, м², м

Бірліктер: СИ; вторая форма верна в любых единицах, если p∗p^*, pˉ\bar p и mm выражены в одних единицах давления

Мына жағдайда орындалады
  • на КВД есть прямая радиального притока (её подтверждают производной давления);
  • однородный пласт, однофазный флюид, область замкнута, скважина в ней одна (или область — часть блока, отделённая от соседних скважин линиями раздела потоков);
  • до остановки дебит был постоянным; на участке прямой Δt≪tp\Delta t \ll t_p.

Отсюда рабочая формула: pˉ=p∗−(m/2,303) pDMBH(tpDA)\bar p = p^* - (m/2{,}303)\,p_{DMBH}(t_{pDA}). Значение pDMBHp_{DMBH} берут с графика MBH для подходящей формы области и положения скважины.

Басқа бірліктер жүйелерінде (2)
Метрические промысловые единицы
p∗−pˉ=0,933 q B μk h pDMBHp^* - \bar p = 0{,}933\,\frac{q\,B\,\mu}{k\,h}\,p_{DMBH}

qq — м³/сут в стандартных условиях, μ\mu — мПа·с, kk — мД, hh — м; давление — МПа.

0,933=1/(4π⋅0,0853)=2,149/2,3030{,}933 = 1/(4\pi\cdot0{,}0853) = 2{,}149/2{,}303 — наклон на единицу натурального логарифма (откуда постоянные).

Промысловые единицы (oilfield)
p∗−pˉ=70,6 q B μk h pDMBHp^* - \bar p = 70{,}6\,\frac{q\,B\,\mu}{k\,h}\,p_{DMBH}

qq — STB/d, BB — rb/STB, μ\mu — сП, kk — мД, hh — ft; давление — psi.

70,6=162,6/2,30370{,}6 = 162{,}6/2{,}303.

Откуда берётся эта функция, показывает суперпозиция (вывод — в дополнительных разделах):

Функция MBH через решение для работающей скважины
pDMBH(tpDA)=4πtpDA−2[pD(tpD)−12(ln⁡tpD+0,809)]p_{DMBH}(t_{pDA}) = 4\pi t_{pDA} - 2\Big[p_D(t_{pD}) - \tfrac12\big(\ln t_{pD} + 0{,}809\big)\Big]
мұндағы:
  • pD(tpD)p_D(t_{pD})функция давления скважины в замкнутой области при постоянном дебите (без скина) в момент остановки, —
  • 12(ln⁡tpD+0,809)\tfrac12(\ln t_{pD} + 0{,}809)то же для бесконечного пласта (Решение для линейного источника и функция давления), —
  • tpDt_{pD}время работы в безразмерном виде, ktp/(ϕμctrw2)kt_p/(\phi\mu c_tr_w^2), —

Бірліктер: безразмерная; ln⁡\ln — натуральный логарифм

Мына жағдайда орындалады
  • те же, что для функции MBH; pDp_D рассчитан для конкретной формы области и положения скважины.

Первое слагаемое — истощение по материальному балансу, pi−pˉp_i - \bar p в единицах qBμ/(4πkh)qB\mu/(4\pi kh). Второе — насколько границы успели понизить забойное давление к моменту остановки по сравнению с бесконечным пластом: ровно на столько же p∗p^* ниже pip_i. Если границы ещё не сказались, второе слагаемое равно нулю, p∗=pip^* = p_i и pDMBH=4πtpDAp_{DMBH} = 4\pi t_{pDA} — поправка сводится к материальному балансу. В псевдоустановившемся режиме pDMBH=ln⁡(CAtpDA)p_{DMBH} = \ln(C_At_{pDA}) — прямая на полулогарифмическом графике; при tpDA=1t_{pDA} = 1 её ордината равна ln⁡CA\ln C_A.

Функции MBH для скважины в замкнутом квадрате и прямоугольниках0,010,1110−10123456Время работы tpDApDMBH4πtквадрат, центр2:1, центр4:1, центр2:1, ¼ длины

График: по горизонтали безразмерное время работы от 0,01 до 10 на логарифмической шкале, по вертикали функция MBH от −1 до 6. Штриховая кривая 4πt поднимается круто. Четыре сплошные кривые сначала идут по ней, затем отходят: квадрат с центральной скважиной выходит на прямую около 0,1 и поднимается до 5,7 при 10; прямоугольник 2:1 с центральной скважиной идёт немного ниже; прямоугольник 4:1 и прямоугольник 2:1 со скважиной в четверти длины проходят через плато около 0,3–0,9 при времени 0,05–0,3 и выходят на свои прямые около 0,5–1, достигая 4 и 3,8 при 10.

ЕсептеуРасчётные функции MBH для скважины в замкнутом квадрате и прямоугольниках (модель урока); оранжевая кривая — блок примеров урока (2:1, скважина в четверти длины). Штриховая кривая — 4πtpDA4\pi t_{pDA}: границы ещё не повлияли на работу скважины. Каждая кривая со временем выходит на свою прямую ln⁡(CAtpDA)\ln(C_At_{pDA}). Чем вытянутее область и чем ближе скважина к её границе, тем ниже кривая и тем позже она выходит на прямую.
Функция MBH для некоторых замкнутых областей (расчёт модели)
Форма и положение скважины
Квадрат, скважина в центре
CAC_A
30,88
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
0,11
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
1,15
при tpDAt_{pDA} = 0,3
2,23
Форма и положение скважины
Прямоугольник 2:1, в центре
CAC_A
21,84
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
0,21
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
0,96
при tpDAt_{pDA} = 0,3
1,88
Форма и положение скважины
Прямоугольник 2:1, на оси в 1/4 длины от короткой стороны
CAC_A
4,52
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
0,98
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
0,76
при tpDAt_{pDA} = 0,3
0,88
Форма и положение скважины
Прямоугольник 2:1, на оси в 1/8 длины от короткой стороны
CAC_A
0,58
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
1,09
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
−0,10
при tpDAt_{pDA} = 0,3
−0,75
Форма и положение скважины
Прямоугольник 4:1, в центре
CAC_A
5,38
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
0,49
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
0,34
при tpDAt_{pDA} = 0,3
0,61
Форма и положение скважины
Прямоугольник 5:1, в центре
CAC_A
2,36
ln⁡(CAtpDA)\ln(C_At_{pDA}) годится при tpDA>t_{pDA} >
0,64
pDMBHp_{DMBH} при tpDAt_{pDA} = 0,1
0,04
при tpDAt_{pDA} = 0,3
−0,05

Скважина — линейный источник в однородном замкнутом прямоугольнике; «годится» — по критерию урока: отличие от ln⁡(CAtpDA)\ln(C_At_{pDA}) меньше 0,02 (ошибка pˉ\bar p меньше 0,01 m0{,}01\,m). Классические табличные CAC_A совпадают с этими в пределах 0,4 %. Отрицательные значения означают, что pˉ>p∗\bar p > p^*.

Порядок расчёта по методу MBH
  1. 1
    Прямая Хорнера

    По точкам радиального притока находят наклон mm, проводимость khkh и p∗p^*, как в уроке Анализ кривой восстановления давления; в функцию Хорнера подставляют время работы tpt_p или Np/qN_p/q.

  2. 2
    Площадь и форма

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

  3. 3
    Безразмерное время

    tpDA=3,553⋅10−3 ktp/(ϕμctA)t_{pDA} = 3{,}553\cdot10^{-3}\,kt_p/(\phi\mu c_tA) при tpt_p в часах, kk в мД, μ\mu в мПа·с, ctc_t в 1/МПа и AA в м².

  4. 4
    Поправка

    С графика MBH для подходящей формы снимают pDMBH(tpDA)p_{DMBH}(t_{pDA}); если tpDAt_{pDA} больше времени, с которого кривая совпадает со своей прямой, берут ln⁡(CAtpDA)\ln(C_At_{pDA}).

  5. 5
    Среднее давление

    pˉ=p∗−(m/2,303) pDMBH\bar p = p^* - (m/2{,}303)\,p_{DMBH}. Результат сверяют с материальным балансом и с предыдущими исследованиями.

КВД добывающей скважины в замкнутом блоке (условные данные)
Δt\Delta t, ч
0,25
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
16 001
pwsp_{ws}, МПа
18,592
Δt\Delta t, ч
0,5
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
8001
pwsp_{ws}, МПа
18,737
Δt\Delta t, ч
1
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
4001
pwsp_{ws}, МПа
18,882
Δt\Delta t, ч
2
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
2001
pwsp_{ws}, МПа
19,026
Δt\Delta t, ч
4
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
1001
pwsp_{ws}, МПа
19,169
Δt\Delta t, ч
8
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
501,0
pwsp_{ws}, МПа
19,308
Δt\Delta t, ч
12
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
334,3
pwsp_{ws}, МПа
19,388
Δt\Delta t, ч
24
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
167,7
pwsp_{ws}, МПа
19,515
Δt\Delta t, ч
48
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
84,3
pwsp_{ws}, МПа
19,627
Δt\Delta t, ч
72
Отношение Хорнера (tp+Δt)/Δt(t_p + \Delta t)/\Delta t
56,6
pwsp_{ws}, МПа
19,683

Время работы tpt_p = 4000 ч с дебитом 45 м³/сут; скважину закрыли у забоя, влияние ствола закончилось в первые минуты. Данные рассчитаны по модели однородного замкнутого блока 1000 × 500 м и округлены до 0,001 МПа.

Есептеу мысалыСреднее давление по КВД методом MBH
Бастапқы деректер (шартты)
Скважина и блок
единственная скважина в замкнутом тектоническом блоке 1000 × 500 м (AA = 50 га) на длинной оси в 250 м от короткой стороны; CAC_A = 4,52
Работа до остановки
45 м³/сут нефти; накоплено 7500 м³, tp=Np/qt_p = N_p/q = 4000 ч
BB; μ\mu; hh; ϕ\phi; ctc_t
1,25 м³/м³; 1,2 мПа·с; 10 м; 0,20; 1,6⋅10−31{,}6\cdot10^{-3} МПа⁻¹
Начальное давление
26,0 МПа
КВД
таблица выше; прямая — по замерам 0,5 и 2 ч
Шешім
  1. Абсциссы Хорнера: lg⁡(4000,5/0,5)=3,90314\lg(4000{,}5/0{,}5) = 3{,}90314 и lg⁡(4002/2)=3,30125\lg(4002/2) = 3{,}30125; наклон m=19,026−18,7373,90314−3,30125=0,4802m = \dfrac{19{,}026 - 18{,}737}{3{,}90314 - 3{,}30125} = 0{,}4802 МПа на цикл.
  2. Проводимость и проницаемость: kh=2,149⋅45⋅1,25⋅1,20,4802=302kh = 2{,}149\cdot\dfrac{45\cdot1{,}25\cdot1{,}2}{0{,}4802} = 302 мД·м, k=30,2k = 30{,}2 мД.
  3. Экстраполированное давление: p∗=19,026+0,4802⋅3,30125=20,611p^* = 19{,}026 + 0{,}4802\cdot3{,}30125 = 20{,}611 МПа.
  4. Безразмерное время работы: tpDA=3,553⋅10−3⋅30,2⋅40000,20⋅1,2⋅1,6⋅10−3⋅500 000=2,235t_{pDA} = \dfrac{3{,}553\cdot10^{-3}\cdot30{,}2\cdot4000}{0{,}20\cdot1{,}2\cdot1{,}6\cdot10^{-3}\cdot500\,000} = 2{,}235. Для этой формы функция MBH совпадает со своей прямой при tpDA>0,98t_{pDA} > 0{,}98 (таблица), поэтому pDMBH=ln⁡(CAtpDA)p_{DMBH} = \ln(C_At_{pDA}).
  5. Поправка: CAtpDA=4,52⋅2,235=10,10C_At_{pDA} = 4{,}52\cdot2{,}235 = 10{,}10, lg⁡10,10=1,0043\lg 10{,}10 = 1{,}0043, mlg⁡(CAtpDA)=0,4823m\lg(C_At_{pDA}) = 0{,}4823 МПа.
  6. Среднее давление: pˉ=20,611−0,4823=20,129\bar p = 20{,}611 - 0{,}4823 = 20{,}129 МПа.
  7. Проверка материальным балансом: поровый объём блока 500 000⋅10⋅0,20=106500\,000\cdot10\cdot0{,}20 = 10^6 м³, pˉ=26,0−1,25⋅75001,6⋅10−3⋅106=20,141\bar p = 26{,}0 - \dfrac{1{,}25\cdot7500}{1{,}6\cdot10^{-3}\cdot10^6} = 20{,}141 МПа.
Жауап: Среднее давление в блоке на момент остановки — около 20,129 МПа; p∗p^* выше него почти на полмегапаскаля.

Бұл нені білдіреді. Метод MBH дал среднее давление с ошибкой 0,012 МПа — около 0,2 % от снижения давления в блоке (5,86 МПа). Если бы за среднее приняли p∗p^*, ошибка составила бы 0,47 МПа, и материальный баланс показал бы меньшее истощение, а значит, завышенные запасы. Дождаться выравнивания нельзя: через 72 ч КВД ещё на 0,46 МПа ниже среднего давления, а до разницы 0,05 МПа надо ждать больше месяца (таблица выше). Данные примера условные: они рассчитаны по модели однородного блока.

Если скважина вышла на псевдоустановившийся режим: формулы Дитца и Рэми — Кобба

Добывающие скважины обычно работают перед остановкой месяцами, и их tpDAt_{pDA} намного больше времени выхода на режим. Тогда график MBH не нужен:

Среднее давление после работы в псевдоустановившемся режиме
pˉ=p∗−m lg⁡(CA tpDA)\bar p = p^* - m\,\lg\left(C_A\,t_{pDA}\right)
мұндағы:
  • CAC_Aкоэффициент формы области с данным положением скважины, —
  • tpDAt_{pDA}время работы, отнесённое к площади, —
  • p∗, mp^*,\ mпродолжение прямой Хорнера при отношении 1 и её наклон на десятичный цикл, МПа

Бірліктер: любые согласованные: p∗p^*, pˉ\bar p и mm — в одних единицах давления

Мына жағдайда орындалады
  • до остановки скважина вышла на псевдоустановившийся режим: tpDAt_{pDA} больше времени, с которого функция MBH совпадает с ln⁡(CAtpDA)\ln(C_At_{pDA}) (таблица выше);
  • p∗p^* и tpDAt_{pDA} рассчитаны с одним и тем же временем работы.

При CAtpDA>1C_At_{pDA} > 1 среднее давление ниже p∗p^*, при CAtpDA<1C_At_{pDA} < 1 — выше. Второе бывает у скважины рядом с границей (малый CAC_A) после не слишком долгой работы.

Две графические формы той же формулы избавляют от вычитания поправки:

  • Метод Дитца (метод Дитца). На графике MDH — давление от lg⁡Δt\lg\Delta t — прямую радиального притока продолжают до времени Δtpˉ\Delta t_{\bar p} (формула ниже); давление на прямой в этот момент и есть pˉ\bar p. Время работы знать не нужно, но оно должно быть много больше Δtpˉ\Delta t_{\bar p}.
  • Метод Рэми — Кобба. На графике Хорнера давление читают на прямой при отношении Хорнера (tp+Δt)/Δt=CAtpDA(t_p + \Delta t)/\Delta t = C_At_{pDA}.
Время Дитца
Δtpˉ=ϕμctAk CA\Delta t_{\bar p} = \frac{\phi\mu c_t A}{k\,C_A}
мұндағы:
  • Δtpˉ\Delta t_{\bar p}время на продолжении прямой MDH, при котором она проходит через pˉ\bar p, с
  • A, CAA,\ C_Aплощадь области и коэффициент формы, м², —
  • k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_tпроницаемость, пористость, вязкость, общая сжимаемость, м², —, Па·с, 1/Па

Бірліктер: СИ

Мына жағдайда орындалады
  • до остановки скважина работала в псевдоустановившемся режиме, и время работы много больше Δtpˉ\Delta t_{\bar p};
  • прямая MDH проведена по точкам радиального притока.

В безразмерном виде ΔtDA=1/CA\Delta t_{DA} = 1/C_A. Это не время, через которое скважина действительно выйдет на среднее давление (оно в несколько раз больше), а точка, где продолжение прямой пересекает уровень pˉ\bar p.

Басқа бірліктер жүйелерінде (2)
Метрические промысловые единицы
Δtpˉ [ч]=ϕ  μ [мПа⋅с]  ct [МПа−1]  A [м2]3,553⋅10−3 k [мД]  CA\Delta t_{\bar p}\,[\text{ч}] = \frac{\phi\;\mu\,[\text{мПа·с}]\;c_t\,[\text{МПа}^{-1}]\;A\,[\text{м}^2]}{3{,}553\cdot10^{-3}\,k\,[\text{мД}]\;C_A}

μ\mu — мПа·с, ctc_t — 1/МПа, AA — м², kk — мД; время — часы.

3,553⋅10−33{,}553\cdot10^{-3} — множитель безразмерного времени в этих единицах (время выхода на режим).

Промысловые единицы (oilfield)
Δtpˉ [h]=ϕ  μ [сП]  ct [psi−1]  A [ft2]2,637⋅10−4 k [мД]  CA\Delta t_{\bar p}\,[\text{h}] = \frac{\phi\;\mu\,[\text{сП}]\;c_t\,[\text{psi}^{-1}]\;A\,[\text{ft}^2]}{2{,}637\cdot10^{-4}\,k\,[\text{мД}]\;C_A}

μ\mu — сП, ctc_t — 1/psi, AA — квадратные футы, kk — мД; время — часы.

2,637⋅10−42{,}637\cdot10^{-4} — множитель безразмерного времени в промысловых единицах.

Какое время работы подставлять в функцию Хорнера

В функцию Хорнера подставляют время работы tpt_p, а при меняющемся дебите — эквивалентное время Np/qN_p/q (метод MDH и эквивалентное время). В ограниченной области от этого выбора зависит p∗p^*: если заменить tpt_p на tp′t_p', точки сдвинутся по оси абсцисс, и продолжение прямой изменится на mlg⁡(tp′/tp)m\lg(t_p'/t_p). Но и поправка mlg⁡(CAtpDA)m\lg(C_At_{pDA}) изменится на ту же величину, поэтому среднее давление по формуле предыдущего раздела останется прежним — если скважина вышла на режим и обе величины посчитаны с одним временем.

Время работы влияет на другое — на то, как долго точки лежат на прямой. После долгой работы в псевдоустановившемся режиме к росту давления после остановки примешивается продолжающееся истощение от воображаемой «работающей» скважины суперпозиции, и наклон на графике Хорнера с фактическим tpt_p уменьшается задолго до того, как сказываются границы:

Местный наклон графика Хорнера после работы в псевдоустановившемся режиме
местный наклонm=(1+Δttp′)(1−4π ΔtDA)\frac{\text{местный наклон}}{m} = \left(1 + \frac{\Delta t}{t_p'}\right)\left(1 - 4\pi\,\Delta t_{DA}\right)
мұндағы:
  • tp′t_p'время работы, подставленное в функцию Хорнера (не обязательно фактическое), ч
  • ΔtDA\Delta t_{DA}время остановки, отнесённое к площади: kΔt/(ϕμctA)k\Delta t/(\phi\mu c_tA), —
  • mmнаклон радиального притока, МПа

Бірліктер: безразмерная; Δt\Delta t и tp′t_p' в одних единицах

Мына жағдайда орындалады
  • до остановки скважина работала в псевдоустановившемся режиме;
  • после остановки вокруг скважины ещё радиальный приток: возмущение от остановки не дошло до границ.

С фактическим долгим tpt_p первая скобка почти равна 1, и наклон падает на 5 % уже при ΔtDA≈0,004\Delta t_{DA} \approx 0{,}004. Если же подставить время, при котором tpDA′=1/(4π)≈0,08t_{pDA}' = 1/(4\pi) \approx 0{,}08, искажения первого порядка взаимно гасятся, и прямая держится в несколько раз дольше: до отклонения 5 % — примерно в четыре раза, до 1 % — в десять. Поэтому после долгой работы удобно строить график Хорнера с tpDA′≈0,08t_{pDA}' \approx 0{,}08, а среднее давление считать по той же формуле с CAtpDA′C_At_{pDA}'.

Есептеу мысалыТе же данные: другое время в функции Хорнера и метод Дитца
Бастапқы деректер (шартты)
КВД и пласт
из предыдущего примера; kk = 30,2 мД, CAC_A = 4,52, ϕμctA\phi\mu c_tA = 192 мПа·с·м²/МПа
Время для функции Хорнера
такое, что tpDA′=0,08t_{pDA}' = 0{,}08
Прямая
по замерам 2 ч (19,026 МПа) и 12 ч (19,388 МПа)
Шешім
  1. Время Хорнера: tp′=0,08⋅1923,553⋅10−3⋅30,2=143,1t_p' = \dfrac{0{,}08\cdot192}{3{,}553\cdot10^{-3}\cdot30{,}2} = 143{,}1 ч.
  2. Абсциссы: lg⁡143,1+22=1,86064\lg\dfrac{143{,}1 + 2}{2} = 1{,}86064, lg⁡143,1+1212=1,11143\lg\dfrac{143{,}1 + 12}{12} = 1{,}11143; наклон m′=19,388−19,0261,86064−1,11143=0,4832m' = \dfrac{19{,}388 - 19{,}026}{1{,}86064 - 1{,}11143} = 0{,}4832 МПа на цикл.
  3. Продолжение прямой: p∗′=19,388+0,4832⋅1,11143=19,925p^{*\prime} = 19{,}388 + 0{,}4832\cdot1{,}11143 = 19{,}925 МПа — на 0,686 МПа ниже, чем с фактическим временем.
  4. Среднее давление: CAtpDA′=4,52⋅0,08=0,3616C_At_{pDA}' = 4{,}52\cdot0{,}08 = 0{,}3616, lg⁡0,3616=−0,4418\lg 0{,}3616 = -0{,}4418; pˉ=19,925−0,4832⋅(−0,4418)=20,138\bar p = 19{,}925 - 0{,}4832\cdot(-0{,}4418) = 20{,}138 МПа.
  5. Метод Дитца: Δtpˉ=1923,553⋅10−3⋅30,2⋅4,52=396\Delta t_{\bar p} = \dfrac{192}{3{,}553\cdot10^{-3}\cdot30{,}2\cdot4{,}52} = 396 ч. Прямая MDH через замеры 0,5 и 2 ч: наклон 0,289/lg⁡4=0,48000{,}289/\lg 4 = 0{,}4800 МПа на цикл; на ней при 396 ч давление 19,026+0,4800⋅lg⁡3962=19,026+0,4800⋅2,2967=20,12819{,}026 + 0{,}4800\cdot\lg\dfrac{396}{2} = 19{,}026 + 0{,}4800\cdot2{,}2967 = 20{,}128 МПа.
Жауап: Короткое время Хорнера даёт 20,138 МПа, метод Дитца — 20,128 МПа, метод MBH в предыдущем примере — 20,129 МПа; материальный баланс — 20,141 МПа.

Бұл нені білдіреді. Расхождения — сотые доли мегапаскаля. Точнее всех вариант с коротким временем Хорнера: прямая по замерам 2–12 ч лежит на участке, где наклон искажён меньше чем на 0,7 %, тогда как на графике с фактическим временем работы к 12 ч наклон уже на 8 % меньше. Поэтому с фактическим временем прямую пришлось строить по первым двум часам. Время Дитца — около 17 сут — не момент, когда скважина действительно выйдет на среднее давление: по модели через 396 ч КВД ещё на 0,18 МПа ниже него. Данные примера условные.

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

Поздние точки КВД: подход к среднему давлению

Если КВД всё-таки длинная — скважину остановили надолго по другим причинам, это наблюдательная скважина или пласт очень проницаемый, — среднее давление можно найти и без формы области. В конце восстановления остаток pˉ−pws\bar p - p_{ws} затухает по экспоненте. Давление выравнивается одновременно по разным направлениям и масштабам, но медленнее всего — вдоль самого длинного размера области; этот «тон» и остаётся последним.

Поздняя КВД в замкнутой области
pˉ−pws(Δt)≈a e−λ Δt,λ≈π2ηL2\bar p - p_{ws}(\Delta t) \approx a\,\e^{-\lambda\,\Delta t},\qquad \lambda \approx \frac{\pi^2\eta}{L^2}
мұндағы:
  • aaпостоянная (зависит от формы области и истории работы), Па
  • λ\lambdaпоказатель затухания, 1/с
  • η=k/(ϕμct)\eta = k/(\phi\mu c_t)пьезопроводность пласта, м²/с
  • LLдлина прямоугольной области (длинная сторона), м

Бірліктер: СИ; в метрических единицах η=3,553⋅10−3 k/(ϕμct)\eta = 3{,}553\cdot10^{-3}\,k/(\phi\mu c_t) м²/ч, и λ\lambda получается в 1/ч

Мына жағдайда орындалады
  • поздний период: остальные, более быстрые «тоны» уже затухли — для прямоугольника это ΔtDA\Delta t_{DA} порядка 0,2–0,5 и больше;
  • формула для λ\lambda — для прямоугольной области и скважины не посередине её длины; посередине этот тон не возбуждается, и затухание быстрее. Для скважины в центре круга λ=14,7 η/re2\lambda = 14{,}7\,\eta/r_e^2;
  • одна скважина в замкнутой области, соседи не работают.

По нескольким поздним замерам находят λ\lambda и aa, а с ними и pˉ\bar p. На этом основаны два давних приёма: подобрать pˉ\bar p так, чтобы lg⁡(pˉ−pws)\lg(\bar p - p_{ws}) от Δt\Delta t легло на прямую (метод Маскета для КВД; не путать с одноимённым методом прогноза режима растворённого газа, Прогноз работы залежи с режимом растворённого газа), или построить скорость роста давления  dpws/ dΔt\dd p_{ws}/\dd\Delta t от самого давления — прямую, которая пересекает ось при pws=pˉp_{ws} = \bar p (метод Арпса — Смита).

Есептеу мысалыСреднее давление по трём поздним замерам
Бастапқы деректер (шартты)
Скважина
та же, но остановленная надолго (продолжение КВД по модели)
Замеры
400 ч — 19,962 МПа; 600 ч — 20,037 МПа; 800 ч — 20,081 МПа
Пьезопроводность; длина блока
η\eta = 277,6 м²/ч; 1000 м
Шешім
  1. Приращения за равные промежутки 200 ч: d1=20,037−19,962=0,075d_1 = 20{,}037 - 19{,}962 = 0{,}075 МПа, d2=20,081−20,037=0,044d_2 = 20{,}081 - 20{,}037 = 0{,}044 МПа; их отношение qr=d2/d1=0,5867q_r = d_2/d_1 = 0{,}5867.
  2. Показатель затухания: λ=−ln⁡0,5867/200=0,00267\lambda = -\ln 0{,}5867/200 = 0{,}00267 ч⁻¹; по формуле π2η/L2=9,87⋅277,6/10002=0,00274\pi^2\eta/L^2 = 9{,}87\cdot277{,}6/1000^2 = 0{,}00274 ч⁻¹ — близко.
  3. Остаток роста после последнего замера — сумма убывающей геометрической прогрессии: d2 qr1−qr=0,044⋅0,58671−0,5867=0,062d_2\,\dfrac{q_r}{1 - q_r} = 0{,}044\cdot\dfrac{0{,}5867}{1 - 0{,}5867} = 0{,}062 МПа.
  4. Среднее давление: pˉ=20,081+0,062=20,143\bar p = 20{,}081 + 0{,}062 = 20{,}143 МПа.
Жауап: По поздним точкам pˉ≈20,143\bar p \approx 20{,}143 МПа (материальный баланс — 20,141 МПа).

Бұл нені білдіреді. Приём работает, только когда остался один тон. По замерам 100, 200 и 300 ч (19,727, 19,831 и 19,905 МПа) та же арифметика дала бы 20,088 МПа — на 0,05 МПа меньше. И он требует остановки больше чем на месяц, что для добывающей скважины почти всегда неприемлемо; для наблюдательной или длительно простаивающей скважины — вполне. Данные примера условные.

Повторные исследования и материальный баланс

Первое исследование новой скважины стоит сделать длинным: по позднему участку записи находят границы, площадь и форму области (Исследование на ограниченность пласта, Непроницаемый разлом у скважины: метод изображений), а значит, и CAC_A. Дальше достаточно коротких остановок — только до появления прямой радиального притока: p∗p^* и поправка MBH дают среднее давление на каждую дату, а скважина простаивает часы, а не недели.

Скин по коротким КВД сравнивают при одном и том же наклоне — проводимости из первого, надёжного исследования (если она заметно не меняется): короткая прямая определяет наклон неточно, и его разброс превратился бы в мнимые изменения скина. А средние давления повторных исследований вместе с накопленной добычей дают материальный баланс области:

Среднее давление замкнутой области и накопленная добыча
pˉ=pi−BctVp Np\bar p = p_i - \frac{B}{c_t V_p}\,N_p
мұндағы:
  • NpN_pнакопленная добыча нефти в стандартных условиях к моменту остановки, м³
  • VpV_pпоровый объём области, связанный со скважиной, м³
  • B, ctB,\ c_tобъёмный коэффициент нефти и общая сжимаемость, м³/м³, 1/МПа
  • pip_iначальное пластовое давление, МПа

Бірліктер: любые согласованные (например, м³, 1/МПа и МПа)

Мына жағдайда орындалады
  • замкнутая область без притока воды; однофазная нефть выше давления насыщения; BB и ctc_t постоянны;
  • в области одна скважина, или NpN_p и VpV_p относятся ко всему замкнутому блоку.

Средние давления в зависимости от накопленной добычи ложатся на прямую. Её наклон даёт ctVpc_tV_p, а продолжение к Np=0N_p = 0 должно вернуть начальное давление — это проверка. Запасы нефти в области равны Vp(1−Swi)/BV_p(1 - S_{wi})/B, где SwiS_{wi} — начальная водонасыщенность, а BB взят при начальном давлении (Общее уравнение материального баланса нефтяной залежи).

Есептеу мысалыПоровый объём по трём исследованиям
Бастапқы деректер (шартты)
Исследования той же скважины
после 2000, 3000 и 4000 ч работы с дебитом 45 м³/сут (NpN_p = 3750, 5625 и 7500 м³)
Прямые Хорнера
mm = 0,480 МПа на цикл во всех трёх; p∗p^* = 23,399; 22,016 и 20,611 МПа
CAC_A; BB; ctc_t
4,52; 1,25; 1,6⋅10−31{,}6\cdot10^{-3} МПа⁻¹; tpDAt_{pDA} при 4000 ч — 2,235 (первый пример)
Шешім
  1. Безразмерное время при 2000 и 3000 ч: 2,235⋅2000/4000=1,1182{,}235\cdot2000/4000 = 1{,}118 и 2,235⋅3000/4000=1,6762{,}235\cdot3000/4000 = 1{,}676 — оба больше 0,98, поэтому годится ln⁡(CAtpDA)\ln(C_At_{pDA}).
  2. Средние давления: 23,399−0,480⋅lg⁡(4,52⋅1,118)=23,399−0,480⋅0,7036=23,06123{,}399 - 0{,}480\cdot\lg(4{,}52\cdot1{,}118) = 23{,}399 - 0{,}480\cdot0{,}7036 = 23{,}061 МПа; 22,016−0,480⋅0,8794=21,59422{,}016 - 0{,}480\cdot0{,}8794 = 21{,}594 МПа; 20,611−0,480⋅1,0044=20,12920{,}611 - 0{,}480\cdot1{,}0044 = 20{,}129 МПа.
  3. Наклон по первому и последнему: при отборе 1,25⋅(7500−3750)=4687,51{,}25\cdot(7500 - 3750) = 4687{,}5 м³ в пластовых условиях среднее давление снизилось на 23,061−20,129=2,93223{,}061 - 20{,}129 = 2{,}932 МПа; ctVp=4687,5/2,932=1599c_tV_p = 4687{,}5/2{,}932 = 1599 м³/МПа, Vp=1599/(1,6⋅10−3)=999V_p = 1599/(1{,}6\cdot10^{-3}) = 999 тыс. м³.
  4. Проверка: продолжение к Np=0N_p = 0 — 23,061+2,932=25,99323{,}061 + 2{,}932 = 25{,}993 МПа — совпадает с начальным давлением 26,0 МПа.
  5. Если вместо средних давлений взять сами p∗p^*: ctVp=4687,5/(23,399−20,611)=1681c_tV_p = 4687{,}5/(23{,}399 - 20{,}611) = 1681 м³/МПа — объём на 5 % больше, а продолжение к Np=0N_p = 0 даёт 26,19 МПа — выше начального давления.
Жауап: Поровый объём, связанный со скважиной, — около 999 тыс. м³ (в модели — 1000 тыс. м³); без поправки MBH он завышен на 5 %.

Бұл нені білдіреді. Разность p∗−pˉ=mlg⁡(CAtpDA)p^* - \bar p = m\lg(C_At_{pDA}) растёт со временем работы, поэтому ряд p∗p^* снижается медленнее среднего давления — отсюда завышенный объём и «начальное давление» выше настоящего. Не менее сильно на результат влияет сжимаемость: VpV_p обратно пропорционален ctc_t, а сжимаемость пор известна хуже всего (Сжимаемость пород и влияние напряжений на свойства керна). При притоке воды точки отклоняются от прямой вверх, и простой баланс неприменим (Модели притока законтурной воды). Данные примера условные.

Пределы методов и практика

Способы найти среднее давление по КВД
Способ
MBH по графику
Что нужно знать
площадь, форму, положение скважины, tpt_p
Когда применим
при любом времени работы; нужна прямая Хорнера
Главный риск
неверная площадь или форма
Способ
Формулы Дитца, Рэми — Кобба
Что нужно знать
площадь и CAC_A
Когда применим
скважина до остановки вышла на режим
Главный риск
работа короче времени выхода на режим
Способ
Поздние точки (экспонента)
Что нужно знать
только сами замеры
Когда применим
остановка на ΔtDA\Delta t_{DA} порядка 0,3 и больше
Главный риск
длительность; работа соседей
Способ
Материальный баланс
Что нужно знать
VpV_p, ctc_t, накопленную добычу
Когда применим
замкнутая область, однофазный флюид
Главный риск
приток воды, неверная ctc_t
  • Слоистый пласт. Слои разной проницаемости истощаются по-разному, и в остановленной скважине между ними идёт переток через ствол. Давление такой скважины — не среднее ни одного слоя; пластовое давление по слоям измеряют опробователем пластов (Замеры давления в разрабатываемых пластах: барьеры и точность).
  • Газ. Сжимаемость газа зависит от давления, поэтому методы применяют в псевдодавлениях (псевдодавление), а материальный баланс строят по p/Zp/Z (Материальный баланс газовой залежи).
  • Давление ниже давления насыщения. В пласте свободный газ, сжимаемость и подвижность меняются с давлением; поправки становятся приближёнными.
  • Подпитка. Если к области подходит вода или в ней есть газовая шапка, область не замкнута, и давление после остановки выходит не на баланс замкнутой области (Модели притока законтурной воды).
  • Соседние скважины. При долгой остановке их работа меняет и рост давления, и уровень, к которому оно стремится (Исследования в разбуренной залежи: метод Слайдера).

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

Среднее пластовое давление pˉ\bar p
среднее по объёму пор давление области дренирования; к нему в замкнутой области выходит давление после бесконечно долгой остановки.
Функция MBH pDMBHp_{DMBH}
безразмерная разность p∗−pˉp^* - \bar p в единицах qBμ/(4πkh)qB\mu/(4\pi kh); зависит от формы области, положения скважины и tpDAt_{pDA}.
tpDAt_{pDA}
время работы до остановки, отнесённое к площади области: ktp/(ϕμctA)kt_p/(\phi\mu c_tA).
Время Дитца Δtpˉ\Delta t_{\bar p}
время на продолжении прямой MDH, при котором она проходит через pˉ\bar p; в безразмерном виде ΔtDA=1/CA\Delta t_{DA} = 1/C_A.
Время Хорнера tp′t_p'
время работы, подставленное в функцию Хорнера; может отличаться от фактического, если tpDAt_{pDA} в поправке посчитано с ним же.
Повторные исследования
короткие КВД для контроля пластового давления и скина после первого, длинного исследования, которое установило границы области.

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

Қате

«Давление за последние часы почти не растёт — значит, КВД закончилась, и последний замер равен среднему давлению».

Шын мәнінде

На полулогарифмической шкале рост может идти ещё недели: в примере урока через 72 ч давление на 0,46 МПа ниже среднего. Последний замер — нижняя граница среднего давления, а не оно само.

Қате

«p∗p^* всегда выше среднего давления».

Шын мәнінде

Обычно да, но у скважины рядом с границей (малый CAC_A) после не слишком долгой работы функция MBH отрицательна, и pˉ>p∗\bar p > p^*.

Қате

«Для поправки достаточно взять CAC_A = 31,6».

Шын мәнінде

31,6 — круг со скважиной в центре. В примере урока CAC_A = 4,52, и с 31,6 среднее давление получилось бы на 0,48⋅lg⁡(31,6/4,52)≈0,410{,}48\cdot\lg(31{,}6/4{,}52) \approx 0{,}41 МПа ниже.

Қате

«Раз время работы известно неточно, метод MBH бесполезен».

Шын мәнінде

После выхода на режим результат от выбранного времени не зависит, если p∗p^* и tpDAt_{pDA} посчитаны с одним и тем же временем.

Қате

«Время Дитца — это сколько нужно держать скважину закрытой».

Шын мәнінде

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

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

Вывод функции MBH и её предельных форм

Запишем КВД суперпозицией в безразмерном виде, pDs=2πkh(pi−pws)/(qBμ)=pD(tp+Δt)−pD(Δt)p_{Ds} = 2\pi kh(p_i - p_{ws})/(qB\mu) = p_D(t_p + \Delta t) - p_D(\Delta t). На участке прямой Хорнера после остановки вокруг скважины радиальный приток, pD(Δt)=12(ln⁡ΔtD+0,809)p_D(\Delta t) = \tfrac12(\ln\Delta t_D + 0{,}809), а Δt≪tp\Delta t \ll t_p, так что pD(tp+Δt)≈pD(tp)p_D(t_p + \Delta t) \approx p_D(t_p). Прибавив и вычтя 12ln⁡tpD\tfrac12\ln t_{pD}, получим pDs=12ln⁡(tp/Δt)+pD(tpD)−12(ln⁡tpD+0,809)p_{Ds} = \tfrac12\ln(t_p/\Delta t) + p_D(t_{pD}) - \tfrac12(\ln t_{pD} + 0{,}809): прямая Хорнера с наклоном qBμ/(4πkh)qB\mu/(4\pi kh) на единицу натурального логарифма, продолжение которой (tp/Δt→1t_p/\Delta t \to 1) равно p∗=pi−qBμ2πkh[pD(tpD)−12(ln⁡tpD+0,809)]p^* = p_i - \frac{qB\mu}{2\pi kh}\big[p_D(t_{pD}) - \tfrac12(\ln t_{pD} + 0{,}809)\big].

Среднее давление в момент остановки — материальный баланс: 2πkhqBμ(pi−pˉ)=2πtpDA\frac{2\pi kh}{qB\mu}(p_i - \bar p) = 2\pi t_{pDA}. Вычитая, получаем 4πkhqBμ(p∗−pˉ)=4πtpDA−2[pD(tpD)−12(ln⁡tpD+0,809)]\frac{4\pi kh}{qB\mu}(p^* - \bar p) = 4\pi t_{pDA} - 2\big[p_D(t_{pD}) - \tfrac12(\ln t_{pD} + 0{,}809)\big] — формулу основного раздела. Радиус скважины и скин в неё не входят: они одинаково сдвигают pDp_D и его логарифмическое приближение.

В псевдоустановившемся режиме pD(tpD)=2πtpDA+12ln⁡4AγCArw2p_D(t_{pD}) = 2\pi t_{pDA} + \tfrac12\ln\frac{4A}{\gamma C_Ar_w^2}, а 0,809=ln⁡(4/γ)0{,}809 = \ln(4/\gamma) и tpD rw2/A=tpDAt_{pD}\,r_w^2/A = t_{pDA}. Подстановка даёт pDMBH=ln⁡(CAtpDA)p_{DMBH} = \ln(C_At_{pDA}). Отсюда и формула Рэми — Кобба (на прямой Хорнера ln⁡[(tp+Δt)/Δt]=ln⁡(CAtpDA)\ln[(t_p + \Delta t)/\Delta t] = \ln(C_At_{pDA})), и время Дитца: на прямой MDH p=p∗−qBμ4πkhln⁡(tp/Δt)p = p^* - \frac{qB\mu}{4\pi kh}\ln(t_p/\Delta t) равно pˉ\bar p при tp/Δt=CAtpDAt_p/\Delta t = C_At_{pDA}, то есть при ΔtDA=1/CA\Delta t_{DA} = 1/C_A.

Та же арифметика показывает, сколько ждать выравнивания. Если работа шла в псевдоустановившемся режиме, то pˉ−pws(Δt)=qBμ4πkh[pDMBH(ΔtDA)−ln⁡(CAΔtDA)]\bar p - p_{ws}(\Delta t) = \frac{qB\mu}{4\pi kh}\big[p_{DMBH}(\Delta t_{DA}) - \ln(C_A\Delta t_{DA})\big]: КВД подходит к среднему давлению так же, как функция MBH, рассчитанная для времени остановки, — к своей прямой. В примере урока при 400 ч остановки (ΔtDA=0,22\Delta t_{DA} = 0{,}22) эта разность около 0,85, то есть 0,18 МПа.

Почему прямая Хорнера после долгой работы короче

Пусть до остановки скважина работала в псевдоустановившемся режиме. Тогда pD(tp+Δt)=pD(tp)+2πΔtDAp_D(t_p + \Delta t) = p_D(t_p) + 2\pi\Delta t_{DA} — пока скважина «продолжает работать» в суперпозиции, область истощается дальше. После остановки вокруг скважины радиальный приток: pD(Δt)=12ln⁡ΔtD+constp_D(\Delta t) = \tfrac12\ln\Delta t_D + \text{const}. Значит, давление в остановленной скважине pws=const+meln⁡Δt−2me⋅2πΔtDAp_{ws} = \text{const} + m_e\ln\Delta t - 2m_e\cdot2\pi\Delta t_{DA}, где me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh).

Абсцисса Хорнера с временем tp′t_p' — X=ln⁡[(tp′+Δt)/Δt]X = \ln[(t_p' + \Delta t)/\Delta t], и  dX/ dΔt=−tp′/[Δt(tp′+Δt)]\dd X/\dd\Delta t = -t_p'/[\Delta t(t_p' + \Delta t)]. Производная давления по Δt\Delta t равна me(1−4πΔtDA)/Δtm_e(1 - 4\pi\Delta t_{DA})/\Delta t. Делим одно на другое и получаем наклон − dpws/ dX=me(1+Δt/tp′)(1−4πΔtDA)-\dd p_{ws}/\dd X = m_e(1 + \Delta t/t_p')(1 - 4\pi\Delta t_{DA}). Раскрыв скобки, видим поправку первого порядка ΔtDA(1/tpDA′−4π)\Delta t_{DA}(1/t_{pDA}' - 4\pi): она исчезает при tpDA′=1/(4π)=0,0796t_{pDA}' = 1/(4\pi) = 0{,}0796. Расчёт по точной модели урока подтверждает формулу с точностью до третьего знака, пока возмущение от остановки не дошло до границ.

Как рассчитаны кривые урока

Скважина — линейный источник в однородном замкнутом прямоугольнике L×WL\times W. Функция Грина такой области — произведение функций Грина двух «полос» с непроницаемыми краями: по xx шириной LL и по yy шириной WW. Каждую при малых временах удобно считать как сумму отражённых источников (скважины-отражения), а при больших — как ряд косинусов. Безразмерное давление на стенке скважины — интеграл по времени от их произведения; вклад самой скважины даёт привычное 12E1(rD2/(4tD))\tfrac12E_1(r_D^2/(4t_D)), вклад границ — всё остальное.

В псевдоустановившемся режиме отношение наклона в обычных координатах к предельному равно ΠxΠy\Pi_x\Pi_y, где Πx=1+2∑ne−n2π2tDAW/Lcos⁡2(nπxD)\Pi_x = 1 + 2\sum_n \e^{-n^2\pi^2 t_{DA}W/L}\cos^2(n\pi x_D) и Πy\Pi_y — то же с L/WL/W и yDy_D. Отсюда видно, что выравнивание — сумма затухающих экспонент, и самая медленная из них — вдоль длинной стороны. Расчёт сверен двумя независимыми способами: прямым суммированием экспоненциальных интегралов от отражённых скважин (совпадение до 10−1310^{-13}) и классическими таблицами CAC_A (в пределах 0,4 %). Среднее по площади давление в модели в любой момент равно 2πtDA2\pi t_{DA} — материальный баланс выполняется.

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

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

  1. 1

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

    1. А.К экстраполированному давлению p∗p^*
    2. Б.К среднему давлению блока в момент остановки
    3. В.К давлению на границе блока в момент остановки
    4. Г.К начальному пластовому давлению
    Жауабын көрсету

    Дұрыс жауап: Б. К среднему давлению блока в момент остановки

    После остановки флюид только перераспределяется внутри блока; при постоянной сжимаемости давление выравнивается на среднем по объёму пор значении — том, что даёт материальный баланс к моменту остановки. Начальное давление недостижимо: часть флюида уже добыта. p∗p^* — продолжение прямой, а не уровень выравнивания. Давление на границе в момент остановки выше среднего.

  2. 2

    Скважина в центре замкнутой прямоугольной области 2:1 (CAC_A = 21,8) работала до остановки до tpDAt_{pDA} = 0,9. Прямая Хорнера: mm = 0,30 МПа на цикл, p∗p^* = 18,40 МПа. Найдите среднее давление.

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

    Для центральной скважины в области 2:1 функция MBH совпадает с ln⁡(CAtpDA)\ln(C_At_{pDA}) уже с tpDA≈0,2t_{pDA} \approx 0{,}2, поэтому pˉ=18,40−0,30⋅lg⁡(21,8⋅0,9)=18,40−0,30⋅lg⁡19,6=18,40−0,30⋅1,292≈18,01\bar p = 18{,}40 - 0{,}30\cdot\lg(21{,}8\cdot0{,}9) = 18{,}40 - 0{,}30\cdot\lg 19{,}6 = 18{,}40 - 0{,}30\cdot1{,}292 \approx 18{,}01 МПа.

  3. 3

    Что обязательно нужно знать, кроме самой КВД, чтобы применить метод MBH?

    1. А.Давление насыщения нефти
    2. Б.Площадь области дренирования, её форму и положение скважины
    3. В.Коэффициент влияния ствола
    4. Г.Только наклон прямой и p∗p^*
    Жауабын көрсету

    Дұрыс жауап: Б. Площадь области дренирования, её форму и положение скважины

    Поправка pDMBH(tpDA)p_{DMBH}(t_{pDA}) зависит от формы области и положения скважины, а безразмерное время — от площади. Наклон и p∗p^* дают сами точки КВД; давление насыщения и влияние ствола важны для выбора участка прямой, но не для поправки.

  4. 4

    Почему после долгой работы прямую Хорнера лучше строить с временем, при котором tpDA′≈0,08t_{pDA}' \approx 0{,}08, и почему среднее давление от такой замены не меняется?

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

    С фактическим долгим временем местный наклон равен m(1+Δt/tp)(1−4πΔtDA)≈m(1−4πΔtDA)m(1 + \Delta t/t_p)(1 - 4\pi\Delta t_{DA}) \approx m(1 - 4\pi\Delta t_{DA}) и падает на 5 % уже при ΔtDA≈0,004\Delta t_{DA} \approx 0{,}004. При tpDA′=1/(4π)t_{pDA}' = 1/(4\pi) поправки первого порядка гасятся, и прямая держится в несколько раз дольше. Замена времени сдвигает p∗p^* на mlg⁡(tp′/tp)m\lg(t_p'/t_p), а mlg⁡(CAtpDA)m\lg(C_At_{pDA}) — на ту же величину, поэтому pˉ=p∗−mlg⁡(CAtpDA)\bar p = p^* - m\lg(C_At_{pDA}) не меняется, если скважина до остановки вышла на режим.

  5. 5

    Скважина давно работает в центре квадратной области площадью 40 га (CAC_A = 30,9); kk = 80 мД, ϕμct=4⋅10−4\phi\mu c_t = 4\cdot10^{-4} МПа⁻¹·мПа·с. При каком времени остановки нужно прочитать давление на продолжении прямой MDH, чтобы получить среднее давление?

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

    По формуле Дитца Δtpˉ=ϕμctA/(3,553⋅10−3 kCA)=4⋅10−4⋅400 000/(3,553⋅10−3⋅80⋅30,9)=160/8,78≈18\Delta t_{\bar p} = \phi\mu c_tA/(3{,}553\cdot10^{-3}\,kC_A) = 4\cdot10^{-4}\cdot400\,000/(3{,}553\cdot10^{-3}\cdot80\cdot30{,}9) = 160/8{,}78 \approx 18 ч. Сама скважина к этому времени ещё не выйдет на среднее давление — это лишь точка на продолжении прямой.

  6. 6

    При подсчёте запасов по материальному балансу вместо средних давлений повторных КВД взяли значения p∗p^*. Как это скажется на найденном поровом объёме?

    1. А.Изменится только начальное давление
    2. Б.Объём будет завышен
    3. В.Никак
    4. Г.Объём будет занижен
    Жауабын көрсету

    Дұрыс жауап: Б. Объём будет завышен

    p∗−pˉ=mlg⁡(CAtpDA)p^* - \bar p = m\lg(C_At_{pDA}) растёт со временем работы, поэтому ряд p∗p^* снижается медленнее среднего давления. Меньшее снижение давления на ту же добычу означает больший ctVpc_tV_p: в примере урока объём завышен на 5 %, а продолжение к нулевой добыче даёт давление выше начального.

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